USGS Science⌕ Search

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 145 records · Page 8Linked to original sources
USGS2001

Relocation of Wyoming mine production blasts using calibration explosions

An important requirement for a comprehensive seismic monitoring system is the capability to accurately locate small seismic events worldwide. Accurate event location can improve the probability of determining whether or not a small event, recorded predominantly by local and regional stations, is a nuclear explosion. For those portions of the earth where crustal velocities are not well established, reference event calibration techniques offer a method of increased locational accuracy and reduced locational bias. In this study, data from a set of mining events with good ground-truth data in the Powder River Basin region of eastern Wyoming are used to investigate the potential of event calibration techniques in the area. Results of this study are compared with locations published in the prototype International Data Center’s Reviewed Event Bulletin (REB). A Joint Hypocenter Determination (JHD) method was applied to a set of 23 events. Four of those events with superior ground-truth control (mining company report or Global Positioning System data) were used as JHD reference events, Nineteen (83%) of the solutions converged and the resulting set of station-phase travel-time corrections from the JHD results was then tested. When those travel-time corrections were applied individually to the four events with good ground-truth control, the average locational error reduced the original REB location error from 16.1 km to 5.7 km (65% improvement). The JHD locations indicated reduced locational bias and all of the individual error ellipses enclosed the actual known event locations. Given a set of well-recorded calibration events, it appears that the JHD methodology is a viable technique for improving locational accuracy of future small events where the location depends on arrival times from predominantly local and/or regional stations. In this specific case, the International Association of Seismology and the Physics of the Earth’s Interior (IASPEI) travel-time tables, coupled with JHDderived travel-time corrections, may obviate the need for an accurately known regional velocity structure in the Powder River Basin region.

Montana, Wyoming↗
USGS2014

Tracking changes in volcanic systems with seismic Interferometry

The detection and evaluation of time-dependent changes at volcanoes form the foundation upon which successful volcano monitoring is built. Temporal changes at volcanoes occur over all time scales and may be obvious (e.g., earthquake swarms) or subtle (e.g., a slow, steady increase in the level of tremor). Some of the most challenging types of time-dependent change to detect are subtle variations in material properties beneath active volcanoes. Although difficult to measure, such changes carry important information about stresses and fluids present within hydrothermal and magmatic systems. These changes are imprinted on seismic waves that propagate through volcanoes. In recent years, there has been a quantum leap in the ability to detect subtle structural changes systematically at volcanoes with seismic waves. The new methodology is based on the idea that useful seismic signals can be generated “at will” from seismic noise. This means signals can be measured any time, in contrast to the often irregular and unpredictable times of earthquakes. With seismic noise in the frequency band 0.1–1 Hz arising from the interaction of the ocean with the solid Earth known as microseisms, researchers have demonstrated that cross-correlations of passive seismic recordings between pairs of seismometers yield coherent signals (Campillo and Paul 2003; Shapiro and Campillo 2004). Based on this principle, coherent signals have been reconstructed from noise recordings in such diverse fields as helioseismology (Rickett and Claerbout 2000), ultrasound (Weaver and Lobkis 2001), ocean acoustic waves (Roux and Kuperman 2004), regional (Shapiro et al. 2005; Sabra et al. 2005; Bensen et al. 2007) and exploration (Draganov et al. 2007) seismology, atmospheric infrasound (Haney 2009), and studies of the cryosphere (Marsan et al. 2012). Initial applications of ambient seismic noise were to regional surface wave tomography (Shapiro et al. 2005). Brenguier et al. (2007) were the first to use ambient noise tomography (ANT) to map the 3D structure of a volcanic interior (at Piton de la Fournaise). Subsequent studies have imaged volcanoes with ANT at Okmok (Masterlark et al. 2010), Toba (Stankiewicz et al. 2010), Katmai (Thurber et al. 2012), Asama (Nagaoka et al. 2012), Uturuncu (Jay et al. 2012), and Kilauea (Ballmer et al. 2013b). In addition, Ma et al. (2013) have imaged a scatterer in the volcanic region of southern Peru by applying array techniques to ambient noise correlations. Prior to and in tandem with the development of ANT, researchers discovered that repeating earthquakes, which often occur at volcanoes, could be used to monitor subtle time-dependent changes with a technique known as the doublet method or coda wave interferometry (CWI) (Poupinet et al. 1984; Roberts et al. 1992; Ratdomopurbo and Poupinet 1995; Snieder et al. 2002; Pandolfi et al. 2006; Wegler et al. 2006; Martini et al. 2009; Haney et al. 2009; De Angelis 2009; Nagaoka et al. 2010; Battaglia et al. 2012; Erdem and Waite 2005; Hotovec-Ellis et al. 2014). Chaput et al. (2012) have also used scattered waves from Strombolian eruption coda at Erebus volcano to image the reflectivity of the volcanic interior with body wave interferometry. However, CWI in its original form was limited in that repeating earthquakes, or doublets, were not always guaranteed to occur. With the widespread use of noise correlations in seismology following the groundbreaking work by Campillo and Paul (2003) and Shapiro et al. (2005), it became evident that the nature of the ambient seismic field, due to its oceanic origin, enabled the continuous monitoring of subtle, time-dependent changes at both fault zones (Wegler and Sens-Schönfelder 2007; Brenguier et al. 2008b; Wegler et al. 2009; Sawazaki et al. 2009; Tatagi et al. 2012) and volcanoes (Sens-Schönfelder and Wegler 2006; Brenguier et al. 2008a) without the need for repeating earthquakes. Seismic precursors to eruptions based on ambient noise we

Book chapter↗
USGS1996

Landslides triggered by the 1994 Northridge, California, earthquake

The 17 January 1994 Northridge, California, earthquake ( M w = 6.7) triggered more than 11,000 landslides over an area of about 10,000 km 2 . Most of the landslides were concentrated in a 1000-km 2 area that included the Santa Susana Mountains and the mountains north of the Santa Clara River valley. We mapped landslides triggered by the earthquake in the field and from 1:60,000-nominal-scale aerial photography provided by the U.S. Air Force and taken the morning of the earthquake; these mapped landslides were subsequently digitized and plotted in a GIS-based format. Most of the triggered landslides were shallow (1- to 5-m thick), highly disrupted falls and slides within weakly cemented Tertiary to Pleistocene clastic sediment. Average volumes of these types of landslides were less than 1000 m 3 , but many had volumes exceeding 100,000 m 3 . The larger disrupted slides commonly had runout paths of more than 50 m, and a few traveled as far as 200 m from the bases of steep parent slopes. Deeper (>5-m thick) rotational slumps and block slides numbered in the tens to perhaps hundreds, a few of which exceeded 100,000 m 3 in volume. Most of these were reactivations of previously existing landslides. The largest single landslide triggered by the earthquake was a rotational slump/block slide having a volume of 8 × 10 6 m 3 . Analysis of the mapped landslide distribution with respect to variations in (1) landslide susceptibility and (2) strong shaking recorded by hundreds of instruments will form the basis of a seismic landslide hazard analysis of the Los Angeles area.

California↗
USGS2016

Noise reduction in long‐period seismograms by way of array summing

Long‐period (>100 s period) seismic data can often be dominated by instrumental noise as well as local site noise. When multiple collocated sensors are installed at a single site, it is possible to improve the overall station noise levels by applying stacking methods to their traces. We look at the noise reduction in long‐period seismic data by applying the time–frequency phase‐weighted stacking method of Schimmel and Gallart (2007) as well as the phase‐weighted stacking (PWS) method of Schimmel and Paulssen (1997) to four collocated broadband sensors installed in the quiet Albuquerque Seismological Laboratory underground vault. We show that such stacking methods can improve vertical noise levels by as much as 10 dB over the mean background noise levels at 400 s period, suggesting that greater improvements could be achieved with an array involving multiple sensors. We also apply this method to reduce local incoherent noise on horizontal seismic records of the 2 March 2016 M w 7.8 Sumatra earthquake, where the incoherent noise levels at very long periods are similar in amplitude to the earthquake signal. To maximize the coherency, we apply the PWS method to horizontal data where relative azimuths between collocated sensors are estimated and compared with a simpler linear stack with no azimuthal rotation. Such methods could help reduce noise levels at various seismic stations where multiple high‐quality sensors have been deployed. Such small arrays may also provide a solution to improving long‐period noise levels at Global Seismographic Network stations.

Bulletin of the Seismological Society of America↗
USGS2000

Data for Quaternary faults, liquefaction features, and possible tectonic features in the Central and Eastern United States, east of the Rocky Mountain Front

The USGS is currently leading an effort to compile published geological information on Quaternary faults, folds, and earthquake-induced liquefaction in order to develop an internally consistent database on the locations, ages, and activity rates of major earthquake-related features throughout the United States. This report is the compilation for such features in the Central and Eastern United States (CEUS), which for the purposes of the compilation, is defined as the region extending from the Rocky Mountain Front eastward to the Atlantic seaboard. A key objective of this national compilation is to provide a comprehensive database of Quaternary features that might generate strong ground motion and therefore, should be considered in assessing the seismic hazard throughout the country. In addition to printed versions of regional and individual state compilations, the database will be available on the World-Wide Web, where it will be readily available to everyone. The primary purpose of these compilations and the derivative database is to provide a comprehensive, uniform source of geological information that can by used to complement the other types of data that are used in seismic-hazard assessments. Within our CEUS study area, which encompasses more than 60 percent of the continuous U.S., we summarize the geological information on 69 features that are categorized into four classes (Class A, B, C, and D) based on what is known about the feature's Quaternary activity. The CEUS contains only 13 features of tectonic origin for which there is convincing evidence of Quaternary activity (Class A features). Of the remaining 56 features, 11 require further study in order to confidently define their potential as possible sources of earthquake-induced ground motion (Class B), whereas the remaining features either lack convincing geologic evidence of Quaternary tectonic faulting or have been studied carefully enough to determine that they do not pose a significant seismic hazard (Classes C and D). The correlation between historical seismicity and Quaternary faults and liquefaction features in the CEUS is generally poor, which probably reflects the long return times between successive movements on individual structures. Some Quaternary faults and liquefaction features are located in aseismic areas or where historical seismicity is sparse. These relations indicate that the record of historical seismicity does not identify all potential seismic sources in the CEUS. Furthermore, geological studies of some currently aseismic faults have shown that the faults have generated strong earthquakes in the geologically recent past. Thus, the combination of geological information and seismological data can provide better insight into potential earthquake sources and thereby, contribute to better, more comprehensive seismic-hazard assessments.

Open-File Report↗
USGS2021

Why do my squiggles look funny? A gallery of compromised seismic signals

Seismic instruments are highly sensitive and capable of recording a large range of different Earth signals. The high sensitivity of these instruments also makes them prone to various failures. Although many failures are very obvious, such as a dead channel, there are other more subtle failures that easily go unnoticed by both network operators and data users. This work documents several different types of failure modes in which the instrument is no longer faithfully recording ground‐motion data. Although some of these failure modes make the data completely unusable, there are also a number of failures in which the data can still be used for certain applications. Of course, the ideal situation is to identify as soon as possible when data become compromised and to have the network operator fix the station. However, knowing how the data became compromised can also help data users to identify if the data can still be used for their particular application. This work in no way attempts to exhaustively document recording failures but rather to communicate examples and equip the reader with ways of identifying failure modes.

Seismological Research Letters↗
USGS2025

Nodal seismic deployment on Mauna Loa volcano, Hawaii: Dataset and preliminary insights

Mauna Loa is the largest active volcano on Earth, comprising ∼51% of the Island of Hawai‘i’s landmass and posing significant risks to the island’s communities, infrastructure, and natural environment. Historical eruptions have produced lava flows that have reached the ocean in as little as 3 hr. The timing and location of such lava flows in the past 200 yr underscore how critical determination of the location and geometry of magma storage and structure is for volcanic hazard assessment and eruption forecasting. Now, after nearly 38 yr of volcanic quiescence, Mauna Loa has erupted again. On 27 November 2022, fissures initiated within the summit caldera and then migrated to the northeast rift zone, where they generated a large lava flow that threatened a major highway. To improve our understanding of the geometry of this magma system, we deployed a temporary nodal array on Mauna Loa in the summer of 2024. This increased our seismic coverage sufficiently to image this magmatic system. This nodal array consists of 33 seismometers distributed on and around the volcano and was deployed for over three months to record seismic signals. The primary objective of this project is to resolve the high‐resolution seismic velocity structure and characterize seismic features associated with magma storage and ascent pathways. In this article, we present an overview of the deployment, evaluate the quality of the data, and show example recordings to evaluate the suitability of the data set for future seismic investigations, including earthquake relocation, seismic tomography, and receiver function analysis. Comparisons with nearby permanent broadband and short‐period seismic stations demonstrate that the nodal array recorded high‐quality waveforms, making it a valuable resource for constraining the magmatic system beneath Mauna Loa at multiple scales.

Hawaii↗
USGS1989

Three-dimensional records of surface displacement on the Superstition Hills fault zone associated with the earthquakes of 24 November 1987

Seven quadrilaterals, constructed at broadly distributed points on surface breaks within the Superstition Hills fault zone, were repeatedly remeasured after the pair of 24 November 1987 earthquakes to monitor the growing surface displacement. Changes in the dimensions of the quadrilaterals are recalculated to right-lateral and extensional components at millimeter resolution, and vertical components of change are resolved at 0.2 mm precision. The displacement component data for four of the seven quadrilaterals record the complete fault movement with respect to an October 1986 base. These data fit with remarkable agreement the power law

California↗
USGS1978

Mineral resources of the Cabinet Mountains Wilderness, Lincoln and Sanders Counties, Montana

This report describes the differential array, of seismometers recently installed at the Hollister, California, Municipal Airport. Such an array of relatively closely spaced seismometers has already been installed in El Centro and provided useful information for both engineering and seismological applications from the 1979 Imperial Valley earthquake. Differential ground motions, principally due to horizontally propagating surface waves, are important in determining the stresses in such extended structures as large mat foundations for nuclear power stations, dams, bridges and pipelines. Further, analyses of the records of the 1979 Imperial Valley earthquake from the differential array have demonstrated the utility of short-baseline array data in tracking the progress of the rupture wave front of an earthquake.

Open-File Report↗
USGS2023

Remote characterization of the 12 January 2020 eruption of Taal Volcano, Philippines, using seismo-acoustic, volcanic lightning, and satellite observations

On 12 January 2020, an eruption began on the shores of the Main Crater Lake (MCL) of Taal Volcano—a caldera system on the southern end of Luzon Island in the Philippines. Taal, one of the most active volcanoes in the Philippines, is located 30 km south of Manila—a major metropolitan area with a population of 13.5 million people. Eruptive activity intensified throughout the day on 12 January, producing prolific volcanic lightning, ashfall, and a sustained plume that reached 16–17 km altitude. The chronology of events was well documented by the Philippine Institute of Volcanology and Seismology and the Tokyo Volcanic Ash Advisory Center. The wealth of data collected during the eruption provides a unique opportunity to investigate how the combination of different remote sensing methods may complement local observations and monitoring. Remote systems tend to provide lower resolution data but are also less likely to be compromised by the eruptive activity, thus providing continuous records of eruptive processes. Here, we present a postevent analysis of the 12 January activity, including data from long‐range lightning, infrasound, and seismic arrays located at distances up to several thousands of kilometers from the volcano. By combining these datasets, we distinguish five phases of activity and infer a major shift in eruption behavior around 12:00 on 12 January (UTC). The remote observations suggest that the most of the water within the MCL ( ⁠∼42 million m 3⁠ ) was vaporized and incorporated into the volcanic plume within the first 12 hr of the eruption.

Taal Volcano↗
USGS2005

Development of ground-motion prediction equations relevant to shallow-mining-induced seismicity in the Trial Mountain area, Emery County, Utah

To provide a basis for assessing the seismic hazard to the Joes Valley Dam due to future coal mining in the nearby Cottonwood Tract, central Utah, we developed ground-motion prediction relations using data recorded by a seismic network, established and operated by the University of Utah Seismograph Stations. The network was centered on the Trail Mountain coal mine, located adjacent to the Cottonwood Tract. From late 2000 until early 2001, this network recorded numerous mining-induced events with magnitudes as large as 2.17. The ground motion from these events, recorded at hypocentral distances ranging from about 500 m to approximately 10 km, were well suited to developing new ground-motion prediction relations, especially when augmented by data from a M 4.2 earthquake in the Willow Creek mine, about 50 km north of Trail Mountain. Using a two-stage regression analysis, we determined prediction relations for peak acceleration, peak velocity, and pseudovelocity response spectra, at 5% damping, for periods of 0.1, 0.2, 0.5, 1.0, and 2.0 s. To illustrate the potential seismic hazard at the Joes Valley dam, we used these ground-motion relations to predict a peak velocity of 6.8 cm/s due to an earthquake with the probable maximum magnitude of 3.9, at a hypocentral distance of 1 km, recorded at a rock site typical for this region. This result does not take into account the site response at the dam.

Bulletin of the Seismological Society of America↗
USGS2017

Recurrent Holocene movement on the Susitna Glacier Thrust Fault: The structure that initiated the Mw 7.9 Denali Fault earthquake, central Alaska

We conducted a trench investigation and analyzed pre‐ and postearthquake topography to determine the timing and size of prehistoric surface ruptures on the Susitna Glacier fault (SGF), the thrust fault that initiated the 2002 M w 7.9 Denali fault earthquake sequence in central Alaska. In two of our three hand‐excavated trenches, we found clear evidence for a single pre‐2002 earthquake (penultimate earthquake [PE]) and determined an age of 2210±420 cal. B.P. (2 σ ) for this event. We used structure‐from‐motion software to create a pre‐2002‐earthquake digital surface model (DSM) from 1:62,800‐scale aerial photography taken in 1980 and compared this DSM with postearthquake 5‐m/pixel Interferometric Synthetic Aperature Radar topography taken in 2010. Topographic profiles measured from the pre‐earthquake DSM show features that we interpret as fault and fold scarps. These landforms were about the same size as those formed in 2002, so we infer that the PE was similar in size to the initial ( M w 7.2) subevent of the 2002 sequence. A recurrence interval of 2270 yrs and dip slip of ∼4.8 m yield a single‐interval slip rate of ∼1.8 mm/yr. The lack of evidence for pre‐PE deformation indicates probable episodic (clustering) behavior on the SGF that may be related to strain migration among other similarly oriented thrust faults that together accommodate shortening south of the Denali fault. We suspect that slip‐partitioned thrust‐triggered earthquakes may be a common occurrence on the Denali fault system, but documenting the frequency of such events will be very difficult, given the lack of long‐term paleoseismic records, the number of potential thrust‐earthquake sources, and the pervasive glacial erosion in the region.

Alaska↗
USGS2024

Limited evidence of late Quaternary tectonic surface deformation in the eastern Tennessee seismic zone, USA

The ~300-km-long eastern Tennessee seismic zone (ETSZ), USA, is the second-most seismically active region east of the Rocky Mountains. Seismicity generally occurs below the Paleozoic fold-and-thrust belt within the Mesoproterozoic basement, at depths of 5–26 km, and earthquake magnitudes during the instrumental record have been moment magnitude ( M w )≤4.8. Evidence of surface deformation may not exist or be difficult to detect because of the vegetated and soil-mantled landscape, landslides, locally steep topography, anthropogenic landscape modification, or long, irregular recurrence intervals between surface-rupturing earthquakes. Despite the deep seismicity, analog models indicate that accumulation of strike-slip or oblique-slip displacement at depth could be expected to propagate upward through the Paleozoic section, producing a detectable surficial signal of distributed faulting. To identify potential surface deformation, we interrogated the landscape at different spatial scales. We evaluated morphotectonic and channel metrics, such as channel sinuosity and catchment-scale hypsometry. Additionally, we mapped possible fault-related topographic features on 1-m lidar. Finally, we integrated our observations with available bedrock and Quaternary surficial mapping and subsurface geophysical data. At a regional scale, most morphotectonic and channel metrics have a strong lithologic control. Within smaller regions of similar lithology, we observe changes in landscape metrics like channel sinuosity and catchment-scale hypsometry that spatially correlate with new lineaments identified in this study and previously mapped east–west Cenozoic faults. These faults have apparent left-lateral offsets, are optimally oriented to slip in the current stress field, and match kinematics from recent focal mechanisms, but do not clearly preserve evidence of late Pleistocene or Holocene tectonic surface deformation. Most newly mapped lineaments might be explained by either tectonic or non-tectonic origins, such as fluvial or karst processes. We also re-evaluated a previously described paleoseismic site and interpret that the exposure does not record evidence of late Pleistocene faulting but instead is explained by fluvial stratigraphy.

Tennessee↗
USGS1968

A method for estimating the uncertainty of seismic velocities measured by refraction techniques

Time residuals from 7 5-km segments of 18 crustal seismic-refraction profiles in the Basin and Range province are used to investigate the validity of the linear-regression model and to make large sample estimates of the variance in the travel time distributions. A formula for unbiased estimates of velocity uncertainty is derived, assuming a linear trend with distance for the variances of the travel-time distributions. If the recording units are symmetric about the center of the recording interval, this formula is equivalent to the one derived assuming the variances are equal. At the 95-per cent confidence level the chi-squared test implied 84 per cent of the time-residual samples were inconsistent with the hypothesis that their parent populations had Gaussian distributions. If the number of recording locations expceeds 8, confidence limits computed without the Gaussian assumption suggest the departures from normality are not significant for velocity uncertainty estimates. The large sample estimates of the time-residual populations may be applicable to other areas. This evidence motivated the development of a method, requiring very little numerical calculation, for estimating uncertainties in velocities. The method requires, in addition to the large sample estimates of the travel time variances, information on the quality of the data, the location of the recording interval, and the number of recording units. The method is useful for the design of new experiments and independent estimates of uncertainty reported in the literature.

Arizona, California, Idaho, Nevada, Oregon, Utah↗
USGS2011

PQLX: A seismic data quality control system description, applications, and users manual

We present a detailed description and users manual for a new tool to evaluate seismic station performance and characteristics by providing quick and easy transitions between visualizations of the frequency and time domains. The software is based on the probability density functions (PDF) of power spectral densities (PSD) (McNamara and Buland, 2004) and builds on the original development of the PDF stand-alone software system (McNamara and Boaz, 2005) and the seismological data viewer application PQL (IRIS-PASSCAL Quick Look) and PQLII (available through the IRIS PASSCAL program: http://www.passcal.nmt.edu/content/pql-ii-program-viewing-data). With PQLX (PQL eXtended), computed PSDs are stored in a MySQL database, allowing a user to access specific time periods of PSDs (PDF subsets) and time series segments through a GUI-driven interface. The power of the method and software lies in the fact that there is no need to screen the data for system transients, earthquakes, or general data artifacts, because they map into a background probability level. In fact, examination of artifacts related to station operation and episodic cultural noise allow us to estimate both the overall station quality and a baseline level of Earth noise at each site. The output of this analysis tool is useful for both operational and scientific applications. Operationally, it is useful for characterizing the current and past performance of existing broadband stations, for conducting tests on potential new seismic station locations, for evaluating station baseline noise levels (McNamara and others, 2009), for detecting problems with the recording system or sensors, and for evaluating the overall quality of data and metadata. Scientifically, the tool allows for mining of PSDs for investigations on the evolution of seismic noise (for example, Aster and others, 2008; and Aster and others, 2010) and other phenomena. Currently, PQLX is operational at several organizations including the USGS National Earthquake Information Center (NEIC), the USGS Albuquerque Seismological Laboratory (ASL), and the Incorporated Research Institutions in Seismology (IRIS) Data Management Center (DMC) for station monitoring and instrument response quality control. The PQLX system is available to the community at large through the U.S. Geological Survey (USGS) (http://ehpm-earthquake.wr.usgs.gov/research/software/pqlx.php) and IRIS (http://www.iris.edu/software/pqlx). Also provided is a fully searchable website for bug reporting and enhancement requests (http://wush.net/bugzilla/PQLX). The first part of this document aims to describe and illustrate some of the features and capabilities of the software. The second part of this document is a detailed users manual that covers installation procedures, system requirements, operations, bug reporting, and software components (Appendix).

Open-File Report↗
USGS1976

Physical constraints on sounds generated by very small earthquakes

Vertical vibrations of the ground surface due to elastic waves from an earthquake will generate acoustic pressure waves in the overlying atmosphere. Sufficiently intense vibrations of the ground at frequencies greater than about 20 HZ can produce pressure waves that are within the audible range of many animals, including humans. People in the epicentral region of moderate to large earthquakes, for instance, commonly relate hearing sounds accompanying individual earthquakes (here we exclude sounds associated with mechanical resonances in structures such as houses). Descriptions of these sounds range from the report of distant guns to the rumble of thunder or the rushing of a wind (Davison, 1938; Lanchow Seismological Brigade, Appendix II). In many instances, the sounds are described as starting several seconds before the earthquake is felt. Experimental and theoretical verification that such earthquake sounds are generated by high-frequency ground motion from local earthquakes is based on simultaneous acoustic and seismic recordings recently obtained during a swarm of earthquakes in the Imperial Valley, California (Hill and others, 1976). These results show that for local earthquakes in the magnitude range 2 to-3, audible sounds are generated by the first arriving compressional (P) wave while perceptible shaking begins with the larger but slower shear (S) wave, explaining the reports that "earthquake sounds" are frequently heard several seconds before the earthquake is felt. A commonly proposed explanation for alleged anomalous animal behavior hours to days prior to large earthquakes suggests that animals with acute, high-frequency hearing respond to weak '''ultrasonic" sounds generated by small earthquakes occurring in the epicentral region of the impending event. Such earthquakes are presumably too small to be either felt or heard by humans in the region and too small to be routinely recognized as earthquakes on seismograms recorded on a local seismograph network. Evidence that such "popping and cracking" may occur in the region of an impending earthquake comes from laboratory experiments in which the rate of microfracturing in a rock sample subjected to large stress differences increases dramatically before catastrophic failure of the sample (Scholz, 1968; Stesky, 1975). In this paper we investigate the following question: What are the constraints on earthquake size (fault dimension and displacement) and hypocentral distance (distance from the earthquake focus to a point on the Earth's surface) such that animals noted for their acute hearing might react to the sound generated by the event while the same event would go undetected by humans (either by feeling, hearing, or routine identification on a local seismograph network)? Factors to be considered include: 1) the nature of the displacement spectral amplitudes of elastic waves radiated from an earthquake, 2) the effects of attenuation on the spectral components of elastic (seismic) waves over propagation paths from the focus to the Earth's surface, 3) the appropriate scaling of spectral amplitudes from moderate to very small earthquakes, 4) the coupling of elastic waves in the Earth to acoustic waves in the atmosphere, and 5) the threshold-of-hearing characteristics for various animals, including man. In considering these factors, we will find that, except for extremely small, shallow "earthquakes" (fracture dimension on the order of 10 cm and focal depths on the order of 10 m), acoustic emissions from premonitory fracturing ought to be equally audible to both people and animals in the immediate area of the fracturing.

Open-File Report↗
USGS2019

Establishing high-frequency noise baselines to 100 Hz based on millions of power spectra from IRIS MUSTANG

Advances in seismic instrumentation have enabled data to be recorded at increasing sample rates. This has in turn created a need to establish higher-frequency baselines for assessing data quality, as the widely-used New High (NHNM) and Low Noise Models (NLNM) of Peterson (1993) do not extend to frequencies above 10 Hz. To provide a baseline for higher frequencies (10-100 Hz), we examine power spectral density probability density functions (PSDPDFs) for high-sample-rate stations available from the Incorporated Research Institutions for Seismology Data Services (IRIS DS) MUSTANG quality control system. We compute high-frequency high and low noise baselines by matching the appropriate composite PSDPDF percentile points to NHNM and NLNM power levels at overlapping frequencies (1-10 Hz) and then extending to higher frequencies (10-100 Hz) with piecewise linear fits to the matching PSDPDF percentile. We find that the Peterson NLNM remains an accurate representation of the lower bound of global ambient Earth noise since it is matched by only 0.1% of Global Seismographic Network (GSN) PSDs. We present high-frequency high and low noise baselines intended primarily for use by temporary networks targeting high-frequency signals (e.g. monitoring of aftershocks or induced seismicity) based on statistics of PSDPDFs from all publicly available high-sample-rate data. Most publicly-available high-sample-rate data is recorded by temporary deployments, and the experiment design and scientific targets of these deployments strongly influence the observed statistical distribution of high-frequency noise. We anticipate that the noise baselines presented here will be useful in automated quality control of high-sample-rate seismic data. However, we note that establishing a low noise model that accurately represents the lowest possible ambient Earth noise at frequencies up to 100 Hz will require additional continuous high-sample-rate data from high-quality permanent stations in low-noise environments.

Bulletin of the Seismological Society of America↗
USGS2006

Site response and attenuation in the Puget Lowland, Washington State

Simple spectral ratio (SSR) and horizontal-to-vertical (HN) site-response estimates at 47 sites in the Puget Lowland of Washington State document significant attenuation of 1.5- to 20-Hz shear waves within sedimentary basins there. Amplitudes of the horizontal components of shear-wave arrivals from three local earthquakes were used to compute SSRs with respect to the average of two bedrock sites and H/V spectral ratios with respect to the vertical component of the shear-wave arrivals at each site. SSR site-response curves at thick basin sites show peak amplifications of 2 to 6 at frequencies of 3 to 6 Hz, and decreasing spectra amplification with increasing frequency above 6 Hz. SSRs at nonbasin sites show a variety of shapes and larger resonance peaks. We attribute the spectral decay at frequencies above the amplification peak at basin sites to attenuation within the basin strata. Computing the frequency-independent, depth-dependent attenuation factor (Qs,int) from the SSR spectral decay between 2 and 20 Hz gives values of 5 to 40 for shallow sedimentary deposits and about 250 for the deepest sedimentary strata (7 km depth). H/V site responses show less spectral decay than the SSR responses but contain many of the same resonance peaks. We hypothesize that the H/V method yields a flatter response across the frequency spectrum than SSRs because the H/V reference signal (vertical component of the shear-wave arrivals) has undergone a degree of attenuation similar to the horizontal component recordings. Correcting the SSR site responses for attenuation within the basins by removing the spectral decay improves agreement between SSR and H/V estimates.

Bulletin of the Seismological Society of America↗