ALBERT

Description: 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.

Description: 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.

Description: 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.

Description: After the 2010 M w 8.8 Maule earthquake, an international collaboration involving teams and instruments from Chile, the US, the UK, France and Germany established the International Maule Aftershock Deployment temporary network over the source region of the event to facilitate detailed, open-access studies of the aftershock sequence. Using data from the first 9-months of this deployment, we have analyzed the detailed spatial distribution of over 2500 well-recorded aftershocks. All earthquakes have been relocated using a hypocentral decomposition algorithm to study the details of and uncertainties in both their relative and absolute locations. We have computed regional moment tensor solutions for the largest of these events to produce a catalogue of 465 mechanisms, and have used all of these data to study the spatial distribution of the aftershock sequence with respect to the Chilean megathrust. We refine models of co-seismic slip distribution of the Maule earthquake, and show how small changes in fault geometries assumed in teleseismic finite fault modelling significantly improve fits to regional GPS data, implying that the accuracy of rapid teleseismic fault models can be substantially improved by consideration of existing fault geometry model databases. We interpret all of these data in an integrated seismotectonic framework for the Maule earthquake rupture and its aftershock sequence, and discuss the relationships between co-seismic rupture and aftershock distributions. While the majority of aftershocks are interplate thrust events located away from regions of maximum co-seismic slip, interesting clusters of aftershocks are identified in the lower plate at both ends of the main shock rupture, implying internal deformation of the slab in response to large slip on the plate boundary interface. We also perform Coulomb stress transfer calculations to compare aftershock locations and mechanisms to static stress changes following the Maule rupture. Without the incorporation of uncertainties in earthquake locations, just 55 per cent of aftershock nodal planes align with faults promoted towards failure by co-seismic slip. When epicentral uncertainties are considered (on the order of just ±2–3 km), 90 per cent of aftershocks are consistent with occurring along faults demonstrating positive stress transfer. These results imply large sensitivities of Coulomb stress transfer calculations to uncertainties in both earthquake locations and models of slip distributions, particularly when applied to aftershocks close to a heterogeneous fault rupture; such uncertainties should therefore be considered in similar studies used to argue for or against models of static stress triggering.

Description: We investigate the ongoing seismicity in the Raton Basin and find that the deep injection of wastewater from the coal-bed methane field is responsible for inducing the majority of the seismicity since 2001. Many lines of evidence indicate that this earthquake sequence was induced by wastewater injection. First, there was a marked increase in seismicity shortly after major fluid injection began in the Raton Basin in 1999. From 1972 through July 2001, there was one M≥4 earthquake in the Raton Basin, whereas 12 occurred between August 2001 and 2013. The statistical likelihood that such a rate change would occur if earthquakes behaved randomly in time is 3.0%. Moreover, this rate change is limited to the area of industrial activity. Earthquake rates remain low in the surrounding area. Second, the vast majority of the seismicity is within 5 km of active disposal wells and is shallow, ranging between 2 and 8 km depth. The two most carefully studied earthquake sequences in 2001 and 2011 have earthquakes within 2 km of high-volume, high-injection-rate wells. Third, injection wells in the area are commonly very high volume and high rate. Two wells adjacent to the August 2011 M  5.3 earthquake injected about 4.9 million cubic meters of wastewater before the earthquake, more than seven times the amount injected at the Rocky Mountain Arsenal well that caused damaging earthquakes near Denver, Colorado, in the 1960s. The August 2011 M  5.3 event is the second-largest earthquake to date for which there is clear evidence that the earthquake sequence was induced by fluid injection. Online Material: Gutenberg–Richter plots for varying decade-long catalogs.

Description: The U.S. Geological Survey National Earthquake Information Center (NEIC) uses a variety of classical network-averaged magnitudes (e.g., m b and M s ) and waveform modeling procedures to determine the moment magnitude ( M w ) of an earthquake from teleseismic observations. Initial magnitude estimates are often inaccurate because of poor azimuthal control (sampling of the focal sphere) and/or intrinsic limitation of each method to a specific range of event size. To provide faster and more accurate estimates of the moment magnitude, source duration, and source complexity, NEIC is exploring the use of a variation of the empirical Green’s function (EGF) deconvolution procedure. This approach uses a predicted focal mechanism derived from the Global Centroid Moment Tensor Catalog to compute teleseismic P -wave synthetic seismograms, which are then deconvolved from observed P and SH waveforms to determine station-specific M w , source time function, and a network-averaged M w . Our EGF approach is validated using broadband waveforms from 246 earthquakes in the magnitude range M w  6.0–9.1. Within approximately 13 min of earthquake origin time, our procedure using teleseismic P waves only computes an M w that lies within ±0.25 of the final W -phase M w in the magnitude range 6–8. Using later arriving teleseismic SH phases results in an M w that lies within ±0.12 of the W -phase M w . For magnitude 8 or larger earthquakes, we underestimated the moment magnitude by up to 0.8 magnitude units, primarily due to the initial P phase not containing the total seismic moment release. Long-period phases such as the W -phase and surface waves that better characterize total moment release can also be incorporated in the processing.

Description: 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).

Description: The 2008 Wells, Nevada, earthquake represents the largest domestic event in the contiguous United States 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 Preliminary Determination of Epicenters catalog with 1928 well-located events. 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 Interferometric Synthetic Aperture Radar (InSAR) observations. By combining the relocated sequence with the finite-fault analysis, we show that the aftershocks occur primarily up-dip and along the southwestern edge of the zone of maximum slip. The aftershock locations illuminate areas of postmainshock 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. Electronic Supplement: Description of methods and figures of relocated aftershocks, interferograms, residuals, transformations, and seismic station locations.

Description: The M w  5.8 earthquake of 23 August 2011 (17:51:04 UTC) (moment, M 0 5.7 x 10 17 N·m) occurred near Mineral, Virginia, within the central Virginia seismic zone and was felt by more people than any other earthquake in United States history. The U.S. Geological Survey (USGS) received 148,638 felt reports from 31 states and 4 Canadian provinces. The USGS PAGER system estimates as many as 120,000 people were exposed to shaking intensity levels of IV and greater, with approximately 10,000 exposed to shaking as high as intensity VIII. Both regional and teleseismic moment tensor solutions characterize the earthquake as a northeast-striking reverse fault that nucleated at a depth of approximately 7±2 km. The distribution of reported macroseismic intensities is roughly ten times the area of a similarly sized earthquake in the western United States ( Horton and Williams, 2012 ). Near-source and far-field damage reports, which extend as far away as Washington, D.C., (135 km away) and Baltimore, Maryland, (200 km away) are consistent with an earthquake of this size and depth in the eastern United States (EUS). Within the first few days following the earthquake, several government and academic institutions installed 36 portable seismograph stations in the epicentral region, making this among the best-recorded aftershock sequences in the EUS. Based on modeling of these data, we provide a detailed description of the source parameters of the mainshock and analysis of the subsequent aftershock sequence for defining the fault geometry, area of rupture, and observations of the aftershock sequence magnitude–frequency and temporal distribution. The observed slope of the magnitude–frequency curve or b -value for the aftershock sequence is consistent with previous EUS studies ( b =0.75), suggesting that most of the accumulated strain was released by the mainshock. The aftershocks define a rupture that extends between approximately 2–8 km in depth and 8–10 km along the strike of the fault plane. Best-fit modeling of the geometry of the aftershock sequence defines a rupture plane that strikes N36°E and dips to the east-southeast at 49.5°. Moment tensor solutions of the mainshock and larger aftershocks are consistent with the distribution of aftershock locations, both indicating reverse slip along a northeast–southwest striking southeast-dipping fault plane. Online Material: Tables of regional moment tensor source parameters and aftershock location.