Geology topics

Harley M. Benz

Publications and source records attributed to Harley M. Benz.

At least 37 records · Page 2Linked to original sources

USGS2018

Seismicity in the Challis, Idaho region, January 2014 - May 2017: Late aftershocks of the 1983 Ms 7.3 Borah Peak earthquake

In April 2014, after about 20 yrs of relatively low seismicity, an energetic earthquake sequence (maximum M L "> M L 4.8) began 25–30 km northwest of the 1983 M s "> M s 7.3 Borah Peak earthquake rupture area near the town of Challis, Idaho. This sequence ended in the fall of 2014, but in January 2015, a second energetic sequence (maximum M L "> M L 5.0) began about 20 km to the southeast. Modest seismicity has continued in both regions with ∼ 1000 "> ∼ 1000 earthquakes detected and located through May 2017. To better characterize the seismicity in the area, we deployed a seven‐station local seismometer network during April–October 2014; one of the stations remained active until July 2015. Here, we report updated locations for earthquakes in the Challis area for 1 January 2014–31 May 2017. Using a combination of absolute and differential arrival times, we generated a catalog of high‐accuracy relocations. The earthquakes clustered into four primary groups, three of them with strikes similar to the Lost River fault—the fault responsible for the 1983 Borah Peak event. We used a modified cut‐and‐paste method to determine moment tensors for 15 of the largest events. All of the moment tensors showed normal faulting with nodal plane strikes consistent with the trend of the relocated seismicity and the regional stress field. We suggest that the recent seismicity near Challis is best interpreted as a continuation of the 1983 M s "> M s Ms 7.3 Borah Peak aftershock sequence, which is unusually long compared to plate boundary aftershock sequences because of the lower regional strain rate.

Idaho↗

USGS2017

2017 Valparaíso earthquake sequence and the megathrust patchwork of central Chile

In April 2017, a sequence of earthquakes offshore Valparaíso, Chile, raised concerns of a potential megathrust earthquake in the near future. The largest event in the 2017 sequence was a M 6.9 on 24 April, seemingly colocated with the last great-sized earthquake in the region—a M 8.0 in March 1985. The history of large earthquakes in this region shows significant variation in rupture size and extent, typically highlighted by a juxtaposition of large ruptures interspersed with smaller magnitude sequences. We show that the 2017 sequence ruptured an area between the two main slip patches during the 1985 earthquake, rerupturing a patch that had previously slipped during the October 1973 M 6.5 earthquake sequence. A significant gap in historic ruptures exists directly to the south of the 2017 sequence, with large enough moment deficit to host a great-sized earthquake in the near future, if it is locked.

Valparaiso↗

USGS2017

Spatio-temporal evolution of the 2011 Prague, Oklahoma aftershock sequence revealed using subspace detection and relocation

The 6 November 2011 M w 5.7 earthquake near Prague, Oklahoma is the second largest earthquake ever recorded in the state. A M w 4.8 foreshock and the M w 5.7 mainshock triggered a prolific aftershock sequence. Utilizing a subspace detection method, we increase by fivefold the number of precisely located events between 4 November and 5 December 2011. We find that while most aftershock energy is released in the crystalline basement, a significant number of the events occur in the overlying Arbuckle Group, indicating that active Meeker-Prague faulting extends into the sedimentary zone of wastewater disposal. Although the number of aftershocks in the Arbuckle Group is large, comprising ~40% of the aftershock catalog, the moment contribution of Arbuckle Group earthquakes is much less than 1% of the total aftershock moment budget. Aftershock locations are sparse in patches that experienced large slip during the mainshock.

Geophysical Research Letters↗

USGS2017

Integrated geophysical characteristics of the 2015 Illapel, Chile, earthquake

On 16 September 2015, a Mw 8.3 earthquake ruptured the subduction zone offshore of Illapel, Chile, generating an aftershock sequence with 14 Mw 6.0–7.0 events. A double source W phase moment tensor inversion consists of a Mw 7.2 subevent and the main Mw 8.2 phase. We determine two slip models for the mainshock, one using teleseismic broadband waveforms and the other using static GPS and InSAR surface displacements, which indicate high slip north of the epicenter and west-northwest of the epicenter near the oceanic trench. These models and slip distributions published in other studies suggest spatial slip uncertainties of ~25 km and have peak slip values that vary by a factor of 2. We relocate aftershock hypocenters using a Bayesian multiple-event relocation algorithm, revealing a cluster of aftershocks under the Chilean coast associated with deep (20–45 km depth) mainshock slip. Less vigorous aftershock activity also occurred near the trench and along strike of the main aftershock region. Most aftershocks are thrust-faulting events, except for normal-faulting events near the trench. Coulomb failure stress change amplitudes and signs are uncertain for aftershocks collocated with deeper mainshock slip; other aftershocks are more clearly associated with loading from the mainshock. These observations reveal a frictionally heterogeneous interface that ruptured in patches at seismogenic depths (associated with many aftershocks) and with homogeneous slip (and few aftershocks) up to the trench. This event likely triggered seismicity separate from the main slip region, including along-strike events on the megathrust and intraplate extensional events.

Illapel↗

USGS2017

The 2008 Wells, Nevada earthquake sequence: Source constraints using calibrated multiple event relocation and InSAR

The 2008 Wells, NV earthquake represents the largest domestic event in the conterminous U.S. outside of California since the October 1983 Borah Peak earthquake in southern Idaho. We present an improved catalog, magnitude complete to 1.6, of the foreshock-aftershock sequence, supplementing the current U.S. Geological Survey (USGS) Preliminary Determination of Epicenters (PDE) catalog with 1,928 well-located events. In order to create this catalog, both subspace and kurtosis detectors are used to obtain an initial set of earthquakes and associated locations. The latter are then calibrated through the implementation of the hypocentroidal decomposition method and relocated using the BayesLoc relocation technique. We additionally perform a finite fault slip analysis of the mainshock using InSAR observations. By combining the relocated sequence with the finite fault analysis, we show that the aftershocks occur primarily updip and along the southwestern edge of the zone of maximum slip. The aftershock locations illuminate areas of post-mainshock strain increase; aftershock depths, ranging from 5 to 16 km, are consistent with InSAR imaging, which shows that the Wells earthquake was a buried source with no observable near-surface offset.

Nevada↗

USGS2017

Oklahoma experiences largest earthquake during ongoing regional wastewater injection hazard mitigation efforts

The 3 September 2016, M w 5.8 Pawnee earthquake was the largest recorded earthquake in the state of Oklahoma. Seismic and geodetic observations of the Pawnee sequence, including precise hypocenter locations and moment tensor modeling, shows that the Pawnee earthquake occurred on a previously unknown left-lateral strike-slip basement fault that intersects the mapped right-lateral Labette fault zone. The Pawnee earthquake is part of an unprecedented increase in the earthquake rate in Oklahoma that is largely considered the result of the deep injection of waste fluids from oil and gas production. If this is, indeed, the case for the M 5.8 Pawnee earthquake, then this would be the largest event to have been induced by fluid injection. Since 2015, Oklahoma has undergone wide-scale mitigation efforts primarily aimed at reducing injection volumes. Thus far in 2016, the rate of M 3 and greater earthquakes has decreased as compared to 2015, while the cumulative moment—or energy released from earthquakes—has increased. This highlights the difficulty in earthquake hazard mitigation efforts given the poorly understood long-term diffusive effects of wastewater injection and their connection to seismicity.

Oklahoma↗

USGS2017

Tectonic summaries of magnitude 7 and greater earthquakes from 2000 to 2015

This paper describes the tectonic summaries for all magnitude 7 and larger earthquakes in the period 2000–2015, as produced by the U.S. Geological Survey National Earthquake Information Center during their routine response operations to global earthquakes. The goal of such summaries is to provide important event-specific information to the public rapidly and concisely, such that recent earthquakes can be understood within a global and regional seismotectonic framework. We compile these summaries here to provide a long-term archive for this information, and so that the variability in tectonic setting and earthquake history from region to region, and sometimes within a given region, can be more clearly understood.

Open-File Report↗

USGS2017

Source modeling of the 2015 Mw 7.8 Nepal (Gorkha) earthquake sequence: Implications for geodynamics and earthquake hazards

The Gorkha earthquake on April 25th, 2015 was a long anticipated, low-angle thrust-faulting event on the shallow décollement between the India and Eurasia plates. We present a detailed multiple-event hypocenter relocation analysis of the Mw 7.8 Gorkha Nepal earthquake sequence, constrained by local seismic stations, and a geodetic rupture model based on InSAR and GPS data. We integrate these observations to place the Gorkha earthquake sequence into a seismotectonic context and evaluate potential earthquake hazard. Major results from this study include (1) a comprehensive catalog of calibrated hypocenters for the Gorkha earthquake sequence; (2) the Gorkha earthquake ruptured a ~ 150 × 60 km patch of the Main Himalayan Thrust (MHT), the décollement defining the plate boundary at depth, over an area surrounding but predominantly north of the capital city of Kathmandu (3) the distribution of aftershock seismicity surrounds the mainshock maximum slip patch; (4) aftershocks occur at or below the mainshock rupture plane with depths generally increasing to the north beneath the higher Himalaya, possibly outlining a 10–15 km thick subduction channel between the overriding Eurasian and subducting Indian plates; (5) the largest Mw 7.3 aftershock and the highest concentration of aftershocks occurred to the southeast the mainshock rupture, on a segment of the MHT décollement that was positively stressed towards failure; (6) the near surface portion of the MHT south of Kathmandu shows no aftershocks or slip during the mainshock. Results from this study characterize the details of the Gorkha earthquake sequence and provide constraints on where earthquake hazard remains high, and thus where future, damaging earthquakes may occur in this densely populated region. Up-dip segments of the MHT should be considered to be high hazard for future damaging earthquakes.

Gorkha↗

USGS2016

Far-field pressurization likely caused one of the largest injection induced earthquakes by reactivating a large pre-existing basement fault structure

The M w 5.1 Fairview, Oklahoma, earthquake on 13 February 2016 and its associated seismicity produced the largest moment release in the central and eastern United States since the 2011 M w 5.7 Prague, Oklahoma, earthquake sequence and is one of the largest earthquakes potentially linked to wastewater injection. This energetic sequence has produced five earthquakes with M w 4.4 or larger. Almost all of these earthquakes occur in Precambrian basement on a partially unmapped 14 km long fault. Regional injection into the Arbuckle Group increased approximately sevenfold in the 36 months prior to the start of the sequence (January 2015). We suggest far-field pressurization from clustered, high-rate wells greater than 12 km from this sequence induced these earthquakes. As compared to the Fairview sequence, seismicity is diffuse near high-rate wells, where pressure changes are expected to be largest. This points to the critical role that preexisting faults play in the occurrence of large induced earthquakes.

Geophysical Research Letters↗

USGS2016

Hydra—The National Earthquake Information Center’s 24/7 seismic monitoring, analysis, catalog production, quality analysis, and special studies tool suite

This report provides an overview of the capabilities and design of Hydra, the global seismic monitoring and analysis system used for earthquake response and catalog production at the U.S. Geological Survey National Earthquake Information Center (NEIC). Hydra supports the NEIC’s worldwide earthquake monitoring mission in areas such as seismic event detection, seismic data insertion and storage, seismic data processing and analysis, and seismic data output. The Hydra system automatically identifies seismic phase arrival times and detects the occurrence of earthquakes in near-real time. The system integrates and inserts parametric and waveform seismic data into discrete events in a database for analysis. Hydra computes seismic event parameters, including locations, multiple magnitudes, moment tensors, and depth estimates. Hydra supports the NEIC’s 24/7 analyst staff with a suite of seismic analysis graphical user interfaces. In addition to the NEIC’s monitoring needs, the system supports the processing of aftershock and temporary deployment data, and supports the NEIC’s quality assurance procedures. The Hydra system continues to be developed to expand its seismic analysis and monitoring capabilities.

Open-File Report↗

USGS2016

Rapid response, monitoring, and mitigation of induced seismicity near Greeley, Colorado

On 1 June 2014 (03:35 UTC), an M w 3.2 earthquake occurred in Weld County, Colorado, a historically aseismic area of the Denver–Julesburg basin. Weld County is a prominent area of oil and gas development, including many high‐rate class II wastewater injection wells. In the days following the earthquake, the University of Colorado, with support from the U.S. Geological Survey and Incorporated Research Institutions for Seismology–Portable Array Seismic Studies of the Continental Lithosphere, rapidly deployed six seismic stations to characterize the seismicity associated with the 1 June earthquake (the Greeley sequence) and to investigate its possible connection to wastewater disposal. The spatial and temporal proximity of earthquakes to a high‐rate wastewater disposal well strongly suggests these earthquakes were induced. Scientific communication between the university, state agencies, and the energy industry led to rapid mitigation strategies to reduce the occurrence of further earthquakes. Mitigation efforts included implementing a temporary moratorium on injection at the well, cementing the bottom portion of the disposal well to minimize hydrologic connectivity between the disposal formation and the underlying crystalline basement, and subsequently allowing injection to resume at lower rates. Following the resumption of wastewater disposal, microseismicity was closely monitored for both increases in earthquake rate and magnitude. Following mitigation efforts, between 13 August 2014 and 29 December 2015, no earthquakes larger than M 1.5 occurred near the Greeley sequence. This study demonstrates that a detailed and rapid characterization of a seismic sequence in space and time relative to disposal, combined with collaboration and communication between scientists, regulators, and industry, can lead to objective and actionable mitigation efforts that potentially reduced the rate of earthquakes and the possible generation of larger earthquakes.

Colorado↗

USGS2016

RMT focal plane sensitivity to seismic network geometry and faulting style

Modern tectonic studies often use regional moment tensors (RMTs) to interpret the seismotectonic framework of an earthquake or earthquake sequence; however, despite extensive use, little existing work addresses RMT parameter uncertainty. Here, we quantify how network geometry and faulting style affect RMT sensitivity. We examine how data-model fits change with fault plane geometry (strike and dip) for varying station configurations. We calculate the relative data fit for incrementally varying geometries about a best-fitting solution, applying our workflow to real and synthetic seismograms for both real and hypothetical station distributions and earthquakes. Initially, we conduct purely observational tests, computing RMTs from synthetic seismograms for hypothetical earthquakes and a series of well-behaved network geometries. We then incorporate real data and station distributions from the International Maule Aftershock Deployment (IMAD), which recorded aftershocks of the 2010 M W 8.8 Maule earthquake, and a set of regional stations capturing the ongoing earthquake sequence in Oklahoma and southern Kansas. We consider RMTs computed under three scenarios: (1) real seismic records selected for high data quality; (2) synthetic seismic records with noise computed for the observed source-station pairings and (3) synthetic seismic records with noise computed for all possible station-source pairings. To assess RMT sensitivity for each test, we observe the ‘fit falloff’, which portrays how relative fit changes when strike or dip varies incrementally; we then derive the ranges of acceptable strikes and dips by identifying the span of solutions with relative fits larger than 90 per cent of the best fit. For the azimuthally incomplete IMAD network, Scenario 3 best constrains fault geometry, with average ranges of 45° and 31° for strike and dip, respectively. In Oklahoma, Scenario 3 best constrains fault dip with an average range of 46°; however, strike is best constrained by Scenario 1, with a range of 26°. We draw two main conclusions from this study. (1) Station distribution impacts our ability to constrain RMTs using waveform time-series; however, in some tectonic settings, faulting style also plays a significant role and (2) increasing station density and data quantity (both the number of stations and the number of individual channels) does not necessarily improve RMT constraint. These results may be useful when organizing future seismic deployments (e.g. by concentrating stations in alignment with anticipated nodal planes), and in computing RMTs, either by guiding a more rigorous data selection process for input data or informing variable weighting among the selected data (e.g. by eliminating the transverse component when strike-slip mechanisms are expected).

Geophysical Journal International↗

USGS2015

Rapid characterization of the 2015 M w 7.8 Gorkha, Nepal, earthquake sequence and its seismotectonic context

Earthquake response and related information products are important for placing recent seismic events into context and particularly for understanding the impact earthquakes can have on the regional community and its infrastructure. These tools are even more useful if they are available quickly, ahead of detailed information from the areas affected by such earthquakes. Here we provide an overview of the response activities and related information products generated and provided by the U.S. Geological Survey National Earthquake Information Center in association with the 2015 M 7.8 Gorkha, Nepal, earthquake. This group monitors global earthquakes 24  hrs/day and 7  days/week to provide rapid information on the location and size of recent events and to characterize the source properties, tectonic setting, and potential fatalities and economic losses associated with significant earthquakes. We present the timeline over which these products became available, discuss what they tell us about the seismotectonics of the Gorkha earthquake and its aftershocks, and examine how their information is used today, and might be used in the future, to help mitigate the impact of such natural disasters.

Seismological Research Letters↗

USGS2015

Reactivated faulting near Cushing, Oklahoma: Increased potential for a triggered earthquake in an area of United States strategic infrastructure

In October 2014 two moderate-sized earthquakes ( Mw 4.0 and 4.3) struck south of Cushing, Oklahoma, below the largest crude oil storage facility in the world. Combined analysis of the spatial distribution of earthquakes and regional moment tensor focal mechanisms indicate reactivation of a subsurface unnamed and unmapped left-lateral strike-slip fault. Coulomb failure stress change calculations using the relocated seismicity and slip distribution determined from regional moment tensors, allow for the possibility that the Wilzetta-Whitetail fault zone south of Cushing, Oklahoma, could produce a large, damaging earthquake comparable to the 2011 Prague event. Resultant very strong shaking levels (MMI VII) in the epicentral region present the possibility of this potential earthquake causing moderate to heavy damage to national strategic infrastructure and local communities.

Oklahoma↗

USGS2015

Seismicity of the Earth 1900‒2013 Mediterranean Sea and vicinity

The Mediterranean region is seismically active due to the convergence of the Africa Plate with the Eurasia plate. Present day Africa-Eurasia motion ranges from ~4 millimeters per year (mm/yr) in a northwest-southeast direction in the western Mediterranean to ~10 mm/yr (north-south) in the eastern Mediterranean. The Africa-Eurasia plate boundary is complex, and includes extensional and translational zones in addition to the dominant convergent regimes characterized by subduction and continental collision. This convergence began at approximately 50 million years ago and was associated with the closure of the Tethys Sea; the Mediterranean Sea is all that remains of the Tethys. The highest rates of seismicity in the Mediterranean region are found along the Hellenic subduction zone of southern Greece and the North Anatolian Fault Zone of northwestern Turkey, but significant rates of current seismicity and large historical earthquakes have occurred throughout the region spanning the Mediterranean Sea.

Open-File Report↗

USGS2015

Hundreds of earthquakes per day: The 2014 Guthrie, Oklahoma, Earthquake Sequence

A remarkable increase in seismic activity in Oklahoma since 2009 has been shown to correlate closely with enhanced hydrocarbon extraction and associated wastewater disposal; 99% of this recent Oklahoma earthquake activity has occurred within 15 km of a call II injection well (Ellsworth, 2013). In response to this increase in seismic activity, the U.S. Geological Survey (USGS) partnered with the Oklahoma Geological Survey (OGS) to exchange waveform data from permanent and temporary seismic stations to improve the cataloging of earthquake source parameters for a broad region of north-central Oklahoma. For a particularly persistent earthquake sequence near Guthrie, Oklahoma, a subspace detection method is applied to data from nearby seismic stations. This approach documents the occurrence of hundreds of readily detectable, highly similar, earthquakes per day, with rates occasionally exceeding 1000 earthquakes per day. Time-varying changes in b-value appear episodic, suggesting a correlation with periods of reversible fault weakening and associated failure.

Oklahoman↗

USGS2015

Earthquake hypocenters and focal mechanisms in central Oklahoma reveal a complex system of reactivated subsurface strike-slip faulting

The sharp increase in seismicity over a broad region of central Oklahoma has raised concern regarding the source of the activity and its potential hazard to local communities and energy industry infrastructure. Since early 2010, numerous organizations have deployed temporary portable seismic stations in central Oklahoma in order to record the evolving seismicity. In this study, we apply a multiple-event relocation method to produce a catalog of 3,639 central Oklahoma earthquakes from late 2009 through 2014. RMT source parameters were determined for 195 of the largest and best-recorded earthquakes. Combining RMT results with relocated seismicity enabled us to determine the length, depth and style-of-faulting occurring on reactivated subsurface fault systems. Results show that the majority of earthquakes occur on near vertical, optimally oriented (NE-SW and NW-SE), strike-slip faults in the shallow crystalline basement. These are necessary first order observations required to assess the potential hazards of individual faults in Oklahoma.

Oklahoma↗

USGS2015

Efforts to monitor and characterize the recent increasing seismicity in central Oklahoma

The sharp increase in seismicity over a broad region of central Oklahoma has raised concerns regarding the source of the activity and its potential hazard to local communities and energy-industry infrastructure. Efforts to monitor and characterize the earthquake sequences in central Oklahoma are reviewed. Since early 2010, numerous organizations have deployed temporary portable seismic stations in central Oklahoma to record the evolving seismicity. A multiple-event relocation method is applied to produce a catalog of central Oklahoma earthquakes from late 2009 into early 2015. Regional moment tensor (RMT) source parameters were determined for the largest and best-recorded earthquakes. Combining RMT results with relocated seismicity enabled determination of the length, depth, and style of faulting occurring on reactivated subsurface fault systems. It was found that the majority of earthquakes occur on near-vertical, optimally oriented (northeast-southwest and northwest-southeast) strike-slip faults in the shallow crystalline basement. In 2014, 17 earthquakes occurred with magnitudes of 4 or larger. It is suggested that these recently reactivated fault systems pose the greatest potential hazard to the region.

Oklahoma↗