Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Geologic Modeling”

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 1,081 records · Page 60
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

BioSiting Tool (BioSiting) v2

The BioSiting Tool provides a geospatial interface for analyzing bioeconomy resources and infrastructure across the continental U.S. The tool integrates empirical and modeled data from a broad range of sources. Bioeconomy resources mapped in the tool include agricultural residues, forest residues, municipal solid waste streams, food waste, manure, fats, oils and greases and potential yields of energy crops. Infrastructure mapped in the tool includes biorefineries, material recovery facilities, anaerobic digesters, wastewater treatment plants, combustion plants, district energy systems, crude oil pipelines, petroleum pipelines, natural gas pipelines, railways and freight terminals. Additional data layers include environmental justice indicators at the census tract level and carbon dioxide geologic storage potential. Users can select a location on the map, define a buffer radius in kilometers and generate an inventory of all bioecomony resources within the buffer zone. Data from the tool can be downloaded from individual buffer zones, or at the state or national level.

Huntington, Tyler↗
DOE OSTI2023

Summary of Thermodynamic Data Integration (M4SF-23LL010301052)

This progress report (Level 4 Milestone Number M4SF-23LL010301052) summarizes research conducted at Lawrence Livermore National Laboratory (LLNL) within the Argillite Activity Number M4SF-23LL010301052. SUPCRTNE is being developed as the primary engine for thermodynamic database development in support of geologic disposal of high-level nuclear waste. As used here, “SUPCRT” refers to both a computer program and its supporting database. Additional letters or numbers refer to specific variants (SUPCRT92, Johnson et al., 1992; SUPCRTBL, Zimmer et al., 2016; and SUPCRTNE). A SUPCRT database contains the data required to calculate the thermodynamic properties of solid, gas, and aqueous species over a wide range of temperature and pressure. Normally SUPCRT is used to create a higher-level data base that directly supports modeling and simulation codes such as EQ3/6, GWB, PHREEQC, and PFLOTRAN.

58 GEOSCIENCES↗
DOE OSTI2024

Impacts of irregularly-distributed acidified brine flow on geo-chemo-mechanical alteration in an artificial shale fracture under differential stress

The efficacy of geological carbon sequestration is reliant on the integrity of the caprock and its resistance to physical and chemical alteration. Caprocks with high abundance of reactive carbonates like calcite are susceptible to acid-promoted dissolution and can result in structural weakening. This work investigates the effect of acidified brine flow through an artificially fractured, high-carbonate (30 % by XRD) shale under differential compressive stress. Cylindrical samples were cut in half vertically and milled to create an artificial fracture with interlocking asperities and open channels. Samples were sheared with a single applied stress in a custom flow cell housed within an industrial CT scanner. Further, either acidic (pH 4) or reservoir-simulated (pH 9.5) brine was flowed through the artificial fracture for 7–8 days under reservoir pressure and room temperature. Model simulations indicate flow mainly occurred in open channels, with limited flow between overlapping asperities. Analysis of fracture surfaces by optical and scanning electron microscopy show increased surface alteration and roughness after exposure to pH 4 versus pH 9.5 brine indicating mineral dissolution/loss, and this effect is greater in areas that receive the highest brine flows. Similarly, surface analysis by scratch testing shows fracture toughness decreases more after exposure to acidic versus reservoir-simulated brine, with the greatest alteration in areas of highest acidic brine flows. Despite weakening, no shear slip occurred. Overall, the results indicate that acidified brine can result in significant physical and geomechanical alteration of irregular fracture surfaces in shale caprock, with greatest effects in preferential flow regions.

58 GEOSCIENCES↗
DOE OSTI2021

Atmospheric River Response to Arctic Sea Ice Loss in the Polar Amplification Model Intercomparison Project

The atmospheric river (AR) response to Arctic sea ice loss in the Northern hemisphere winter is investigated using simulations from the Polar Amplification Model Intercomparison Project. Results have shown that the midlatitude responses are dominated by dynamic effects. Poleward of around 60°N, the dynamic and thermodynamic effects cancel each other, resulting in relatively small responses. The response uncertainty can be characterized by leading uncertainty modes, with the responses over the Pacific and Atlantic projecting onto the northeastward extension and equatorward shift mode, respectively. In addition, the responses seem to be mean state-dependent: under the same forcing, models with more poleward-located climatological ARs tend to show stronger equatorward shifts over the Atlantic; over the Pacific, models with more westward-located climatological AR core tend to show stronger northeastward extensions. Finally, these relationships highlight the importance of improving the AR climatology representation on reducing the response uncertainty to Arctic sea ice loss.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Trajectories for the evolution of bacterial CO 2 -concentrating mechanisms

Cyanobacteria rely on CO 2 -concentrating mechanisms (CCMs) to grow in today’s atmosphere (0.04% CO 2 ). These complex physiological adaptations require ≈15 genes to produce two types of protein complexes: inorganic carbon (Ci) transporters and 100+ nm carboxysome compartments that encapsulate rubisco with a carbonic anhydrase (CA) enzyme. Mutations disrupting any of these genes prohibit growth in ambient air. If any plausible ancestral form—i.e., lacking a single gene—cannot grow, how did the CCM evolve? Here, we test the hypothesis that evolution of the bacterial CCM was “catalyzed” by historically high CO 2 levels that decreased over geologic time. Using an E. coli reconstitution of a bacterial CCM, we constructed strains lacking one or more CCM components and evaluated their growth across CO 2 concentrations. We expected these experiments to demonstrate the importance of the carboxysome. Instead, we found that partial CCMs expressing CA or Ci uptake genes grew better than controls in intermediate CO 2 levels (≈1%) and observed similar phenotypes in two autotrophic bacteria, Halothiobacillus neapolitanus and Cupriavidus necator . To understand how CA and Ci uptake improve growth, we model autotrophy as colimited by CO 2 and HCO 3 − , as both are required to produce biomass. Our experiments and model delineated a viable trajectory for CCM evolution where decreasing atmospheric CO 2 induces an HCO 3 − deficiency that is alleviated by acquisition of CA or Ci uptake, thereby enabling the emergence of a modern CCM. This work underscores the importance of considering physiology and environmental context when studying the evolution of biological complexity.

Flamholz, Avi I.↗
DOE OSTI2020

Using Machine Learning to Identify Optimum Prospects for Offshore Geologic CO 2 Storage and Enhanced Oil Recovery in the Gulf of Mexico

The SECARB Offshore Partnership is a government-industry partnership focused on assembling the knowledge base required for secure, long-term, large-scale CO 2 subsea storage in saline formations and mature oil reservoirs. The SECARB Offshore Partnership is evaluating potential storage opportunities in the Central and Eastern planning areas of the outer continental shelf (OCS) of the Gulf of Mexico. The evaluation focuses on active and depleted oil and gas fields and potentially associated CO 2 -enhanced oil recovery (CO 2 -EOR), as well as, deep saline storage resources. As part of Subtask 3.0 (Offshore Storage Resources Characterization) and Subtask 4.0 (Risk Assessment, Simulation, and Modeling), Oklahoma State University (OSU) is developing a machine learning system that employs the SAS Institute’s Viya platform with visual data mining and machine learning to identify and characterize optimum prospects for enhanced oil recovery and stacked storage in saline reservoirs. This report summarizes this effort to date and discusses the many variables that can be used to characterize these storage objectives, as well as, the design vision of the Viya machine learning system.

02 PETROLEUM↗
DOE OSTI2020

From particle attachment to space-filling coral skeletons

Reef-building corals and their aragonite (CaCO 3 ) skeletons support entire reef ecosystems, yet their formation mechanism is poorly understood. Here we used synchrotron spectromicroscopy to observe the nanoscale mineralogy of fresh, forming skeletons from six species spanning all reef-forming coral morphologies: Branching, encrusting, massive, and table. In all species, hydrated and anhydrous amorphous calcium carbonate nanoparticles were precursors for skeletal growth, as previously observed in a single species. The amorphous precursors here were observed in tissue, between tissue and skeleton, and at growth fronts of the skeleton, within a low-density nano- or microporous layer varying in thickness from 7 to 20 µm. Brunauer-Emmett-Teller measurements, however, indicated that the mature skeletons at the microscale were space-filling, comparable to single crystals of geologic aragonite. Nanoparticles alone can never fill space completely, thus ion-by-ion filling must be invoked to fill interstitial pores. Such ion-by-ion diffusion and attachment may occur from the supersaturated calcifying fluid known to exist in corals, or from a dense liquid precursor, observed in synthetic systems but never in biogenic ones. Concomitant particle attachment and ion-by-ion filling was previously observed in synthetic calcite rhombohedra, but never in aragonite pseudohexagonal prisms, synthetic or biogenic, as observed here. Models for biomineral growth, isotope incorporation, and coral skeletons’ resilience to ocean warming and acidification must take into account the dual formation mechanism, including particle attachment and ion-by-ion space filling.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

Endothermic physiology of extinct megatooth sharks

The evolution of the extinct megatooth shark, Otodus megalodon , and its close phylogenetic relatives remains enigmatic. A central question persists regarding the thermophysiological origins of these large predatory sharks through geologic time, including whether O. megalodon was ectothermic or endothermic (including regional endothermy), and whether its thermophysiology could help to explain the iconic shark’s gigantism and eventual demise during the Pliocene. To address these uncertainties, we present unique geochemical evidence for thermoregulation in O. megalodon from both clumped isotope paleothermometry and phosphate oxygen isotopes. Our results show that O. megalodon had an overall warmer body temperature compared with its ambient environment and other coexisting shark species, providing quantitative and experimental support for recent biophysical modeling studies that suggest endothermy was one of the key drivers for gigantism in O. megalodon and other lamniform sharks. The gigantic body size with high metabolic costs of having high body temperatures may have contributed to the vulnerability of Otodus species to extinction when compared to other sympatric sharks that survived the Pliocene epoch.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2023

Fusion of Multiple Models for Improving Gross Primary Production Estimation With Eddy Covariance Data Based on Machine Learning

Abstract Terrestrial gross primary production (GPP) represents the magnitude of CO 2 uptake through vegetation photosynthesis, and is a key variable for carbon cycles between the biosphere and atmosphere. Light use efficiency (LUE) models have been widely used to estimate GPP for its physiological mechanisms and availability of data acquisition and implementation, yet each individual GPP model has exhibited large uncertainties due to input errors and model structure, and further studies of systematic validation, comparison, and fusion of those models with eddy covariance (EC) site data across diverse ecosystem types are still needed in order to further improve GPP estimation. We here compared and fused five GPP models (VPM, EC‐LUE, GOL‐PEM, CHJ, and C‐Fix) across eight ecosystems based on FLUXNET2015 data set using the ensemble methods of Bayesian Model Averaging (BMA), Support Vector Machine (SVM), and Random Forest (RF) separately. Our results showed that for individual models, EC‐LUE gave a better performance to capture interannual variability of GPP than other models, followed by VPM and GLO‐PEM, while CHJ and C‐Fix were more limited in their estimation performance. We found RF and SVM were superior to BMA on merging individual models at various plant functional types (PFTs) and at the scale of individual sites. On the basis of individual models, the fusion methods of BMA, SVM, and RF were examined by a five‐fold cross validation for each ecosystem type, and each method successfully improved the average accuracy of estimation by 8%, 18%, and 19%, respectively.

Environmental Sciences & Ecology↗
DOE OSTI2020

Global ocean dimethyl sulfide climatology estimated from observations and an artificial neural network

Marine dimethyl sulfide (DMS) is important to climate due to the ability of DMS to alter Earth's radiation budget. Knowledge of the global-scale distribution, seasonal variability, and sea-to-air flux of DMS is needed in order to improve understanding of atmospheric sulfur, aerosol/cloud dynamics, and albedo. Here we examine the use of an artificial neural network (ANN) to extrapolate available DMS measurements to the global ocean and produce a global climatology with monthly temporal resolution. A global database of 82 996 ship-based DMS measurements in surface waters was used along with a suite of environmental parameters consisting of latitude–longitude coordinates, time of day, time of year, solar radiation, mixed layer depth, sea surface temperature, salinity, nitrate, phosphate, and silicate. Linear regressions of DMS against the environmental parameters show that on a global-scale mixed layer depth and solar radiation are the strongest predictors of DMS. These parameters capture ~9 % and ~7 % of the raw DMS data variance, respectively. Multilinear regression can capture more of the raw data variance (~39 %) but strongly underestimates DMS in high-concentration regions. In contrast, the artificial neural network captures ~66 % of the raw data variance in our database. Like prior climatologies our results show a strong seasonal cycle in surface ocean DMS with the highest concentrations and sea-to-air fluxes in the high-latitude summertime oceans. We estimate a lower global sea-to-air DMS flux (20.12±0.43 Tg S yr –1 ) than the prior estimate based on a map interpolation method when the same gas transfer velocity parameterization is used. Our sensitivity test results show that DMS concentration does not change unidirectionally with each of the environmental parameters, which emphasizes the interactions among these parameters. The ANN model suggests that the flux of DMS from the ocean to the atmosphere will increase with global warming. Given that larger DMS fluxes induce greater cloud albedo, this corresponds to a negative climate feedback.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

A comparative analysis of numerical approaches for the description of gas flow in clay-based repository systems: From a laboratory to a large-scale gas injection test

There is nowadays a consensus among many countries that geological disposal is a favourable solution for the long-term management. Although different host formations and different barrier systems are under consideration around the world, clay-based materials form an important component for waste isolation in most national programmes. Hence, a good comprehension of the effect of gas flow on the hydro-mechanical behaviour of clay-based soils is essential, both at laboratory and field scale. Task B under the international cooperative project DECOVALEX-2023 has recently shown that, after some enhancement, models can be employed to reproduce laboratory scale tests, even with different sample geometries 37 . However, further work is required to understand whether they can be applied to simulate a large-scale experiment. Up-scaling of models for the advective transport of gas through clay-based low permeable material presents a number of problems related to the difficulty in obtaining consistent hydrogeological parameters and constitutive relationships at both laboratory and field scale. Based on a unique dataset from a large-scale gas injection test (Lasgit) performed at the Äspö Hard Rock Laboratory (Sweden), Task B within DECOVALEX-2023 has explored the refinement of these numerical strategies applied to the simulation of gas flow. Work performed within the task reveals that codes do not need to be substantially modified from the laboratory models to reproduce full-scale tests: indeed, model parameters calibrated and validated at laboratory scale have been applied to predict field scale gas flow at Lasgit, including peak gas pressure and injected cumulative gas volume. By means of (1) the introduction of interfaces between blocks to reflect the experimental configuration and the (2) adjustment of some parameters (e.g., higher permeability), the updated models are able to represent most of the key features observed in the experimental data, even at a large scale.

Tamayo-Mas, E↗
DOE OSTI2022

Role of copper in nitrous oxide accumulation in wetland soils from Illinois and South Carolina and stream sediments from Tennessee

Denitrification is microbially-mediated through enzymes containing metal cofactors. Laboratory studies of pure cultures have highlighted that the availability of copper (Cu), required for the multicopper enzyme nitrous oxide reductase, can limit nitrous oxide (N2O) reduction. However, in natural aquatic systems, such as wetlands and hyporheic zones in stream beds, the role of Cu in controlling denitrification remains incompletely understood. In this study, we collected soils and sediments from three natural environments -- riparian wetlands, marsh wetlands, and a stream -- to investigate their nitrogen species transformation activity at background Cu levels and different supplemented Cu loadings. All of the systems contained solid-phase associated Cu below or around geological levels (40–280 nmol g-1) and exhibited low dissolved Cu (3–50 nM), which made them appropriate sites for evaluating the effect of limited Cu availability on denitrification.The dataset contains the variation in the nutrient (nitrate (NO3-), nitrite (NO2-), and ammonium (NH4+)) and N2O concentrations during the incubation experiments to evaluate the effect of Cu. We have also reported the variation in metal concentrations (Cu, iron(Fe) and manganese (Mn)) and dissolved organic carbon (DOC) concentrations during the incubations. The parameters obtained using a kinetic model to quantitatively report the effect of Cu on nitrogen species conversion in natural aquatic systems are also added in the dataset. The dataset also contains the labile concentration of Cu estimated using a speciation model in MINTEQ. All data are provided in text-based CSV format with header sections indicating the data contained in each file and the corresponding units.Our study suggests that high concentrations of N2O accumulated in all microcosms lacking Cu amendment except for one stream sediment sample. With Cu added to provide dissolved concentrations at trace levels (10–300 nM), the reduction rate of N2O to N2 in the wetland soils and stream sediments was enhanced. A kinetic model could account for the trends in nitrogen species by combining the reactions for microbial reduction of NO3- to NO2-/N2O/N2 and abiotic reduction of NO2 to nitrogen (N2_. The model revealed that the rate of N2O to N2 conversion increased significantly in the presence of Cu. For riparian wetland soils and stream sediments, the kinetic model also suggested that overall denitrification is driven by abiotic reduction of NO2- in the presence of inorganic electron donors. This study demonstrated that natural aquatic systems containing Cu at concentrations less than or equal to crustal abundances may display incomplete reduction of N2O to N2 that would cause N2O accumulation and release to the atmosphere.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Wellbore injectivity response to step-rate CO 2 injection: Coupled thermo-poro-elastic analysis in a vertically heterogeneous formation

Safe and permanent carbon geological storages require an elaborate analysis of geomechanical stability. Subsurface injection of CO 2 changes the local temperature and pore pressure and, further, alters the stress due to thermo-poro-elastic responses. A permanent CO 2 storage testing was conducted in the water leg of Cranfield reservoir in Mississippi, USA, where hosted CO 2 enhanced-oil recovery activities. During CO 2 injection, the injection rate was ramped up twice but the bottom-hole pressure did not increase with the imposed injection rates as expected. This unexpected field observation suggests the possibility of an open-mode fracture development at the injector. However, the injector response and potential fracture development have not been rigorously interpreted with detailed near-wellbore temperature and pore pressure change upon the CO 2 injection. In this study, we performed history matching of CO 2 injection using coupled thermo-poro-elastic reservoir simulation to examine the possibility of fracturing. We built a reservoir model including vertical heterogeneity in both petrophysical and geomechanical properties estimated from well-logging analysis and laboratory experiments. The simulation results show that CO 2 injection changes stresses and support the hypothesis of development of an open-mode fracture at the injector during the Cranfield test. The near-injector region exhibits a large temperature reduction up to 55 °C with ensuing effective horizontal stress reduction up to 9.1 MPa. However, a caprock integrity issue is unlikely because of the horizontal stress contrast within layers and low hydraulic communication with the injection zone. This work indicates that the injectant temperature should be considered in the design of high-rate CO 2 injector.

58 GEOSCIENCES↗
DOE OSTI2023

Preservation of corneous β-proteins in Mesozoic feathers

Fossil proteins are valuable tools in evolutionary biology. Recent technological advances and better integration of experimental methods have confirmed the feasibility of biomolecular preservation in deep time, yielding new insights into the timing of key evolutionary transitions. Keratins (formerly α-keratins) and corneous β-proteins (CBPs, formerly β-keratins) are of particular interest as they define tissue structures that underpin fundamental physiological and ecological strategies and have the potential to inform on the molecular evolution of the vertebrate integument. Reports of CBPs in Mesozoic fossils, however, appear to conflict with experimental evidence for CBP degradation during fossilization. Further, the recent model for molecular modification of feather chemistry during the dinosaur–bird transition does not consider the relative preservation potential of different feather proteins. Here we use controlled taphonomic experiments coupled with infrared and sulfur X-ray spectroscopy to show that the dominant β-sheet structure of CBPs is progressively altered to α-helices with increasing temperature, suggesting that (α-)keratins and α-helices in fossil feathers are most likely artefacts of fossilization. Our analyses of fossil feathers shows that this process is independent of geological age, as even Cenozoic feathers can comprise primarily α-helices and disordered structures. Critically, our experiments show that feather CBPs can survive moderate thermal maturation. Importantly, as predicted by our experiments, analyses of Mesozoic feathers confirm that evidence of feather CBPs can persist through deep time.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

2D-to-3D image translation of complex nanoporous volumes using generative networks

Image-based characterization offers a powerful approach to studying geological porous media at the nanoscale and images are critical to understanding reactive transport mechanisms in reservoirs relevant to energy and sustainability technologies such as carbon sequestration, subsurface hydrogen storage, and natural gas recovery. Nanoimaging presents a trade off, however, between higher-contrast sample-destructive and lower-contrast sample-preserving imaging modalities. Furthermore, high-contrast imaging modalities often acquire only 2D images, while 3D volumes are needed to characterize fully a source rock sample. In this work, we present deep learning image translation models to predict high-contrast focused ion beam-scanning electron microscopy (FIB-SEM) image volumes from transmission X-ray microscopy (TXM) images when only 2D paired training data is available. We introduce a regularization method for improving 3D volume generation from 2D-to-2D deep learning image models and apply this approach to translate 3D TXM volumes to FIB-SEM fidelity. We then segment a predicted FIB-SEM volume into a flow simulation domain and calculate the sample apparent permeability using a lattice Boltzmann method (LBM) technique. Results show that our image translation approach produces simulation domains suitable for flow visualization and allows for accurate characterization of petrophysical properties from non-destructive imaging data.

58 GEOSCIENCES↗
DOE OSTI2020

Ordovician-Cambrian Units: Hierarchical Evaluation of Geologic Carbon Storage Resource Estimates

The Indiana Geological and Water Survey (IGWS) led subtask 1.1 to assess the regional distribution and estimate the storage capacity of Ordovician-Cambrian stratigraphic units located within the partnership region. A comprehensive data set of wireline logs and petrophysical information was used to generate these interpretations. These data include core analysis for porosity and permeability, mercury injection capillary pressure (MICP), and existing well data including location and stratigraphic information. This report includes storage resource estimates (SREs) for three potential storage reservoirs(limestone and dolostone from the Upper Ordovician Trenton Limestone/Black River Group and equivalent units, the Middle Ordovician St. Peter Sandstone, and primary target reservoir rocks of the Lower Ordovician and Upper Cambrian Knox Supergroup and equivalent units) calculated using six methodologies: (1) a fixed value of porosity of 10 percent in all units evaluated; (2) a unique average porosity (per well) from wireline-derived porosity (neutron, sonic, and/or density porosity for each unit); (3) porosity values from core analysis; (4) a depth-dependent porosity model (Knox Supergroup only); (5) porosity based on a model based on petrophysical facies; and (6) SREs using National Energy Technology Laboratory’s CO2 Storage prospeCtive Resource Estimation Excel aNalysis (CO2-SCREEN beta V2). All methods used the same values for thickness for each unit. However, the areal extent of each assessment was limited by the data available for each method. Estimated volumes were calculated in 1-by-1 kilometer grid cells and summarized as county and total stratigraphic unit volumes. The resultant SREs mass are displayed using boxplots, which allow for comparing data statistics (mean values and variability) between methods. Differences observed in SRE results from the six methods are mainly attributable to differences in the data and conceptual models used to interpret or estimate porosity in each method. Based on this systematic variability between methods, it is inferred that methods 1, 4, and 6 are best used for regional-scale reconnaissance estimates of storage capacity while methods 2, 3, and 5 are more appropriate for local scales where more data is required. All estimates are data-density dependent and different methods require different amounts of data for reasonable assessments. ArcMap 10.5.1 software was used to portray SREs to help visualize spatial variance of estimates for each methodology, and more importantly, to highlight those areas having the greatest total storage potential estimates.

01 COAL, LIGNITE, AND PEAT↗
DOE OSTI2025

2032 and 2036 risk enhancement from NEOs in the Taurid stream: Is there a significant coherent component to impact risk?

Impact risk is normally quantified by summing the product of the probability of an event and some measure of its consequences over the set of all possible events. The probability factor is considered to be more objective and is based on the size frequency distribution of Near Earth objects (NEOs) and an implicit assumption of randomness, which can be described as “stochastic risk”. Impact frequency does change with time, however, and there have been episodes in the deep geological past when the flux has been much higher. The hypothesis of “coherent catastrophism” suggests large variations on shorter timescales. It postulates the existence of a “Taurid resonant swarm” (TRS) of debris associated with Comet Encke that is stabilized by Jupiter and in a 7:2 resonance with it. The hypothetical cluster orbits in the broad Taurid stream, which crosses Earth's orbit twice a year at its nodes. Whereas the most extreme and fanciful versions of coherent catastrophism in recent geologic history have been comprehensively refuted, the possibility of a significant component of coherent (time-dependent) risk, associated with non-random correlated events, remains. This paper incorporates recently published data from observational campaigns associated with the Taurid stream. There is no evidence for objects in the Taurid stream that are above the global catastrophe threshold, but the possibility of a few large objects and a significant population of objects the size of the Chelyabinsk or Tunguska bodies has not been eliminated. Eyewitness accounts and comparison of airburst models to ground truth suggest that the Tunguska object was a Beta Taurid. The Tunguska event is probably an outlier when compared to the size frequency distributions under stochastic assumptions, but if there is a significant coherent component it may be representative of events that take place far more frequently, and risk assessments may have underestimated the contribution from airbursts. If so, the Earth will experience increased risk from objects this size, peaking in November 2032 and June 2036. We recommend targeted survey campaigns during these hypothetical TRS node crossings to quantify the population, search for potentially hazardous objects, and identify imminent impactors.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

Geologic stress modulates fluid mixing at fracture intersections

Fracture intersections are critical links that enable flow and transport in subsurface fracture networks, and their behavior strongly influences fluid mixing in a network. Although all subsurface fractures are subjected to geological stress, we lack a fundamental understanding of how fracture intersection geometry evolves under stress and how these changes influence fluid mixing. Here, we combine 3D printing, 3D X-ray tomographic imaging, and 3D pore-scale numerical simulations to reveal stress-induced changes in intersection geometry and their impact on mixing. Mixing is found to be strongly affected by partial closure of an intersection under stress. As an intersection closes, the void area for fluid flow and diffusion decreases leading to substantial deviations between conventional mixing models and full pore-scale modeling. To address this, we propose a modified mixing model that accounts for intersection deformation, which is essential for accurate modeling of solute transport and mixing through fracture networks.

15 GEOTHERMAL ENERGY↗