SEARCH · Search USGS

Results for “Bulletin Seismological Society America”

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 667 records · Page 37Linked to original sources

USGS2021

Probabilistic fault displacement hazard assessment (PFDHA) for nuclear installations according to IAEA safety standards

In the last 10 yr, the International Atomic Energy Agency (IAEA) revised its safety standards for site evaluations of nuclear installations in response to emerging fault displacement hazard evaluation practices developed in Member States. New amendments in the revised safety guidance (DS507) explicitly recommend fault displacement hazard assessment, including separate approaches for candidate new sites versus existing sites. If there is insufficient basis to conclusively determine that a fault is not capable of surface displacement at an existing site, then a probabilistic fault displacement hazard analysis (PFDHA) is recommended to better characterize the hazard. This new recommendation has generated the need for the IAEA to provide its Member States with guidance on performing PFDHA, including its formulation and implementation. This article provides an overview of current PFDHA state‐of‐practice for nuclear installations that is consistent with the new IAEA safety standards. We also summarize progress in an ongoing international PFDHA benchmark project that will ultimately provide technical guidance to Member States for conducting site‐specific fault displacement hazard assessments.

Bulletin of the Seismological Society of America↗

USGS2021

Improved scaling relationships for seismic moment and average slip of strike-slip earthquakes incorporating fault slip rate, fault width and stress drop

We develop a self‐consistent scaling model relating magnitude M w to surface rupture length ( ⁠L E ⁠ ), surface displacement D E ⁠ , and rupture width W E ⁠ , for strike‐slip faults. Knowledge of the long‐term fault‐slip rate S F improves magnitude estimates. Data are collected for 55 ground‐rupturing strike‐slip earthquakes that have geological estimates of L E ⁠ , D E ⁠ , and S F⁠ , and geophysical estimates of W E ⁠ . We begin with the model of Anderson et al. (2017) , which uses a closed form equation for the seismic moment of a surface‐rupturing strike‐slip fault of arbitrary aspect ratio and given stress drop, Δτ C ⁠ . Using W E estimates does not improve M w estimates. However, measurements of D E plus the relationship between Δτ C and surface slip provide an alternate approach to study W E ⁠ . A grid of plausible stress drop and width pairs were used to predict displacement and earthquake magnitude. A likelihood function was computed from within the uncertainty ranges of the corresponding observed M w and D E values. After maximizing likelihoods over earthquakes in length bins, we found the most likely values of W E for constant stress drop; these depend on the rupture length. The best‐fitting model has the surprising form W E ∝logL E —a gentle increase in width with rupture length. Residuals from this model are convincingly correlated to the fault‐slip rate and also show a weak correlation with the crustal thickness. The resulting model thus supports a constant stress drop for ruptures of all lengths, consistent with teleseismic observation. The approach can be extended to test other observable factors that might improve the predictability of magnitude from a mapped fault for seismic hazard analyses.

Bulletin of the Seismological Society of America↗

USGS2021

A big problem for small earthquakes: Benchmarking routine magnitudes and conversion relationships with coda-envelope-derived Mw in southern Kansas and northern Oklahoma

Earthquake magnitudes are widely relied upon measures of earthquake size. Although moment magnitude ( ⁠ M w "> M w Mw ⁠ ) has become the established standard for moderate and large earthquakes, difficulty in reliably measuring seismic moments for small (generally M w < 4 "> M w < 4 Mw<4 ⁠ ) earthquakes has meant that magnitudes for these events remain plagued by a patchwork of inconsistent measurement scales. Because of this, magnitudes of small earthquakes and statistics derived from them can be biased. Furthermore, because small earthquakes are much more numerous than large ones, many applications, such as seismic hazard modeling, depend critically on analysis of events characterized by magnitudes other than M w "> M w Mw ⁠ . To assess this problem, we apply coda envelope analysis to reliably determine moment magnitudes for a case study of small earthquakes from northern Oklahoma and southern Kansas. Not surprisingly, we find significant differences among M L "> M L ML ⁠ , m b L g "> m b L g mbLg ⁠ , and M w "> M w Mw for M ∼2–4 earthquakes examined here. More troublingly, we find that relations designed to convert other magnitudes to M w "> M w Mw ⁠ , which are relied upon for important applications such as seismic hazard analysis, often increase rather than decrease this bias for our dataset. In our case study, we find that converted magnitudes can result in a systematic bias sometimes exceeding 0.5 magnitude units, a difference that typically corresponds to a factor of ∼3 in seismicity rate. Moreover, we find a correspondingly large bias in Gutenberg–Richter b ‐values, controlled primarily by inaccurate magnitude scaling in the conversion relationships. Although this study focuses on a relatively small geographic area, we can expect that similar issues exist with varying severity in other regions. Therefore, magnitudes of small earthquakes and their associated statistics, including seismicity rates and b ‐values, should be treated with caution.

Kansas, Oklahoma↗

USGS2021

Refinements to the Graves–Pitarka kinematic rupture generator, including a dynamically consistent slip‐rate function, applied to the 2019 Mw 7.1 Ridgecrest earthquake

The main objective of this study is to develop physics‐based constraints on the spatiotemporal variation of the slip‐rate function using a simplified dynamic rupture model. First, we performed dynamic rupture modeling of the 2019 M w 7.1 Ridgecrest, California, earthquake, to analyze the effects of depth‐dependent stress and material friction on slip rate. Then, we used our modeling results to guide refinements to the slip‐rate function that were implemented in the Graves–Pitarka kinematic rupture generation technique. The dynamic ruptures were computed on a surface‐rupturing, planar strike‐slip fault that includes a weak (negative to low‐stress‐drop) zone in the upper 4 km of the crust. Below the weak zone, we placed high‐stress‐drop patches designed to mirror the large‐slip areas seen in various rupture model inversions of the event. The locations of the high‐stress‐drop patches and the hypocenter were varied in multiple realizations to investigate how changing the dynamic conditions affected the resulting rupture kinematics, in particular, the slip rate. From these simulations, we observed a systematic change in the shape of the slip‐rate function from Kostrov type below the weak zone to a predominantly symmetric shape within the weak zone, along with a depth‐dependent reduction of peak slip rate. We generalized these shallow rupture features into a depth‐dependent parametric variation of the slip‐rate function and implemented it in the Graves–Pitarka kinematic rupture model generator. The performance of the updated kinematic approach was then verified in 0–4 Hz simulations of the M w 7.1 Ridgecrest earthquake, which showed that incorporating the depth‐dependent variation in the shape of the slip‐rate function improves the fit to the observed near‐fault ground motions in the 0.5–3 s period range.

California↗

USGS2021

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

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

Bulletin of the Seismological Society of America↗

USGS2022

Active‐source interferometry in marine and terrestrial environments: Importance of directionality and stationary phase

We utilize active‐source seismic interferometry with dense seismic arrays both offshore and onland to explore the utility of this method to create virtual sources and reveal body‐wave reflections in these two different environments. We first utilize data from an ocean‐bottom cable (OBC) array in the Gulf of Mexico with equal numbers of sources (160 airgun shots) and receivers (160 ocean‐bottom four‐component sensors). We next use data from a geophone array across the Bighorn Mountains of Wyoming with many receivers (1300 vertical‐component geophones) but a small number of sources (14 borehole active‐source shots). We find that the OBC virtual source results, which produce strong reflections from sub‐seafloor structures, are far superior to the onland results which lack usable reflections, and we explore reasons for these differences through a set of selective stacking approaches. We present techniques to account for the direction the seismic waves travel (directionality) and stationary phase and show that improvements can be made when incorporating these corrections. Although interferometric methods are based on assumptions of large numbers of widely distributed actual sources, we find that selective exclusion of potentially problematic source–receiver pairs can yield improved results. These geometric adjustments to active‐source interferometry methods have utility for dense‐nodal‐array surveys that are now common in academic studies, but that often suffer from sparse source geometry.

Bulletin of the Seismological Society of America↗

USGS2021

The impact of 3D finite‐fault information on ground‐motion forecasting for earthquake early warning

We identify aspects of finite‐source parameterization that strongly affect the accuracy of estimated ground motion for earthquake early warning (EEW). EEW systems aim to alert users to impending shaking before it reaches them. The U.S. West Coast EEW system, ShakeAlert, currently uses two algorithms based on seismic data to characterize the earthquake’s location, magnitude, and origin time, treating it as a point or line source. From this information, ShakeAlert calculates shaking intensity and alerts locations where shaking estimates exceed a threshold. Several geodetic EEW algorithms under development would provide 3D finite‐fault information. We investigate conditions under which this information produces sufficiently better intensity estimates to potentially improve alerting. Using scenario crustal and subduction interface sources, we (1) identify the most influential source geometry parameters for an EEW algorithm’s shaking forecast, and (2) assess the intensity alert thresholds and magnitude ranges for which more detailed source characterization affects alert accuracy. We find that alert regions determined using 3D‐source representations of correct magnitude and faulting mechanism are generally more accurate than those obtained using line sources. If a line‐source representation is used and magnitude is calculated from the estimated length, then incorrect length estimates significantly degrade alert region accuracy. In detail, the value of 3D‐source characterization depends on the user’s chosen alert threshold, tectonic regime, and faulting style. For the suite of source models we tested, the error in shaking intensity introduced by incorrect geometry could reach levels comparable to the intrinsic uncertainty in ground‐motion calculations (e.g., 0.5–1.3 modified Mercalli intensity [MMI] units for MMI 4.5) but, especially for crustal sources, was often less. For subduction interface sources, 3D representations substantially improved alert area accuracy compared to line sources, and incorrect geometry parameters were more likely to cause error in calculated shaking intensity that exceeded uncertainties.

Bulletin of the Seismological Society of America↗

USGS2021

Post-seismic relaxation following the 2019 Ridgecrest, California, earthquake sequence

The 2019 Ridgecrest, California, earthquake sequence involved predominantly right‐lateral strike slip on a northwest–southeast‐trending subvertical fault in the 6 July M 7.1 mainshock, preceded by left‐lateral strike slip on a northeast–southwest‐trending subvertical fault in the 4 July M 6.4 foreshock. To characterize the postseismic deformation, we assemble displacements measured by Global Positioning System (GPS) and Interferometric Synthetic Aperture Radar. The geodetic measurements illuminate vigorous postseismic deformation for at least 21 months following the earthquake sequence. The postseismic transient deformation is particularly well constrained from survey‐mode GPS (sGPS) in the epicentral region carried out during the weeks after the mainshock. We interpret these observations with mechanical models including afterslip and viscoelastic relaxation of the lower crust and mantle asthenosphere. During the first 21 months, up to several centimeters of horizontal motions are measured at continuous GPS and sGPS sites, with amplitude that diminishes slowly with distance from the mainshock rupture, suggestive of deeper afterslip or viscoelastic relaxation. We find that although afterslip involving right‐lateral strike slip along the mainshock fault traces and their deeper extensions reach a few decimeters, most postseismic deformation is attributable to viscoelastic relaxation of the lower crust and mantle. Within the Basin and Range crust and mantle, we infer a transient lower crust viscosity several times that of the mantle asthenosphere. The transient mantle asthenosphere viscosity is ∼ 1.3 × 10 17    Pa   s ">

California↗

USGS2022

Empirical map-based nonergodic models of site response in the greater Los Angeles area

We develop empirical estimates of site response at seismic stations in the Los Angeles area using recorded ground motions from 414 M 3–7.3 earthquakes in southern California. The data are from a combination of the Next Generation Attenuation‐West2 project, the 2019 Ridgecrest earthquakes, and about 10,000 newly processed records. We estimate site response using an iterative mixed‐effects residuals partitioning approach, accounting for azimuthal variations in anelastic attenuation and potential bias due to spatial clusters of colocated earthquakes. This process yields site response for peak ground acceleration, peak ground velocity, and pseudospectral acceleration relative to a 760 m/s shear‐wave velocity ( ⁠ V S "> V V S "> s ) reference condition. We employ regression kriging to generate a spatially continuous site response model, using the linear site and basin terms from Boore et al. (2014) as the background model, which depend on V S "> V V S "> s V S "> 30 ⁠ and depth to the 1 km/s V S "> V V S "> s isosurface. This is different from past approaches to nonergodic models, in which spatially varying coefficients are regressed. We validate the model using stations in the Community Seismic Network (CSN) that are in the middle of our model spatial domain but were not considered in model development, finding strong agreement between the interpolated model and CSN data for long periods. Our model could be implemented in regional seismic hazard analyses, which would lead to improvements especially at long return periods. Our site response model also has potential to improve both ground‐motion accuracy and warning times for the U.S. Geological Survey ShakeAlert earthquake early warning (EEW) system. For a point‐source EEW simulation of the 1994 M 6.7 Northridge earthquake, our model produces ground motions more consistent with the ground‐truth ShakeMap and would alert areas with high population density such as downtown Los Angeles at lower estimated magnitudes (i.e., sooner) than an ergodic model for a modified Mercalli intensity 4.5 alerting threshold.

California↗

USGS2022

Seismic background noise levels across the continental United States from USArray Transportable Array: The influence of geology and geography

Since 2004, the most complete estimate of background noise levels across the continental U.S. was attained using 61 broadband seismic stations to calculate power spectral density (PSD) probability density functions. To improve seismic noise estimates across the U.S., we examine vertical component seismic data from the EarthScope USArray Transportable Array seismic network that rolled across the U.S. and southeastern Canada between 2004 and 2015 and form a large (10 TB) PSD database from 1679 stations that contains no smoothing or binning of the spectral estimates. Including station outages, our database has a mean of 98.9% data completeness, and we present maps showing the spatial and temporal variability of seismic noise in six bands of interest between 0.2- and 75-s period. At 0.2 s period, seismic noise across the eastern U.S. is predominantly anthropogenically generated and may be subsequently amplified more than 20 decibels in the sandy and water-saturated sediments of the southeastern U.S. Coastal Plain and Mississippi Embayment. In these sediments, 1 s noise shows similar amplification and is generated through a variety of mechanisms including cultural activity throughout Kentucky and the southeastern Appalachian Mountains, lake waves around the Great Lakes, and ocean waves throughout New England, the Pacific Northwest, and Florida. Both 0.2 and 1 s noise levels are the lowest in the Intermountain West portion of the U.S. We attribute this to a combination of installations on crystalline rocks and reduced population density. Finally, we find that sensors emplaced in sandy, water-saturated sediments observe median, diurnal variations in vertical component power at 18 to 75 s period, which we infer arise through local deformation driven by pressure variations. Ultimately, our results underscore that for shallow (<5 m depth) sensor installation, bedrock provides superior broadband noise performance compared to unconsolidated sediments.

Bulletin of the Seismological Society of America↗

USGS2022

Partitioning ground motion uncertainty when conditioned on station data

Rapid estimation of earthquake ground shaking and proper accounting of associated uncertainties in such estimates when conditioned on strong‐motion station data or macroseismic intensity observations are crucial for downstream applications such as ground failure and loss estimation. The U.S. Geological Survey ShakeMap system is called upon to fulfill this objective in light of increased near‐real‐time access to strong‐motion records from around the world. Although the station data provide a direct constraint on shaking estimates at specific locations, these data also heavily influence the uncertainty quantification at other locations. This investigation demonstrates methods to partition the within‐ (phi) and between‐event (tau) uncertainty estimates under the observational constraints, especially when between‐event uncertainties are heteroscedastic. The procedure allows the end users of ShakeMap to create separate between‐ and within‐event realizations of ground‐motion fields for downstream loss modeling applications in a manner that preserves the structure of the underlying random spatial processes.

Bulletin of the Seismological Society of America↗

USGS2022

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

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

Bulletin of the Seismological Society of America↗

USGS2022

An efficient, analytic solution using order statistics for probabilistic seismic‐hazard assessment without the Poisson assumption

Standard approaches to probabilistic seismic‐hazard assessment (PSHA) assume that earthquakes are random, independent events that follow a Poisson distribution of occurrences in a given time period ( Cornell, 1968 ). To overcome the limitations of the Poisson assumption, such as ignoring earthquake clustering, we introduce an analytic method for PSHA that uses order statistics to allow for arbitrary distributions of earthquake occurrence. Cornell (1968) used the Poisson assumption to achieve a computationally efficient method that enables users to explore the impact of parameters used in earthquake occurrence and ground‐motion models. We apply our order statistics method to the highly clustered seismicity associated with caldera collapses at Kīlauea and explore the general implications of non‐Poisson behavior for PSHA. We find that non‐Poisson behavior has the greatest impact for high probabilities of exceedance, low‐mean rates of occurrence, and multiple exceedances. Those conditions can be important for applications such as operating standards for buildings and infrastructure engineering, standards for temporary structures and during construction, the insurance industry, the design of earthquake early warning, and to assess hazards due to clustered processes such as aftershock sequences and earthquake swarms. The commonly used rate of exceedance hides the difference between the hazard due to a non‐Poisson distribution and a Poisson distribution with the same mean rate of earthquakes. Thus, including non‐Poisson behavior in PSHA means that we must plot and discuss PSHA results as the probability and not the rate of exceedance.

Bulletin of the Seismological Society of America↗

USGS2022

Earthquake early warning for estimating floor shaking levels of tall buildings

This article investigates methods to improve earthquake early warning (EEW) predictions of shaking levels for residents of tall buildings. In the current U.S. Geological Survey ShakeAlert EEW system, regions far from an epicenter will not receive alerts due to low predicted ground‐shaking intensities. However, residents of tall buildings in those areas may still experience significant shaking due to the acceleration amplification caused by tall buildings’ dynamic behavior, as recently experienced by residents of the 52‐story building in downtown Los Angeles (DTLA) during the 2019 M 7.1 Ridgecrest earthquake. Using more than 400 recorded response data acquired from 77 instrumented buildings in California, here we compare the Federal Emergency Management Agency (FEMA) P‐58 and American Society of Civil Engineers (ASCE) 7‐16 simplified equations for peak floor acceleration (PFA), finding that the ASCE estimation is close to the median of data recorded in large and long‐distance events, whereas the current FEMA estimation is not suitable. In the second part of this article, four instrumented tall buildings in DTLA are extensively studied, and the performance of the simplified and response spectrum (RS) methods giving both an estimation of the free‐field horizontal peak ground acceleration (PGA) and pseudospectral acceleration is evaluated. The results show that the RS method is as accurate as the response history analysis as long as the ground‐motion RS is accurate, whereas the ASCE 7‐16 prediction is conservative. However, when ground‐motion RS or PGA is estimated for DTLA using a ground‐motion model (GMM), the performance of the RS method significantly degrades due to underestimation by the GMM at long periods. The results of this study imply that a nonergodic GMM, which may give more accurate prediction in Los Angeles, could improve the results for PFA when the building’s behavior is dominated by a few long‐period fundamental modes, as is the case for the 52‐story building in DTLA.

California↗

USGS2022

Alert optimization of the PLUM earthquake early warning algorithm for the western United States

We determine an optimal alerting configuration for the propagation of local undamped motion (PLUM) earthquake early warning (EEW) algorithm for use by the U.S. ShakeAlert system covering California, Oregon, and Washington. All EEW systems should balance the primary goal of providing timely alerts for impactful or potentially damaging shaking while limiting alerts for shaking that is too low to be of concern (precautionary alerts). The PLUM EEW algorithm forward predicts observed ground motions to nearby sites within a defined radius without accounting for attenuation, avoiding the earthquake source parameter estimation step of most EEW algorithms. PLUM was originally developed in Japan where the alert regions and ground motions for which alerts are issued differ from those implemented by ShakeAlert. We compare predicted ground motions from PLUM to ShakeMap‐reported ground motions for a set of 22 U.S. West Coast earthquakes of magnitude 4.4–7.2 and evaluate available warning times. We examine a range of prediction radii (20–100 km), thresholds used to issue an alert (alert threshold), and levels of impactful or potentially damaging shaking (target threshold). We find optimal performance when the alert threshold is close to the target threshold, although higher target ground motions benefit from somewhat lower alert thresholds to ensure timely alerts. We also find that performance, measured as the cost reduction that a user can achieve, depends on the user’s tolerance for precautionary alerts. Users with a low target threshold and high tolerance for precautionary alerts achieve optimal performance when larger prediction radii (60–100 km) are used. In contrast, users with high target thresholds and low tolerance for precautionary alerts achieve better performance for smaller prediction radii (30–60 km). Therefore, setting the PLUM prediction radius to 60 km balances the needs of many users and provides warning times of up to ∼20 s.

California, Oregon, Washington↗

USGS2022

Waveform signatures of earthquakes located close to the subducted Gorda Plate interface

Complex seismic velocity structure near the earthquake source can affect rupture dynamics and strongly modify the seismic waveforms recorded near the fault. Fault‐zone waves are commonly observed in continental crustal settings but are less clear in subduction zones due to the spatial separation between seismic stations and the plate boundary fault. We observed anomalously long duration S waves from earthquake clusters located near the interface of the subducted Gorda plate north of the Mendocino triple junction. In contrast, earthquakes located just a few kilometers below each cluster show impulsive S waves. A nodal array experiment was conducted around the Northern California Seismic Network station KCT for two months to investigate the origin of the complex S waves. Beamforming analysis shows that the S waves contain three arrivals that have different horizontal slownesses, which we term S 1, S 2, and S coda. Similar analysis on P waves also show two arrivals with different horizontal slownesses, which we term P 1 and P 2. P 1 and S 1 have larger horizontal slowness than P 2 and S 2, respectively, indicating that the phase pairs are body waves with different ray paths. Building upon a seismic refraction profile, we construct 1D velocity models and test different thicknesses and V P / V S ratios for the subducted oceanic crust. The arrival times and relative slownesses of P 1/ P 2 and S 1/ S 2 phases indicate that they are the direct and the Moho reflected phases, respectively. Their properties are consistent with a crustal thickness of ∼6 km and a moderate V P / V S ratio (∼1.8). The S coda is more difficult to characterize but has a clear dominant frequency that likely reflects the near‐source velocity and attenuation structure. Our study indicates that waveforms from earthquakes near the interface of the subducted slab can be used to infer detailed structural information about the plate‐boundary zone at seismogenic depths.

California↗

USGS2022

Topographic response to simulated Mw 6.5-7.0 earthquakes on the Seattle Fault

We explore the response of ground motions to topography during large crustal fault earthquakes by simulating several magnitude 6.5–7.0 rupture scenarios on the Seattle fault, Washington State. Kinematic simulations are run using a 3D spectral element code and a detailed seismic velocity model for the Puget Sound region. This model includes realistic surface topography and a near‐surface low‐velocity layer; a mesh spacing of ∼30 m at the surface allows modeling of ground motions up to 3 Hz. We simulate 20 earthquake scenarios using different slip distributions and hypocenter locations on a planar fault surface. Results indicate that average ground motions in simulations with and without topography are similar. However, shaking amplification is common at topographic highs, and more than a quarter of all sites experience short‐period (≤2 s) ground‐motion amplification greater than 25%–35%, compared with models without topography. Comparisons of peak ground velocity at the top and bottom of topographic features demonstrate that amplification is sensitive to period, with the greatest amplifications typically manifesting near a topographic feature’s estimated resonance frequency and along azimuths perpendicular to its primary axis of elongation. However, interevent variability in topographic response can be significant, particularly at shorter periods (<1 s). We do not observe a clear relationship between source centroid‐to‐site azimuths and the strength of topographic amplification. Overall, our results suggest that although topographic resonance does influence the average ground motions, other processes (e.g., localized focusing and scattering) also play a significant role in determining topographic response. However, the amount of consistent, significant amplification due to topography suggests that topographic effects should likely be considered in some capacity during seismic hazard studies.

Washington↗

USGS2022

Seismometer records of ground tilt induced by debris flows

A change in surface loading causes the Earth’s surface to deform. Mass movements, such as debris flows, can cause a tilt large enough to be recorded by nearby instruments, but the signal is strongly dependent on the mass loading and subsurface parameters. Specifically designed sensors for such measurements (tiltmeters) are cumbersome to install. Alternatively, broadband seismometers record translational motion and also tilt signals, often at periods of tens to hundreds of seconds. Their horizontal components are thereby the most sensitive to tilt. In this study, we show how to obtain tilt caused by the passing by of debris flows from seismic measurements recorded within tens of meters of the flow and investigate the usefulness of this signal for flow characterization. We investigate the problem on three scales (1) large‐scale laboratory experiments at the U.S. Geological Survey debris‐flow flume, where broadband seismometers and tiltmeters were installed for six 8 – 10    m 3 "> 8 – 10 m 3 experiments, (2) the Illgraben torrent in Switzerland, one of the most active mass wasting sites in the European Alps, where a broadband seismometer placed within a few meters of the channel recorded 15 debris‐flow events with volumes up to 10 5    m 3 "> 10 5 m 3 ⁠ , and (3) Volcán de Fuego, Guatemala, where a broadband seismometer recorded two lahars. We investigate how the tilt signals compare to debris‐flow parameters such as mean normal stresses, usually measured by expensive force plates, and debris‐flow height. We model the elastic ground deformation as the response of an elastic half‐space to a moving surface load. In addition, we use the model with some simplifications to determine the maximum debris‐flow heights of Volcán de Fuego events, where no force plate measurements are available. Finally, we address how and under what assumptions the relatively affordable and straightforward tilt measurements may be utilized to infer debris‐flow parameters, as opposed to force plates and other complicated instrument setups.

Oregon↗