Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “reservoir engineering”

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 379 records · Page 21

Hydropower Resilience Database for Assessing Microgrid Formation Capability and Enhancing Power Grid Resilience

In the face of increasing frequency of extreme events, enhancing resilience and reliability of energy infrastructure demands innovative solutions. Hydropower, with its inherent generation flexibility and grid-forming capability, offers significant potential for enhancing power grid resilience through the establishment of microgrids. To enable informed decision-making and strategic planning, we present the Hydropower Resilience Database (HRD) that integrates data from various sources such as Oakridge National Laboratory (ORNL) HydroSource, National Inventory of Dams, and US Western Grid database. By integrating relevant information on hydropower plant characteristics, including dams, reservoirs, and connected electrical grid, this resource enables hydropower plant owners, utilities, and community stakeholders to identify and evaluate the feasibility of using hydropower resources in microgrids to support nearby communities and critical infrastructure in various challenging scenarios. Using HRD, a set of metrics are evaluated to quantify the capability of hydropower plants to support essential microgrid functions. An interactive tool is developed using ArcGIS, enabling the visualization and analysis using HRD. Our research contributes to strategic planning efforts for fortifying energy infrastructure, ensuring reliable power supply during disruptions, and advancing the development of robust and resilient energy systems in regions susceptible to grid vulnerabilities.

13 HYDRO ENERGY↗

Hydropower Resilience Database for Assessing Microgrid Formation Capability and Enhancing Power Grid Resilience

In the face of increasing frequency of extreme events, enhancing resilience and reliability of energy infrastructure demands innovative solutions. Hydropower, with its inherent generation flexibility and grid-forming capability, offers significant potential for enhancing power grid resilience through the establishment of microgrids. To enable informed decision-making and strategic planning, we present the Hydropower Resilience Database (HRD) that integrates data from various sources such as Oakridge National Laboratory (ORNL) HydroSource, National Inventory of Dams, and US Western Grid database. By integrating relevant information on hydropower plant characteristics, including dams, reservoirs, and connected electrical grid, this resource enables hydropower plant owners, utilities, and community stakeholders to identify and evaluate the feasibility of using hydropower resources in microgrids to support nearby communities and critical infrastructure in various challenging scenarios. Using HRD, a set of metrics are evaluated to quantify the capability of hydropower plants to support essential microgrid functions. An interactive tool is developed using ArcGIS, enabling the visualization and analysis using HRD. Our research contributes to strategic planning efforts for fortifying energy infrastructure, ensuring reliable power supply during disruptions, and advancing the development of robust and resilient energy systems in regions susceptible to grid vulnerabilities.

13 HYDRO ENERGY↗

Evaluating Production Implications of Pressure Maintenance in Unconventional Oil and Gas Wells using a Machine Learning Modeling Approach: Case Study in the Permian Basin

This study implements the proprietary machine learning-based model (model) developed under the 2022 study titled “Evaluating the Impact of Proprietary Oil & Gas Data on Machine Learning Model Performance Using a Quasi-Experimental Analytical Approach” for forecasting unconventional oil and gas production using well data from the Permian Basin. The model, developed using an exclusive dataset that includes time series production data from an operator in the Permian Basin, is designed to jointly predict daily oil, gas, and water production for horizontal wells as a function of bottomhole pressure drawdown, spatial placement across the study domain, and well-completion attributes. In this study, the model was explicitly applied to explore its utility to evaluate the impact of varying drawdown strategies on the production forecast of a well in the Permian Basin dataset. Managing pressure drawdown has been identified as a way to improve estimated ultimate recovery (EUR) from unconventional shale wells due to the stress-dependent nature of fractures in shale reservoirs. Research has shown that applying a lower pressure drawdown helps to maintain the reservoir conductivity, resulting in higher productivity over the life of a well. Historic bottomhole pressure data from the well over time was used as a benchmark from which to set more and less aggressive pressure decline rates as bounding modeling cases. All pressure decline rates/strategies were forecasted over 5 years, and the model was used to generate oil, water, and gas prediction over the same timeframe. This report presents the results from the production forecast and discusses potential operational and economic implications, with contrasting perspectives between well productivity and profitability given typical oil and gas economics and the volatility in the oil and gas market.

02 PETROLEUM↗

Subsurface Structure in the Martian Polar Layered Deposits: The Deep Space 2 Impact Accelerometry Experiment

While primarily a technology demonstration mission, the New Millenium Mars Microprobes (also known as Deep Space 2, or simply DS2)will also provide the first in situ science measurements of the martian subsurface. The DS2 impact accelerometry experiment will provide both engineering data about the depth of probe emplacement and science data about the physical nature of the subsurface at the probes' landing sites. Little is known about the detailed physical properties or small-scale vertical structure of the subsurface at the DS2 landing site in the southern martian polar layered deposits. Imaging data from the Viking Orbiters and Mars Global Surveyor reveal alternating bands of light and dark material in this region with thicknesses at least as small as the limit of resolution, about 10 m. The overall composition of these layers is poorly constrained, but generally thought to be a mixture of dust and ice with the layers being caused by variations in the dust/ice ratio, or perhaps by dust deposits of different densities. Low thermal inertias in the region suggest that the top few centimeters of the surface are composed of a mantling of fine-grained dust. However, 3.5-cm radar returns indicate that the maximum depth of this dust is not greater than a few tens of centimeters. Thermal models generally agree that, while the layered deposits do provide a potential near-surface reservoir for ice, the uppermost few centimeters to meters in these regions are likely to be ice-free because of sublimation losses. Finally, while it is generally agreed that the layered deposits are the product of variations in the martian climate, no direct correlation has been made between band sequences and specific climate changes. Our intention is to shed light on some of these questions about the martian polar layered deposits by using the DS2 accelerometry experiment to determine the physical nature of the layered deposits, and to detect the presence of any subsurface layering of dust, ice, and/or rock. In the process, we will also determine the final resting depth of the two microprobes, an important parameter in the interpretation of other DS2 experiments.

Moersch, J. E.↗

Maximizing the Proppant Carrying and Viscoelastic Properties of the Bakken Hypersaline-Produced Water with High-Viscosity Friction Reducers for Sustainable Applications

Summary The development and production of unconventional reservoirs, such as the Bakken Formation, have become a resolved mystery for operators in North America since the arrival and advancement of horizontal drilling and hydraulic fracturing technologies. As a result, unconventional reservoir assets became the central focus of the oil and gas industry at the state, national, and global levels. The produced water from these activities in the Bakken Formation have high salt contents (110,000–350,000 ppm) total dissolve solids (TDS) and can pose significant challenges to the environment if not treated. Deep injection into disposal wells is the routine method used to get rid of the Bakken produced water. However, there have been some concerns that unrestrained injections, in addition to polluting the groundwater, could potentially lead to seismic activities either at the time of injection or in the near future. To diminish the environmental impacts that may be associated with induced seismicity, including the reduction of the costs of water acquisition, the produced water can be treated and reused in the hydraulic fracturing processes. Also, the treated water could be used for irrigation purposes, for power generation, and coal mining operations. The issues of waste water and residual oil high in TDS are challenges yet to be effectively addressed despite preceding research and studies on advancing produced water technologies. The goal of this study is to explore all applicable ways by which the produced water from the Bakken Formation can serve as a replacement base fluid for use with polymers like the high-viscosity friction reducers (HVFRs) to create hydraulic fracturing fluids that can be stable with reservoir conditions and also be able to minimize environmental impacts and cost of operations. Experimental investigations using the high-salinity produced water from the Bakken Formation with HVFRs were carried out. The studies included a base case that served as a bench mark for comparing the effectiveness of the other scenarios. The results indicate that the Bakken hyper-saline produced water can withstand effect of heavy metals, salinity, hardness and remain stable through different shear rates (66–330 s−1) when treated with higher dosages [4–8 gal/1,000 gal (gpt)] of HVFRs. Filtration and dilution were the only methods used on the Bakken Formation produced wate for this research

Engineering↗

Benchmarking the Performance of Neuromorphic and Spiking Neural Network Simulators

Software simulators play a critical role in the development of new algorithms and system architectures in any field of engineering. Neuromorphic computing, which has shown potential in building brain-inspired energy-efficient hardware, suffers a slow-down in the development cycle due to a lack of flexible and easy-to-use simulators of either neuromorphic hardware itself or of spiking neural networks (SNNs), the type of neural network computation executed on most neuromorphic systems. While there are several openly available neuromorphic or SNN simulation packages developed by a variety of research groups, they have mostly targeted computational neuroscience simulations, and only a few have targeted small-scale machine learning tasks with SNNs. Evaluations or comparisons of these simulators have often targeted computational neuroscience-style workloads. In this work, we seek to evaluate the performance of several publicly available SNN simulators with respect to non-computational neuroscience workloads, in terms of speed, flexibility, and scalability. We evaluate the performance of the NEST, Brian2, Brian2GeNN, BindsNET and Nengo packages under a common front-end neuromorphic framework. Our evaluation tasks include a variety of different network architectures and workload types to mimic the computation common in different algorithms, including feed-forward network inference, genetic algorithms, and reservoir computing. We also study the scalability of each of these simulators when running on different computing hardware, from single core CPU workstations to multi-node supercomputers. Our results show that the BindsNET simulator has the best speed and scalability for most of the SNN workloads (sparse, dense, and layered SNN architectures) on a single core CPU. However, when comparing the simulators leveraging the GPU capabilities, Brian2GeNN outperforms the others for these workloads in terms of scalability. NEST performs the best for small sparse networks and is also the most flexible simulator in terms of reconfiguration capability NEST shows a speedup of at least 2x compared to the other packages when running evolutionary algorithms for SNNs. The multi-node and multi-thread capabilities of NEST show at least 2x speedup compared to the rest of the simulators (single core CPU or GPU based simulators) for large and sparse networks. We conclude our work by providing a set of recommendations on the suitability of employing these simulators for different tasks and scales of operations. We also present the characteristics for a future generic ideal SNN simulator for different neuromorphic computing workloads.

97 MATHEMATICS AND COMPUTING↗

Thermodynamics of quantum measurement and Maxwell's demon's arrow of time

Here, we discuss the thermodynamic aspects of a single qubit based device, powered by weak quantum measurements, and feedback controlled by a quantum Maxwell's demon. We discuss both discrete and time-continuous operation of the measurement based device at finite temperature of the reservoir. In the discrete example where a demon acquires information via discrete weak measurements, we find that the thermodynamic variables including the heat exchanged, extractable work, and the entropy produced are completely determined by an information theoretic measure of the demon's perceived arrow of time. We also discuss a realistic time-continuous operation of the device where the feedback is applied after a sequence of weak measurements. In the time-continuous limit, we derive the exact finite-time statistics of work, heat and entropy changes along individual quantum trajectories of the quantum measurement process, and relate them to the demon's arrow of time.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Assessment of chemo-mechanical impacts of CO 2 sequestration on the caprock formation in Farnsworth oil field, Texas

This study evaluates the chemo-mechanical influence of injected CO 2 on the Morrow B sandstone reservoir and the upper Morrow shale caprock utilizing data from the inverted 5-spot pattern centered on Well 13-10A within the Farnsworth unit (FWU). This study also seeks to evaluate the integrity of the caprock and the long-term CO 2 storage capability of the FWU. The inverted 5-spot pattern was extracted from the field-scale model and tuned with the available field observed data before the modeling work. Two coupled numerical simulation models were utilized to continue the study. First, a coupled hydro-geochemical model was constructed to simulate the dissolution and precipitation of formation minerals by modeling three intra-aqueous and six mineral reactions. In addition, a coupled hydro-geomechanical model was constructed and employed to study the effects of stress changes on the caprock’s porosity, permeability, and ground displacement. The Mohr–Coulomb circle and failure envelope were used to determine caprock failure. In this work, the CO 2 -WAG injection is followed by the historical field-observed strategy. During the forecasting period, a Water Alternating Gas (WAG) injection ratio of 1:3 was utilized with a baseline bottom-hole pressure constraint of 5500 psi for 20 years. A post-injection period of 1000 years was simulated to monitor the CO 2 plume and its effects on the CO 2 storage reservoir and caprock integrity. The simulation results indicated that the impacts of the geochemical reactions on the porosity of the caprock were insignificant as it experienced a decrease of about 0.0003% at the end of the 1000-year post-injection monitoring. On the other hand, the maximum stress-induced porosity change was about a 1.4% increase, resulting in about 4% in permeability change. It was estimated that about 3.3% of the sequestered CO 2 in the formation interacted with the caprock. Despite these petrophysical property alterations and CO 2 interactions in the caprock, the caprock still maintained its elastic properties and was determined to be far from its failure.

04 OIL SHALES AND TAR SANDS↗

Geochemical evolution in Cacapon member: Fluid-rock interaction experiments and model insights for Appalachian Basin geothermal development

Here, this study combines recirculated flow-through experimental results conducted for 17 days at 90C and 200 PSI with reactive transport modeling to estimate fluid-rock interactions occurring in a sandy mudstone using an interbedded sandstone-shale sample from the Cacapon Member of the upper Tuscarora Sandstone/lower Rose Hill Formation for the purpose of geothermal exploration. Results suggest that the fluid and rock are likely to be in or near partial-equilibrium after approximately one year. In addition, after >400 h of continuous injection at 0.05 mL/min (5*10 −8 m 3 /min) the reactive front is restricted to the first ∼13% of the 4 cm*1.6 cm experimental rock length, whereas after >9100 h (∼1 year), the reactive front extends to 30% total length. The rate of changes in dissolution or precipitation are however, very minimal, with all major rock-forming minerals having rates <10 −11 mol/L porous media/s. Reservoir rock in the presence of dilute brine as may be the case during operation of an enhanced geothermal system would experience little alteration during the shut-in phase, and possibly up to one year. These results have utility in geothermal exploration for reservoirs at similar temperatures as well as general fluid-sandy mudstone rock interaction geochemistry.

Appalachian Basin↗

Response surface modeling of CO 2 dynamic storage efficiency factor in high permeability thick sandstones

In order to utilize as much of the pore space for CO 2 storage in high permeability thick saline aquifers, it is vital to investigate the interactions of injected CO 2 with formation brine and rock. In order to quantify the displacement process, we investigate the dynamic storage efficiency factor (DSEF) for saline aquifers where pressure increase is minimal during the injection phase. Dimensionless numbers are derived from basic governing equations, constitutive equations, initial and boundary conditions using the inspection analysis. Then using the Hammersley sequence sampling, 178 numerical experiments are designed, and a compositional reservoir simulator is used to perform these simulations. In the next step, response surface regression analysis is used to establish a relationship between DSEF obtained from the numerical simulations and the corresponding dimensionless numbers. The simulation results show that for the studied conditions the underground dynamics is mostly influenced by the gravity number, followed by effective aspect ratio and dip numbers. The results from the response surface regression analysis are used to develop a correlation, which can be used to estimate the dynamic CO 2 storage capacity of relevant zones. This study provides quantitative measures for the different competing mechanisms involved in underground displacement of fluids in CO 2 geological storage, which can serve as a useful tool during planning phase of storage projects.

02 PETROLEUM↗

Assessment of CO 2 trapping mechanisms in partially depleted oil-bearing sands

The goal of this work is to evaluate various CO 2 sequestration mechanisms occurring in the Morrow B Sandstone in the Farnsworth Unit. A history-matched numerical simulation model was created using extensive geological, petrophysical, and operational data collected from the field. The numerical model is competent to investigate the impact of residual, structural-stratigraphic, solubility, and mineral trapping mechanisms on the fluid transportation dynamics and petrophysical property variations. The model forecasts the field response of 20 years of WAG injections. Afterward, all wells were shut-in, and the reservoir was allowed to evolve for 1000 years to investigate the fate of injected CO 2 . In this paper, we assess the impacts of various trapping mechanisms on oil recovery and CO 2 storage efficacy. By analyzing the results reported from the numerical simulation model, the in-situ fluid composition and mineralogy changes are also investigated. More importantly, we seek to confirm the petrophysical property variations due to the CO 2 injection with observations from laboratory measurements. The experiences gained from this study provide valuable insights regarding physiochemical storage induced by the CO 2 injection activities and serve as a benchmark case for future CO 2 enhanced oil recovery (EOR) projects involving reactive solute transport.

42 ENGINEERING↗

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↗

The Multistage Compressor Facility

Research and developments of new aerospace technologies is one of Glenn Research Center's specialties. One facility that deals with the research of aerospace technologies is the High-speed Multistage Compressor Facility. This facility will be testing the performance and efficiency of an Ultra Efficient Engine Technology (UEET) two-stage compressor. There is a lot of preparation involved with testing something of this caliber. Before the test article can be installed into the test rig, the facility must be fully operational and ready to run. Meaning all the necessary instrumentation must be calibrated and installed in the facility. The test rig should also be in safe operating condition, and the proper safety permits obtained. In preparation for the test, the Multistage Compressor Facility went through a few changes. For instance the facility will now be utilizing slip rings, the gearbox went through some maintenance, new lubrications systems replaced the old ones, and special instrumentation needs to be fine tuned to achieve the maximum amount of accurate data. Slips rings help gather information off of a rotating device - in this case from a shaft - onto stationary contacts. The contacts (or brushes) need to be cooled to reduce the amount of frictional heat produced between the slip ring and brushes. The coolant being run through the slip ring is AK-225, a material hazardous to the ozone. To abide by the safety regulations the coolant must be run through a closed chiller system. A new chiller system was purchased but the reservoir that holds the coolant was ventilated which doesn t make the system truly closed and sealed. My task was to design and have a new reservoir built for the chiller system that complies with the safety guidelines. The gearbox had some safety issues also. Located in the back of the gearbox an inching drive was set up. When the inching drive is in use the gears and chain are bare and someone can easily get caught up in it. So to prevent anyone from getting hurt in the gears I designed a chain guard. Additional information is included in the original extended abstract.

Flegel, Ashlie↗

Data Arrays for Microearthquake (MEQ) Monitoring using Deep Learning for the Newberry EGS Sites

The 'Machine Learning Approaches to Predicting Induced Seismicity and Imaging Geothermal Reservoir Properties' project looks to apply machine learning (ML) methods to Microearthquake (MEQ) data for imaging geothermal reservoir properties and forecasting seismic events, in order to advance geothermal exploration and safe geothermal energy production. As part of the project, this submission provides data arrays for 149 microearthquakes between the year 2012 and 2013 at the Newberry EGS Site for use with the Deep Learning Algorithm that has been developed. The data provided includes raw waveform data, location data, normalized waveform data, and processed waveform data. Penn State Geothermal Team has shared the following files from the project: - 149 microearthquakes (MEQs) between 2012 and 2013 at Newberry EGS sites, 'Normalized Waveform Inputs.npz' are normalized waveforms. - labels of 149 MEQs: Processed Waveform Inputs.npz - location labels of 149 MEQs: Location Data.npz Note: .npz is the python file format by NumPy that provides storage of array data.

15 GEOTHERMAL ENERGY↗

Feasibility of using in situ deformation to monitor CO 2 storage

Deformation during CO 2 injection can lead to problems, like seismicity or fluid leaks, but small strains have the potential to be a useful signal for monitoring. The objective of this paper is to evaluate the possible evolution of the strain field during injection, and then assess existing and emerging techniques for measuring the strain field. Poroelastic analyses show that normal strains caused by injection into a reservoir are tensile in the vicinity of the well, but everywhere else at least one strain component is compressive. The vertical strain is compressive in the confining unit, and the radial strain decreases and changes sign from tensile to compressive with distance from the well. Tilting is away from the injection well at the ground surface, but it is towards the well overlying the reservoir. Methods for measuring in-situ strain include instruments that are grouted in the annulus between casing and wall rock (~ 0.1 microstrain resolution), portable strain sensors that are temporarily clamped to the borehole wall (~0.01 microstrain resolution), and strainmeters that are grouted in place (~0.001 microstrain resolution). Instruments for measuring in-situ normal strains at the magnitudes and rates expected during injection are emerging, but they have yet to be fully evaluated in applications related to CO 2 storage. Finally, in-situ strain data measured with emerging instruments promises to fill an important gap between the episodes of fast strain rates measured by seismic data, and the slow strains measured over relatively long periods of time by InSAR and GPS.

42 ENGINEERING↗

Dataset: "Widespread Drought-driven Declines in Streamflows and Water quality in the Upper Colorado River Basin (1998-2022)"

This data package contains the associated data and scripts for Nagamoto, E., Ombadi, M., Ciulla, F. et al. Widespread drought-driven declines in streamflows and water quality in the Upper Colorado River Basin during 1998-2022. Commun Earth Environ 7, 734 (2026). https://doi.org/10.1038/s43247-026-03890-5. This purpose of this study was to investigate the impact of the 21st century drought on water quantity and quality at catchments throughout the Upper Colorado River Basin (UCRB). We used stream flow, water temperature, specific conductance, air temperature, precipitation, and catchment attribute data for over 200 sites in the UCRB, collected from the National Water Information System using Basin3D (Varadharajan, 2023), GAGESII (Falcone, 2010), and the Google Earth Engine. We identified years of severe drought between 1998 and 2022 using the Standardized Precipitation Evaporation Index (SPEI), then calculated the relative change percentage of the stream flow, water temperature, and specific conductance from drought versus non-drought years. We used the attribute information from GAGESII to investigate what physical traits of catchments are associated streamflow vulnerability (greater relative change) or resilience to drought. We used land cover data from the National Land Cover Database (USGS, 2024) to assess any changes to physical attributes that may not be represented in the static attributes information in GAGESII. To increase data availability, we modeled stream temperature using methods from Willard, 2023. While the study period is water years 1998 to 2022, the raw water quantity and quality data extends to 1950 and the meteorological data extends to 1980. The data and code can be downloaded via the UCRB_drought.zip. Within the zip, the files are organized as follows: - INPUTS: Contains all input data used in UCRB_Drought_Workflow.ipynb - OUTPUTS: Contains all intermediate data created from UCRB_Drought_Workflow.ipynb as well as final products including the calculated Standardized Evapotranspiration Index (SPEI) - climatic_variables: The code used to collect meteorologic data from Google Earth Engine - feature_importance: The code used for the catchment attributes analysis - preprocessing: Code used in UCRB_Drought_Workflow_Preprocessing.ipynb - pyeto: Code used in UCRB_Drought_Workflow_Preprocessing.ipynb - calculations: Code used in UCRB_Drought_Workflow_Impacts.ipynb - plotting: Code used in UCRB_Drought_Workflow_Impacts.ipynb - README.md - UCRB_Drought_Workflow_Preprocessing.ipynb: The code used to prep raw data for the analysis - UCRB_Drought_Workflow_Impact.ipynb: The code which uses the prepped raw data for analysis, and plots all figures - requirements_ucrb-drought_v2.yml: The requirements file to create a virtual environment and Jupyter Lab kernel to run the code The INPUTS folder is organized into the following major directories and sub-directories. The "RDC_WT_SC_RAW" folder contains raw data for streamflow, water temperature, and specific conductance in a ".h5" file. The "NLCD_RAW" folder contains ".csv" files with annual land cover percentages for counties within the UCRB. The "MET_RAW" folder contains a ".csv" file with monthly meteorological data (air temperature and precipitation) for the sites in the UCRB which was obtained from code in the climatic_variables folder. The "GAGESII" folder contains ".csv" files with physical catchment attribute variables for catchments across the country. The "WT_LSTM_data" folder contains ".csv" files with calculated WT (Willard, 2023) and the associated RMSEs. The "Upper_Colorado_River_Basin_Boundary" folder contains geographic data including a shapefile for plotting in the UCRB_Drought_Workflow.ipynb. The "RESERVOIRS_RAW" folder contains ".csv" files for each reservoir in the UCRB with daily reservoir storage. There are also two files in the INPUTS folder that have combined reservoir storage data and reservoir metadata. The OUTPUTS folder is organized into the following major directories and sub-directories. The "RDC_WT_SC_data" folder contains a folder "Water_year" with the associated cleaned data, metadata, and data availability information in ".csv" files, a folder "Median_Relchange" with the relative change comparing drought to non-drought years in ".csv" files, and a folder "Peak95_Min5_Relchange" that has ".csv" files for the relative change in peak (95th %) and minimum (5th %) variables. The "NLCD_data" folder contains the difference in land cover from the beginning to end of the study period and the percentage of the county that is within UCRB bounds can be found in Nagamoto et al (2025)). The "MET_data" folder contains separated monthly air temperature and precipitation data and the calculated PET in ".csv" files. The "SPEI_data" folder contains ".csv" files with calculated SPEI values (one restricted to the study period and the other with information from the entire MET data period). The "Paper_Tables" folder contains two ".csv" files containing site information and data availability and information about the GAGESII trait aggregated categories. The base directory includes the file “flmd.csv” for a list and description of all files and the file “dd.csv” for data dictionaries. Scripts for preprocessing, analysis, and figure generation are located in the associated GitHub repository found at [https://github.com/iNAIADS/drought-impacts/tree/develop/UCRB-drought]. UPDATE 1: Title and code file updated to match submitted manuscript 10-15-2025. UPDATE 2: Code and data files updated to match revised manuscript 3-4-2026. UPDATE 3: Code and data files updated to match revised manuscript 6-7-2026. ** NOTE: DD and FLMD have not been updated yet. UPDATE 4: Added associated Manuscript information and DD and FLMD have been updated. To cite this code, please use the following BibTeX: @misc{nagamoto2025drought, author = {Emily Nagamoto and Fabio Ciulla and Mohammad Ombadi and Jared Willard and Rosemary Carroll and Charuleka Varadharajan}, title = {Dataset: "Widespread Drought-driven Declines in Streamflows and Water quality in the Upper Colorado River Basin (1998-2022)"}, year = {2025}, doi = {10.15485/2551894}, publisher = {ESS-DIVE Repository}, url = {https://data.ess-dive.lbl.gov/datasets/doi:10.15485/2551894} }

54 ENVIRONMENTAL SCIENCES↗

Factors controlling injection-induced rupture of intersecting faults during geological sequestration of CO 2

This study addresses coupled multiphase fluid flow and geomechanics effects on potential fault activation associated with subsurface CO 2 injection around intersecting faults. An enhanced fault-representation model is used to capture geomechanical responses of two intersecting faults with finite length during CO 2 injection. The faults are embedded in a strike-slip stress regime of a caprock-reservoir-basement system with the faults represented by zero-thickness interfaces with adjacent finite-thickness damage zones. A sensitivity analysis is conducted to study the effect of fault permeability, slip-weakening behavior, well location relative to the orientation of faults, and well placement (the number and location of injection wells). Five metrics (pressure, CO 2 plume, shear state on the fault, as well as shear displacement and stress path at selected fault monitoring points) are selected to assess CO 2 migration and reactivation of intersecting faults. The results show that induced ruptures are favored by low permeability faults due to high pressure buildup and by slip-weakening behavior resulting from fault strength reduction. The location of one injection well relative to fault orientation determines the magnitude of changes in effective normal stress and shear stress, affecting the location of induced ruptures. Well placement (two injection wells used in the paper) dominates pressure diffusion around the intersection and tips of faults. This redistributes changes in effective normal stress caused by each injection well, influencing the spatial distribution of ruptures along faults. A larger injection volume induces far-field ruptures that are controlled by stress transfer within the injection layer. The findings presented here can provide valuable insights into engineering operations for a long-term, safe, and reliable geologic CO 2 storage.

Fault permeability↗

Advanced supercritical water-based process concepts for treatment and beneficial reuse of brine in oil/gas production

Produced water generated from oil/gas reservoirs is a cost burden for oil/gas producers, with over 21 billion barrels of the waste generated in 2019. Average disposal costs range from 4.00–8.00·bbl -1 , with up to 20·bbl -1 spent on transportation alone. Dissolved solids content is a contributing factor which limits economical treatment options for this waste stream. This work proposes a novel technique that can handle high salinity waste, employing favorable properties past the critical point of water. In this study, process simulation of two supercritical water desalination (SCWD) scenarios was completed. Zero liquid discharge (ZLD) and brine concentration operating scenarios were compared, weighing the associated economics and benefits for each case. The results were shown to be economically feasible for brines with a high dissolved solids content, ranging from 3.49 to 17.28·bbl -1 in an expanded sensitivity analysis.

42 ENGINEERING↗