Geology ReportsSearch

USGS · 70203088

Reservoir analysis of the Denver earthquakes: A case of induced seismicity

Abstract

Injection of fluid wastes into the fractured Precambrian crystalline bedrock beneath the Rocky Mountain Arsenal near Denver triggered earthquakes in the 1960's. An analysis, based on the assumption that fluid flow in the fractured reservoir can be approximated by flow in a porous medium, is presented. The configuration and hydrologic properties of the reservoir are determined from two lines of evidence: (1) locations of earthquake hypocenters determined by seismic arrays installed at the Arsenal and (2) observed long-term decline in fluid levels in the injection well. Together these two sets of data indicate that a long, narrow reservoir, aligned in the direction N 60°W, exists. The reservoir is 3.35 km in width, extends 30.5 km to the northwest and infinitely to the southeast, and spans a depth interval from 3.7 to 7.0 km below land surface. It has a transmissivity of 1.08×10-5 m2/s and a storage coefficient of 1.0×10-5. Computed pressure buildup along the length of the reservoir is compared with the spatial distribution of earthquake epicenters. The comparison shows that earthquakes are confined to that part of the reservoir where the pressure buildup exceeds 32 bars. This critical value is interpreted as the pressure buildup above which earthquakes occur. The migration of earthquake epicenters away from the injection well, a phenomenon noted by previous investigators, can be accounted for by the outward propagation of the critical pressure buildup. The analysis is extended to examining the effects of rapid flow in fractures opened by high injection pressure. The results show that the effect is confined to a small region within 1 km of the injection well. The existence of a critical pressure buildup above which earthquakes occur is completely consistent with the theory on the role of fluid pressure in fault movement as presented by Hubbert and Rubey.

Explore related subjects

90° N90° S · 180° W ← longitude → 180° E
Source-reported bounding extent: 39.00851330385611° to 40.5930995321649° latitude; -105.99609375° to -103.73291015625° longitude. This indicates report coverage, not an exact sampling location. View area on OpenStreetMap.

Keep this discovery

Explore connections, maps & timelines

BibTeXRIS

Paul A. Hsieh, John D. Bredehoeft. 2012-09-20. Reservoir analysis of the Denver earthquakes: A case of induced seismicity. https://doi.org/10.1029/jb086ib02p00903

Cite the original work for its findings. Save a collection to share your selection of sources.

KEEP EXPLORING

Related USGS reports

Assessing the seasonal dynamics of nitrate and sulfate aerosols at the South Pole utilizing stable isotopes

Atmospheric nitrate (NO 3 − = particulate NO 3 − + gas‐phase nitric acid [HNO 3 ]) and sulfate (SO 4 2− ) are key molecules that play important roles in numerous atmospheric processes. Here, the seasonal cycles of NO 3 − and total suspended particulate sulfate (SO 4 2− (TSP) ) were evaluated at the South Pole from aerosol samples collected weekly for approximately 10 months (26 January to 25 October) in 2002 and analyzed for their concentration and isotopic compositions. Aerosol NO 3 − was largely affected by snowpack emissions in which [NO 3 − ] and δ 15 N(NO 3 − ) were highest (49.3 ± 21.4 ng/m 3 , n = 8) and lowest (−47.0 ± 11.7‰, n = 5), respectively, during periods of sunlight in the interior of Antarctica. The seasonal cycle of Δ 17 O(NO 3 − ) reflected tropospheric chemistry year‐round with lower values observed during sunlight periods and higher values observed during dark periods, reflecting shifts from HO x ‐ to O 3 ‐dominated oxidation chemistry. SO 4 2− (TSP) concentrations were highest during austral summer and fall (86.7 ± 73.7 ng/m 3 , n = 18) and are indicated to be derived from dimethyl sulfide (DMS) emissions, as δ 34 S(SO 4 2− ) (TSP) values (18.5 ± 1.0‰, n = 10) were similar to literature δ 34 S(DMS) values. The seasonal cycle of Δ 17 O(SO 4 2− ) (TSP) exhibited minima during austral summer (0.9 ± 0.1‰, n = 5) and maxima during austral fall (1.3 ± 0.3‰, n = 6) and austral spring (1.6 ± 0.1‰, n = 5), indicating a shift from HO x ‐ to O 3 ‐dominated chemistry in the atmospheric derived SO 4 2− component. Overall, the budgets of NO 3 − and SO 4 2− (TSP) at the South Pole were complex functions of transport, localized chemistry, biological activity, and meteorological conditions, and these results will be important for interpretations of oxyanions in ice core records in the interior of Antarctica.

Journal of Geophysical Research D: Atmospheres

Influence of land-atmosphere feedbacks on temperature and precipitation extremes in the GLACE-CMIP5 ensemble

We examine how soil moisture variability and trends affect the simulation of temperature and precipitation extremes in six global climate models using the experimental protocol of the Global Land-Atmosphere Coupling Experiment of the Coupled Model Intercomparison Project, Phase 5 (GLACE-CMIP5). This protocol enables separate examinations of the influences of soil moisture variability and trends on the intensity, frequency, and duration of climate extremes by the end of the 21st century under a business-as-usual (Representative Concentration Pathway 8.5) emission scenario. Removing soil moisture variability significantly reduces temperature extremes over most continental surfaces, while wet precipitation extremes are enhanced in the tropics. Projected drying trends in soil moisture lead to increases in intensity, frequency, and duration of temperature extremes by the end of the 21st century. Wet precipitation extremes are decreased in the tropics with soil moisture trends in the simulations, while dry extremes are enhanced in some regions, in particular the Mediterranean and Australia. However, the ensemble results mask considerable differences in the soil moisture trends simulated by the six climate models. We find that the large differences between the models in soil moisture trends, which are related to an unknown combination of differences in atmospheric forcing (precipitation, net radiation), flux partitioning at the land surface, and how soil moisture is parameterized, imply considerable uncertainty in future changes in climate extremes.

Journal of Geophysical Research D: Atmospheres

Testing the accuracy of a 1-D volcanic plume model in estimating mass eruption rate

During volcanic eruptions, empirical relationships are used to estimate mass eruption rate from plume height. Although simple, such relationships can be inaccurate and can underestimate rates in windy conditions. One-dimensional plume models can incorporate atmospheric conditions and give potentially more accurate estimates. Here I present a 1-D model for plumes in crosswind and simulate 25 historical eruptions where plume height H obs was well observed and mass eruption rate M obs could be calculated from mapped deposit mass and observed duration. The simulations considered wind, temperature, and phase changes of water. Atmospheric conditions were obtained from the National Center for Atmospheric Research Reanalysis 2.5° model. Simulations calculate the minimum, maximum, and average values ( M min , M max , and M avg ) that fit the plume height. Eruption rates were also estimated from the empirical formula M empir = 140 H obs 4.14 ( M empir is in kilogram per second, H obs is in kilometer). For these eruptions, the standard error of the residual in log space is about 0.53 for M avg and 0.50 for M empir . Thus, for this data set, the model is slightly less accurate at predicting M obs than the empirical curve. The inability of this model to improve eruption rate estimates may lie in the limited accuracy of even well-observed plume heights, inaccurate model formulation, or the fact that most eruptions examined were not highly influenced by wind. For the low, wind-blown plume of 14–18 April 2010 at Eyjafjallajökull, where an accurate plume height time series is available, modeled rates do agree better with M obs than M empir .

Journal of Geophysical Research D: Atmospheres