USGS ScienceSearch

SEARCH · USGS Science

Results for “The Seismological Record”

Search indexed USGS publications on groundwater, aquifers, geologic maps, mineral resources and earthquakes. Explore source records by subject and place.

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

At least 163 records · Page 9Linked to original sources

Evaluation of intensity prediction equations (IPEs) for small-magnitude earthquakes

This study assesses existing intensity prediction equations (IPEs) for small unspecified magnitude ( M ≤3.5) earthquakes at short hypocentral distances ( ⁠ D h ⁠ ) and explores such earthquakes’ contribution to the felt shaking hazard. In particular, we consider IPEs by Atkinson and Wald (2007) and Atkinson et al. (2014) , and evaluate their performance based on “Did You Feel It” (DYFI) reports and recorded peak ground velocities (PGVs) in the central United States. Both IPEs were developed based on DYFI reports in the central and eastern United States with moment magnitudes above M w 3.0. DYFI reports are often used as the ground truth when evaluating and developing IPEs, but they could be less reliable when there are limited responses for small‐magnitude earthquakes. We first compare the DYFI reports with intensities interpolated from recorded PGVs. Results suggest a minimal discrepancy between the two when the intensity is large enough to be felt (i.e., M >2 and D h < 15 km ⁠ ). We then compare intensities from 31,617 DYFI reports of 3049 earthquakes with the two IPEs. Results suggest that both the IPEs match well with observed intensities for 2.0< M <3.0 and D h < 10 km ⁠ , but the IPE by Atkinson et al. (2014) matches better for larger distances. We also observe that intensities from DYFI reports attenuate faster compared with the two IPEs, especially for distances greater than 10 km. We then group DYFI reports by inferred V S 30 as a proxy for site amplification effects. We observe that intensities at sites with V S 30 around 300 m/s are consistently higher than at sites with V S 30 around 700 m/s and are also closer to the two IPEs. Finally, we conduct hazard disaggregation for earthquakes at close distances ( ⁠ D h = 7.5 km ⁠ ) using the observed records. Results suggest that earthquakes with magnitudes below M 3.0 contribute more than 40% to the occurrence of felt shaking.

Bulletin of the Seismological Society of America

The potential of using fiber optic distributed acoustic sensing (DAS) in earthquake early warning applications

As the seismological community embraces fiber optic distributed acoustic sensing (DAS), DAS arrays are becoming a logical, scalable option to obtain strain and ground‐motion data for which the installation of seismometers is not easy or cheap, such as in dense offshore arrays. The potential of strain data in earthquake early warning (EEW) applications has been recently demonstrated using records from borehole strainmeters (BSMs). However, current BSM networks are sparse, installing more BSMs is expensive and often impractical, and BSMs have the same limitations in offshore environments as other traditional seismic instruments. Here, we aim to provide a road map about how DAS data could be used in existing EEW applications, using the ShakeAlert EEW System for the West Coast of the United States as an example. We review the data requirements for EEW systems, examine ways in which strain‐derived ground‐motion data can be incorporated into such systems without significant modifications, and determine what is still needed for full utilization of DAS data in these applications. Importantly, EEW algorithms require ground‐motion amplitude information for rapid earthquake source characterization; thus, accurate strain amplitude observations, not only phase information, are necessary for deriving these ground‐motion metrics from DAS data. To obtain high‐quality ground‐motion observations, EEW‐compatible DAS arrays need to be multicomponent, well coupled, and low noise. We suggest ways to achieve such data requirements using existing DAS technology and discuss areas in which further research is needed to optimize DAS array performance for EEW.

Bulletin of the Seismological Society of America

Coal-mining seismicity and ground-shaking hazard: A case study in the Trail Mountain area, Emery County, Utah

We describe a multipart study to quantify the potential ground-shaking hazard to Joes Valley Dam, a 58-m-high earthfill dam, posed by mining-induced seismicity (MIS) from future underground coal mining, which could approach as close as ???1 km to the dam. To characterize future MIS close to the dam, we studied MIS located ???3-7 km from the dam at the Trail Mountain coal mine. A 12-station local seismic network (11 stations above ground, one below, combining eight triaxial accelerometers and varied velocity sensors) was operated in the Trail Mountain area from late 2000 through mid-2001 for the dual purpose of (1) continuously monitoring and locating MIS associated with longwall mining at a depth of 0.5-0.6 km and (2) recording high-quality data to develop ground-motion prediction equations for the shallow MIS. (Ground-motion attenuation relationships and moment-tensor results are reported in companion articles.) Utilizing a data set of 1913 earthquakes (M ??? 2.2), we describe space-time-magnitude distributions of the observed MIS and source-mechanism information. The MIS was highly correlated with mining activity both in space and time. Most of the better-located events have depths constrained within ??0.6 km of mine level. For the preponderance (98%) of the 1913 located events, only dilatational P-wave first motions were observed, consistent with other evidence for implosive or collapse-type mechanisms associated with coal mining in this region. We assess a probable maximum magnitude of M 3.9 (84th percentile of a cumulative distribution) for potential MIS close to Joes Valley Dam based on both the worldwide and regional record of coal-mining-related MIS and the local geology and future mining scenarios.

Bulletin of the Seismological Society of America

A preliminary summary of a seismic-refraction survey in the vicinity of the Tonto Forest Observatory, Arizona

The U.S. Geological Survey complete d a seismic-refraction survey in the vicinity of the Tonto Forest Seismological Observatory (T.F.S.O.) in April and May 1964. More than 1200 km of reversed profiles were surveyed to determine the crustal structure and crustal and upper mantle velocities in this area. The purpose of this work was to provide information on wave-propagation paths of seismic events recorded at T.F.S.O. and to improve the performance of the Observatory in locating and identifying these events. First arrivals indicate that the Mohorovicic discontinuity dips to the northeast by as much as 6 degrees under T.F.S.O., and may even be displaced vertically by as much as 5 km immediately north of the Observatory near the boundary of the Basin and Range and the Colorado Plateau Provinces. A preliminary examination of the first arrivals indicates that the crust at T.F.S.O. is at least 30 km thick and is made up of at least two seismic layers. A thin veneer at the surface with a velocity of approximately 4 km/sec is underlain by a layer with a velocity of approximately 5.9 km/sec to 6.1 km/sec. An intermediate layer with velocity of 6.6 to 7.0 km/sec is probably present in the lower crust, but is not revealed by first arrivals. The velocity of seismic waves in the upper mantle is about 7.9 km/sec.

Arizona

Improvements in absolute seismometer sensitivity calibration using local earth gravity measurements

The ability to determine both absolute and relative seismic amplitudes is fundamentally limited by the accuracy and precision with which scientists are able to calibrate seismometer sensitivities and characterize their response. Currently, across the Global Seismic Network (GSN), errors in midband sensitivity exceed 3% at the 95% confidence interval and are the least‐constrained response parameter in seismic recording systems. We explore a new methodology utilizing precise absolute Earth gravity measurements to determine the midband sensitivity of seismic instruments. We first determine the absolute sensitivity of Kinemetrics EpiSensor accelerometers to 0.06% at the 99% confidence interval by inverting them in a known gravity field at the Albuquerque Seismological Laboratory (ASL). After the accelerometer is calibrated, we install it in its normal configuration next to broadband seismometers and subject the sensors to identical ground motions to perform relative calibrations of the broadband sensors. Using this technique, we are able to determine the absolute midband sensitivity of the vertical components of Nanometrics Trillium Compact seismometers to within 0.11% and Streckeisen STS‐2 seismometers to within 0.14% at the 99% confidence interval. The technique enables absolute calibrations from first principles that are traceable to National Institute of Standards and Technology (NIST) measurements while providing nearly an order of magnitude more precision than step‐table calibrations.

Bulletin of the Seismological Society of America

Introduction to the special issue on rotational seismology and engineering applications

Because of a lack of suitable instruments, rotational ground motions have not been observed until the last decade. Rotational measurements in the near field of earthquakes in Japan ( Takeo, 1998 ) indicate that rotational ground motions are many times larger than expected from the classical elasticity theory. After failing to obtain useful rotational ground motions (using similar rotational sensors as Takeo did), we deployed a far more sensitive rotational velocity sensor (R-1) at the HGSD station in eastern Taiwan. From 7 December 2004 to 12 November 2006, several hundreds of earthquakes were recorded during our Phase 1 operation. This was mostly a learning exercise to solve field operation problems; Phase 1 operations ended when our two R-1 sensors ceased to operate. A K2+R1 instrument was deployed in the spring of 2007 to start our Phase 2 operation. From 8 May 2007 to 17 February 2008, we observed 52 local earthquakes with good rotational velocity signals (with signal-to-noise ratio >∼5), together with excellent translational acceleration signals (with signal-to-noise ratio >∼10). Unfortunately, field operation was interrupted due to flooding of the HGSD station site in mid-February 2008; we just resumed normal operation in June 2008.

Bulletin of the Seismological Society of America

Earthquake recordings from the 2002 Seattle Seismic Hazard Investigation of Puget Sound (SHIPS), Washington state

This report describes seismic data obtained during the fourth Seismic Hazard Investigation of Puget Sound (SHIPS) experiment, termed Seattle SHIPS . The experiment was designed to study the influence of the Seattle sedimentary basin on ground shaking during earthquakes. To accomplish this, we deployed seismometers over the basin to record local earthquakes, quarry blasts, and teleseisms during the period of January 26 to May 27, 2002. We plan to analyze the recordings to compute spectral amplitudes at each site, to determine the variability of ground motions over the basin. During the Seattle SHIPS experiment, seismometers were deployed at 87 sites in a 110-km-long east-west line, three north-south lines, and a grid throughout the Seattle urban area (Figure 1). At each of these sites, an L-22, 2-Hz velocity transducer was installed and connected to a REF TEK Digital Acquisition System (DAS), both provided by the Program for Array Seismic Studies of the Continental Lithosphere (PASSCAL) of the Incorporated Research Institutes for Seismology (IRIS). The instruments were installed on January 26 and 27, and were retrieved gradually between April 18 and May 27. All instruments continuously sampled all three components of motion (velocity) at a sample rate of 50 samples/sec. To ensure accurate computations of amplitude, we calibrated the geophones in situ to obtain the instrument responses. In this report, we discuss the acquisition of these data, we describe the processing and merging of these data into 1-hour long traces and into windowed events, we discuss the geophone calibration process and its results, and we display some of the earthquake recordings.

Washington

Debris avalanches in the northern California Coast Range triggered by plate boundary earthquakes

Determining the timing and cause for ancient hillslope failures proves difficult in the western United States, yet critical as it ties directly into groundmotion estimates for hazardous events. This knowledge gap is important to confront as hillslope failures are candidates to be triggered by earthquakes along active plate boundaries. We identify two prehistoric, i.e., preinstrumental history, debris avalanches (3 10^6 and 6 10^6 m^3) in the Coast Ranges of northern California. These debris avalanches are well suited for studying the timing (to the exact year) and trigger as the densely forested landscape enables effective use of dendrochronology and high-resolution radiocarbon to pinpoint time of failure, and the legacy of landslide studies in the region provides context for evaluating a seismic trigger as the most likely triggering mechanism. Neither debris avalanche was triggered by the largest precipitation event of the instrumental record in 1964 CE, and the two debris avalanche sites are physiographically suited to accommodate topographic amplification of seismic shaking. Through a suite of geologic and dendrochronologic evidence, we establish the likely time of failure of the two debris avalanches as 1906 CE and 933 CE. The failure of the younger landslide coincides with the 1906 CE San Francisco earthquake on the Pacific/North American plate boundary, while the older landslide dates to 933 CE falling within the broad age range (850-966 CE) of a likely megathrust earthquake on the Cascadia subduction zone as recorded in coastal marshes in Humboldt Bay, California. The precise age on the 933 CE debris avalanche could trim the uncertainty on a Cascadia subduction zone event from 116 years to one year. Utilizing debris avalanche records from sites suitable to record seismic shaking improves understanding of plate-boundary earthquake timing and extent of shaking.

Bulletin of the Seismological Society of America

Are the horizontal-to-vertical spectral ratios of earthquakes and microtremors the same?

We consider the similarities and differences between earthquake and microtremor horizontal‐to‐vertical spectral ratios (eHVSR and mHVSR, respectively) using a dataset of 161 sites in southern California. Quantitative comparisons are made in terms of the eHVSR and mHVSR lognormal median curves, as well as the frequencies and amplitudes associated with the fundamental‐ and higher‐mode resonances where present. The results show only 58% of the eHVSR–mHVSR pairs agree in terms of their median curve and only 25% of the eHVSR–mHVSR pairs agree in terms of shared resonances, which increases to 68% if flat HVSRs are considered equivalent. Furthermore, while the shared resonances match very well in terms of frequency (root mean square error, RMSE, <0.11 Hz), the amplitudes of those resonances do not agree (RMSE >1.6). These findings demonstrate that while eHVSR and mHVSR agree at some sites, they are not equivalent at all sites. To investigate if the agreement between eHVSR and mHVSR could be related to features of the microtremor data, earthquake recordings, and/or the site conditions, three machine learning (ML) models at varying levels of interpretability are presented. The ML models—which include multivariate logistic regression, gradient‐boosted trees, and support vector machines—show only partial success at using site‐specific data to predict whether eHVSR and mHVSR will likely agree in terms of their median curve (accuracy of 78%) and number of resonances (accuracy of 84%). Therefore, we conclude that while eHVSR and mHVSR can be quite similar in terms of resonant frequencies at some sites, they are not identical at all sites. Furthermore, preliminary evidence shows that the agreement of eHVSR and mHVSR can be predicted a priori given features of the microtremor measurements, earthquake recordings, and site conditions, although a larger dataset will be necessary for developing a robust predictive model.

California

The sedimentary record of the 2018 Anchorage Earthquake in Eklutna Lake, Alaska: Calibrating the lacustrine seismograph

The 30 November 2018 M w "> M M w "> w 7.1 Anchorage earthquake caused modified Mercalli intensities of V¼ to V½ at Eklutna Lake (south central Alaska). A few hours after the earthquake, a “dirt streak” was observed on the lake surface, followed by a peak in sediment turbidity values ( ⁠ &#x223C; 80 "> ∼ 80 times normal) at a drinking water facility, which receives water from the lake through a pipe. These observations hint toward turbidity currents triggered by the earthquake in Eklutna Lake. Here, we study 32 short sediment cores retrieved from across Eklutna Lake and observe a millimeter‐to‐centimeter scale turbidite that can be confidently attributed to the 2018 earthquake in all coring locations. X‐ray computed tomography, grain‐size, and color‐spectral analyses of the turbidite show that it shares physical characteristics with the turbidite generated by the 1964 M w "> M w 9.2 Great Alaska earthquake, while it is considerably different from turbidites caused by historical floods. The 2018 turbidite reaches its largest thickness in the inflow‐proximal basin, but when compared to the 1964 turbidite and thereby canceling out local site effects, it is relatively thick in the inflow‐distal sub‐basin. The latter was exposed to stronger shaking during the 2018 earthquake, and this relative thickness trend may therefore be attributed to shaking intensity and gives an indication of the location of the earthquake epicenter relative to the basin axis. Furthermore, in contrast to the 1964 turbidite, which was sourced from both deltas and hemipelagic slopes, the 2018 turbidite was sourced from deltas only, as evidenced by its distribution. These results confirm that while it is generally accepted that shaking intensities of &#x2265; VI "> ≥ VI are needed to trigger turbidity currents from hemipelagic slopes, intensities as low as V¼ can be sufficient to trigger turbidity currents from deltaic slopes. Our results show that proglacial lakes can sensitively record differences in shaking intensity and that investigating deposits from recent earthquakes is crucial to calibrate the lacustrine seismograph.

Alaska

Seismological mapping of fine structure near the base of the Earth's mantle

The Earth's core–mantle boundary (CMB) juxtaposes liquid iron and crystalline silicates, and is a region of large vertical thermal gradients. The D" region, which extends up to 200–300 km above the CMB, often has elevated shear-wave velocity and suggestions of lateral variations in structure 1 . Recent improvements in our ability to assemble and analyse records from regional seismic networks have allowed us to examine long profiles of travel times, amplitudes and waveforms from more than a thousand short-period seismometers 2 . We observe, across Canada and the United States, P waves that have grazed the CMB from the powerful nuclear test in Lop Nor, China, on 21 May 1992. First-arrival travel times and large secondary arrivals are consistent with a 1.5% compressional velocity increase with depth ∼ 130 km above the CMB—about half the thickness of D" in this locality 3 . Our observations, together with evidence for the absence of such a thin, fast layer in neighbouring regions, suggest the presence of lateral heterogeneity in composition or phase at the base of the mantle.

Nature

Proceedings of the Guidelines for Seismometer Testing Workshop, Albuquerque, New Mexico, 9-10 May 2005 ("GST2")

Testing and specification of seismic and earthquake-engineering sensors and recorders has been marked by significant variations in procedures and selected parameters. These variations cause difficulty in comparing such specifications and test results. In July 1989, and again in May 2005, the U.S. Geological Survey hosted international public/private workshops with the goal of defining widely accepted guidelines for the testing of seismological inertial sensors, seismometers, and accelerometers. This document reports the Proceedings of the 2005 workshop and includes as Appendix 6 the report of the 1989 workshop. In a future document, we will attempt to collate and rationalize a single set of formal guidelines for testing and specifying seismic sensors, supplementing Advanced National Seismic System (ANSS) guidelines on instrumentation likely used by ANSS as its standard for verification, acceptance, and intermittent testing, as well as for responses to ANSS instrument requisitions.

Open-File Report

What do data used to develop ground-motion prediction equations tell us about motions near faults?

A large database of ground motions from shallow earthquakes occurring in active tectonic regions around the world, recently developed in the Pacific Earthquake Engineering Center&rsquo;s NGA-West2 project, has been used to investigate what such a database can say about the properties and processes of crustal fault zones. There are a relatively small number of near-rupture records, implying that few recordings in the database are within crustal fault zones, but the records that do exist emphasize the complexity of ground-motion amplitudes and polarization close to individual faults. On average over the whole data set, however, the scaling of ground motions with magnitude at a fixed distance, and the distance dependence of the ground motions, seem to be largely consistent with simple seismological models of source scaling, path propagation effects, and local site amplification. The data show that ground motions close to large faults, as measured by elastic response spectra, tend to saturate and become essentially constant for short periods. This saturation seems to be primarily a geometrical effect, due to the increasing size of the rupture surface with magnitude, and not due to a breakdown in self similarity.

Pure and Applied Geophysics

Kinematic rupture and 3D wave propagation simulations of the 2019 Mw 7.1 Ridgecrest, California, earthquake

We model the kinematic rupture process of the 2019 M w "> M w Mw 7.1 Ridgecrest, California, earthquake using numerical simulations to reproduce the elastodynamic wave field observed by inertial seismometers, high‐rate Global Navigation Satellite System stations, and borehole strainmeters. This was the largest earthquake in Southern California in 20 yr and was widely felt throughout the region. The M w "> M w Mw 7.1 mainshock was part of a large sequence of &#x223C; 30 , 000 "> ∼ 30 , 000 ∼30,000 aftershocks and was notably preceded by an M w "> M w Mw 6.4 foreshock by 34 hr on fault structures that were once poorly understood. A large number of seismic and geodetic instruments measured the rupture process for both events, with many stations located in the near field. Hence, this is a rare opportunity to better understand complex earthquake processes that arise in an immature fault zone using advanced computing. Of the kinematic rupture models that we tested, our preferred is the simplest one that reproduces signals recorded by the three different geophysical datasets; it is composed of four distinct ruptures that progressively migrate to the southeast with delayed initiation times, and typical rupture speeds. This type of model does a better job at matching the recorded ground motions and deformations than does one composed of a continuous rupture with very low‐rupture velocity, as proposed in other studies of this earthquake.

California

Interpreting the tectonic evolution of Pacific Rim margins using plate kinematics and slab window volcanism

The possibility that slab windows might form in the wake of ridge subduction was first conceptualized in the late 1970s as earth scientists explored the implications of plate tectonic theory in three dimensions. Acceptance of slab-window occurrence, however, has been a long time coming (McCrory and Wilson, this issue). With records of both current and past slab-window transients now firmly established through detailed studies of volcanism and tectonism tied to kinematic plate reconstructions, we are able to exploit these spatially and temporally restricted events as paleogeographic markers. Slab-window research encompasses and integrates a broad range of disciplines including earthquake seismology, plate kinematics, volcanic geochemistry, lithospheric dynamics, seismic tomography, and structural geology. In this special issue we use this powerful new tool to constrain and interpret the Cenozoic evolution of Pacific Rim subduction margins from Alaska, to Patagonia, to Antarctica, to Tonga.

Tectonophysics

Some observations on colocated and closely spaced strong ground-motion records of the 1999 Chi-Chi, Taiwan, earthquake

The digital accelerograph network installed in Taiwan produced a rich set of records from the 20 September 1999 Chi-Chi, Taiwan earthquake (Mw 7.6). Teledyne Geotech model A-800 and A-900A* digital accelerographs were colocated at 22 stations that recorded this event. Comparisons of the amplitudes, frequency content, and baseline offsets show that records from several of the A-800 accelerographs are considerably different than those from the colocated A-900A accelerographs. On this basis, and in view of the more thorough predeployment testing of the newer A-900A instruments, we recommend that the records from the A-800 instruments be used with caution in analyses of the mainshock and aftershocks. At the Hualien seismic station two A-900A and one A-800 instruments were colocated, along with a Global Positioning System instrument. Although the records from the two A-900A instruments are much more similar than those from a colocated A-800 instrument, both three-component records contain unpredictable baseline offsets, which produced completely unrealistic ground displacements derived from the accelerations by double integration, as do many of the strong-motion data from this event; the details of the baseline offsets differ considerably on the two three-component records. There are probably numerous sources of the baseline offsets, including sources external to the instruments, such as tilting or rotation of the ground, and sources internal to the instruments, such as electrical or mechanical hysteresis in the sensors. For the two colocated A-900A records at the Hualien seismic station, however, the differences in the baseline offsets suggest that the principal source is some transient disturbance within the instrument. The baseline offsets generally manifest themselves in the acceleration time series as pulses or steps, either singly or in combination. We find a 0.015-Hz low-cut filter can almost completely eliminate the effects of the baseline offsets, but then information regarding the permanent displacements is lost. The causative mechanisms of the baseline offsets are unknown presently. Hence, it is very difficult to recover the permanent displacements from the modern digital records, although for records close to large earthquakes, the signal-to-noise ratio should theoretically be adequate to obtain ground motions with periods of hundreds of seconds. This study reinforces our conclusion from previous studies that the sources of baseline offsets occurring in digital strong-motion records are very complex and often unpredictable, and that, therefore, it is difficult to remove the baseline effects to maximize the information content of the record. The baseline offsets only affect very long period motions (e.g., >20 sec), however, and therefore are of little or no engineering concern.

Bulletin of the Seismological Society of America

Performance of the radial semblance method for the location of very long period volcanic signals

We investigate the performance of a source location method that combines multichannel semblance and particle motions and is being increasingly used to obtain estimates of the source locations of very long period (VLP) seismic signals recorded on volcanoes. The method makes use of the radial particle motions and large wavelengths that characterize the VLP events. To assess the capabilities of this radial semblance method , and to better understand its limitations, we quantify the effects of window length, noise contents of the signal , inaccurate velocity models, receiver coverage, and orientation errors in the horizontal components of the receivers. Our results show that the semblance method performs best when (1) the noise level is low enough to allow a good characterization of the waveforms, (2) the sources are located at distances between one half of the average receiver spacing and about two times the network aperture, and (3) the orientations of the horizontal components of the seismometers are known with relative accuracy. When these requirements are met, the radial semblance method constitutes an adequate tool to obtain preliminary locations of VLP volcanic signals recorded by broadband networks. Moreover, we provide a formula to determine the radial semblance level that should be used to define error regions associated to the estimated source locations.

Bulletin of the Seismological Society of America

The LArge-n Seismic Survey in Oklahoma (LASSO) experiment

In 2016, the U.S. Geological Survey deployed more than 1,800 vertical-component nodal seismometers in Grant County, Oklahoma to study induced seismic activity associated with production of the Mississippi Limestone Play. The LArge-n Seismic Survey in Oklahoma (LASSO) array operated for approximately one month, covering a 25-km-by-32-km region with a nominal station spacing of ~400 m. Primary goals of the deployment were to detect microseismic events not captured by the sparser regional network stations and to provide nearly unaliased records of the seismic wavefield. A more complete record of earthquakes allows us to map the spatiotemporal evolution of induced event sequences and illuminates the structures on which the events occur. Dense records of the seismic wavefield also provide improved measurements of the earthquake source, including focal mechanisms and stress drops. Taken together, we can use these findings to glean insights into the processes that induce earthquakes. Here, we describe the array layout, features of the nodal sensors, data recording configurations, and the field deployment. We also provide examples of earthquake waveforms recorded by the array to illustrate data quality and initial observations. LASSO array data provide a significant resource for understanding the occurrence of earthquakes induced by wastewater disposal.

Kansas, Oklahoma