SEARCH · USGS Science
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.
The project for the International Deployment of Accelerometers (IDA) is intended to provide very low frequency seismic data for the study of Earth structure and source mechanisms. This article outlines preliminary results on attenuation and splitting derived from IDA recordings of a recent large earthquake. ?? 1979 Nature Publishing Group.
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.
The primary purpose of this short note is to formally record the liquefaction prediction (Holzer et al., 1986) made in connection with this predicted earthquake. In addition, this note serves to alert the seismic engineering community to special instrumentation being installed at the prediction site. The instrumentation will consist of 4 downhole accelerometers at depths ranging from 3-30 m, a surface accelerometer, 7 dynamic piezometers distributed in the sand strata between depths of 5 and 15 m, and a network of bench marks for measuring permanent ground deformation.
A reliable method for collection, display, and analysis of low-frequency geophysical data from isolated sites, which can be throughout North and South America and the Pacific Rim, has been developed for use with the Geostationary Operational Environmental Satellite (GOES) system. Geophysical data primarily intended for earthquake hazard and crustal deformation monitoring are digitized with either 12-bit or 16-bit resolution and transmitted every 10 min through a satellite link to a bank of UNIX-based computers in Menlo Park, California. There the data are available for analysis and display within a few seconds of their transmit time. This system provides real-time monitoring of crustal deformation parameters such as tilt, strain, fault displacement, local magnetic field, crustal geochemistry, and water levels, as well as meteorological and other parameters, along faults in California and Alaska, and in volcanic regions in the western United States, Rabaul, and other locations in the New Britain region of the South Pacific. Various mathematical, statistical, and graphical algorithms process the incoming data to detect changes in crustal deformation and fault slip that may indicate the first stages of catastrophic fault failure. Alert trigger levels based on physical models, signal resolution, and previous history have been defined for particular instrument types. Computer-driven remote paging and mail systems are used to notify appropriate personnel when alarm status is reached. The system supports continuous historical records of low-frequency geophysical data, software for extensive analysis of these data, and programs for modeling fault rupture with and without seismic radiation, as well as providing an environment for real-time attempts at earthquake prediction.
In mountain headwater streams, the quality and resilience of summer cold-water habitat is generally regulated by stream discharge, longitudinal stream channel connectivity and groundwater exchange. These critical hydrologic processes are thought to be influenced by the stream corridor bedrock contact depth (sediment thickness), a parameter often inferred from sparse hillslope borehole information, piezometer refusal and remotely sensed data. To investigate how local bedrock depth might control summer stream temperature and channel disconnection (dewatering) patterns, we measured stream corridor bedrock depth by collecting and interpreting 191 passive seismic datasets along eight headwater streams in Shenandoah National Park (Virginia, USA). In addition, we used multi-year stream temperature and streamflow records to calculate several baseflow-related metrics along and among the study streams. Finally, comprehensive visual surveys of stream channel dewatering were conducted in 2016, 2019 and 2021 during summer low flow conditions (124 total km of stream length). We found that measured bedrock depths along the study streams were not well-characterized by soils maps or an existing global-scale geologic dataset where the latter overpredicted measured depths by 12.2 m (mean) or approximately four times the average bedrock depth of 2.9 m. Half of the eight study stream corridors had an average bedrock depth of less than 2 m. Of the eight study streams, Staunton River had the deepest average bedrock depth (3.4 m), the coldest summer temperature profiles and substantially higher summer baseflow indices compared to the other study steams. Staunton River also exhibited paired air and water annual temperature signals suggesting deeper groundwater influence, and the stream channel did not dewater in lower sections during any baseflow survey. In contrast, Paine Run and Piney River did show pronounced, patchy channel dewatering, with Paine Run having dozens of discrete dry channel sections ranging from 1 to greater than 300 m in length. Stream dewatering patterns were apparently influenced by a combination of discrete deep bedrock (20+ m) features and more subtle sediment thickness variation (1–4 m) depending on local stream valley hydrogeology. In combination, these unique datasets show the first large-scale empirical support for existing conceptual models of headwater stream disconnection based on spatially variable underflow capacity and shallow groundwater supply.
On December 11, 1989, geologists working in the crater at Mount St. Helens discovered two thin layers of ash separated by fresh snow-clear evidence that at least two small explosions had occurred recently. The explosions were neither seen nor heard, but on December 7 scientists suspected that a small ash-producing explosion had occurred when seismometers near the volcano recorded a long explosion-like signal, and titlt and displacement meters showed minor deformation of the dome. There were no other large seismic signals to account for the second ash layer, which was most likely associated with one of several smaller signals in early December. The December ash-producing explosions were the first eruptive activity at Mount St. Helens since October 1986. There have been at least five more ash-producing explosions since December 1989, all without recognized seismic or other geophysical precursors. The ash from these explosions appears to be pulverized pieces of dacite dome. The absence of glass shards in the ash suggests that no new magmatic material was ejected. Several of the explosions were accompanied by snow and rock avalanches, pyroclastic flows, ballistic showers, and debris flows. These ash-producing explosions are part of a series of at least 28 explosion-like seismic events that began on August 24, 989. Seismic signals from these events resemble those associated with confirmed ash-producing explosions in April-May 1986. Yet not all of the 1989-1991 events produced ash plumes. Excellent visual observations during four of the events indicated that neither a steam nor ash plume was generated. There is little information about the other events because they occurred when the mountain was not visible, nor was there physical evidence of ashfall or surface changes when scientists visited the crater days to weeks alter. Considerable deformation of the north side of the dome occurred during the series of explosion-like seismic events. Sections of the dome slumped northward and two new vents were formed. However, monitoring the changes associated with individual events was often impossible because several key electronic-distance-meter (EDM) targets and tiltmeters were destroyed by the series of events.
A relational database was developed over a five-year period to support ground motion model (GMM) development for the Next Generation Attenuation-Subduction (NGA-Sub) project. The relational database has components that interact according to a database schema, including a source and path component used to describe attributes of seismic sources in global subduction regions and to compute source-to-site distances, a site component that describes attributes of sites where recordings have been made, and a ground motion component. The source component of the database has information for 1880 earthquakes, mainly from the following regions: the Pacific Northwest region of North America, Alaska and the Aleutian Islands, Japan, Taiwan, New Zealand, South America, Central America, and Mexico. Of the 1880 earthquakes, 88 have finite fault models (FFMs) from the literature that were systematically reviewed, distilled to one more rectangular shapes, and trimmed according to procedures based on percentage of total slip. For earthquakes without FFMs, a simulation routine is used to represent finite fault effects required for distance calculations. This simulation routine was adjusted and made more uniform in its application than in prior NGA projects. All earthquakes are classified as interface, intraslab, shallow crustal, or outer rise, using uniform protocols developed for this project. All earthquakes are also assigned class designations adapted from a prior NGA project for active regions, that allows foreshock, mainshock, and aftershock events to be distinguished. The site component of the database is described in a companion paper (Ahdi et al. 2020 [1]). The ground motion component of the database consists of median – and maximum – horizontal component peak parameters (peak ground acceleration, PGA and peak ground velocity, PGV) and pseudo-spectral accelerations (PSa) at 111 oscillator periods and 11 damping ratios. Response spectra were also computed for the vertical component. Fourier amplitude spectra (FAS) and duration metrics were also computed. The ground motion recordings were obtained from collaborating organizations world-wide as uncorrected (Vol 1) digital recordings, that were corrected (componentspecific low – and high – pass filters and baseline correction, as needed) following Pacific Earthquake Engineering Research Center (PEER)/NGA protocols. The relational database operates on each of these (and other) database components to dynamically draw relevant parameters into a single file, known as a flatfile, that is used by researchers engaged in GMM development. The flatfiles used in model development are being published with the NGA-Sub GMMs as products of the NGA-Sub project.
Shallow shear‐wave velocities ( V s ) sometimes are estimated from joint inversions of horizontal‐to‐vertical (H/V) spectral ratios and surface‐wave dispersion curves derived from ambient noise or small active sources. Here, we evaluate carrying out these inversions using Rayleigh‐wave dispersion curves computed from crustal‐scale P ‐wave seismic refraction data. We use data from the 2014–2015 Eastern North American Margin (ENAM) experiment in Virginia and North Carolina, but similar seismic refraction data sets have been acquired over sedimentary basins of interest for seismic hazard studies, including in major urban areas. The ENAM project deployed a pair of ∼215 km long, northwest–southeast linear arrays with ∼300 m receiver spacing to record 11 dynamite shots, and 80 continuously recording seismometers with 5–6 km spacing along the same arrays to record offshore airguns. The arrays crossed the onland portion of the Atlantic Coastal Plain sediments, which are a seaward‐thickening wedge of Cretaceous and younger sediments deposited mostly on crystalline bedrock. We compute Rayleigh‐wave dispersion curves from 3 to 9 km long portions of the receiver arrays on each side of the dynamite shots, and we compute ambient‐noise H/V ratios from the continuously recording seismometers. We use a genetic inversion algorithm in which forward velocity models in each “generation” are evaluated for misfits compared to the observed data, with subsequent generations constructed from the models with the smallest misfits. Velocities to depths of 500 m are defined well, as shown by a narrow range of velocities in the best‐fit models, by the consistency between multiple inversion runs at a site, and by forward modeling of site responses. The resulting velocity cross‐section of the Coastal Plain strata has seaward‐dipping contours in the thinner portions of the Coastal Plain but smaller dips in the deeper portions. We interpret these results as showing that velocity contours in the ACP strata are influenced by a combination of lithology and overburden pressure. Results demonstrate that existing seismic refraction data have the potential for determining detailed shallow shear‐wave velocity profiles.
Micro-tsunami waves with a maximum amplitude of 4-6 mm were detected with the ocean-bottom pressure gauges on a cabled deep seafloor observatory south of Hokkaido, Japan, following the January 28, 2000 earthquake (M w 6.8) in the southern Kuril subduction zone. We model the observed micro-tsunami and estimate the focal depth and other source parameters such as fault length and amount of slip using grid searching with the least-squares method. The source depth and stress drop for the January 2000 earthquake are estimated to be 50 km and 7 MPa, respectively, with possible ranges of 45-55 km and 4-13 MPa. Focal depth of typical inter-plate earthquakes in this region ranges from 10 to 20 km and stress drop of inter-plate earthquakes generally is around 3 MPa. The source depth and stress drop estimates suggest that the earthquake was an intra-slab event in the subducting Pacific plate, rather than an inter-plate event. In addition, for a prescribed fault width of 30 km, the fault length is estimated to be 15 km, with possible ranges of 10-20 km, which is the same as the previously determined aftershock distribution. The corresponding estimate for seismic moment is 2.7x10 19 Nm with possible ranges of 2.3x10 19 -3.2x10 19 Nm. Standard tide gauges along the nearby coast did not record any tsunami signal. High-precision ocean-bottom pressure measurements offshore thus make it possible to determine fault parameters of moderate-sized earthquakes in subduction zones using open-ocean tsunami waveforms. Published by Elsevier Science B. V.
Although seismicity near Koyna Reservoir (India) has persisted for ~50 years and includes the largest induced earthquake ( M 6.3) reported worldwide, the seismotectonic framework of the area is not well understood. We recorded ~1800 earthquakes from 6 January 2010 to 28 May 2010 and located a subset of 343 of the highest-quality earthquakes using the tomoDD code of Zhang and Thurber (2003) to better understand the framework. We also inverted first arrivals for 3-D Vp , Vs , and Vp / Vs and Poisson's ratio tomography models of the upper 12 km of the crust. Epicenters for the recorded earthquakes are located south of the Koyna River, including a high-density cluster that coincides with a shallow depth (<1.5 km) zone of relatively high Vp and low Vs (also high Vp / Vs and Poisson's ratios) near Warna Reservoir. This anomalous zone, which extends near vertically to at least 8 km depth and laterally northward at least 15 km, is likely a water-saturated zone of faults under high pore pressures. Because many of the earthquakes occur on the periphery of the fault zone, rather than near its center, the observed seismicity-velocity correlations are consistent with the concept that many of the earthquakes nucleate in fractures adjacent to the main fault zone due to high pore pressure. We interpret our velocity images as showing a series of northwest trending faults locally near the central part of Warna Reservoir and a major northward trending fault zone north of Warna Reservoir.
An increase in seismic stations also having microbarographs has led to increased interest in the field of seismoacoustics. A review of the recent advances in this field can be found in Dannemann Dugick et al. (2023) . The goal of this note is to draw the attention of the readers of Dannemann Dugick et al. (2023) to several additional interactions between the solid Earth and atmosphere that have not been classically considered in the field of seismoacoustics. The 15 January 2022 Hunga Tonga–Hunga Ha‘api eruption produced acoustic gravity waves that were recorded globally. For example, the Lamb wave from this eruption produced early‐arriving and long‐lasting tsunami waves. This eruption also provided globally recorded coupling of atmospheric modes with solid Earth modes, providing another example of the complex interactions that can occur at the boundary between the atmosphere and the solid Earth. Even in the absence of large atmospheric signals, collocated pressure sensors at seismic stations can be a useful tool for estimating the local substructure, such at V S 30 , the average shear velocity of the upper 30 m. Finally, at low frequencies, it is possible to use pressure records to correct out atmospheric disturbances recorded on seismometers. We briefly review the aforementioned, nontraditional seismoacoustic topics that we feel are important to consider as part of the full suite of interactions occurring between the solid Earth and atmosphere.
In the seismic design of structures, estimates of design forces are usually provided to the engineer in the form of elastic response spectra. Predictive equations for elastic response spectra are derived from empirical recordings of ground motion. The geometric mean of the two orthogonal horizontal components of motion is often used as the response value in these predictive equations, although it is not necessarily the most relevant estimate of forces within the structure. For some applications it is desirable to estimate the response value on a randomly chosen single component of ground motion, and in other applications the maximum response in a single direction is required. We give adjustment factors that allow converting the predictions of geometric-mean ground-motion predictions into either of these other two measures of seismic ground-motion intensity. In addition, we investigate the relation of the strike-normal component of ground motion to the maximum response values. We show that the strike-normal component of ground motion seldom corresponds to the maximum horizontal-component response value (in particular, at distances greater than about 3 km from faults), and that focusing on this case in exclusion of others can result in the underestimation of the maximum component. This research provides estimates of the maximum response value of a single component for all cases, not just near-fault strike-normal components. We provide modification factors that can be used to convert predictions of ground motions in terms of the geometric mean to the maximum spectral acceleration (SaMaxRot) and the random component of spectral acceleration (SaArb). Included are modification factors for both the mean and the aleatory standard deviation of the logarithm of the motions.
Seismic soundings of Lake Vostok have been performed by the Polar Marine Geological Research Expedition in collaboration with the Russian Antarctic Expedition since the early 1990s. The seismograms recorded show at least two relatively closely spaced reflections associated with the lake bottom. These were initially interpreted as boundaries of a layer of unconsolidated sediments at the bottom of the lake. A more recent interpretation suggests that the observed reflections are side echoes from the rough lake bottom, and that there are no unconsolidated sediments at the bottom of the lake. The major goal of this paper is to reveal the nature of those reflections by testing three hypotheses of their origin. The results show that some of the reflections, but not all of them, are consistent with the hypothesis of a non-flat lake bottom along the source-receiver line (2D case). The reflections were also evaluated as side echoes from an adjacent sloping interface, but these tests implied unreasonably steep slopes (at least 8 degrees) at the lake bottom. The hypothesis that is the most compatible with seismic data is the presence of a widespread layer of unconsolidated sediments at the bottom of Lake Vostok. The modeling suggests the presence of a two hundred meter thick sedimentary layer with a seismic velocity of 1700 -1900 m/sec in the southern and middle parts of the lake. The sedimentary layer thickens in the northern basin to ~350 m
The Hawaiian Volcano Observatory (HVO) summary presents seismic data gathered during the year and a chronological narrative describing the volcanic events. The seismic summary is offered without interpretation as a source of preliminary data. It is complete in the sense that all data for events of M≥1.5 routinely gathered by the Observatory are included. The emphasis in collection of tilt and deformation data has shifted from quarterly measurements at a few water-tube tilt stations (“wet” tilt) to a larger number of continuously recording borehole tiltmeters, repeated measurements at numerous spirit-level tilt stations (“dry” tilt), and surveying of level and trilateration networks. Because of the large quantity of deformation data now gathered and differing schedules of data reduction, the seismic and deformation summaries are published separately. The HVO summaries have been published in various forms since 1956. Summaries prior to 1974 were issued quarterly, but cost, convenience of preparation and distribution, and the large quantities of data dictated an annual publication beginning with Summary 74 for the year 1974. Summary 86 (the introduction of CUSP at HVO) includes a description of the seismic instrumentation, calibration, and processing used in recent years. The present summary includes enough background information on the seismic network and processing to allow use of the data and to provide an understanding of how they were gathered.
The Hawaiian Volcano Observatory (HVO) summary presents seismic data gathered during the year and a chronological narrative describing the volcanic events. The seismic summary is offered without interpretation as a source of preliminary data. It is complete in the sense that all data for events of M>1.5 routinely gathered by the Observatory are included. The emphasis in collection of tilt and deformation data has shifted from quarterly measurements at a few water-tube tilt stations ("wet" tilt) to a larger number of continuously recording borehole tiltmeters, repeated measurements at numerous spirit-level tilt stations ("dry" tilt), and surveying of level and trilateration networks. Because of the large quantity of deformation data now gathered and differing schedules of data reduction, the seismic and deformation summaries are published separately. The HVO summaries have been published in various forms since 1956. Summaries prior to 1974 were issued quarterly, but cost, convenience of preparation and distribution, and the large quantities of data dictated an annual publication beginning with Summary 74 for the year 1974. Summary 86 (the introduction of CUSP at HVO) includes a description of the seismic instrumentation, calibration, and processing used in recent years. The present summary includes enough background information on the seismic network and processing to allow use of the data and to provide an understanding of how they were gathered.
The Hawaiian Volcano Observatory (HVO) summary presents seismic data gathered during the year and a chronological narrative describing the volcanic events. The seismic summary is offered without interpretation as a source of preliminary data. It is complete in the sense that most data for events of M= 1.5 routinely gathered by the Observatory are included. The emphasis in collection of tilt and deformation data has shifted from quarterly measurements at a few water-tube tilt stations ('wet' tilt) to a larger number of continuously recording borehole tiltmeters, repeated measurements at numerous spirit-level tilt stations ('dry' tilt), and surveying of level and trilateration networks. Because of the large quantity of deformation data now gathered and differing schedules of data reduction, the seismic and deformation summaries are published separately. The HVO summaries have been published in various forms since 1956. Summaries prior to 1974 were issued quarterly, but cost, convenience of preparation and distribution, and the large quantities of data dictated an annual publication beginning with Summary 74 for the year 1974. Summary 86 (the introduction of CUSP at HVO) includes a description of the seismic instrumentation, calibration, and processing used in recent years. The present summary includes background information on the seismic network and processing to allow use of the data and to provide an understanding of how they were gathered.
The Integrated Ocean Drilling Program Expedition 318 to the Wilkes Land margin of Antarctica recovered a sedimentary succession ranging in age from lower Eocene to the Holocene. Excellent stratigraphic control is key to understanding the timing of paleoceanographic events through critical climate intervals. Drill sites recovered the lower and middle Eocene, nearly the entire Oligocene, the Miocene from about 17 Ma, the entire Pliocene and much of the Pleistocene. The paleomagnetic properties are generally suitable for magnetostratigraphic interpretation, with well-behaved demagnetization diagrams, uniform distribution of declinations, and a clear separation into two inclination modes. Although the sequences were discontinuously recovered with many gaps due to coring, and there are hiatuses from sedimentary and tectonic processes, the magnetostratigraphic patterns are in general readily interpretable. Our interpretations are integrated with the diatom, radiolarian, calcareous nannofossils and dinoflagellate cyst (dinocyst) biostratigraphy. The magnetostratigraphy significantly improves the resolution of the chronostratigraphy, particularly in intervals with poor biostratigraphic control. However, Southern Ocean records with reliable magnetostratigraphies are notably scarce, and the data reported here provide an opportunity for improved calibration of the biostratigraphic records. In particular, we provide a rare magnetostratigraphic calibration for dinocyst biostratigraphy in the Paleogene and a substantially improved diatom calibration for the Pliocene. This paper presents the stratigraphic framework for future paleoceanographic proxy records which are being developed for the Wilkes Land margin cores. It further provides tight constraints on the duration of regional hiatuses inferred from seismic surveys of the region.
The giant tsunami that swept the Pacific from Alaska to Antarctica in 1946 was generated along one of three Alaska Trench instrumentally recorded aftershock areas following great and giant earthquakes. Aftershock areas were investigated during the past decade with multibeam bathymetry, ocean bottom seismograph wide‐angle seismic, reprocessed legacy, and new seismic reflection images. Summarized and updated here are previous papers and additional data. Tectonic structures collocated with aftershock area boundaries indicate possible lengths of rupture in future great earthquakes. NE aftershock area boundaries relate to subducted lower plate structures whereas the SW zone upper plate retains Beringian structural relicts. The lower to middle slope transition separating a stronger continental framework rock from a weaker accreted prism occurs along splay fault zones previously interpreted as backstops in seismic images. Damage zones along splay faults are generally 1‐km‐wide dipping typically 21°. Splays form slip paths from the plate interface to the seafloor much shorter than the 3–4° dipping plate interface beneath the frontal prism. Associated seafloor vent structures indicate overpressured fluids at depth. Splay fault dip and its rigid hanging wall impart greater seafloor uplift than the accreted prism per unit of slip making them effective tsunami generators. Backstop splay fault zones (BSFZs) run along the entire Alaska Trench. Beneath the frontal prism, active bend faults add rugosity to the plate interface and km high relief is commonly imaged in reprocessed legacy and new seismic data. The 1946 Unimak great (M8.6) earthquake epicenter is located near the BSFZ.