Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “waveform”

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 1,063 records · Page 59

Characteristics of Airplanes and Helicopters Recorded by a Dense Seismic Array Near Anza California

We observe frequent air traffic events in continuous seismic waveforms recorded for about 30 days by 1,108 vertical geophones in a tight array on the San Jacinto fault zone. The waveforms of the air traffic events resemble tremor or collections of small earthquakes. However, time-frequency analysis reveals clear Doppler effects that can be modeled with basic equations and fitted well with parameters corresponding to airplanes and helicopters. The flying traces can be inverted by fitting the parameters at each station across the entire array. About 31 air traffic events are detected per day in the relatively remote study area, with peak activity from about 8:00 to 18:00 and significantly fewer events between 23:00 and 4:00 local time. The average event duration is about 200 s, so they cover together>7% of the day. To estimate the total time covered by earthquakes, we derive a scaling relation logτ(M) = 0.41M + 0.89 between earthquake duration and magnitudes using data of 266 earthquakes in the magnitude range 0 ≤ M ≤ 3. The results indicate that in most places, the duration of air traffic events is likely to exceed considerably the total time covered by earthquakes.

58 GEOSCIENCES↗

Computational Seismology Workshop Trains Early-Career Scientists

The Computational Infrastructure for Geodynamics (CIG), in collaboration with Lawrence Livermore National Laboratory (LLNL), held a workshop last year that focused on computational seismology. The workshop provided training in seismic waveform processing, visualization, and high-performance computing (HPC) waveform simulation. Fifty-five predominantly early-career participants—graduate students and postdocs—from the United States and 16 other countries attended. The workshop was the first of its kind to feature full access to HPC resources for research-grade example problems.

58 GEOSCIENCES↗

Seismoacoustic Analysis of Chemical Explosions at the Nevada National Security Site

In recent years, two sets of chemical explosions have been conducted at the Nevada National Security Site (NNSS): 6 explosions from Phase I of the Source Physics Experiment (SPE) and 4 explosions from the Forensics Surface Events (FSE). We have estimated the explosive yield and depth-of-burial/height-of-burst of both from the synthesis of seismic and low-frequency acoustic data. Analysis of the seismic signal is performed using a waveform envelope technique. Using seismic data only, however, there is a tradeoff between the yield and depth or height of the explosion, which becomes especially acute for near-surface events. This tradeoff between yield and depth/height can be broken by including complementary analysis of the acoustic signal using the acoustic waveform inversion method. The acoustic signal also has a strong tradeoff between a yield and depth/height, but the slope is opposite that of the seismic for near-surface events. We use a likelihood function to combine the seismic and acoustic analysis and improve our estimates of yield and depth/height for these sets of explosions. The results are generally consistent with the known parameters, but can be improved upon by higher-frequency seismic analysis and improved acoustic impulse scaling for large scaled depths-of-burial. These events serve as prime examples of data fusion from multiple data streams. The challenge will be to apply the method to events at longer regional distances where the seismic signals will include phases traveling in the upper mantle, and the acoustic signals will be recorded as infrasound signals, which are subject to atmospheric variability.

58 GEOSCIENCES↗

Peak Ground Displacement Saturates Exactly When Expected: Implications for Earthquake Early Warning

The scaling of rupture properties with magnitude is of critical importance to earthquake early warning (EEW) systems that rely on source characterization using limited snapshots of waveform data. ShakeAlert, a prototype EEW system that is being developed for the western United States, provides real-time estimates of earthquake magnitude based on P-wave peak ground displacements calculated at stations triggered by the event. The algorithms used in ShakeAlert assume that the displacement measurements at each station are statistically independent and that there exists a linear and time-independent relation between log peak ground displacement and earthquake magnitude. Here we challenge this basic assumption using the largest dataset assembled for this purpose to date: a comprehensive database of more than 140,000 vertical-component waveforms from M4.5–M9 earthquakes occurring near Japan from 1997 through 2018 and recorded by the K-NET and KiK-net strong-motion networks. By analyzing the time-evolution of P-wave peak ground displacements for these earthquakes, we show that there is a break, or saturation, in the magnitude-displacement scaling that depends on the length of the measurement time window. We demonstrate that the magnitude at which this saturation occurs is well-explained by a simple and non-deterministic model of earthquake rupture growth. We then utilize the predictions of this saturation model to develop a Bayesian framework for estimating posterior uncertainties in real-time magnitude estimates.

58 GEOSCIENCES↗

A Case Study of Nonresonant Mode 3-s ULF Waves Observed by MMS

The nature of the 3-s ultralow frequency (ULF) wave in the Earth's foreshock region and the associated wave-particle interaction are not yet well understood. In this work, we investigate the 3-s ULF waves using Magnetospheric Multiscale (MMS) observations. By combining the plasma rest frame wave properties obtained from multiple methods with the instability analysis based on the velocity distribution in the linear wave stage, the ULF wave is determined to be due to the ion/ion nonresonant mode instability. The interaction between the wave and ions is analyzed using the phase relationship between the transverse wave fields and ion velocities and using the longitudinal momentum equation. During the stage when ULF waves have sinusoidal waveforms up to |dB|/|B 0 | ~ 3, where dB is the wave magnetic field and B0 is the background magnetic field, the wave electric fields perpendicular to B 0 do negative work to solar wind ions; along B 0 , a longitudinal electric field develops, but the V × B force is stronger and leads to solar wind ion deceleration. During the same wave stage, the backstreaming beam ions gain energy from the transverse wave fields and get deceleration along B 0 by the longitudinal electric field. The ULF wave leads to electron heating, preferentially in the direction perpendicular to the local magnetic field. Secondary waves are generated within the ULF waveforms, including whistler waves near half of the electron cyclotron frequency, high-frequency electrostatic waves, and magnetosonic whistler waves. The work improves the understanding of the nature of 3-s ULF waves and the associated wave-particle interaction.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Regional Moment Tensor Inversion Using Rotational Observations

There are benefits from the addition of rotational motions to translational displacements for moment tensor (MT) inversions. The rotational radiation pattern is orthogonal to the shear radiation pattern, thus incorporating rotations is equivalent to gaining another observation point on the focal sphere. We demonstrated this by simulating curl and displacement wavefields for a regional distance station. Thus, one 6-C station (3-Component translational + 3-Component rotational) gathers the same information on radiation pattern as two 3-C stations at 90° azimuth from one another along the focal plane axis, which is sometimes difficult to obtain when restricted to surface sensors. We added rotational Green's functions to a regional MT inversion scheme (long-period, time-domain, and linear inversion) by computing spatial gradients from f-κ reflectivity synthetics. For rotational data, we used Array Derived Rotations (ADR) from Piñon Flats Observatory Array and Golay array deployed during IRIS Community Wavefield Demonstration Experiment. The hope is to ultimately use compact and field deployable broadband rotational seismometers instead of ADRs. Rotational motions were predicted from well-constrained deviatoric MT solutions of nine earthquakes recorded by the arrays and other seismic networks. Additionally, we formed three station sparse datasets with (3-C and 6-C) and without rotational ground motions (3-C only) for eight earthquakes testing the benefits of including rotational waveforms in MT inversions when they are equally weighted with translational waveforms. Adding rotational motions improved the double-couple components and reduced compensated-linear-vector-dipole and isotropic components in full-MT solutions using sparse datasets with poor coverage.

58 GEOSCIENCES↗

Seismic Imaging of the Mw 7.1 Ridgecrest Earthquake Rupture Zone From Data Recorded by Dense Linear Arrays

Abstract We analyze seismograms recorded by four arrays (B1–B4) with 100 m station spacing and apertures of 4–8 km that cross the surface rupture of the 2019 Mw 7.1 Ridgecrest earthquake. The arrays extend from B1 in the northwest to B4 in the southeast of the surface rupture. Delay times between P wave arrivals associated with ∼1,200 local earthquakes and four teleseismic events are used to estimate local velocity variations beneath the arrays. Both teleseismic and local P waves travel faster on the northeast than the southwest side of the fault beneath arrays B1 and B4, but the velocity contrast is less reliably resolved at arrays B2 and B3. We identify several 1–2 km wide low‐velocity zones with much slower inner cores that amplify S waveforms, inferred as damage zones, beneath each array. The damage zones at arrays B2 and B4 also generate fault‐zone head and trapped waves. An automated detector, based on peak ground velocities and durations of high‐amplitude waves, identifies candidate fault‐zone trapped waves (FZTWs) in a localized zone for ∼600 earthquakes at array B4. Synthetic waveform modeling of averaged FZTWs, generated by ∼30 events with high‐quality signals, indicates that the trapping structure at array B4 has a width of ∼300 m, depth of 3–5 km, S wave velocity reduction of ∼20% with respect to the surrounding rock, Q ‐value of ∼30, and S wave velocity contrast of ∼4% across the fault (faster on the northeast side). The results show complex fault‐zone internal structures (velocity contrasts and low‐velocity zones) that vary along fault strike.

Qiu, Hongrui↗

An End-To-End Earthquake Detection Method for Joint Phase Picking and Association Using Deep Learning

Earthquake monitoring by seismic networks typically involves a workflow consisting of phase detection/picking, association, and location tasks. In recent years, the accuracy of these individual stages has been improved through the use of machine learning techniques. Here, in this study, we introduce a new end-to-end approach that improves overall earthquake detection accuracy by jointly optimizing each stage of the detection pipeline. We propose a neural network architecture for the task of multi-station processing of seismic waveforms recorded over a seismic network. This end-to-end architecture consists of three sub-networks: a backbone network that extracts features from raw waveforms, a phase picking sub-network that picks P- and S-wave arrivals based on these features, and an event detection sub-network that aggregates the features from multiple stations to associate and detect earthquakes across a seismic network. We use these sub-networks together with a shift-and-stack module based on back-projection that introduces kinematic constraints on arrival times, allowing the neural network model to generalize to different velocity models and to variable station geometry in seismic networks. We evaluate our proposed method on the STanford EArthquake Dataset (STEAD) and on the 2019 Ridgecrest, CA earthquake sequence. The results demonstrate that our end-to-end approach can effectively pick P- and S-wave arrivals and achieve earthquake detection accuracy rivaling that of other state-of-the-art approaches. Because our approach preserves information across tasks in the detection pipeline, it has the potential to outperform approaches that do not.

58 GEOSCIENCES↗

Identification of Acoustic Wave Signatures in the Ionosphere From Conventional Surface Explosions Using MF/HF Doppler Sounding

Here, we present an experiment to detect one ton TNT-equivalent chemical explosions using pulsed Doppler radar observations of isodensity layers in the ionospheric E region during two campaigns. The first campaign, conducted on 15 October 2019, produced potential detections of all three shots. The detections closely resemble the temporal and spectral properties predicted using the InfraGA ray tracing and weakly nonlinear waveform propagation model. Here the model predicts that within 6.5–7.25 min of each shot a waveform peaking between 0.9 and 0.4 Hz will impact the ionosphere at 100 km. As the waves pass through this region, they will imprint their signal on an isodensity layer, which is detectable using a Doppler radar operating at the plasma frequency of the isodensity. Within the time windows of each of the three shots in the first campaign, we detect enhanced wave activity peaking near 0.5 Hz. These waves were imprinted on the Doppler signal probing an isodensity layer at 2.785 MHz near 100 km altitude. Despite these detections, the method appears to be unreliable as none of the six shots from the second campaign, conducted on 10 July 2020 were detected. The observations from this campaign were characterized by an increased acoustic noise environment in the microbarom band and persistent scintillation on the radar returns. These effects obscured any detectable signal from these shots and the baseline noise was well above the detection levels of the first campaign.

47 OTHER INSTRUMENTATION↗

Combining Deep Learning With Physics Based Features in Explosion-Earthquake Discrimination

This paper combines the power of deep-learning with the generalizability of physics-based features, to present an advanced method for seismic discrimination between earthquakes and explosions. The proposed method contains two branches: a deep learning branch operating directly on seismic waveforms or spectrograms, and a second branch operating on physics-based parametric features. These features are high-frequency P/S amplitude ratios and the difference between local magnitude (ML) and coda duration magnitude (MC). The combination achieves better generalization performance when applied to new regions than models that are developed solely with deep learning. Further, we also examined which parts of the waveform data dominate deep learning decisions (i.e., via Grad-CAM). Such visualization provides a window into the black-box nature of the machine-learning models and offers new insight into how the deep learning derived models use data to make decisions.

58 GEOSCIENCES↗

Spatiotemporal Graph Convolutional Networks for Earthquake Source Characterization

Abstract Accurate earthquake location and magnitude estimation play critical roles in seismology. Recent deep learning frameworks have produced encouraging results on various seismological tasks (e.g., earthquake detection, phase picking, seismic classification, and earthquake early warning). Many existing machine learning earthquake location methods utilize waveform information from a single station. However, multiple stations contain more complete information for earthquake source characterization. Inspired by recent successes in applying graph neural networks (GNNs) in graph‐structured data, we develop a Spatiotemporal Graph Neural Network (STGNN) for estimating earthquake locations and magnitudes. Our graph neural network leverages geographical and waveform information from multiple stations to construct graphs automatically and dynamically by adaptive message passing based on graphs' edges. Using a recent graph neural network and a fully convolutional neural network as baselines, we apply STGNN to earthquakes recorded by the Southern California Seismic Network from 2000 to 2019 and earthquakes collected in Oklahoma from 2014 to 2015. STGNN yields more accurate earthquake locations than those obtained by the baseline models and performs comparably in terms of depth and magnitude prediction, though the ability to predict depth and magnitude remains weak for all tested models. Our work demonstrates the potential of using GNNs and multiple stations for better automatic estimation of earthquake epicenters.

58 GEOSCIENCES↗

Analysis of Narrow Bipolar Events Using Mode Decomposition Methods

Multi-resolution analysis methods can reveal the underlying physical dynamics of nonstationary signals, such as those from lightning. In this paper we demonstrate the application of two multi-resolution analysis methods: Ensemble Empirical Mode Decomposition (EEMD) and Variational Mode Decomposition (VMD) in a comparative way in the analysis of electric field change waveforms from lightning. EEMD and VMD decompose signals into a set of Intrinsic Mode Functions (IMFs). The IMFs can be combined using distance and divergence metrics to obtain noise reduction or to obtain new waveforms that isolate the physical processes of interest while removing irrelevant components of the original signal. We apply the EEMD and VMD methods to the observations of three close Narrow Bipolar Events (NBEs) that were reported by Rison et al. The ΔE observations reveal the occurrence of complex oscillatory processes after the main NBE sferic. We show that both EEMD and VMD are able to isolate the oscillations from the main NBE, with VMD being more effective of the two methods since it requires the least user supervision. The oscillations are found to begin at the end of the NBEs' downward fast positive breakdown and appear to be produced by a half-wavelength standing wave within a weakly-conducting resonant ionization cavity left behind in the wake of the streamer-based NBE event. In conclusion, additional analysis shows that one of the NBEs was likely initiated by an energetic cosmic ray shower, and also corrects a misinterpretation in the literature that fast breakdown is an artifact of NBE-like events in interferometer observations.

54 ENVIRONMENTAL SCIENCES↗

Parameters of the Lightning Attachment Processes in a Negative Cloud-To-Ground Stroke Observed on a Microsecond Timescale

We present time-correlated ultra-high-speed video camera and electromagnetic field measurements of the attachment processes in a natural negative cloud-to-ground stroke. The video camera frame exposure time and pixel resolution were 740 ns and 0.91 m/pixel, respectively. The common streamer zone (CSZ) was first observed 2.52 µs preceding the first frame showing the return stroke (RS) in progress, when the upward and downward leader-tips were 9.8 m apart. In the next frame, the two leaders were observed to have propagated toward each other within the CSZ, with their tips being 0.91 m apart. Our observations show with unprecedented precision/clarity that (a) the slow front in the field waveform is associated with the CSZ, and (b) the “proper” start of the RS is marked by the onset of the fast transition in the field waveform which occurs at the completion of the attachment processes (when the upward and downward leaders have merged).

54 ENVIRONMENTAL SCIENCES↗

Enhancing Regional Seismic Velocity Models With Higher‐Resolution Local Results Using Sparse Dictionary Learning

Abstract We use sparse dictionary learning to develop transformations between seismic velocity models of different resolution and spatial extent. Starting with data in the common region of both models, the method can enhance a regional lower‐resolution model to match the style and resolution of local higher‐resolution results while preserving its regional coverage. The method is demonstrated by applying it to two‐dimensional V S and three‐dimensional V P and V S regional and local velocity models in southern California. The enhanced reconstructed regional results exhibit clear visual improvements, especially in the reconstructed V P / V S ratios, and better correlations with geological features. Moreover, the reconstructed regional V P , V S models outperform the original ones in comparison of simulated earthquake waveforms to observations. The improved fitting to observed waveforms extends beyond the domain of the overlapping region. The developed dictionary learning approach provides physically interpretable results and offers a powerful tool for additional applications of data enhancement in earth sciences.

Geochemistry & Geophysics↗

A Frequency Domain Methodology for Quantitative Evaluation of Diffuse Wavefield With Applications to Seismic Imaging

Abstract Ambient Noise Imaging (ANI) of subsurface structures relies on seismic interferometry of diffuse seismic wavefields. However, the lack of effective methods to quantify and identify highly diffuse waves hampers applications of ANI, particularly in evaluating seismic attenuation and monitoring structural changes with high temporal resolution. Conventional ANI approaches require data normalization, which effectively suppresses the non‐diffuse component with large amplitude but also results in significant loss of amplitude and phase information in the continuous seismic records. In this study, we propose a frequency domain method to quantitatively evaluate the degree of diffuseness of seismic wavefields by analyzing their statistical characteristics of modal amplitudes for stationarity and randomness. Tests on synthetic waveform and field nodal records show that the proposed method can effectively distinguish between diffuse and non‐diffuse waveforms for either single‐ or three‐component data. As an application, we identify a 60‐s‐long diffuse coda of a local M 2.2 earthquake recorded by a dense nodal array on the San Jacinto Fault Zone, and successfully extract high‐quality dispersion curve andQ‐value without performing data normalization. These results are consistent with those obtained by conventional methods that assess the correlation between coherency and the Green's function, and by modeling ballistic waves generated by road traffic. Our proposed method can advance the imaging of subsurface velocity and attenuation structures as well as monitoring temporal changes for scientific studies and engineering applications.

Geochemistry & Geophysics↗

Electro-optic characterization of synthesized infrared-visible light fields

The measurement and control of light field oscillations enable the study of ultrafast phenomena on sub-cycle time scales. Electro-optic sampling (EOS) is a powerful field characterization approach, in terms of both sensitivity and dynamic range, but it has not reached beyond infrared frequencies. Here, we show the synthesis of a sub-cycle infrared-visible pulse and subsequent complete electric field characterization using EOS. The sampled bandwidth spans from 700 nm to 2700 nm (428 to 110 THz). Tailored electric-field waveforms are generated with a two-channel field synthesizer in the infrared-visible range, with a full-width at half-maximum duration as short as 3.8 fs at a central wavelength of 1.7 µm (176 THz). EOS detection of the complete bandwidth of these waveforms extends it into the visible spectral range. To demonstrate the power of our approach, we use the sub-cycle transients to inject carriers in a thin quartz sample for nonlinear photoconductive field sampling with sub-femtosecond resolution.

77 NANOSCIENCE AND NANOTECHNOLOGY↗

Automatic speech recognition predicts contemporaneous earthquake fault displacement

Abstract Significant progress has been made in probing the state of an earthquake fault by applying machine learning to continuous seismic waveforms. The breakthroughs were originally obtained from laboratory shear experiments and numerical simulations of fault shear, then successfully extended to slow-slipping faults. Here we apply the Wav2Vec-2.0 self-supervised framework for automatic speech recognition to continuous seismic signals emanating from a sequence of moderate magnitude earthquakes during the 2018 caldera collapse at the Kīlauea volcano on the island of Hawai’i. We pre-train the Wav2Vec-2.0 model using caldera seismic waveforms and augment the model architecture to predict contemporaneous surface displacement during the caldera collapse sequence, a proxy for fault displacement. We find the model displacement predictions to be excellent. The model is adapted for near-future prediction information and found hints of prediction capability, but the results are not robust. The results demonstrate that earthquake faults emit seismic signatures in a similar manner to laboratory and numerical simulation faults, and artificial intelligence models developed for encoding audio of speech may have important applications in studying active fault zones.

58 GEOSCIENCES↗

Protecting qubit coherence by spectrally engineered driving of the spin environment

Abstract Modern quantum technologies rely crucially on techniques to mitigate quantum decoherence; these techniques can be either passive, achieved for example via materials engineering, or active, typically achieved via pulsed monochromatic driving fields applied to the qubit. Using a solid-state defect spin coupled to a microwave-driven spin bath, we experimentally demonstrate a decoherence mitigation method based on spectral engineering of the environmental noise with a polychromatic drive waveform, and show that it outperforms monochromatic techniques. Results are in agreement with quantitative modeling, and open the path to active decoherence protection using custom-designed waveforms applied to the environment rather than the qubit.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗