Geology topics

Justin L. Rubinstein

Publications and source records attributed to Justin L. Rubinstein.

48 records · Page 3Linked to original sources

USGS2015

Preface to the focus section on injection-induced seismicity

The ongoing, dramatic increase in seismicity in the central United States that began in 2009 is believed to be the result of injection‐induced seismicity ( Ellsworth, 2013 ). Although the basic mechanism for activation of slip on a fault by subsurface fluid injection is well established ( Healy et al. , 1968 ; Raleigh et al. , 1976 ; Nicholson and Wesson, 1992 ; McGarr et al. , 2002 ; Ellsworth, 2013 ), the occurrence of damaging M ≥5 earthquakes and the dramatic increase in seismicity in the central United States has brought heightened attention to this issue. The elevated seismicity is confined to a limited number of areas, and accumulating evidence indicates that the seismicity in these locations is directly linked to nearby industrial operations. This Seismological Research Letters ( SRL ) focus section presents a selected set of seven technical papers that cover various aspects of this topic, including basic seismological and ground‐motion observations, case studies, numerical simulation of fault activation, and risk mitigation.

Seismological Research Letters↗

USGS2015

Near-surface versus fault zone damage following the 1999 Chi-Chi earthquake: Observation and simulation of repeating earthquakes

We observe crustal damage and its subsequent recovery caused by the 1999  M 7.6 Chi-Chi earthquake in central Taiwan. Analysis of repeating earthquakes in Hualien region, ~70 km east of the Chi-Chi earthquake, shows a remarkable change in wave propagation beginning in the year 2000, revealing damage within the fault zone and distributed across the near surface. We use moving window cross correlation to identify a dramatic decrease in the waveform similarity and delays in the S wave coda. The maximum delay is up to 59 ms, corresponding to a 7.6% velocity decrease averaged over the wave propagation path. The waveform changes on either side of the fault are distinct. They occur in different parts of the waveforms, affect different frequencies, and the size of the velocity reductions is different. Using a finite difference method, we simulate the effect of postseismic changes in the wavefield by introducing S wave velocity anomaly in the fault zone and near the surface. The models that best fit the observations point to pervasive damage in the near surface and deep, along-fault damage at the time of the Chi-Chi earthquake. The footwall stations show the combined effect of near-surface and the fault zone damage, where the velocity reduction (2–7%) is twofold to threefold greater than the fault zone damage observed in the hanging wall stations. The physical models obtained here allow us to monitor the temporal evolution and recovering process of the Chi-Chi fault zone damage.

Journal of Geophysical Research B: Solid Earth↗

USGS2015

Coping with earthquakes induced by fluid injection

Large areas of the United States long considered geologically stable with little or no detected seismicity have recently become seismically active. The increase in earthquake activity began in the mid-continent starting in 2001 ( 1 ) and has continued to rise. In 2014, the rate of occurrence of earthquakes with magnitudes ( M ) of 3 and greater in Oklahoma exceeded that in California (see the figure). This elevated activity includes larger earthquakes, several with M > 5, that have caused significant damage ( 2 , 3 ). To a large extent, the increasing rate of earthquakes in the mid-continent is due to fluid-injection activities used in modern energy production ( 1 , 4 , 5 ). We explore potential avenues for mitigating effects of induced seismicity. Although the United States is our focus here, Canada, China, the UK, and others confront similar problems associated with oil and gas production, whereas quakes induced by geothermal activities affect Switzerland, Germany, and others.

Science↗

USGS2015

High-rate injection is associated with the increase in U.S. mid-continent seismicity

An unprecedented increase in earthquakes in the U.S. mid-continent began in 2009. Many of these earthquakes have been documented as induced by wastewater injection. We examine the relationship between wastewater injection and U.S. mid-continent seismicity using a newly assembled injection well database for the central and eastern United States. We find that the entire increase in earthquake rate is associated with fluid injection wells. High-rate injection wells (>300,000 barrels per month) are much more likely to be associated with earthquakes than lower-rate wells. At the scale of our study, a well’s cumulative injected volume, monthly wellhead pressure, depth, and proximity to crystalline basement do not strongly correlate with earthquake association. Managing injection rates may be a useful tool to minimize the likelihood of induced earthquakes.

Science↗

USGS2015

Increasing seismicity in the U. S. midcontinent: Implications for earthquake hazard

Earthquake activity in parts of the central United States has increased dramatically in recent years. The space-time distribution of the increased seismicity, as well as numerous published case studies, indicates that the increase is of anthropogenic origin, principally driven by injection of wastewater coproduced with oil and gas from tight formations. Enhanced oil recovery and long-term production also contribute to seismicity at a few locations. Preliminary hazard models indicate that areas experiencing the highest rate of earthquakes in 2014 have a short-term (one-year) hazard comparable to or higher than the hazard in the source region of tectonic earthquakes in the New Madrid and Charleston seismic zones.

The Leading Edge↗

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↗

USGS2015

Incorporating induced seismicity in the 2014 United States National Seismic Hazard Model: results of the 2014 workshop and sensitivity studies

The U.S. Geological Survey National Seismic Hazard Model for the conterminous United States was updated in 2014 to account for new methods, input models, and data necessary for assessing the seismic ground shaking hazard from natural (tectonic) earthquakes. The U.S. Geological Survey National Seismic Hazard Model project uses probabilistic seismic hazard analysis to quantify the rate of exceedance for earthquake ground shaking (ground motion). For the 2014 National Seismic Hazard Model assessment, the seismic hazard from potentially induced earthquakes was intentionally not considered because we had not determined how to properly treat these earthquakes for the seismic hazard analysis. The phrases “potentially induced” and “induced” are used interchangeably in this report, however it is acknowledged that this classification is based on circumstantial evidence and scientific judgment. For the 2014 National Seismic Hazard Model update, the potentially induced earthquakes were removed from the NSHM’s earthquake catalog, and the documentation states that we would consider alternative models for including induced seismicity in a future version of the National Seismic Hazard Model. As part of the process of incorporating induced seismicity into the seismic hazard model, we evaluate the sensitivity of the seismic hazard from induced seismicity to five parts of the hazard model: (1) the earthquake catalog, (2) earthquake rates, (3) earthquake locations, (4) earthquake Mmax (maximum magnitude), and (5) earthquake ground motions. We describe alternative input models for each of the five parts that represent differences in scientific opinions on induced seismicity characteristics. In this report, however, we do not weight these input models to come up with a preferred final model. Instead, we present a sensitivity study showing uniform seismic hazard maps obtained by applying the alternative input models for induced seismicity. The final model will be released after further consideration of the reliability and scientific acceptability of each alternative input model. Forecasting the seismic hazard from induced earthquakes is fundamentally different from forecasting the seismic hazard for natural, tectonic earthquakes. This is because the spatio-temporal patterns of induced earthquakes are reliant on economic forces and public policy decisions regarding extraction and injection of fluids. As such, the rates of induced earthquakes are inherently variable and nonstationary. Therefore, we only make maps based on an annual rate of exceedance rather than the 50-year rates calculated for previous U.S. Geological Survey hazard maps.

Open-File Report↗

USGS2014

Seismological and geodetic constraints on the 2011 M w 5.3 Trinidad, Colorado earthquake and induced deformation in the Raton Basin

The Raton Basin of southern Colorado and northern New Mexico is an actively produced hydrocarbon basin that has experienced increased seismicity since 2001, including the August 2011 M w 5.3 Trinidad normal faulting event. Following the 2011 earthquake, regional seismic observations were used to relocate 21 events, including the 2011 main shock, two foreshocks, and 13 aftershocks. Additionally, interferometric synthetic aperture radar (InSAR) observations of both the 2011 event and preevent basin deformation place constraint on the spatial kinematics of the 2011 event and localized basin subsidence due to ground water or gas withdrawal. We find that the 2011 earthquake ruptured an 8–10 km long segment of a normal fault at depths of 1.5–6.0 km within the crystalline Precambrian basement underlying the Raton Basin sedimentary rocks. The earthquake also nucleated within the crystalline basement in the vicinity of an active wastewater disposal site. The ensuing aftershock sequence demonstrated statistical properties expected for intraplate earthquakes, though the length of the 2011 earthquake is unexpectedly long for an M w 5.3 event, suggesting that wastewater disposal may have triggered a low stress drop, otherwise natural earthquake. Additionally, preevent and postevent seismicity in the Raton Basin spatially correlates to regions of subsidence observed in InSAR time series analysis. While these observations cannot discern a causal link between hydrocarbon production and seismicity, they constrain spatial relationships between active basin deformation and geological and anthropogenic features. Furthermore, the InSAR observations highlight the utility of space-based geodetic observations for monitoring and assessing anthropogenically induced and triggered deformation.

Journal of Geophysical Research B: Solid Earth↗

USGS2014

The 2001-present induced earthquake sequence in the Raton Basin of northern New Mexico and southern Colorado

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.

Colorado;New Mexico↗

USGS2012

Fixed recurrence and slip models better predict earthquake behavior than the time- and slip-predictable models 1: repeating earthquakes

The behavior of individual events in repeating earthquake sequences in California, Taiwan and Japan is better predicted by a model with fixed inter-event time or fixed slip than it is by the time- and slip-predictable models for earthquake occurrence. Given that repeating earthquakes are highly regular in both inter-event time and seismic moment, the time- and slip-predictable models seem ideally suited to explain their behavior. Taken together with evidence from the companion manuscript that shows similar results for laboratory experiments we conclude that the short-term predictions of the time- and slip-predictable models should be rejected in favor of earthquake models that assume either fixed slip or fixed recurrence interval. This implies that the elastic rebound model underlying the time- and slip-predictable models offers no additional value in describing earthquake behavior in an event-to-event sense, but its value in a long-term sense cannot be determined. These models likely fail because they rely on assumptions that oversimplify the earthquake cycle. We note that the time and slip of these events is predicted quite well by fixed slip and fixed recurrence models, so in some sense they are time- and slip-predictable. While fixed recurrence and slip models better predict repeating earthquake behavior than the time- and slip-predictable models, we observe a correlation between slip and the preceding recurrence time for many repeating earthquake sequences in Parkfield, California. This correlation is not found in other regions, and the sequences with the correlative slip-predictable behavior are not distinguishable from nearby earthquake sequences that do not exhibit this behavior.

Journal of Geophysical Research B: Solid Earth↗

USGS2011

Nonlinear site response in medium magnitude earthquakes near Parkfield, California

Careful analysis of strong-motion recordings of 13 medium magnitude earthquakes (3.7 ≤ M ≤ 6.5) in the Parkfield, California, area shows that very modest levels of shaking (approximately 3.5% of the acceleration of gravity) can produce observable changes in site response. Specifically, I observe a drop and subsequent recovery of the resonant frequency at sites that are part of the USGS Parkfield dense seismograph array (UPSAR) and Turkey Flat array. While further work is necessary to fully eliminate other models, given that these frequency shifts correlate with the strength of shaking at the Turkey Flat array and only appear for the strongest shaking levels at UPSAR, the most plausible explanation for them is that they are a result of nonlinear site response. Assuming this to be true, the observation of nonlinear site response in small ( M <5) earthquakes implies that nonlinear site response can occur at much lower levels of shaking than previously believed. Below I present observations of a resonant frequency shift during five M ≤ 5 earthquakes near Parkfield, California, strongly contrasting with previous studies that have only identified nonlinear site effects for much larger events. In addition to the nonlinear effects seen for the smaller events, nonlinear site response is also observed for two largest earthquakes in the region during the study period (the 2003 M 6.5 San Simeon earthquake and the 2004 M 6 Parkfield earthquake).

Bulletin of the Seismological Society of America↗

USGS2009

Seismic wave triggering of nonvolcanic tremor, episodic tremor and slip, and earthquakes on Vancouver Island

[1] We explore the physical conditions that enable triggering of nonvolcanic tremor and earthquakes by considering local seismic activity on Vancouver Island, British Columbia during and immediately after the arrival of large-amplitude seismic waves from 30 teleseismic and 17 regional or local earthquakes. We identify tremor triggered by four of the teleseismic earthquakes. The close temporal and spatial proximity of triggered tremor to ambient tremor and aseismic slip indicates that when a fault is close to or undergoing failure, it is particularly susceptible to triggering of further events. The amplitude of the triggering waves also influences the likelihood of triggering both tremor and earthquakes such that large amplitude waves triggered tremor in the absence of detectable aseismic slip or ambient tremor. Tremor and energy radiated from regional/local earthquakes share the same frequency passband so that tremor cannot be identified during these smaller, more frequent events. We confidently identify triggered local earthquakes following only one teleseism, that with the largest amplitude, and four regional or local events that generated vigorous aftershock sequences in their immediate vicinity. Earthquakes tend to be triggered in regions different from tremor and with high ambient seismicity rates. We also note an interesting possible correlation between large teleseismic events and episodic tremor and slip (ETS) episodes, whereby ETS events that are “late” and have built up more stress than normal are susceptible to triggering by the slight nudge of the shaking from a large, distant event, while ETS events that are “early” or “on time” are not.

Vancouver Island↗