Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “seismic”

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 145 records · Page 8
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Mid-height seismic isolation of equipment in nuclear power plants: Numerical simulations and design recommendations

Seismic isolation can mitigate earthquake demands on safety-class equipment and is being considered for application to some advanced nuclear reactors. For tall, slender vessels that could represent an advanced reactor, a steam generator or a heat exchanger, mid-height seismic isolation has been shown to be beneficial and practical. Here, this paper focuses on numerical modeling of a mid-height isolated tall, slender vessel with two primary goals 1) investigate and quantify the benefits of mid-height seismic isolation for a range of support structure stiffnesses, isolation systems, and seismic inputs, and 2) provide analysis and design recommendations for isolation of safety-class equipment. Results of response-history analysis were in good agreement with experimental measurements and demonstrated that mid-height isolation can substantially reduce seismic demands on tall, slender vessels for a range of support structure stiffnesses, isolation systems, and seismic inputs. Importantly, the reductions in horizontal spectral accelerations in the mid-height isolated vessel from the non-isolated condition were not affected by the stiffness of the support structure. Recommendations are made for analysis and design of isolated equipment and testing of isolators used for equipment protection in nuclear facilities, which are also applicable to other industries.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
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 OSTI2022

Review: Induced Seismicity During Geoenergy Development—A Hydromechanical Perspective

The basic triggering mechanism underlying induced seismicity traces back to the mid-1960s that relied on the process of pore-fluid pressure diffusion. The last decade has experienced a renaissance of induced seismicity research and data proliferation. An unprecedent opportunity is presented to us to synthesize the robust growth in knowledge. The objective of this article is to provide a concise review of the triggering mechanisms of induced earthquakes with a focus on hydro-mechanical processes. Four mechanisms are reviewed: pore-fluid pressure diffusion, poroelastic stress, Coulomb static stress transfer, and aseismic slip. For each, an introduction of the concept is presented, followed by case studies. Diving into these mechanisms sheds light on several outstanding questions. For example, why did some earthquakes occur far from fluid injection or after injection stopped? Our review converges on the following conclusions: (a) Pore-fluid pressure diffusion remains a basic mechanism for initiating inducing seismicity in the near-field. (b) Poroelastic stresses and aseismic slip play an important role in inducing seismicity in regions beyond the influence of pore-fluid pressure diffusion. (c) Coulomb static stress transfer from earlier seismicity is shown to be a viable mechanism for increasing stresses on mainshock faults. (d) Multiple mechanisms have operated concurrently or consecutively at most induced seismicity sites. (e) Carbon dioxide injection is succeeding without inducing earthquakes and much can be learned from its success. Future research opportunities exist in deepening the understanding of physical and chemical processes in the nexus of geoenergy development and fluid motion in the Earth’s crust.

58 GEOSCIENCES↗
DOE OSTI2025

Analysis and optimization of seismic monitoring networks with Bayesian optimal experimental design

SUMMARY Monitoring networks increasingly aim to assimilate data from a large number of diverse sensors covering many sensing modalities. Bayesian optimal experimental design (OED) seeks to identify data, sensor configurations or experiments which can optimally reduce uncertainty and hence increase the performance of a monitoring network. Information theory guides OED by formulating the choice of experiment or sensor placement as an optimization problem that maximizes the expected information gain (EIG) about quantities of interest given prior knowledge and models of expected observation data. Therefore, within the context of seismo-acoustic monitoring, we can use Bayesian OED to configure sensor networks by choosing sensor locations, types and fidelity in order to improve our ability to identify and locate seismic sources. In this work, we develop the framework necessary to use Bayesian OED to optimize a sensor network’s ability to locate seismic events from arrival time data of detected seismic phases at the regional-scale. This framework requires five elements: (i) A likelihood function that describes the distribution of detection and traveltime data from the sensor network, (ii) A prior distribution that describes a priori belief about seismic events, (iii) A Bayesian solver that uses a prior and likelihood to identify the posterior distribution of seismic events given the data, (iv) An algorithm to compute EIG about seismic events over a data set of hypothetical prior events, (v) An optimizer that finds a sensor network which maximizes EIG. Once we have developed this framework, we explore many relevant questions to monitoring such as: how to trade off sensor fidelity and earth model uncertainty; how sensor types, number and locations influence uncertainty; and how prior models and constraints influence sensor placement.

58 GEOSCIENCES↗
DOE OSTI2021

High-Quality Revision of the Israeli Seismic Bulletin

Seismic bulletins, with trustworthy phase picks, origin times, and source locations are key for regional seismic studies, such as travel-time (TT) tomography, attenuation tomography, and anisotropy studies. To lay the groundwork for such studies in Israel, we revised the seismic bulletin of Israel and the surrounding area and obtained a trustworthy TT data set. From the earthquake and explosion bulletins of the Geophysical Institute of Israel, we compiled a starting data set of about 123,000 earthquakes and explosions that occurred during the past 40 yr. After screening out the poorly recorded events, we were left with a data set of ~38,000 well-recorded events. We then revised the remaining data set in two consecutive steps. In the first, we reviewed and updated station metadata, including changes in station metadata parameters over time. In the second step, we jointly relocated a list of selected seismic events, using the Bayesian hierarchical location software package (BayesLoc) of Myers et al. (2007) that performs joint relocation of multiple events. We observed striking dissimilarities between the spatial distributions of the newly relocated catalog and the initial locations. Although the depth distribution of the starting catalog is trimodal with peaks at 0, 5, and 10 km, the distribution in this study is unimodal, with a broad peak between 7.5 and 12.5 km. By differencing the observed arrival times and the origin times obtained through relocation with BayesLoc, we obtained a revised TT database that consists of 261,336 Pg, 132,876 Pn, 114,816 Sg, and 60,394 Sn arrivals, from a set of 30,458 jointly relocated seismic sources. In this work, we compared prerevision and postrevision TTs as a function of epicentral distance and concluded that the revised data set contains far fewer outliers and inconsistencies than the original data set. The revised TT data set may be used for seismic studies, such as TT tomography, attenuation tomography, and anisotropy studies.

58 GEOSCIENCES↗
DOE OSTI2022

Reducing the Overnight Capital Cost of Advanced Reactors Using Equipment-Level Seismic Protective Systems

Consideration of the effects of earthquake shaking on the design and construction of nuclear power plants adds substantially to the overnight capital cost, with anecdotal estimates as high as 35+%, attributed to additional construction materials, need for one-off and sub-optimal designs of equipment due to conflicting design choices, the high cost of seismic qualification of equipment, and regulatory review. Safety-critical equipment in large light water reactors is designed and qualified for seismic demands imposed by the supporting reactor building, optimal mechanical designs are not possible, and designs of a given piece of equipment may vary with height above grade. Similar negative impacts are expected for advanced reactors unless the seismic design paradigm is changed. The overarching goal of this transformational MEITNER project, which involved a multidisciplinary engineering team and designers of three fundamentally different advanced reactors, was to adapt proven seismic isolation and damping technologies to operationalize modular protective systems for safety-class equipment inside advanced reactor buildings. Such seismic protective systems would be tightly integrated into design development for reactor support systems and balance-of-plant construction. The adoption of the technology, which is widely used in non-nuclear sectors, would simplify plant design, enable the use of standardized equipment and buildings, optimized for operational performance, and reduce plant size and weight. The need for site-specific equipment would be eliminated, enabling identical equipment to be used across multiple plants sited across the US and economies of scale, and catalyzing new interest and investment. The equipment-based protective systems would allow siting of advanced reactors in regions of high seismic hazard.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2023

Refining Principal Stress Measurements in Reservoir Underburden in Regions of Induced Seismicity through Seismological Tools, Laboratory Experiments - Final Technical Report

This project developed methodologies to measure the in-situ principal stress in the deep subsurface through use of multiple independent, but complementary, seismic methods, laboratory verification, and development of theoretical frameworks. By leveraging existing regional and local datasets we developed, tested, and refined a set of diagnostic tools for determining the in-situ stress state with reduced uncertainty at and below reservoir depths (1.5-6 km). A set of novel tools was produced that are scale independent, such that their utility is equivalent on regional, field scale, and near borehole monitoring of principal stresses in reservoir underburden for carbon storage projects. During a 4-year Department of Energy (DOE) and Southern Company funded project, carried out by the Electric Power Research Institute (EPRI), Lawrence Livermore National Laboratory (LLNL), the University of Oklahoma (OU), and the U.S. Geological Survey (USGS), the project team developed methodologies to measure the far-field in-situ principal stress in the deep subsurface, leveraging induced seismicity data from waste-water disposal projects. These methodologies consisted in the use of well-established and technically advanced seismic processing methods, such as virtual seismometer method-moment tensor (VSM-MT) and shear wave splitting (SWS), that are adept at recovering the stress orientation and certain components of the stress tensor. These methods were applied to robust seismicity catalogs created with matched filter techniques near sites of active fluid disposal—a proxy for carbon storage sites where such datasets are more limited. Estimates of the stress orientation made with seismic processing tools were considered along with laboratory acoustic emission experiments conducted on rock samples from the region of interest. Stress orientations in the studied region do not vary significantly across distances of ~100 km, nor are they found to rotate through time as a consequence of local wastewater disposal, as previously speculated. Finally, the project team investigated the trade-offs among the different seismic methods and evaluated the range of uncertainty that is generated with these methodologies, which led to a practical use and refinement of the VSM-MT technique when it is applied to field datasets. Understanding the trade-offs between these different methods highlighted the potential benefits of improved quantification of uncertainties on stress field estimations.

58 GEOSCIENCES↗
DOE OSTI2022

Porewater Geochemical Assessment of Seismic Indications for Gas Hydrate Presence and Absence: Mahia Slope, East of New Zealand’s North Island

We compare sediment vertical methane flux off the Mahia Peninsula, on the Hikurangi Margin, east of New Zealand’s North Island, with a combination of geochemical, multichannel seismic and sub-bottom profiler data. Stable carbon isotope data provided an overview of methane contributions to shallow sediment carbon pools. Methane varied considerably in concentration and vertical flux across stations in close proximities. At two Mahia transects, methane profiles correlated well with integrated seismic and TOPAS data for predicting vertical methane migration rates from deep to shallow sediment. However, at our “control site”, where no seismic blanking or indications of vertical gas migration were observed, geochemical data were similar to the two Mahia transect lines. This apparent mismatch between seismic and geochemistry data suggests a potential to underestimate gas hydrate volumes based on standard seismic data interpretations. To accurately assess global gas hydrate deposits, multiple approaches for initial assessment, e.g., seismic data interpretation, heatflow profiling and controlled-source electromagnetics, should be compared to geochemical sediment and porewater profiles. A more thorough data matrix will provide better accuracy in gas hydrate volume for modeling climate change and potential available energy content.

03 NATURAL GAS↗
DOE OSTI2023

Seismic Characterization of the Blue Mountain Geothermal Field

Subsurface characterization is crucial for geothermal energy exploration and production. Yet hydrothermal reservoirs usually reside in highly fractured and faulted zones where accurate characterization is very challenging because of low signal-to-noise ratios of land seismic data and lack of coherent reflection signals. We perform an active-source seismic characterization for the Blue Mountain geothermal field in Nevada using active seismic data to reveal the elastic medium property complexity and fault distribution at this field. We first employ an unsupervised machine learning method to attenuate groundroll and near-surface guided-wave noise and enhance coherent reflection and scattering signals from noisy seismic data. We then build a smooth initial P-wave velocity model based on an existing magnetotellurics survey result, and use 3D first-arrival traveltime tomography to refine the initial velocity model. We then derive a set of elastic wave velocities and anisotropic parameters using elastic full-waveform inversion, and obtain PP and PS images using elastic reverse-time migration. We identify major faults by analyzing the variations of seismic velocities and anisotropy parameters, and reveal mid- to small-scale faults by applying a supervised machine learning method to the seismic migration images. Our characterization reveals complex velocity heterogeneities and anisotropies, as well as faults, with a high spatial resolution. These results can provide valuable information for optimal placement of future injection and production wells to increase geothermal energy production at the Blue Mountain geothermal power plant.

58 GEOSCIENCES↗
DOE OSTI2020

EGS Collab Experiment 1: 3D Seismic Velocity Model and Updated Microseismic Catalog Using Transfer-Learning Aided Double-Difference Tomography

This package contains a 3D Seismic velocity model and an updated microseismic catalog associated with a proceedings paper (Chai et al., 2020) published in the 45th Workshop on Geothermal Reservoir Engineering. The 3D_seismic_velocity_model text file contains x (m), y(m), z(m), P-wave velocity (km/s), P-wave velocity quality indicator (1 for well-constrained; 0 for poorly constrained), S-wave velocity (km/s), and S-wave velocity quality indicator (1 for well-constrained; 0 for poorly constrained). The Updated_MEQ_catalog text file contains event origin time, x(m), y(m), z(m), error in x (m), error in y (m), error in z (m), and RMS misfit (millisecond). The 3D_seismic_P-wave_velocity_model animation file shows slices of the 3D P-wave velocity model. The 3D_seismic_S-wave_velocity_model animation file shows slices of the 3D S-wave velocity model. The Interactive_MEQ_locations API file is an interactive visualization of the updated microseismic event locations. The visualization allows users to view the event locations by dragging, rotating, and zooming in. References: Chai, C., Maceira, M., Santos-Villalobos, H. J., Venkatakrishnan, S. V., Schoenball, M., and EGS Collab Team, 2020, Automatic Seismic Phase Picking Using Deep Learning for the EGS Collab Project, in PROCEEDINGS, 45th Workshop on Geothermal Reservoir Engineering, edited, Stanford University, Stanford, California, 45, 1266-1276.

15 GEOTHERMAL ENERGY↗
DOE OSTI2021

Seismic savanna: machine learning for classifying wildlife and behaviours using ground‐based vibration field recordings

Abstract We develop a machine learning approach to detect and discriminate elephants from other species, and to recognise important behaviours such as running and rumbling, based only on seismic data generated by the animals. We demonstrate our approach using data acquired in the Kenyan savanna, consisting of 8000 h seismic recordings and 250 k camera trap pictures. Our classifiers, different convolutional neural networks trained on seismograms and spectrograms, achieved 80%–90% balanced accuracy in detecting elephants up to 100 m away, and over 90% balanced accuracy in recognising running and rumbling behaviours from the seismic data. We release the dataset used in this study: SeisSavanna represents a unique collection of seismic signals with the associated wildlife species and behaviour. Our results suggest that seismic data offer substantial benefits for monitoring wildlife, and we propose to further develop our methods using dense arrays that could result in a seismic shift for wildlife monitoring.

Szenicer, Alexandre↗
DOE OSTI2020

Geophysical monitoring using active seismic techniques at the Citronelle Alabama CO 2 storage demonstration site

Between August 2012 and September 2014, about 114,000 metric tonnes of CO 2 was captured from the coal-fired Plant Barry Power Station at Bucks Alabama and injected into the Paluxy Formation above the oil pool in the southeast unit of the Citronelle Oilfield. Various monitoring methods were deployed at land surface and in project wells to measure system performance, comply with permit requirements and test new and innovative monitoring tools. The monitoring program relied heavily on active seismic methods for subsurface imaging of geologic structure and time-lapse seismic techniques to track the CO 2 migration in the injection interval. Both conventional geophone/hydrophone and fiber-optic based Distributed Acoustic Sensing (DAS) arrays were deployed and tested, allowing a side by side comparison of the equipment and techniques. Geophysical imaging of the subsurface was successful using DAS in the offset vertical seismic profile (OVSP) survey configuration. A high resolution OVSP image of the subsurface was obtained in 2014 with DAS, which exceeded project expectations in comparison to a lower resolution image obtained in 2012 using a conventional 80-level geophone array. A time-lapse image of the redistribution of CO 2 after injection ended in September 2014 was obtained with two DAS OVSP surveys from June 2014 and December 2015, thus successfully demonstrating its proof-of-concept. Unfortunately, a pre-injection baseline survey with DAS, which was in its initial stage of technology development in 2012, did not have sufficient quality for use, making it difficult to interpret the acquired DAS time-lapse difference. Additional research in this area has since demonstrated the utility of time-lapse DAS OVSP. DAS data were also acquired during a cross-well seismic survey conducted in 2014. Unfortunately, the DAS technique was not success in the cross-well survey configuration because the system noise level was too high in the crosswell frequency output range (100–1200 Hz) of the piezoelectric source (increasing by a factor of ten compared to VSP frequency band). Additionally, the cross-well geometry causes sub-horizontal (broadside) incidence on the vertical DAS fiber cable, which is known to be problematic. Current research is focused on improving the DAS cable response to broadside acoustic energy. Time-lapse seismic surveys using commercially available conventional arrays were also acquired. In contrast to the DAS acquired data, the cross-well seismic results obtained with the conventional array was highly successful and clearly showed the CO 2 remained in zone at the end of injection. Time-lapse differencing of the OSVP surveys acquired with the conventional arrays proved to be inconclusive. Finally, changes in wellbore conditions between surveys and unavoidable changes in equipment (the array used for the baseline survey was retired) affected data quality, making it difficult to interpret the OVSP results.

58 GEOSCIENCES↗
DOE OSTI2020

Using Ambient Noise Seismic Interferometry and Local and Teleseismic Earthquakes to Determine Crustal Thickness and Moho Structure of the Northwestern Gulf of Mexico Margin

Abstract The northwestern part of the Gulf of Mexico has undergone two episodes of continental rifting and collision and produced structural artifacts that are now buried under many kilometers of sediments, complicating investigations of the region. The deep sedimentary package precludes outcrops and points to a need for the application of seismic techniques, but low rates of seismicity in the region and sparse seismic monitoring limit the utility of traditional seismic methods. We therefore use diverse data to perform two‐dimensional seismic tomography across the dry land portion of the margin. Data are gleaned from teleseismic and regional earthquakes, postcritical SsPmp arrivals, and direct P wave energy identified with seismic interferometry that were recorded by a broadband, three‐component array and partially overlapping short‐period, vertical‐component array. The Pn and postcritical SsPmp phase help constrain the Moho discontinuity, which a previous receiver function study suggested was absent beneath the seaward portion of this transect. A high‐velocity body is observed in the crust at the same location as the Houston Magnetic Anomaly, possibly marking rocks from the Alleghenian continental assembly. The crust thins from NW to SE, indicating that extension occurred mostly to the south of the Ouachita orogeny. Our model indicates that the margin's sediment package reaches a maximum thickness of ~15 km at the coast and becomes unresolvably thin near the Llano Uplift.

58 GEOSCIENCES↗
DOE OSTI2020

Identifying Different Classes of Seismic Noise Signals Using Unsupervised Learning

Abstract Proper classification of nontectonic seismic signals is critical for detecting microearthquakes and developing an improved understanding of ongoing weak ground motions. We use unsupervised machine learning to label five classes of nonstationary seismic noise common in continuous waveforms. Temporal and spectral features describing the data are clustered to identify separable types of emergent and impulsive waveforms. The trained clustering model is used to classify every 1 s of continuous seismic records from a dense seismic array with 10–30 m station spacing. We show that dominate noise signals can be highly localized and vary on length scales of hundreds of meters. The methodology demonstrates the complexity of weak ground motions and improves the standard of analyzing seismic waveforms with a low signal‐to‐noise ratio. Application of this technique will improve the ability to detect genuine microseismic events in noisy environments where seismic sensors record earthquake‐like signals originating from nontectonic sources.

Johnson, Christopher W.↗
DOE OSTI2022

Modeling Seismic Recordings of High-Frequency Guided Infrasound on Mars

NASA's Interior Exploration using Seismic Investigations, Geodesy and Heat Transport (InSight) mission records several high-frequency (>0.5 Hz) dispersive seismic signals on Mars. These signals are due to the acoustic-to-seismic coupling of infrasound generated by the entry and impact of meteorites. This dispersion property is due to infrasound propagating in a structured atmosphere, and we refer to this dispersive infrasound as guided infrasound. We propose to model the propagation of guided infrasound and the seismic coupling to the ground analytically; we use a 1D layered atmosphere on a three-layer solid subsurface medium. The synthetic ground movements fit the observed dispersive seismic signals well and the fitting indicates that the regolith beneath InSight is about 40-m in thickness. We also examine and validate the previously-published subsurface models derived from InSight ambient seismic vibration data.

58 GEOSCIENCES↗
DOE OSTI2024

Probing the evolution of fault properties during the seismic cycle with deep learning

We use seismic waves that pass through the hypocentral region of the 2016 M6.5 Norcia earthquake together with Deep Learning (DL) to distinguish between foreshocks, aftershocks and time-to-failure (TTF). Binary and N-class models defined by TTF correctly identify seismograms in test with > 90% accuracy. We use raw seismic records as input to a 7 layer CNN model to perform the classification. Here we show that DL models successfully distinguish seismic waves pre/post mainshock in accord with lab and theoretical expectations of progressive changes in crack density prior to abrupt change at failure and gradual postseismic recovery. Performance is lower for band-pass filtered seismograms (below 10 Hz) suggesting that DL models learn from the evolution of subtle changes in elastic wave attenuation. Tests to verify that our results indeed provide a proxy for fault properties included DL models trained with the wrong mainshock time and those using seismic waves far from the Norcia mainshock; both show degraded performance. Our results demonstrate that DL models have the potential to track the evolution of fault zone properties during the seismic cycle. If this result is generalizable it could improve earthquake early warning and seismic hazard analysis.

58 GEOSCIENCES↗
DOE OSTI2021

Explosion Discrimination Using Seismic Gradiometry and Spectral Filtering of Data

Here, we present a new method to discriminate between earthquakes and buried explosions using observed seismic data. The method is different from previous seismic discrimination algorithms in two main ways. First, we use seismic spatial gradients, as well as the wave attributes estimated from them (referred to as gradiometric attributes), rather than the conventional three-component seismograms recorded on a distributed array. The primary advantage of this is that a gradiometer is only a fraction of a wavelength in aperture compared with a conventional seismic array or network. Second, we use the gradiometric attributes as input data into a machine learning algorithm. The resulting discrimination algorithm uses the norms of truncated principal components obtained from the gradiometric data to distinguish the two classes of seismic events. Using high-fidelity synthetic data, we show that the data and gradiometric attributes recorded by a single seismic gradiometer performs as well as a conventional distributed array at the event type discrimination task.

58 GEOSCIENCES↗
DOE OSTI2023

Relocation of the Seismicity of the Caucasus Region

Our objective is to improve the view of the seismicity in the Caucasus region using instrumental data between 1951 and 2019. To create a comprehensive catalog, we combine the bulletins of local agencies and the International Seismological Centre, and use an advanced single-event location algorithm, iLoc, to obtain better locations. We show that relocations with iLoc, using travel-time predictions from the 3D upper mantle velocity model, Regional Seismic Travel Time, improve the locations. Then, using the iLoc results as initial locations and the ground-truth events identified in the iLoc results as fix points, we apply Bayesloc, a multiple-event location algorithm, to simultaneously relocate the entire seismicity of the Caucasus region. Here we demonstrate that the simultaneous relocation of the seismicity with Bayesloc clarifies the location and geometry of major active structures accommodating ongoing convergence between the Arabian and Eurasian continents between the Black and Caspian Seas. Among our major findings is the confirmation of widespread seismicity in the mantle beneath the northern flank of the Greater Caucasus and central Caspian, resulting from north-dipping subduction of the Kura and South Caspian basins and the identification of a discrete band of crustal seismicity beneath the southern flank of the Greater Caucasus.

58 GEOSCIENCES↗