Search NASA⌕ Search

SEARCH · Search NASA

Results for “seismic data processing”

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 415 records · Page 23

Seismoelectric Effects for Geothermal Resources Assessment and Monitoring (SEE4GEO)

The seismoelectric effects technique (SEE) is a new and innovative approach for geothermal subsurface imaging and monitoring at reservoir scale. The objective of this project is to assess SEE in terms of data acquisition, cost and quality, and to determine its capability in comparison with classical imaging and monitoring techniques, particularly decoupled seismic and electromagnetic methods. This will be achieved by (1) development of a fast, true 3D numerical package, handling SEE imaging and subsurface properties characterization, including resistivity and permeability, (2) laboratory experiments performed in a controlled environment to define optimal deployment design, data quality, and inform field deployment, and (3) field surveys to ultimately test and draw lessons for practical use of SEE technology. There is a relatively extensive body of work in the literature on SEE, and members of this consortium have been involved in theoretical and numerical development of SEE modeling as well as laboratory experiments. Nevertheless, to our knowledge few, if any, documented efforts have been specifically targeting the use of SEE for geothermal subsurface imaging and monitoring. The strength and originality of our proposal rely on an integrated approach leveraging numerical, laboratory and field experiments, to properly document the practical use of SEE. Through this process, SEE in-hand technology for the geothermal industry will be able to progress from a TRL 1 to TRL 3.

58 GEOSCIENCES↗

Data for Monitoring Fracture Saturation with Internal Transportable Seismic Sources and Twin Neural Networks

Seismic coda-wave analysis is a well-developed method for detecting subtle physical changes in complex media by measuring arrival times in the late-arriving energy from multiply-scattered or reflected waves. However, a challenge arises when multiply-scattered waves are not sufficiently separated in time from the direct arrivals to provide a clear coda wave train. Additional complications for monitoring changes in fracture systems arise when the signals originate from unsynchronized internal sources, such as natural induced seismicity, from acoustic emission, or from transportable intra-fracture sources (chattering dust), that generate uncontrolled signals that vary in arrival time, amplitude and frequency content. Here, we use a twin neural network (also known as a Siamese neural network) for dimensionality reduction to analyze signals from chattering dust to classify the fluid saturation state of a synthetic fracture system. The twin neural network with shared weights generates a low-dimensional representation of the data input by minimizing contrastive loss, serving as the input to a multiclass classifier that accurately classifies whether multiple fractures in a fracture system are fully saturated or partially saturated, or whether a change in saturation has occurred in different fractures in the system. These results show that information buried in unresolved codas from uncontrolled sources can be extracted using machine learning to monitor the evolution of fracture systems caused by physical and chemical processes even when the scattered and direct wave fields overlap.

58 GEOSCIENCES↗

The evolution of seafloor venting from hydrate-sealed gas reservoirs

Here we use 3D seismic data to show that three rows of seafloor gas mounds can be traced downward to leak points that lie at the hydrate-gas contact within three individual dipping coarse-grained sand bodies in the Terrebonne Basin, Walker Ridge block 313, northern Gulf of Mexico. We predict the overpressure within the sand bodies by assuming that the gas pressure at the vent points equal the least principal stress. We interpret that free gas accumulates at the base of the hydrate stability zone, trapped by the overlying sand which has a high methane hydrate saturation. The free gas accumulates until the gas pressure at the base of the hydrate stability zone reaches the least principal stress in the overlying mudrock, whereupon hydraulic fractures form and fluids are vented to the surface. The warm rising fluids and perhaps localized exothermic formation of hydrate raise the local salinity and temperature. This process progressively shifts the base of the hydrate stability zone to shallower depths and dissociates the hydrate seal within the sand, which creates new leak points and results in the observed migration of the seafloor vents. Within the southwest corner of the Terrebonne Basin, this process has repeated multiple times within the Blue, Orange, and Green sands. Here this study shows how the hydrate stability zone can influence the location of fluid expulsion and in turn be affected by the warm, rising, saline fluids sealed by hydrate.

58 GEOSCIENCES↗

Structural algorithm to reservoir reconstruction using passive seismic data (synthetic example)

Using of passive seismic observations to detect a reservoir is a new direction of prospecting and exploration of hydrocarbons. In order to identify thin reservoir model we applied the modification of Gaussian elimination method in conditions of incomplete synthetic data. Because of the singularity of a matrix conventional method does not work. Therefore structural algorithm has been developed by analyzing the given model as a complex model. Numerical results demonstrate of its advantage compared with usual way of solution. We conclude that the gas reservoir is reconstructed by retrieving of the image of encasing shale beneath it.

58 GEOSCIENCES↗

Fluoroethylenepropylene ferroelectret films with cross-tunnel structure for piezoelectric transducers and micro energy harvesters

Layered fluoroethylenepropylene (FEP) ferroelectret films with cross-tunnel structure were fabricated from sheets of FEP films by template-patterning followed by a fusion-bonding process and contact charging. The typical piezoelectric d{sub 33} coefficients, measured by a quasi-static method of samples not annealed, are in the range of 1000–3700 pC/N. The resonance behavior of the samples is analyzed by dielectric spectroscopy which also yields Young's modulus. Microphones built with such films exhibit a somewhat decreasing frequency response up to 1 kHz, an increase of the responses due to diffraction effects at higher frequencies, and eventually a peak probably due to a thickness resonance at about 40 kHz. Annealing at 125 °C indicates that the sensitivity stabilizes at about 40% of the original value. From this data, stable dynamic d{sub 33} coefficients of up to 300 pC/N can be calculated. A micro energy harvesting generator utilizing these films and based on the excitation of thickness vibrations is also described. With an active area of 4.3 cm{sup 2} and a seismic mass of 69.5 g, power up to 0.5 μW referred to an acceleration of 1 g can be generated at a frequency of 120 Hz.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Robust inverse scattering full waveform seismic tomography for imaging complex structure

Seismic tomography becomes important tool recently for imaging complex subsurface. It is well known that imaging complex rich fault zone is difficult. In this paper, The application of time domain inverse scattering wave tomography to image the complex fault zone would be shown on this paper, especially an efficient time domain inverse scattering tomography and their run in cluster parallel computer which has been developed. This algorithm is purely based on scattering theory through solving Lippmann Schwienger integral by using Born's approximation. In this paper, it is shown the robustness of this algorithm especially in avoiding the inversion trapped in local minimum to reach global minimum. A large data are solved by windowing and blocking technique of memory as well as computation. Parameter of windowing computation is based on shot gather's aperture. This windowing technique reduces memory as well as computation significantly. This parallel algorithm is done by means cluster system of 120 processors from 20 nodes of AMD Phenom II. Benchmarking of this algorithm is done by means Marmoussi model which can be representative of complex rich fault area. It is shown that the proposed method can image clearly the rich fault and complex zone in Marmoussi model even though the initial model is quite far from the true model. Therefore, this method can be as one of solution to image the very complex mode.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Hydraulic fracture mapping using wavelet-based fusion of wave transmission and emission measurements

Fracture characterization is essential to hydrocarbon and geothermal exploration and production as well as carbon geo-sequestration. Active and passive measurements have been extensively applied to map fractures across various scales. Active and passive methods have inherent advantages and limitations that complement each other. Here, in this study, to improve the mapping of embedded fractures and the surrounding geomechanically altered regions, we integrate active shear-wave transmission measurement with passive acoustic-emission measurement collected during a lab-scale hydraulic-fracturing experiment. The proposed data fusion leverages the advantages of the two modalities of measurements, while minimizing their individual limitations. Polarized shear-wave transmission (active) measurement was collected before and after the hydraulic fracturing, and the acoustic emission (passive) measurement was collected during hydraulic fracturing. Two sets of two-dimensional maps of fracture and fracture-induced damage in axial, median and frontal planes were obtained by separately processing active and passive measurements. These 2D maps are then fed to wavelet-based image fusion technique to integrate the two sources of information for the reliable mapping/imaging of the embedded fractures and the surrounding geomechanically altered regions.

58 GEOSCIENCES↗

Preliminary result of P-wave speed tomography beneath North Sumatera region

The structure of P-wave speed beneath the North Sumatra region was determined using P-wave arrival times compiled by MCGA from time periods of January 2009 to December 2012 combining with PASSCAL data for February to May 1995. In total, there are 2,246 local earthquake events with 10,666 P-wave phases from 63 stations seismic around the study area. Ray tracing to estimate travel time from source to receiver in this study by applying pseudo-bending method while the damped LSQR method was used for the tomographic inversion. Based on assessment of ray coverage, earthquakes and stations distribution, horizontal grid nodes was set up of 30×30 km2 for inside the study area and 80×80 km2 for outside the study area. The tomographic inversion results show low Vp anomaly beneath Toba caldera complex region and around the Sumatra Fault Zones (SFZ). These features are consistent with previous study. The low Vp anomaly beneath Toba caldera complex are observed around Mt. Pusuk Bukit at depths of 5 km down to 100 km. The interpretation is these anomalies may be associated with ascending hot materials from subduction processes at depths of 80 km down to 100 km. The obtained Vp structure from local tomography will give valuable information to enhance understanding of tectonic and volcanic in this study area.

58 GEOSCIENCES↗

EGS Collab Project: Status, Tests, and Data

The EGS (Enhanced Geothermal Systems) Collab project is performing stimulation and flow experiments in highly-monitored and well-characterized intermediate-scale (approximately10 to 20 meter) field test beds at a depth of approximately 1,500 meters in the Sanford Underground Research Facility (SURF) in the Black Hills of South Dakota. Our fracture stimulation and interwell flow tests are performed to better understand processes that control formation of effective subsurface heat exchangers that are critical to the development and success of EGS. Different EGS Collab stimulations will be performed under dissimilar stress conditions produce data for model comparisons that better differentiate stimulation mechanisms and the evolution of permeability enhancement in crystalline rock. EGS Collab experiments provide a means of testing tools, concepts and strategies that could later be employed under geothermal reservoir conditions at DOE’s Frontier Observatory for Research in Geothermal Energy (FORGE) and other enhanced geothermal systems. Key to the project is using numerical simulations in the experiment design and interpretation of results, and the collection of high-quality stimulation/flow, thermal, and geologic/geophysical data. A wide range of geologic and geophysical data have been acquired and interpreted for pre-stimulation testbed characterization. Multiple stimulation and flow tests have been performed at SURF, and multiple data sets have been collected and analyzed during these tests including electrical resistance tomography, microseismic, continuous active source seismic monitoring, distributed temperature, and the SIMFIP tool has been used to identify fracture formation and characteristics. An overview of some of the numerical simulations and experiments performed at SURF are presented.

15 GEOTHERMAL ENERGY↗

Preliminary results of local earthquake tomography around Bali, Lombok, and Sumbawa regions

Bali, Sumbawa, and Lombok regions are located in active tectonic influence by Indo-Australia plate subducts beneath Sunda plate in southern part and local back-arc thrust in northern part the region. Some active volcanoes also lie from eastern part of Java, Bali, Lombok and Sumbawa regions. Previous studies have conducted subsurface seismic velocity imaging using regional and global earthquake data around the region. In this study, we used P-arrival time from local earthquake networks compiled by MCGA, Indonesia within time periods of 2009 up to 2013 to determine seismic velocity structure and simultaneously hypocenter adjustment by applying seismic tomography inversion method. For the tomographic inversion procedure, we started from 1-D initial velocity structure. We evaluated the resolution of tomography inversion results through checkerboard test and calculating derivative weigh sum. The preliminary results of tomography inversion show fairly clearly high seismic velocity subducting Indo-Australian and low velocity anomaly around volcano regions. The relocated hypocenters seem to cluster around the local fault system such as back-arc thrust fault in northern part of the region and around local fault in Sumbawa regions. Our local earthquake tomography results demonstrated consistent with previous studies and improved the resolution. For future works, we will determine S-wave velocity structure using S-wave arrival time to enhance our understanding of geological processes and for much better interpretation.

58 GEOSCIENCES↗

Integrated pre-feasibility study for CO 2 geological storage in the Cascadia Basin, offshore Washington State, British Columbia

The Cascadia CarbonSAFE project conducted a pre-feasibility study to evaluate technical and nontechnical aspects of storing 50 million metric tons (MMT) of carbon dioxide in a safe, ocean basalt reservoir offshore Washington State and British Columbia. Sub-seafloor basalts are very common on Earth and enable geological mineralization as a long-term storage mechanism, permanently sequestering the carbon in solid rock form. This project evaluated the offshore storage complex, developed potential industrial source/transport scenarios, built an inventory of existing geophysical/geological data and environmental monitoring capabilities, conducted laboratory studies and reservoir modeling to determine storage capacity, and analyzed economic factors, regulatory requirements, and project management risks. Our team included researchers at Columbia University in New York, Pacific Northwest National Laboratory (PNNL) in Washington, University of Victoria in British Columbia, University of California at Santa Cruz, and GHG Underground in Maine, as well as collaborators at University of Iceland in Reykjavik and University of Washington in Saint Louis. In the study region, large emitters (>100,000 MT/year) generate a total of approximately 40 MMT annually from stationary sources, such as power plants, refineries, ammonia production, and mineral processing plants. Five potential industry source/transport scenarios were identified (three in U.SA S. and two in Canada) to potentially provide 50 MMT of carbon dioxide to the offshore reservoir. One scenario included net carbon-negative sources that would reduce atmospheric levels. The inventory of existing data illustrated that basalt properties in the region are potentially beneficial for long-term storage. Permeability is on the order of 0.1 to 1 Darcies within the uppermost 600 m of ocean crust and density-derived porosity reaches 10–20% in thin layers that provide act as flow channels and provide access to porous and permeable basalt. Fine-grained sediments overlying basalt provide a low-permeability barrier separating the storage reservoir from the overlying ocean. Subsurface monitoring of the storage complex benefits from a unique opportunity close proximity of to the nearby NEPTUNE cabled seafloor observatory, which actively transmits data to shore and could be used to monitor reservoir properties. Injection pressure changes, reservoir chemistry, seismicity, well head corrosion, and potential leakage monitoring would be feasible in real time. Using existing data and published models, the STOMP-CO 2 modeling code developed at PNNL was used to simulate carbon dioxide injection, storage, and mineralization in the reservoir. Preliminary simulations indicate that a 50 MMT carbon dioxide plume injected over a 20-year period remains well within the reservoir area, both laterally and vertically, for at least a 50-year period of post-injection monitoring. Laboratory studies on core samples from Cascadia measured dissolution rates in CO 2 –basalt–seawater mixtures and indicated large variability in rates with the potential for complete mineralization. Compared to published data, the dissolution rate measured in these samples is faster, decreasing the time until full conversion to carbonate minerals would occur. Adding geochemical reactivity in the STOMP-CO 2 model indicates that carbon dioxide would be fully converted to carbonate minerals in 135 years or less after injection ceases. The project assessed the viability and costs for the five scenarios using either dedicated offshore pipelines or shipping vessels. Carbon dioxide capture costs were not considered. Shipping incurred lower upfront CAPEX costs, and because most viable sources were located near the shore and waterways, shipping is preferred; no dedicated carbon dioxide pipelines exist in the region. The project also considered potential risks from mixed-use impacts on surrounding lands and ocean areas, such as active use by commercial and recreational fisherman, and other community considerations in planning for public engagement in the future. Other project risks were identified, including project management, permitting and compliance, and financing, as well as technical issues such as long-term instrumentation, plume migration, and model development. Geological concerns like reservoir permeability, seal, and induced seismicity held relatively low risk. The regulatory framework for offshore storage in the U.S. and Canada may be covered by general legal/regulatory programs, but is not yet fully defined. Existing programs appear to restrict the transportation of carbon dioxide across national boundaries (e.g., Canada to U.S.) and its sub-seabed disposal in some circumstances. Our assessment suggests that developing specific and separate pilot projects will help to advance the technical feasibility and regulatory requirements in each country, and lessons learned from such activities may be transferable elsewhere around the globe.

42 ENGINEERING↗

Poromechanics of stick-slip frictional sliding and strength recovery on tectonic faults

Pore fluids influence many aspects of tectonic faulting including frictional strength aseismic creep and effective stress during the seismic cycle. But, the role of pore fluid pressure during earthquake nucleation and dynamic rupture remains poorly understood. Here we report on the evolution of pore fluid pressure and porosity during laboratory stick-slip events as an analog for the seismic cycle. We sheared layers of simulated fault gouge consisting of glass beads in a double-direct shear configuration under true triaxial stresses using drained and undrained fluid conditions and effective normal stress of 5–10 MPa. Shear stress was applied via a constant displacement rate, which we varied in velocity step tests from 0.1 to 30 µm/s. Here, we observe net pore pressure increases, or compaction, during dynamic failure and pore pressure decreases, or dilation, during the interseismic period, depending on fluid boundary conditions. In some cases, a brief period of dilation is attendant with the onset of dynamic stick slip. Our data show that time-dependent strengthening and dynamic stress drop increase with effective normal stress and vary with fluid conditions. For undrained conditions, dilation and preseismic slip are directly related to pore fluid depressurization; they increase with effective normal stress and recurrence time. Microstructural observations confirm the role of water-activated contact growth and shear-driven elastoplastic processes at grain junctions. These results indicate that physicochemical processes acting at grain junctions together with fluid pressure changes dictate stick-slip stress drop and interseismic creep rates and thus play a key role in earthquake nucleation and rupture propagation.

58 GEOSCIENCES↗

Unsupervised learning from three-component accelerometer data to monitor the spatiotemporal evolution of meso-scale hydraulic fractures

Enhanced geothermal systems can provide a substantial share of the global energy demand. There exist several hurdles in the engineering implementations of such geothermal systems. One such hurdle is the accurate monitoring of the fracture networks created in subsurface through hydraulic stimulation of these systems. Micro seismicity associated with the stimulation is the primary means to locate the event hypocenters for estimating the stimulated rock volume. Existing methods for location the hypocenters are restricted to only the highest amplitude impulsive signals that are simultaneously detected on several sensors. Consequently, a large portion (usually ~99%) of the measurements are left unused. In this paper, an unsupervised manifold-approximation followed by clustering of 3-component accelerometer data is used to analyze the seismicity recorded on a monitoring well. With this method, a larger portion of the measured signal is used for the monitoring of the hydraulic fracture network. We analyze the EGS Collab experiment 1 microseismic data, recorded at the Sanford Underground Research Facility, South Dakota. Using the data from a single three-component accelerometer, the polarization features viz. Azimuth, incidence, rectilinearity, and planarity are used as inputs for the unsupervised manifold approximation followed by clustering. Our study shows that density-based clusters in the projected 3D space correspond to distinct types of hydraulically fractured zones around the injection point. Finally, we show that the temporal evolution of these clusters can be used to track fracture creation and propagation.

58 GEOSCIENCES↗

Quantitative interpretation of time-lapse seismic data at Farnsworth field unit: Rock physics modeling, and calibration of simulated time-lapse velocity responses

Here, this study investigates the contribution of fluid saturation variation to the time-lapse velocity response by performing fluid substitution modeling. The methodology is exemplified by the time-lapse seismic monitoring of carbon dioxide at Farnsworth field unit (FWU). In order to evaluate the fluid distribution in a matured oil reservoir, the Southwest Regional Partnership (SWP) acquired multiple vertical seismic profile (VSP) surveys at different times during the CO 2 –water alternatinggas (WAG) injection period. In this work, we present a thorough methodology for computing the elastic response of the saturated rock for different fluid saturations using a site-specific petro-elastic model (PEM). The output from the PEM was combined with results from a fluid compositional model to compute the seismic velocities at times corresponding to each VSP survey. To produce a calibrated simulated response, the measured time-lapse seismic velocities were integrated into the numerical simulation model. The mismatches between the predicted and measured time-lapse velocities were minimized through an iterative calibration process using a trained artificial neural network proxy (ANN) coupled with a particle swarm optimizer (PSO). Our study indicates that the hybrid optimization workflow can effectively perform the history matching. With an accurate prediction of the hydrodynamic properties, the migration of CO 2 within the subsurface was modeled by predicting the spatial velocity distribution for a radius of 305 m around the injection well. The technology demonstrated and the expertise gained from this study can guide similar CO 2 -WAG projects.

58 GEOSCIENCES↗

The Local Seismoacoustic Wavefield of a Research Nuclear Reactor and Its Response to Reactor Power Level

In this article, we describe the seismoacoustic wavefield recorded outdoors but inside the facility fence of the High Flux Isotope Reactor (HFIR) at Oak Ridge National Laboratory (Tennessee). HFIR is a research nuclear reactor that generates neutrons for scattering, irradiation research, and isotope production. This reactor operates at a nominal power of 85 MW, with a full-power period between 24 and 26 days. This study uses data from a single seismoacoustic station that operated for 60 days and sampled a full operating reactor cycle, that is, full-power operation and end-of-cycle outage. The analysis presented here is based on identifying signals that characterize the steady, that is, full-power operation and end-of-cycle outage, and transitional, that is, start-up and shutdown, states of the reactor. We found that the overall seismoacoustic energy closely follows the main power cycle of the reactor and identified spectral regions excited by specific reactor operational conditions. In particular, we identified a tonal noise sequence with a fundamental frequency around 21.4 Hz and multiple harmonics that emerge as the reactor reaches 90% of nominal power in both seismic and acoustic channels. We also utilized temperature measurements from the monitoring system of the reactor to suggest links between the operation of reactor’s subsystems and seismoacoustic signals. We demonstrate that seismoacoustic monitoring of an industrial facility can identify and track some industrial processes and detect events related to operations that involve energy transport.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Frequency domain analysis of errors in cross-correlations of ambient seismic noise

We analyse random errors (variances) in cross-correlations of ambient seismic noise in the frequency domain, which differ from previous time domain methods. Extending previous theoretical results on ensemble averaged cross-spectrum, we estimate confidence interval of stacked cross-spectrum of finite amount of data at each frequency using non-overlapping windows with fixed length. The extended theory also connects amplitude and phase variances with the variance of each complex spectrum value. Analysis of synthetic stationary ambient noise is used to estimate the confidence interval of stacked cross-spectrum obtained with different length of noise data corresponding to different number of evenly spaced windows of the same duration. This method allows estimating Signal/Noise Ratio (SNR) of noise cross-correlation in the frequency domain, without specifying filter bandwidth or signal/noise windows that are needed for time domain SNR estimations. Based on synthetic ambient noise data, we also compare the probability distributions, causal part amplitude and SNR of stacked cross-spectrum function using one-bit normalization or pre-whitening with those obtained without these pre-processing steps. Natural continuous noise records contain both ambient noise and small earthquakes that are inseparable from the noise with the existing pre-processing steps. Using probability distributions of random cross-spectrum values based on the theoretical results provides an effective way to exclude such small earthquakes, and additional data segments (outliers) contaminated by signals of different statistics (e.g. rain, cultural noise), from continuous noise waveforms. This technique is applied to constrain values and uncertainties of amplitude and phase velocity of stacked noise cross-spectrum at different frequencies, using data from southern California at both regional scale (~35 km) and dense linear array (~20 m) across the plate-boundary faults. A block bootstrap resampling method is used to account for temporal correlation of noise cross-spectrum at low frequencies (0.05–0.2 Hz) near the ocean microseismic peaks.

58 GEOSCIENCES↗

Central and Eastern North America Ground-Motion Characterization (NGA-East)

This document is the final project report of the Next Generation Attenuation for Central and Eastern North America (CENA) project (NGA-East). The NGA-East objective was to develop a new ground-motion characterization (GMC) model for the CENA region. The GMC model consists of a set of new ground-motion models (GMMs) for median and standard deviation of ground motions and their associated weights to be used with logic-trees in probabilistic seismic hazard analyses (PSHA). NGA-East is a large multidisciplinary project coordinated by the Pacific Earthquake Engineering Research Center (PEER), at the University of California. The project has two components: (1) a set of scientific research tasks, and (2) a model-building component following the framework of the “Seismic Senior Hazard Analysis Committee (SSHAC) Level 3” (Budnitz et al. 1997; NRC 2012). Component (2) is built on the scientific results of component (1) of the NGA-East project. This report documents the tasks under component (2) of the project. Under component (1) of NGA-East, several scientific issues were addressed, including: (a) development of a new database of ground motion data recorded in CENA; (b) development of a regionalized ground-motion map for CENA, (c) definition of the reference site condition; (d) simulations of ground motions based on different methodologies; and (e) development of numerous GMMs for CENA. The scientific tasks of NGA-East were all documented as a series of PEER reports. The scope of component (2) of NGA-East was to develop the complete GMC. This component was designed as a SSHAC Level 3 study with the goal of capturing the ground motions’ center, body, and range of the technically defensible interpretations in light of the available data and models. The SSHAC process involves four key tasks: evaluation, integration, formal review by the Participatory Peer Review Panel (PPRP), and documentation (this report). Key tasks documented in this report include review and evaluation of the empirical groundmotion database, the regionalization of ground motions, and screening sets of candidate GMMs. These are followed by the development of new median and standard deviation GMMs, the development of new analyses tools for quantifying the epistemic uncertainty in ground motions, and the documentation of implementation guidelines of the complete GMC for PSHA computations. Appendices include further documentation of the relevant SSHAC process and additional supporting technical documentation of numerous sensitivity analyses results. The PEER reports documenting component (1) of NGA-East are also considered “attachments” to the current report and are all available online on the PEER website (https://peer.berkeley.edu/). The final NGA-East GMC model includes a set of 17 GMMs defined for 24 ground-motion intensity measures, applicable to CENA in the moment magnitude range of 4.0 to 8.2 and covering distances up to 1500 km. Standard deviation models are also provided for site-specific analysis (single-station standard deviation) and for general PSHA applications (ergodic standard deviation). Adjustment factors are provided for consideration of source-depth effects and hanging-wall effects, as well as for hazard computations at sites in the Gulf Coast region.

58 GEOSCIENCES↗

Accelerating full waveform inversion via source stacking and cross-correlations

SUMMARY Accurate synthetic seismic wavefields can now be computed in 3-D earth models using the spectral element method (SEM), which helps improve resolution in full waveform global tomography. However, computational costs are still a challenge. These costs can be reduced by implementing a source stacking method, in which multiple earthquake sources are simultaneously triggered in only one teleseismic SEM simulation. One drawback of this approach is the perceived loss of resolution at depth, in particular because high-amplitude fundamental mode surface waves dominate the summed waveforms, without the possibility of windowing and weighting as in conventional waveform tomography. This can be addressed by redefining the cost-function and computing the cross-correlation wavefield between pairs of stations before each inversion iteration. While the Green’s function between the two stations is not reconstructed as well as in the case of ambient noise tomography, where sources are distributed more uniformly around the globe, this is not a drawback, since the same processing is applied to the 3-D synthetics and to the data, and the source parameters are known to a good approximation. By doing so, we can separate time windows with large energy arrivals corresponding to fundamental mode surface waves. This opens the possibility of designing a weighting scheme to bring out the contribution of overtones and body waves. It also makes it possible to balance the contributions of frequently sampled paths versus rarely sampled ones, as in more conventional tomography. Here we present the results of proof of concept testing of such an approach for a synthetic 3-component long period waveform data set (periods longer than 60 s), computed for 273 globally distributed events in a simple toy 3-D radially anisotropic upper mantle model which contains shear wave anomalies at different scales. We compare the results of inversion of 10 000 s long stacked time-series, starting from a 1-D model, using source stacked waveforms and station-pair cross-correlations of these stacked waveforms in the definition of the cost function. We compute the gradient and the Hessian using normal mode perturbation theory, which avoids the problem of cross-talk encountered when forming the gradient using an adjoint approach. We perform inversions with and without realistic noise added and show that the model can be recovered equally well using one or the other cost function. The proposed approach is computationally very efficient. While application to more realistic synthetic data sets is beyond the scope of this paper, as well as to real data, since that requires additional steps to account for such issues as missing data, we illustrate how this methodology can help inform first order questions such as model resolution in the presence of noise, and trade-offs between different physical parameters (anisotropy, attenuation, crustal structure, etc.) that would be computationally very costly to address adequately, when using conventional full waveform tomography based on single-event wavefield computations.

Romanowicz, Barbara↗