Geology topics

Harley M. Benz

Publications and source records attributed to Harley M. Benz.

At least 19 recordsLinked to original sources

USGS2025

Estimating earthquake source depth using teleseismic broadband waveform modeling at the USGS National Earthquake Information Center

The U.S. Geologic Survey National Earthquake Information Center (NEIC) monitors global seismicity, producing a catalog of earthquake source parameters in near-real-time to provide information that can help mitigate the societal impact of earthquakes. The NEIC commonly relies on teleseismic observations to constrain earthquake source parameters (e.g., location, depth, magnitude, and mechanism) due to a lack of local and regional observations. For these ‘teleseismic’ events, depth phase (i.e., pP , sP ) arrival time observations provide the best estimate on source depth. However, depth phases are often difficult to accurately identify and/or pick. Therefore, NEIC relies on waveform modeling, such as those determined from W-phase ( M ww ), body wave ( M wb ), and regional ( M wr ) moment tensor estimations, to provide constraints on source depth. While depth estimates from these approaches are informative, higher frequency observations provide more precise estimates because depth phases are more prominently observed at higher frequencies. Here, we present NEIC’s relatively high-frequency (~0.04 to 1 Hz) teleseismic waveform modeling approach, termed Synthetic Depth Phase Modeling (SynDepth), for determining source depth. SynDepth was developed to provide NEIC with a tool that enables rapid, accurate, and quantifiable estimates of earthquake source depth in cases where locator depths are not reliable. This relatively simple and fast procedure searches over 1 km-incremented source depths and an expanding triangular source-time function to find the best-fitting solution. We compare automatic SynDepth solutions for a dataset of 1,216 earthquakes (M5.5-M7.6) between 2017 and 2021 to NEIC-derived depth estimates from other methods. Our approach provides a robust depth estimate for earthquakes lacking local arrival time data, and it minimizes the need for analyst review of depth-phase picks ( pP , sP ) or using predefined ‘fixed’ depths.

Seismological Research Letters↗

USGS2024

Uncertainty and spatial correlation in station measurements for mb magnitude estimation

The body‐wave magnitude (⁠⁠) is a long‐standing network‐averaged, amplitude‐based magnitude used to estimate the magnitude of seismic sources from teleseismic observations. The U.S. Geological Survey National Earthquake Information Center (NEIC) relies on in its global real‐time earthquake monitoring mission. Although waveform modeling‐based moment magnitudes are the modern standard to characterize earthquake size, is important because (1) in many cases, waveform modeling is not possible (e.g., low signal‐to‐noise events), (2) is applicable over a broad range of magnitudes, ∼M 4–7, and (3) there is a many decades‐long history of estimating magnitudes. We use the NEIC Preliminary Determination of Epicenters earthquake catalog to investigate the uncertainty in NEIC station measurements. We show that measurements are spatially correlated, which can bias event ⁠, and we describe an empirical relation between this spatial correlation and station‐to‐station distance. We further describe an approach to mitigate bias from the spatial correlation. Accounting for the spatial covariance of observations can change the event from −0.15 to 0.07 units (10th to 90th percentile) for smaller events (⁠⁠). These smaller events have the largest standard deviations ranging from 0.05 to 0.15 units (10th to 90th percentile).

The Seismic Record↗

USGS2024

Formation and evolution of the Pacific-North American (San Andreas) plate boundary: Constraints from the crustal architecture of northern California

The northward migration of the Mendocino triple junction (MTJ) drives a fundamental plate boundary transformation from convergence to translation; producing a series of strike-slip faults, that become the San Andreas plate boundary. We find that the 3-D structure of the Pacific plate lithosphere in the vicinity of the MTJ controls the location of San Andreas plate boundary formation. At the time of initiation of the Pacific-North America plate boundary (∼30 Ma), the sequential interaction with the western margin of North America of the Pioneer Fracture Zone, soon followed by the Mendocino Fracture Zone, led to the capture of a small segment of partially subducted Farallon lithosphere by the Pacific plate, termed the Pioneer Fragment (PF). Since that time, the PF has translated with the Pacific Plate along the western margin of North America. Recently developed, high-resolution seismic-tomographic imagery of northern California indicates that (a) the PF is extant, occupying the western half of the slab window, immediately south of the MTJ; (b) the eastern edge of the PF lies beneath the newly forming Maacama fault system, which develops to become the locus for the primary plate boundary structure after approximately 6–10 Ma; and (c) the location of the translating PF adjacent to the asthenosphere of the slab window generates a shear zone within and below the crust that develops into the plate boundary faults. As a result, the San Andreas plate boundary forms interior to the western margin of North America, rather than at its western edge.

California↗

USGS2023

MLAAPDE: A machine learning dataset for determining global earthquake source parameters

The Machine Learning Asset Aggregation of the Preliminary Determination of Epicenters (MLAAPDE) dataset is a labeled waveform archive designed to enable rapid development of machine learning (ML) models used in seismic monitoring operations. MLAAPDE consists of more than 5.1 million recordings of 120 s long three‐component broadband waveform data (raw counts) for P , Pn , Pg , S , Sn , and Sg arrivals. The labeled catalog is collected from the U.S. Geological Survey National Earthquake Information Center’s (NEIC) Preliminary Determination of Epicenters bulletin, which includes local to teleseismic observations for earthquakes ∼ M 2.5 and larger. Each arrival in the labeled dataset has been manually reviewed by NEIC staff. An accompanying Python module enables users to develop customized training datasets, which includes different time‐series lengths, distance ranges, sampling rates, and/or phase lists. MLAAPDE is distinct from other publicly available datasets in containing local (14%), regional (36%), and teleseismic (50%) observations, in which local, regional, and teleseismic distance are 0°–3°, 3°–30°, and 30°+, respectively. A recent version of the dataset is publicly available (see Data and Resources ), and user‐specific versions can be generated locally with the accompanying software. MLAAPDE is an NEIC supported, curated, and periodically updated dataset that can contribute to seismological ML research and development.

Seismological Research Letters↗

USGS2023

Rapid characterization of the February 2023 Kahramanmaraş, Turkey, earthquake sequence

The 6 February 2023 M w "> M w 7.8 Pazarcık and subsequent M w "> M w 7.5 Elbistan earthquakes generated strong ground shaking that resulted in catastrophic human and economic loss across south‐central Türkiye and northwest Syria. The rapid characterization of the earthquakes, including their location, size, fault geometries, and slip kinematics, is critical to estimate the impact of significant seismic events. The U.S. Geological Survey National Earthquake Information Center (NEIC) provides real‐time monitoring of earthquakes globally, including rapid source characterization and impact estimates. Here, we describe the seismic characterization products generated and made available by the NEIC over the two weeks following the start of the earthquake sequence in southeast Türkiye, their evolution, and how they inform our understanding of regional seismotectonics and hazards. The kinematics of rupture for the two earthquakes was complex, involving multiple fault segments. Therefore, incorporating observations from rupture mapping was critical for characterizing these events. Dense local datasets facilitated robust source characterization and impact assessment once these observations were obtained and converted to NEIC product input formats. We discuss how we may improve the timeliness of NEIC products for rapid assessment of future seismic hazards, particularly in the case of complex ruptures.

The Seismic Record↗

USGS2023

Substantial upper plate faulting above a shallow subduction megathrust earthquake: Mechanics and implications of the surface faulting during the 2016 Kaikoura, New Zealand, earthquake

The 2016 moment magnitude 7.8 Kaikoura, New Zealand, earthquake occurred at the southern end of the Hikurangi subduction zone where the upper plate above the shallow megathrust is exposed sub-aerially. As a result, the substantial co-seismic deformation in the upper plate above the megathrust rupture was observed geologically and geodetically. We explore the relationship between this surface faulting and the subduction megathrust rupture and find that the greatest upper plate fault slip occurred coincident (in time and location) with the megathrust rupture. Models of Coulomb stress change demonstrate that these surface faults become positively loaded as the upper plate rebounds during the megathrust event, favoring fault slip. In addition, during the megathrust rupture these faults terminate against an uncoupled subduction plate interface. We simulate the effects of decoupling at the base of these faults and find that very large fault slip is an expected consequence of this decoupling, allowing near-complete strain release. In contrast, typical strike-slip faults, pinned at their base, would have lower amounts of fault slip. These two conditions—increased Coulomb stress and basal decoupling—combine to produce the extreme co-seismic upper plate faulting observed above the shallow Kaikoura megathrust earthquake. Similar conditions occur in other global subduction zones, but in most subduction zones the region above the coupled megathrust is underwater and poorly observed. Our analysis of the Kaikoura earthquake indicates a need to reevaluate patterns of strain accumulation and release in these regions, rather than assuming simple models of elastic rebound.

Tectonics↗

USGS2022

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↗

USGS2022

Achievements and prospects of global broadband seismographic networks after 30 years of continuous geophysical observations

Global seismographic networks (GSNs) emerged during the late nineteenth and early twentieth centuries, facilitated by seminal international developments in theory, technology, instrumentation, and data exchange. The mid- to late-twentieth century saw the creation of the World-Wide Standardized Seismographic Network (1961) and International Deployment of Accelerometers (1976), which advanced global geographic coverage as seismometer bandwidth increased greatly allowing for the recording of the Earth's principal seismic spectrum. The modern era of global observations and rapid data access began during the 1980s, and notably included the inception of the GEOSCOPE initiative (1982) and GSN (1988). Through continual improvements, GEOSCOPE and the GSN have realized near-real time recording of ground motion with state-of-art data quality, dynamic range, and timing precision to encompass 180 seismic stations, many in very remote locations. Data from GSNs are increasingly integrated with other geophysical data (e.g., space geodesy, infrasound and Interferometric Synthetic Aperture Radar). Globally distributed seismic data are critical to resolving crust, mantle, and core structure; illuminating features of the plate tectonic and mantle convection system; rapid characterization of earthquakes; identification of potential tsunamis; global nuclear test verification; and provide sensitive proxies for environmental changes. As the global geosciences community continues to advance our understanding of Earth structure and processes controlling elastic wave propagation, GSN infrastructure offers a springboard to realize increasingly multi-instrument geophysical observatories. Here, we review the historical, scientific, and monitoring heritage of GSNs, summarize key discoveries, and discuss future associated opportunities for Earth Science.

Reviews of Geophysics↗

USGS2021

Seismotectonic analysis of the 2019–2020 Puerto Rico sequence: The value of absolute earthquake relocations in improved interpretations of active tectonics

We present a new catalog of calibrated earthquake relocations from the 2019–2020 Puerto Rico earthquake sequence related to the 7 January 2020 M w 6.4 earthquake that occurred offshore of southwest Puerto Rico at a depth of 15.9 km. Utilizing these relocated earthquakes and associated moment tensor solutions, we can delineate several distinct fault systems that were activated during the sequence and show that the M w 6.4 mainshock may have resulted from positive changes in Coulomb stress from earlier events. Seismicity and mechanisms define (1) a west–southwest (∼260°) zone of seismicity comprised of largely sinistral strike‐slip and oblique‐slip earthquakes that mostly occurs later in the sequence and to the west of the mainshock, (2) an area of extensional faulting that includes the mainshock and occurs largely within the mainshock’s rupture area, and (3) an north–northeast (∼30°)‐striking zone of seismicity, consisting primarily of dextral strike‐slip events that occurs before and following the mainshock and generally above (shallower than) the normal‐faulting events. These linear features intersect within the M w 6.4 mainshock’s fault plane in southwest Puerto Rico. In addition, we show that earthquake relocations for M 4+ normal‐faulting events, when traced along their fault planes, daylight along east–west‐trending bathymetric features offshore of southwest Puerto Rico. Correlation of these normal‐faulting events with bathymetric features suggests an active fault system that may be a contributor to previously uncharacterized seismic hazards in southwest Puerto Rico.

Puerto Rico↗

USGS2021

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↗

USGS2021

Short‐period surface‐wave tomography in the continental United States— A resource for research

The variation of phase and group velocity dispersion of Love and Rayleigh waves was determined for the continental United States and adjacent Canada. By processing ambient noise from the broadband channels of the Transportable Array (TA) of USArray and several Program for the Array Seismic Studies of the Continental Lithosphere experiments and using some earthquake recordings, the effort was focused on determining dispersion down to periods as short as 2 s. The relatively short distances between TA stations permitted the use of a 25    km × 25    km "> 25 km × 25 km grid for the four independent tomographic inversions (Love and Rayleigh and phase and group velocity). One reason for trying to obtain short‐period dispersion was to have a data set capable of constraining upper crust velocity models for use in determining regional moment tensors. The benefit of focusing on short‐period dispersion is apparent in the tomography maps—shallow geologic structures such as the Mid‐Continent Rift, and the Michigan, Illinois, Anadarko, Arkoma, and Appalachian basins are imaged. In our processing, we noted that the phase velocities were more robustly determined than the group velocities. We also noted that the inability to obtain dispersion at short periods shows distinct regional patterns that may be related to the local upper crust structure.

Seismological Research Letters↗

USGS2020

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↗

USGS2020

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↗

USGS2019

The August 2018 Kaktovik earthquakes: Active tectonics in northeastern Alaska revealed With InSAR and seismology

The largest earthquakes recorded in northern Alaska (M w 6.4 and M w 6.0) occurred ~6 hours apart on August 12, 2018 in the northeastern Brooks Range. The earthquakes were captured by Sentinel-1 InSAR satellites and Earthscope Transportable Array seismic data, giving insight into the little-known active tectonic processes of Arctic Alaska, obscured until recently by sparse data availability. In this study, InSAR modelling, teleseismic back projections, calibrated hypocentral relocations and regional moment tensor solutions resolve two previously unknown, SSW-dipping right-lateral fault segments. These are the first active faults identified as conjugate to the NE-trending sinistral Canning Displacement Zone directly to the west, which is therefore a more complex zone of diffuse faulting than previously thought. The northeastern Brooks Range has been characterized as an area of low to moderate seismic hazard, but these earthquakes illustrate the potential for larger, possibly destructive events in a region earmarked for rapid resource development.

Alaska↗

USGS2019

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↗

USGS2019

GLASS3: A standalone multi-scale seismic detection associator

The automated global real-time association of phase picks into seismic sources comes with unique challenges when simultaneously monitoring at local, regional and global scales. High spatial variability in seismic station density, transitory seismic data availability, and time-varying noise characteristics of individual stations must be considered in the design of an associator that is fast and accurate with a low false association rate. These challenges are particularly apparent at the U.S. Geological Survey (USGS) National Earthquake Information Center (NEIC), which monitors seismicity in near-real time on local, regional, and global scales using seismic data from roughly 2,100 real-time seismic stations. In order to fully leverage this large dataset, NEIC developed a stand-alone, self-configuring seismic phase associator, GLASS3 (GLobal ASSociator 3) that simultaneously processes variably scaled 3D association webs, each with a unique set of nucleation criteria (e.g., nucleation stack threshold). GLASS3 has many useful features for real-time monitoring including its computational efficiency, instantaneous pick processing, and on-the-fly configurability such as the creation and removal of targeted association webs and updates to supporting station metadata. GLASS3 runs both as part of a real-time event processing system, and as a configurable standalone associator that can be applied to a large variety of seismic problems. Here we describe the GLASS3 algorithm and demonstrate (including input data and configuration files) its use in associating phase-ambiguous picks on multiple scales.

BSSA↗

USGS2019

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↗

USGS2018

Spatiotemporal analysis of the Foreshock-Mainshock-Aftershock sequence of the 6 July 2017 M5.8 Lincoln, Montana, earthquake

A MW 5.8 earthquake occurred on 6 July 2017 at 12.2 km depth, 11 km southeast of Lincoln in west central Montana. No major damage or injuries were reported; however, the widely felt mainshock generated a prolific aftershock sequence with more than 1200 located events through the end of 2017. The Lincoln event is the latest in a series of moderate-to-large earthquakes that have affected western Montana. We characterize the spatiotemporal evolution of the sequence using matched filter detection and multiple-event relocation techniques. Moment tensor solutions and aftershock locations indicate faulting occurred on a 9-km-long NNE-striking, near-vertical, strike-slip fault antithetic to the Lewis and Clark Line, the main through-going fault system. Seismicity primarily occurs between 6 and 16 km depth, consistent with seismicity in the Intermountain Seismic Belt. We estimate a fault rupture area of ~64 km2 and ~30 cm of average fault displacement. We identified four foreshocks in the three days prior to, and 3005 aftershocks in the three weeks following the mainshock. The supplemented catalog frequency-magnitude distribution has a b-value of 0.79 and a minimum magnitude of completeness of 0.7. The overall decay rate is consistent with a modified Omori decay law p-value of 0.76 and c-value of 0.32. This event demonstrates that unmapped faults antithetic to major geologic structures play a role in accommodating regional strain in Western Montana and can host significant earthquakes

Montana↗