SEARCH · Search NASA

Results for “full 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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 235 records · Page 13

All sources NASA NTRS DOE OSTI

DOE OSTI2024

Adjoint waveform tomography of East Asia for improved waveform prediction

We present a preliminary version of the East Asia Tomography (EAT) model, an adjoint waveform tomography model of East and Southeast Asia. We used SPiRaL (Simmons et al., 2021) as our starting model and source parameters for 238 earthquakes from the Global Centroid Moment Tensor catalogue (Ekström et al., 2012). After 50 iterations on Lawrence Livermore National Laboratory’s Lassen supercomputer, we converge on a model with a minimum period of 50 seconds. The preliminary model shows improved slab structure compared to SPiRaL and shows significantly reduced misfit. In later versions of the model, we aim to harness techniques proposed in other studies to improve waveform predictions that travel through the ocean (e.g., Wehner et al., 2022) and use relative amplitude-based misfit functions (e.g., Tao et al., 2018) to better constrain Earth structure. We hope to iterate the current extent of the model to 25 seconds minimum period before iterating to shorter periods for a smaller subregion of the full model.

58 GEOSCIENCES↗

NASA NTRS1979

An acoustical study of the XV-15 Tilt Rotor Research Aircraft

Acoustic data were obtained during a full-scale test of the XV-15 Tilt Rotor Research Aircraft in the Ames 40- by 80-Foot Wind Tunnel. The XV-15 has two 25-ft-diameter, three-bladed rotors at the tips of a 32-ft span wing. The rotors are used as lifting rotors, as propellers, and in various intermediate stages. Acoustic waveforms, dBA, and spectra as functions of different rotor conditions are presented and discussed. The noise level was found to be sensitive to change in rotor-disk angle of attack in the helicopter configuration. Much higher noise levels and harmonic contents were found in the helicopter mode than in the airplane mode. The measurements are compared with existing theoretical predictions which include wind-tunnel wall reflections.

Lee, A.↗

DOE OSTI2020

Effect of net direct current on the properties of radio frequency sheaths: simulation and cross-code comparison

In order to understand, predict and control ion cyclotron range of frequency (ICRF) interactions with tokamak scrape-off layer plasmas, computational tools which can model radio frequency (RF) sheaths are needed. In particular, models for the effective surface impedance and DC rectified sheath potentials may be coupled with full wave RF simulation codes to predict self-consistent wave fields near surfaces and the resulting power dissipation and plasma-material interactions from ion sputtering. In this study, previous work assuming zero net DC current flow through the sheath is generalized to allow the surface to collect net positive or negative current, as is often observed in experiments. The waveforms, DC potential and RF admittance are investigated by means of analytical theory, nonlinear fluid and particle-in-cell (PIC) codes. Cross-code comparisons provide detailed model verification and elucidate the roles of ion and electron kinetics. When the sheath draws negative (positive) DC current, the voltage rectification is reduced (increased) compared with the zero-current case, and the magnitude of both the real and imaginary parts of the admittance are increased (reduced). In conclusion, a previous four-input parametrization of the sheath rectification and admittance properties is generalized to include a fifth parameter describing the DC sheath current.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2025

Improving ICARUS Track Reconstruction Algorithms

The ICARUS experiment is part of the Short-Baseline Neutrino (SBN) program at Fermilab. The main goal of the experiment is to investigate the possibility of sterile neutrinos in the O(1 eV) mass region and provide clarification of the anomaly detected from the Liquid Scintillator Neutrino Detector (LSND) and MiniBooNE experiments. The ICARUS-T600 detector is a Liquid Argon Time Projection Chamber (LAr-TPC), that can provide excellent 3D imaging and calorimetric reconstruction of any ionizing particles. This detection technique allows a detailed study of neutrino interactions, spanning a wide energy spectrum (from a few keV to several hundreds of GeV). The detector consists of two identical adjacent modules, filled with a total of 760 tons of ultra-pure liquid argon. Each module houses two LAr-TPCs separated by a common cathode with a maximum drift distance of 1.5 m, equivalent to about 1 ms drift time for the nominal $500$ V/m electric drift field. The anode is made of three parallel wire planes positioned 3 mm apart, where the stainless-steel wires are oriented on each plane at a different angle with respect to the horizontal direction ($+60^\degree$,$-60^\degree$,$0^\degree$). The first two planes (Induction 1 and Induction 2) provide a non-destructive charge measurement, whereas the ionization charge is fully collected by the last collection plane. In total, 53248 wires with a 3 mm pitch and length up to 9 m are installed in the detector. In the first stage of the reconstruction, segments of waveforms corresponding to physical signals (hits) are searched for in the deconvolved wire waveform with a threshold-based hit-finding algorithm. Each hit is then fitted with a Gaussian, whose area is proportional to the number of drift electrons generating the signal. In the second stage of the reconstruction, hits are passed as input to Pandora, a framework software composed of different pattern recognition algorithms, that performs a 3D reconstruction of the full image recorded in the collected event, including the identification of interaction vertices and tracks and showers inside the TPC. These are organized into a hierarchical structure (called slice) of particles generated starting from a primary interaction vertex. In some cases, related to the inefficiencies in the hit detection or excessive deflection of the particle trajectory, Pandora breaks the particle's track into two or more smaller pieces and considers each piece as an independent track. We studied this phenomenon focusing on primary muons from ν_μ CC interactions contained in a single module with a track at least 20 cm long, to exclude delta rays. The study determined that about $7-8\%$ of the muon tracks are broken. Approximately $80\%$ of the times, Pandora assigns all segments of the track to the same slice (intra-slice track split), while in the remaining $20\%$ of the cases, one of the segments is associated with another slice (extra-slice track split). To mitigate this phenomenon, we designed an algorithm that detects and stitches the tracks broken by Pandora for the intra-slice split. In Monte Carlo simulations, the algorithm showed an efficiency exceeding $80\%$ and a purity exceeding $93\%$.

Ricci, Alessandro Maria [Pisa U.; INFN, Pisa] (ORC↗

DOE OSTI2026

Constructing a High‐Resolution Aftershock Catalog for the 2017 Mw 8.2 Tehuantepec Earthquake Sequence Using a Machine Learning–Based Workflow

The 8 September 2017 Mw 8.2 Tehuantepec earthquake was the largest instrumentally recorded normal‐faulting earthquake in Mexico. The mainshock occurred offshore within the Tehuantepec seismic gap, generating >30,000 aftershocks in the following year. We applied an open‐source, machine learning (ML)–assisted workflow to construct a high‐resolution aftershock catalog using data from temporary and permanent seismic networks in southern Mexico. The workflow integrates PhaseNet for phase detection; GaMMA for phase association; and VELEST, HypoInverse, and HypoDD for velocity modeling and relocation. We processed seven months of continuous waveform data from 29 broadband stations, including a temporary rapid‐response deployment that improved station coverage of the offshore rupture zone. To evaluate performance, we compared our results against analyst‐reviewed picks and event locations from the Servicio Sismológico Nacional catalog. The resulting catalog contains 11,374 relocated earthquakes and represents the most comprehensive published dataset for this sequence, incorporating the first full use of the temporary network. Relocated hypocenters show improved depth control and align well with the Slab2.0 subduction geometry, revealing clearer separation between offshore slab events and onshore crustal seismicity. This study demonstrates that combining ML‐based detection with established methods provides a scalable and reproducible approach for constructing high‐quality earthquake catalogs in tectonically complex environments and offers practical guidance for adapting similar workflows to other earthquake sequences.

Garcia, Marc [The University of Texas at El Paso, ↗

NASA NTRS1975

The high current transient generator, theory and operation for simulating lightning induced voltages into aerospace electrical circuits

Because of the difficulty, hazard, and cost of simulating full-scale lightning currents on flight vehicles, a nondestructive test technique using reduced-scale lightning currents was used to determine if threat level voltages will be induced into aircraft electrical circuits in the event the aircraft is actually struck by lightning. The reasoning and theory of selecting a particular simulated lightning waveshape and magnitude is given, showing that extrapolation of induced waveforms over several orders of magnitude should be avoided. The operation of high-current transient generators and hardware are described along with typical applications. The means by which voltages were generated, induced, and measured, as well as possible sources of error, are discussed.

Schulte, E. H.↗

DOE OSTI2025

Abiotic influences on continuous conifer forest structure across a subalpine watershed

Understanding the abiotic drivers of high-elevation forest physiognomy is essential for forecasting how mountain ecosystems will respond to emerging environmental pressures. Most prior studies of these relationships have relied on small samples of the full landscape, resulting in limited power to detect dominant covariates and their interactions. Here we report the first evaluation of abiotic influences on a complement of accurate, wall-to-wall estimates of conifer forest structure and composition at the watershed scale. In a subalpine conifer domain in the Colorado Rocky Mountains (USA), we developed a novel method for deriving stand structure metrics from waveform LiDAR data, which showed high fidelity with field inventory. We quantified the relationships between structural and compositional metrics and climate, topographic, edaphic, and geologic factors. Our results showed that peak snow water equivalent (SWE), snow disappearance rate, and elevation explained most of the variation in forest structure. The highest stand density, basal area, maximum canopy height, and quadratic mean diameter occurred in sites with SWE around one standard deviation below mean, but with long snow residence times. Stand density decreased linearly with elevation, while other metrics peaked between 3000 m and 3200 m. Substrate properties had weaker influence. Continuous mapping of through-canopy forest structure enabled our novel findings of the dominant role of snowpack in explaining structural and compositional variation, and of elevation thresholds. Our reproducible approach facilitates assessment of forest-topoclimate relationships in other conifer-dominated landscapes and improves understanding of the baseline patterns controlling forest structure, which is needed for predicting long-term ecological change.

forest ecology↗

NASA NTRS2006

Getting a Kick Out of Numerical Relativity

Operating ground-based gravitational wave detectors and a planned instrument in space are bringing about the new field of gravitational wave astronomy. A prime source for any of these observatories is the merger of a system of two black holes. Brought together by copious losses of gravitational-wave energy, these systems merge in a burst of energy with a peak power exceeding any electromagnetic source. Observations of these sources will generate a wealth of astrophysical information, and may provide an unparalleled probe of strong-field gravitational physics, but a full interpretation of the observations will require detailed predictions from General Relativity. I will discuss recent advances in numerical simulations of binary black hole systems which are generating dramatic progress in understanding binary black hole mergers. Recent achievements include the first simulations of binary black hole systems through several orbits and merger, leading to detailed predictions for the final portion of the gravitational radiation waveforms from equal-mass mergers. For unequal-mass mergers, it has recently become possible to measure the impulsive kick imparted to the final black hole, by the asymmetry of the merger radiation. These first results announce an accelerating wave of progress soon to come from the energetic field of numerical relativity.

Source record↗

NASA NTRS2009

Pixelized Device Control Actuators for Large Adaptive Optics

A fully integrated, compact, adaptive space optic mirror assembly has been developed, incorporating new advances in ultralight, high-performance composite mirrors. The composite mirrors use Q-switch matrix architecture-based pixelized control (PMN-PT) actuators, which achieve high-performance, large adaptive optic capability, while reducing the weight of present adaptive optic systems. The self-contained, fully assembled, 11x11x4-in. (approx.= 28x28x10-cm) unit integrates a very-high-performance 8-in. (approx.=20-cm) optic, and has 8-kHz true bandwidth. The assembled unit weighs less than 15 pounds (=6.8 kg), including all mechanical assemblies, power electronics, control electronics, drive electronics, face sheet, wiring, and cabling. It requires just three wires to be attached (power, ground, and signal) for full-function systems integration, and uses a steel-frame and epoxied electronics. The three main innovations are: 1. Ultralightweight composite optics: A new replication method for fabrication of very thin composite 20-cm-diameter laminate face sheets with good as-fabricated optical figure was developed. The approach is a new mandrel resin surface deposition onto previously fabricated thin composite laminates. 2. Matrix (regenerative) power topology: Waveform correction can be achieved across an entire face sheet at 6 kHz, even for large actuator counts. In practice, it was found to be better to develop a quadrant drive, that is, four quadrants of 169 actuators behind the face sheet. Each quadrant has a single, small, regenerative power supply driving all 169 actuators at 8 kHz in effective parallel. 3. Q-switch drive architecture: The Q-switch innovation is at the heart of the matrix architecture, and allows for a very fast current draw into a desired actuator element in 120 counts of a MHz clock without any actuator coupling.

Knowles, Gareth J.↗

NASA NTRS1995

Shuttle Laser Altimeter (SLA): A pathfinder for space-based laser altimetry and lidar

The Shuttle Laser Altimeter (SLA) is a Hitchhiker experiment now being integrated for first flight on STS-72 in November 1995. Four Shuttle flights of the SLA are planned at a rate of about a flight every 18 months. They are aimed at the transition of the Goddard Space Flight Center airborne laser altimeter and lidar technology to low Earth orbit as a pathfinder for operational space-based laser remote sensing devices. Future alser altimeter sensors such as the Geoscience Laser Altimeter System (GLAS), an Earth Observing System facility instrument, and the Multi-Beam Laser Altimeter (MBLA), the land and vegetation laser altimeter for the NASA TOPSAT (Topography Satellite) Mission, will utilize systems and approaches being tested with SLA. The SLA Instrument measures the distance from the Space Shuttle to the Earth's surface by timing the two-way propagation of short (approximately 10 na noseconds) laser pulses. laser pulses at 1064 nm wavelength are generated in a laser transmitter and are detected by a telescope equipped with a silicon avalanche photodiode detector. The SLA data system makes the pulse time interval measurement to a precision of about 10 nsec and also records the temporal shape of the laser echo from the Earth's surface for interpretation of surface height distribution within the 100 m diam. sensor footprint. For example, tree height can be determined by measuring the characteristic double-pulse signature that results from a separation in time of laser backscatter from tree canopies and the underlying ground. This is accomplished with a pulse waveform digitizer that samples the detector output with an adjustable resolution of 2 nanoseconds or wider intervals in a 100 sample window centered on the return pulse echo. The digitizer makes the SLA into a high resolution surface lidar sensor. It can also be used for cloud and atmospheric aerosol lidar measurements by lengthening the sampling window and degrading the waveform resolution. Detailed test objectives for the STS-72 mission center on the acquisition of sample data sets for land topography and vegetation height, waveform digitizer performance, and verification of data acquisition algorithms. The operational concept of SLA is illustrated in Fig. 1 where a series of 100 m footprints stretch in a profile of Earth surface topography along the nadir track of the Space Shuttle. The location of SLA as a dual canister payload on the Hitchhiker Bridge Assembly in Bay 12 of the Space Shuttle Endeavor can also be noted in this figure. Full interpretation of the SLA range measurement data set requires a 1 m knowledge of the Orbiter trajectory and better than 0.1 deg knowledge of Orbiter pointing angle. These ancillary data sets will be acquired during the STS-72 mission with an on-board Global Positioning System (GPS) receiver, K-band range and range-rate tracking of the Orbiter through TDRSS, and use of on-board inertial measurement units and star trackers. Integration and interpretation of all these different data sets as a pathfinder investigation for accurate determination of Earth surface elevation is the overall science of the SLA investigation.

Bufton, Jack↗

NASA NTRS1972

Study of mathematical modeling of communication systems transponders and receivers

The modeling of communication receivers is described at both the circuit detail level and at the block level. The largest effort was devoted to developing new models at the block modeling level. The available effort did not permit full development of all of the block modeling concepts envisioned, but idealized blocks were developed for signal sources, a variety of filters, limiters, amplifiers, mixers, and demodulators. These blocks were organized into an operational computer simulation of communications receiver circuits identified as the frequency and time circuit analysis technique (FATCAT). The simulation operates in both the time and frequency domains, and permits output plots or listings of either frequency spectra or time waveforms from any model block. Transfer between domains is handled with a fast Fourier transform algorithm.

Walsh, J. R.↗

NASA NTRS1993

Streak camera based SLR receiver for two color atmospheric measurements

To realize accurate two-color differential measurements, an image digitizing system with variable spatial resolution was designed, built, and integrated to a photon-counting picosecond streak camera, yielding a temporal scan resolution better than 300 femtosecond/pixel. The streak camera is configured to operate with 3 spatial channels; two of these support green (532 nm) and uv (355 nm) while the third accommodates reference pulses (764 nm) for real-time calibration. Critical parameters affecting differential timing accuracy such as pulse width and shape, number of received photons, streak camera/imaging system nonlinearities, dynamic range, and noise characteristics were investigated to optimize the system for accurate differential delay measurements. The streak camera output image consists of three image fields, each field is 1024 pixels along the time axis and 16 pixels across the spatial axis. Each of the image fields may be independently positioned across the spatial axis. Two of the image fields are used for the two wavelengths used in the experiment; the third window measures the temporal separation of a pair of diode laser pulses which verify the streak camera sweep speed for each data frame. The sum of the 16 pixel intensities across each of the 1024 temporal positions for the three data windows is used to extract the three waveforms. The waveform data is processed using an iterative three-point running average filter (10 to 30 iterations are used) to remove high-frequency structure. The pulse pair separations are determined using the half-max and centroid type analysis. Rigorous experimental verification has demonstrated that this simplified process provides the best measurement accuracy. To calibrate the receiver system sweep, two laser pulses with precisely known temporal separation are scanned along the full length of the sweep axis. The experimental measurements are then modeled using polynomial regression to obtain a best fit to the data. Data aggregation using normal point approach has provided accurate data fitting techniques and is found to be much more convenient than using the full rate single shot data. The systematic errors from this model have been found to be less than 3 ps for normal points.

Varghese, Thomas K.↗

DOE OSTI2020

Time Resolution Measurements of EJ-232Q With Single- and Dual-Sided Readouts

In view of its importance to applications in high energy physics, nuclear security, and time-of-flight PET, the time resolution of 0.5% benzophenone-quenched EJ-232Q plastic scintillators of various lengths was studied using a fast waveform digitizer and fast PMTs. The scintillator sizes were 5 mm × 5 mm × 10, 25, 50, 75, or 100 mm. First, the timing measurements were performed by using a single PMT reading out one end. The degradation in time resolution with increasing scintillator length was studied. The irradiation geometry was also observed to have a non-negligible impact on the time resolution due to the variation of depth-of-interaction positions, especially for the scintillators with a high aspect ratio. Second, two PMTs were employed for a dual-sided readout. This readout together with the mean time method was demonstrated to be an effective approach to improve the timing performance since the depth-of-interaction-related effect was considerably mitigated. A time resolution of 79 ± 1 ps (full-width at half-maximum) has been achieved for the 5 × 5 × 10 mm 3 scintillator. The 5 × 5 × 100 mm 3 scintillator can still provide a time resolution as good as 172 ± 1 ps. Here, the energy range used for each was between 340 keVee and 450 keVee. In addition, the time information from individual PMTs in a dual-sided readout was able to show a time resolution almost identical to that determined with the mean-time-based technique after correcting for the depth-of-interaction variation for the 100 mm long scintillator.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

NASA NTRS1978

Acoustical effects of blade tip shape changes on a full scale helicopter rotor in a wind tunnel

Four tip shapes were tested. They were rectangular, swept, tapered, and swept-tapered. The measured data covered a wide range of operating conditions. The range of advancing tip Mach numbers were between 0.72 to 0.96, and the advance ratios were from 0.2 to 0.375. At low and moderate advancing tip Mach number, the data in the dbA scale appear to indicate the swept tip is the quietest, swept tapered the second, tapered third and rectangular the most noisy. Above an advancing tip Mach number of about 0.89, a distinct acoustical pulse can be observed, which dominates the acoustical waveform. The pulse shape is symmetric at moderate tip Mach number, changing to a sawtooth shape at high advancing tip Mach numbers. Based on the amplitude of the impulsive noise, it appears the swept-tapered tip is the quietest, tapered tip the second, swept tip third and square tip the most noisy. The data presented in this report should be useful as data bases for modeling and evaluating helicopter impulsive noise.

Lee, A.↗

DOE OSTI2023

A New Plasma Radar Concept for Simultaneous Magnetic and Density Measurements

An innovative, compact 288GHz interferometer has been fabricated, tested, installed and successfully demonstrated on the LAPD-U magnetized plasma at UCLA. The system takes advantage of frequency modulated (FM) radar techniques to deliver a compact heterodyne system. In addition, the reflected power from the source is taken advantage of to eliminate the need for additional quasi-optical components. Electron density in LAPD-U plasma has recently been increased substantially thereby requiring a higher frequency/shorter wavelength interferometer to avoid deleterious refractive effects. This system satisfies those needs. The system uses a 96GHz varactor tuned Gunn oscillator which passes to a passive tripler. This tripler has ~3% conversion efficiency. The 288GHz radiation is then coupled to free space using a so-called dual-mode or Pickett horn. The output 288GHz beam is then coupled to an aspheric lens manufactured from low-loss, high-density polyethylene. This lens is employed to collimate the emerging beam. Small axial adjustment of the lens position can also be used to create a slowly focusing beam so as to optimize the measured signal. In addition, up-down or side-to-side adjustment of the lens can be utilized to steer the beam vertically or horizontally – again to optimize alignment. The propagating beam passes through a beam splitter and then through a water-free, bubble-free fused quartz window into the LAPD-U vacuum vessel. The beam-splitter is a thin sheet of G10 which reflects a small fraction of the incident power (~5 %) towards a zero-bias detector optimized for the frequency range from 220 to 300GHz. Note that waveguides at this frequency have dimensions of ~0.9mm x 0.45mm and so have very large conductive losses. This drives the use of quasi-optical propagation. The detector requires no DC bias and is very responsive (> 1V/mW into 1MΩ). Radiation is coupled to the detector via a similar lens-horn arrangement used for the launch. This reflected beam acts as the local oscillator or reference millimeter-wave beam for the detector. The remainder of the launched source beam then enters the LAPD-U vacuum vessel and passes through the plasma at the mid-plane until reaching the opposing port which is closed off with an aluminum flange. This flange is used as a mirror to retroreflect the incident 288GHz beam back along its path. The retroreflected beam exits the input port but does NOT couple directly into the zero-bias detector. Instead, the majority of the return power continues towards the 288GHz source. As mentioned above the transmitted beam enters the source a second time. This would appear undesirable. However, at these frequencies multipliers are highly non-linear elements which results in a significant portion of the return beam (~20%) re-emerging from the multiplier and horn and then coupling via the G10 beam-splitter to the zero-bias detector. This approach eliminated the need for a second quasi-optical beam-splitter. The system is extremely compact measuring approximately 28 inches x 20 inches. The above did not explain how heterodyne operation was achieved. As mentioned above the Gunn oscillator is able to be varactor tuned. This allows a low voltage to be applied to control the operating frequency of the Gunn oscillator. During heterodyne operation a sawtooth shaped voltage is applied to the varactor at 750kHz using an 80MHz Arbitrary Waveform Generator (AWG). This voltage changes the Gunn frequency linearly during the up-sweep which is then reset abruptly at the sawtooth crash to be immediately followed by another linear sweep. Passage through the 288GHz multiplier triples the frequency change experienced by the electromagnetic wave. These frequency changes are small – tens of megahertz. This FM radar approach results in the launched electromagnetic wave frequencies at the detector for the reference and plasma wave to be different. The approximately 10 ns delay propagation delay for the plasma beam results in the local oscillator and plasma beams NOT having an identical frequency – there is in fact a fixed difference frequency. The frequency tuning level of the Gunn oscillator is then adjusted so that there is ONE cycle of this difference frequency during each linear ramp. During the sawtooth crash or downward re-sweep this one cycle replays in reverse but on a very fast timescale. The process then repeats. Low-pass filtering eliminates the fast re-sweep to leave a pure sine wave heterodyne signal. When the plasma is present it introduces a phase delay in the sine wave (caused by the extremely small Doppler shift resulting from the optical path length change). Of course, to measure this phase change we need a reference. This is simply obtained from an arbitrary waveform generator which provides a synchronized output pulse train which again is low pass filtered to obtain a 750kHz sinusoidal voltage reference for the interferometer. The interferometer was installed on LAPD-U where it has worked reliably and has established that electron densities exceeding 1x10 13 cm -3 are routinely achieved. In addition, the system sensitivity was able to easily observe density fluctuation at frequencies up to 50kHz. FM Radar techniques have enabled a full demonstration of a compact, sensitive, high frequency (288GHz/1mm) heterodyne interferometer.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2021

Seismo-Acoustic Wave Simulation for Earthquake-Generated Infrasound

It is well known that underground explosions or earthquakes can generate infrasound in the atmosphere. Seismoacoustic coupling between the ground motions and atmosphere is a key mechanism for infrasound generation. Epicentral infrasound which is direct response to epicentral ground motions are extensively studied for event detection and discrimination. Diffracted infrasound generated from Rayleigh waves are often observed and reported. Recent studies suggested the surface topography is also important for infrasound generation, and mountains or topographic peaks can act as effective secondary infrasound sources. In this study, we perform full 3-D seismoacoustic simulations to understand a coupled seismic and acoustic wave generation by an earthquake and energy partitioning across the solid Earthocean-atmosphere system. By including the effects of realistic topography, bathymetry, and earthquake focal mechanism, we provide quantitative information about seismoacoustic wavefields and secondary infrasound sources across the solid-water, solid-air, and water-air interfaces. Acoustic wave transmission across water-air interface is also investigated by using synthetic waveforms. Our study suggests that shallow seas near the coastline or over elevated seafloors can transmit significant amount of energy from acoustic waves in water and can be effective secondary infrasound sources.

58 GEOSCIENCES↗

DOE OSTI2024

Collaborative Research: Properties and Dynamics of the Shallow Crust (Final Report)

Ground motions recorded at one location are often extrapolated to nearby regions within a given radius. Here, we explore the appropriateness of spatial extrapolation using data from seven small aperture seismic network deployments in southern California. Six of these deployments are linear arrays of 4-13 stations, and one is a 2D array of 13 stations at Pinyon Flats Observatory. The spatial footprint array diameters are 3 km or less, and each array was operational for a year or more. From our base catalog (M2.5+ earthquakes; 4038 events; September 2010 - June 2023), automated methods remove temporally overprinted waveforms from nearby events (< 5 km) in quick succession (< 5 min) and data with nonviable waveforms. These 200 samples per second data are filtered at 0.5-25 Hz and must have signal-tonoise ratios (SNRs) of 2.5+. Peak ground acceleration (PGA) and peak ground velocity (PGV) are derived individually from the maximum absolute values of each of the 3-component waveforms (vertical, northsouth, and east-west). Five of the seven arrays traverse the San Jacinto fault, and two do not. Ground motion observations are compared with theoretical estimates from Abrahamson et al., 2014. On average, arrays deployed within and across fault zones consistently record ground motions above theoretical expectations, whereas off-fault arrays record ground motions at or slightly below theoretical expectations. We attribute these differences to site conditions because these trends prevail for the full data suites. For each network and each individual channel, the coefficient of variation indicates that the standard deviations are ~30±6% of the mean. Exploring relative ground motion contributions from all three channels (ternary plots), as expected, most data show that vertical ground motions are attenuated compared to horizontal ones. However, this is not always the case for RA array data, where vertical motions can be ~2-3 times larger than horizontal motions for select events near Cahuilla, CA. These anomalously high vertical motions are focal mechanism-related. These results suggest that ternary plots created using only a small amount of data can be used as a data quality metric and a tool to find anomalous features in three-component data.

58 GEOSCIENCES↗

DOE OSTI2021

Seismo-Acoustic Wave Simulation for the 2016 Kumamoto Earthquake

It is well known that underground explosions or earthquakes can generate infrasound in the atmosphere. Seismoacoustic coupling between the ground motions and atmosphere is a key mechanism for infrasound generation. Epicentral infrasound which is direct response to epicentral ground motions are extensively studied for event detection and discrimination. Diffracted infrasound generated from Rayleigh waves are often observed and reported. Recent studies suggested the surface topography is also important for infrasound generation, and mountains or topographic peaks can act as effective secondary infrasound sources. In this study, we perform full 3-D seismoacoustic simulations to understand a coupled seismic and acoustic wave generation by an earthquake and energy partitioning across the solid Earth-ocean-atmosphere system. By including the effects of realistic topography, bathymetry, and earthquake focal mechanism, we provide quantitative information about seismoacoustic wavefields and secondary infrasound sources across the solid-water, solid-air, and water-air interfaces. Acoustic wave transmission across water-air interface is also investigated by using synthetic waveforms. Our study suggests that shallow seas near the coastline or over elevated seafloors can transmit significant amount of energy from acoustic waves in water and can be effective secondary infrasound sources.

58 GEOSCIENCES↗