Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “seismic wave”

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

The Observed Inefficiency of Explosions to Produce Large Aftershocks: Båth’s Law for Explosions is 2.5

Underground explosions are observed to produce fewer and smaller aftershocks than similar size earthquakes. The seismic magnitude difference $Δm_x$ between an explosion and its largest aftershock is an expression of Båth’s law for explosions. Based on an analysis of a compilation of aftershock studies from Soviet testing at the Semipalatinsk test site in Kazakhstan and observations from American testing at the Nevada National Security Site (NNSS), here we find that the average magnitude difference for explosions $\overline{Δm_x}$ is about 2.5. Based on the NNSS data, two standard deviations of $Δm_x$ is about 1.5. In all the cases studied, from ton to megaton yield, from shallow to overburied depth, and chemical or nuclear source, no explosion aftershock has been larger than the explosion that preceded it. In fact, the two events at the NNSS with the largest aftershock magnitudes relative to the explosion are associated with the collapse of the cavity created by the explosion. This is similar to observations from North Korean testing at the Punggye-ri Test Site, where the largest seismic event following the test is attributed to the collapse after the 2017 explosion and is from 0.8 to 2 magnitude units less than the mainshock.

58 GEOSCIENCES↗
DOE OSTI2021

Unusual Seismic Signals in the Sevier Desert, Utah Possibly Related to the Black Rock Volcanic Field

Abstract On September 12, 2018 and April 14, 2019 shallow (<2.5 km from the surface, <1 km relative to sea level) earthquakes ( M w 4.0 and 4.1, respectively) struck a relatively aseismic area in south‐central Utah, within the Black Rock Desert volcanic field. Seismic waveforms are dominated by long‐period energy, no clear direct S‐wave arrivals, and signal durations of ∼100 s. Using co‐seismic interferograms, we detect clear and distinguishable surface deformation confirming the shallowness of the two mainshocks, while full moment tensor solutions rule out an isotropic source, such as an explosion or collapse. For the April 2019 sequence, we build an earthquake catalog using a dense temporary geophone array. The enhanced earthquake catalog shows a rapidly decaying aftershock sequence, with 85% of aftershocks occurring during the first three hours. Combining these observations, we argue that these earthquake sequences are related to the Black Rock Desert volcanic field.

Mesimeri, Maria↗
NASA NTRS1975

Lunar shear velocity structure at Apollo sites 12, 14, and 15

Spectral amplitude ratios of horizontal-to-vertical motion produced on seismograms of meteoroid impacts at the Apollo 12, 14, and 15 sites, the Apollo 14 and 15 lunar-module impacts, and the Apollo 15 S4B impact show consistent differences among the recording sites. On the assumption that the motion represents predominantly fundamental-mode Reyleigh waves and that the compressional-wave velocity structure is similar to that derived in other investigations involving the Apollo seismic experiments, estimates are made of the shear-wave velocity structure under the three sites. Near-surface velocities are about 35 m/s at the three sites. The results for site 14 indicate an increase to about 100 m/s near the 8-m depth and to 200 m/s at the 38-m depth. Results for sites 12 and 15 show a smoother gradient and generally a greater velocity at a given depth than that indicated at site 14, reaching velocities of about 400 m/s near the 120-m depth. If the assumed P velocity structures are correct and if changes in S velocity coincide with changes in P velocity, then the ratio of these velocities decreases from 2.9 to 2.0 in the upper 19 meters at site 12, in the upper 38 meters at site 14, and in the upper 21 meters at site 15.

Mark, N.↗
DOE OSTI2024

Evaluating Physics-Informed Neural Network Performance for Seismic Discrimination between Earthquakes and Explosions

In this article, we evaluate adding a weak physics constraint, that is, a physics‐based empirical relationship, to the loss function with a physics‐informed manner in local distance explosion discrimination in the hope of improving the generalization capability of the machine learning (ML) model. We compare the proposed model with the two‐branch model we previously developed, as well as with a pure data‐driven model. Unexpectedly, the proposed model did not consistently outperform the pure data‐driven model. By varying the level of inconsistency in the training data, we find this approach is modulated by the strength of the physics relationship. In conclusion, this result has important implications for how to best incorporate physical constraints in ML models.

58 GEOSCIENCES↗
NASA NTRS2011

Towards Simulating a Realistic Planetary Seismic Wavefield: The Contribution of the Megaregolith and Low-Velocity Waveguides

Lunar seismograms are distinctly different from their terrestrial counterparts. The Apollo lunar seismometers recorded moonquakes without distinct P- or S-wave arrivals; instead waves arrive as a diffuse coda that decays over several hours making the identification of body waves difficult. The unusual character of the lunar seismic wavefield is generally tied to properties of the megaregolith: it consists of highly fractured and broken crustal rock, the result of extensive bombardment of the Moon. The megaregolith extends several kilometers into the lunar crust, possibly into the mantle in some regions, and is covered by a thin coating of fine-scale dust. These materials possess very low seismic velocities that strongly scatter the seismic wavefield at high frequencies. Directly modeling the effects of the megaregolith to simulate an accurate lunar seismic wavefield is a challenging computational problem, owing to the inherent 3-D nature of the problem and the high frequencies (greater than 1 Hz) required. Here we focus on modeling the long duration code, studying the effects of the low velocities found in the megaregolith. We produce synthetic seismograms using 1-D slowness integration methodologies, GEMINI and reflectivity, and a 3-D Cartesian finite difference code, Wave Propagation Program, to study the effect of thin layers of low velocity on the surface of a planet. These codes allow us generate seismograms with dominant frequencies of approximately 1 Hz. For background lunar seismic structure we explore several models, including the recent model of Weber et al., Science, 2011. We also investigate variations in megaregolithic thickness, velocity, attenuation, and seismogram frequency content. Our results are compared to the Apollo seismic dataset, using both a cross correlation technique and integrated envelope approach to investigate coda decay. We find our new high frequency results strongly support the hypothesis that the long duration of the lunar seismic codes is generated by the presence of the low velocity megaregolith, and that the diffuse arrivals are a combination of scattered energy and multiple reverberations within this layer. The 3-D modeling indicates the extreme surface topography of the Moon adds only a small contribution to scattering effects, though local geology may play a larger role. We also study the effects of the megaregolith on core reflected and converted phases and other body waves. Our analysis indicates detection of core interacting arrivals with a polarization filter technique is robust and lends the possibility of detecting other body waves from the Moon.

Schmerr, Nicholas C.↗
DOE OSTI2022

Development, verification, and validation of comprehensive acoustic fluid-structure interaction capabilities in an open-source computational platform

The acoustic fluid-structure interaction (FSI) formulation is a practical numerical approach for the seismic analysis of fluid-filled tanks. However, there are no verification and validation studies reported in the literature that demonstrate the ability of an acoustic FSI numerical model to predict responses important to structural and mechanical design for intense translational and rotational earthquake inputs. Herein, an acoustic FSI formulation is implemented in the open-source Multiphysics Object-Oriented Simulation Environment (MOOSE), and is formally verified and validated using analytical solutions and code-to-code verification, and experimental data, respectively. The analytical solutions are for small amplitude, unidirectional seismic inputs. The code-to-code verification utilizes a previously verified and validated Arbitrary Lagrangian-Eulerian (ALE) numerical model in the commercial finite element code LS-DYNA. The validation studies utilize a comprehensive data set assembled from results of 3D earthquake-simulator tests of a fluid-filled vessel. The acoustic numerical model in MOOSE is verified and validated for hydrodynamic pressures and support reactions except for cases that involve significant convective response. For small amplitude inputs, numerically predicted wave heights match those of the analytical solutions. The numerical model is not verified and validated for wave height calculations under intense 3D seismic inputs. The run times for the acoustic FSI simulations in MOOSE are an order of magnitude, or more, shorter than for the corresponding ALE simulations in LS-DYNA. The utility of the MOOSE acoustic FSI implementation is demonstrated by seismic analysis of a building equipped with a fluid-filled, advanced nuclear reactor.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2021

PRISM3D: a 3-D reference seismic model for Iberia and adjacent areas

SUMMARY We present PRISM3D, a 3-D reference seismic model of P- and S-wave velocities for Iberia and adjacent areas. PRISM3D results from the combination of the most up-to-date earth models available for the region. It extends horizontally from 15°W to 5°E in longitude, 34°N to 46°N in latitude and vertically from 3.5 km above to 200 km below sea level, and is modelled on a regular grid with 10 and 0.5 km of grid node spacing in the horizontal and vertical directions, respectively. It was designed using models inferred from local and teleseismic body-wave tomography, earthquake and ambient noise surface wave tomography, receiver function analysis and active source experiments. It includes two interfaces, namely the topography/bathymetry and the Mohorovičić (Moho) discontinuity. The Moho was modelled from previously published receiver function analysis and deep seismic sounding results. To that end we used a probabilistic surface reconstruction algorithm that allowed to extract the mean of the Moho depth surface along with its associated standard deviation, which provides a depth uncertainty estimate. The Moho depth model is in good agreement with previously published models, although it presents slightly sharper gardients in orogenic areas such as the Pyrenees or the Betic-Rif system. Crustal and mantle P- and S-wave wave speed grids were built separately on each side of the Moho depth surface by weighted average of existing models, thus allowing to realistically render the speed gradients across that interface. The associated weighted standard deviation was also calculated, which provides an uncertainty estimation on the average wave speed values at any point of the grid. At shallow depths (<10 km), low P and S wave speeds and high VP/VS are observed in offshore basins, while the Iberian Massif, which covers a large part of western Iberia, appears characterized by a rather flat Moho, higher than average VP and VS and low VP/VS. Conversely, the Betic-Rif system seems to be associated with low VP and VS, combined with high VP/VS in comparison to the rest of the study area. The most prominent feature of the mantle is the well known high wave speed anomaly related to the Alboran slab imaged in various mantle tomography studies. The consistency of PRISM3D with previous work is verified by comparing it with two recent studies, with which it shows a good general agreement.The impact of the new 3-D model is illustrated through a simple synthetic experiment, which shows that the lateral variations of the wave speed can produce traveltime differences ranging from –1.5 and 1.5 s for P waves and from –2.5 and 2.5 s for S waves at local to regional distances. Such values are far larger than phase picking uncertainties and would likely affect earthquake hypocentral parameter estimations. The new 3-D model thus provides a basis for regional studies including earthquake source studies, Earth structure investigations and geodynamic modelling of Iberia and its surroundings.

Arroucau, P.↗
DOE OSTI2022

Seeking Repeating Anthropogenic Seismic Sources: Implications for Seismic Velocity Monitoring at Fault Zones

Abstract Seismic velocities in rocks are highly sensitive to changes in permanent deformation and fluid content. The temporal variation of seismic velocity during the preparation phase of earthquakes has been well documented in laboratories but rarely observed in nature. It has been recently found that some anthropogenic, high‐frequency (>1 Hz) seismic sources are powerful enough to generate body waves that travel down to a few kilometers and can be used to monitor fault zones at seismogenic depth. Anthropogenic seismic sources typically have fixed spatial distribution and provide new perspectives for velocity monitoring. In this work, we propose a systematic workflow to seek such powerful seismic sources in a rapid and straightforward manner. We tackle the problem from a statistical point of view, considering that persistent, powerful seismic sources yield highly coherent correlation functions (CFs) between pairs of seismic sensors. The algorithm is tested in California and Japan. Multiple sites close to fault zones show high‐frequency CFs stable for an extended period of time. These findings have great potential for monitoring fault zones, including the San Jacinto Fault and the Ridgecrest area in Southern California, Napa in Northern California, and faults in central Japan. However, extra steps, such as beamforming or polarization analysis, are required to determine the dominant seismic sources and study the source characteristics, which are crucial to interpreting the velocity monitoring results. Train tremors identified by the present approach have been successfully used for seismic velocity monitoring of the San Jacinto Fault in previous studies.

58 GEOSCIENCES↗
DOE OSTI2021

3D Seismic Anatomy of a Watershed Reveals Climate–Topography Coupling That Drives Water Flowpaths and Bedrock Weathering

To investigate how bedrock transforms to soil, we mapped the topography of the interface demarcating onset of weathering under an east-west trending shale watershed in the Valley and Ridge province in the U.S.A. Using wave-equation travel-time tomography from a seismic array of >4000 geophones, we obtained a 3D P-wave velocity (Vp) model that resolves structures ~20 m below land surface (mbls). The depth of mobile soil and the onset of dissolution of chlorite roughly match Vp = 600 m/s and Vp = 2700 m/s, respectively. Chlorite dissolution initiates porosity growth in the shale matrix. Depth to the 2700 m/s contour is greater under the N- as compared to S-facing hillslopes and under sub-planar as compared to concave-up land surfaces. Broadly, the geometries of the ‘soil’ and ‘chlorite’ Vp contours are consistent with the calculated potential for shear fracture opening under weak regional compression. However, this calculated fracture potential does not consistently explain observations related to N- vs. S-facing aspect nor fracture density observed by borehole televiewer. Apparently, regional compression is only a secondary influence on Vp: the primary driver of P-wave slowing in the upper layers of this catchment is topographic control of reactive water flowpaths and their integrated effects on weathering. Here, the Vp result is best explained as the long-term integrated effect of groundwater flow-induced geochemical weathering of shale in response to climate-driven patterns of micro- and macro-topography.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Fast Seismic Anomalies Under Continents Explained by the Delaminated Lower Continental Crust—Implications From High Pressure‐Temperature Elasticity of Jadeite

Seismic tomography has shown that the shear wave velocities (Vs) under continents, especially under cratons, are extremely fast at 100–200 km depth, which is difficult to explain by low temperatures or high Mg#. Alternatively, delaminated eclogitic lower continental crust has been proposed to account for these fast seismic anomalies. However, the thermoelastic properties of jadeite which constitutes up to 60–80 mol% of clinopyroxene in the potentially delaminated lower continental crust are not well constrained. In this study, we measured the single-crystal elasticity of jadeite by Brillouin spectroscopy under simultaneous high pressure and temperature conditions for the first time. We found that the temperature dependence of Vs of jadeite is extremely small if not negligible. The seismic velocities of the potentially delaminated lower continental crusts were subsequently modeled and found to match the widely observed fast seismic anomalies under cratons between 100 and 200 km depth.

Geology↗
DOE OSTI2023

Crustal and uppermost mantle structure near the Gloria Fault, North Atlantic, from ocean bottom seismometer surface wave observations

SUMMARY In this work, we present both 1-D and 3-D shear wave velocity (Vs) models of the oceanic crust and uppermost mantle below the Deep OCean Test ARray area, located ∼ 70 km north of the central section of the Gloria Fault, in the eastern North-Atlantic Ocean. The velocity models are inferred from the dispersion of surface waves recorded on ocean bottom seismometers. Dispersion measurements are obtained from the analysis of ambient seismic noise at short periods (< 14 s) and teleseismic surface waves at long periods (> 14 s) using the two-station method. The 1-D Vs model is inferred from the joint inversion of Rayleigh wave phase velocities and Love wave group and phase velocities. The 3-D tomographic model is obtained by inversion of 2-D Love wave group velocity maps as a function of depth, further constrained by the average of Love wave phase velocities obtained from ambient noise (4–9 s) and the average Rayleigh and Love wave phase velocities calculated from teleseismic data (14–44 s). The 1-D Vs model shows a sediment layer with a low velocity of 1.05 km s−1, similar to previous studies in the region. Below the sediments, we find an oceanic crust with velocities ranging from 3.3 to 4.5 km s−1. The model reaches an unusually high velocity of 4.9 km s−1 in a 20 km thick layer at depths between 16 and 36 km. We interpret this fast velocity layer as indicative of the presence of harzburgite, a residue of enhanced melting that might have been formed by the proximity between the Mid-Atlantic Ridge and the Azores mantle plume. At greater depths the velocity decreases, forming a low-velocity zone that reaches a minimum at ∼ 70 km depth, which we interpret as the maximum depth for the lithosphere–asthenosphere boundary. The 3-D model shows a structure that is mostly horizontally layered, with Vs isocontours at 3.5–4.5 km s−1 highlighting oscillations of the crustal structure with wavelengths of ∼25–30 km. These oscillations may be due to changes in the rate of mantle upwelling and magma supply rate.

Pinzón, Juan I. (ORCID:0000000338078068)↗
DOE OSTI2024

Numerical modelling of impact seismic sources using the stress glut theory

SUMMARY Meteorite impacts have proved to be a significant source of seismic signal on the Moon, and have now been recorded on Mars by InSight seismometers. Understanding how impacts produce seismic signal is key to the interpretation of this unique data, and to improve their identification in continuous seismic records. Here, we use the seismic Representation Theorem, and particularly the stress glut theory, to model the seismic motion resulting from impact cratering. The source is described by equivalent forces, some resulting from the impactor momentum transfer, and others from the stress glut, which represents the mechanical effect of plasticity and non linear processes in the source region. We condense these equivalent forces into a point-source with a time-varying single force and nine-component moment tensor. This analytical representation bridges the gap between the complex dynamics of crater formation, and the linear point-source representation classically used in seismology. Using the multiphysics modelling software HOSS, we develop a method to compute the stress glut of an impact, and the associated point-source from hypervelocity impact simulations. For a vertical and an oblique impact at 1000 m s−1, we show that the moment tensor presents a significant deviatoric component. Hence, the source is not an ideal isotropic explosion contrary to previous assumptions, and draws closer to a double couple for the oblique impact. The contribution of the point force to the seismic signal appears negligible. We verify this model by comparing two signals: (1) HOSS is coupled to SPECFEM3D to propagate the near-source signal elastically to remote seismic stations; (2) the point-source model derived from the stress-glut theory is used to generate displacements at the same distance. The comparison shows that the point-source model is accurately simulating the low-frequency impact seismic waveform, and its seismic moment is in trend with Lunar and Martian impact data. High-frequencies discrepancies exist, which are partly related to finite-source effects, but might be further explained by the difference in mathematical framework between classical seismology and HOSS’ numerical modelling.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Seismic Imaging of a Shale Landscape Under Compression Shows Limited Influence of Topography–Induced Fracturing

We used seismic refraction to image the P-wave velocity structure of a shale watershed experiencing regional compression in the Valley and Ridge Province (USA). From estimates showing strong compressional stress, we expected the depth to unweathered bedrock to mirror the hill-valley-hill topography (“bowtie pattern”) by analogy to seismic velocity patterns in crystalline bedrock in the North American Piedmont that also experience compression. Previous researchers used failure potentials calculated for strong compression in the Piedmont to suggest fractures are open deeper under hills than valleys to explain the “bowtie” pattern. Seismic images of the shale watershed, however, show little evidence of such a “bowtie.” Instead, they are consistent with weak (not strong) compression. Here, this contradiction could be explained by the greater importance of infiltration-driven weathering than fracturing in determining seismic velocities in shale compared to crystalline bedrock, or to local perturbations of the regional stress field due to lithology or structures.

58 GEOSCIENCES↗
DOE OSTI2023

Seismoacoustic Analysis of the Large Surface Explosion Coupling Experiment Using a Large- N Seismic Array

Seismoacoustic signals at local distance (<~10 km) are widely used as important constraints on source parameters for near-surface events, yet the seismoacoustic wave generation and energy partitioning are not fully understood. Spatially dense sensors could provide observations in high resolution to capture the full wavefield for better understanding wave propagation and improving source estimation. Recently, spatially dense observations of the local seismoacoustic wavefield produced by a pair of 1-ton surface explosions have been recorded using a large-N seismic array. This large-N array consists of 446 geophones and covers an area of approximately 2×2.5 km 2 . The two surface explosions occur at the same location but at different times with different atmospheric conditions. Both seismic and air–ground coupled acoustic waves from the two surface explosions are well observed. Analyses of signals recorded by the large-N seismic array show different acoustic wave speed and amplitude for the two explosions. A strong spatial variability in acoustic wave speed and amplitude for each explosion is also observed. Finally, the observations suggest the important role of local atmosphere state on wave propagation and source estimation and demonstrate how the use of a large-N capability can improve characterization of the propagation medium and source.

58 GEOSCIENCES↗
DOE OSTI2022

Atmospheric waves and global seismoacoustic observations of the January 2022 Hunga eruption, Tonga

The 15 January 2022 climactic eruption of Hunga volcano, Tonga, produced an explosion in the atmosphere of a size that has not been documented in the modern geophysical record. The event generated a broad range of atmospheric waves observed globally by various ground-based and spaceborne instrumentation networks. Most prominent was the surface-guided Lamb wave (≲0.01 hertz), which we observed propagating for four (plus three antipodal) passages around Earth over 6 days. As measured by the Lamb wave amplitudes, the climactic Hunga explosion was comparable in size to that of the 1883 Krakatau eruption. The Hunga eruption produced remarkable globally detected infrasound (0.01 to 20 hertz), long-range (~10,000 kilometers) audible sound, and ionospheric perturbations. Seismometers worldwide recorded pure seismic and air-to-ground coupled waves. Air-to-sea coupling likely contributed to fast-arriving tsunamis. In this article, we highlight exceptional observations of the atmospheric waves.

58 GEOSCIENCES↗
NASA NTRS2022

Atmospheric Waves and Global Seismoacoustic Observations of the January 2022 Hunga Eruption, Tonga

The 15 January 2022 climactic eruption of Hunga volcano, Tonga, produced an explosion in the atmosphere of a size that has not been documented in the modern geophysical record. The event generated a broad range of atmospheric waves observed globally by various ground-based and spaceborne instrumentation networks. Most prominent was the surface-guided Lamb wave (≲0.01 hertz), which we observed propagating for four (plus three antipodal) passages around Earth over 6 days. As measured by the Lamb wave amplitudes, the climactic Hunga explosion was comparable in size to that of the 1883 Krakatau eruption. The Hunga eruption produced remarkable globally detected infrasound (0.01 to 20 hertz), long-range (~10,000 kilometers) audible sound, and ionospheric perturbations. Seismometers worldwide recorded pure seismic and air-to-ground coupled waves. Air-to-sea coupling likely contributed to fast-arriving tsunamis. Here, we highlight exceptional observations of the atmospheric waves.

Robin S Matoza↗
DOE OSTI2022

Wintertime Brine Discharge at the Surface of a Cold Polar Glacier and the Unexpected Absence of Associated Seismicity

A subglacial groundwater system beneath Taylor Glacier, Antarctica, discharges hypersaline, iron-rich brine episodically at the glacier surface to create Blood Falls. However, the triggering mechanism for these brine release events is not yet understood. Identifying which fracture processes are observed seismically can help us better characterize the hydrological system at Taylor Glacier, and more generally, provide us with a broader understanding of englacial hydrologic activity in cold glaciers. We document wintertime brine discharge using time-lapse photography. Subfreezing air temperatures during the brine discharge indicate that surface melt-induced hydrofracture is an unlikely trigger for brine release. Further, we analyze local seismic data to test a hypothesis that fracturing generates elevated surface wave energy preceding and/or coinciding with brine release events. Our results show no discernible elevated Rayleigh wave activity prior to or during Blood Falls brine release. Instead, we find a pattern of seismic events dominated by a seasonal signal, with more Rayleigh events occurring in the summer than the winter from the Blood Falls source area. We calculate that the volumetric opening of cracks that would generate Rayleigh waves at our detection limits are of similar size to myriad cracks in glacier ice, lake ice, and frozen sediment in the terminus area. We therefore propose that any fracturing coincident with brine release activity likely consists of a series of smaller opening events that are masked by other seismicity in the local environment.

58 GEOSCIENCES↗
DOE OSTI2021

Local-Distance Seismic Event Relocation and Relative Magnitude Estimation, Applications to Mining Related Seismicity in the Powder River Basin, Wyoming

Recent efforts to characterize small ( M w < 3 ) seismic events at local distances have become more important because of the increased observation of human-triggered and induced seismicity and the need to advance nuclear explosion monitoring capabilities. The signals generated by low-magnitude seismic sources necessitate the use of nearby short-period observations, which are sensitive to local geological heterogeneity. Local to near-regional distance ( < 300 km ) surface and shear waves can dominate short-period observations from small, shallow seismic sources. In this work, we utilize these observations to estimate precise, relative locations and magnitudes of ~ 700 industrial mining events in Wyoming, using nearly 360,000 observations. The precise, relative location estimates (with formal location uncertainty estimates of less than 1 km) collapse a diffuse collection of mining events into discrete clusters associated with individual blasting operations. We also invert the cross-correlation amplitudes to estimate precise, relative moment magnitude estimates, which help validate and identify disparities in the event sizes reported by regional network catalogs. Joint use of multiple phases allows for the inclusion of more seismic events due to the increase in the number of observations. In some cases, using a single phase allowed us to relocate only 50% of the original reported seismic events within a cluster. Combining shear- and surface-wave phases increased the number of events to above 90% of the original events, allowing us to characterize a broader range of event sizes, source to station distances, and event distributions. This analysis takes a step toward making a fuller characterization of small industrial seismic events observed at local distances.

58 GEOSCIENCES↗