USGS Science⌕ Search

SEARCH · USGS Science

Results for “Rock Mechanics and Rock Engineering”

Search indexed USGS publications on groundwater, aquifers, geologic maps, mineral resources and earthquakes. Explore source records by subject and place.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

42 records · Page 3Linked to original sources

Mechanical and hydraulic properties of rocks related to induced seismicity

Witherspoon, P.A. and Gale, J.E., 1977. Mechanical and hydraulic properties of rocks related to induced seismicity. Eng. Geol., 11(1): 23-55. The mechanical and hydraulic properties of fractured rocks are considered with regard to the role they play in induced seismicity. In many cases, the mechanical properties of fractures determine the stability of a rock mass. The problems of sampling and testing these rock discontinuities and interpreting their non-linear behavior are reviewed. Stick slip has been proposed as the failure mechanism in earthquake events. Because of the complex interactions that are inherent in the mechanical behavior of fractured rocks, there seems to be no simple way to combine the deformation characteristics of several sets of fractures when there are significant perturbations of existing conditions. Thus, the more important fractures must be treated as individual components in the rock mass. In considering the hydraulic properties, it has been customary to treat a fracture as a parallel-plate conduit and a number of mathematical models of fracture systems have adopted this approach. Non-steady flow in fractured systems has usually been based on a two-porosity model, which assumes the primary (intergranular) porosity contributes only to storage and the secondary (fracture) porosity contributes only to the overall conductivity. Using such a model, it has been found that the time required to achieve quasi-steady state flow in a fractured reservoir is one or two orders of magnitude greater than it is in a homogeneous system. In essentially all of this work, the assumption has generally been made that the fractures are rigid. However, it is clear from a review of the mechanical and hydraulic properties that not only are fractures easily deformed but they constitute the main flow paths in many rock masses. This means that one must consider the interaction of mechanical and hydraulic effects. A considerable amount of laboratory and field data is now available that clearly demonstrates this stress-flow behavior. Two approaches have been used in attempting to numerically model such behavior: (1) continuum models, and (2) discrete models. The continuum approach only needs information as to average values of fracture spacing and material properties. But because of the inherent complexity of fractured rock masses and the corresponding decrease in symmetry, it is difficult to develop an equivalent continuum that will simulate the behavior of the entire system. The discrete approach, on the other hand, requires details of the fracture geometry and material properties of both fractures and rock matrix. The difficulty in obtaining such information has been considered a serious limitation of discrete models, but improved borehole techniques can enable one to obtain the necessary data, at least in shallow systems. The possibility of extending these methods to deeper fracture systems needs more investigation. Such data must be considered when deciding whether to use a continuum or discrete model to represent the interaction of rock and fluid forces in a fractured rock system, especially with regard to the problem of induced seismicity. When one is attempting to alter the pressure distribution in a fault zone by injection or withdrawal of fluids, the extent to which this can be achieved will be controlled in large measure by the behavior of the fractures that communicate with the borehole. Since this is essentially a point phenomenon, i.e., the changes will propagate from a relatively small region around the borehole, the use of a discrete model would appear to be preferable.

Engineering Geology↗

Observations of debris flows at Chalk Cliffs, Colorado, USA: Part 1, in-situ measurements of flow dynamics, tracer particle movement and video imagery from the summer of 2009

Debris flows initiated by surface-water runoff during short duration, moderate- to high-intensity rainfall are common in steep, rocky, and sparsely vegetated terrain. Yet large uncertainties remain about the potential for a flow to grow through entrainment of loose debris, which make formulation of accurate mechanical models of debris-flow routing difficult. Using a combination of in situ measurements of debris flow dynamics, video imagery, tracer rocks implanted with passive integrated transponders (PIT) and pre- and post-flow 2-cm resolution digital terrain models (terrain data presented in a companion paper by STALEY et alii, 2011), we investigated the entrainment and transport response of debris flows at Chalk Cliffs, CO, USA. Four monitored events during the summer of 2009 all initiated from surface-water runoff, generally less than an hour after the first measurable rain. Despite reach-scale morphology that remained relatively constant, the four flow events displayed a range of responses, from long-runout flows that entrained significant amounts of channel sediment and dammed the main-stem river, to smaller, short-runout flows that were primarily depositional in the upper basin. Tracer-rock travel-distance distributions for these events were bimodal; particles either remained immobile or they travelled the entire length of the catchment. The long-runout, large-entrainment flow differed from the other smaller flows by the following controlling factors: peak 10-minute rain intensity; duration of significant flow in the channel; and to a lesser extent, peak surge depth and velocity. Our growing database of natural debris-flow events can be used to develop linkages between observed debris-flow transport and entrainment responses and the controlling rainstorm characteristics and flow properties.

Colorado↗

Field evidence and indicators of rockfall fragmentation and implications for mobility

Rockfall fragmentation can play an important role in hazard studies and the design of protective measures. However, the current lack of modeling tools that incorporate rock fragmentation mechanics is a limitation to enhancing studies and design. This research investigates the fragmentation patterns of rockfalls and analyzes the resulting distribution of fragment sizes within corresponding rockfall deposits. We focus on small rock fragments, which provide insights into the dynamics of the rockfall event and can be used as input for numerical modeling. We analyzed multiple rockfall events from locations worldwide, each exhibiting different degrees of fragmentation. Using image analysis techniques, we mapped all visible blocks, determined their volumes, and measured the distances they travelled from the initial point of impact. A key finding is the identification of three indicators of fragmentation. First, in cases where fragmentation was largely absent, we observed a trend of increasing block size with distance from the impact point or source area, which aligns with previously published findings. However, for energetic rockfall events characterized by intense fragmentation, we observed that small fragments exhibited longer travel distances compared to larger fragments. This distinction allowed us to differentiate blocks primarily resulting from the disaggregation process from those primarily resulting from dynamic fragmentation, with implications for rockfall mobility. Second, although the size distribution of rockfall deposits exhibits a power-law scaling for volumes larger than a minimum size threshold corresponding to a rollover of the distribution, in some case studies a deviation from power-law scaling is observed, indicating a process of larger block comminution due to fragmentation. Third, we found that rockfalls with fragmentation experience reduced mobility, indicated by higher reach angles, and higher lateral dispersion showing a wider distribution of trajectories. We interpret these findings as being directly related to the energy-consuming nature of fragmentation, which prevents farther deposition of fragmented rock blocks.

Albacete province, Lombardy and Aosta Valley, Yose↗

Hydrologic evaluation methodology for estimating water movement through the unsaturated zone at commercial low-level radioactive waste disposal site

This report identifies key technical issues related to hydrologic assessment of water flow in the unsaturated zone at low-level radioactive waste (LLW) disposal facilities. In addition, a methodology for incorporating these issues in the performance assessment of proposed LLW disposal facilities is identified and evaluated. The issues discussed fall into four areas: Estimating the water balance at a site (i.e., infiltration, runoff, water storage, evapotranspiration, and recharge); Analyzing the hydrologic performance of engineered components of a facility; Evaluating the application of models to the prediction of facility performance; and Estimating the uncertainty in predicted facility performance. An estimate of recharge at a LLW site is important since recharge is a principal factor in controlling the release of contaminants via the groundwater pathway. The most common methods for estimating recharge are discussed in Chapter 2. Many factors affect recharge; the natural recharge at an undisturbed site is not necessarily representative either of the recharge that will occur after the site has been disturbed or of the flow of water into a disposal facility at the site. Factors affecting recharge are discussed in Chapter 2. At many sites engineered components are required for a LLW facility to meet performance requirements. Chapter 3 discusses the use of engineered barriers to control the flow of water in a LLW facility, with a particular emphasis on cover systems. Design options and the potential performance and degradation mechanisms of engineered components are also discussed. Water flow in a LLW disposal facility must be evaluated before construction of the facility. In addition, hydrologic performance must be predicted over a very long time frame. For these reasons, the hydrologic evaluation relies on the use of predictive modeling. In Chapter 4, the evaluation of unsaturated water flow modeling is discussed. A checklist of items is presented to guide the evaluation. Several computer simulation codes that were used in the examples (Chapter 6) are discussed with respect to this checklist. The codes used include HELP, UNSAT-H, and VAM3DCG. To provide a defensible estimate of water flow in a LLW disposal facility, the uncertainty associated with model predictions must be considered. Uncertainty arises because of the highly heterogeneous nature of most subsurface environments and the long time frame required in the analysis. Sources of uncertainty in hydrologic evaluation of the unsaturated zone and several approaches for analysis are discussed in Chapter 5. The methods of analysis discussed include a bounding approach, sensitivity analysis, and Monte Carlo simulation. To illustrate the application of the discussion in Chapters 2 through 5, two examples are presented in Chapter 6. The first example is of a below ground vault located in a humid environment. The second example looks at a shallow land burial facility located in an arid environment. The examples utilize actual site-specific data and realistic facility designs. The two examples illustrate the issues unique to humid and arid sites as well as the issues common to all LLW sites. Strategies for addressing the analytical difficulties arising in any complex hydrologic evaluation of the unsaturated zone are demonstrated. The report concludes with some final observations and recommendations.

Report↗

Shorter contributions to stratigraphy and structural geology, 1979

PART A: A system of anticlines lies along the trend of the sinuous course of the Colorado River for a distance of 97 km in the central Grand Canyon. Similar anticlines occur in some perennially wet side canyons. The anticlines are most abundant and well developed along northeast-trending reaches of the main canyon where it is floored by the Cambrian Muav Limestone. Dips of the folded strata are as great as 60?, and the folding locally extends more than 250 m from the river. Low-angle thrust faults in the limbs of the anticlines parallel the river and have formed in response to folding of the comparatively brittle carbonate strata. High-angle reverse kink bands, along which rocks are displaced up toward the river, also parallel the anticlines and have develop2d in response to the upward bulging of the canyon floor. The river anticlines are an unloading phenomenon. They result from lateral squeezing toward the river of saturated shaly parts of the Muav Limestone and underlying Bright Angel Shale. The driving mechanism for the deformation is a stress gradient that results from a difference in lithostatic load between the heavily loaded rocks under the 650-m-high canyon walls and the unloaded canyon floor. Saturation appears to weaken the shaly rocks sufficiently to allow deformation to take place. River anticlines are not present in the eastern Grand Canyon, where the Cambrian rocks also occur at river level. Their absence is explained by a lack of shaly rocks that could flow when saturated. PART B: The current interest in contemporary tectonic processes in the Eastern United States is turning up abundant evidence of crustal movements in late geologic time. Topographic analysis of the highland areas from the southern Blue Ridge to the Adirondack Mountains indicates that most of the landforms owe their origin to erosion of rocks of different resistance rather than to tectonic processes. Most areas of high relief and high altitude have been formed on resistant rocks. The Cambrian and Ordovician belt, containing mostly shale and carbonate rock, on the other hand, forms an extensive lowland from Alabama to the Canadian border and girdles the Adirondack Mountains. Differences in altitude can be explained by the presence of resistant rocks outside the belt; these resistant rocks form local base levels on the streams that drain the belt. A few areas may have undergone local uplift at a higher rate than areas nearby--for example, the Piedmont region northwest of Chesapeake Bay. Most estimates of erosion rates, based on the load transported by streams and of uplift rates, based on removal during a known period of time, are of the same order of magnitude, averaging almost 4x 10^-2 millimeters per year. Rates of uplift, based on study of tilted Pleistocene beaches and repeated geodetic traverses, are at least an order of magnitude higher for comparable areas. Tectonic uplift of the highlands has been slow and involves mostly warping or tilting on a large scale. Erosion rates keep up with or exceed the rate of uplift and have been sufficient to mask evidence of faulting or other differential movements. The high rates of uplift that are inferred on tilted water planes in the glaciated regions or that are measured by differences in repeated geodetic traverses cannot have been sustained for long periods of time. PART C: The Hanson Creek Formation southwest of Eureka, Nev., in the Bellevue Peak Quadrangle is composed of three lithostratigraphic members: (1) a basal dark-gray dolomite, (2) a middle silty thin- to thick-bedded, locally nodular, dark-gray, light-yellow-mottled limestone topped by light-gray dolomite, and (3) an upper dark-gray dolomite, which is herein named the Combs Canyon Dolomite Member. Detailed geologic mapping and accompanying fossil collecting prove that the same lithostratigraphic and biostratigraphic sequence is present in the Mountain Boy Range and 11 km to the south near Wood Cone Peak. Minor differences in

Professional Paper↗

Lateral and vertical channel movement and potential for bed-material movement on the Madison River downstream from Earthquake Lake, Montana

The 1959 Hebgen Lake earthquake caused a massive landslide (Madison Slide) that dammed the Madison River and formed Earthquake Lake. The U.S. Army Corps of Engineers excavated a spillway through the Madison Slide to permit outflow from Earthquake Lake. In June 1970, high streamflows on the Madison River severely eroded the spillway channel and damaged the roadway embankment along U.S. Highway 287 downstream from the Madison Slide. Investigations undertaken following the 1970 flood events concluded that substantial erosion through and downstream from the spillway could be expected for streamflows greater than 3,500 cubic feet per second (ft 3 /s). Accordingly, the owners of Hebgen Dam, upstream from Earthquake Lake, have tried to manage releases from Hebgen Lake to prevent streamflows from exceeding 3,500 ft 3 /s measured at the U.S. Geological Survey (USGS) gaging station 0638800 Madison River at Kirby Ranch, near Cameron, Montana. Management of flow releases from Hebgen Lake to avoid exceeding the threshold streamflow at USGS gaging station 06038800 is difficult, and has been questioned for two reasons. First, no road damage was reported downstream from the Earthquake Lake outlet in 1993, 1996, and 1997 when streamflows exceeded the 3,500-ft 3 /s threshold. Second, the 3,500-ft 3 /s threshold generally precludes releases of higher flows that could be beneficial to the blue-ribbon trout fishery downstream in the Madison River. In response to concerns about minimizing streamflow downstream from Earthquake Lake and the possible armoring of the spillway, the USGS, in cooperation with the Madison River Fisheries Technical Advisory Committee (MADTAC; Bureau of Land Management; Montana Department of Environmental Quality; Montana Fish, Wildlife and Parks; PPL-Montana; U.S. Department of Agriculture Forest Service - Gallatin National Forest; and U.S. Fish and Wildlife Service), conducted a study to determine movement of the Madison River channel downstream from Earthquake Lake and to investigate the potential for bed material movement along the same reach. The purpose of this report is to present information about the lateral and vertical movement of the Madison River from 1970 to 2006 for a 1-mile reach downstream from Earthquake Lake and for Raynolds Pass Bridge, and to provide an analysis of the potential for bed-material movement so that MADTAC can evaluate the applicability of the previously determined threshold streamflow for initiation of damaging erosion. As part of this study channel cross sections originally surveyed by the USGS in 1971 were resurveyed in 2006. Incremental channel-movement distances were determined by comparing the stream centerlines from 14 aerial photographs taken between 1970 and 2006. Depths of channel incision and aggregation were determined by comparing the 2006 and 1971 cross-section and water-surface data. Particle sizes of bed and bank materials were measured in 2006 and 2008 using the pebble-count method and sieve analyses. A one-dimensional hydraulic-flow model (HEC-RAS) was used to calculate mean boundary-shear stresses for various streamflows; these calculated boundary-shear stresses were compared to calculated critical-shear stresses for the bed materials to determine the potential for bed-material movement. A comparison of lateral channel movement distances with annual peak streamflows shows that streamflows higher than the 3,500-ft 3 /s threshold were followed by lateral channel movement except from 1991 to 1992 and possibly from 1996 to 1997. However, it was not possible to discern whether the channel moved gradually or suddenly, or in response to one peak flow, to several peak flows, or to sustained flows. The channel moved between 2002 and 2005 even when streamflows were less than the threshold streamflow of 3,500 ft 3 /s. Comparisons of cross sections and aerial photographs show that the channel has moved laterally and incised and aggraded to varying degrees. The channel has developed meander bends and has incised as much as 5–12 feet (ft) through the upstream part of the Madison Slide (cross sections 1400–800). Near cross section 800, the stream has eroded into the steep right bank between the stream and the road where fill was mechanically placed after 1970. Channel movement also was noted downstream from the Madison Slide. Near Raynolds Pass Bridge, about 3 miles (mi) downstream from Earthquake Lake, elevations across the channel have changed by -1.4 ft to +1.9 ft, but these changes were local in nature and could represent a few rocks or depressions in the bed. Overall, it does not appear that the materials eroded from the Madison Slide are causing aggradation in the subreach near the Raynolds Pass Bridge. Comparisons of critical shear stresses to mean boundary-shear stresses indicate that the D50 particle sizes (median size) along the right side of the bed between cross sections 400 and 500 and along the right side of the bed between cross sections 1300 and 1400 could move at the threshold streamflow. In contrast, most of the D84 particle sizes at those two locations probably will not move at the threshold streamflow. This lack of movement for the larger particles at the threshold streamflow could lead to further armoring of the bed as the D50 and smaller-sized particles are removed from the bed and transported downstream. The Shields parameter values from 0.04 to 0.08 that were used to calculate critical shear stresses could be conservative for a high-gradient stream such as the Madison. A higher, less conservative, Shields parameter would result in higher critical shear stresses, meaning that higher streamflows would be required to move material than those reported herein. In addition, because materials in the channel thalweg are exposed to higher boundary-shear stresses than the materials along the sides of the channel, larger, more erosion-resistant materials likely exist in the deeper parts of the channel where high-flow depths and velocities prevented sediment sampling. Movement of these materials might require higher critical shear stresses than estimated in this report. Characterization of sediment sizes in the center of the stream and observation of bed-material movement for a range of streamflows could provide information to help refine the Shields parameter and critical-shear stress estimates for bed materials in the Madison River downstream from Earthquake Lake. Furthermore, resurveying cross sections and water-surface elevations more frequently (either annually or after high streamflows) could better define the relation between streamflow and lateral and vertical channel movement.

Montana↗