Geology Reports⌕ Search

Geology topics

Brad T. Aagaard

Publications and source records attributed to Brad T. Aagaard.

At least 37 records · Page 2Linked to original sources

2019 M7.1 Ridgecrest earthquake slip distribution controlled by fault geometry inherited from Independence dike swarm

Faults often form through reactivation of pre-existing structures, developing geometries and mechanical properties specific to the system’s geologic inheritance. Competition between fault geometry and other factors (e.g., lithology) to control slip at Earth’s surface is an open question that is central to our knowledge of fault processes and seismic hazards. Here we use remote sensing data and field observations to investigate the origin of the 2019 M 7.1 Ridgecrest, California, earthquake rupture geometry and test its impact on the slip distribution observed at Earth’s surface. Common geometries suggest the fault system evolved through reactivation of structures within the surrounding Independence dike swarm (IDS). Mechanical models testing a range of fault geometries and stress fields indicate that the inherited rupture geometry strongly controlled the M 7.1 earthquake slip distribution. These results motivate revisiting the development of other large-magnitude earthquake ruptures (1992 M 7.3 Landers, 1999 M 7.1 Hector Mine) and tectonic provinces within the IDS.

California↗

2022 Crustal Deformation Modeling Workshop Report

The 2022 Crustal Deformation Modeling Workshop was held June 20–24 at the Colorado School of Mines in Golden, Colorado. The workshop included two days of tutorials on the use of the open-source software PyLith for crustal deformation modeling followed by three days of science talks and discussions. The workshop focused on three primary themes: (1) Earthquake cycle modeling; (2) Inversions for fault slip; and (3) Faulting, fluids, and surface loading. The talks highlighted how computational models provide insight into intriguing observations of Earth and planetary behavior. These included (1) earthquake synchronization of rupture patches due to their close proximity to each other, (2) the influence of fault geometry and damage zones on hypocenter depth and rupture propagation, (3) a lack of steady-state faulting behavior due to long time scales for grain size evolution in the mid-crust, (4) crustal deformation due to tidal, hydrological, and atmospheric loads, and (5) plumes of gas and icy particles due to tidal driven faulting on Enceladus (one of Saturn’s moons). The talks also described new computational modeling capabilities for incorporating complex geologic structure into Bayesian inversions for fault slip and efficient implementation of earthquake cycle models using a symmetric interior discontinuous Galerkin method. The complete agenda is available on the Computational Infrastructure for Geodynamic (CIG) website.

Report↗

2023 PyLith Hackathon report

The 3rd Pylith Hackathon was held June 12–17, 2023, at the Colorado School of Mines in Golden, Colorado with funding from the Computational Infrastructure for Geodynamics (CIG). The hackathon involved 17 participants working on 5 different projects to implement new features and create new examples for the PyLith crustal deformation modeling software. The projects included (1) spontaneous rupture using fault friction, (2) extending the poroelasticity implementation, (3) developing 2D and 3D examples involving strike-slip faults, (4) integrating PyLith with the cascading adaptive transitional metropolis in parallel (CATMIP) Bayesian inversion framework for use in studies inverting for static fault slip, and (5) adding self-gravitation using the current multiphysics formulation in PyLith. Participants learned how to navigate the PyLith code base, implement point-wise functions for governing equations and bulk and fault rheologies using the finite-element method, extend the code using the modular, object-oriented design, write examples that demonstrate how to use the new features in PyLith simulations, and implement method of manufactured solutions tests and full-scale tests. The PyLith development team benefitted from discussions with the other participants (contributors) about the technical aspects of the various projects as well as general discussions about PyLith design. The in-person format and 6-day duration allowed the groups to make significant progress. Participants appreciated the project-based organization of the hackathon and recommended that future hackathons include online meetings of the various projects before the in-person gathering to self-organize and prepare. Sarah Minson (remote) provided technical advice on the use of the CATMIP Bayesian inversion framework, and this type of participation could be expanded to allow additional technical presentations and advice on various topics in future hackathons.

Report↗

Implementation of basin models and sediment depth terms in the 2023 update of the U.S. National Seismic Hazard Model: Example from Reno, Nevada

We present a framework to evaluate the inclusion of candidate basin depth models in the U.S. Geological Survey National Seismic Hazard Model. We compute intensity measures (peak and spectral amplitudes) from uniformly processed earthquake ground motions in and around the basin of interest and compare these to ground-motion model (GMM) estimates over a range of oscillator periods. The GMMs use depth to specific shear-wave velocity isosurfaces (Z x ) as a proxy for basin depth. We quantify whether the GMM estimates using Z x from the candidate basin depth model outperform the default estimates based on V S30 (the time-averaged shear wave velocity in the upper 30 m). We partition GMM residuals into event and site terms and compare site terms for stations within the basin to non-basin sites. We apply this framework to the greater Reno, Nevada, region, which has shallow basin depths (less than 450 m) and empirical amplifications of up to a factor of 2.2 at 1.0 s. There are no strong trends in site terms with basin depth, and the use of Z x values from the new candidate basin model shows marginal improvement of GMM total residuals compared to the default V S30 -based model (the condition of no relative basin amplification).

Nevada↗

Improving the Development Pipelines for USGS Earthquake Hazards Program Real-Time and Scenario Products

The real-time and scenario products of the U.S. Geological Survey (USGS) Earthquake Hazards Program, such as the ComCat catalog, Did You Feel It?, ShakeMap, ShakeCast, and PAGER, are highly visible and used by a wide variety of stakeholders. We propose two significant enhancements to the development pipelines for the Earthquake Hazards Program real-time and scenario products that have far-reaching benefits. First, we propose incorporating processed and archived ground-motion records into the data streams for real-time products. This increases reproducibility and transparency for ShakeMap and downstream products that serve critical functions in earthquake response and long-term research. It will also provide comprehensive, open access databases of ground-motion metrics (for example, peak ground acceleration, peak ground velocity, and acceleration response spectra) and ground-motion time histories that are fundamental tools in most engineering seismology studies. Second, we propose extending the pipeline for scenario products to provide a full set of complementary products to the real-time pipeline. This would define a comprehensive set of standards for archiving scenarios, including three-dimensional ground-motion simulations, and allow the suite of scenario products to be disseminated in the same way as real-time products. Ultimately, these enhancements would increase the value of some of the most important Earthquake Hazards Program products and transform the way USGS scientists and the engineering seismology community conduct ground-motion research.

Conference Paper↗

COSMOS Ground-Motion Simulation Working Group workshops #1 and #2

These 2 workshops were held in response to interest generated from sessions on the use of simulated earthquake ground motions at the 2020 and 2021 Consortium of Organizations for Strong Motion Observation Systems (COSMOS) Technical Sessions. The discussions at the Technical Sessions highlighted desires to promote the use of simulated earthquake ground motions for engineering applications and the need to coordinate efforts to validate and disseminate the data. Our first workshop focused on curating and disseminating simulated earthquake ground-motion data. Our second workshop focused on validating simulated ground-motions for engineering applications. Both workshops included a few invited presentations from providers (ground-motion simulators) and engineering user perspectives followed by open discussions. About 100 people participated in each of the workshops.

Report↗

2021 Computational Infrastructure for Geodynamics Developers Workshop

The CIG Developers Workshop resulted in a number of recommendations that we think will help expand the CIG developer community, make software more accessible to new users, and increase developer productivity through use of common infrastructure and best practices for software development. This includes building a broad user base with sufficient support through documentation, tutorials, user forums, hackathons, scientific workshops, and mentoring to maintain a healthy suite of software developers and maintainers. Communities also need to offer opportunities, like this workshop, for developer teams to interact with each other to exchange ideas, identify common infrastructure, and interact with users to discuss modeling workflows and development priorities.

Conference Paper↗

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↗

Science plan for improving three-dimensional seismic velocity models in the San Francisco Bay region, 2019–24

This five-year science plan outlines short-term and long-term goals for improving three-dimensional seismic velocity models in the greater San Francisco Bay region as well as how to foster a community effort in reaching those goals. The short-term goals focus on improving the current U.S. Geological Survey San Francisco Bay region geologic and seismic velocity model using existing data. The long-term goals focus on acquiring new data and leveraging better analytic tools to improve the model and characterize the uncertainty. The plan describes opportunities for contributions by members of the community to develop these seismic velocity models, provides current and potential users with general information on where efforts will likely be focused to improve these models and how new versions of the models will be released, and outlines funding needs and obstacles for improving and maintaining such models. Several aspects of this plan, including how to foster a community effort, are independent of the geographic region and apply to other similar efforts.

California↗

The 2019 Ridgecrest, California, earthquake sequence ground motions: Processed records and derived intensity metrics

Following the 2019 Ridgecrest, California, earthquake sequence, we compiled ground‐motion records from multiple data centers and processed these records using newly developed ground‐motion processing software that performs quality assurance checks, performs standard time series processing steps, and computes a wide range of ground‐motion metrics. In addition, we compute station and waveform metrics such as the time‐averaged shear‐wave velocity to 30 m depth ( ⁠ V S 30 "> V S 30 ⁠ ), finite‐rupture distances, and spectral accelerations. This data set includes 22,708 records from 133 events from 4 July 2019 (UTC) to 18 October 2019 with a magnitude range from 3.6 to 7.1. We expect that the rapid collection and dissemination of this information will facilitate detailed studies of these ground motions. In this article, we describe the data selection, processing steps, and how to access the data.

California↗

Research to improve ShakeAlert earthquake early warning products and their utility

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

Open-File Report↗

A suite of exercises for verifying dynamic earthquake rupture codes

We describe a set of benchmark exercises that are designed to test if computer codes that simulate dynamic earthquake rupture are working as intended. These types of computer codes are often used to understand how earthquakes operate, and they produce simulation results that include earthquake size, amounts of fault slip, and the patterns of ground shaking and crustal deformation. The benchmark exercises examine a range of features that scientists incorporate in their dynamic earthquake rupture simulations. These include implementations of simple or complex fault geometry, off‐fault rock response to an earthquake, stress conditions, and a variety of formulations for fault friction. Many of the benchmarks were designed to investigate scientific problems at the forefronts of earthquake physics and strong ground motions research. The exercises are freely available on our website for use by the scientific community.

Seismological Research Letters↗

Integrate urban‐scale seismic hazard analyses with the U.S. National Seismic Hazard Model

For more than 20 yrs, damage patterns and instrumental recordings have highlighted the influence of the local 3D geologic structure on earthquake ground motions (e.g., M "> M M 6.7 Northridge, California, Gao et al. , 1996 ; M "> M M 6.9 Kobe, Japan, Kawase, 1996 ; M "> M M 6.8 Nisqually, Washington, Frankel, Carver, and Williams, 2002 ). Although this and other local‐scale features are critical to improving seismic hazard forecasts, historically they have not been explicitly incorporated into the U.S. National Seismic Hazard Model (NSHM, national model and maps), primarily because the necessary basin maps and methodologies were not available at the national scale. Instead, the U.S. Geological Survey (USGS), its partners, and external groups developed urban seismic hazard maps (urban models and maps) that consider detailed site effects in local areas (e.g., Wong et al. , 2002 ; Cramer et al. , 2006 ; Frankel et al. , 2007 ; Graves et al. , 2011 ). The disconnect between the urban and national hazard models, however, means that the national models, which underlie U.S. building codes and other applications, do not make use of all of the scientific results informing earthquake ground‐shaking hazards. We recommend that future U.S. national seismic hazard assessment make use of all available regional information, including that in urban models. In this column, we describe the roles of and differences between the urban and national models, and discuss the obstacles to and benefits of integrating the urban models with the national model. Future progress on this issue will require further coordination and implementation efforts between the USGS and external partners.

Seismological Research Letters↗

U.S. Geological Survey National Strong-Motion Project strategic plan, 2017–22

The mission of the National Strong-Motion Project is to provide measurements of how the ground and built environment behave during earthquake shaking to the earthquake engineering community, the scientific community, emergency managers, public agencies, industry, media, and other users for the following purposes: Improving engineering evaluations and design methods for facilities and systems; Providing timely information for earthquake early warning, damage assessment, and emergency response action; and Contributing to a greater understanding of the mechanics of earthquake rupture, groundmotion characteristics, and earthquake effects.

Open-File Report↗

Subsurface geometry of the San Andreas fault in southern California: Results from the Salton Seismic Imaging Project (SSIP) and strong ground motion expectations

The San Andreas fault (SAF) is one of the most studied strike‐slip faults in the world; yet its subsurface geometry is still uncertain in most locations. The Salton Seismic Imaging Project (SSIP) was undertaken to image the structure surrounding the SAF and also its subsurface geometry. We present SSIP studies at two locations in the Coachella Valley of the northern Salton trough. On our line 4, a fault‐crossing profile just north of the Salton Sea, sedimentary basin depth reaches 4 km southwest of the SAF. On our line 6, a fault‐crossing profile at the north end of the Coachella Valley, sedimentary basin depth is &#x223C; 2 &#x2013; 3 &#x2009;&#x2009; km "> ∼ 2 – 3 km and centered on the central, most active trace of the SAF. Subsurface geometry of the SAF and nearby faults along these two lines is determined using a new method of seismic‐reflection imaging, combined with potential‐field studies and earthquakes. Below a 6–9 km depth range, the SAF dips &#x223C; 50 &#xB0; &#x2013; 60 &#xB0; "> ∼ 50 ° – 60 ° NE, and above this depth range it dips more steeply. Nearby faults are also imaged in the upper 10 km, many of which dip steeply and project to mapped surface fault traces. These secondary faults may join the SAF at depths below about 10 km to form a flower‐like structure. In Appendix D , we show that rupture on a northeast‐dipping SAF, using a single plane that approximates the two dips seen in our study, produces shaking that differs from shaking calculated for the Great California ShakeOut, for which the southern SAF was modeled as vertical in most places: shorter‐period ( T &lt; 1 &#x2009;&#x2009; s "> T < 1 s ) shaking is increased locally by up to a factor of 2 on the hanging wall and is decreased locally by up to a factor of 2 on the footwall, compared to shaking calculated for a vertical fault.

California↗

Earthquake outlook for the San Francisco Bay region 2014–2043

Using information from recent earthquakes, improved mapping of active faults, and a new model for estimating earthquake probabilities, the 2014 Working Group on California Earthquake Probabilities updated the 30-year earthquake forecast for California. They concluded that there is a 72 percent probability (or likelihood) of at least one earthquake of magnitude 6.7 or greater striking somewhere in the San Francisco Bay region before 2043. Earthquakes this large are capable of causing widespread damage; therefore, communities in the region should take simple steps to help reduce injuries, damage, and disruption, as well as accelerate recovery from these earthquakes.

California↗

The Mw6.0 24 August 2014 South Napa earthquake

The M w 6.0 South Napa earthquake, which occurred at 10:20 UTC 24 August 2014 was the largest earthquake to strike the greater San Francisco Bay area since the M w 6.9 1989 Loma Prieta earthquake. The rupture from this right‐lateral earthquake propagated mostly unilaterally to the north and up‐dip, directing the strongest shaking toward the city of Napa, where peak ground accelerations (PGAs) between 45% g and 61% g were recorded and modified Mercalli intensities (MMIs) of VII&ndash;VIII were reported. Tectonic surface rupture with dextral slip of up to 46 cm was observed on a 12.5 km long segment, some of which was along a previously mapped strand of the West Napa fault system, although the rupture extended to the north of the mapped Quaternary strand. Modeling of seismic and geodetic data suggests an average coseismic slip of 50 cm, with a maximum slip of about 1 m at depths of 10&ndash;11 km. We observed up to 35 cm of afterslip along the surface trace in the week following the mainshock, primarily along the southern half of the surface rupture that experienced relatively little coseismic offset. Relocation of the sparse aftershock sequence suggests en echelon southwest‐ and northeast‐dipping fault planes, reflective of the complex fault geometry in this region. The Napa basin and historic and late Holocene alluvial flood deposits in downtown Napa amplified the ground motions there. Few ground failures were mapped, reflecting the dry season (as well as a persistent drought that had lowered the groundwater table) and the short duration of strong shaking in the epicentral area.

California↗

Key recovery factors for the August 24, 2014, South Napa Earthquake

Through discussions between the Federal Emergency Management Agency (FEMA) and the U.S. Geological Survey (USGS) following the South Napa earthquake, it was determined that several key decision points would be faced by FEMA for which additional information should be sought and provided by USGS and its partners. This report addresses the four tasks that were agreed to. These tasks are (1) assessment of ongoing fault movement (called afterslip) especially in the Browns Valley residential neighborhood, (2) assessment of the shaking pattern in the downtown area of the City of Napa, (3) improvement of information on the fault hazards posed by the West Napa Fault System (record of past earthquakes and slip rate, for example), and (4) imagery acquisition and data processing to provide overall geospatial information support to FEMA.

California↗