SEARCH · Engineering Papers

Results for “array seismology”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

SourceResearch connectionFrom yearThrough year

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

35 records · Page 2

All sources NASA NTRS DOE OSTI

NASA NTRS2020

Cluster Analysis of Thermal Icequakes Using the Seismometer to Investigate Ice and Ocean Structure (SIIOS): Implications for Ocean World Seismology

Ocean Worlds are of high interest to the planetary community due to the potential habitability of their subsurface oceans. Over the next few decades several missions will be sent to ocean worlds including the Europa Clipper, Dragonfly, and possibly a Europa lander. The Dragonfly and Europa lander missions will carry seismic payloads tasked with detecting and locating seismic sources. The Seismometer to Investigate Ice and Ocean Structure (SIIOS) is a NASA PSTAR funded project that investigates ocean world seismology using terrestrial analogs. The goals of the SIIOS experiment include quantitatively comparing flight-candidate seismometers to traditional instruments, comparing single-station approaches to a small-aperture array, and characterizing the local seismic environment of our field sites. Here we present an analysis of detected local events at our field sites at Gulkana Glacier in Alaska and in Northwest Greenland approximately 80 km North of Qaanaaq, Greenland. Both field sites passively recorded data for about two weeks. We deployed our experiment on Gulkana Glacier in September 2017 and in Greenland in June 2018. At Gulkana there was a nearby USGS weather station which recorded wind data. Temperature data was collected using the MERRA satellite. In Greenland we deployed our own weather station to collect temperature and wind data. Gulkana represents a noisier and more active environment. Temperatures fluctuated around 0°C, allowing for surface runoff to occur during the day. The glacier had several moulins, and during deployment we heard several rockfalls from nearby mountains. In addition to the local environment, Gulkana is located close to an active plate boundary (relative to Greenland). This meant that there were more regional events recorded over two weeks, than in Greenland. Greenland’s local environment was also quieter, and less active. Temperatures remained below freezing. The Greenland ice was much thicker than Gulkana (~850 m versus ~100 m) and our stations were above a subglacial lake. Both conditions can reduce event detections from basal motion. Lastly, we encased our Greenland array in an aluminum vault and buried it beneath the surface unlike our array in Gulkana where the instruments were at the surface and covered with plastic bins. The vault further insulated the array from thermal and atmospheric events.

Marusiak, A. G.↗

DOE OSTI2020

Urban Seismic Site Characterization by Fiber-Optic Seismology

ccurate ground motion prediction requires detailed site effect assessment, but in urban areas where such assessments are most important, geotechnical surveys are difficult to perform, limiting their availability. Distributed acoustic sensing (DAS) offers an appealing alternative by repurposing existing fiber-optic cables, normally employed for telecommunication, as an array of seismic sensors. We present a proof-of-concept demonstration by using DAS to produce high-resolution maps of the shallow subsurface with the Stanford DAS array, California. We describe new methods and their assumptions to assess H/V spectral ratio—a technique widely used to estimate the natural frequency of the soil—and to extract Rayleigh wave dispersion curves from ambient seismic field. These measurements are jointly inverted to provide models of shallow seismic velocities and sediment thicknesses above bedrock in central campus. The good agreement with an independent survey validates the methodology and demonstrates the power of DAS for microzonation.

Geochemistry & Geophysics↗

NASA NTRS2011

Imaging the Moon's Core with Seismology

Constraining the structure of the lunar core is necessary to improve our understanding of the present-day thermal structure of the interior and the history of a lunar dynamo, as well as the origin and thermal and compositional evolution of the Moon. We analyze Apollo deep moonquake seismograms using terrestrial array processing methods to search for the presence of reflected and converted energy from the lunar core. Although moonquake fault parameters are not constrained, we first explore a suite of theoretical focal spheres to verify that fault planes exist that can produce favorable core reflection amplitudes relative to direct up-going energy at the Apollo stations. Beginning with stacks of event seismograms from the known distribution of deep moonquake clusters, we apply a polarization filter to account for the effects of seismic scattering that (a) partitions energy away from expected components of ground motion, and (b) obscures all but the main P- and S-wave arrivals. The filtered traces are then shifted to the predicted arrival time of a core phase (e.g. PcP) and stacked to enhance subtle arrivals associated with the Moon s core. This combination of filtering and array processing is well suited for detecting deep lunar seismic reflections, since we do not expect scattered wave energy from near surface (or deeper) structure recorded at varying epicentral distances and stations from varying moonquakes at varying depths to stack coherently. Our results indicate the presence of a solid inner and fluid outer core, overlain by a partial-melt-containing boundary layer (Table 1). These layers are consistently observed among stacks from four classes of reflections: P-to-P, S-to-P, P-to-S, and S-to-S, and are consistent with current indirect geophysical estimates of core and deep mantle properties, including mass, moment of inertia, lunar laser ranging, and electromagnetic induction. Future refinements are expected following the successful launch of the GRAIL lunar orbiter and SELENE 2 lunar lander missions.

Weber, Renee C.↗

DOE OSTI2025

Manipulation of Geographic Information in Global Seismology

Geographic data, such as seismic event locations, station locations, etc., are generally given in geographic latitude Φ ’, longitude θ , and depth below sea level, ζ , using the WGS84 ellipsoid as a reference. In software systems that use this type of geographic data, it is necessary to manipulate the data mathematically in order to perform such tasks as finding the angular distance or azimuth from one point to another, to find an array of points along a great circle, to rotate a point about a pole of rotation, to move a point some angular distance in a specified direction, to find the intersections of two great circles or to find the intersections of a great circle and a small circle. In this paper, equations are presented that convert geographic locations first to geocentric coordinates and then to Earth-centered Cartesian coordinates where many mathematical manipulations can be performed conveniently and efficiently.

58 GEOSCIENCES↗

DOE OSTI2022

On the measurement of S diff splitting caused by lowermost mantle anisotropy

Seismic anisotropy has been detected at many depths of the Earth, including its upper layers, the lowermost mantle and the inner core. While upper mantle seismic anisotropy is relatively straightforward to resolve, lowermost mantle anisotropy has proven to be more complicated to measure. Due to their long, horizontal ray paths along the core–mantle boundary (CMB), S waves diffracted along the CMB (S diff ) are potentially strongly influenced by lowermost mantle anisotropy. S diff waves can be recorded over a large epicentral distance range and thus sample the lowermost mantle everywhere around the globe. S diff therefore represents a promising phase for studying lowermost mantle anisotropy; however, previous studies have pointed out some difficulties with the interpretation of differential SH diff –SV diff traveltimes in terms of seismic anisotropy. Here, we provide a new, comprehensive assessment of the usability of S diff waves to infer lowermost mantle anisotropy. Using both axisymmetric and fully 3-D global wavefield simulations, we show that there are cases in which S diff can reliably detect and characterize deep mantle anisotropy when measuring traditional splitting parameters (as opposed to differential traveltimes). First, we analyze isotropic effects on S diff polarizations, including the influence of realistic velocity structure (such as 3-D velocity heterogeneity and ultra-low velocity zones), the character of the lowermost mantle velocity gradient, mantle attenuation structure, and Earth’s Coriolis force. Secondly, we evaluate effects of seismic anisotropy in both the upper and the lowermost mantle on SH diff waves. In particular, we investigate how SH diff waves are split by seismic anisotropy in the upper mantle near the source and how this anisotropic signature propagates to the receiver for a variety of lowermost mantle models. We demonstrate that, in particular and predictable cases, anisotropy leads to S diff splitting that can be clearly distinguished from other waveform effects. Here, these results enable us to lay out a strategy for the analysis of S diff splitting due to anisotropy at the base of the mantle, which includes steps to help avoid potential pitfalls, with attention paid to the initial polarization of S diff and the influence of source-side anisotropy. We demonstrate our S diff splitting method using three earthquakes that occurred beneath the Celebes Sea, measured at many transportable array stations at a suitable epicentral distance. We resolve consistent and well-constrained S diff splitting parameters due to lowermost mantle anisotropy beneath the northeastern Pacific Ocean.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2021

Data Release Report for the Source Physics Experiment Phase II: Dry Alluvium Geology Experiments (DAG-1 through DAG-4), Nevada National Security Site

The Dry Alluvium Geology (DAG) project was Phase II of the Source Physics Experiment and consisted of a series of four chemical explosive tests conducted in the same source hole on the Nevada National Security Site. This hole is located at 37.1146°N and -116.0693°W, with a surface elevation of 1,285.2 meters (m) (4,216.5 feet [ft]) above sea level. The first test (DAG-1) was conducted on July 20, 2018, at 16:51:52.67838 Coordinated Universal Time (UTC). The explosive source for DAG-1 was nitromethane initiated by a small plastic-bonded explosive (PBX) charge, detonated at the depth of 385.0 m (1,263.2 ft) below ground surface. DAG-1 had a trinitrotoluene (TNT) equivalent yield of 0.908 metric tons (2,002 pounds [lbs]). DAG-2 was conducted on December 19, 2018, at 18:45:56.92115 UTC. This test was the largest in the series, with a TNT equivalent yield of 50.997 metric tons (112,429 lbs). The explosive source for DAG-2 was nitromethane initiated by a small PBX charge, detonated at the depth of 299.8 m (983.6 ft) below ground surface. DAG-3 was conducted on April 27, 2019, at 15:49:01.84183 UTC. The explosive source for this test was nitromethane initiated by a small PBX charge, detonated at the depth of 149.9 m (492.0 ft) below ground surface. DAG-3 had a TNT equivalent yield of 0.908 metric tons (2,002 lbs). The final DAG test (DAG-4) was conducted on June 22, 2019, at 21:06:19.87632 UTC. The explosive source for DAG-4 was nitromethane initiated by a small PBX charge, detonated at the depth of 51.6 m (169.3 ft) below ground surface. DAG-4 had a TNT equivalent yield of 10.357 metric tons (22,833 lbs). The four tests were recorded by an extensive set of instrumentation that included sensors both at near-field (less than 200 m) and far-field (200 m or greater) distances. The near-field instruments consisted of three-component (3C) accelerometers installed at various depths ranging from 51.6 to 385 m (169.3 to 1,263.1 ft) below ground surface in boreholes positioned around the source hole, and arrays of single-component and 3C accelerometers on the surface. The far-field network comprised a variety of seismic and acoustic sensors, including short-period geophones, broadband seismometers, and 3C accelerometers at distances of 200 m to 400 kilometers. In addition, the DAG-2, DAG-3, and DAG-4 explosions were recorded by a temporary array of 496 geophones arranged in a densely spaced grid pattern known as “Large N.” This report coincides with the release of these data for analysts and organizations that are not participants in this program. This report describes the four DAG tests and the various types of near-field, far-field, and other data that are available. Assembled data sets are accessible through: Incorporated Research Institutions for Seismology, Data Management Center 1408 NE 45th Street, Suite 201, Seattle, Washington 98105 USA. www.iris.washington.edu

58 GEOSCIENCES↗

NASA NTRS2003

Kepler Mission: A Mission to Find Earth-size Planets in the Habitable Zone

The Kepler Mission is a Discovery-class mission designed to continuously monitor the brightness of 100,000 solar-like stars to detect the transits of Earth-size and larger planets. It is a wide field of view photometer Schmidt-type telescope with an array of 42 CCDs. It has a 0.95 m aperture and 1.4 m primary and is designed to attain a photometric precision of 2 parts in 10(exp 5) for 12th magnitude solar-like stars for a 6 hr transit duration. It will continuously observe 100,000 main-sequence stars from 9th to 14th magnitude in the Cygnus constellation for a period of four years with a cadence of 4/hour. An additional 250 stars can be monitored at a cadence of l/minute to do astro-seismology of stars brighter than 11.5 mv. The photometer is scheduled to be launched into heliocentric orbit in 2007. When combined with ground-based spectrometric observations of these stars, the positions of the planets relative to the habitable zone can be found. The spectra of the stars are also used to determine the relationships between the characteristics of terrestrial planets and the characteristics of the stars they orbit. In particular, the association of planet size and occurrence frequency with stellar mass and metallicity will be investigated. Based on the results of the current Doppler-velocity discoveries, over a thousand giant planets will also be found. Information on the albedos and densities of those giants showing transits will be obtained. At the end of the four year mission, hundreds of Earth-size planets should be discovered in and near the HZ of their stars if such planets are common. A null result would imply that terrestrial planets in the HZ are very rare and that life might also be quite rare.

Borucki, W. J.↗

NASA NTRS2003

Kepler Mission: A Technical Overview

Borucki, W. J.↗

DOE OSTI2021

In-situ observation of pre-, co- and post-seismic shear slip at 1.5 km depth

Understanding the initiation and arrest of earthquakes is one of the long-standing challenges of seismology. Here we report on direct observations of borehole displacement by a meter-sized shear rupture induced by pressurization of metamorphic rock at 1.5 km depth. We observed the acceleration of sliding, followed by fast co-seismic slip and transient afterslip. Total displacements were about 7, 5.5 and 9.5 micrometers, respectively for the observed pre-slip, co-seismic slip and afterslip. The observed pre-slip lasted about 0.4 seconds. Co-seismic slip was recorded by the 1 kHz displacement recording and a 12-component array of 3-C accelerometers sampled at 100 kHz. The observed afterslip is consistent with analytical models of arrest in a velocity-strengthening region and subsequent stress relaxation. < p >The observed slip vector agrees with the activation of a bedding plane within the phyllite, which is corroborated by relocated seismic events that were observed during the later stages of the injection experiment.

Schoenball, Martin↗

NASA NTRS1994

Permanent GPS Geodetic Array in Southern California (PGGA) and GPS observations in Indonesia

The Permanent GPS Geodetic Array (PGGA) is a network of permanent monitoring GPS stations in southern California devoted to the continuous measurement of crustal deformation in near real-time. The PGGA plays a unique role in studies of the kinematics of crustal deformation and the earthquake cycle in southern California because it is also providing temporally dense geodetic measurements of crustal motion over periods of minutes to variations in regional crustal strain. As it expands and matures the PGGA will play an increasingly important role in the study of active tectonics of southern California by bridging the frequency range between seismology, observatory geodesy, paleoseismology, and geology. In Indonesia GPS data is used for measurement of a large scale crustal deformation, extending from north China to the Indonesian archipelago. Indonesia offers a tremendous laboratory to study some of the processes that build continents, and mountains are active there. We began GPS observations in August 1989 on mainland Sumatra and the Mentawai Islands to study the phenomena of oblique plate convergence. We have analyzed the Indonesian data in conjunction with data collected on Christmas and Cocos Islands and at Darwin, Australia, and with the triangulation data in Sumatra.

Bock, Yehuds↗

DOE OSTI2022

Near-surface seismic anisotropy in Antarctic glacial snow and ice revealed by high-frequency ambient noise

Ambient seismic recordings taken at broad locations across Ross Ice Shelf and a dense array near West Antarctic Ice Sheet (WAIS) Divide, Antarctica, show pervasive temporally variable resonance peaks associated with trapped seismic waves in near-surface firn layers. These resonance peaks feature splitting on the horizontal components, here interpreted as frequency-dependent anisotropy in the firn and underlying ice due to several overlapping mechanisms driven by ice flow. Frequency peak splitting magnitudes and fast/slow axes were systematically estimated at single stations using a novel algorithm and compared with good agreement with active source anisotropy measurements at WAIS Divide determined via active sources recorded on a 1 km circular array. The approach was further applied to the broad Ross Ice Shelf (RIS) array, where anisotropy axes were directly compared with visible surface features and ice shelf flow lines. The near-surface firn, depicted by anisotropy above 30 Hz, was shown to exhibit a novel plastic stretching mechanism of anisotropy, whereby the fast direction in snow aligns with accelerating ice shelf flow.

58 GEOSCIENCES↗

DOE OSTI2025

Bayesian event categorization matrix approach for explosion monitoring

Current efforts to correctly categorize natural events from suspected explosion sources with data that is collected by ground- or space-based sensors presents historical challenges that remain unaddressed by the Event Categorization Matrix (ECM) model. Smaller historical events (lower yield explosions) may have data available from fewer measurement techniques than are available today, and therefore, a historical event record can lack a complete set of discriminants. The covariance structures can also differ between such observations of event (source-type) categories. Both obstacles are problematic for the classic ECM model. Our work addresses this gap and presents a Bayesian update to the previous ECM model, termed the Bayesian Event Categorization Matrix model, which can be trained on partial observations and does not rely on a pooled covariance structure. We further augment the ECM model with Bayesian Decision Theory so that false negative or false positive rates of an event categorization can be reduced in an intuitive manner. To demonstrate improved categorization rates for the Bayesian Event Categorization Matrix model, we compare an array of Bayesian and classic models with multiple performance metrics using Monte Carlo experiments. We use both synthetic and real data. Our Bayesian models show consistent gains in overall accuracy and lower false negative rates relative to the classic ECM model. Here, we propose future avenues to improve Bayesian Event Categorization Matrix models’ decision making and predictive capability.

58 GEOSCIENCES↗

NASA NTRS2011

The Effect of Sonic Booms on Earthquake Warning Systems

Several aerospace companies are designing quiet supersonic business jets for service over the United States. These aircraft have the potential to increase the occurrence of mild sonic booms across the country. This leads to interest among earthquake warning (EQW) developers and the general seismological community in characterizing the effect of sonic booms on seismic sensors in the field, their potential impact on EQW systems, and means of discriminating their signatures from those of earthquakes. The SonicBREWS project (Sonic Boom Resistant Earthquake Warning Systems) is a collaborative effort between Seismic Warning Systems, Inc. (SWS) and NASA Dryden Flight Research Center. This project aims to evaluate the effects of sonic booms on EQW sensors. The study consists of exposing high-sample-rate (1000 sps) triaxial accelerometers to sonic booms with overpressures ranging from 10 to 600 Pa in the free field and the built environment. The accelerometers record the coupling of the sonic boom to the ground and surrounding structures, while microphones record the acoustic wave above ground near the sensor. Sonic booms are broadband signals with more high-frequency content than earthquakes. Even a 1000 sps accelerometer will produce a significantly aliased record. Thus the observed peak ground velocity is strongly dependent on the sampling rate, and increases as the sampling rate is reduced. At 1000 sps we observe ground velocities that exceed those of P-waves from M L 3 earthquakes at local distances, suggesting that sonic booms are not negligible for EQW applications. We present the results of several experiments conducted under SonicBREWS showing the effects of typical-case low amplitude sonic booms and worst-case high amplitude booms. We show the effects of various sensor placements and sensor array geometries. Finally, we suggest possible avenues for discriminating sonic booms from earthquakes for the purposes of EQW.

Wurman, Gilead↗

DOE OSTI2024

Localization of infrasonic sources via Bayesian back projection

SUMMARY A Bayesian framework is investigated for event-specific localization of infrasonic sources using back projection ray tracing. Direction-of-arrival information from array-based detection analysis is used to initialize a back projection ray path originating from the detecting array location and quantifying propagation characteristics from hypothetical source locations. The Fisher statistic, computed from the array’s beam coherence, is mapped into uncertainty in the launch angles of the ray path. Auxiliary parameters previously introduced for solving the Transport equation to compute geometric spreading along ray paths are used to map uncertainty in the ray launch angles into spatial and temporal uncertainties in the ray path. An atmospheric ensemble approach is applied to account for atmospheric uncertainty, and the relation between uncertainties in the atmospheric state and confidence in estimated localization are evaluated using several ensembles with specified variances. The method is evaluated using a synthetic event in the western United States constructed via forward propagation simulations as well as a single-station, multi-arrival detection from a surface explosion in the western United States. Localization results using this event-specific approach are more accurate and exhibit improved precision than existing Bayesian localization methods that leverage generalized, pre-computed propagation statistics.

58 GEOSCIENCES↗

DOE OSTI2026

Reply to “Comment on ‘The Reduced Detection Rate of Signals That Are Hidden by Earthquakes: Case Studies with Spotlight Detectors That Operate at Seismic Arrays,’ by Joshua D. Carmichael, Brent G. Delbridge, and Richard Alfaro-Diaz” by Paul G. Richards

We reply to a comment that Paul G. Richards (Richards, 2026) directed to an article by Carmichael et al. (2025). Richards expressed concern over two issues. The first issue relates to the size of the explosions discussed in the article. The second relates to the probability that earthquakes and explosions sourced near the same location coincide in time, by chance. We do not address the first issue because the original article does not assign significance to the size of the explosions from the point of view of an experimental party. We do respond to Richards’ comment on the second issue and claim that it can be explained through a clarification about conditional versus joint probabilities. Our present discussion, therefore, describes the difference between the conditional probability that a correlator fails, given that an explosion and earthquake coincidently occur in the same place, and the joint probability that a correlator fails while an earthquake and explosion also coincidentally occur in the same place. This latter, joint probability appears to be rare. We clarify that our original article did not treat the joint probability, only the conditional probability.

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗

DOE OSTI2023

The influence of physical and algorithmic factors on simulated far-field waveforms and source–time functions of underground explosions using unsupervised machine learning

SUMMARY Characterizing explosion sources and differentiating between earthquake and underground explosions using distributed seismic networks becomes non-trivial when explosions are detonated in cavities or heterogeneous ground material. Moreover, there is little understanding of how changes in subsurface physical properties affect the far-field waveforms we record and use to infer information about the source. Simulations of underground explosions and the resultant ground motions can be a powerful tool to systematically explore how different subsurface properties affect far-field waveform features, but there are added variables that arise from how we choose to model the explosions that can confound interpretation. To assess how both subsurface properties and algorithmic choices affect the seismic wavefield and the estimated source functions, we ran a series of 2-D axisymmetric non-linear numerical explosion experiments and wave propagation simulations that explore a wide array of parameters. We then inverted the synthetic far-field waveform data using a linear inversion scheme to estimate source–time functions (STFs) for each simulation case. We applied principal component analysis (PCA), an unsupervised machine learning method, to both the far-field waveforms and STFs to identify the most important factors that control variance in the waveform data and differences between cases. For the far-field waveforms, the largest variance occurs in the shallower radial receiver channels in the 0–50 Hz frequency band. For the STFs, both peak amplitude and rise times across different frequencies contribute to the variance. We find that the ground equation of state (i.e. lithology and rheology) and the explosion emplacement conditions (i.e. tamped versus cavity) have the greatest effect on the variance of the far-field waveforms and STFs, with the ground yield strength and fracture pressure being secondary factors. Differences in the PCA results between the far-field waveforms and STFs could possibly be due to near-field non-linearities of the source that are not accounted for in the estimation of STFs and could be associated with yield strength, fracture pressure, cavity radius and cavity shape parameters. Other algorithmic parameters are found to be less important and cause less variance in both the far-field waveforms and STFs, meaning algorithmic choices in how we model explosions are less important, which is encouraging for the further use of explosion simulations to study how physical Earth properties affect seismic waveform features and estimated STFs.

58 GEOSCIENCES↗

DOE OSTI2025

Observations and seismoacoustic simulations of earthquake-generated infrasound waves in non-epicentral regions

We analysed infrasound waves associated with the Gyeongju earthquake (ML 5.8) that occurred on 2016 September 12, in the southeastern Korean Peninsula. For infrasound wave detection, the Progressive Multichannel Correlation method was applied to the infrasound data set recorded at seven arrays operating in South Korea at epicentral distances ranging from 178 to 472 km. Based on the back-projection method constrained by array-dependent celerity and azimuth deviation models, the source regions were identified in both the epicentral and non-epicentral regions. Remarkably, the non-epicentral secondary sources of this earthquake were located in regions with shallow water depths: (i) the western coastal area in the Yellow Sea and (ii) the shallow ocean basin and bank in the East Sea. The location results obtained from the earthquake could be corroborated through its foreshock (ML 5.1), yielding location results consistent with those of the main shock. The generation of infrasound waves over shallow water depths was fortuitously validated by direct recordings of dominant single-frequency (~0.3 Hz) infrasound waves at close range via temporary sensors near the ocean basin and bank. We interpreted that low-frequency infrasound signals could be generated from interactions among the ocean floor, shallow seawater and atmosphere. We performed numerical simulations of seismoacoustic fields to predict ground motions on the seafloor and acoustic transmission efficiency between the water and air interface. The simulations quantified the energy transfer through different media and clarified our observational results. We found that because this solid Earth‒water‒atmosphere coupled air wave has a relatively low frequency (~0.3 Hz), it can survive propagation over long distances compared with high-frequency infrasound waves generated in inland and mountain regions. In this study, we extend our understanding of water‒atmosphere coupling and the monitoring framework for earthquake-associated non-epicentral infrasound waves, encompassing not only inland ground shaking but also shallow sea regions located far from the epicentre.

Earthquake monitoring and test-ban treaty verifica↗