SEARCH · Search USGS

Results for “The Seismic Record”

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

SourceResearch connectionFrom yearThrough year

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

At least 1,009 records · Page 56Linked to original sources

USGS2026

Syn-magmatic subsidence during the early stages of continental rifting in the Mesoproterozoic—A reanalysis of legacy data for the Midcontinent Rift, western Lake Superior

The Midcontinent Rift system (ca. 1.1 Ga) is a 2000-km-long series of elongated volcanic and sedimentary troughs and associated intrusive centers exposed chiefly in the Lake Superior region of North America. The rift system represents a long history of intense magmatism and subsequent sedimentation that was arrested by far-field tectonic events before sea-floor spreading was established. The premature cessation preserved a record of processes related to the beginning of continental rifting. The rift system under Lake Superior has been long studied using seismic-reflection data collected as part of the Great Lakes International Multidisciplinary Program on Crustal Evolution (GLIMPCE). We reexamine GLIMPCE Line C by developing a detailed velocity model for time to depth conversion constrained by other legacy data. We corroborate the model and develop a geologic interpretation using gravity and magnetic modeling and ties to geology mapped onshore. We recognize superposed subsiding sedimentary and volcanic basins for the southern half of the Line C depth section. This interpretation differs from previous paradigms that show major crustal faults that bound half-grabens or full grabens. We conclude that high-velocity (6.9 km/s) intrusive zones rather than major crustal faults border the sides of the basins. We speculate that the volcanic basin represents the initiation of seaward dipping reflectors. The syn-magmatic subsidence can be explained by dike injection and volcanic loading. Discrete lava basins throughout the region likely subsided at different times in a disorganized manner along the rift trend, raising questions about the long-term role of lithospheric thinning and melt generation.

Michigan, Minnesota, Ontario, Wisconsin↗

USGS2016

A rare moderate‐sized (Mw 4.9) earthquake in Kansas: Rupture process of the Milan, Kansas, earthquake of 12 November 2014 and its relationship to fluid injection

The largest recorded earthquake in Kansas occurred northeast of Milan on 12 November 2014 ( M w 4.9) in a region previously devoid of significant seismic activity. Applying multistation processing to data from local stations, we are able to detail the rupture process and rupture geometry of the mainshock, identify the causative fault plane, and delineate the expansion and extent of the subsequent seismic activity. The earthquake followed rapid increases of fluid injection by multiple wastewater injection wells in the vicinity of the fault. The source parameters and behavior of the Milan earthquake and foreshock–aftershock sequence are similar to characteristics of other earthquakes induced by wastewater injection into permeable formations overlying crystalline basement. This earthquake also provides an opportunity to test the empirical relation that uses felt area to estimate moment magnitude for historical earthquakes for Kansas.

Seismological Research Letters↗

USGS2014

Three-dimensional ground-motion simulations of earthquakes for the Hanford area, Washington

This report describes the results of ground-motion simulations of earthquakes using three-dimensional (3D) and one-dimensional (1D) crustal models conducted for the probabilistic seismic hazard assessment (PSHA) of the Hanford facility, Washington, under the Senior Seismic Hazard Analysis Committee (SSHAC) guidelines. The first portion of this report demonstrates that the 3D seismic velocity model for the area produces synthetic seismograms with characteristics (spectral response values, duration) that better match those of the observed recordings of local earthquakes, compared to a 1D model with horizontal layers. The second part of the report compares the response spectra of synthetics from 3D and 1D models for moment magnitude (M) 6.6–6.8 earthquakes on three nearby faults and for a dipping plane wave source meant to approximate regional S-waves from a Cascadia great earthquake. The 1D models are specific to each site used for the PSHA. The use of the 3D model produces spectral response accelerations at periods of 0.5–2.0 seconds as much as a factor of 4.5 greater than those from the 1D models for the crustal fault sources. The spectral accelerations of the 3D synthetics for the Cascadia plane-wave source are as much as a factor of 9 greater than those from the 1D models. The differences between the spectral accelerations for the 3D and 1D models are most pronounced for sites with thicker supra-basalt sediments and for stations with earthquakes on the Rattlesnake Hills fault and for the Cascadia plane-wave source.

Washington↗

USGS2016

Comparisons between vs30 and spectral response for 30 sites in Newcastle, Australia from collocated seismic cone penetrometer, active- and passive-source vs data

Although the time‐averaged shear‐wave velocity down to 30 m depth ( V S 30 ) can be a proxy for estimating earthquake ground‐motion amplification, significant controversy exists about its limitations when used as a single parameter for the prediction of amplification. To examine this question in absence of relevant strong‐motion records, we use a range of different methods to measure the shear‐wave velocity profiles and the resulting theoretical site amplification factors (AFs) for 30 sites in the Newcastle area, Australia, in a series of blind comparison studies. The multimethod approach used here combines past seismic cone penetrometer and spectral analysis of surface‐wave data, with newly acquired horizontal‐to‐vertical spectral ratio, passive‐source surface‐wave spatial autocorrelation (SPAC), refraction microtremor (ReMi), and multichannel analysis of surface‐wave data. The various measurement techniques predicted a range of different AFs. The SPAC and ReMi techniques have the smallest overall deviation from the median AF for the majority of sites. We show that V S 30 can be related to spectral response above a period T of 0.5 s but not necessarily with the maximum amplification according to the modeling done based on the measured shear‐wave velocity profiles. Both V S 30 and AF values are influenced by the velocity ratio between bedrock and overlying sediments and the presence of surficial thin low‐velocity layers (<2 m thick and <150 m/s), but the velocity ratio is what mostly affects the AF. At 0.2< T <0.4 s, the AFs are largely controlled by the surficial geology of a particular site. AF maxima are the highest in the hard classes, which is the inverse of the findings used in the Australian Building Code. Only for T >0.5 s do the amplification curves consistently show higher values for soft site classes and lower for hard classes.

Bulletin of the Seismological Society of America↗

USGS2012

Aftershock seismicity of the 2010 Maule Mw=8.8 Chile, earthquake: Correlation between co-seismic slip models and aftershock distribution?

The 27 February 2010 Maule, Chile (Mw=8.8) earthquake is one of the best instrumentally observed subduction zone megathrust events. Here we present locations, magnitudes and cumulative equivalent moment of the first -2 months of aftershocks, recorded on a temporary network deployed within 2 weeks of the occurrence of the mainshock. Using automatically-determined onset times and a back projection approach for event association, we are able to detect over 30,000 events in the time period analyzed. To further increase the location accuracy, we systematically searched for potential S-wave arrivals and events were located in a regional 2D velocity model. Additionally, we calculated regional moment tensors to gain insight into the deformation history of the aftershock sequence. We find that the aftershock seismicity is concentrated between 40 and 140 km distance from the trench over a depth range of 10 to 35 km. Focal mechanisms indicate a predominance of thrust faulting, with occasional normal faulting events. Increased activity is seen in the outer-rise region of the Nazca plate, predominantly in the northern part of the rupture area. Further down-dip, a second band of clustered seismicity, showing mainly thrust motion, is located at depths of 40–45 km. By comparing recent published mainshock source inversions with our aftershock distribution, we discriminate slip models based on the assumption that aftershocks occur in areas of rapid transition between high and low slip, surrounding high-slip regions of the mainshock.

Geophysical Research Letters↗

USGS1989

Seismic methods for determining earthquake source parameters and lithospheric structure

The seismologic methods most commonly used in studies of earthquakes and the structure of the continental lithosphere are reviewed in three main sections: earthquake source parameter determinations, the determination of earth structure using natural sources, and controlled-source seismology. The emphasis in each section is on a description of data, the principles behind the analysis techniques, and the assumptions and uncertainties in interpretation. Rather than focusing on future directions in seismology, the goal here is to summarize past and current practice as a companion to the review papers in this volume. Reliable earthquake hypocenters and focal mechanisms require seismograph locations with a broad distribution in azimuth and distance from the earthquakes; a recording within one focal depth of the epicenter provides excellent hypocentral depth control. For earthquakes of magnitude greater than 4.5, waveform modeling methods may be used to determine source parameters. The seismic moment tensor provides the most complete and accurate measure of earthquake source parameters, and offers a dynamic picture of the faulting process. Methods for determining the Earth’s structure from natural sources exist for local, regional, and teleseismic sources. One-dimensional models of structure are obtained from body and surface waves using both forward and inverse modeling. Forward-modeling methods include consideration of seismic amplitudes and waveforms, but lack the formal resolution estimates obtained with inverse methods. Two- and three-dimensional lithospheric models are derived using various inverse methods, but at present most of these methods consider only traveltimes of body waves. Controlled-source studies of the Earth’s structure are generally divided by method into seismic refraction/wide-angle reflection and seismic reflection studies. Seismic refraction profiles are usually interpreted in terms of two-dimensional structure by forward modeling of traveltimes and amplitudes. The refraction method gives excellent estimates of seismic velocities, but relatively low resolution of structure. Formal resolution estimates are not possible for models derived from forward modeling, but informal estimates can be obtained by perturbing the best-fitting model. Inversion methods for seismic refraction data for one-dimensional models are well established, and two- and three-dimensional methods, including tomography, have recently been developed. Seismic reflection data provide the highest resolution of crustal structure, and have provided many important geological insights in the past decade. The acquisition and processing of these data have been greatly advanced by the hydrocarbon exploration industry. However, reliable crustal velocity control is generally lacking, and the origin of deep crustal reflections remains unclear, resulting in nonunique interpretations. A new form of lithospheric seismology has recently emerged that combines the advantages of seismic refraction and seismic reflection profiles, and the distinction between the two methods is steadily diminishing. Major challenges for future work will be the collection of data that are more densely sampled in space, and the development of interpretation methods that provide quantitative estimates of the uncertainties in the calculated models.

GSA Memoirs↗

USGS2018

Research to improve ShakeAlert earthquake early warning products and their utility

Earthquake early warning (EEW) is the rapid detection of an earthquake and issuance of an alert or notification to people and vulnerable systems likely to experience potentially damaging ground shaking. The level of ground shaking that is considered damaging is defined by the specific application; for example, manufacturing equipment may experience damage at a lower intensity ground shaking than would cause damage to a building. Along the West Coast of the United States, the warning times for ground shaking could range as high as tens of seconds for moderate levels of ground shaking, or potentially longer, if a lower ground-shaking threshold is used to issue alerts. However, it is not always possible to provide advance warning of ground shaking, particularly for locations close to an earthquake that are most likely to experience very strong ground shaking. EEW alerts may be useful to individuals who can use a few seconds to move to a safe zone and to electromechanical systems that can take automatic actions to reduce damage and injuries. An EEW system, ShakeAlert, has been under development in the United States since 2006. Federal and State governments, as well as the private sector, are now investing in the ShakeAlert prototype system that will, when completed, become an operational public system for the West Coast of the United States. While the current prototype is delivering alerts to test users, improvements to the accuracy, timeliness, and utility of the alerts are needed. For this reason, it is essential that the ShakeAlert system be continuously improved through targeted research, involving not only the current ShakeAlert partner organizations, but also the broader scientific, engineering, and emergencyresponse communities. To this end, this report describes the opportunities for improvement that can be addressed through research and development over the next 5 years. Our recommendations are organized into four areas: (1) understand EEW capabilities and user needs, (2) make alerts as fast and accurate as possible, (3) ensure reliability when it counts, and (4) explore the use of new instrumentation. The first challenge is to understand EEW capabilities and user needs. EEW must deliver actionable information to people and to automated systems to mitigate short- and longterm impacts of damaging ground shaking, so development of EEW must be motivated by the needs of users. Within this challenge, we must study the technical capabilities and limitations of EEW in general, and the ShakeAlert system specifically. This includes development of performance metrics that assess the timeliness and accuracy of alerts to understand the value and utility of the ShakeAlert EEW product(s) for various user groups, including different industry sectors, emergency-management agencies, and the public. Research is needed to define the alerting choices that maximize the utility of the system for users and to determine what the available communication pathways are for providing timely alert information. Additionally, we engage users to assess how alerts will be used by different sectors to mitigate losses and to inform EEW product design. Further, social-science research is needed to develop alert messaging, including what relevant prior and follow-up information are required, to ensure effective use of alerts. The second challenge is to make alerts as fast and as accurate as possible. The timeliness and accuracy of an EEW alert is important because it will set in motion a series of actions and downstream products. An EEW alert will trigger notification across emergency-alert systems and across multiple communication channels to populations in impacted regions. The EEW alert region may grow as the earthquake fault-rupture length increases, and the EEW system’s characterization of it, evolves. We must continue research into new or improved seismic and geodetic waveform-processing methods necessary to rapidly characterize the expected ground shaking and associated uncertainties. It is important to thoroughly evaluate whether new methods improve alerts through more accurate ground-motion estimates and (or) reduced latencies (that is, longer warning times). New methods could include tracking the extent of a large rupture in real time (known as finite-fault algorithms) and ground-motionbased EEW algorithms. Additionally, ground motion predictions could be optimized for each earthquake as the earthquake fault rupture progresses by using, for example, event terms to shift ground-motion curves for more (or less) energetic ruptures. The third challenge is to ensure reliability when it counts. This challenge requires us to explore approaches that assess the expected performance of ShakeAlert across the range of earthquake magnitudes, locations, and depths that may occur within the alerting region. Large, damaging earthquakes and their associated aftershock sequences matter most for hazard and for EEW, but these large-earthquake sequences occur infrequently. We expect ShakeAlert to respond robustly to these large-earthquake sequences despite potentially long periods of relative seismic quiescence in the intervening years, and in spite of inevitable communication challenges that arise during and after a large earthquake. We must develop methods to utilize the broadest available datasets to test EEW performance, including ground-motion data recorded in other parts of the world. The observational period for large, damaging earthquakes in any particular region has been short in comparison to estimated large-earthquake recurrence times. Ground-motion records for very large, damaging western United States events and major aftershock sequences do not yet exist, nor do data exist for all potential sources of noise and spurious signals that ShakeAlert must be “tuned” to reject. In addition, robust synthetic data could provide the flexibility to test a wider range of earthquake magnitude, tectonic-setting, and noise scenarios than are covered by existing observational data. Synthetic ground-motion data must be thoroughly vetted against records of smaller magnitude earthquakes to ensure that they accurately capture both the onset and the amplitude of the ground shaking. The final challenge is to explore the use of new instrumentation. The development of EEW around the world to date has focused on the use of high-quality, scientific-grade seismic and geodetic instrumentation. The use of additional types of instrumentation or information may also improve EEW products by filling gaps in sensor coverage in countries that already have dense seismic networks or enable EEW in countries without such networks. We must keep up with these developments and continuously assess their value in supplementing existing EEW systems, such as ShakeAlert, or enabling EEW where such systems do not exist. Such developments include low-cost instrumentation with microelectromechanical system (MEMS) sensors and global positioning system (GPS)/global navigation satellite system (GNSS) antennas embedded in low-cost consumer electronics, sea-floor seismometers, geodetic instrumentation deployed along the Cascadia and Alaska megathrust margins of western North America, and borehole strainmeters that are already deployed across the region.

Open-File Report↗

USGS2017

Rapid Holocene coastal change revealed by high-resolution micropaleontological analysis, Pamlico Sound, North Carolina, USA

Foraminiferal analyses of 404 contiguous samples, supported by diatom, lithologic, geochronologic and seismic data, reveal both rapid and gradual Holocene paleoenvironmental changes in an 8.21-m vibracore taken from southern Pamlico Sound, North Carolina. Data record initial flooding of a latest Pleistocene river drainage and the formation of an estuary 9000. yr ago. Estuarine conditions were punctuated by two intervals of marine influence from approximately 4100 to 3700 and 1150 to 500. cal. yr BP. Foraminiferal assemblages in the muddy sand facies that accumulated during these intervals contain many well-preserved benthic foraminiferal species, which occur today in open marine settings as deep as the mid shelf, and significant numbers of well-preserved planktonic foraminifera, some typical of Gulf Stream waters. We postulate that these marine-influenced units resulted from temporary destruction of the southern Outer Banks barrier islands by hurricanes. The second increase in marine influence is coeval with increased rate of sea-level rise and a peak in Atlantic tropical cyclone activity during the Medieval Climate Anomaly. This high-resolution analysis demonstrates the range of environmental variability and the rapidity of coastal change that can result from the interplay of changing climate, sea level and geomorphology in an estuarine setting. ?? 2011 University of Washington.

North Carolina↗

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↗

USGS1988

Parkfield, California, liquefaction prediction (USA)

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.

Califonia↗

USGS1989

A satellite-based digital data system for low-frequency geophysical data

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.

Bulletin of the Seismological Society of America↗

USGS2022

Bedrock depth influences spatial patterns of summer baseflow, temperature and flow disconnection for mountainous headwater streams

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.

Virginia↗

USGS1992

Small explosions interrupt 3-year quiescence at Mount St. Helens, Washington

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.

Washington↗

USGS2022

Data resources for NGA-subduction project

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.

Extramural-Authored Publication Paper↗

USGS2024

Using crustal-scale refraction data of joint inversions of Rayleigh-wave dispersion curves and H/V spectral ratios for Atlantic Coastal Plain velocity structure, eastern U.S.

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.

eastern United States↗

USGS2003

Source depth dependence of micro-tsunamis recorded with ocean-bottom pressure gauges: The January 28, 2000 Mw 6.8 earthquake off Nemuro Peninsula, Japan

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.

Earth and Planetary Science Letters↗

USGS2014

Seismicity, faulting, and structure of the Koyna-Warna seismic region, Western India from local earthquake tomography and hypocenter locations

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.

Maharashtra↗

USGS2023

Comment on “A new decade in seismoacoustics (2010–2022)” by Fransiska Dannemann Dugick, Clinton Koch, Elizabeth Berg, Stephen Arrowsmith, and Sarah Albert

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.

Bulletin of the Seismological Society of America↗