Search USGS⌕ Search

Geology topics

Gavin P. Hayes

Publications and source records attributed to Gavin P. Hayes.

At least 19 recordsLinked to original sources

Subduction zone earthquake catalog separation tool: Implementation in the USGS 2025 Puerto Rico and U.S. Virgin Islands National Seismic Hazard Model

The U.S. Geological Survey (USGS) periodically releases updates to National Seismic Hazard Model (NSHM) for the United States and its territories leveraging current scientific knowledge and methodologies to guide public policy, building codes, and risk assessments regarding potential ground shaking due to earthquakes that may result in infrastructure damage. In subduction zones, there is a need to separate the earthquake catalog into tectonic regimes to create specific seismicity models for which the most appropriate ground‐motion models are then applied. Here, we describe newly developed methods and software, called CatSep, that classifies subduction zone events into three primary tectonic regimes: crustal, interface, and intraslab. This method incorporates information about the location of the earthquake relative to the subducting slab, the depth of the Mohorovičić discontinuity, and the earthquake’s moment tensor. Applying this method is a first step in the NSHM workflow for regions covering U.S. subduction zones. Results using this subduction zone earthquake catalog separation tool for the 2025 Puerto Rico and U.S. Virgin Islands NSHM earthquake catalog are presented and analyzed.

Puerto Rico, U.S. Virgin Islands↗

U.S. Geological Survey Earthquake Hazards Program decadal science strategy, 2024–33

Executive Summary Earthquakes represent one of our Nation’s most significant and costly natural hazards, with estimated annual loses from earthquakes close to $15 billion in 2023. Over the past two centuries, 37 U.S. States have experienced an earthquake exceeding a magnitude of 5, and 50 percent of States have a significant potential for future damaging shaking; these statistics speak to the need for nationwide interest and investment in earthquake hazard characterization and risk reduction. Authorized under the Earthquake Hazards Reduction Authorization Act, the U.S. Geological Survey (USGS) Earthquake Hazards Program (EHP) provides the scientific information, situational awareness, and knowledge necessary to reduce deaths, injuries, and economic losses from earthquakes and earthquake-induced tsunamis, landslides, and soil liquefaction. The EHP supports activities in three focused topical areas: (1) earthquake monitoring, (2) hazard assessment, and (3) applied research, using the results of each—and the coordination among them—to further support risk translation and communication in regions at risk nationwide. For earthquake monitoring, the Advanced National Seismic System (ANSS), a cooperative effort of USGS networks, university partner regional seismic networks, and real-time geodetic networks, collects and analyzes data on earthquakes; issues timely, reliable notifications of their occurrence and impacts; and provides data for earthquake research, hazard, and risk assessment as a foundation for building an earthquake-resilient Nation. The USGS-operated ShakeAlert Earthquake Early Warning system is a recent addition to EHP’s ANSS infrastructure. In the realm of earthquake hazard assessment, the EHP contributes to earthquake risk mitigation strategies by developing the National Seismic Hazard Model and maps, and other related products, that describe the likelihood and potential effects of earthquakes nationwide, especially in the urban areas of highest risk. The EHP also conducts research on the causes, characteristics, and effects of earthquakes and prioritizes work that directly increases the accuracy and precision of earthquake hazards assessments, earthquake forecasts, and earthquake monitoring and situational-awareness products and that supports the Nation’s earthquake mitigation practices. Bridging the EHP’s efforts across research, hazard assessments, and earthquake monitoring is a broad and comprehensive collection of earthquake information products, including the National Seismic Hazard Model, ShakeAlert, and other products describing impact, such as ShakeMap and PAGER (Prompt Assessment of Global Earthquakes for Response), which have been developed and integrated into EHP’s real-time monitoring systems. EHP funds external partners to carry out many important collaborative activities through an active external grants program—one of the largest in the USGS—and through cooperative agreements with other partners such as the university-operated regional seismic networks, funded as part of the ANSS. To continue its support of earthquake hazard characterization and risk reduction, the EHP aims to strengthen its foundational products and practices while positioning itself to respond to the evolving needs of the Nation and follow best practices of the scientific community. This document describes a strategy for the program to ensure it can meet these demands. The foundational priorities outlined in this strategy represent those activities that remain critical to the core functionality of the program and those that can be supported under current fiscal year 2024-level appropriations. Priorities described as aspirational are important for future growth, and to maintain the program’s position as a leading global resource in earthquake science, but would require increases in appropriated funding to be fully realized. Across the program’s portfolio of activities, several major themes have been identified as the most critical activities to advance EHP science over the coming decade. Together, these activities provide the framework necessary to integrate critical hazard characterization and risk reduction activities across the program. They provide the structure for research to advance the understanding of where, when, and why earthquakes occur and how we can use improved knowledge to drive short-term and actionable forecasts of seismic activity. They expand the usefulness of critical earthquake products and advance the sophistication of those products to keep pace with the rapidly evolving needs of an ever-expanding user base while maintaining the position of the USGS as a global leader in earthquake science. Focus on system-level science. Establish an automated earthquake-processing pipeline. Enhance the accuracy and reliability of the ShakeAlert Earthquake Early Warning system and plan for extension to other regions. Implement time-dependent earthquake forecasting. Develop physically realistic models. Expand computational capacity. This science strategy is organized into three primary sections. The first section provides an overview of the EHP and its budget, governance, and program council. Readers familiar with the program may wish to focus on the second section, which describes the core of the science strategy, including priorities across each of the EHP’s major program activities in monitoring, hazard assessment, and targeted research. The third section outlines science priorities that cut across program activities, including those involving collaborations external to the EHP.

American Samoa, Guam, Northern Marianas Islands, U↗

A global catalog of calibrated earthquake locations

We produced a globally distributed catalog of earthquakes and nuclear explosions with calibrated hypocenters, referred to as the Global Catalog of Calibrated Earthquake Locations (GCCEL). This dataset currently contains 18,782 events in 289 clusters with >3.2 million arrival times observed at 19,258 stations. The term “calibrated” refers to the property that the hypocenters are minimally biased by unknown Earth structure. In addition, we calculate uncertainties using empirically determined variability of the arrival‐time data itself, specific to each calibrated cluster of hypocenters. Outliers in the arrival‐time dataset are removed based on measured variability of the data. In each cluster, we estimate the empirically determined uncertainty for each set of station‐phase arrival times. We use a version of the hypocentroidal decomposition multiple event relocation algorithm specifically adapted for calibrated relocations of clusters of seismic events. Most clusters are calibrated by fitting the subset of direct crustal first arrivals ( Pg and Sg ) with a locally appropriate travel‐time model to estimate the cluster hypocentroid. A few clusters are calibrated by aligning the pattern of relative locations in space and time with one or more events for which a ground‐truth hypocenter is available from an independent source with known uncertainty, such as a nuclear explosion. Epicentral uncertainties in GCCEL typically range from 1 to 5 km with a 90% confidence interval. Most events have depth constraint from one or more sources, usually with an uncertainty of ≤5 km. GCCEL is a significant resource for research at local, regional, and global scales because it provides minimally biased absolute hypocenters, meaningful associated error estimates, and curated arrival times as a reference dataset that can be used as prior constraints in the development of new regional, national, and global earthquake catalogs; validation of new location techniques; and the generation of advanced Earth models.

Seismological Research Letters↗

A ground motion model for GNSS peak ground displacement

We present an updated ground‐motion model (GMM) for M w "> M w Mw 6–9 earthquakes using Global Navigation Satellite Systems (GNSS) observations of the peak ground displacement (PGD). Earthquake GMMs inform a range of Earth science and engineering applications, including source characterization, seismic hazard evaluations, loss estimates, and seismic design standards. A typical GMM is characterized by simplified metrics describing the earthquake source (magnitude), observation distance, and site terms. Most often, GMMs are derived from broadband seismometer and accelerometer observations, yet during strong shaking, these traditional seismic instruments are affected by baseline offsets, leading to inaccurate recordings of low‐frequency ground motions such as displacement. The incorporation of geodetic data sources, particularly for characterizing the unsaturated ground displacement of large‐magnitude events, has proven valuable as a complement to traditional seismic approaches and led to the development of an initial point‐source GMM based on PGD estimated from high‐rate GNSS data. Here, we improve the existing GMM to more effectively account for fault finiteness, slip heterogeneity, and observation distance. We evaluate the limitations of the currently available GNSS earthquake data set to calibrate the GMM. In particular, the observed earthquake data set is lacking in observations within 100 km of large‐magnitude events ( ⁠ M w > 8 "> M w > 8 Mw>8 ⁠ ), inhibiting evaluation of fault dimensions for earthquakes too large to be represented as point sources in the near field. To that end, we separately consider previously validated synthetic GNSS waveforms within 10–1000 km of M w "> M w Mw 7.8–9.3 Cascadia subduction zone scenario ruptures. The synthetic data highlight the importance of fault distance rather than point‐source metrics and improve our preparedness for large‐magnitude earthquakes with spatiotemporal qualities unlike those in our existing data set.

Bulletin of the Seismological Society of America↗

Seismic monitoring during crises at the NEIC in support of the ANSS

Over the past two decades, the U.S. Geological Survey (USGS) National Earthquake Information Center (NEIC) has overcome many operational challenges. These range from minor disruptions, such as power outages, to significant operational changes, including system reconfiguration to handle unique earthquake sequences and the need to handle distributed work during a pandemic. Our ability to overcome crises is built on the development and implementation of a continuity of operations plan, well‐designed infrastructure, adaptive software systems, experienced staff, and extensive collaboration. The NEIC does not operate in a vacuum but benefits from contributions of United States and international seismic networks. Similarly, the overall resilience of earthquake monitoring in the United States and around the globe benefits from the NEIC’s role as the national center for the Advanced National Seismic System (ANSS). Here, we highlight significant adaptations the NEIC has made in the face of crises. We discuss the COVID‐19 pandemic, which represents the most significant operational crisis to impact the NEIC. The NEIC has maintained continuous operations during the ongoing COVID‐19 pandemic by shifting from a fully onsite operations center to a distributed hybrid of onsite and telework staffing. We then discuss cases in which the NEIC has supported regional monitoring in the face of significant crises. In 2018, the NEIC assisted the Hawaiian Volcano Observatory with the Kīlauea volcano eruption by responding to large events, implementing contingency monitoring procedures, and calculating moment magnitudes for the low‐frequency caldera collapses. Impacts of a crisis extend beyond the immediate response and often require a significant postevent assessment and a rebuilding phase. After the 2017 Hurricane Maria, the NEIC, the USGS National Strong‐Motion Program, and the USGS Albuquerque Seismological Laboratory worked with the Puerto Rico Seismic Network and the Puerto Rico Strong‐Motion program to assess, plan, and implement upgrades at sites that experienced storm damage.

Seismological Research Letters↗

Leveraging deep learning in global 24/7 real-time earthquake monitoring at the National Earthquake Information Center

Machine‐learning algorithms continue to show promise in their application to seismic processing. The U.S. Geological Survey National Earthquake Information Center (NEIC) is exploring the adoption of these tools to aid in simultaneous local, regional, and global real‐time earthquake monitoring. As a first step, we describe a simple framework to incorporate deep‐learning tools into NEIC operations. Automatic seismic arrival detections made from standard picking methods (e.g., short‐term average/long‐term average [STA/LTA]) are fed to trained neural network models to improve automatic seismic‐arrival (pick) timing and estimate seismic‐arrival phase type and source‐station distances. These additional data are used to improve the capabilities of the NEIC associator. We compile a dataset of 1.3 million seismic‐phase arrivals that represent a globally distributed set of source‐station paths covering a range of phase types, magnitudes, and source distances. We train three separate convolutional neural network models to predict arrival time onset, phase type, and distance. We validate the performance of the trained networks on a subset of our existing dataset and further extend validation by exploring the model performance when applied to NEIC automatic pick data feeds. We show that the information provided by these models can be useful in downstream event processing, specifically in seismic‐phase association, resulting in reduced false associations and improved location estimates.

Seismological Research Letters↗

Geometric controls on megathrust earthquakes

The role of subduction zone geometry in the nucleation and propagation of great-sized earthquake ruptures is an important topic for earthquake hazard, since knowing how big an earthquake can be on a given fault is fundamentally important. Past studies have shown subducting bathymetric features (e.g. ridges, fracture zones, seamount chains) may arrest a propagating rupture. Other studies have correlated the occurrence of great-sized earthquakes with flat megathrusts and homogenous stresses over large distances. It remains unclear, however, how subduction zone geometry and the potential for great-sized earthquakes ( M 8+) are quantifiably linked—or indeed whether they can be. Here, we examine the potential role of subduction zone geometry in limiting earthquake rupture by mapping the planarity of seismogenic zones in the Slab2 subduction zone geometry database. We build from the observation that historical great-sized earthquakes have preferentially occurred where the surrounding megathrust is broadly planar, and we use this relationship to search for geometrically similar features elsewhere in subduction zones worldwide. Assuming geometry exerts a primary control on earthquake propagation and termination, we estimate the potential size distribution of large ( M 7+) earthquakes and the maximum earthquake magnitude along global subduction faults based on geometrical features alone. Our results suggest that most subduction zones are capable of hosting great-sized earthquakes over much of their area. Many bathymetric features previously identified as barriers are indistinguishable from the surrounding megathrust from the perspective of slab curvature, meaning that they either do not play an important role in arresting earthquake rupture or that their influence on slab geometry at depth is not resolvable at the spatial scale of our subduction zone geometry models.

Geophysical Journal International↗

Seismicity of the Earth 1900–2018

This map illustrates 119 years of global seismicity in the context of global plate tectonics and the Earth’s physiography. Primarily designed for use by earth scientists, engineers, and educators, this map provides a comprehensive overview of strong (magnitude [ M ] 5.5 and larger) earthquakes since 1900. The map clearly identifies the locations of the “great” earthquakes ( M 8.0 and larger) and the aftershock or rupture area (green fill), if known, of the M 8.3 or larger earthquakes. The circular earthquake symbols are scaled to be proportional to the moment magnitude and therefore to the area of faulting, thus providing a better understanding of the relative sizes and distribution of earthquakes in the magnitude range 5.5 to 9.5. Plotting the known rupture or aftershock areas (which are closely related) of the largest earthquakes also provides a better appreciation of the faulting extent of some of the most famous and damaging instrumentally recorded earthquakes in modern history.

Scientific Investigations Map↗

Structural control on megathrust rupture and slip behavior: Insights from the 2016 Mw 7.8 Pedernales Ecuador earthquake

The heterogeneous seafloor topography of the Nazca Plate as it enters the Ecuador subduction zone provides an opportunity to document the influence of seafloor roughness on slip behavior and megathrust rupture. The 2016 Mw 7.8 Pedernales Ecuador earthquake was followed by a rich and active postseismic sequence. An internationally coordinated rapid response effort installed a temporary seismic network to densify coastal stations of the permanent Ecuadorian national seismic network. A combination of 82 onshore short and intermediate period and broadband seismic stations and six ocean bottom seismometers recorded the postseismic Pedernales sequence for over a year after the mainshock. A robust earthquake catalog combined with calibrated relocations for a subset of magnitude ≥4 earthquakes shows pronounced spatial and temporal clustering. A range of slip behavior accommodates postseismic deformation including earthquakes, slow slip events, and earthquake swarms. Models of plate coupling and the consistency of earthquake clustering and slip behavior through multiple seismic cycles reveal a segmented subduction zone primarily controlled by subducted seafloor topography, accreted terranes, and inherited structure. The 2016 Pedernales mainshock triggered moderate to strong earthquakes (5 ≤ M ≤ 7) and earthquake swarms north of the mainshock rupture close to the epicenter of the 1906 Mw 8.8 earthquake and in the segment of the subduction zone that ruptured in 1958 in a Mw 7.7 earthquake.

Journal of Geophysical Research: Earth Surface↗

An update of USGS bear-real-time earthquake shaking and impact products

We report on advancements in both hazard and consequence modeling that form the core of the U.S. Geological Survey’s (USGS) strategy to improve rapid earthquake shaking and loss estimates. Whereas our primary goal is to improve our operational capabilities of the USGS National Earthquake Information Center, the science, software, and datasets behind these systems continue to advance uses and studies of earthquake shaking and impact by the seismological, engineering, financial, and risk modeling communities. Several important updates to our integrated shaking and impact products are outlined and we introduce new earthquake information products that have recently been brought online, including rapid ground failure estimates and more spatially refined loss estimates domestically (in the U.S). We continue to compile, develop, and refine key openly available models and datasets that contribute to calibrating these systems and report on the collection and storage of new inventories. We also describe some of the basic operational considerations in the current generation of these shaking and loss-estimation systems. A key aspect of the product integration and development is leveraging earthquake-hazard and loss-modeling science done internally (within the USGS) and by external researchers and collaborators. Lastly, we outline new opportunities for further research and development by emphasizing scientific, data, and application gaps and challenges that must be solved in order to improve our shaking and impact information tools.

Conference Paper↗

USGS near-real-time products-and their use-for the 2018 Anchorage earthquake

In the minutes to hours after a major earthquake, such as the recent 2018 M w "> M w 7.1 Anchorage event, the U.S. Geological Survey (USGS) produces a suite of interconnected earthquake products that provides diverse information ranging from basic earthquake source parameters to loss estimates. The 2018 Anchorage earthquake is the first major domestic earthquake to occur since several new USGS products have been developed, thus providing an opportunity to discuss the newly expanded USGS earthquake product suite, its timeliness, performance, and reception. Overall, the products were relatively timely, accurate, well received, and widely used, including by the media, who used information and visualizations from many products to frame their early reporting. One downside of the codependence of multiple products is that reasonable updates to upstream products (e.g., magnitude and source characterization) can result in significant changes to downstream products; this was the case for the Anchorage earthquake. However, the coverage of strong‐motion stations and felt reports was so dense that the ShakeMap and downstream products were relatively insensitive to changes in magnitude or fault‐plane orientation once the ground‐motion data were available. Shaking and loss indicators initially fluctuated in the first hour or two after the earthquake, but they stabilized quickly. To understand how the products are being used and how effectively they are being communicated, we analyze the media coverage of USGS earthquake products. Most references to USGS products occurred within the first 48 hr after the event. The lack of coverage after 48 hr could indicate that longer‐term products addressing what actions the USGS is taking or what early reconnaissance has revealed might be useful for those people wanting additional information about the earthquake.

Alaska↗

National earthquake information center strategic plan, 2019–23

Executive Summary Damaging earthquakes occur regularly around the world; since the turn of the 20th century, hundreds of earthquakes have caused significant loss of life and (or) millions of dollars or more in economic losses. While most of these did not directly affect the United States and its Territories, by studying worldwide seismicity we can better understand how to mitigate the effects of earthquakes when they do occur within U.S. borders. Within the U.S. Government, this mandate falls on the U.S. Geological Survey (USGS) National Earthquake Information Center (NEIC), which has the statutory responsibility for monitoring and reporting on earthquakes domestically and globally. The NEIC has been operating since 1966, and throughout its history has been recognized as a world leader for earthquake information. For much of this time, NEIC has been cooperating with a number of regional seismic networks (RSNs) which operate in areas of heightened seismicity in the United States. In 2000, the Advanced National Seismic System (ANSS) was founded as a cooperative umbrella for earthquake-related data collection, analysis, and dissemination in the United States, thereby promoting advanced interoperability between the NEIC and RSN partners. The NEIC also cooperates and coordinates with dozens of global seismic networks. At present (2019), NEIC acquires real-time waveform data from more than 2,000 seismic stations worldwide, contributed from more than 130 seismic networks. Since 2006, the NEIC has operated on a 24-hour, 7-days per week (24/7) basis, and reports on about 30,000 earthquakes per year. Soon after the occurrence of a significant global earthquake, notifications are issued to government representatives, aid agencies, the press, and members of the general public by the Earthquake Notification Service (ENS), electronic feeds, and through the USGS Earthquake Hazards Program (EHP) website. Event-specific web pages provide detailed source parameter information outlining the location and magnitude of the earthquake, including more detailed source characteristics like moment magnitude and focal mechanisms and finite fault solutions. Further, NEIC produces a suite of real-time situational awareness products, including ShakeMap, ShakeCast, Did-You-Feel-It? (DYFI?), and Prompt Assessment of Global Earthquakes for Response (PAGER), to characterize the shaking resulting from the earthquake and the impact it is likely to have on nearby populations and infrastructure. All of these products are ultimately archived in the ANSS Comprehensive Catalog (ComCat), hosted and served by the NEIC. The NEIC also pursues an active research program to improve its ability to characterize earthquakes and understand their hazards. These efforts are all aimed at mitigating the risks of earthquakes to humankind. To maintain its prominent position in earthquake monitoring, the NEIC must continue to evolve, concurrently improving its operations and 24/7 robustness, streamlining services and infrastructure, and keeping pace with research and innovation in the field of seismology. This document outlines how the NEIC might best achieve such goals, by describing specific avenues and opportunities for development in the next five years (2019–23). Several key areas of operational and research focus are identified in this plan as being of the highest importance. First, NEIC must finalize improvements to its regional monitoring capabilities, including the implementation of a variety of improved earthquake detection and association algorithms. One of the most exciting avenues of recent research expansion in earthquake monitoring has involved the use of machine learning; NEIC must explore the benefits of machine learning for improved earthquake detection and source characterization. NEIC also needs to address issues related to the timeliness of earthquake information, exploring the benefits of distributing information as it becomes available, rather than when certain quality criteria are met. To that end, the incorporation of real-time Global Positioning System (GPS) data into the NEIC operational workflow will help improve the speed and accuracy of information for moderate-to-large earthquakes. Finally, NEIC should explore how to further expand and improve the quality and content of the products served during earthquake response efforts, including the generation of new earthquake sequence-specific products, adding an evolutionary component to earthquake information, and continued improvements to earthquake impact products.

Circular↗

Global earthquake response with imaging geodesy: recent examples from the USGS NEIC

The U.S. Geological Survey National Earthquake Information Center leads real-time efforts to provide rapid and accurate assessments of the impacts of global earthquakes, including estimates of ground shaking, ground failure, and the resulting human impacts. These efforts primarily rely on analysis of the seismic wavefield to characterize the source of the earthquake, which in turn informs a suite of disaster response products such as ShakeMap and PAGER. In recent years, the proliferation of rapidly acquired and openly available in-situ and remotely sensed geodetic observations has opened new avenues for responding to earthquakes around the world in the days following significant events. Geodetic observations, particularly from interferometric synthetic aperture radar (InSAR) and satellite optical imagery, provide a means to robustly constrain the dimensions and spatial complexity of earthquakes beyond what is typically possible with seismic observations alone. Here, we document recent cases where geodetic observations contributed important information to earthquake response efforts – from informing and validating seismically-derived source models to independently constraining earthquake impact products – and the conditions under which geodetic observations improve earthquake response products. We use examples from the 2013 Mw7.7 Baluchistan, Pakistan, 2014 Mw6.0 Napa, California, 2015 Mw7.8 Gorkha, Nepal, and 2018 Mw7.5 Palu, Indonesia earthquakes to highlight the varying ways geodetic observations have contributed to earthquake response efforts at the NEIC. We additionally provide a synopsis of the workflows implemented for geodetic earthquake response. As remote sensing geodetic observations become increasingly available and the frequency of satellite acquisitions continues to increase, operational earthquake geodetic imaging stands to make critical contributions to natural disaster response efforts around the world.

Remote Sensing↗

The 12 November 2017 Mw 7.3 Ezgeleh–Sarpolzahab (Iran) earthquake and active tectonics of the Lurestan arc

The 12 November 2017 M w 7.3 Ezgeleh‐Sarpolzahab earthquake is the largest instrumentally recorded earthquake in the Zagros Simply Folded Belt by a factor of ∼10 in seismic moment. Exploiting local, regional, and teleseismic data and synthetic aperture radar interferometry imagery, we characterize the rupture, its aftershock sequence, background seismicity, and regional tectonics. The mainshock ruptured slowly (∼2 km/s), unilaterally southward, for ∼40 km along an oblique (dextral‐thrust) fault that dips ∼14°E beneath the northwestern Lurestan arc. Slip is confined to basement depths of ∼12–18 km, resolvably beneath the sedimentary cover which is ∼8 km thick in this area. The gentle dip angle and basement location allow for a broad slip area, explaining the large magnitude relative to earthquakes in the main Fars arc of the Zagros, where shallower, steeper faults are limited in rupture extent by weak sedimentary layers. Early aftershocks concentrate around the southern and western edges of the mainshock slip area and therefore cluster in the direction of rupture propagation, implying a contribution from dynamic triggering. A cluster of events ∼100 km to the south near Mandali (Iraq) reactivated the ∼50° dipping Zagros Foredeep Fault. The basement fault responsible for the Ezgeleh‐Sarpolzahab earthquake probably accounts for the ∼1 km elevation contrast between the Lurestan arc and the Kirkuk embayment but is distinct from sections of the Mountain Front Fault that define frontal escarpments elsewhere in the Zagros. It may be related to a seismic interface underlying the central and southern Lurestan arc, and a key concern is whether or not the more extensive regional structure is also seismogenic.

Journal of Geophysical Research: Solid Earth↗

Characterizing large earthquakes before rupture is complete

Whether large and very large earthquakes are distinguishable from each other early on in the rupture process has been a subject often debated over the past several decades. Studies have shown that the frequency content of radiated seismic energy in the first few seconds of an earthquake scales with the final magnitude of the event, implying determinism. Other studies have shown that the recordings of ground displacement from small-to-moderate sized earthquakes are indistinguishable, and thus earthquakes share a universal early rupture process. Regardless of how earthquakes start, however, at some point in the rupture process events of different sizes must be distinguishable from one another. If that difference occurs early - i.e., before the rupture duration of the smaller event - this implies that earthquakes demonstrate some level of determinism. Here we show through analysis of a large database of source time functions and near-source displacement records that after an initiation phase, ruptures of M7-9 earthquakes organize into a slip pulse, the kinematic properties of which scale with magnitude. As such, early in the rupture process - after about 10s - large and very large earthquakes demonstrate different properties and can thus be distinguished.

Science Advances↗

Incorporating teleseismic tomography data into models of upper mantle slab geometry

Earthquake-based models of slab geometry are limited by the distribution of earthquakes within a subducting slab, which is often heterogeneous. The fast seismic velocity signature of slabs in tomography studies is independent of the distribution of earthquakes within the slab, providing a critical constraint on slab geometry when earthquakes are absent. In order to utilize this constraint, researchers typically hand-contour images of subducting slabs in tomography models, leading to a subjective final slab model. With this paper, we present an automated procedure for extracting slab geometry from teleseismic tomography volumes that limits this subjectivity and provides constraints on the structure of aseismic segments of slabs. This procedure is designed as a complement to earthquake-based slab models rather than as a replacement, which can help to broaden the extent of existing subduction zone geometry databases.

Geophysical Journal International↗

Slab2, a comprehensive subduction zone geometry model

Subduction zones are home to the most seismically active faults on the planet. The shallow megathrust interfaces of subduction zones host our largest earthquakes and are likely the only faults capable of magnitude 9+ ruptures. Despite these facts, our knowledge of subduction zone geometry—which likely plays a key role in determining the spatial extent and ultimately the size of subduction zone earthquakes—is incomplete. We calculated the three-dimensional geometries of all seismically active global subduction zones. The resulting model, called Slab2, provides a uniform geometrical analysis of all currently subducting slabs.

Science↗

Systematic observations of the slip pulse properties of large earthquake ruptures

In earthquake dynamics there are two end member models of rupture: propagating cracks and self-healing pulses. These arise due to different properties of faults and have implications for seismic hazard; rupture mode controls near-field strong ground motions. Past studies favor the pulse-like mode of rupture; however, due to a variety of limitations, it has proven difficult to systematically establish their kinematic properties. Here we synthesize observations from a database of >150 rupture models of earthquakes spanning M 7– M 9 processed in a uniform manner and show the magnitude scaling properties of these slip pulses indicates self-similarity. Further, we find that large and very large events are statistically distinguishable relatively early (at ~15 s) in the rupture process. This suggests that with dense regional geophysical networks strong ground motions from a large rupture can be identified before their onset across the source region.

Geophysical Research Letters↗