Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “fracturing”

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 199 records · Page 11
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Experimentally validated multiphysics modeling of fracture induced by thermal shocks in sintered UO 2 pellets

Uranium Dioxide (UO 2 ) fuel powers almost all commercial Nuclear Power Plants (NPPs) worldwide, generating carbon-free energy and contributing to the fight against climate change. UO 2 fuel incurs damage and fractures due to large thermal gradients that develop across the fuel pellet during normal and transient operating conditions. A comprehensive understanding of the underlying mechanisms by which these processes take place is still lacking. A combined experimental and computational approach is utilized here to quantify the behavior of UO 2 fuel fracture induced by thermal shock. Here, this work introduces both (1) an experimental study to understand the fuel fracturing behavior of sintered UO 2 pellets when exposed to thermal shock, and (2) a Multiphysics phase-field fracture model capable of simulating this process. Parametric studies were conducted to evaluate the effects of uncertainties in fracture properties on the fracture behavior of UO 2 due to thermal shocking. A set of energy release rate (or equivalently fracture toughness) and contract area (the part of the fuel pellet in direct contact with the cold bath) were able to capture the overall fracture trends of the corresponding experimental data. Our combined approach presents a new method for accounting for the effects of microstructure and sample size on the energy release rate/fracture toughness. The experimental data were collected from multiple experiments that exposed UO 2 pellets to high-temperature conditions (589–676 °C) followed by a quench in sub-zero water. This work demonstrates that joint experimental and computational efforts are able to advance the understanding of thermal fracture in the primary fuel source for existing and future NPPs.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Multiscale spatial analysis of fracture arrangement and pattern reconstruction using Ripley's K-function

This work presents novel multiscale spatial data analytics using Ripley's K-function, as a measure of spatial interaction, to study one-dimensional arrangement of fractures. Fracture spatial arrangements are classified into clustered, anticlustered, or indistinguishable from random by testing statistical significance of the calculated Ripley's K-function. Characterizations of fracture arrangements are performed as a function of length scale and position. Analysis of the K-function along the study interval identifies where fracture clustering and anticlustering occur. Here, a simulation technique is also introduced here to statistically reconstruct spatial arrangements and to generate fracture realizations that are spatially similar to the fractures observed in the field. With this simulation technique, one can also fill spatial gaps in fracture measurements where data are absent, unreliable, or unused. Synthetic as well as field-measured 1D fracture datasets are used for testing and demonstration. Methods introduced in this work can be readily applied to fracture datasets observed in outcrops, borehole image logs, and cores.

58 GEOSCIENCES↗
DOE OSTI2025

Monitoring Fracture Hydromechanical Evolution in the Lab and Field Using Unsupervised Metric Learning

Fractures evolve in time through thermal‐hydraulic‐mechanical‐chemical (THMC) processes that alter their long‐range hydraulic transport properties and modify subsurface behavior and activities. The location of subsurface fractures makes it necessary to use remote sensing techniques such as passive or active seismic monitoring for fracture characterization. In this paper, we develop a machine learning approach to monitor the evolution of fracture properties using passive seismic sources in a laboratory setting and using active seismic monitoring from the Sanford Underground Research Facility in Lead, South Dakota, at a depth of 1.25 km in amphibolite rock during stimulation of natural fractures as well as during induced fracturing. The unsupervised metric learning technique applies tandem neural networks (twin (Siamese) or triplet) with contrastive loss and adaptive margins to track slowly varying systems for which class or similarity labels are not available. The approach adopts locality‐sensitive hashing to divide time‐ordered contiguous data into an arbitrary number of pseudo‐classes. Contrastive‐loss training with many hash bins generates an evolving latent‐space trajectory. This approach enables unsupervised metric learning for seismic data stacks under the condition of contiguous state sampling and slowly varying fracture properties. The displacement discontinuity theory provides a mechanistic foundation for the fracture‐dependent trajectories that are related to relaxation of fractures with time‐dependent specific stiffness responding to changes in stress or fluid saturation.

02 PETROLEUM↗
DOE OSTI2026

Extended Rice–Thomson analysis and atomistic simulations revealing grain boundary effects on fracture in refractory high-entropy alloys

Significance This work serves to extend the fundamental ductile vs. brittle fracture theory, specifically the Rice–Thomson criterion, by introducing a grain boundary ahead of an initiating crack which propagates at an oblique angle to impinge the boundary. Atomistic fracture simulations on two refractory complex concentrated alloys, the brittle NbMoTaW and the ductile Nb 45 Ta 25 Ti 15 Hf 15 , demonstrate qualitative correspondence with the extended Rice–Thomson criterion and experimental observations. Abstract Understanding how grain boundaries mediate fracture remains a critical challenge in designing ductile, high-performance refractory alloys. Here, we extend the Rice–Thomson criterion to account for the angle between cracks and the impinging grain boundaries (GBs), capturing the competition between intergranular fracture and dislocation-mediated plasticity. Using machine learning interatomic potentials, we performed molecular statics simulations to probe fracture mechanisms in nanocrystalline NbMoTaW and Nb 45 Ta 25 Ti 15 Hf 15 , each with two different grain sizes, revealing trends consistent with experimental observations and the extended Rice model. Comparison with averaged R-curves for bulk samples demonstrates that GBs enhance ductility in Nb 45 Ta 25 Ti 15 Hf 15 in both grain sizes investigated. In contrast, GBs only locally improve fracture resistance in NbMoTaW when cracks are temporarily pinned at GBs inclined at high angles from the crack, but generally promote brittle intergranular fracture. These contrasting behaviors are attributed to differences in GB cohesion, reflecting clear alloying trends that align with ab-initio calculations and trends observed experimentally. Our results bridge classical fracture theory, atomistic simulations, and experimental observations, providing a comprehensive understanding of the fracture mechanisms in nanocrystalline refractory complex concentrated alloys.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Transforming microseismic clouds into near real-time visualization of the growing hydraulic fracture

SUMMARY Microseismic observations during unconventional reservoir stimulation are typically seen as a proxy for clusters of hydraulic fractures and the extent of the stimulated reservoir. Such straightforward interpretation is often misleading and fails to provide a physically reasonable image of the fracturing process. This paper demonstrates the application of a physics-based machine learning algorithm which enables a rapid and accurate fracture mapping from the microseismic data. Our training and validation data set relies on a history-matched geomechanical modelling workflow implemented in GEOS software for the Hydraulic Fracturing Test Site 1 (HFTS-1) project. For this study we augmented the simulated fracture growth through geostatistical modelling of induced seismicity, so that the synthetic microseismic catalogue matches the main statistical properties of the field observations. We formulated the problem of mapping the actual fracture in the clutter of events to parallel common video segmentation workflows: several past video frames (microseismic density snapshots) are passed through a deep convolutional network to classify whether a given voxel is associated with a fracture or intact rock. We found that for accurate fracture mapping, the network’s input and architecture must be augmented to incorporate the fluid injection parameters (pressure, rate, concentration of proppant, and location of the perforation within the cluster). The error rate for the network reached as little as 10 per cent of the fracture area, while a conventional microseismic interpretation approach yielded ∼300 per cent. Our approach also yields must faster predictions than conventional methods (minutes instead of weeks), and could enable engineers to make rapid decisions regarding engineering parameters (pumping rate, viscosity) in real time during stimulation.

58 GEOSCIENCES↗
DOE OSTI2021

Impact Fracture and Fragmentation of Glass via the 3D Combined Finite-Discrete Element Method

A driving technical concern for the automobile industry is their assurance that developed windshield products meet Federal safety standards. Besides conducting innumerable glass breakage experiments, product developers also have the option of utilizing numerical approaches that can provide further insight into glass impact breakage, fracture, and fragmentation. The combined finite-discrete element method (FDEM) is one such tool and was used in this study to investigate 3D impact glass fracture processes. To enable this analysis, a generalized traction-separation model, which defines the constitutive relationship between the traction and separation in FDEM cohesive zone models, was introduced. The mechanical responses of a laminated glass and a glass plate under impact were then analyzed. For laminated glass, an impact fracture process was investigated and results were compared against corresponding experiments. Correspondingly, two glass plate impact fracture patterns, i.e., concentric fractures and radial fractures, were simulated. The results show that for both cases, FDEM simulated fracture processes and fracture patterns are in good agreement with the experimental observations. The work demonstrates that FDEM is an effective tool for modeling of fracture and fragmentation in glass.

42 ENGINEERING↗
DOE OSTI2022

Postmortem and Antemortem Forensic Assessment of Pediatric Fracture Healing from Radiographs and Machine Learning Classification

A timeline of pediatric bone healing using fracture healing characteristics that can be assessed solely using radiographs would be practical for forensic casework, where the fracture event may precede death by days, months, or years. However, the dating of fractures from radiographs is difficult, imprecise, and lacks consensus, as only a few aspects of the healing process are visible on radiographs. Multiple studies in both the clinical and forensic literature have attempted to develop a usable scale to assess pediatric bone healing on radiographs using various healing characteristics. In contrast to the orthopedic definition, a fracture in forensic casework is only considered to be healed when the area around the fracture has been remodeled to the point that the fracture is difficult to detect on a radiograph or on the surface of the bone itself, a process that can take several years. We subjectively assessed visible characteristics of healing in radiograms of fractures occurring in 942 living children and adolescents. By dividing these assessments into learning and test (validation) sets, the accuracy of a newly proposed fracture healing scale was compared to a previous study. Two machine learning models were used to test predictions of the new scale. All three models produced similar estimates with substantial imprecision. Results corroborate the Malone model with an independent dataset and support the efficacy of using less complex models to estimate fracture age in children.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2022

Fractal Characterization of Multimodal, Multiscale Images of Shale Rock Fracture Networks

An array of multimodal and multiscale images of fractured shale rock samples and analogs was collected with the aim of improving the numerical representation of fracture networks. 2D and 3D reconstructions of fractures and matrices span from 10 –6 to 100 m. The origin of the fracture networks ranged from natural to thermal maturation to hydraulically induced maturation. Images were segmented to improve fracture identification. Then, the fractal dimension and length distribution of the fracture networks were calculated for each image dataset. The resulting network connectivity and scaling associations are discussed at length on the basis of scale, sample and order of magnitude. Fracture network origin plays an important role in the resulting fracture systems and their scaling. Rock analogs are also evaluated using these descriptive tools and are found to be faithful depictions of maturation-induced fracture networks.

58 GEOSCIENCES↗
DOE OSTI2024

Fracture Network Influence on Rock Damage and Gas Transport following an Underground Explosion

Simulations of rock damage and gas transport following underground explosions that omit preexisting fracture networks in the subsurface cannot fully characterize the influence of geo-structural variability on gas transport. Previous studies do not consider the impact that fracture network structure and variability have on gas seepage. In this study, we develop a sequentially coupled, axi-symmetric model to look at the damage pattern and resulting gas breakthrough curves following an underground explosion given different fracture network realizations. We simulate 0.327 and 0.164 kT chemical explosives with burial depths of 100 m for 90 stochastically generated fracture networks. Gases quickly reach the surface in 30% of the higher yield simulations and 5% of the lower yield simulations. The fast breakthrough can be attributed to the formation of connected pathways between fractures to the surface. The formation of a connected damage pathway to the surface is not clearly correlated with the fracture intensity (P32) in our simulations. Breakthrough curves with slower transport are highly variable depending on the fracture network sample. The variability in the breakthrough behavior indicates that ignoring the influence of fracture networks on rock damage, which strongly influences the hydraulic properties following an underground explosion, will likely lead to a large underestimation of the uncertainty in the gas transport to the surface. This work highlights the need for incorporation of fracture networks into models for accurately predicting gas seepage following underground explosions.

58 GEOSCIENCES↗
DOE OSTI2023

Using Low-Field Nuclear Magnetic Resonance and X-Ray Computed Microtomography Imaging to Explore Potential of Microbially-Induced Calcium Carbonate Precipitation Treatment to Seal Shale Fractures

Microbially-induced calcium carbonate precipitation (MICP) is a biological process in which microbially produced urease enzymes convert urea and calcium into solid calcium carbonate (CaCO 3 ) deposits. Studies have shown that MICP can be used to seal fractures in shale, raising the possibility of applying this technology to restimulate fracking wells by plugging underperforming fractures. For this and other applications to become a reality, non-invasive tools are needed to determine how effectively MICP seals shale fractures under subsurface conditions. In this study, a 2.54 cm wide and 5.08 cm long Marcellus shale core with a single, ~1 mm wide fracture held open by sand "proppant" underwent MICP-treatment at 60°C until reaching three orders of magnitude permeability reduction. Low-field nuclear magnetic resonance (LF-NMR) and X-Ray computed microtomography (μ-CT) techniques were used to assess the extent of biomineralization within the fracture. These tools revealed that while CaCO 3 precipitation occurred throughout the fracture, there was preferential precipitation around proppant, and the core sealed at the effluent end before filling most of the fracture. Both tools were able to independently calculate of the amount of solid biomineral formed inside the fracture. Furthermore, this study found that the distribution of proppant within the shale fracture was an important parameter controlling the degree of biomineralization.

clastic rock↗
NASA NTRSDate not supplied

Macro-scale Testing and Micromechanics Modeling of Fracture Behaviors for Boron Carbide Composites with Hierarchical Microstructures

With a unique combination of properties including high hardness, low density, wear and corrosion resistance, thermal stability, high neutron absorption, and semi-conductivity, boron carbide (B4C) is a candidate material for various engineering applications that involve extreme environment. The current applications of boron carbide, however, are limited by its intrinsic brittleness due to its strong covalent bonding. To toughen boron carbide, in this work hierarchical microstructure designs was used to provide multiple toughening mechanisms including crack deflection/bridging, micro-crack toughening, etc. Using field assisted sintering, B4C composites with hierarchical microstructure features including graphite platelets, micron and sub-micron sized TiB2 reinforcements were fabricated. The fracture toughness of fabricated B4C composites were previously measured at micro-scale using micro-indentation followed by post-testing microstructure inspection. However, questions including whether the fracture toughness enhancement measured at micro-scale can translate to macro-scale mechanical properties, and what are the fundamental mechanisms behind observed fracture toughness enhancement, remain to be answered. In this study, the fabricated B4C composites were tested using standardized four-point bending method to obtain fracture toughness at macro-scale. In addition, micromechanics modeling was conducted using MAC/GMC code and crack-band model to study the effect of residual stress and weak interphases on fracture behaviors of B4C composites reinforced with TiB2 particles. Through standardized four-point bending tests, fracture toughness enhancements up to 2.85, 3.32, and 3.65 MPa∙m(exp 1/2) (from 2.38 MPa∙m(exp 1/2)) were achieved for B4C composites with graphite platelets addition (micro/nano B4C), with TiB2 formation (micro B4C-TiB2), and with both graphite and TiB2 addition (micro/nano B4C-TiB2) respectively. Micromechanics modeling indicated that introduction of thermal residual stress and weak interphases caused enhanced micro-cracking behavior and resulted in the observed fracture toughness enhancement. These results furthered understanding about the mechanical behaviors at macro-scale and the mechanisms behind observed fracture toughness enhancement for B4C composites with hierarchical microstructures and can provide reference data for the future design of B4C composites with optimized microstructures for further fracture toughness enhancement.

Boron Carbide↗
DOE OSTI2025

A chain stretch-based gradient-enhanced model for damage and fracture in elastomers

Similar to quasi-brittle materials, it has been recently shown that elastomers can exhibit a macroscopically diffuse damage zone that accompanies the fracture process. In this study, we introduce a stretch-based gradient-enhanced damage (GED) model that allows the fracture to localize and also captures the development of a physically diffuse damage zone. This capability contrasts with the paradigm of the phase field method for fracture, where a sharp crack is numerically approximated in a diffuse manner. Capturing fracture localization and diffuse damage in our approach is achieved by considering nonlocal effects that encompass network topology, heterogeneity, and imperfections. These considerations motivate the use of a statistical damage function dependent upon the nonlocal deformation state. From this model, fracture toughness is realized as an output. While GED models have been classically utilized for damage modeling of structural engineering materials (e.g., concrete), they face challenges when trying to capture the cascade from damage to fracture, often leading to damage zone broadening (de Borst and Verhoosel, 2016). This deficiency contributed to the popularity of the phase-field method over the GED model for elastomers and other quasi-brittle materials. Other groups have proceeded with damage-based GED formulations that prove identical to the phase-field method (Lorentz et al., 2012), but these inherit the aforementioned limitations. To address this issue in a thermodynamically consistent framework, we implement two modeling features (a nonlocal driving force bound and a simple relaxation function) specifically designed to capture the evolution of a physically meaningful damage field and the simultaneous localization of fracture, thereby overcoming a longstanding obstacle in the development of these nonlocal strain- or stretch-based approaches. Here, we discuss several numerical examples to understand the features of the approach at the limit of incompressibility, and compare them to the phase-field method as a benchmark for the macroscopic response and fracture energy predictions.

Elastomers↗
DOE OSTI2021

Hydro-mechanical measurements and x-ray imaging of sheared crystalline rock fractures from EGS Collab experiments 1 & 2

Measuring hydro-mechanical properties of natural fractures is a prerequisite for optimizing hydraulic stimulation design and well placement. We completed experiments to characterize shear on natural fractures in schist, amphibolite, and rhyolite specimens drilled from EGS Collab Project’s field sites at the Sanford Underground Research Facility (SURF) in South Dakota. A triaxial direct shear method and coupled x-ray imaging were used to perform hydroshearing and mechanical shearing at the site’s in-situ stress conditions. This produced simultaneous measurements of fracture and matrix strength, permeability, stress-dependent aperture, dilation, and friction strength. Our results identified that only a subset of the natural fractures was weak enough for hydroshearing. Generally, hydroshearing increases fracture permeability by a factor of 10 or more and the enhancement is retainable over time. However, the shear slip does not always result in permeability enhancement. High content of phyllosilicates was found to associate with exceptionally weak fractures that also exhibited poor or even negative enhancement after stimulation. Combining our measurements with site data, we can predict that most observable fractures at the two EGS Collab sites do not meet the criteria for hydroshearing before tensile opening. In some cases, the visible fractures are low permeability and as strong as the adjacent rock. To induce hydroshearing before tensile opening, injection must target known weak and favorably oriented fractures with confirmed pre-existing permeability.

58 GEOSCIENCES↗
DOE OSTI2020

Seismicity and Stress Associated With a Fluid-Driven Fracture: Estimating the Evolving Geometry

A coupled approach, combining the theory of rate- and state-dependent friction and methods from poroelasticity, forms the basis for a quantitative relationship between displacements and fluid leak-off from a growing fracture and changes in the rate of seismic events in the region surrounding the fracture. Poroelastic Green's functions link fracture aperture changes and fluid flow from the fracture to changes in the stress field and pore pressure in the adjacent formation. The theory of rate- and state-dependent friction provides a connection between Coulomb stress changes and variations in the rate of seismic events. Numerical modeling indicates that the Coulomb stress changes can vary significantly between formations with differing properties. The relationship between the seismicity rate changes and the changes in the formation stresses and fluid pressure is nonlinear, but a transformation produces a quantity that is linearly related to the aperture changes and fluid leak-off from the fracture. The methodology provides a means for mapping changes in seismicity into fracture aperture changes and to image an evolving fracture. An application to observed microseismicity associated with a hydrofracture reveals asymmetric fracture propagation within two main zones, with extended propagation in the upper zone. Furthermore, the time-varying volume of the fracture agrees with the injected volume, given by the integration of rate changes at the injection well, providing validation of the estimated aperture changes.

58 GEOSCIENCES↗
DOE OSTI2023

Characterizing the Impacts of Multi-Scale Heterogeneity on Solute Transport in Fracture Networks

We model flow and transport in fracture networks with varying degrees of fracture-to-fracture aperture heterogeneity and network intensity to show how changes in these properties can cause the emergence of anomalous flow and transport behavior. If fracture-to-fracture aperture heterogeneity is increased in sparse networks, velocity fluctuations inhibit high flow rates and solute transport can be delayed. Surprisingly, transport can be slowed even in cases where hydraulic aperture is monotonically increased. As the intensity of the networks is increased, more connected pathways allow for particles to bypass these effects. There exists transition behavior where with relatively few connected pathways in a network, first arrival times of particles are not heavily affected by fracture-to-fracture aperture heterogeneity, but the scaling behavior of the tails is strongly influenced. These results reinforce the importance of considering multi-scale effects in fractured systems and can inform flow and transport processes in both natural and engineered fractured systems.

58 GEOSCIENCES↗
DOE OSTI2023

Emergence of Unstable Focused Flow Induced by Variable‐Density Flows in Vertical Fractures

Abstract Fluids with different densities often coexist in subsurface fractures and lead to variable‐density flows that control subsurface processes such as seawater intrusion, contaminant transport, and geologic carbon sequestration. In nature, fractures have dip angles relative to gravity, and density effects are maximized in vertical fractures. However, most studies on flow and transport through fractures are often limited to horizontal fractures. Here, we study the mixing and transport of variable‐density fluids in vertical fractures by combining three‐dimensional (3D) pore‐scale numerical simulations and visual laboratory experiments. Two miscible fluids with different densities are injected through two inlets at the bottom of a fracture and exit from an outlet at the top of the fracture. Laboratory experiments show the emergence of an unstable focused flow path, which we term a “runlet.” We successfully reproduce the unstable runlet using 3D numerical simulations and elucidate the underlying mechanisms triggering the runlet. Dimensionless number analysis shows that the runlet instability arises due to the Rayleigh‐Taylor instability (RTI), and flow topology analysis is applied to identify 3D vortices that are caused by the RTI. Even under laminar flow regimes, fluid inertia is shown to control the runlet instability by affecting the size and movement of vortices. Finally, we confirm the emergence of a runlet in rough‐walled fractures. Since a runlet dramatically affects fluid distribution, residence time, and mixing, the findings in this study have direct implications for the management of groundwater resources and subsurface applications.

Environmental Sciences & Ecology↗
DOE OSTI2024

Coupled Transport, Reactivity, and Mechanics in Fractured Shale Caprocks

Abstract Shales are low‐permeability caprocks that confine fluid, such as CO 2 , nuclear waste, and hydrogen, in storage formations. Stress‐induced fractures in shale caprocks provide pathways for fluid to leak and potentially contaminate fresh water aquifers. Fractured shales are also increasingly considered as resources for CO 2 sequestration, enhanced geothermal, and unconventional energy recovery. Injecting reactive fluids into shales introduces chemical disequilibrium, causing an onset of a series of dissolution, precipitation, and fines mobilization mechanisms. The reactions have rapid kinetics and significant impact on porosity and permeability; consequently, flow and storage properties of caprocks. While previous research has explored the separate effects of these reactions, this study aims to uncover their simultaneous occurrence and collective influence. This study unveils these highly coupled transport and reactivity mechanisms by tracking and visualizing the reaction‐induced alterations in the matrix, microcracks, and fractures of shales over time. We conducted brine injection experiments sequentially at pH 4 and 2 in a naturally fractured Wolfcamp shale sample while simultaneously imaging the dynamic processes using X‐ray computed tomography (CT). CT images are validated by finer resolution images obtained using micro‐CT and scanning electron microscopy. We also tracked the sample permeability and fluid chemistry using brine permeability and inductively coupled plasma mass spectrometry, respectively. Findings show that fluid primarily flowed through fractures, dissolving reactive minerals and mobilizing fines on fracture surfaces. Dissolution of fracture asperities under confining stress resulted in the closing of fractures. Clogging in narrow fracture pathways, caused by fines accumulation, diverted fluid flow into matrix pores.

58 GEOSCIENCES↗
DOE OSTI2022

Microbial colonization and persistence in deep fractured shales is guided by metabolic exchanges and viral predation

Background: Microbial colonization of subsurface shales following hydraulic fracturing offers the opportunity to study coupled biotic and abiotic factors that impact microbial persistence in engineered deep subsurface ecosystems. Shale formations underly much of the continental USA and display geographically distinct gradients in temperature and salinity. Complementing studies performed in eastern USA shales that contain brine-like fluids, here we coupled metagenomic and metabolomic approaches to develop the first genome-level insights into ecosystem colonization and microbial community interactions in a lower-salinity, but high-temperature western USA shale formation. Results: We collected materials used during the hydraulic fracturing process (i.e., chemicals, drill muds) paired with temporal sampling of water produced from three different hydraulically fractured wells in the STACK (Sooner Trend Anadarko Basin, Canadian and Kingfisher) shale play in OK, USA. Relative to other shale formations, our metagenomic and metabolomic analyses revealed an expanded taxonomic and metabolic diversity of microorganisms that colonize and persist in fractured shales. Importantly, temporal sampling across all three hydraulic fracturing wells traced the degradation of complex polymers from the hydraulic fracturing process to the production and consumption of organic acids that support sulfate- and thiosulfate-reducing bacteria. Furthermore, we identified 5587 viral genomes and linked many of these to the dominant, colonizing microorganisms, demonstrating the key role that viral predation plays in community dynamics within this closed, engineered system. Lastly, top-side audit sampling of different source materials enabled genome-resolved source tracking, revealing the likely sources of many key colonizing and persisting taxa in these ecosystems. Conclusions: These findings highlight the importance of resource utilization and resistance to viral predation as key traits that enable specific microbial taxa to persist across fractured shale ecosystems. We also demonstrate the importance of materials used in the hydraulic fracturing process as both a source of persisting shale microorganisms and organic substrates that likely aid in sustaining the microbial community. Moreover, we showed that different physicochemical conditions (i.e., salinity, temperature) can influence the composition and functional potential of persisting microbial communities in shale ecosystems. Together, these results expand our knowledge of microbial life in deep subsurface shales and have important ramifications for management and treatment of microbial biomass in hydraulically fractured wells.

59 BASIC BIOLOGICAL SCIENCES↗