ABSTRACT Fragile geologic features (FGFs) can be used to estimate an upper bound to the strength of earthquake ground shaking in the past thousands to tens of thousands of years. FGFs could be especially useful in the eastern United States, where few active faults are available for paleoseismic studies and where low seismic attenuation means that an earthquake could damage FGFs at large distances. This article describes a pilot study of one form of FGF, precariously balanced rocks (PBRs), in the Blue Ridge of Virginia, ∼160 km from the central Virginia seismic zone (CVSZ) and ∼100 km from Virginia’s Giles County seismic zone (GCSZ). The main goal is to develop and refine methods for studying PBRs in the eastern United States. Photogrammetric 3D models of the rocks were used to locate their centers of mass and delineate their basal contacts, with the latter corroborated by wrapping a wire around the bases to verify the circumferences. Relations between hill geometry and amplification were used to estimate topographic amplification. Seismometers on the PBRs provided the frequencies, amplitudes, and directions of rocking after giving the rocks gentle pushes. The levels of ground shaking likely to topple the rocks were compared with the median hazard curves from the U.S. Geological Survey’s National Seismic Hazard Model and were used to map magnitudes of earthquakes in the nearby region likely to topple the rocks. Results indicate that an Mw≥7.5 earthquake is needed in the CVSZ or near the 2020 Sparta, North Carolina, earthquake to topple the two most fragile PBRs and that an earthquake of ∼Mw 7–7.5 in the GCSZ is needed. The fragilities and estimated ages of the PBRs (20–40 ka) are thus consistent with the absence of such events in the historical record in the nearby region.

Semantic mapping in 3D is fundamental to a wide range of geoscientific studies and applications, including geomorphology, hazard assessment, and environmental monitoring. However, automatically segmenting geological features from large-scale photogrammetric datasets remains a significant challenge. We present a methodology to address this gap. Using overlapping images collected over environments of interest, Structure-from-Motion (SfM) produces georeferenced point clouds and estimates camera poses. Existing large vision models, such as Segment Anything Model, segment objects in the images, generating pixel-segmentation associations. To produce pixel-point associations, we project the points back onto the camera image planes. As objects are independently segmented across multiple images with different perspectives, we develop a segmentation mosaicking algorithm to build probabilistic point-segmentation associations that combines the pixel-segmentation associations and pixel-point associations. Our methodology is validated using both synthetic data generated by Kubric and real-world UAV-SfM data. The implementation is designed to be compatible with existing SfM software, including Agisoft and OpenDroneMap, for photogrammetry mapping in geoscience studies. As a case study, we apply our method to the semantic mapping of precariously balanced rocks (PBRs), which provide upper-bound constraints on historical ground motion shaking intensity. To support object-level identification of PBRs, we additionally integrated Grounding DINO, enabling text-prompted segmentation of features of interest within UAV imagery. This case study demonstrates the effectiveness of our method in generating a 3D semantic map of PBRs, enabling spatial distribution of PBR fragility for earthquake hazard analysis.

AbstractStrong earthquakes along subduction zones are often devastating events, but sparse records along some tectonic margins limit our understanding of seismic hazards. Constraining shaking intensities is critical, especially in subduction zones with infrequent but large‐magnitude earthquakes like the Cascadia Subduction Zone (CSZ), where the lack of recorded ground motions has led to uncertainty in the severity and potential impacts of future earthquakes. Here we fill this observational gap with a novel inventory of quantitative estimates of past shaking intensities from geotechnical modeling of coastal landforms. One hundred fifty‐four deep‐seated landslides and 65 fragile geologic features constrain minimum and maximum peak ground accelerations, respectively. These estimates are broadly consistent with model predictions of M9 ruptures, suggesting strong shaking of 0.4–0.8 g during past CSZ earthquakes. Local discrepancies between our geologic shaking constraints and earthquake simulations may inform past rupture behavior, leading to better predictions of shaking intensity for future earthquakes.

We analyze the seismic response of a class of fragile geologic features (FGFs), referred to as rock towers (RTs) at the Trona Pinnacles, a group of RTs in southern California that suffered strong shaking during the 2019 Ridgecrest earthquake sequence. FGFs, including RTs, may provide maximum constraints on past earthquake shaking intensity, and thereby support probabilistic seismic hazard assessments (PSHAs). In a rare case study, we explore the hypothesis that RT structural integrity is time dependent, as damage accumulates progressively. We develop finite-element method (FEM) models of the RTs using photogrammetric shape models. We validate the models by comparing numerical simulations of their response to broadband ground shaking with low-intensity seismic recordings obtained at the Pinnacles. Results of our simulations are in good agreement with the seismic recordings of actual earthquake aftershocks. We next use the results of the FEM models to analyze the response and evolution of RTs. Our analyses elucidate the influence of geometry over their seismic response, providing a rationale that may explain the rarity of slender RTs at Trona: high-aspect-ratio structures that respond in bending develop detrimental tensile stresses that crack the rock, whereas low-aspect-ratio ones' response also includes shearing, which does not compromise material integrity as much as tension. Field measurements with a rebound hammer support this finding, suggesting that the material around the base of slender rocks has been weakened relative to other parts of the RT. We also study how to define simplified mechanical models ("archetypes") to predict the natural frequencies of RTs. Results from our work illuminate the fundamental mechanisms of seismic response and progressive failure of RTs, and open new avenues of research to potentially incorporate these geologic features as long-return period constraints on PSHA, in ways analogous to those of the widely used precariously balanced rocks.

We articulate a scientific vision and roadmap for the development of improved Earthquake Rupture Forecast (ERF) models, which are one of the two main modeling components used in modern seismic hazard and risk analysis. One primary future objective is to provide fully time-dependent models that include both elastic rebound and spatiotemporal clustering nationwide, which can be particularly important for shorter-term hazard and risk considerations (e.g., earthquake insurance products). We also discuss the importance and perennial challenges associated with quantifying epistemic uncertainties, including those associated with deformation-model slip rates, unquantified sampling errors with respect to off-fault seismicity, and any spatial covariances. The need for more physics-based approaches is also emphasized, as is the benefit of adding model valuation (quantifying usefulness) to our verification and validation protocols. Given the multidisciplinary and system-level nature of this activity, modular design is critical. Future updates will also draw from best-available science by both the U.S. Geological Survey and the external community. The primary goal of this article is to highlight plans that guide research and facilitate community engagement with model development, especially with respect to lowering the entry barrier for early-career scientists and engineers. The article is written so readers can focus on the sections that interest them most (see Manuscript Organization section), with the Introduction and the Discussion section providing a stand-alone overview and summary, respectively.

We present the 2023 U.S. Geological Survey time -independent earthquake rupture forecast for the conterminous United States, which gives authoritative estimates of the magnitude, location, and time -averaged frequency of potentially damaging earthquakes throughout the region. In addition to updating virtually all model components, a major focus has been to provide a better representation of epistemic uncertainties. For example, we have improved the representation of multifault ruptures, both in terms of allowing more and less fault connectivity than in the previous models, and in sweeping over a broader range of viable models. An unprecedented level of diagnostic information has been provided for assessing the model, and the development was overseen by a 19 -member participatory review panel. Although we believe the new model embodies significant improvements and represents the best available science, we also discuss potential model limitations, including the applicability of logic tree branch weights with respect different types of hazard and risk metrics. Future improvements are also discussed, with deformation model enhancements being particularly worthy of pursuit, as well as better representation of sampling errors in the gridded seismicity components. We also plan to add time -dependent components, and assess implications with a wider range of hazard and risk metrics.

ABSTRACT Precariously balanced rocks (PBRs) and other fragile geologic features have the potential to constrain the maximum intensity of earthquake ground shaking over millennia. Such constraints may be particularly useful in the eastern United States (U.S.), where few earthquake-source faults are reliably identified, and moderate earthquakes can be felt at great distances due to low seismic attenuation. We describe five PBRs in northern New York and Vermont—a region of elevated seismic hazard associated with historical seismicity. These boulders appear to be among the most fragile PBRs in the region, based on reports from hobbyists. The PBRs are glacial erratics, best evidenced by glacial striations on bedrock pedestals. The pedestals themselves are locally high knobs, often situated on regionally high topography; this setting limits soil development and indicates that any outwash deposits were likely ephemeral. As a result, PBR ages can be reliably established by the retreat of the last continental ice sheet, ∼15–13 ka. To quantify the fragility of the PBRs, we surveyed them with ground-based light detection and ranging and calculated geometric parameters from the point clouds, field observations, and seismic responses. Preliminary validation of the 2023 time-independent U.S. National Seismic Hazard Model (NSHM) shows that the existence of PBRs is generally consistent with the median site-specific hazard curves. Only the Blue Ridge Road site suggests a modest reduction in hazard. To visualize the ensemble of data, we mapped the minimum permissible distance to potential source faults around each PBR site as a function of source magnitude by using the ground-motion models from the 2023 NSHM. Viewed in this manner, our data are consistent with potential M∼6.5 earthquake-source faults in many parts of the Lake Champlain Valley and northern Adirondack Mountains. Our work illustrates a potential pathway for better constraining earthquake-source faults in regions of cryptic faults.

The US National Seismic Hazard Model (NSHM) was updated in 2023 for all 50 states using new science on seismicity, fault ruptures, ground motions, and probabilistic techniques to produce a standard of practice for public policy and other engineering applications (defined for return periods greater than ∼475 or less than ∼10,000 years). Changes in 2023 time-independent seismic hazard (both increases and decreases compared to previous NSHMs) are substantial because the new model considers more data and updated earthquake rupture forecasts and ground-motion components. In developing the 2023 model, we tried to apply best available or applicable science based on advice of co-authors, more than 50 reviewers, and hundreds of hazard scientists and end-users, who attended public workshops and provided technical inputs. The hazard assessment incorporates new catalogs, declustering algorithms, gridded seismicity models, magnitude-scaling equations, fault-based structural and deformation models, multi-fault earthquake rupture forecast models, semi-empirical and simulation-based ground-motion models, and site amplification models conditioned on shear-wave velocities of the upper 30 m of soil and deeper sedimentary basin structures. Seismic hazard calculations yield hazard curves at hundreds of thousands of sites, ground-motion maps, uniform-hazard response spectra, and disaggregations developed for pseudo-spectral accelerations at 21 oscillator periods and two peak parameters, Modified Mercalli Intensity, and 8 site classes required by building codes and other public policy applications. Tests show the new model is consistent with past ShakeMap intensity observations. Sensitivity and uncertainty assessments ensure resulting ground motions are compatible with known hazard information and highlight the range and causes of variability in ground motions. We produce several impact products including building seismic design criteria, intensity maps, planning scenarios, and engineering risk assessments showing the potential physical and social impacts. These applications provide a basis for assessing, planning, and mitigating the effects of future earthquakes.

Radiocarbon dates from alluvial sections provide maximum deposit ages because of the time lag between formation of the dated material and deposition at the sample site, potentially producing decade- to century-long biases in the dates of historic events, paleoclimatic change, fire histories, and paleoearthquakes. This bias, called the inherited age, combines the inbuilt age distribution, which reflects the age composition of the vegetation of the source area, and the residence time distribution, which includes transport and interim storage prior to final deposition. We tackle inherited age and its components by comparing charcoal dates from two modern fires in southern California, the 2020 Bobcat Fire and the 2013 Grand Fire, with a well-dated late Holocene terrace deposit in the Pallett Creek watershed. Fifty-six radiocarbon dates from the modern fires provide an inbuilt age distribution with a median of 25 years pre-fire (320-year 95% range). An inherited age distribution calculated from 175 terrace deposit dates is older, with a median age of similar to 90 years (850-year 95% range). Comparing inherited ages calculated from organic-rich versus clastic terrace deposits reveals a slight facies dependence suggesting longer residence times in clastic deposits. We develop a modeled inherited age that incorporates larger calibration uncertainties in pre-1950s samples by combining the modern fire sample distribution with the pre-bomb portion of the calibration curve. The modeled inherited age is younger than the terrace deposit inherited age by only 21 years, indicating inbuilt age, not long residence times, dominates inherited age in this setting. The results imply that paleoearthquakes and climatic event age estimates in the Western Transverse Ranges are up to a century too old. More broadly, dating charcoal from modern fires can constrain inherited age and the resulting distributions can improve the accuracy of dates of past environmental and tectonic events.

To constrain fault processes and hazard, fault slip rates may be extrapolated over different fault lengths or time intervals. Here, we investigate slip rates for the Cucamonga Fault (CF). The CF is located at the junction of the Transverse Range fault system with the San Andreas and San Jacinto Faults, and it is hypothesized to connect with these faults, promoting the propagation of large, multi‐fault earthquakes. Previous work has shown that CF displacements on late Quaternary alluvial fan surfaces are highly variable along strike. We present two new 10 Be surface exposure ages from depth profiles on the alluvial fans. Slip rates are consistent with a rate of 1.4 ± 0.3 m/kyr over time intervals of ∼20, ∼30, and ∼40 kyr. If the CF participates in multi‐fault ruptures, then these earthquakes occur either rarely or with sufficient regularity to maintain apparently steady rates over multiple intervals. We also explore along‐strike fault displacement variability using a calibrated morphological model. The model successfully reproduces scarp profiles and indicates that fault displacement variability can be explained in part by scarp age but not uplift rate. We infer that both erosion by ephemeral gullying and distributed deformation contribute to fault displacement variability, although both are difficult to detect confidently without excavations across the scarp. These investigations show that better characterization of cumulative‐slip variability along strike may improve accuracy and precision of slip rates. Slip rates that do not consider epistemic uncertainties may not be suitable for extrapolation over longer fault sections.

ABSTRACT We present the 2023 U.S. Geological Survey time-independent earthquake rupture forecast for the conterminous United States, which gives authoritative estimates of the magnitude, location, and time-averaged frequency of potentially damaging earthquakes throughout the region. In addition to updating virtually all model components, a major focus has been to provide a better representation of epistemic uncertainties. For example, we have improved the representation of multifault ruptures, both in terms of allowing more and less fault connectivity than in the previous models, and in sweeping over a broader range of viable models. An unprecedented level of diagnostic information has been provided for assessing the model, and the development was overseen by a 19-member participatory review panel. Although we believe the new model embodies significant improvements and represents the best available science, we also discuss potential model limitations, including the applicability of logic tree branch weights with respect different types of hazard and risk metrics. Future improvements are also discussed, with deformation model enhancements being particularly worthy of pursuit, as well as better representation of sampling errors in the gridded seismicity components. We also plan to add time-dependent components, and assess implications with a wider range of hazard and risk metrics.

The 2019 M-w 7.1 Ridgecrest California earthquake rupture passed within 4 km of the Trona Pinnacles, a large group of tufa rock pillars. Reconnaissance following the Ridgecrest mainshock documented fresh damage to several of the Pinnacles. Repeated aerial photogrammetric surveys also documented damage during subsequent aftershocks. Here, we describe the photogrammetric data with emphasis on a specific rock spire that toppled during an M-w 5.5 aftershock. We calculate the volumes of the intact spire and of its subsequent debris. To explore the utility of the pinnacles as fragile geologic features for constraining past earthquake shaking intensity, we calculate the quasi static, horizontal acceleration required to break the spire at its base. We also examine the response of this feature to observed shaking using a dynamic model of the spire. In this case, we find that the quasi-static estimate provides a conservative maximum constraint on fragility. The dynamic model of the spire suggests that shaking during the M-w 7.1 mainshock likely generated tensile stresses in excess of the spire's bulk strength, thereby making it vulnerable to collapse in subsequent aftershocks.

Recurrence intervals for ground rupturing earthquakes are critical data for assessing seismic hazard. Recurrence intervals are presented here for 38 paleoseismic sites in California. Eleven of these include new or updated data; the remainder use data pre-viously included in the Unified California Earthquake Rupture Forecast Version 3 (UCERF3). The methods and results are consistent with UCERF3. In addition, revised recurrence intervals are presented at every site. The revised recurrence intervals incor-porate uncertainty in the interpretation of paleoseismic evidence, which is expressed as event likelihood. Event likelihood is the probability that the evidence has been correctly interpreted as a unique earthquake. Event likelihoods are estimated here for 85 inferred past earthquakes at eight paleoseismic sites in California, using a single, consistent methodology. The average event likelihood is 0.85. The revised recurrence intervals are 16% longer, on average, than conventional estimates, and their confidence intervals are disproportionately wider. These recurrence intervals are suitable for inclusion in a "grand inversion" rupture forecast, and they may be important for addressing a system-atic misfit in the UCERF3 grand inversion. The revised recurrence intervals may also be important for assessing the unusually long earthquake hiatus in California. Other appli-cations may not need to consider event likelihoods because the effects are small relative to typical uncertainties.

ABSTRACT Fragile geologic features (FGFs), which are extant on the landscape but vulnerable to earthquake ground shaking, may provide geological constraints on the intensity of prior shaking. These empirical constraints are particularly important in regions such as the Pacific Northwest that have not experienced a megathrust earthquake in written history. Here, we describe our field survey of FGFs in southern Oregon. We documented 58 features with fragile geometric characteristics, as determined from field measurements of size and strength, historical photographs, and light detection and ranging point clouds. Among the surveyed FGFs, sea stacks have particular advantages for use as ground-motion constraints: (1) they are frequently tall and thin; (2) they are widely distributed parallel to the coast, proximal to the trench and the likely megathrust rupture surface; and (3) they are formed by sea cliff retreat, meaning that their ages may be coarsely estimated as a function of distance from the coast. About 40% of the surveyed sea stacks appear to have survived multiple Cascadia megathrust earthquakes. Using a quasi-static analysis, we estimate the minimum horizontal ground accelerations that could fracture the rock pillars. We provide context for the quasi-static results by comparing them with predictions from kinematic simulations and ground-motion prediction equations. Among the sea stacks old enough to have survived multiple megathrust earthquakes (n = 16), eight yield breaking accelerations lower than the predictions, although they generally overlap within uncertainty. FGFs with the lowest breaking accelerations are distributed uniformly over 130 km of coastline. Results for inland features, such as speleothems, are in close agreement with the predictions. We conclude that FGFs show promise for investigating both past earthquake shaking and its spatial variability along the coasts of Oregon and Washington, where sea stacks are often prevalent. Future work can refine our understanding of FGF age and evolution.

In tectonically active regions, sedimentary records are overprinted by landscape response to climate, fire, and local earthquakes. We explore this issue using a new paleoclimate record developed at the Pallett Creek paleoseismic site in southern California USA, a recently incised distal fan located along the San Andreas Fault at the base of a 35 km2 catchment in the San Gabriel Mountains. To date, we have analyzed 6 m of section, spanning the last 1300 yr, for grain size, total organic material (TOM), carbon/nitrogen (C/N) ratios, magnetic susceptibility, and charcoal count. Existing C-14 dates (Scharer et al., 2011) inform rates of sediment deposition and charcoal accumulation (CHAR). Additional dating and macrofossil analysis is ongoing. Sedimentological variability within the section is dominated by two general units. Unit 1 is characterized by high % clay, % silt, and % TOM, while Unit 2 is distinctly coarser with higher % sand and lower % TOM. Pulses of high CHAR occur from 1150-1260 yr BP and during the Little Ice Age (100-500 yr BP) and are associated with high sedimentation rates (0.3-2 cm/yr), while only a few relatively weak fire episodes are recorded in the Medieval Climate Anomaly (700-1000 yr BP), despite similarly high sedimentation rates (0.6 cm/yr). Ten earthquakes documented at the site (Sieh et al., 1989) occurred about every 135 years and impart no obvious short-term impact on sedimentation rates, perhaps reflecting the distance between the site and steeper portions of the drainage network (>4 km) likely to produce mass wasting. Overall, the landscape response of this large, integrated catchment appears to reflect a stronger influence of fire and climate than earthquakes. Future work will focus on the impact of the fire episodes on sediment delivery and resultant paleoearthquake ages.

ABSTRACT Induced seismicity predominantly occurs along faults that are optimally oriented to the local principal compressive stress direction, and the characterization of these stress orientations is an important component of understanding seismic hazards. The seismicity rate in southern Kansas rapidly increased in 2013 primarily due to the disposal of large volumes of wastewater into the Arbuckle Group. Previously, local stress orientations in this area were poorly constrained, which limited our understanding of the complex faulting and diverse earthquake mechanisms in this region. We use shear-wave splitting and focal mechanism inversion techniques to create multiple, independent estimates of maximum horizontal stress directions (SHmax) and their spatial variation across the study area. We then create an integrated model of stress orientations for southern Kansas and northern Oklahoma using our local results in conjunction with previous, regional stress orientation estimates. We find that SHmax in both southern Kansas and central Oklahoma exhibits an east-northeast (∼N78° E) orientation, and these regions bound a northeast (∼N59° E) rotation within a ∼20 km area in northern Oklahoma near the Nemaha ridge.

We provide an overview of a 2019 workshop on the use of fragile geologic features (FGFs) to evaluate seismic hazard models. FGFs have been scarcely utilized in the evaluation of seismic hazard models, despite nearly 30 yr having passed since the first recognition of their potential value. Recently, several studies have begun to focus on the implementation of FGFs in seismic hazard modeling. The workshop was held to capture a “snapshot” of the state-of-the-art in FGF work and to define key research areas that would increase confidence in FGF-based evaluation of seismic hazard models. It was held at the annual meeting of the Southern California Earthquake Center on 8 September 2019, and the conveners were Mark Stirling (University of Otago, New Zealand) and Michael Oskin (University of California, Davis). The workshop attracted 44 participants from a wide range of disciplines. The main topics of discussion were FGF fragility age estimation (age at which an FGF achieved its current fragile geometry), fragility estimation, FGF-based evaluation of seismic hazard models, and ethical considerations relating to documentation and preservation of FGFs. There are now many scientists working on, or motivated to work on, FGFs, and more types of FGFs are being worked on than just the precariously balanced rock (PBR) variety. One of the ideas presented at the workshop is that fragility ages for FGFs should be treated stochastically rather than assuming that all share a common age. In a similar vein, new studies propose more comprehensive methods of fragility assessment beyond peak ground acceleration and peak ground velocity-based approaches. Two recent studies that apply PBRs to evaluate probabilistic seismic hazard models use significantly different methods of evaluation. Key research needs identified from the workshop will guide future, focused efforts that will ultimately facilitate the uptake of FGFs in seismic hazard analysis.