SEARCH · Geology Reports
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.
Fractal geometry is used to characterize the roughness of cracked concrete surfaces through a specially built profilometer, and the fractal dimension is subsequently correlated to the fracture toughness and direction of crack propagation. Preliminary results indicate that the fracture surface is indeed fractal over two orders of magnitudes with a dimension of approximately 1.20.
Geologists and petroleum engineers have struggled to identify the mechanisms that drive productivity in horizontal hydraulically fractured oil wells. The machine learning algorithms of Random Forest (RF), gradient boosting trees (GBT) and extreme gradient boosting (XGBoost) were applied to a dataset containing 7311 horizontal hydraulically fractured wells drilled into the middle member of the Bakken Formation from 2010 through 2017. The initial goal is to use these data‐driven machine learning algorithms to identify the most important explanatory predictors of well productivity within nine subareas and the composite area. Predictor variables representing initial gas production, the initial 180‐day water cut, and vertical depth vary spatially and are identified with geologically favorable areas. Well‐completion predictors include the well lateral length, number of fracture stages, volume of proppant per stage, and the volume of injected fluids per stage. The performance of methods is compared based on a common test sample. The analysis then examines the comparative predictive performance of the three algorithms for 1330 wells that had initiated production after the initial 7311 well sample had been producing. The computations of predictor importance identified the initial 180‐day water cut and the 30‐day initial gas production predictors as having a dominant influence in most subareas and for the composite area. The relative importance of well completion predictor variables, that is, the number of fracture stages per well, volume of injected proppant per stage, volume of injected fluids per stage, and lateral length, varied considerably across the subareas. For the common test or holdout sample, the models calibrated with the XGBoost algorithm had superior predictive power. The predictive power of all the algorithms trained on the data from the original sample suffered some loss when tested with a sample of wells that had started production after the end of that period. Implications of the empirical findings and strategies to mitigate loss of predictive power are discussed in the concluding section.
Large-scale gravitational spreading and movement along fractures of steep-sided ridges in the mountainous areas of the western United States are characterized by linear fissures, trenches, and uphill-facing scarps on tops and sides of ridges. Spreading appears to take place by movement along disconnected planes and/or by slow plastic deformation of a rock mass. In some places, valleyward squeezing out of weak shales overlain by rigid rocks causes extensional fracturing and outward movement of the rigid layers, as illustrated by extension of two laccoliths overlying Mancos Shale, one at Dolores Peak and another at Crested Butte in western Colorado. Gravitational forces acting on a ridge of more homogeneous material causes tensional spreading of the ridge parallel to its long axis, for example in fractured granitic rock north of Mt. Massive in central Colorado, where a survey course has been established to monitor the movement. Recognition and understanding of these large-scale gravitational features and the mechanism that causes them are pertinent to site selection and design of engineering structures in high mountains. If fractures extend to considerable depth and if movement is continuing, engineering structures in valleys or tunnels through the spreading ridges could be damaged.
In recent years, geologists and petroleum engineers have struggled to clearly identify the mechanisms that drive productivity in horizontal, hydraulically-fractured oil wells producing from the middle member of the Bakken formation. This paper fills a gap in the literature by showing how this play’s heterogeneity affects factors that drive well productivity. It is important because understanding the relative strength of productivity drivers and how predictors vary spatially facilitates best-practices for well site selection and well completion design. The paper describes an application of the Random Forest (RF) machine learning technique to identify these mechanisms and to evaluate their importance across 9 subareas of the North Dakota portion of the Bakken play. The study examined productivity of 7311 wells initiating production from 2010 through 2017. Well productivity varied considerably across the 9 subareas within the play, so it was not surprising that the dominant predictors, the initial 180-day water cut and the 30-day initial gas production, vary spatially to mirror local conditions that strongly affect well productivity. The relative importance of well completion predictor variables, that is, the numbers of fractures stages per well, volume of injected proppant per stage, volume of injected fluids per stage, and lateral length, varied considerably across the subareas. Statistical permutation tests are presented that generally confirm the importance rankings. Subarea Random Forest models explained from 50 percent to 82 percent of the variation in productivity test samples while the play-wide model explained 73 percent of the test sample well productivity. Weakness in the predictive ability of the Random Forest models are traced to the limited variability in the training data. Implications of the empirical findings regarding the Bakken play for operators and for research and government institutions are discussed in the concluding section.
Rock environments both underground and on Earth’s surface show indications of energetic macroscale fracture. In tunnels and excavations, these manifest as rockbursts—energetic explosions of rock that can damage engineering projects, and may pose ongoing financial and safety risk as rock stresses adjust during post-failure relaxation. In natural settings at the surface, evidence for rockbursts exist in the form of tent-like structures of ruptured exfoliation sheets, but few direct observations of such events exist, precluding the analysis of how natural rock formations may evolve after rupture. Here we investigate the post-failure evolution of a granitic rock dome following rapid fracture events (i.e., surficial rockbursts) that occurred in California, USA during 2014–2016. Building upon previous work that showed a thermal stress origin for the observed fracturing, we investigate the return to background stress conditions (i.e., stress relaxation) observed in both short- (week, month) and long-term (multi-year) rock deformation trends. Acoustic emissions, deformation, and environmental monitoring data indicate that partially detached rock sheets forming the surface of the dome undergo fracture aperture closing during cooling periods, concurrent with reduction of rock stress by the source of forcing (i.e., thermal stress). However, with sufficient critical and/or subcritical fracture, our observations also show that rock sheets can become decoupled from the source of stress, resulting in a long-term return to background stress conditions. Our results provide insight into the cyclic and likely ephemeral nature of rock fracture in surficial rock domes, as well as in underground rockburst environments.
Landslides can cause seismic disturbances; landslides can also result from seismic disturbances, and earthquake-induced slides have caused loss of life in many countries. Slides can cause disastrous flooding, particularly when landslide dams across streams are breached, and flooding may trigger slides. Slope movement in general is a major process of the geologic environment that places constraints on engineering development. In order to understand and foresee both the causes and effects of slope movement, studies must be made on a regional scale, at individual sites, and in the laboratory. Areal studies — some embracing entire countries — have shown that certain geologic conditions on slopes facilitate landsliding; these conditions include intensely sheared rocks; poorly consolidated, fine-grained clastic rocks; hard fractured rocks underlain by less resistant rocks; or loose accumulations of fine-grained surface debris. Field investigations as well as mathematical- and physical-model studies are increasing our understanding of the mechanism of slope movement in fractured rock, and assist in arriving at practical solutions to landslide problems related to all kinds of land development for human use. Progressive failure of slopes has been studied in both soil and rock mechanics. New procedures have been developed to evaluate earthquake response of embankments and slopes. The finite element method of analysis is being extensively used in the calculation of slope stability in rock broken by joints, faults, and other discontinuities.
This paper examines the environmental processes and mechanisms that govern differential heaving in steeply dipping claystone bedrock near Denver, Colorado. Three potential heave mechanisms and causal processes were evaluated: (1) rebound expansion, from reduced overburden stress; (2) expansive gypsum-crystal precipitation, from oxidation of pyrite; and (3) swelling of clay minerals, from increased ground moisture. First, we documented the effect of short-term changes in overburden stress, atmospheric exposure, and ground moisture on bedrock at various field sites and in laboratory samples. Second, we documented differential heaving episodes in outcrops and at construction and developed sites. We found that unloading and exposure of the bedrock in construction-cut areas are essentially one-time processes that result in drying and desiccation of the near-surface bedrock, with no visible heaving response. In contrast, wetting produces a distinct swelling response in the claystone strata, and it may occur repeatedly as natural precipitation or from lawn irrigation. We documented 2.5 to 7.5 cm (1 to 3 in.) of differential heaving in 24 hours triggered by sudden infiltration of water at the exposed ground surface in outcrops and at construction sites. From these results, we interpret that rebound and pyrite weathering, both of which figure strongly into the long-term geologic evolution of the geologic framework, do not appear to be major heave mechanisms at these excavation depths. Heaving of the claystone takes two forms: (1) hydration swelling of dipping bentonitic beds or zones, and (2) hydration swelling within bedrock blocks accommodated by lateral, thrust-shear movements, along pre-existing bedding and fracture planes.
A geotechnical investigation of the Pierre Shale near Hayes, South Dakota, was conducted by the U.S. Geological Survey as a basis for evaluating problems in deep excavations into that formation. The physical and mechanical properties of the shale were determined through use of core holes drilled to a maximum depth of 184 m. In situ borehole determinations included a gravimeter survey, pressuremeter testing, thermal profile measurements, and borehole velocity measurements. Onsite and offsite laboratory measurements included rebound measurements, sonic velocity measurements of shear and primary waves, X-ray mineralogy and major element determinations, size analyses, fracture analyses, fabric analyses, and determination of thermal properties.Below 15–22 m, the shale is an unweathered, saturated, overconsolidated, underpressured clay shale with a clay-mineral content ranging between 50 and 100%, dominantly composed of mixed-layer illitic smectites. The physical and mechanical properties vary widely. The variation is related to the clay mineral content (especially in bentonite zones), a large transverse mechanical anisotropy, and zones of fractures and microfractures, which may result from rebound caused by erosion. These may contribute to slope instability over large areas. The thermal and mechanical properties change markedly if the shale is permitted to dry out. The state of stress and overconsolidation appear to be functions of the depositional and erosional history of the deposit. Both are markedly affected by the large fracture zones. The properties of the clay shale indicate problems that may be encountered in excavation and use of deep underground facilities. Key words : anisotropy, characterization, clay shale, consolidation state, physical properties, rebound, relaxation, stress state, thermal properties.
Slow gravitational failures of mountain peaks and ridges are poorly understood. Herein, we report on 50 years of studies at a slowly spreading castellate ridge in the Sawatch Range in central Colorado. The orientations of geomorphic-structural features indicate that the fractured Precambrian granitic rock underlying the ridge has extended and spread northwestward toward the formerly glacier-covered Busk Creek valley. Results from surveying, field-based geomorphic-structural mapping using lidar, rock mass quality measurements, a passive seismic survey, and satellite radar provide a major update to research started by U.S. Geological Survey researchers in the 1970s and 1980s. New insights include a recognition that the entire ridge has slowly moved by concurrent sliding along an inferred northwest dipping, compound basal-slip surface (or zone), and through the formation of multiple grabens by normal faulting and flexural toppling along sets of pre-existing fractures that dip moderately (∼45°) to the southeast. We were unable to distinguish the presence of a sudden and strong contrast in seismic velocity across the inferred slip surface. Movement during the 50-year study period has been episodic and gradually decreasing, in correspondence with decreasing cumulative annual precipitation and increasing mean annual air temperatures. The fastest moving area, just upslope from the glacier trimline, had an average horizontal velocity of 3 – 4 mm/yr. Evidence suggests that movement started as a paraglacial response mechanism, but because of the site’s proximity to the Rio Grande Rift, we cannot exclude earthquake shaking as a mechanism for initiation or enhancement of slope movement. An estimate of longer-term horizontal movement from the exposed basal-slip surface at the uphill side of the ridgetop graben is ∼1.1 mm/yr for the 13 – 14 ky post-glacial period. Broad implications of our work are that: (1) long-term measurements (decades or longer) of slope movement can add insights into how sackungen form and evolve through time; (2) the identification of thrust faults and toes in zones of compression near valley bottoms can be crucial for interpreting sackung failure mechanisms, and (3) the use of passive-seismic techniques to identify the depth to a slip surface may not be successful in granitic terrain dominated by planar fractures and subtle changes in rock-mass characteristics.
Tertiary and/or Quaternary tectonic faulting is documented in three areas of southernmost Illinois: the Fluorspar Area Fault Complex (FAFC) in Pope and Massac Counties, the Ste. Genevieve Fault Zone (SGFZ) in Alexander and Union Counties, and the Commerce Fault Zone (CFZ) in Alexander County. In the FAFC, faults that strike NE and NNE displace Mounds Gravel (late Miocene to early Pleistocene) and, locally, the Metropolis terrace gravel (Pleistocene; pre-Woodfordian). No Woodfordian or younger deposits are deformed. Faults typically outline narrow, linear grabens that formed under tension with a component of strike slip. North-south to NW-trending vertical faults near the southeast end of the SGFZ displace Eocene sediments. Again, faults outline narrow grabens and show indications of strike slip. Deformed Quaternary sediments have not been observed. The CFZ, which trends northeast, displaces Mounds Gravel in Illinois and units as young as Peoria Silt (Woodfordian) in Missouri. Quaternary movement has been interpreted as right-lateral strike-slip. The CFZ coincides with a subtle gravity and magnetic lineament and seems to reflect a major feature in the basement. Surface expression in Illinois is subtle, but mafic and ultramafic intrusions, hydrothermal alteration and small faults align with the Commerce geophysical lineament. Earthquake foci in Missouri and Illinois lie on or close to the CFZ; some focal mechanisms fit the fault trend. Among these structures, only the CFZ exhibits slip that conforms to the current stress field (principal compressive stress axis E-W to ENE-WSW). Possibly, the stress field changed during Neogene time. Alternatively, high fluid pressures or local stress concentrations may have induced slip on less favorably oriented fractures. Tighter constraints are needed on timing, magnitude, and direction of Neogene displacement.
Samples of sedimentary material from interbeds between basalt flows and from fractures in the flows, taken from two drill cores at the Radioactive Waste Management Complex at the Idaho National Engineering Laboratory were analyzed for (1) particle-size dribution, (2) bulk mineralogy, (3) clay mineralogy, (4) cation-exchange capacity, and (5) carbonate content. Thin sections of selected sediment material were made for petrographic examination. Preliminary interpretations indicate that (1) it may be possible to distinguish the various sediment interbeds on the basis of their mineralogy, (2) the presence of carbonate horizons in sedimentary interbeds may be utilized to approximate the time of exposure and the climate while the surface was exposed (which affected the hydrogeologic character of the sediment), and the type and orientation of fracture-filling material may be utilized to determine the mechanism by which fractures were filled. (USGS)
Cycles of stress build-up and release are inherent to tectonically active planets. Such stress oscillations impart strain and damage, prompting mechanically loaded rocks and materials to fail. Here, we investigate, under uniaxial conditions, damage accumulation and weakening caused by time-dependent creep (at 60, 65, and 70% of the rocks’ expected failure stress) and repeating stress oscillations (of ± 2.5, 5.0 or 7.5% of the creep load), simulating earthquakes at a shaking frequency of ~ 1.3 Hz in volcanic rocks. The results show that stress oscillations impart more damage than constant loads, occasionally prompting sample failure. The magnitudes of the creep stresses and stress oscillations correlate with the mechanical responses of our porphyritic andesites, implicating progressive microcracking as the cause of permanent inelastic strain. Microstructural investigation reveals longer fractures and higher fracture density in the post-experimental rock. We deconvolve the inelastic strain signal caused by creep deformation to quantify the amount of damage imparted by each individual oscillation event, showing that the magnitude of strain is generally largest with the first few oscillations; in instances where pre-existing damage and/or the oscillations’ amplitude favour the coalescence of micro-cracks towards system scale failure, the strain signal recorded shows a sharp increase as the number of oscillations increases, regardless of the creep condition. We conclude that repetitive stress oscillations during earthquakes can amplify the amount of damage in otherwise mechanically loaded materials, thus accentuating their weakening, a process that may affect natural or engineered structures. We specifically discuss volcanic scenarios without wholesale failure, where stress oscillations may generate damage, which could, for example, alter pore fluid pathways, modify stress distribution and affect future vulnerability to rupture and associated hazards.
The typical initial piezometric response to longwall undermining is a decline in head due to the opening of fractures and bedding planes during early subsidence. However, in studies over two active longwall mines in southern Illinois, temporary rises in head just before subsidence were observed in piezometers constructed in low-permeability units. Although the initial phase of subsidence is considered dilational, these head increases indicate compressional effects that raise pore-water pressures. Possible mechanisms are shear stresses at the leading edge of subsidence or transmission of stress related to dewatering of underlying permeable units.
Landsliding is a common problem in eastern Fairfax County, an area underlain by Potomac Group (Lower Cretaceous) Coastal Plain deposits of silt and clay interbedded and interfingered with sand and gravel. The slides commonly are present in clay and silt that, on the basis of laboratory tests, appear to be much too strong to have failed. However, the very plastic silt and clay deposits are commonly cut by long continuous to short discontinuous high-angle and subhorizontal joints, shears, and faults. These fractures can contribute to sliding. Failure along faults and shears takes place because the relative movement has greatly weakened and softened the clay and silt. The mechanism for failure along the joints is less obvious but may be related to infilling along joints, slight movement due to erosional unloading, swelling of clays, and softening along joints. Field investigations suggest that the three best developed high-angle joint sets in eastern Fairfax County have trends in the northeast quadrant. Lineaments plotted on a lineament map prepared from Landsat imagery also trend northeast parallel to the regional joint orientation. This regional orientation suggests that some of the fractures are tectonic in origin. Therefore, they should be anticipated in a broad zone many miles wide in the vicinity of the Fall Line. In addition, fractures of unknown origin contributing to landsliding, have been observed in massive Potomac Group clay and silt deposits. Irrespective of origin, these fractures require careful investigation and engineering to prevent slope failure at critical locations during and after major excavation.
A new high-resolution directional borehole radar-logging tool (DBOR tool) was used to log three wells at the Idaho National Engineering and Environmental Laboratory (INEEL). The radar system uses identical directional cavity-backed monopole transmitting and receiving antennas that can be mechanically rotated while the tool is stationary or moving slowly in a borehole. Faster reconnaissance logging with no antenna rotation was also done to find zones of interest. The microprocessor-controlled motor/encoder in the tool can rotate the antennas azimuthally, to a commanded angle, accurate to a within few degrees. The three logged wells in the unsaturated zone at the INEEL had been cored with good core recovery through most zones. After coring, PVC casing was installed in the wells. The unsaturated zone consists of layered basalt flows that are interbedded with thin layers of coarse-to-fine grained sediments. Several zones were found that show distinctive signatures consistent with fractures in the basalt. These zones may correspond to suspected preferential flow paths. The DBOR data were compared to core, and other borehole log information to help provide better understanding of hydraulic flow and transport in preferential flow paths in the unsaturated zone basalts at the INEEL.
Shiveluch volcano (Kamchatka, Russia) is an active andesitic volcano with a history of explosive activity, dome extrusion, and structural collapse during the Holocene. The most recent major (> 1 km 3 ) dome collapse occurred in November 1964, producing a ~ 1.5 km 3 debris avalanche that traveled over 15 km from the vent and triggered a phreatic explosion followed by a voluminous (~ 0.8 km 3 ) eruption of juvenile pyroclastic material. Seismic records suggest that the collapse was likely triggered by a magnitude 5.1 earthquake associated with the ascent of magma into the edifice. The geomechanical properties of the pre-1964 dome are unknown; accordingly, the mechanics of the collapse are poorly understood. This project employs numerical slope stability modeling using the finite element method to constrain probable ranges of geomechanical properties for the materials involved in the collapse, considering earthquake loading as the most likely triggering mechanism. Model results show good agreement with the 1964 collapse geometry considering Geological Strength Index and horizontal pseudo-static seismic coefficient ranges of 30–60 and 0.05–0.15 g, respectively, representing variably fractured and altered dome rocks under moderate earthquake loading, confirming that ground acceleration alone could have triggered the dome collapse. Deep-seated rotational sliding is the dominant failure mode, but local extension within the dome during failure appears to play an important role in the development of the collapse. The findings of this work allow for better forward modeling of potential future collapses, the results of which can be incorporated into regional hazard and risk assessments.
Much has been published about double-diffusive convection as a mechanism for explaining variations in composition and temperature within all-liquid natural systems. However, relatively little is known about the applicability of this phenomenon within the heterogeneous rocks of currently active geothermal systems where primary porosity may control fluid flow in some places and fractures may control it in others. The main appeal of double-diffusive convection within hydrothermal systems is-that it is a mechanism that may allow efficient transfer of heat mainly by convection, while at the same time maintaining vertical and lateral salinity gradients. The Salton Sea geothermal system exhibits the following reservoir characteristics: (1) decreasing salinity and temperature from bottom to top and center toward the sides, (2) a very high heat flow from the top of the system that seems to require a major component of convective transfer of heat within the chemically stratified main reservoir, and (3) a relatively uniform density of the reservoir fluid throughout the system at all combinations of subsurface temperature, pressure, and salinity. Double-diffusive convection can account for these characteristics very nicely whereas other previously suggested models appear to account either for the thermal structure or for the salinity variations, but not both. Hydrologists, reservoir engineers, and particularly geochemists should consider the possibility and consequences of double-diffusive convection when formulating models of hydrothermal processes, and of the response of reservoirs to testing and production.
Characterization of rock discontinuities and rock bridges is required to define stability conditions of fractured rock masses in both natural and engineered environments. Although remote sensing methods for mapping discontinuities have improved in recent years, remote detection of intact rock bridges on cliff faces remains challenging, with their existence typically confirmed only after failure. In steep exfoliating cliffs, such as El Capitan in Yosemite Valley (California, USA), rockfalls mainly occur along cliff-parallel exfoliation joints, with rock bridges playing a key role in the stability of partially detached exfoliation sheets. We employed infrared thermal imaging (i.e., thermography) as a new means of detecting intact rock bridges prior to failure. An infrared thermal panorama of El Capitan revealed cold thermal signatures for the surfaces of two granitic exfoliation sheets, consistent with the expectation that air circulation cools the back of the partially detached sheets. However, we also noted small areas of warm thermal anomalies on these same sheets, even during periods of nocturnal rock cooling. Rock attachment via rock bridges is the likely cause for the warm anomalies in the thermal data. 2-D model simulations of the thermal behavior of one of the monitored sheets reproduce the observed anomalies and explain the temperature differences detected in the rock bridge area. Based on combined thermal and ground-based lidar imaging, and using geometric and rock fracture mechanics analysis, we are able to quantify the stability of both sheets. Our analysis demonstrates that thermography can remotely detect intact rock bridges and thereby greatly improve rockfall hazard assessment.