Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “rock deformation”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

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

At least 163 records · Page 9

A Model of the Chicxulub Impact Basin Based on Evaluation of Geophysical Data, Well Logs, and Drill Core Samples

Abundant evidence now shows that the buried Chicxulub structure in northern Yucatan, Mexico, is indeed the intensely sought-after source of the ejecta found world-wide at the Cretaceous-Tertiary (K/T) boundary. In addition to large-scale concentric patterns in gravity and magnetic data over the structure, recent analyses of drill-core samples reveal a lithological assemblage similar to that observed at other terrestrial craters. This assemblage comprises suevite breccias, ejecta deposit breccias (Bunte Breccia equivalents), fine-grained impact melt rocks, and melt-matrix breccias. All these impact-produced lithologies contain diagnostic evidence of shock metamorphism, including planar deformation features in quartz, feldspar, and zircons; diaplectic glasses of quartz and feldspar; and fused mineral melts and whole-rock melts. In addition, elevated concentrations of Ir, Re, and Os, in meteoritic relative proportions, have been detected in some melt-rock samples from the center of the structure. Isotopic analyses, magnetization of melt-rock samples, and local stratigraphic constraints identify this crater as the source of K/T boundary deposits.

Sharpton, Virgil L.↗

Gravity-driven deformation of the crust on Venus

Intense horizontal deformation features are observed on the surface of Venus despite the apparent absence of terrestrial-style plate tectonics. High surface temperature, low erosion rate, and regions of elevated terrain suggest that gravity spreading may be responsible for some of these features. Gravity spreading on the Earth is known to cause extensive deformation in certain thin-skinned tectonic regions, where the uppermost layers become detached along a decollement. While the high pore pressures and soft sedimentary rocks frequently observed along decollements on the Earth are not plausible for Venus, temperature induced regions of ductile deformation may exist. Simple viscous models of gravity spreading above a ductile region of the crust are developed to examine the effects of rheology, thermal gradient, surface slope and crustal thickness on deformation. For a reasonable range of parameters, the models show that surface-slope-driven ductile deformation at relatively shallow depths in the crust is possible and may be capable of detaching the uppermost rigid portion of the crust, resulting in horizontal translation.

Smrekar, Suzanne↗

Impactite and pseudotachylite from Roter Kamm Crater, Namibia

Pseudotachylite is known to occur in a variety of geologic settings including thrust belts (e.g., the Alps and the Himalayas) and impact craters such as Roter Kamm, Namibia. Controversy exists, however, as to whether pseudotachylite can be produced by shock brecciation as well as by tectonic frictional melting. Also open to debate is the question of whether pseudotachylites form by frictional fusion or by cataclasis. It was speculated that the pseudotachylite at Roter Kamm was formed by extensional settling and adjustment of basement blocks during 'late modification stage' of impact. The occurrence of pseudotachylite in association with rocks resembling quenched glass bombs and melt breccias in a relatively young crater of known impact origin offers a rare opportunity to compare features of these materials. Petrographic, x-ray diffraction, and electron microprobe analyses of the impactites and pseudotachylites are being employed to determine the modes of deformation and to assess the role of frictional melting and comminution of adjacent target rocks.

Degenhardt, J. J., Jr.↗

Structure of Southwestern Michipicoten Greenstone Belt, Ontario: Evidence for Archean accretion

The results of detailed mapping of a well exposed (by a large fume kill) Archean Michipicoten greenstone belt of Ontario is discussed. Numerous structural features are described including soft sediment deformation, thrust faults, isoclinal folds, and sill intrusion, which apparently formed prior to the earliest flattening cleavage. It was inferred that all of these events may have taken place before the rocks were completely lithified, and the question was raised as to whether early soft sediment and tectonic deformation were coeval. If so, the overall picture would permit (but not compel) recognition of a convergent oceanic environment such as an accretionary wedge or forearc basin.

Mcgill, George E.↗

A preliminary assessment of the fluid erosional processes that shaped the Martian outflow channels

Catastrophic floods may be the mechanism which created Martian outflow channels and the distinctive assemblage of channel landforms, including regional and local anastomosing patterns, expanding and contracting reaches correlated with flow constrictions, streamlined uplands, inner channels with recessional headcuts, pendant forms on the downcurrent sides of flow obstacles, and scour marks around the obstacles. A terrestrial landscape which contains analogs to all the Martian outflow features is the Channeled Scabland of eastern Washington, formed by catastrophic flooding across differentially resistant basalt, loess and sedimentary rock. Both hydrodynamic theory and the properties of the boundary undergoing deformation need to be considered in studying the forms that result from catastrophic flooding.

Baker, V. R.↗

Modeling of fluid injection-induced fault reactivation using coupled fluid flow and mechanical interface model

The present study is aimed at developing a numerical model to reproduce coupled hydro-mechanical processes associated with fault reactivation by fluid injection in low permeability rock, as part of the DECOVALEX-2019 project Task B. We proposed a modeling approach for simulating the processes using the TOUGH-FLAC simulator, and modeled a fault reactivation experiment conducted at Mont Terri Rock Laboratory in Switzerland. The first step of the study involved benchmark calculations considering a simplified fault plane and geometry. Fluid flow along a fault was modeled using elements of aperture-sized thickness on the basis of Darcy's law and the cubic law in TOUGH2, whereas the mechanical behavior of a single fault was represented by zero-thickness interface elements in FLAC3D upon which a slip and/or separation is allowed. A methodology to connect a TOUGH2 volume element to a FLAC3D interface element was developed for handling the hydro-mechanical interactions on the fault during fluid injection. Two different fault models for describing the evolutions of hydraulic aperture by elastic fracture opening and failure-induced aperture increase were considered in the benchmark calculations. In the coupling process, the changes in geometrical features and hydrological properties induced by mechanical deformation were continuously updated. The transient responses of the fault and host rock to stepwise pressurization were examined during the simulation. The hydro-mechanical behavior, including the injection flow rate, pressure distribution around the borehole, stress conditions, and displacements in normal and shear directions were monitored in the surrounding rock and along the fault. The results of benchmark calculations suggest that the developed model reasonably represents the hydro-mechanical behavior of a fault and the surrounding rock. This modeling approach was applied to the fault reactivation experiment of the Mont Terri Rock Laboratory. In this interpretive modeling, a parametric study was conducted to examine the effects of input parameters regarding in situ stress and fault properties on the hydro-mechanical responses of the fault to water injection. Then, an optimal parameter set to reproduce the field experiment results was chosen by trial-and-error. The injection flow rate and pressure response during fault reactivation closely matched those obtained at the site, which indicates the capability of the model to appropriately capture the progressive pathway evolution during fault reactivation tests at the site. Finally, the anchor displacements were overestimated by the model, but a fair agreement was obtained in terms of the order of magnitude and the variation tendency.

58 GEOSCIENCES↗

Microscale mechanical modeling of deformable geomaterials with dynamic contacts based on the numerical manifold method

Abstract Micromechanical modeling of geomaterials is challenging because of the complex geometry of discontinuities and potentially large number of deformable material bodies that contact each other dynamically. In this study, we have developed a numerical approach for micromechanical analysis of deformable geomaterials with dynamic contacts. In our approach, we detect contacts among multiple blocks with arbitrary shapes, enforce different contact constraints for three different contact states of separated, bonded, and sliding, and iterate within each time step to ensure convergence of contact states. With these features, we are able to simulate the dynamic contact evolution at the microscale for realistic geomaterials having arbitrary shapes of grains and interfaces. We demonstrate the capability with several examples, including a rough fracture with different geometric surface asperity characteristics, settling of clay aggregates, compaction of a loosely packed sand, and failure of an intact marble sample. With our model, we are able to accurately analyze (1) large displacements and/or deformation, (2) the process of high stress accumulated at contact areas, (3) the failure of a mineral cemented rock samples under high stress, and (4) post-failure fragmentation. The analysis highlights the importance of accurately capturing (1) the sequential evolution of geomaterials responding to stress as motion, deformation, and high stress; (2) large geometric features outside the norms (such as large asperities and sharp corners) as such features can dominate the micromechanical behavior; and (3) different mechanical behavior between loosely packed and tightly packed granular systems.

58 GEOSCIENCES↗

Deciphering Extreme Mineral Records; Microstructural Phase Heritage of Shocked Materials

High-pressure minerals such as coesite, stishovite, and reidite are often used as indicators of shock metamorphism. However, high-pressure and high-temperature phases are metastable at ambient conditions and, therefore, often revert to a more energetically favorable phase. During solid-state transformations, the crystallographic orientations of the stable polymorph are controlled by the transformation pathway from the parent. Microstructural orientation analysis using electron backscatter diffraction (EBSD) of the stable phase can reveal systematic intercrystalline orientation relationships (OR) diagnostic of solid-state transformations from these high-pressure or high-temperature polymorph phases (Cayron et al., 2006). Knowledge of polymorph stability fields for a given system can thus be used to infer minimum shock P and/or T conditions that rocks have experienced, which are otherwise unavailable using traditional thermobarometers such as element partition coefficients. Using orientation relationships to infer the former presence of mineral phases has been termed microstructural “phase heritage”, and has successfully been used to probe for evidence of extreme pressures and temperatures resulting from hypervelocity impacts. The reconstruction of polymorphs of ZrSiO4 (zircon, reidite) and ZrO2 (baddeleyite) have further elucidated the behavior of these phases under hypervelocity impact conditions (Timms et al., 2017a; White et al., 2018), and OR analysis of shocked monazite has identified a previously unrecognized tetragonal high-P polymorph (Erickson et al., 2019). Zircon will convert to the low tetragonal polymorph reidite at shock stresses above ~21 GPa. However, upon decompression at temperatures above 1200 °C reidite will revert to zircon (Kusaba et al., 1985). Recrystallized zircon neoblasts encased in impact melt often show systematic misorientation relationships of 90° about <110> indicative of reversion from the reidite (Fig. 1), and preserving cryptic evidence of the high-P history that is often erased by the post shock thermal spike (Cavosie et al., 2016). Microstructural OR analysis of a dissociated zircon corona composed of baddeleyite (monoclinic-ZrO2) has identified the former presence of cubic zirconia within impact melt glass from the Mistastin Lake impact structure, Labrador, CA. Based on stability fields in the ZrO2-SiO2 binary phase diagram, cubic zirconia records thermal conditions in excess of 2370 °C at ambient P, indicating that the Mistastin impact melt is the hottest rock identified on Earth’s surface (Timms et al., 2017b). Shock-deformed monazite (monoclinic La,Ce,ThPO4) containing lamellae comprised of interlocking laths has been identified from the Haughton Dome impact structure, Nunavut, Canada and Nördlinger-Ries Crater, southern Germany. The lath-structured lamellae are composed of four systematic orientation variants, which OR analyses indicate originate from a tetragonal parent whereby the [010]monoclinic aligns with either [100] or [010]tetragonal (Erickson et al., 2019). These results highlight utility of microstructural phase heritage analyses for the identification of unstable high-P or high-T polymorphs that uniquely record the extremely transient shock conditions produced by hypervelocity impacts.

Timmons Erickson↗

Topographic Response to the Yakutat Block Collision

The principal objective of this grant and this research were to investigate the topographic development of an active glaciated orogenic belt in southern Alaska as that development relates to patterns of erosion and crustal deformation. A specific objective of the research was to investigate feedbacks between mountain building, orographic affects on climate, and patterns of exhumation and rock uplift. To that end, an orogen-scale analysis of topography was conducted with the aid of digital elevation models, magnitudes and patterns of crustal deformation were compiled from existing literature, present and past climate patterns were constrained using the modern and past distribution of glaciers, and styles, magnitudes, and extent of erosion were constrained with observations from the 1998 field season.

Stock, Joann M.↗

Crystallographic Texture, Structure, and Stress Transmission in Nugget Sandstone Examined With X‐Ray Tomography and Diffraction Microscopy

Subsurface processes in sandstones are controlled by porosity, permeability, and deformation mechanisms, all of which are controlled by a complex interplay of crystallographic rock texture, structure, and micromechanics. Texture, structure, and micromechanics have historically been studied using optical and electron microscopy of thin-sections. Here, we employed a new combination of in situ X-ray tomography and ray diffraction microscopy to study crystallographic texture, structure, and grain stresses in 3D. We examined these features in a sample of Nugget sandstone, a sandstone constituting hydrocarbon reservoirs across the American West. Our aims are threefold. First, we demonstrate the utility of X-ray diffraction microscopy probes for revealing texture, structure, and stress transmission in 3D. Second, we apply these techniques to Nugget sandstone and discuss findings in the context of prior work. Third, we study grain stress tensor evolution during mechanical compression to examine whether their heterogeneity and orientation evolution reflect that of inter-particle forces in granular materials. Our results show: (a) larger grains featured higher intra-granular misorientations, possibly from an increased prevalence of cements; (b) pores closed parallel to the loading direction and opened normal to loading; (c) grain stresses featured heterogeneity and orientations similar to inter-particle forces in non-cohesive granular materials; (d) grains featured compressive stresses in the loading direction and tensile stresses orthogonal to the loading direction, the latter resisting sample dilation and grain separation. Our work demonstrates the first known application of multi-modal X-ray tomography and diffraction microscopy to sandstone, providing new 3D insight into the nature of quartz cement and stress evolution.

high energy diffraction microscopy↗

A reconnaissance space sensing investigation of the crustal structure for a strip from the eastern Sierra Nevada to the Colorado Plateau: April 1971

The author has identified the following significant results. An area of anomalous linear topographic grain and color expressions was recognized in Apollo 9 and ERTS-1 imagery along the Colorado River of northwestern Arizona and southern Nevada. Field reconnaissance and analysis of U-2 photography has shown the anomaly to be a zone of north to north-northwest trending dike swarms and associated granitic plutons. The dikes vary in composition from rhyolite to diabase, with an average composition nearer rhyolite. Shearing and displacement of host rocks along dikes suggest dike emplacement along active fault zones. Post-dike deformation has resulted in shearing and complex normal faulting along a similar north-south trend. The epizonal plutonism and volcanism of this north-south belt appears to represent a structurally controlled volcanogenic province which ends abruptly in the vicinity of Lake Mead at a probable eastern extension of the Las Vegas Shear Zone. The magnitude and chronology of extensional faulting and plutonism recognized in the north-south zone, support the hypothesis that the Las Vegas Shear Zone is a transform fault separating two areas of crustal spreading.

Bechtold, I. C.↗

Shocked quartz and more: Impact signatures in K-T boundary clays and claystones

Quartz grains displaying multiple sets of planar features are described from numerous Cretaceous-Tertiary (K-T) boundary clays and claystones at both marine and nonmarine depositional sites around the world. All these sites also show anomalously high amounts of iridium and enrichments of other siderophile elements in cosmic ratios within these boundary units. This combination of mineralogical and geochemical features are used in support of an impact hypothesis for the end-Cretaceous event. Recently, it was suggested that some combination of explosive and nonexplosive volcanism associated with the formation of the Deccan traps in India could be responsible for the mineralogy and geochemistry seen in the K-T boundary units. Besides the obvious contradition of simultaneous explosive and nonexplosive volcanism from one locality during an instant of geologic time, there remains the difficulty of spreading both iridium (and trace elements in cosmic proportions) and quartz grains around the world by volcanic (atmospheric) transport. In addition, the ability of volcanism to produce the type of shock metamorphism seen in minerals at the K-T boundary was not demonstrated. Multiple sets of shock lamellae in quartz are considered characteristic of shock metamorphism in rocks at the sites of known impact craters and are the type of deformation seen in quartz from K-T boundary clays and claystones. Single sets of poorly defined lamellae described from rare quartz grains in certain volcanic deposits are characteristic of tectonic deformation and do not correspond to the shock lamellae in quartz from K-T sediments and impact structures. So-called shock mosaicism in quartz and feldspar grains described from volcanic deposits can result from many processes other than shock metamorphism, and therefore is not considered to be an effect characteristic solely of shock. The mineralogy of shock-metamorphosed grains at the K-T boundary also argues against a volcanic origin.

Bohor, Bruce F.↗

Lithospheric-scale buckling and thrust structures on Mars: The Coprates rise and south Tharsis ridge belt

Detailed photogeologic mapping documents for the first time the deformation sequence of the Coprates rise. Materials of Middle Noachian through perhaps Early Hesperian age were involved in the horizontal shortening of crustal and lithospheric rocks that produced uplift and asymmetric, east-vergent folding of these layers during the Late Noachian and perhaps Early Hesperian. Previous suggestions of normal faulting as the deforming mechanism are not supported by the new mapping or by dislocation models of rise topography. The most likely mechanism for producing the Coprates rise is folding due to buckling instability and perhaps thrust faulting. An extensive ridge belt concentric to south Tharsis is defined by a newly recognized system of approximately coeval long-wavelength topogtaphic ridges and swales that include the Coprates rise. Horizontal shortening of crust and lithosphere is manifested by periodically spaced structures at two different scales: wrinkle ridges spaced tens of kilometers apart and the south Tharsis system of ridges spaced hundreds of kilometers apart. Both types of structures may be associated with early volcanotectonic activity and crustal thickening in south Tharsis. These and other large ridges indicate that lithospheric buckling and thrust faulting may have been common on early Mars.

Schultz, Richard A.↗

Viscous Behavior of Clay‐Rich Rocks and Its Role in Focused Fluid Flow

Abstract Focused fluid flow is common in sedimentary basins worldwide, where flow structures often penetrate through sandy reservoir rocks, and clay‐rich caprocks. To better understand the mechanisms forming such structures, the impacts of the viscoelastic deformation and strongly nonlinear porosity‐dependent permeability of clay‐rich materials are assessed from an experimental and numerical modeling perspective. The experimental methods to measure the poroviscoelastic and transport properties of intact and remolded shale have been developed, and the experimental data is used to constrain the numerical simulations. It is demonstrated that viscoelastic deformation combined with nonlinear porosity‐dependent permeability triggers the development of localized flow channels, often imaged as seismic chimneys. The permeability inside a channel increases by several orders of magnitude compared to the background values. In addition, the propagation time scale and the channel size strongly depend on the material properties of the fluid and the rock. The time‐dependent behavior of the clay‐rich rock may play a key role in the long‐term integrity of the subsurface formations.

58 GEOSCIENCES↗

A scale-invariant bonded particle model for simulating large deformation and failure of continua

In this paper, a scale-invariant bond model for simulating two-dimensional large deformation of continua is presented. An accurate torque form caused by relative rotation between bonded particles is derived based on beam theory. Validations of the proposed bond model are carried out by large deformation simulations of a simply supported beam subjected to central loading and a cantilever beam with end loading. Comparison of simulations with analytical solutions show a very good agreement. Later, Brazilian tests are simulated to demonstrate the ability of the framework to model failure and fracture process of rock.

58 GEOSCIENCES↗

Global Assessment of Volcanic Debris Hazards from Space

Hazard (slope stability) assessment for different sectors of volcano edifices was successfully obtained from volcanoes in North and South America. The assessment entailed Hyperion images to locate portions of the volcano that were hydrothermally altered to clay rich rocks with zones that were also rich in alunite and other minerals. The identified altered rock zones were field checked and sampled. The rock strength of these zones was calculated from the field and laboratory measurements. Volcano modeling utilizing the distinct element method and limit equilibrium technique, with the calculated strength data was used to assess stability and deformation of the edifice. Modeling results give indications of possible failure volumes, velocities and direction. The models show the crucial role hydrothermally weak rock plays in reducing the strength o the volcano edifice and the rapid identification of weak rock through remote sensing techniques. Volcanoes were assessed in the Cascade Range (USA), Mexico, and Chile (ongoing).

Watters, Robert J.↗

New Constraints on the Slate Islands Impact Structure, Ontario, Canada

The Slate Islands in northern Lake Superior represent the eroded remains of a complex impact crater, originally approximately 32 km in diameter. New field studies there reveal allogenic crater fill deposits along the eastern and northern portions of the islands indicating that this 500-800 Ma impact structure is not as heavily eroded as previously thought. Near the crater center, on the western side or Patterson Island, massive blocks of target rocks, enclosed within a matrix of fine-grained polymict breccia, record the extensive deformation associated with the central uplift. Shatter cones are a common structural feature on the islands and range from less than 3 cm to over 10 m in length. Although shatter cones are powerful tools for recognizing and analyzing eroded impact craters, their origin remains poorly constrained.

Sharpton, Virgil L.↗

Diaplectic Glass Content in Experimentally Shock-loaded Quartz Determined by X-Ray Powder Diffraction

Quartz is the most common mineral of terrestrial crustal rocks and thus a widespread indicator for impact cratering and associated shock metamorphism. Planar deformation features (PDFs) are among the most prominent and diagnostic shock features in quartz and they represent thin lamellae of glass that formed via solid-state transformations. This socalled 'diaplectic' glass becomes pervasive at higher pressures and results in optically isotropic and X-ray-amorphous phases that resemble texturally the original quartz grains (without evidence of melt flow). In the past, it has been shown that the amount of this amorphous material in experimentally shock-loaded quartz correlates with peak shock pressure. Both reports derive the glass content from density measurements of individual crystals employing the equation X(%) = (rho(sub x) - rho(sub 0))/(rho(sub x) - rho(sub gl)), where x and 0 stands for X-ray and average (optical) density, respectively. The density of glass, rho(sub gl), was adopted as 2.2 g/cu cm. Though the same procedures had been applied, the resulting glass content differs significantly among the above studies. In the present study, we are using a new approach based solely on the integral intensity of a single, carefully selected reflection in the XRD pattern, and we will compare our data to those reported in the literature.

Skala, R.↗