Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “USG”

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

Data From: "Warming and snow loss increase reliance on old groundwater in a Colorado River headwater"

This repository contains the data and code associated with the paper titled "Warming and snow loss increase reliance on old groundwater in a Colorado River headwater," published in Nature Geoscience, 2026. This study seeks to answer how various ages of groundwater interact with mountainous streamflow in mountainous headwaters such as the East River. It includes various model-data processing scripts, primarily for ParFlow-CLM analysis of simulated water years 2015-2021, and two numerical warming experiments (+2.5 and +4.0 degrees C), including run scripts, forcing scripts, and post-processing, as well as comparison to observation datasets, detailed below. This data requires the use of R (.r, .rmd), Python (.py), Jupyter Notebook or Jupyter Lab (.ipynb), ParFLOW-CLM, EcoSLIM. Further information on the use of all file formats mentioned below (e.g. .tff. .nc) are provided within the associated scripts and directory where the files are located. Contents & Usage ASO/: ​​Contains the bash and python scripts used to convert airborne snow observatory (ASO) data (ASO, 2023) in various data formats (georeferenced tiff file, NetCDF, UTM, and to latitude/longitude) then regrided to the ParFlow equivalent grid. Output data are in regrid_regll_data.zip and subsequently visualized and analyzed in plot_and_compare.py for Supplementary Figures A14 and A15. The wksht_ASO_comparison.xlsx spreadsheet is used to calculate the data for Supplementary Figure A16. EcoSLIM/: Contains the scripts and input files to run the EcoSLIM particle tracking simulations (/run_scripts) and the post-processing python script (/plot_scripts/eco_agedist_plots.ipynb). Jasechko et al./: Contains the jupyter notebook (Extract_Elevation.ipynb) to determine the outlet elevations of the 260 watersheds used in Jasechko et al. (2016), and the corresponding table, Table_S1_Watersheds_alt.csv. Used to create Supplementary Information Figure A2. PLM_Wells/: Contains the QA/QC-ed groundwater level time series of the PLM-1 and PLM-6 Monitoring Wells from Faybishenko et al. (2023), reformatted to water years used for Supplementary Figures A19 and and A20. ParFlow/: Contains the input files and run scripts to run ParFlow-CLM (/run_scripts), the python and tool command language (Tcl) scripts to create and distribute the ParFlow forcing simulation files (/forcing), and various scripts and intermediary files to analyze the model outputs (/post_process). SQUIRE/: Contains the processing scripts and intermediary files for the Surface QUantitatIve pRecipitation Estimation (SQUIRE) data (Grover, 2023) used to generate Supplementary Figure A18. USGS_Streamflow/: Contains the raw and gap-filled United States Geological Survey streamflow data (U.S. Geological Survey, 2026) used at the Almont station (site number 09112500). Gap-filling is performed in the R script with data from the Taylor station (site number 09110000). (/USGS_09112500_EAST_RIVER_AT_ALMONT_GAP_FILLED/code_almont_streamflow_gap_fill.Rmd). discharge/: Contains the gap-filled discharge data at the Watershed Function SFA East River pumphouse site (Newcomer et al., 2022) used to generate Supplementary Figure A13 and to compute hourly Nash-Sutcliffe model efficiency coefficients (NSE) in Table A4. snotel_and_flux_tower/: Contains the snow telemetry data (U.S. Department of Agriculture, 2024) from the Butte (site ID 380) and Schofield (site ID 737) stations, reformatted by water year, accessed with the snotelr R package. Used to create Supplementary Figure A17. Also contains the flux tower observational data (FluxTower_Pumphouse_ESS-DIVE.ET_only.h.txt) from Ryken et al. (2022) and sap flux transpiration data (MaxB_Transpiration_5Sites.daily_sums.h.txt) from Ryken (2021), used to create Supplementary Figures A22 and A23, respectively. Raw EcoSLIM model outputs are in excess of 24TB, and are stored on National Energy Research Scientific Computing Center (NERSC) and publicly available via the external link provided in the paper.

atmospheric warming↗

A bespoke model of Arctic river basins based on hillslope delineation: Model Archive

This dataset is a model archive of the paper A bespoke model of Arctic river basins based on hillslope delineation (in prep), which introduces a watershed decomposition and parameterization method for large scale permafrost hydrology simulation. With this dataset, this study aims to address the research question: whether a computationally efficient hillslope-based modeling framework can reliably simulate discharge at Arctic river-basin scales. This dataset contains model input and output data for five modeling scenarios at a study site located in the Sagavanirktok River basin. The five modeling scenarios include three modeling cases under temperate conditions using full 3D, decomposed 3D, and decomposed 2D modeling strategies; and two modeling cases under actual Arctic conditions with permafrost using full 3D and decomposed 2D modeling strategies. Simulations were performed using the Advanced Terrestrial Simulator (ATS, v1.6 for three temperate scenarios and v1.5 for two Arctic scenarios), a physics-rich integrated surface–subsurface hydrologic model with cryo-hydrology features. For the three temperate models, simulations were conducted for the period of 10/01/1993 - 09/30/2002; and for the two Arctic models, simulations were conducted for the period of 01/01/1994 - 12/31/2002. To facilitate reproducibility of simulations, all datasets are organized hierarchically. The dataset contains: (1) Mesh files (.exo) for full 3D model, decomposed 3D models, and decomposed 2D models, located in huc/190604020802_gauge15906000/mesh/. Mesh files can be visualized through Paraview or read by Python. (2) Climate forcings (.h5) for full 3D model and decomposed 3D/2D models are located in huc/190604020802_gauge15906000/daymet_onePiece/, and huc/190604020802_gauge15906000/vp_pr_revised_daymet_1980_2006_with_wind/ separately. Accessible by Python. (3) Raw measured gage discharge (.csv) from USGS, located in huc/190604020802_gauge15906000/gaged_basin15906000_discharge_usgs/. Accessible by Python. (4) Delineated subdomain raster (.tif) and shape files (.shp), and the final parameterized results (.npy) for decomposed models, located in huc/190604020802_gauge15906000/data_preprocessed-meshing. Accessible by Python. (5) Temperate models are located in nonpermaf_huc190604020802_gauge15906000/, which includes three cases: decomposed 2D models (inside model_0*-hillslope_*), decomposed 3D models (inside model_1*-subcatchment_*), and full 3D model (inside model_2*-onepiece_*). Two step spin-up results (checkpoint_final.h5) are located in model_*1-*_spinup_steadystate and model_*2-*_spinup_cycle, separately, which are used to initialize real transient models. The input files (.xml) and output results (.dat) of the real transient models are located in model_*3-*_transient/. Especially, for two example hillslope models (ID=-11 and 11), additional h5py files are included in model_03-hillslope_transient/hillslope-11/, model_03-hillslope_transient/hillslope11, model_13-subcatchment_transient/subcatchment-11/, model_13-subcatchment_transient/subcatchment/11, respectively, which are used to plot the saturation figure (Figure 5) in the manuscript. Accessible by Python. (6) Arctic models are located in huc190604020802_gauge15906000/, which includes two cases: decomposed 2D models (inside model_04-hillslope_transient), and full 3D model (inside model_05-onepiece_transient_mannp1_ra). Three step spin-up results (checkpoint_final.h5) are located in model_01-column_freezeup/, model_02-column_spinup/, model_03-hillslope_spinup/, respectively, which are used to initialize real 2D transient hillslope models. The input files (.xml) and output results (.dat) of transient 2D hillslope models are located in model_04-hillslope_transient/. The input files (.xml) and output results (.dat) of the full 3D transient model is located in model_05-onepiece_transient_mannp1_ra/. The full 3D transient model is initialized by model_02-column_spinup/. Accessible by Python. (7) The MOSART routed discharge results (.csv) under Arctic conditions is located in huc190604020802_gauge15906000/MOSART/. Accessible by Python. (8) All Python codes (.py) used to parameterize full 3D model to decomposed 2D models are located in script/. These codes fit with watershed workflow (a watershed delineation tool) v1.4 under the branch gaob/v1.4 from https://github.com/gaobhub/watershed-workflow.git.

EARTH SCIENCE > CRYOSPHERE↗

AmeriFlux FLUXNET-1F US-Tw1 Twitchell Wetland West Pond

This is the AmeriFlux Management Project (AMP) created FLUXNET-1F version of the carbon flux data for the site US-Tw1 Twitchell Wetland West Pond. This is the FLUXNET version of the carbon flux data for the site US-Tw1 Twitchell Wetland West Pond produced by applying the standard ONEFlux (1F) software. Site Description - The Twitchell Wetland site is a 7.4-acre restored wetland on Twitchell Island, that is managed by the California Department of Water Resources (DWR) and the U.S. Geological Survey (USGS). In the fall of 1997, the site was permanently flooded to a depth of approximately 25 cm. The wetland was almost completely covered by cattails and tules by the third growing season. A flux tower equipped to analyze energy, H2O, CO2, and CH4 fluxes was installed on May 17, 2012.

Valach, Alex↗

AmeriFlux FLUXNET-1F US-Tw5 East Pond Wetland

This is the AmeriFlux Management Project (AMP) created FLUXNET-1F version of the carbon flux data for the site US-Tw5 East Pond Wetland. This is the FLUXNET version of the carbon flux data for the site US-Tw5 East Pond Wetland produced by applying the standard ONEFlux (1F) software. Site Description - The Twitchell Wetland site is a 6.5 acre restored wetland on Twitchell Island, that is managed by the California Department of Water Resources (DWR) and the U.S. Geological Survey (USGS). In the fall of 1997, the site was permanently flooded to a depth of approximately 55 cm. The wetland remained fairly unvegetated in patches increasing in size towards the east. The site underwent a major disturbance in 2013 when the vegetation was removed to seed a nearby restored wetland. A flux tower equipped to analyze energy, H2O, CO2, and CH4 fluxes was installed on April 17, 2018.

Valach, Alex↗

Crowdsourcing Felt Reports Using the MyShake Smartphone App

MyShake is a free citizen science smartphone app that provides a range of features related to earthquakes. Features available globally include rapid postearthquake notifications, live maps of earthquake damage as reported by MyShake users, safety tips, and various educational features. The app also uses the accelerometer in the mobile device to detect earthquake shaking, and to record and submit waveforms to a central archive. In addition, MyShake delivers earthquake early warning alerts in California, Oregon, and Washington. Here, in this study, we compare the felt shaking reports provided by MyShake users in California with the U.S. Geological Survey’s (USGSs) “Did You Feel It?” intensity reports. The MyShake app simply asks, “What strength of shaking did you feel?” and users report on a five-level scale. When the MyShake reports are averaged in spatial or time bins, we find strong correlation with the Modified Mercalli Intensity scale values reported by the USGS based on the DYFI surveys. The MyShake felt reports can therefore contribute to the creation of shaking intensity maps.

58 GEOSCIENCES↗

National Critical Minerals Data Dashboard

The Critical Minerals (CM) Data Dashboard provides data-driven insights to help unlock the potential of conventional and unconventional CM resources in the United States. The dashboard focuses on 12 CM essential for producing zero-emission transportation and clean power technology. More information about these CM and potential resources is available in our CM Story Map. Dashboard data sets can be filtered geographically using the options provided along the top of the dashboard. Data sets can also be filtered by various attributes using the options in the expandable side panel on the left side of the screen. The data sets in this dashboard were compiled from multiple sources including the USGS, USEPA, OSMRE, and DHS., Users are encouraged to provide feedback/suggestions on the function and content of the dashboard, as well as any issues that may arise.

Conventional Resources↗

Constituent Data Replacement Tool

The purpose of this tool is to estimate key parameters that may be missing in public wastewater composition datasets. The tool can be applied to develop complete treatment and critical mineral extraction profiles for leachate, produced water and other aqueous waste streams. The tool applies machine learning algorithms to replace missing data in a user’s water data set that are adjusted based on user preferences for options including algorithm type, number of features, and classification variables. The tool can use the user’s data alone or combine user data with the NEWTS USGS Produced Water Database for more robust training. This research was funded by the U.S. Department of Energy’s Office Fossil Energy and Carbon Management (FECM) through National Energy Technology Laboratory’s ongoing research under the Water Management for Power System Field Work Proposal, DE-FECM 1022428 and Critical Minerals Field Work Proposal, DE-FECM 1022420.

Aqueous Chemistry↗

Datasets and U-Net Model for "A Deep Learning Based Framework to Identify Undocumented Orphaned Oil and Gas Wells from Historical Maps: a Case Study for California and Oklahoma"

This dataset has results and the model associated with the publication Ciulla et al., (2024). It contains a U-Net semantic segmentation model (unet_model.h5) and associated code implemented in tensorflow 2.0 for the model training and identification of oil and gas well symbols in USGS historical topographic maps (HTMC). Given a quadrangle map (7.5 minutes), downloadable at this url: https://ngmdb.usgs.gov/topoview/, and a list of coordinates of the documented wells present in the area, the model returns the coordinates of oil and gas symbols in the HTMC maps. For reproducibility of our workflow, we provide a sample map in California and the documented well locations for the entire State of California (CalGEM_AllWells_20231128.csv) downloaded from https://www.conservation.ca.gov/calgem/maps/Pages/GISMapping2.aspx. Additionally, the locations of 1,301 potential undocumented orphaned wells identified using our deep learning framework or the counties of Los Angeles and Kern in California, and Osage and Oklahoma in Oklahoma are provided in the file found_potential_UOWs.zip. The results of the visual inspection of satellite imagery in Osage County is in the file visible_potential_UOWs.zip. The dataset also includes a custom tool to validate the detected symbols in the HTMC maps (vetting_tool.py). More details about the methodology can be found in the associated paper: Ciulla, F., Santos, A., Jordan, P., Kneafsey, T., Biraud, S.C., and Varadharajan, C. (2024) A Deep Learning Based Framework to Identify Undocumented Orphaned Oil and Gas Wells from Historical Maps: a Case Study for California and Oklahoma. Accepted for publication in Environmental Science and Technology. The geographical coordinates provided correspond to the locations of potential undocumented orphaned oil and gas wells (UOWs) extracted from historical maps. The actual presence of wells need to be confirmed with on-the-ground investigations. For your safety, do not attempt to visit or investigate these sites without appropriate safety training, proper equipment, and authorization from local authorities. Approaching these well sites without proper personal protective equipment (PPE) may pose significant health and safety risks. Oil and gas wells can emit hazardous gasses including methane, which is flammable, odorless and colorless, as well as hydrogen sulfide, which can be fatal even at low concentrations. Additionally, there may be unstable ground near the wellhead that may collapse around the wellbore. This dataset was prepared as an account of work sponsored by the United States Government. While this document is believed to contain correct information, neither the United States Government nor any agency thereof, nor the Regents of the University of California, nor any of their employees, makes any warranty, express or implied, or assumes any legal responsibility for the accuracy, completeness, or usefulness of any information, apparatus, product, or process disclosed, or represents that its use would not infringe privately owned rights. Reference herein to any specific commercial product, process, or service by its trade name, trademark, manufacturer, or otherwise, does not necessarily constitute or imply its endorsement, recommendation, or favoring by the United States Government or any agency thereof, or the Regents of the University of California. The views and opinions of authors expressed herein do not necessarily state or reflect those of the United States Government or any agency thereof or the Regents of the University of California.

Artificial Intelligence↗

Evolution and Trends of Industrial Control System Cyber Incidents since 2017

The industrial control systems (ICSs) that manage our critical infrastructure are increasingly converging with corporate networks and the Internet as technology and businesses prioritize digital connectivity. These connections make them more vulnerable and available to malicious cyber actors who traditionally targeted the companies’ more public-facing information technology (IT) networks. This paper will review select publicly reported cyber incidents to highlight the continued and growing threat to ICS devices and operational technology (OT) environments. It will summarize the incident and when available, will provide information on the cyber actors, the vulnerabilities they exploited, and any publications the U.S. Government (USG) provided in response. Data belonging to the Department of Homeland Security (DHS) will be used to highlight quantitative trends concerning ICS incidents. This paper builds on “History of Industrial Control System Cyber Incidents” (Hemsley & Fisher 2018), a paper that highlighted select noteworthy threats and incidents to ICS systems up to 2017. This paper will similarly review select incidents occurring after the last previously reviewed incident, Triton/HatMan, December 2017, and will note ICS incident trends including IT/OT convergence and advances in cyber-threat actors’ capabilities in observed in the examined incidents.

99 GENERAL AND MISCELLANEOUS↗

Southwest Wind Power R&D Test Site Development at the National Wind Technology Center (CRADA CRD-12-00472 Final Report)

Southwest Windpower, Inc. (SWWP) has been designing and distributing small wind turbines for more than 22 years and is the recognized global leader in the design, manufacturing and distribution of small wind systems (400-3000 watts). The company has been a pioneer in the development of wind technology and has built and shipped more than 170,000 wind turbines to over 120 countries worldwide. Headquartered in Flagstaff AZ, SWWP has sales representatives in over 88 countries. Applications for SWWP systems include residential homes, commercial properties, micro grids, remote cabins, telecom transmitters, offshore platforms, water pumping and sailboats. In addition, U.S. Department of Defense uses SWWP’s products in Forward Operating Bases (FOB’s), USGS uses them for remote monitoring of glacier movements and the CIA uses them to provide power seismic monitoring. Today, SWWP plans development of new technologies to address distributed energy market needs and become even more aggressive in international markets. This strategy includes extensive Research, Development, Demonstration & Deployment (RDD&D) activities to improve existing production line into more simple, reliable, and cost-effective wind turbine systems. Important part of such activities is the field testing of the new systems, which should provide: (i) validation of the reliable operation of hardware and software; (ii) validation of numerical models of the system and its components; (iii) field comparison of various control and optimization strategies; (iv) customized reliability testing; (v) customized monitoring and analysis testing; and (vi) IEC certification testing (vii) filed demonstration of system performance and reliability. Such field testing facility should encompass multiple towers and appropriate data acquisition systems. SWWP is exploring to move part of the engineering department to Colorado, and execute its RDD&D strategies in this new location. Partnership with NREL will provide great insights into SWWP’s RDD&D processes and more efficient commercialization by utilizing NWTC as the field testing site, and working close with NREL personnel. NWTC is ideally suited to strain the design. NREL is very experienced with customized and IEC testing. During this project, Southwest Windpower aims to improve the value proposition of distributed renewable energy, and make it a competitive choice in the energy markets. As of today, more then 170,000 small wind turbines, manufactured in the U.S., are installed in more then 180 countries. With the improved value proposition, Southwest Windpower wants to create and address new market needs, substantially increase its sales domestically and internationally, which should have positive impact on the creation of American jobs, competitiveness of U.S. economy, export or U.S. products, and the use of renewable energy.

17 WIND ENERGY↗

Optimized Microwave Digestion and Quantification Procedure for Boron Carbide samples

A microwave-assisted HNO3-H2SO4-HF digestion system was explored/optimized for the total dissolution of Boron Carbide samples followed by multi-elemental determination using ICP-MS and ICP OES, in order to improve the methods used in a previous try using a microwave-assisted H2O2-HNO3-HF digestion system. The samples were microwave digested to accomplish complete dissolution needed to perform quantitative analysis of their metal content. Two Boron Carbide samples were provided to us: JM10 and JM11. The Boron Carbide samples were completely dissolved at the end of the optimized microwave digestion procedure that was developed. The digested samples were analyzed for metals mostly by ICP MS; few elements were analyzed by ICP-OES. Microwave digestion reactions were performed in a Titan MPS instrument (Perkin Elmer, USA). Titan MPS instrument is equipped with Temperature and Pressure regulations and controlled by software. Each microwave digestion batch contained the samples of interest along with a reagent blank and one certified standard (1632e) in order to (i) compensate for the contaminations present in the reagents, (ii) ensure that the digestion was complete and (iii) that there was a good recovery of all the constituents. The certified standard used (1632e) was provided with USGS certificates of analysis.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

CO2 Capture Strategies via Mineralization with Industrial Waste Brines

Large coal-fired power plants (>500 MW) account for 30% of global CO2 emissions, and long-term management of this CO2 to is urgently needed mitigate global temperature increases. Sequestration of CO2 within stable mineral carbonates (e.g., CaCO3) represents an attractive emission reduction strategy because it offers a leakage-free alternative to geological storage of CO2 in an environmentally friendly form. We have previously described a mineralization process in which divalent cations are sourced from various waste streams (e.g., produced water and brackish water) and alkalinity is induced via regenerable ion-exchange materials (Bustillos et. al. Frontiers in Energy Research. 2020, 8, 352). In our process, aqueous carbonate-bearing streams with pH > 8 are produced by contacting fresh water and carbon dioxide with various ion-exchange materials (e.g., Na form zeolites or ion exchange resins). These streams are mixed with produced water containing varying concentrations (~0.01 – 1.0 M) of Ca2+ leading to the precipitation of solid calcium carbonate (PCC). This process has the advantages of using regenerable solids in a simple and continuous process to increase the pH of water by ion exchange instead of relying on the consumption of costly and unsustainable sources of alkalinity (e.g., sodium hydroxide). While once-through column experiments showed the above benefits, the same were yet to established in a steady-state process with recycle streams. In this work, we set up a process simulation to quantify the energy requirements and CO2 emissions associated with the process and seek optimal produced water compositions and CO2 concentrations (5 – 20 vol%). The process simulation was set up in ASPEN Plus using eRNTL as the thermodynamic property method and sequential modular strategy. Ion exchange alkaline solution was simulated using sodium hydroxide and validated against the experimental data obtained from once-through kinetic experiments. Nanofiltration and reverse osmosis membrane steps were also implemented for the separation of divalent cations and production of fresh water and a regeneration stream following mineralization. Sensitivity analysis was carried out using a range of produced water compositions (0.01 – 1.0 M Ca2+, 0.001 – 0.15 M Mg2+, 0.5 – 3.5 M Na+ and 0.0004 – 0.002 M Fe2+) according to the United States Geological Survey (USGS) database. Calcium carbonate yields increased with increasing CO2 concentrations and were maximized using produced water compositions with larger Ca2+ concentrations. Maximum calcium carbonate yields produced at 5 vol%, 12 vol% and 20 vol% CO2 were 2.3 mmol/L, 5.5 mmol/L, and 9.3 mmol/L, respectively, with the formation of brucite (a magnesium hydroxide phase, Mg(OH)2) and goethite (an iron hydroxide phase, FeOOH) as the primary contaminant phases (99% calcite, 0.6% brucite, 0.4% goethite), which agree with phases detected by XRD experimentally. These results indicate high purity calcium carbonate can be precipitated using industrial waste streams. Consequentially, energy consumption and net CO2 emissions were minimized where precipitated calcium carbonate was maximized for all produced water compositions and CO2 concentrations. Minimum energy consumptions were 0.21 kWh/ton CO2 processed, with 98% of the energy input required coming from the membrane filtration steps. Produced water compositions with large Na+ concentrations (> 0.5 M) were effective at reducing energy consumptions due to faster regeneration time of ion exchange materials. Additionally, calculated net CO2 emissions were negative for the process and ranged from -0.02 kg/ton CO2 to -0.15 kg/ton CO2 processed, indicating a low emission process. We will also present techno-economic assessment showing the economic benefits of the current process as an alternative to the addition of stoichiometric bases to induce alkalinity for the precipitation of CaCO3.

Simonetti, Dante↗

Middle Kittanning Coal Waste and Underclay as an Alternative Rare Earth Elements Feedstock

In order to secure domestic sources of rare earth elements (REE) from coal related materials, there must be validation of representative feedstocks. Actively producing coal mines that target the Middle Kittanning coal seam in the Appalachian Basin were compiled. These mines were cross-referenced with publicly available geochemical data such as the U.S. Geological Survey (USGS) Earth Mapping Resources Initiative geochemical data along with samples evaluated and characterized by the National Energy Technology Laboratory (NETL). This work evaluated the extent of elevated Middle Kittanning underclay concentrations of REE in comparison to other underclay formations. Therefore, further up-scaling of research associated with the separation and extraction of REE and other critical minerals can be beneficial to utilizing domestic REE supplies for various technology sectors, to include energy, biomedical, and defense. Numerous current active mines targeting the Middle Kittanning coal seam represent a geographically significant opportunity for shared feedstocks and collaborations to further understand the role of Middle Kittanning underclay as a critical mineral feedstock.

01 COAL, LIGNITE, AND PEAT↗

Financial Analysis of Experimental Releases Conducted at Glen Canyon Dam during Water Year 2021

This report was prepared by Argonne National Laboratory (Argonne) in support of a financial analysis of two Glen Canyon Dam (GCD) flow experiments that were conducted in Water Year 2021. The first experiment, called “spring disturbance flow”, was conducted from March 15 to 26, 2021. This disturbance flow experiment is designed to simulate a spring-timed runoff event so scientists can study its effects on the Colorado River ecosystem. The second experiment was a reduced and steady discharge from May 29 to June 4, 2021. This reduced steady flow experiment was conducted for the United States Geological Survey (USGS) to collect high-resolution aerial imagery over Grand Canyon National Park. This analysis was funded by the Colorado River Storage Project (CRSP) Office of the U.S. Department of Energy’s Western Area Power Administration (WAPA). CRSP markets electricity produced by hydroelectric facilities collectively known as the Salt Lake City Area Integrated Projects including dams equipped for power generation on the Colorado, Green, Gunnison, and Rio Grande Rivers and on Plateau Creek in the states of Arizona, Colorado, New Mexico, Utah, and Wyoming. Staff members in Argonne’s Energy Systems Division prepared this technical memorandum with assistance from WAPA’s CRSP and Energy Marketing and Management Offices (EMMO).

58 GEOSCIENCES↗

Development of a molecularly informed biogeochemical framework for reactive transport modeling of subsurface carbon inventories, transformations and fluxes (Final Report)

The overall objective of the project was to combine new molecular-level characterization strategies with soil carbon flux measurements to develop and evaluate model representations of subsurface carbon cycling. We expanded on existing studies in the East River watershed, Colorado, in collaboration with the Berkeley Lab Watershed Function Scientific Focus Area (SFA), SLAC Groundwater Quality SFA, and Rocky Mountain Biological Laboratory (RMBL) to develop an elevation and vegetation gradient that is now the subject of long-term monitoring by the USGS. To achieve the overall objective, we combined field studies of soil respiration with laboratory analyses, ranging from spectroscopy to incubation studies. This combination of techniques enabled us to develop a new understanding of the drivers of high-elevation soil respiration. In the process, we developed three new modeling approaches that improve our ability to conceptualize and ultimately to represent soil respiration in numerical models. The first approach is a plot-scale transient inverse model that can be used to determine in situ CO 2 production rates from measured concentration profiles and surface fluxes. Application of this method revealed the importance of plant phenology and deep CO 2 production in moderating CO 2 fluxes to the atmosphere. The second modeling approach is a molecular-scale tool that enables spectroscopic and elemental data for carbon speciation to be transformed into functional group abundances, or the ‘SOC-fga model’. This method uniquely enables carbon speciation to be tracked within a reactive transport framework to partition carbon among different pathways and storage zones within the soil. The third modeling tool builds strongly on the previous approaches and captures the microbial processes driving heterotrophic respiration. This ‘dormancy model’ allows the native soil microbial population to respond transiently to the presence or absence of water in order to catalyze carbon respiration. This approach was also compared to the simpler and more widely used first-order model using two experimental datasets with different temporal and spatial resolutions. Our results illustrate that the simpler first order model provides a robust and efficient representation of deep (>1 meter) soil respiration, but that shallow soils, where most respiration occurs, require explicit representation of moisture-dependent activation and dormancy rates. By assessing soil organic carbon turnover at multiple scales, we see a complex array of controls emerge. Atte scale of a hillslope, spatial heterogeneity in soil respiration rates dominates and is uncorrelated with instantaneous soil moisture and plant community. At the profile scale, the balance between plant inputs and water availability is the dominant control. At the microbial to molecular scale, physiological processes associated with carbon use and carbon speciation are important controls. These scale-dependent controls emphasize the need for new modeling approaches that examine their interactions and hierarchies.

54 ENVIRONMENTAL SCIENCES↗

Refining Principal Stress Measurements in Reservoir Underburden in Regions of Induced Seismicity through Seismological Tools, Laboratory Experiments - Final Technical Report

This project developed methodologies to measure the in-situ principal stress in the deep subsurface through use of multiple independent, but complementary, seismic methods, laboratory verification, and development of theoretical frameworks. By leveraging existing regional and local datasets we developed, tested, and refined a set of diagnostic tools for determining the in-situ stress state with reduced uncertainty at and below reservoir depths (1.5-6 km). A set of novel tools was produced that are scale independent, such that their utility is equivalent on regional, field scale, and near borehole monitoring of principal stresses in reservoir underburden for carbon storage projects. During a 4-year Department of Energy (DOE) and Southern Company funded project, carried out by the Electric Power Research Institute (EPRI), Lawrence Livermore National Laboratory (LLNL), the University of Oklahoma (OU), and the U.S. Geological Survey (USGS), the project team developed methodologies to measure the far-field in-situ principal stress in the deep subsurface, leveraging induced seismicity data from waste-water disposal projects. These methodologies consisted in the use of well-established and technically advanced seismic processing methods, such as virtual seismometer method-moment tensor (VSM-MT) and shear wave splitting (SWS), that are adept at recovering the stress orientation and certain components of the stress tensor. These methods were applied to robust seismicity catalogs created with matched filter techniques near sites of active fluid disposal—a proxy for carbon storage sites where such datasets are more limited. Estimates of the stress orientation made with seismic processing tools were considered along with laboratory acoustic emission experiments conducted on rock samples from the region of interest. Stress orientations in the studied region do not vary significantly across distances of ~100 km, nor are they found to rotate through time as a consequence of local wastewater disposal, as previously speculated. Finally, the project team investigated the trade-offs among the different seismic methods and evaluated the range of uncertainty that is generated with these methodologies, which led to a practical use and refinement of the VSM-MT technique when it is applied to field datasets. Understanding the trade-offs between these different methods highlighted the potential benefits of improved quantification of uncertainties on stress field estimations.

58 GEOSCIENCES↗

Rock Valley Dense Gravity Acquisition

Characterizing the shallow structure of the Rock Valley region of the Nevada National Security Site is a critical component of the Rock Valley Direct Comparison project. Geophysical data of the region is needed for operational decisions, to constrain geologic models used for simulation, and to facilitate the analysis of future explosive source data. Local measurements of gravity are a key piece of geophysical information that helps to resolve the underlying geologic composition, fault structure, and density characteristics, yet, in the Rock Valley region these measurements are sparse on the scale of the testbed. In this report, we present the details of a recent gravity data acquisition survey designed to collect a dense dataset in the region of interest that complements the existing gravity work but greatly enhances our resolution. This dataset will be integrated with a complementary Los Alamos National Laboratory gravity collection and combined with the existing seismic data in a joint inversion. These measurements were conducted over two weeks with a portable gravimeter and high-resolution GPS and include repeat measurements at a USGS base station as well as reoccupation of gravity sites in the regional dataset. This collection of over 100 new dense gravity measurements will facilitate refinement of the existing Geologic Framework Model and directly complement newly acquired dense seismic data, ultimately improving the project’s ability to investigate the direct comparison of shallow earthquake and explosive sources.

58 GEOSCIENCES↗

EQSIM and RAJA: Enabling Exascale Predictions of Earthquake Effects on Critical Infrastructure

Nearly 120 years ago, the great “San Francisco” earthquake of 1906 provided a stark and sobering view of the havoc that can be caused by the sudden and violent movement of Earth’s tectonic plates. According to USGS, the rupture along the San Andreas fault extended 296 miles (447 kilometers) and shook so violently that the motion could be felt as far north as Oregon and east into Nevada. The estimated 7.9-magnitude quake and subsequent fires decimated the major metropolis and surrounding areas: buildings turned to ruins, hundreds of thousands of people left homeless, and a death toll exceeding 3,000. Today, as evidenced by the catastrophic 7.8-magnitude earthquake that struck Turkey in February of 2023, large earthquakes still present a significant potential danger to life and economic security as researchers work to develop ways to better understand earthquake phenomena and quantify associated risks.

42 ENGINEERING↗