Search NASA⌕ Search

SEARCH · Search NASA

Results for “seismic reflection”

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 55 records · Page 3
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Neural net detection of seismic features related to gas hydrates and free gas accumulations on the northern U.S. Atlantic margin

Bottom-simulating reflections (BSRs) that sometimes mark the base of the gas hydrate stability zone in marine sediments are often identified based on the reverse polarity reflections that cut across stratigraphic layering in seismic amplitude data. On the northern U.S. Atlantic margin (USAM) between Cape Hatteras and Hudson Canyon, legacy seismic data have revealed pronounced BSRs south of the deepwater extension of Hudson Canyon and more subtle ones from offshore Delaware south to Cape Hatteras, where the reflections sometimes follow stratigraphic layering. Using high-resolution seismic data acquired during the 2018 Mid-Atlantic Resource Imaging Experiment and a supervised neural net, we identify seismic features associated with gas hydrates and/or the top of gas between Hudson Canyon and Cape Hatteras. Using seismic attributes especially sensitive to the presence of gas, we train a neural network algorithm on seismic data from an area with strong BSRs and then apply the model to the rest of the data set. The results indicate that gas hydrate and/or shallow free gas are significantly more widespread on the northern part of the USAM than previously known. Seismic indicators of gas extend landward from the 2000 m isobath to the upper continental slope in sectors with (offshore Virginia) and, to a lesser extent, without (offshore New Jersey) pervasive upper slope methane seeps. Higher sand content and intermediate sediment thickness, factors related to the container size and gas charge in a petroleum systems framework, are associated with more robust gas indicators.

Geochemistry & Geophysics↗
DOE OSTI2021

Determining fast-S and slow-S propagation directions with SV-P data produced by buried explosives and recorded with vertical geophones

We have evaluated the concept of practicing S-wave reflection seismology with legacy 3D seismic data generated by a P-wave source and recorded with only vertical geophones. This type of S-wave imaging is based on the principle that seismic P-wave sources not only produce a downgoing illuminating P wavefield, but they also simultaneously produce a downgoing illuminating SV wavefield that, in almost all cases, is suitable for S-wave reflection imaging. The S-mode used in this study is the SV-P, or converted-P, mode. This mode involves a downgoing illuminating SV wavefield and an upgoing reflected P-mode that is recorded by vertical geophones. In flat-layered stratigraphy, the lengths of the SV and P raypaths in SV-P imaging are identical to the lengths of the SV and P raypaths in P-SV imaging with P-sources and 3C geophones. P-SV imaging of deep rocks has been practiced for more than two decades; SV-P imaging is a new concept. SV-P data should provide the same options for investigating deep rocks as do P-SV data. We have determined one of the equivalences between SV-P data extracted from vertical-geophone data and P-SV data extracted from horizontal geophones: that both modes react to azimuth-dependent variations in the S velocity in anisotropic rocks. Azimuthal variations in the SV-P traveltime can be used to define the polarization direction of the fast-S-wave mode, which is also the azimuth of the maximum horizontal stress (SHmax). Our investigation demonstrates a noninvasive method for monitoring changes in the SHmax azimuth across a CO 2 storage reservoir, or any targeted porous rock, as fluids are cycled into, and then out of, that rock’s pore space.

Geochemistry & Geophysics↗
DOE OSTI2021

Joint impedance and facies inversion of time-lapse seismic data for improving monitoring of CO 2 incidentally stored from CO 2 EOR

Time-lapse seismic monitoring is an effective and proven technology for mapping the distribution of CO 2 in a subsurface reservoir. When injected CO 2 displaces other reservoir fluids, porous-medium properties are changed and thus the seismic impedance changes, causing time-lapse seismic amplitude differences in the injection zones. The analysis and interpretation of images created from these amplitude differences can provide information about reservoir architecture and the CO 2 migration within the reservoir. Incorporating seismic inversion and rock physics into the interpretation of time-lapse seismic data can considerably improve the modeling and monitoring to detect and assess the location of CO 2 over time. The joint inversion method presented in this paper has an integral representation of the geology in the inversion algorithm using elastic facies, which provides information about the spatial distribution of the geologic heterogeneities controlling the movement of fluids in the reservoir. The method was successfully applied to time-lapse seismic data from a mature oil field undergoing CO 2 enhanced oil recovery. The estimated seismic acoustic impedances and facies reflect the characteristics of individual geologic facies and fluid conditions of the reservoir subject to CO 2 injection. Finally, the probabilities estimated by the joint impedance and facies inversion for the reservoir's litho-fluid facies can be used for forecasting CO 2 saturation and pressure changes within the target reservoir.

4D seismic↗
DOE OSTI2022

Mineralogy, fabric and deformation domains in D" across the southwestern border of the African LLSVP

Recent advances in seismic anisotropy studies that jointly use reflections and shear wave splitting have proven to place tight constraints on the plausible anisotropic and deformation scenarios in the D" region. We apply this novel methodology to a large area of the D" region beneath the South Atlantic, in proximity to and within the African large low seismic velocity province (LLSVP). This area of the mantle is characterized by a transition from fast to slow seismic velocity anomalies and it is thought to be the location of deep-seated plumes responsible for hotspot volcanism. Attempting to probe mantle composition and deformation along the LLSVP borders may provide key information on mantle dynamics. By analysing seismic phases sampling this region, we detect a D" discontinuity over a large area beneath the South Atlantic, with inferred depth ranges ~170 to ~240 km above the core–mantle boundary. We find evidence for a D" reflector within the area of the LLSVP. Shear wave splitting observations suggest that anisotropy is present in this region of the mantle, in agreement with previous studies that partially sampled this region. We model the observations considering lattice- and shape-preferred orientation of materials expected in the D" region. A regional variation of mineralogy, phase transition boundaries, and deformation direction is required to explain the data. We infer two distinct domains of mineralogy and deformation: aligned post-perovskite outside the LLSVP and aligned bridgmanite within the LLSVP. Finally, the scenario depicted by this study agrees well with the current hypotheses for the composition of the LLSVP and with the prevalence of vertical deformation directions expected to occur along the LLSVPs borders.

58 GEOSCIENCES↗
DOE OSTI2025

Exploring Data Set Bias and Decision Support with Predictive Uncertainty Through Bayesian Approximations and Convolutional Neural Networks

Individual seismic catalogs can contain multiscale observations from fault level to global scales and associated waveforms from discrete events reflect crustal structure across many different scales and locations. Seismic network aperture, geographic location, and observation distance may not provide informative guidance or intuition on how different catalogs will behave across models trained under different conditions. We rely on uncertainty to provide guardrails for when to trust model decisions, but understanding when our uncertainty is trustworthy is an open challenge. Here, in this work, we explore Bayesian approximation methods for assigning predictive uncertainty in seismic event classification problems. We find that computationally expensive Bayesian approximations do not outperform simple ensemble methods. We also find that when exploiting multiple seismic event catalogs, joint training with data from all the catalogs combined with Bayesian approximations and supervised training for classification can obscure bias and result in less robust uncertainty while also not providing substantial performance benefits compared to training individual models for each catalog.

58 GEOSCIENCES↗
DOE OSTI2023

Dating submarine landslides using the transient response of gas hydrate stability

Abstract Submarine landslides are prevalent on the modern-day seafloor, yet an elusive problem is constraining the timing of past slope failure. We present a novel age-dating technique based on perturbations to underlying gas hydrate stability caused by slide-impacted seafloor changes. Using three-dimensional (3-D) seismic data, we mapped an irregular bottom simulating reflection (BSR) underneath a submarine landslide in the Orca Basin, Gulf of Mexico. The irregular BSR mimics the pre-slide seafloor geometry rather than the modern bathymetry. Therefore, we suggest that the gas hydrate stability zone (GHSZ) is still adjusting to the post-slide sediment temperature. We applied transient conductive heat-flow modeling to constrain the response of the GHSZ to the slope failure, which yielded a most likely age of ca. 8 ka, demonstrating that gas hydrate can respond to landslides even on multimillennial time scales. We further provide a generalized analytical solution that can be used to remotely date submarine slides in the absence of traditional dating techniques.

Geology↗
DOE OSTI2020

Stochastic Inversion of Gaussian Random Media Using Transverse Coherence Functions for Reflected Waves: Theory and Method

The transverse coherence functions (TCFs) of phase and amplitude fluctuations of a seismic wave are powerful to estimate the spatial distribution, length scales, and strength of random heterogeneities. However, TCFs have been formulated for transmitted waves only, not for reflected waves. In this paper, we derive reflection TCFs for Gaussian random media. Furthermore, we propose to invert for Gaussian random media using the reflection TCFs based on the grid search. We validate the new reflection TCF formulas using 2D finite-difference numerical experiments. The numerical example also illustrates the feasibility and efficiency of the inversion. The stochastic inversion using reflected waves can be used in both exploration and global seismology.

58 GEOSCIENCES↗
DOE OSTI2021

Listening for the Landing: Seismic Detections of Perseverance's Arrival at Mars With InSight

The entry, descent, and landing (EDL) sequence of NASA's Mars 2020 Perseverance Rover will act as a seismic source of known temporal and spatial localization. We evaluate whether the signals produced by this event will be detectable by the InSight lander (3,452 km away), comparing expected signal amplitudes to noise levels at the instrument. Modeling is undertaken to predict the propagation of the acoustic signal (purely in the atmosphere), the seismoacoustic signal (atmosphere-to-ground coupled), and the elastodynamic seismic signal (in the ground only). Our results suggest that the acoustic and seismoacoustic signals, produced by the atmospheric shock wave from the EDL, are unlikely to be detectable due to the pattern of winds in the martian atmosphere and the weak air-to-ground coupling, respectively. However, the elastodynamic seismic signal produced by the impact of the spacecraft's cruise balance masses on the surface may be detected by InSight. The upper and lower bounds on predicted ground velocity at InSight are 2.0 × 10 -14 and 1.3 × 10 -10 m s -1 . The upper value is above the noise floor at the time of landing 40% of the time on average. The large range of possible values reflects uncertainties in the current understanding of impact-generated seismic waves and their subsequent propagation and attenuation through Mars. Uncertainty in the detectability also stems from the indeterminate instrument noise level at the time of this future event. A positive detection would be of enormous value in constraining the seismic properties of Mars, and in improving our understanding of impact-generated seismic waves.

58 GEOSCIENCES↗
DOE OSTI2020

The Role of Premagmatic Rifting in Shaping a Volcanic Continental Margin: An Example From the Eastern North American Margin

Both magmatic and tectonic processes contribute to the formation of volcanic continental margins. Such margins are thought to undergo extension across a narrow zone of lithospheric thinning (~100 km). New observations based on existing and reprocessed data from the Eastern North American Margin contradict this hypothesis. With ~64,000 km of 2-D seismic data tied to 40 wells combined with published refraction, deep reflection, receiver function, and onshore drilling efforts, we quantified along-strike variations in the distribution of rift structures, magmatism, crustal thickness, and early post-rift sedimentation under the shelf of Baltimore Canyon Trough (BCT), Long Island Platform, and Georges Bank Basin (GBB). Results indicate that BCT is narrow (80–120 km) with a sharp basement hinge and few rift basins. The seaward dipping reflectors (SDR) there extend ~50 km seaward of the hinge line. In contrast, the GBB is wide (~200 km), has many syn-rift structures, and the SDR there extend ~200 km seaward of the hinge line. Early post-rift depocenters at the GBB coincide with thinner crust suggesting “uniform” thinning of the entire lithosphere. Models for the formation of volcanic margins do not explain the wide structure of the GBB. We argue that crustal thinning of the BCT was closely associated with late syn-rift magmatism, whereas the broad thinning of the GBB segment predated magmatism. Correlation of these variations to crustal terranes of different compositions suggests that the inherited rheology determined the premagmatic response of the lithosphere to extension.

58 GEOSCIENCES↗
DOE OSTI2022

Estimating SHmax azimuth with P sources and vertical geophones: Use P-P reflection amplitudes or use SV-P reflection times?

We compared two methods for extracting the azimuth of maximum horizontal stress (SHmax) from 3D land-based seismic data generated by a P source and recorded with vertical geophones. In the first method, we used the direct-SV mode that is produced by all land-based P sources. P sources generate SV illumination that radiates in all azimuth directions from a source station and creates SV-P reflections that are recorded by vertical geophones. Unless stratigraphy has steep dip, SV-P raypaths recorded by vertical geophones are the reverse of P-SV raypaths recorded by horizontal geophones. Thus, SV-P data provide the same S-wave sensitivity to stress fields as popular P-SV data do. In the second method, we retrieved P-P reflections and then performed an amplitude-variation-with-azimuth (AVA) analysis of the amplitude-gradient behavior of P-P reflection wavelets. We did this analysis in narrow azimuth corridors to determine the gradient of reflection-wavelet amplitudes as a function of azimuth. This P-P AVA amplitude-gradient method has been of great interest in the reflection seismology community since it was introduced in the late 1990s. Each of these methods, AVA analysis of the gradient of P-P reflection amplitudes and azimuth-dependent arrival times of SV-P reflections, can be used to determine the azimuth of SHmax stress. We compare the results of the two methods with ground truth measurements of SHmax azimuth at a CO 2 sequestration site in the Michigan Basin. SHmax azimuths were determined from P-P and SV-P data at three major boundaries at depths of approximately 3500 ft (1067 m), 5500 ft (1676 m), and 7500 ft (2286 m). Two estimates of SHmax azimuth (one using SV-P data and one using P-P data) were made at each stacking bin inside a 24 mi 2 (62 km 2 ) image space. The result was approximately 98,000 estimates of SHmax azimuth across each of these three boundaries for each of these two prediction strategies. Histogram displays of PP AVA gradient estimates had peaks at correct azimuths of SHmax at all three depths, but the spread of the distributions widened with depth and split into two peaks at the deepest boundary. In contrast, each histogram of SHmax azimuth predicted by azimuth-dependent SV-P traveltimes had a single, definitive peak that was positioned at the correct SHmax azimuth at all three boundary depths.

Geochemistry & Geophysics↗
DOE OSTI2025

Evaluation of Station Performance of the Idaho National Laboratory Seismic Monitoring Network Using Network Detection Thresholds

The Idaho National Laboratory (INL) Seismic Monitoring Network is located in eastern Idaho and monitors a portion of the intermountain seismic belt. It has been in place for 50 yr and has undergone several major changes, the most recent of which has been the transition to the Antelope real‐time acquisition system and the implementation of automatic phase picking algorithms to aid in analysis. This study discusses the efforts to evaluate the performance of the INL seismic monitoring network (and other surrounding stations) using the new real‐time acquisition system. The method outlined by Wilson et al. (2021) is used to develop an empirical relationship between the observability of local earthquakes as a function of magnitude and distance. This relationship is used to produce detection thresholds for Pwaves for all stations of interest. The INL seismic network has two main goals: monitor tectonic‐and volcanic‐related events and measure ground motions for input into seismic hazard analysis. Because of these two overall objectives, several seismic stations have been installed near critical facilities and, therefore, are not as quiet as stations that are used primarily for earthquake detection. This is reflected in their detection thresholds, which are much smaller for stations away from facilities. This study shows that the INL Seismic Monitoring Network is able to detect earthquakes near INL facilities with M L > 1.2, with redundancies built in to ensure this sensitivity even if data became unavailable from some stations. This study also shows “holes” in the monitoring network where the detection of smaller earthquakes is highly dependent on sparsely placed seismic stations. In conclusion, the results of this study will be used to govern plans for expansion of earthquake monitoring in Idaho and the surrounding region and to fine‐tune the detection thresholds for individual stations.

58 - GEOSCIENCES↗
DOE OSTI2024

Revealing complex subsurface dynamics with continuous seismic monitoring: Observations using distributed acoustic sensing and surface orbital vibrators during hydraulic fracturing

Understanding hydraulic fracturing is crucial to improving the stimulation of unconventional reservoirs and increasing fluid production. This study develops a novel seismic monitoring technology using distributed acoustic sensing (DAS) and surface orbital vibrators (SOV) to capture fracture seismic response and mechanical properties at high temporal intervals. We analyze continuous time-lapse vertical seismic profiling (VSP) data acquired every hour during the first nine days of treatment of an unconventional reservoir in the Austin Chalk/Eagle Field Laboratory. The VSP data contain clear seismic signals scattered from the activated fractures. The spatiotemporal changes of the fracture reflectivity revealed by the SOV/DAS data correlate well with the observations of fracture locations inferred from low-frequency DAS data. These results capture the fracture opening and closure processes, as well as highlight potential prestage activations of the fractures due to hydraulic connectivity with preexisting fracture systems. Therefore, analysis of the presented data set provides a unique opportunity to understand fracture initiation and subsequent evolution, not only in the context of unconventional resources but also in enhanced geothermal systems.

Correa, Julia↗
DOE OSTI2022

Seismic monitoring of well integrity

Seismic tube waves, produced by flexure of the well boundary, pressure changes in the fluid in the well, and deformation of the material immediately surrounding the well, are particularly sensitive to variations in the state of the well. In this work, we evaluate a direct approach for generating and observing tube waves as a means of detecting well damage. While we find that it can be difficult to reliably excite observable tube waves without a very strong surface source, time-frequency techniques can be employed to increase the detectability of tube wave reflections. New technologies, particularly distributed acoustic sensing, hold great promise for evaluating well integrity by monitoring tube waves, temperature changes, and seismic noise due to well deformation and fluid leakage.

58 GEOSCIENCES↗
DOE OSTI2020

The geometry of signal space: a case study of direct mapping between seismic signals and event distribution

SUMMARY Under favourable circumstances, seismic waveforms corresponding to an ensemble of events related by a common, spatially distributed process collectively exhibit a regular, signal-space geometry. When events in the ensemble have a common, or nearly common, source mechanism, this geometry is a distorted image of the distribution of events in the source region. The signal-space image can be visualized using a relatively simple waveform alignment and projection operation. Ensemble waveform correlation measurements can be inverted to estimate the distribution of the events in the source region, up to an arbitrary rotation, reflection and scaling, with residual distortion. We demonstrate these concepts with synthetic waveforms and with observations of long-wall mining induced seismicity for which substantial ground truth information is available. Our experience with these data has implications for location, correlation detection and machine learning and possible application to studies of repeating events in induced, volcanic and glacial seismicity. Our results place limits on the widely held assumption that waveform correlation is a useful measure of event separation. We suggest that the constraints on event separation need to be evaluated in the context of a population of related events, whose waveforms sample the signal space image of the source region. A better indicator of event separation is the length of the shortest path in signal space along the image.

58 GEOSCIENCES↗
DOE OSTI2021

LL18-Signal Analysis-NDD2Ad (Final Report)

Under favorable circumstances, seismic waveforms corresponding to an ensemble of events related by a common, spatially-distributed process collectively exhibit a regular, signal-space geometry. When events in the ensemble have a common, or nearly common, source mechanism, this geometry is a distorted image of the distribution of events in the source region. The signal space image can be visualized using a relatively simple waveform alignment and projection operation. Ensemble waveform correlation measurements can be inverted to estimate the distribution of the events in the source region, up to an arbitrary rotation, reflection, and scaling with residual distortion. We demonstrate these concepts with synthetic waveforms and with observations of long-wall mining induced seismicity for which substantial ground truth information is available. Our experience with these data has implications for location, correlation detection and machine learning, and possible application to studies of repeating events in induced, volcanic, and glacial seismicity. Our results place limits on the widely-held assumption that waveform correlation is a useful measure of event separation. We suggest that the constraints on event separation need to be evaluated in the context of a population of related events, whose waveforms sample the signal space image of the source region. A better indicator of event separation is the length of the shortest path in signal space along the image.

58 GEOSCIENCES↗
DOE OSTI2024

Heterogeneous mantle effects on the behavior of SmKS waves and outermost core imaging

Seismic traveltime anomalies of waves that traverse the uppermost 100–200 km of the outer core have been interpreted as evidence of reduced seismic velocities (relative to radial reference models) just below the core–mantle boundary (CMB). These studies typically investigate differential traveltimes of SmKS waves, which propagate as P waves through the shallowest outer core and reflect from the underside of the CMB m times. The use of SmKS and S(m-1)KS differential traveltimes for core imaging are often assumed to suppress contributions from earthquake location errors and unknown and unmodelled seismic velocity heterogeneity in the mantle. The goal of this study is to understand the extent to which differential SmKS traveltimes are, in fact, affected by anomalous mantle structure, potentially including both velocity heterogeneity and anisotropy. Velocity variations affect not only a wave's traveltime, but also the path of a wave, which can be observed in deviations of the wave's incoming direction. Since radial velocity variations in the outer core will only minimally affect the wave path, in contrast to other potential effects, measuring the incoming direction of SmKS waves provides an additional diagnostic as to the origin of traveltime anomalies. Here we use arrays of seismometers to measure traveltime and direction anomalies of SmKS waves that sample the uppermost outer core. We form subarrays of EarthScope's regional Transportable Array stations, thus measuring local variations in traveltime and direction. We observe systematic lateral variations in both traveltime and incoming wave direction, which cannot be explained by changes to the radial seismic velocity profile of the outer core. Moreover, we find a correlation between incoming wave direction and traveltime anomaly, suggesting that observed traveltime anomalies may be caused, at least in part, by changes to the wave path and not solely by perturbations in outer core velocity. Modelling of 1-D ray and 3-D wave propagation in global 3-D tomographic models of mantle velocity anomalies match the trend of the observed traveltime anomalies. Overall, we demonstrate that observed SmKS traveltime anomalies may have a significant contribution from 3-D mantle structure, and not solely from outer core structure.

58 GEOSCIENCES↗
DOE OSTI2024

Heterogeneous mantle effects on the behavior of SmKS waves and outermost core imaging

Seismic traveltime anomalies of waves that traverse the uppermost 100–200 km of the outer core have been interpreted as evidence of reduced seismic velocities (relative to radial reference models) just below the core–mantle boundary (CMB). These studies typically investigate differential traveltimes of SmKS waves, which propagate as P waves through the shallowest outer core and reflect from the underside of the CMB m times. The use of SmKS and S(m- 1)KS differential traveltimes for core imaging are often assumed to suppress contributions from earthquake location errors and unknown and unmodelled seismic velocity heterogeneity in the mantle. The goal of this study is to understand the extent to which differential SmKS traveltimes are, in fact, affected by anomalous mantle structure, potentially including both velocity heterogeneity and anisotropy. Velocity variations affect not only a wave’s traveltime, but also the path of a wave, which can be observed in deviations of the wave’s incoming direction. Since radial velocity variations in the outer core will only minimally affect the wave path, in contrast to other potential effects, measuring the incoming direction of SmKS waves provides an additional diagnostic as to the origin of traveltime anomalies. Here we use arrays of seismometers to measure traveltime and direction anomalies of SmKS waves that sample the uppermost outer core. We form subarrays of EarthScope’s regional Transportable Array stations, thus measuring local variations in traveltime and direction. We observe systematic lateral variations in both traveltime and incoming wave direction, which cannot be explained by changes to the radial seismic velocity profile of the outer core. Moreover, we find a correlation between incoming wave direction and traveltime anomaly, suggesting that observed traveltime anomalies may be caused, at least in part, by changes to the wave path and not solely by perturbations in outer core velocity. Modelling of 1-D ray and 3-D wave propagation in global 3-D tomographic models of mantle velocity anomalies match the trend of the observed traveltime anomalies. Overall, we demonstrate that observed SmKS traveltime anomalies may have a significant contribution from 3-D mantle structure, and not solely from outer core structure.

58 GEOSCIENCES↗
DOE OSTI2022

Observations of Mantle Seismic Anisotropy Using Array Techniques: Shear-Wave Splitting of Beamformed SmKS Phases

Shear-wave splitting measurements are commonly used to resolve seismic anisotropy in both the upper and lowermost mantle. Typically, such techniques are applied to SmKS phases that have reflected (m-1) times off the underside of the core-mantle boundary before being recorded. Practical constraints for shear-wave splitting studies include the limited number of suitable phases as well as the large fraction of available data discarded because of poor signal-to-noise ratios (SNRs) or large measurement uncertainties. Array techniques such as beamforming are commonly used in observational seismology to enhance SNRs, but have not been applied before to improve SmKS signal strength and coherency for shear wave splitting studies. Here, we investigate how a beamforming methodology, based on slowness and backazimuth vespagrams to determine the most coherent incoming wave direction, can improve shear-wave splitting measurement confidence intervals. Through the analysis of real and synthetic seismograms, we show that (a) the splitting measurements obtained from the beamformed seismograms (beams) reflect an average of the single-station splitting parameters that contribute to the beam; (b) the beams have (on average) more than twice as large SNRs than the single-station seismograms that contribute to the beam; (c) the increased SNRs allow the reliable measurement of shear wave splitting parameters from beams down to average single-station SNRs of 1.3. Beamforming may thus be helpful to more reliably measure splitting due to upper mantle anisotropy. Moreover, here we show that beamforming holds potential to greatly improve detection of lowermost mantle anisotropy by demonstrating differential SKS–SKKS splitting analysis using beamformed USArray data.

54 ENVIRONMENTAL SCIENCES↗