Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “seismic location”

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 127 records · Page 7

A Methodological Overview of Seismic Analysis for Nuclear Event Detection

Underground explosions generate potentially detectable signatures, including energy waves that travel through the Earth’s subsurface (i.e., seismic waves), low-frequency sound waves (i.e., infrasound and hydroacoustic waves), and radioactive gases and/or particles that might leak from the test cavity (if the event was nuclear). There can also be intelligence indicators of a test, such as observations of modified patterns of life and activity at a suspected test site. If all of these detectable signatures and intelligence indicators are present and self-consistent, then analysts have high confidence in classifying a signature generating event as an explosion. However, because only partial information about an event is likely to be available, determining whether an event was natural (e.g., an earthquake or landslide) or manmade (e.g., a chemical or nuclear explosion) is much more challenging. This primer describes how one category of event signatures—seismic signatures—can augment event analyses. While universities and government organizations have generated detailed technical descriptions of seismic analytic techniques, we seek to translate seismic event analysis for a broad, non-technical audience. When the geologic conditions near an event are well-characterized, seismic data can be used to calculate critical information, such as event location and depth, with relatively high accuracy. Moreover, specific features within seismic datasets can help determine whether an event was an explosion. However, a key challenge in seismic analysis is that geologic site conditions are often poorly characterized, complicating the ability to discern the true nature of the event. To overcome this challenge, geologists answer a series of questions (discussed in section 1) to guide seismic event analysis and determine the most probable nature of an event. As more information is gathered during each analytic step, confidence grows regarding the nature of the event. Section 2 addresses uncertainties in seismic analysis and the vital nature of high-fidelity geologic data for accurate seismic event analysis.

58 GEOSCIENCES↗

Pore Pressure Diffusion and Onset of Induced Seismicity

The Raton Basin of Southern Colorado and Northern New Mexico experienced a significant increase in seismicity over the last several decades, including two M5+ earthquakes. The increase in seismicity started soon after an increase in wastewater disposal associated with coal-bed methane production. Twenty-nine wastewater injection wells have operated in the Basin, with several high-rate injectors located within a kilometer of the most seismically active regions. We developed a pore pressure model to investigate the relationship between pore pressure diffusion and onset of induced seismicity, with an emphasis on the Trinidad fault zones. Our model indicates that sufficient pore pressure increase had accumulated on the northern extent of the Trinidad fault zone to induce seismicity at the time of early-stage seismic activity. The hydrogeologic architecture of the Basin played a key role in controlling pore pressure perturbation. The basal sedimentary interval with moderately low diffusivity restricted vertical pore pressure diffusion, creating a substantial lag time between wells becoming active and the expression of perturbation in the fault zones. The spatiotemporal relationship between early-stage seismicity and wastewater injection in conjunction with modeled pore pressures suggests that the early-stage seismicity in the Raton Basin was likely triggered by pore pressure diffusion.

58 GEOSCIENCES↗

Updates to the Regional Seismic Travel Time (RSTT) Model: 2. Path-dependent Travel-time Uncertainty

Abstract The regional seismic travel time (RSTT) model and software were developed to improve travel-time prediction accuracy by accounting for three-dimensional crust and upper mantle structure. Travel-time uncertainty estimates are used in the process of associating seismic phases to events and to accurately calculate location uncertainty bounds (i.e. event location error ellipses). We improve on the current distance-dependent uncertainty parameterization for RSTT using a random effects model to estimate slowness (inverse velocity) uncertainty as a mean squared error for each model parameter. The random effects model separates the error between observed slowness and model predicted slowness into bias and random components. The path-specific travel-time uncertainty is calculated by integrating these mean squared errors along a seismic-phase ray path. We demonstrate that event location error ellipses computed for a 90% coverage ellipse metric (used by the Comprehensive Nuclear-Test-Ban Treaty Organization International Data Centre (IDC)), and using the path-specific travel-time uncertainty approach, are more representative (median 82.5% ellipse percentage) of true location error than error ellipses computed using distance-dependent travel-time uncertainties (median 70.1%). We also demonstrate measurable improvement in location uncertainties using the RSTT method compared to the current station correction approach used at the IDC (median 74.3% coverage ellipse).

58 GEOSCIENCES↗

Lunar seismicity and tectonics

Results are presented for an analysis of all moonquake data obtained by the Apollo seismic stations during the period from November 1969 to May 1974 and a preliminary analysis of critical data obtained in the interval from May 1974 to May 1975. More accurate locations are found for previously located moonquakes, and additional sources are located. Consideration is given to the sources of natural seismic signals, lunar seismic activity, moonquake periodicities, tidal periodicities in moonquake activity, hypocentral locations and occurrence characteristics of deep and shallow moonquakes, lunar tidal control over moonquakes, lunar tectonism, the locations of moonquake belts, and the dynamics of the lunar interior. It is concluded that: (1) moonquakes are distributed in several major belts of global extent that coincide with regions of the youngest and most intense volcanic and tectonic activity; (2) lunar tides control both the small quakes occurring at great depth and the larger quakes occurring near the surface; (3) the moon has a much thicker lithosphere than earth; (4) a single tectonic mechanism may account for all lunar seismic activity; and (5) lunar tidal stresses are an efficient triggering mechanism for moonquakes.

Lammlein, D. R.↗

On component isolation of conceptual advanced reactors

Implementation of component isolation in nuclear industry is challenging due to gaps in research and the lack of specific guidelines. In this study, parameters affecting component-level isolation of advanced reactor vessels are identified based on a representative numerical model with explicit consideration of nonlinear soil-structure interaction (SSI). The objective of this study is to evaluate the effectiveness of, and to identify potential limitations of using conventional friction pendulum and lead-rubber bearings to seismically isolate vessels. It is found that slender vessels or components are particularly vulnerable to rotational accelerations at the isolation interface, which are caused by rotation of the sub-structure and by excitation of higher modes in the horizontal direction of the seismically isolated system. Component isolation is found to be more effective for relatively stiffer vessels and at sites with stiff soil properties. Considering that conventional isolators are deficient in resisting axial tension, it is observed that the optimum location for supporting a component to achieve seismic isolation, is at a cross-sectional plane passing through the center of mass of the vessel. These findings are corroborated by numerous simulations of the response of seismically isolated reactor vessels at different nuclear sites subject to a variety of ground motions.

99 GENERAL AND MISCELLANEOUS↗

Emulation of seismic-phase traveltimes with machine learning

SUMMARY We present a machine learning (ML) method for emulating seismic-phase traveltimes that are computed using a global-scale 3-D earth model and physics-based ray tracing. Accurate traveltime predictions based on 3-D earth models are known to reduce the bias of event location estimates, increase our ability to assign phase labels to seismic detections and associate detections to events. However, practical use of 3-D models is challenged by slow computational speed and the unwieldiness of pre-computed lookup tables that are often large and have prescribed computational grids. In this work, we train a ML emulator using pre-computed traveltimes, resulting in a compact and computationally fast way to approximate traveltimes that are based on a 3-D earth model. Our model is trained using approximately 850 million P-wave traveltimes that are based on the global LLNL-G3D-JPS model, which was developed for more accurate event location. The training-set consists of traveltimes between 10 393 global seismic stations and randomly sampled event locations that provide a prescribed, distance-dependent geographic sample density for each station. Prediction accuracy is dependent on event-station distance and whether the station was included in the training set. For stations included in the training set the mean absolute deviation (MAD) of the difference between traveltimes computed using ray tracing through the 3-D model and the ML emulator for local, regional, and teleseismic distances are 0.090, 0.125 and 0.121 s, respectively. For tested station locations not included in the training set, MAD values for the three distance ranges increase to 0.173, 0.219 and 0.210 s, respectively. Empirical traveltime residuals for a global reference data are indistinguishable when ML emulation or the 3-D model is used to compute traveltimes. This result holds regardless of whether the recording station is used in ML training or not.

58 GEOSCIENCES↗

Fourier-DeepONet: Fourier-enhanced deep operator networks for full waveform inversion with improved accuracy, generalizability, and robustness

In this article, full waveform inversion (FWI) infers the subsurface structure information from seismic waveform data by solving a non-convex optimization problem. Data-driven FWI has been increasingly studied with various neural network architectures to improve accuracy and computational efficiency. Nevertheless, the applicability of pre-trained neural networks is severely restricted by potential discrepancies between the source function used in the field survey and the one utilized during training. Here, we develop a Fourier-enhanced deep operator network (Fourier-DeepONet) for FWI with the generalization of seismic sources, including the frequencies and locations of sources. Specifically, we employ the Fourier neural operator as the decoder of DeepONet, and we utilize source parameters as one input of Fourier-DeepONet, facilitating the resolution of FWI with variable sources. To test Fourier-DeepONet, we develop three new and realistic FWI benchmark datasets (FWI-F, FWI-L, and FWI-FL) with varying source frequencies, locations, or both. Our experiments demonstrate that compared with existing data-driven FWI methods, Fourier-DeepONet obtains more accurate predictions of subsurface structures in a wide range of source parameters. Moreover, the proposed Fourier-DeepONet exhibits superior robustness when handling data with Gaussian noise or missing traces and sources with Gaussian noise, paving the way for more reliable and accurate subsurface imaging across diverse real conditions.

42 ENGINEERING↗

Microseismicity Modulation Due To Changes in Geothermal Production at San Emidio, Nevada, USA

Brief cessations of geothermal production can induce seismicity, a phenomenon that has drawn increasing attention in recent years. Such observations are rare, and the underlying mechanism requires careful analysis. In April 2022, a dense seismic and hydrologic monitoring system was deployed at the San Emidio geothermal field, Nevada, to accompany a planned power plant shutdown. Using the dense seismic array data, we detected and located ∼1,800 microseismic events (MSEs) and developed a high-resolution tomographic P-wave velocity model. We observed substantially increased microseismicity during shutdown. Most MSEs occurred on pre-existing normal faults, which are contained within extremely low-velocity zones that are likely damaged, fluid-filled, and hydraulically connected to nearby production wells. Hydrologic data show rapid fluid pressure increases of <60 kPa following the shutdown. We suggest that the cessation of production rapidly increased fluid pressures along pre-existing fault zones, activating critically stressed fault patches and fractures and producing microseismicity.

Guo, Hao [University of Wisconsin-Madison, WI (Uni↗

Prairie State Generating Company 2D Seismic Interpretation

The objectives of the Illinois Storage Corridor (ISC) project are to accelerate commercial deployment of carbon capture utilization and storage at two individual sites and receive approvals for Underground Injection Control (UIC) Class VI permits for construction at each site (ISC Project Narrative, 2020). As part of this project, and as part of the subsurface geologic characterization, two-dimensional (2D) seismic data was acquired at both sites. This report summarizes the findings from two phases of 2D acquisition and seismic interpretation at the Prairie State Generating Company site near Marissa, Illinois. The seismic data indicates the presence of three faults that completely transect the storage and confining units. The initial 2021 2D seismic acquisition revealed the presence of a feature two miles east of the Lively Grove #1 characterization well. The feature is a polyphase fault with a component of strike-slip motion, forming a small positive flower structure. The second phase of 2D seismic acquisition in 2022 constrained the maximum extent of the fault, indicating it has a relatively limited length. Other 2022 seismic lines revealed the presence of two other faults that transect the storage and confining units in the southwest and northeast portions of the 2D seismic acquisition area. The 2021 and 2022 2D seismic surveys have identified three specific locations within the project area that may have an elevated risk of out of zone CO 2 migration due to faulting. Further technical work will be needed to quantify this risk and determine how and if this risk will impact the placement of CO 2 injection wells across the project area.

01 COAL, LIGNITE, AND PEAT↗

Lunar seismology - The internal structure of the moon

It is pointed out that seismology has provided the most detailed information concerning the structure and state of the earth's interior. Beginning in 1969, seismometers were landed on the moon by the Apollo missions, providing the first opportunity to attempt similar studies on another planetary body. In September 1977 the operation of these instruments was terminated. A description is presented of the internal structure of the moon, as determined from the obtained lunar seismic data. The analysis of the lunar data is approached in a systematic fashion, using appropriate techniques to minimize the number of necessary assumptions, extract the maximum amount of structural information, and determine its reliability. The completed lunar seismic network consists of four stations located at the landing sites of Apollo missions 12, 14, 15, and 16. Attention is given to crustal structure, the structure of the lunar mantle, the attenuating region, and the core.

Goins, N. R.↗

Characterization of the Marquez Dome buried impact crater using gravity and magnetic data

The buried impact crater, Marquez Dome, located in Leon County in east central Texas, is an approximately 15 km diameter structure whose central uplift is now partially exposed due to headward erosion of the post-impact cover. The central uplift is approximately 3 km in diameter and the rocks within it have been uplifted more than 1200 m above their regional level. The crater rim remains buried and previous attempts to determine its location have had to rely on seismic reflection data and geologic well logs. These attempts have been somewhat successful in mapping the extent of the disturbed zone around Marquez Dome, but more limited in their ability to image the shallow buried rim. In an attempt to define accurately the whole Marquez Dome structure and assist in the selection of drilling sites, a geophysical investigation involving gravity and magnetic data over the central uplift and the surrounding area has been undertaken.

Wong, A. M.↗

Coupling geophysical, geological, geochemical and mineralogical assessments to examine preferential contaminant transport pathways in interbedded fractured bedrock

This study shows that a multi-faceted approach, combining borehole geophysical logging and surface seismic P-wave first-arrival tomography with confirmatory coring, well installation, and chemical and mineralogical analysis, is effective for identifying difficult-to-locate preferential contaminant transport pathways in deeper fractured bedrock. Seismic tomography detected porous 10–20 m wide elongated fractured conduits that allow acidic groundwater contaminated with uranium (U) and nitrate (NO 3 - ) to migrate within interbedded shale-limestone bedrock over 1000 m from a former disposal facility (S-3 Ponds site) located at the DOE Y-12 National Security Complex in Tennessee (USA). Conventional drilling techniques would easily miss these conduits because they are oriented parallel with fractured bedding planes. Synchrotron analysis of aquifer solids revealed that > 95 % of the U is hexavalent (U VI ). This uranyl (UO 2 2+ ) species is coordinated with carbonate, iron oxide, silicate and phosphate minerals within cemented white to yellow precipitates, which contain U concentrations as high as ∼21.6 % by weight fraction. Identifying the presence of these mineral phases, enables a further understanding of the potential effectiveness of remediation actions. The combination of methodologies presented here can also be applied to other explorations, such as the detection of water supply.

Fate and transport↗

Gulf of Mexico Miocene CO2 Site Characterization Mega Transect

The project Gulf of Mexico Miocene CO2 Site Characterization Mega Transect is one of 9 site characterization projects that were implemented as part of ARRA (American Recovery and Reinvestment Act). Data from this project was used to improve resolution of data in NATCARB in the area of study. Data related to this study has already been incorporated in NATCARB Atlas. The University of Texas at Austin study focused on (1) collecting high-resolution, 3-D (HR3D) seismic datasets of the shallow geologic section to evaluate the capacity of potential reservoirs and the sealing capacity of the confining system, and (2) producing a regional CO2 atlas of Miocene age units of the upper Texas coast submerged lands (southeast Texas). The goal of the study was to evaluate the capability of the Miocene age geologic section of Texas submerged lands to permanently store large volumes of anthropogenic CO2. The geographic locations of three HR3D, “P-Cable” seismic surveys collected in the near-shore waters of southeast Texas offers a vertical and horizontal view (data “cube”) example from the second survey. The dataset's high quality allowed for defining recognizable geologic morphologies (shapes) and possible fluid content, which are important for determining an area’s prospects for geologic storage. The regional CO2 atlas of Miocene age units highlights several topics, including: (1) a regional analyses of petroleum systems as analogs for CO2 storage; (2) petrography, petrology, and extent of confining systems of the Miocene section; (3) static capacity estimate; and (4) examples of characterization methodologies of prospect areas. The HR3D datasets and regional analyses within the CO2 atlas of the offshore Texas Miocene provide a sound basis for a future generation of specific CO2 prospects in the study area.

.las↗

Geothermal Play-Fairway Analysis of Washington State Prospects: Final Report

The Washington State Geothermal Play-Fairway Analysis overcomes the exploration challenges posed by dense vegetation, glacial deposits, and extreme precipitation. The geothermal play-fairways we target are locations where heat, permeability, and saturated porosity are present in sufficient volume to provide adequate heat exchange at depths accessible by modern drilling technology. The three study areas lie along the Cascade Range magmatic arc and are near Mount Baker, Mount St. Helens, and the Wind River Valley. The seven-year project is divided into three phases. In Phase 1 we build on a previous statewide assessment of geothermal resources and develop an initial modeling approach. The results are a series of favorability, uncertainty, and risk maps for three targeted study areas. Based on these initial results, we collect new geologic and geophysical data to further refine our modeling and reduce exploration uncertainty in Phase 2. We improve the modeling method to handle the new data and update the favorability, uncertainty, and risk maps. We also update the conceptual geothermal resource models. In Phase 3 we validate our modeling approach by drilling two temperature-gradient holes and collecting and analyzing core, image logs, and new geochemistry. Our modeling approach improves on an earlier statewide method through a more-rigorous and detailed assessment of heat and permeability. Permeability potential is assessed through geomechanical modeling of the deformation that can generate and maintain reservoir porosity and permeability. Metrics to inform heat potential include temperature-gradient wells, which are sparse in Washington; proximity of Quaternary volcanic vents and young intrusive rock; spring temperature; and reservoir temperature inferred from geothermometry. We weight the individual components using an expert-guided approach known as the Analytical Hierarchy Process. During Phase 2 we also develop a fluid-filled fracture model, and an infrastructure model that helps to delineate areas which are more favorable for geothermal development based on proximity to transmission lines, elevation, land ownership and use restrictions, and availability of process water. New geologic and geophysical data is collected during Phase 2 in each of our three main study areas. At Mount Baker and north of Mount St. Helens we conduct 1:24,000-scale geologic mapping and lidar analysis to better constrain the location and character of surface faults; detailed mapping in the Wind River Valley was completed just prior to the start of this project. Ages of intrusive rocks are determined with 40 Ar/ 39 Ar geochronology, though all of our samples are Miocene or older. We collect ground based gravity observations (a total of 1,580 new stations) in all of our study areas and ground-based magnetic lines (a total of 93 km) at Mount Baker. These data are combined with existing gravity and aeromagnetic data and used to constrain fault locations and geometry. Two to three cross sections are constructed at each study area using the mapped surface geology and forward-modeling of the gravity and magnetic data; these cross sections form the basis for our updated conceptual models. We collect magnetotelluric surveys at Mount Baker and Mount St. Helens and these data are inverted to form a resistivity model from the surface to about 10 km depth; each model shows conductive zones that can be interpreted as upwelling geothermal fluids. At Mount St. Helens we deploy a passive seismic array and use the newly detected events to refine the location of the Saint Helens seismic zone. We also employ ambient-noise tomography to develop a detailed seismic-velocity model for the study area and use this model to help constrain our cross sections and conceptual model. Based on the new data collected during Phase 2—and our updated models—we develop a campaign of temperature-gradient holes and core analysis to validate our modeling in Phase 3. Drill hole MB76-31 is located near Little Park Creek, 11 km west-southwest of the summit of Mount Baker, and is 1,471 ft deep. About 410 ft of core from the lower portion of the hole—and image logs from ~175 ft below ground surface to the bottom—are collected and analyzed. Water samples are collected and processed for geothermometry. Drill hole MSH17-24 is located along upper Schultz Creek, 16 km north-northeast of Mount St. Helens and has core from 470 ft to the bottom at 1,053 ft. We did not collect image logs due to borehole stability concerns, but water samples are collected and analyzed for geothermometry. Repeat temperature-gradient measurements are made at both sites and thermal conductivity is measured from core samples. At MB76-31, the equilibrated temperature gradient of 64°C/km and calculated heat flow of 141–159 mW/m 2 is more than twice the regional average. Detailed mapping and analysis of the core, coupled with correlation to the image logs, indicates a history of permeability generation consistent with our predictions of high permeability. Because the site has high favorability in the Phase 2 model, we consider the results a positive validation of the modeling. At site MSH17-24, the equilibrated temperature gradient of ~15°C/km and calculated heat flow of 41–43 mW/m 2 are similar to regional. Geochemical analysis of the water samples indicates a meteoric source without any geothermal component. Detailed outcrop-based mapping of fault exposures near the drill site and analysis of image logs from nearby boreholes indicates a history of permeability generation consistent with our predictions. Because the site has low favorability in the Phase 2 model, we consider the results a positive validation of the modeling. Together, the two sites provide a reasonably positive validation of the Phase 2 modeling and should encourage future use of this modeling approach.

15 GEOTHERMAL ENERGY↗

Lessons Learned from Recent Rapid Imaging of Earthquake Ruptures with Satellite Geodesy

Rapid determination of the location and extent of earthquake ruptures at the surface and at depth is helpful for disaster response, as it allows prediction of the likely area of major damage from the earthquake and can help with rescue and recovery planning. The Caltech-Jet Propulsion Laboratory (JPL) Advanced Rapid Imaging and Analysis (ARIA) project has responded to many recent large earthquakes to process geodetic data and help determine the location and extent of ruptures. With the increasing availability of near real-time data from the Global Positioning System (GPS) and other global navigation satellite system receivers in active tectonic regions, and with the shorter repeat times of many recent and newly launched radar satellites, geodetic data can now be obtained quickly after earthquakes or other disasters. We have been building an ARIA data system that can ingest, catalog, and process geodetic data and combine it with seismic analysis to estimate the fault rupture locations and slip distributions for large earthquakes that are on or near land.

Simons, Mark↗

Crustal Deformation and the Seismic Cycle across the Kodiak Islands, Alaska

The Kodiak Islands are located approximately 130 to 250 km from the Alaska-Aleutian Trench where the Pacific plate is underthrusting the North American plate at a rate of about 57 mm/yr. The southern extent of the 1964 Prince William Sound (${M-w}$ = 9.2) earthquake rupture occurred offshore and beneath the eastern portion of the Kodiak Islands. Here we report GPS results (1993-2001) from northern Kodiak Island that span the transition between the 1964 uplift region along the eastern coast and the region of coseismic subsidence further inland. The horizontal velocity vectors range from 22.9 $\pm$ 2.2 mm/yr at N26.3$\deg$W $\pm$ 2.5$\deg$, about 150 km from the trench, to 5.9 $\pm$ 1.3 mm/yr at N65.9$\deg$W $\pm$ 6.6$\deg$, about 190 km from the trench. Near the northeastern coast of Kodiak the velocity vector above the shallow, locked main thrust zone is between the orientation of PCFC-NOAM plate motion (N22$/deg$W) and the trench-normal (N3O$\deg$W). Further west, our geodetic results suggest the accumulation of shear strain that will be released eventually as left-lateral motion on upper plate faults such as the Kodiak Island fault. These results are consistent with the hypothesis that the difference between the Pacific-North American plate motion and the orientation of the down going slab would lead to 4-8 mm/yr of left-lateral slip. Short-term geodetic uplift rates range from 2 - 14 mm/yr, with the maximum uplift located near the axis of maximum subsidence during the 1964 earthquake. We evaluated alternate interseismic models for Kodiak to test the importance of various mechanisms responsible for crustal deformation rates. These models are based on the plate interface slip history inferred from earlier modeling of coseismic and post-seismic geodetic results. The horizontal (trench perpendicular) and vertical deformation rates across Kodiak are consistent with a model that includes the viscoelastic response to : (1) a downgoing Pacific plate interface that is locked at shallow depths,(2) coseismic slip in the 1964 and (3) interseismic creep below the seismogenic zone. The change in orientation of the horizontal velocity vector occurs down-dip from the locked main thrust zone. In southern Kodiak, the coseismic slip in the 1964 earthquake was smaller than in the northern Kodiak region; yet, the horizontal, interseismic velocities as a function of distance from the trench are comparable to those in northern Kodiak. Based on the earthquake history prior to, and following the 1964 earthquake, we hypothesize that the plate interface in southern Kodiak slips in more frequent large earthquakes than in northern Kodiak.

Sauber, Jeanne↗

Local Explosion Detection and Infrasound Localization by Reverse Time Migration Using 3-D Finite-Difference Wave Propagation

Infrasound data are routinely used to detect and locate volcanic and other explosions, using both arrays and single sensor networks. However, at local distances (< km) topography often complicates acoustic propagation, resulting in inaccurate acoustic travel times leading to biased source locations when assuming straight-line propagation. Here we present a new method, termed Reverse Time Migration-Finite-Difference Time Domain (RTM-FDTD), that integrates numerical modeling into the standard RTM back-projection process. Travel time information is computed across the entire potential source grid via FDTD modeling to incorporate the effects of topography. The waveforms are then back-projected and stacked at each grid point, with the stack maximum corresponding to the likely source. We apply our method to three volcanoes with different network configurations, source-receiver distances, and topography. At Yasur Volcano, Vanuatu, RTM-FDTD locates explosions within ~20 m of the source and differentiates between multiple vents. RTM-FDTD produces a more accurate location for the two Yasur subcraters than standard RTM and doubles the number of detected events. At Sakurajima Volcano, Japan, RTM-FDTD locates the source within 50 m of the active vent despite notable topographic blocking. The RTM-FDTD location is similar to that from the Time Reversal Mirror method, but is more computationally efficient. Lastly, at Shishaldin Volcano, Alaska, RTM and RTM-FDTD both produce realistic source locations (<50 m) for ground-coupled airwaves recorded on a four-station seismic network. We show that RTM is an effective method to detect and locate infrasonic sources across a variety of scenarios, and by integrating numerical modeling, RTM-FDTD produces more accurate source locations and increases the detection capability.

58 GEOSCIENCES↗

Spatiotemporal Graph Convolutional Networks for Earthquake Source Characterization

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

58 GEOSCIENCES↗