Geology ReportsSearch

SEARCH · Geology Reports

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.

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

At least 541 records · Page 30Linked to original sources

Magnetic field variations in Alaska: Recording space weather events on seismic stations in Alaska

Seismometers are highly sensitive instruments to not only ground motion but also many other nonseismic noise sources (e.g., temperature, pressure, and magnetic field variations). We show that the Alaska component of the Transportable Array is particularly susceptible to recording magnetic storms and other space weather events because the sensors used in this network are unshielded and magnetic flux variations are stronger at higher latitudes. We also show that vertical‐component seismic records across Alaska are directly recording magnetic field variations between 40 and 800 s period as opposed to actual ground motion during geomagnetic events with sensitivities ranging from 0.004 to 0.48    ( m / s 2 ) / T "> 0.48 ( m / s 2 ) / T 0.48 (m/s2)/T ⁠ . These sensitivities were found on a day where the root mean square variation in the magnetic field was 225 nT. Using a method developed by Forbriger (2007 , his section 3.1), we show that improving vertical seismic resolution of an unshielded sensor by as much as 10 dB in the 100–400 s period band using magnetic data from a collocated three‐component magnetometer is possible. However, due to large spatial variations in Earth’s magnetic field, this methodology becomes increasingly ineffective as the distance between the seismometer and magnetometer increases (no more than 200 km separation). A potential solution to this issue may be to incorporate relatively low‐cost magnetometers as an additional environmental data stream at high‐latitude seismic stations. We demonstrate that the Bartington Mag‐690 sensors currently deployed at Global Seismographic Network sites are not only acceptable for performing corrections to seismic data, but are also capable of recording many magnetic field signals with similar signal‐to‐noise ratios, in the 20–1000 s period band, as the observatory grade magnetometers operated by the U.S. Geological Survey Geomagnetism Program. This approach would densify magnetic field observations and could also contribute to space weather monitoring by supplementing highly calibrated magnetometers with additional sensors.

Alaska

Real-time performance of the PLUM earthquake early warning method during the 2019 M6.4 and M7.1 Ridgecrest, California, Earthquakes

We evaluate the timeliness and accuracy of ground‐motion‐based earthquake early warning (EEW) during the July 2019 M 6.4 and 7.1 Ridgecrest earthquakes. In 2018, we began retrospective and internal real‐time testing of the propagation of local undamped motion (PLUM) method for earthquake warning in California, Oregon, and Washington, with the potential that PLUM might one day be included in the ShakeAlert EEW system. A real‐time version of PLUM was running on one of the ShakeAlert EEW system’s development servers at the time of the 2019 Ridgecrest sequence, allowing us to evaluate the timeliness and accuracy of PLUM’s warnings for the M 6.4 and 7.1 mainshocks in real time with the actual data availability and latencies of the operational ShakeAlert EEW system. The latter is especially important because high‐data latencies during the M 7.1 earthquake degraded ShakeAlert’s performance. PLUM proved to be largely immune to these latencies. In this article, we present a retrospective analysis of PLUM performance and explore three potential regional alerting strategies ranging from spatially large regions (counties), to moderate‐size regions (National Weather Service public forecast zones), to high‐spatial specificity (50 km regular geographic grid). PLUM generated initial shaking forecasts for the two mainshocks 5 and 6 s after their respective origin times, and faster than the ShakeAlert system’s first alerts. PLUM was also able to accurately forecast shaking across southern California for all three alerting strategies studied. As would be expected, a cost‐benefit analysis of each approach illustrates trade‐offs between increasing warning time and minimizing the area receiving unneeded alerts. Choosing an optimal alerting strategy requires knowledge of users’ false alarm tolerance and minimum required warning time for taking protective action, as well as the time required to distribute alerts to users.

California

Optimizing earthquake early warning alert distance strategies using the July 2019 Mw6.4 and Mw7.1 Ridgecrest, California, earthquakes

The ShakeAlert earthquake early warning system aims to alert people who experience modified Mercalli intensity (MMI) IV+ shaking during an earthquake using source estimates (magnitude and location) to estimate median‐expected peak ground motions with distance, then using these ground motions to determine median‐expected MMI and thus the extent of MMI IV shaking. Because median ground motions are used, even if magnitude and location are correct, there will be people outside the alert region who experience MMI IV shaking but do not receive an alert (missed alerts). We use 91,000 “Did You Feel It?” survey responses to the July 2019 Mw 6.4 and Mw 7.1 Ridgecrest, California, earthquakes to determine which ground‐motion to intensity conversion equation (GMICE) best fits median MMI with distance. We then explore how incorporating uncertainty from the ground‐motion prediction equation and the GMICE in the alert distance calculation can produce more accurate MMI IV alert regions for a desired alerting strategy (e.g., aiming to alert 95% of people who experience MMI IV+ shaking), assuming accurate source characterization. Without incorporating ground‐motion uncertainties, we find MMI IV alert regions using median‐expected ground motions alert fewer than 20% of the population that experiences MMI IV+ shaking. In contrast, we find >94% of the people who experience MMI IV+ shaking can be included in the MMI IV alert region when two standard deviations of ground‐motion uncertainty are included in the alert distance computation. The optimal alerting strategy depends on the false alert tolerance of the community due to the trade‐off between minimizing missed and false alerts. This is especially the case for situations like the Mw 6.4 earthquake when alerting 95% of the 5 million people who experience MMI IV+ also results in alerting 14 million people who experience shaking below this level and do not need to take protective action.

California

Liquefaction and related ground failure from July 2019 Ridgecrest earthquake sequence

The 2019 Ridgecrest earthquake sequence produced a 4 July M 6.5 foreshock and a 5 July M 7.1 mainshock, along with 23 events with magnitudes greater than 4.5 in the 24 hr period following the mainshock. The epicenters of the two principal events were located in the Indian Wells Valley, northwest of Searles Valley near the towns of Ridgecrest, Trona, and Argus. We describe observed liquefaction manifestations including sand boils, fissures, and lateral spreading features, as well as proximate non‐ground failure zones that resulted from the sequence. Expanding upon results initially presented in a report of the Geotechnical Extreme Events Reconnaissance Association, we synthesize results of field mapping, aerial imagery, and inferences of ground deformations from Synthetic Aperture Radar‐based damage proxy maps (DPMs). We document incidents of liquefaction, settlement, and lateral spreading in the Naval Air Weapons Station China Lake US military base and compare locations of these observations to pre‐ and postevent mapping of liquefaction hazards. We describe liquefaction and ground‐failure features in Trona and Argus, which produced lateral deformations and impacts on several single‐story masonry and wood frame buildings. Detailed maps showing zones with and without ground failure are provided for these towns, along with mapped ground deformations along transects. Finally, we describe incidents of massive liquefaction with related ground failures and proximate areas of similar geologic origin without ground failure in the Searles Lakebed. Observations in this region are consistent with surface change predicted by the DPM. In the same region, geospatial liquefaction hazard maps are effective at identifying broad percentages of land with liquefaction‐related damage. We anticipate that data presented in this article will be useful for future liquefaction susceptibility, triggering, and consequence studies being undertaken as part of the Next Generation Liquefaction project.

California

Pseudo-prospective evaluation of UCERF3-ETAS forecasts during the 2019 Ridgecrest sequence

The 2019 Ridgecrest sequence provides the first opportunity to evaluate Uniform California Earthquake Rupture Forecast v.3 with epidemic‐type aftershock sequences (UCERF3‐ETAS) in a pseudoprospective sense. For comparison, we include a version of the model without explicit faults more closely mimicking traditional ETAS models (UCERF3‐NoFaults). We evaluate the forecasts with new metrics developed within the Collaboratory for the Study of Earthquake Predictability (CSEP). The metrics consider synthetic catalogs simulated by the models rather than synoptic probability maps, thereby relaxing the Poisson assumption of previous CSEP tests. Our approach compares statistics from the synthetic catalogs directly against observations, providing a flexible approach that can account for dependencies and uncertainties encoded in the models. We find that, to the first order, both UCERF3‐ETAS and UCERF3‐NoFaults approximately capture the spatiotemporal evolution of the Ridgecrest sequence, adding to the growing body of evidence that ETAS models can be informative forecasting tools. However, we also find that both models mildly overpredict the seismicity rate, on average, aggregated over the evaluation period. More severe testing indicates the overpredictions occur too often for observations to be statistically indistinguishable from the model. Magnitude tests indicate that the models do not include enough variability in forecasted magnitude‐number distributions to match the data. Spatial tests highlight discrepancies between the forecasts and observations, but the greatest differences between the two models appear when aftershocks occur on modeled UCERF3‐ETAS faults. Therefore, any predictability associated with embedding earthquake triggering on the (modeled) fault network may only crystalize during the presumably rare sequences with aftershocks on these faults. Accounting for uncertainty in the model parameters could improve test results during future experiments.

California

The predictive skills of elastic Coulomb rate-and-state aftershock forecasts during the 2019 Ridgecrest, California, earthquake sequence

Operational earthquake forecasting protocols commonly use statistical models for their recognized ease of implementation and robustness in describing the short-term spatiotemporal patterns of triggered seismicity. However, recent advances on physics-based aftershock forecasting reveal comparable performance to the standard statistical counterparts with significantly improved predictive skills when fault and stress field heterogeneities are considered. Here, we perform a pseudo-prospective forecasting experiment during the first month of the 2019 Ridgecrest (California) earthquake sequence. We develop seven Coulomb rate-and-state models that couple static stress change estimates with continuum mechanics expressed by the rate-and-state friction laws. Our model parametrization supports a gradually increasing complexity; we start from a preliminary model implementation with simplified slip distributions and spatially homogeneous receiver faults to reach an enhanced one featuring optimized fault constitutive parameters, finite-fault slip models, secondary triggering effects, and spatially heterogenous planes informed by pre-existing ruptures. The data-rich environment of Southern California allows us to test whether incorporating data collected in near real-time during an unfolding earthquake sequence boosts our predictive power. We assess the absolute and relative performance of the forecasts by means of statistical tests used within the Collaboratory for the Study of Earthquake Predictability (CSEP) and compare their skills against a standard benchmark ETAS model for the short (24 hours after the two Ridgecrest mainshocks) and intermediate-term (one month). Stress-based forecasts expect heightened rates along the whole near-fault region and increased expected seismicity rates in Central Garlock Fault. Our comparative model evaluation supports that faulting heterogeneities coupled with secondary triggering effects are the most critical success components behind physics-based forecasts, but also underlines the importance of model updates incorporating near real-time available aftershock data reaching better performance than ETAS models. We explore the physical basis behind our results by investigating the localized shut down of pre-existing normal faults in the Ridgecrest near-source area.

California

Mild displacements of boulders during the 2019 Ridgecrest Earthquakes

Strong seismic waves from the July 2019 Ridgecrest, California, earthquakes displaced rocks in proximity to the M 7.1 mainshock fault trace at several locations. In this report, we document large boulders that were displaced at the Wagon Wheel Staging Area (WWSA), approximately 4.5 km southeast of the southern terminus of the large M 6.4 foreshock rupture (hereafter “the large foreshock”) and 9 km southwest of the nearest approach of the M 7.1 mainshock surface rupture. Some boulders appear to have slid along essentially flat surfaces, which suggest that dynamic stresses overcame the coefficient of friction. Other boulders appear to have rocked within their sockets. In both cases, we use simple mechanical models to estimate total peak dynamic accelerations between 0.5 and 1g⁠, commensurate with modified Mercalli intensity 9. It is unclear if the strongest shaking at this location occurred during the large foreshock or the M 7.1 mainshock. The inferred accelerations are higher than predicted mainshock ground motions at WWSA, although local high accelerations could have been generated by path, site, or source effects. Gaps between boulders and their sockets are easily visible in the immediate aftermath of earthquakes and provide a quick indication of strong shaking. More importantly, the gaps quickly fill with surficial organic debris, including seeds and leaves of the year, that quickly become entombed. Boulders may thus potentially be extracted to examine gap fillings associated with past earthquakes, providing a new datable paleoseismic method.

California

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

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

California

ShakeAlert Earthquake Early Warning System Performance During the 2019 Ridgecrest Earthquake Sequence

During July 2019, a sequence of earthquakes including a Mw6.4 foreshock and a Mw7.1 mainshock occurred near Ridgecrest, California. ShakeAlert, the U.S. Geological Survey (USGS) ShakeAlert public Earthquake Early Warning (EEW) system being developed for the U.S. West Coast, was operational during this time, though public alerting was only available within LA County. ShakeAlert created alert messages for many of the earthquakes, including the two largest events, and for many of the larger aftershocks. In this study, we dissect log files and replay data through the system to reconstruct the sequence of events and analyze the performance of the system during that time period. While the system performed reasonably well overall, the sequence also revealed various issues and short comings that will be addressed in impending and future system upgrades, with most parts of the system working as they should. ShakeAlert correctly detected and rapidly characterized both the Mw6.4 and Mw7.1 earthquakes within 6.9 s of their origin times and created alert messages that were available to ShakeAlert’s pilot users. No public alerts were sent out by the ShakeAlertLA cellphone app (the only publicly available alerting method at the time) because the predicted shaking for LA County was below the app’s alerting threshold of MMI 4.0. For the Mw6.4 event this was accurate. For the Mw7.1 event, public alerts for LA County were warranted, but ShakeAlert underpredicted the shaking levels because both the point-source and the finite-fault algorithms underestimated the magnitude of the earthquake by 0.8 units. A number of software and hardware issues that were responsible for the magnitude underestimation of the mainshock have been identified and will be addressed in future ShakeAlert releases. We also analyze the hypothetical alerting performance of ShakeAlert had public alerting been available throughout southern California with a lower alerting threshold of 〖MMI〗_alert=2.5MMI 2.5. We find that, despite the magnitude underestimation, ShakeAlert could have provided timely warnings to a large fraction of affected sites, including some of the near-epicentral sites with high ground motion intensities.

California

Evidence of previous faulting along the 2019 Ridgecrest, California earthquake ruptures

The July 2019 Ridgecrest earthquake sequence in southeastern California was characterized as surprising because only ~35% of the rupture occurred on previously mapped faults. Employing more detailed inspection of pre-event high-resolution topography and imagery in combination with field observations, we document evidence of active faulting in the landscape along the entire fault system. Scarps, deflected drainages, and lineaments and contrasts in topography, vegetation, and ground color demonstrate previous slip on a dense network of orthogonal faults, consistent with patterns of surface rupture observed in 2019. Not all of these newly mapped fault strands ruptured in 2019. Outcrop-scale field observations additionally reveal tufa lineaments and sheared Quaternary deposits. Neotectonic features are commonly short (<2 km), discontinuous, and display en echelon patterns along both the M 6.4 and M 7.1 ruptures. These features are generally more prominent and better preserved outside the late Pleistocene lake basins. Fault expression may also be related to deformation style: scarps and topographic lineaments are more prevalent in areas where substantial vertical motion occurred in 2019. Where strike-slip displacement dominated in 2019, the faults are mainly expressed by less prominent tonal and vegetation features. Both the NE- and NW-trending active fault systems are subparallel to regional bedrock fabrics that were established as early as ~150 Ma, and may be reactivating these older structures. Overall, we estimate that 50-70% (i.e., an additional 15-35%) of the 2019 surface ruptures could have been recognized as active faults with detailed inspection of pre-event data. Similar detailed mapping of potential neotectonic features could help improve seismic hazard analyses in other regions of eastern California and elsewhere that have distributed faulting or incompletely mapped faults. In areas where faults cannot be resolved as single thoroughgoing structures, a zone of potential faulting should be used as a hazard model input.

California

Surface displacement distributions for the July 2019 Ridgecrest, California earthquake ruptures

Surface rupture in the 2019 Ridgecrest, California, earthquake sequence occurred along two orthogonal cross faults and includes dominantly left‐lateral and northeast‐striking rupture in the M w "> M w 6.4 foreshock and dominantly right‐lateral and northwest‐striking rupture in the M w "> M w 7.1 mainshock. We present &gt; 650 "> > 650 field‐based, surface‐displacement observations for these ruptures and synthesize our results into cumulative along‐strike displacement distributions. Using these data, we calculate displacement gradients and compare our results with historical strike‐slip ruptures in the eastern California shear zone. For the M w "> M w 6.4 rupture, we report 96 displacements measured along 18 km of northeast‐striking rupture. Cumulative displacement curves for the rupture yield a mean left‐lateral displacement of 0.3–0.5 m and maximum of 0.7–1.6 m. Net mean vertical displacement based on the difference of down‐to‐the‐west (DTW) and down‐to‐the‐east (DTE) displacement curves is close to zero (0.02 m DTW). The M w "> M w 6.4 displacement distribution shows that the majority of displacement occurred southwest of the intersection with the M w "> M w 7.1 rupture. The M w "> M w 7.1 rupture is northwest‐striking and 50 km long based on 576 field measurements. Displacement curves indicate a mean right‐lateral displacement of 1.2–1.7 m and a maximum of 4.3–7.0 m. Net vertical displacement in the rupture averages 0.3 m DTW. The M w "> M w 7.1 displacement distributions demonstrate that maximum displacement occurred along a 12‐km‐long portion of the fault near the M w "> M w 7.1 epicenter, releasing 66% of the geologically based seismic moment along 24% of the total rupture length. Using our displacement distributions, we calculate kilometer‐scale displacement gradients for the M w "> M w 7.1 rupture. The steepest gradients ( ⁠ &#x223C; 1 &#x2013; 3 &#x2009;&#x2009; m / km "> ∼ 1 – 3 m / km ⁠ ) flank the 12‐km‐long region of maximum displacement. In contrast, gradients for the 1992 M w "> M w 7.3 Landers and 1999 M w "> M w 7.1 Hector Mine earthquakes are &lt; 0.6 &#x2009;&#x2009; m / km "> < 0.6 m / km . Our displacement distributions are important for understanding the influence of cross‐fault rupture on M w "> M w 6.4 and 7.1 rupture length and displacement and will facilitate comparisons with distributions generated remotely and at broader scales.

California

Minimal clustering of injection-induced earthquakes observed with a large-n seismic array

The clustering behavior of injection‐induced earthquakes is examined using one month of data recorded by the LArge‐ n Seismic Survey in Oklahoma (LASSO) array. The 1829‐node seismic array was deployed in a 25 km×32 km area of active saltwater disposal in northern Oklahoma between 14 April and 10 May 2016. Injection rates in the study area are nearly constant around the time of the deployment. We develop a local magnitude ( ⁠ M L ⁠ ) equation for the region and estimate magnitudes for 1104 earthquakes recorded by the deployment. The determined earthquake magnitudes range from M L 0.01 to 3.0. The majority of earthquakes occurred between 1.5 and 5.5 km depth, and the shallowest earthquake depths overlap with the base of injection wells at depths between 1.5 and 2.5 km. We compute focal mechanisms of the largest events ( ⁠ M L > 2.0 ⁠ ), and find a mix of normal‐ and strike‐slip‐faulting types. Earthquakes occur regularly in time during the deployment, but are not evenly distributed in space across the study area, that is, they are spatially clustered. Analysis of the nearest‐neighbor distances in the space–time–magnitude domain shows the seismicity is dominated by single‐event clusters (i.e., independent events). This high proportion of single‐event clusters compared with multievent clusters has been previously noted for induced events at geothermal sites. When clustering occurs, the number of events in a cluster is typically small. We observe only four clusters with 10 or more events. For these larger clusters, we find equivalent numbers of foreshocks and aftershocks; however, the foreshock sequences are significantly longer in duration lasting days to tens of days, while aftershock sequences are observed only on the order of one day. The minimal clustering observed for events in the LASSO array suggests that the majority of events are being directly driven by stress changes due to local saltwater disposal.

Oklahoma

Factors influencing the probability of hydraulic fracturing induced seismicity in Oklahoma

Injection‐induced seismicity became an important issue over the past decade, and although much of the rise in seismicity is attributed to wastewater disposal, a growing number of cases have identified hydraulic fracturing (HF) as the cause. A recent study identified regions in Oklahoma where ≥75% of seismicity from 2010 to 2016 correlated with nearly 300 HF wells. To identify factors associated with increased probability of induced seismicity, we gathered publicly available information about the HF operations in these regions including: injected volume, number of wells on a pad, injected fluid (gel vs. slickwater), vertical depth of the well, proximity of the well to basement rock, and the formation into which the injection occurred. To determine the statistical strength of the trends, we applied logistic regression, bootstrapping, and odds ratios. We see no trend with total injected volume in our Oklahoma dataset, in contrast to strong trends observed in Alberta and Texas, but we note those regions have many more multiwell pads leading to larger cumulative volumes within a localized area. We found a > ∼50% lower probability of seismicity with the use of gel compared to slickwater. We found that HF wells targeting older formations had a higher probability of seismicity; however, these wells also tend to be deeper, and we found the trend with well depth to be stronger than the trend with age of formation. When isolated to the Woodford formation, well depth produced the strongest relationship, increasing from ∼5%> to ∼50% probability from 1.5 to 5.5 km. However, no trend was seen in the proximity to basement parameter. Based on previously measured pore pressure gradients, we interpret the strong absolute depth relationship to be a result of the increasing formation overpressure measured in deeper portions of the basin that lower the stress change needed to induce seismicity.

Oklahoma

More fault connectivity Is needed in seismic hazard analysis

Did the third Uniform California Earthquake Rupture Forecast (UCERF3) go overboard with multifault ruptures? Schwartz (2018) argues that there are too many long ruptures in the model. Here, I address his concern and show that the UCERF3 rupture‐length distribution matches empirical data. I also present evidence that, if anything, the UCERF3 model could be improved by adding more connectivity to the fault system. Adding more connectivity would improve model misfits with data, particularly with paleoseismic data on the southern San Andreas fault; make the model less characteristic on the faults; potentially improve aftershock forecasts; and reduce model sensitivity to inadequacies and unknowns in the modeled fault system.

Bulletin of the Seismological Society of America

VS30 and Dominant Site Frequency (⁠fd⁠) as Provisional Station ML Corrections (⁠dML⁠) in California

New seismic stations added to a regional seismic network cannot be used to calculate local magnitude ( ⁠ M L "> M L ⁠ ) until a revised regionwide amplitude decay function is developed. Each station must record a minimum number of local and regional earthquakes that meet specific amplitude requirements prior to recalibration of the amplitude decay function. Station component adjustments ( ⁠ d M L "> d M L ⁠ ; Uhrhammer et al. , 2011 ) are then calculated after inverting for a new regional amplitude decay function, constrained by the sum of d M L "> d M L ⁠ for long‐running stations. Therefore, there can be significant delay between when a new station starts contributing real‐time waveform packets and when data can be included in magnitude determinations. We propose the use of known estimates of seismic site conditions such as the time‐averaged shear‐wave velocity ( ⁠ V S "> V S ⁠ ) of the upper 30 m ( ⁠ V S 30 "> V S 30 ⁠ ) and the site dominant frequency ( ⁠ f d "> f d ⁠ ) to calculate d M L "> d M L ⁠ ⁠ . Previously established d M L "> d M L ⁠ ⁠ , measured V S 30 "> V S 30 ⁠ ⁠ , and f d "> f d data are available for between 126 and 458 horizontal components (east–west and north–south) at 81 seismic stations in the California Integrated Seismic Network; d M L "> d M L ⁠ d M L "> data range from −1.10 to 0.39, V S 30 "> V S 30 ⁠ values range from 202 to 1464 &#x2009;&#x2009; m / s "> 1464 m / s ⁠ , and 440 f d "> f d values are compiled from earthquake and microseismic records that range from 0.13 to 21 Hz. We find V S 30 "> V S 30 ⁠ and d M L "> d M L ⁠ exhibit a positive coefficient of determination ( ⁠ R = 0.59 "> R = 0.59 ⁠ ), indicating that as V S 30 "> V S 30 ⁠ increases, d M L "> d M L ⁠ increases. This implies that greater site amplification (lower V S 30 "> V S 30 ⁠ ⁠ ) results in smaller d M L "> d M L ⁠ ⁠ . f d "> f d and d M L "> d M L ⁠ also generally exhibit a positive correlation ( ⁠ R 2 &lt; 0.56 "> R 2 < 0.56 ⁠ ), which implies lower d M L "> d M L values are related to site resonance at depth‐dependent frequencies. Using the developed relationships, V S 30 "> V S 30 ⁠ or f d "> f d measurements can be used to establish a provisional d M L "> d M L ⁠ for newly established stations. This procedure allows new stations to contribute to regional network M L "> M L determinations immediately without the need to wait until a minimum set of earthquake data has been recorded.

California

Strain-estimated ground motions associated with recent earthquakes in California

Peak ground velocity (PGV) is a commonly used parameter in earthquake ground‐motion models (GMMs) and hazard analyses, because it is closely related to structural damage and felt ground shaking, and is typically measured on broadband seismometers. Here, we demonstrate that strainmeters, which directly measure in situ strain in the bulk rock, can easily be related to ground velocity by a factor of bulk shear‐wave velocity and, thus, can be used to measure strain‐estimated PGV. We demonstrate the parity of velocity to strain utilizing data from borehole strainmeters deployed along the plate boundaries of the west coast of the United States for nine recent M "> M M 4.4–7.1 earthquakes in California, including the largest two events of the July 2019 Ridgecrest earthquake sequence. PGVs derived from maximum horizontal shear strains fall within the range of seismic‐estimated values recorded at the same distances. We compare the strain‐estimated data with GMMs based on seismic PGVs and find consistency in residual polarity (positive vs. negative; the sign of the difference between observed and modeled data) for certain earthquake–station paths, where some paths indicate an overestimation and others indicate an underestimation of strain‐derived PGVs, as compared with the GMMs. We surmise that this may be indicative of over or underestimation of shear‐wave velocity along those paths, as compared with the average velocity used to derive PGV from strain measurements, or indicative of repeatable site and path effects that are not accounted for in our analyses. This direct comparison of strain with velocity can highlight physical path effects, as well as improve the density and capability of ground‐motion recordings.

California

San Andreas fault exploration using refraction tomography and S-wave-type and Fϕ-mode guided waves

Surface ruptures from the 18 April 1906 M∼7.9 San Francisco earthquake were distributed over an ∼35 ‐meter‐wide zone at San Andreas Lake on the San Francisco Peninsula in California ( Schussler, 1906 ). Since ∼1906⁠ , the surface ruptures have been largely covered by water, but with water levels at near‐historic low levels in 2008–2011, we observed that the 1906 surface ruptures were no longer visible. As a fault imaging test, we acquired refraction tomography and guided‐wave data across the 1906 surface ruptures in 2011. We found that individual fault traces, as mapped by Schussler (1906) , can be identified on the basis of discrete low‐velocity zones ( ⁠V S and V P⁠ , reduced ∼40% and ∼34%⁠ , respectively) and high‐amplitude guided waves. Guided waves have traditionally been observed as large‐amplitude waveforms over wide (hundreds of meters to kilometers) zones of faulting, but we demonstrate that by evaluating guided waves (including Rayleigh/Love‐ and P / SV ‐types) in terms of peak ground velocity (PGV), individual near‐surface fault traces within a fault zone can be precisely located, even more than 100 yr after the surface ruptures. Such precise exploration can be used to focus paleoseismic trenching efforts and to identify or exclude faulting at specific sites. We evaluated PGV of both S ‐wave‐type and Fϕ ‐mode‐type guided waves and found that both wave types can be used to identify subsurface fault traces. At San Andreas Lake (main fault), S ‐wave‐type guided waves travel up to 18% slower than S body waves, and Fϕ ‐mode guided waves travel ∼60% slower than P body waves but ∼15% faster than S body waves. We found that guided‐wave amplitudes vary with frequency but are up to five times higher than those of body waves, including the S wave. Our data are consistent with the concept that guided waves can be a strong‐shaking hazard during large‐magnitude earthquakes.

Calfornia

Human behavioral response in the Ridgecrest earthquakes: Assessing immediate actions based on data from “Did You Feel It?”

Human behavioral response to earthquake ground motion has long been a subject of multidisciplinary interest and research. In most versions of seismic intensity scales, human perceptions and behavior are one component of the assignment of intensity. Public health research has shown that actions taken during earthquakes have a significant impact on the incidence of injury or the maintenance of safety. Based on this research, emergency managers and organizations promoting emergency preparedness have advocated strategies such as drop, cover, and hold on (DCHO) and promoted this safety measure through public education and annual drills. The “Did You Feel It?” (DYFI) mapping system (see Data and Resources) based on an online questionnaire developed and maintained by the U.S. Geological Survey has provided opportunities for those who have experienced an earthquake to report this experience worldwide since 2004. The DYFI questionnaire, although designed to assign intensity, also contains questions regarding the behavior in which one has engaged during the earthquake. The questionnaire includes other important information that may elucidate behavioral response to earthquakes, including assigned intensity, emotional reaction, and whether damage occurred at the location where the earthquake was experienced. The very large number of people who completed DYFI questionnaires following the July 2019 Ridgecrest, California, earthquakes provides a robust dataset for analysis and suggests that as intensity and levels of fear increase, behavior becomes more active in terms of physical movement to locations of presumed safety. Among active responses including DCHO, going to a doorway, and running outside, DCHO was the least likely to be implemented. The study provides possible explanations for low participation in DCHO despite active campaigns to promote this strategy.

California