Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “NEST”

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

Improved Representation of Horizontal Variability and Turbulence in Mesoscale Simulations of an Extended Cold-Air Pool Event

Abstract Cold-air pools (CAPs), or stable atmospheric boundary layers that form within topographic basins, are associated with poor air quality, hazardous weather, and low wind energy output. Accurate prediction of CAP dynamics presents a challenge for mesoscale forecast models in part because CAPs occur in regions of complex terrain, where traditional turbulence parameterizations may not be appropriate. This study examines the effects of the planetary boundary layer (PBL) scheme and horizontal diffusion treatment on CAP prediction in the Weather Research and Forecasting (WRF) Model. Model runs with a one-dimensional (1D) PBL scheme and Smagorinsky-like horizontal diffusion are compared with runs that use a new three-dimensional (3D) PBL scheme to calculate turbulent fluxes. Simulations are completed in a nested configuration with 3-km/750-m horizontal grid spacing over a 10-day case study in the Columbia River basin, and results are compared with observations from the Second Wind Forecast Improvement Project. Using event-averaged error metrics, potential temperature and wind speed errors are shown to decrease both with increased horizontal grid resolution and with improved treatment of horizontal diffusion over steep terrain. The 3D PBL scheme further reduces errors relative to a standard 1D PBL approach. Error reduction is accentuated during CAP erosion, when turbulent mixing plays a more dominant role in the dynamics. Last, the 3D PBL scheme is shown to reduce near-surface overestimates of turbulence kinetic energy during the CAP event. The sensitivity of turbulence predictions to the master length-scale formulation in the 3D PBL parameterization is also explored. Significance Statement In this article, we demonstrate how a new framework for modeling atmospheric turbulence improves cold pool predictions, using a case study from January 2017 in the Columbia River basin (U.S. Pacific Northwest). Cold pools are regions of cold, stagnant air that form within valleys or basins, and improved forecasts could help to mitigate the risks they pose to air quality, transportation, and wind energy production. For the chosen case study, our tests show a reduction in temperature and wind speed errors by up to a factor of 2–3 relative to standard model options. These results strongly motivate continued development of the framework as well as its application to other complex weather events.

17 WIND ENERGY↗
DOE OSTI2023

Projecting Future Energy Production from Operating Wind Farms in North America. Part I: Dynamical Downscaling

Abstract New simulations at 12-km grid spacing with the Weather and Research Forecasting (WRF) Model nested in the MPI Earth System Model (ESM) are used to quantify possible changes in wind power generation potential as a result of global warming. Annual capacity factors (CF; measures of electrical power production) computed by applying a power curve to hourly wind speeds at wind turbine hub height from this simulation are also used to illustrate the pitfalls in seeking to infer changes in wind power generation directly from low-spatial-resolution and time-averaged ESM output. WRF-derived CF are evaluated using observed daily CF from operating wind farms. The spatial correlation coefficient between modeled and observed mean CF is 0.65, and the root-mean-square error is 5.4 percentage points. Output from the MPI-WRF Model chain also captures some of the seasonal variability and the probability distribution of daily CF at operating wind farms. Projections of mean annual CF (CF A ) indicate no change to 2050 in the southern Great Plains and Northeast. Interannual variability of CF A increases in the Midwest, and CF A declines by up to 2 percentage points in the northern Great Plains. The probability of wind droughts (extended periods with anomalously low production) and wind bonus periods (high production) remains unchanged over most of the eastern United States. The probability of wind bonus periods exhibits some evidence of higher values over the Midwest in the 2040s, whereas the converse is true over the northern Great Plains. Significance Statement Wind energy is playing an increasingly important role in low-carbon-emission electricity generation. It is a “weather dependent” renewable energy source, and thus changes in the global atmosphere may cause changes in regional wind power production (PP) potential. We use PP data from operating wind farms to demonstrate that regional simulations exhibit skill in capturing actual power production. Projections to the middle of this century indicate that over most of North America east of the Rocky Mountains annual expected PP is largely unchanged, as is the probability of extended periods of anomalously high or low production. Any small declines in annual PP are of much smaller magnitude than changes due to technological innovation over the last two decades.

Meteorology & Atmospheric Sciences↗
DOE OSTI2022

Effects of Horizontal Resolution, Domain Size, Boundary Conditions, and Surface Heterogeneity on Coarse LES of a Convective Boundary Layer

Atmospheric properties in a convective boundary layer vary over a wide range of spatial scales and are commonly studied using large-eddy simulations (LES) in various configurations. We examine how the boundary layer depth and distribution of variability across scales are affected by LES grid spacing, domain size, inhomogeneity of surface properties, and external forcing. Two different setups of the Weather Research and Forecasting (WRF) model are analyzed. A semi-idealized configuration uses a periodic domain, flat surface, prescribed homogeneous surface heat fluxes, and horizontally uniform profiles of large-scale advective tendencies. A nested LES setup employs a larger domain and realistic initial and boundary conditions, including an interactive land surface model with representative topography and vegetation and soil types. Subdomains of identical size are analyzed for all simulations. Characteristic structure sizes are quantified using the variability scales L 50 and L 95 , defined such that features smaller than that contain 50% and 95% of the total variance, respectively. Progressive increase in L 50 from vertical velocity to temperature and moisture structures is systematically reproduced in all simulation configurations. This dependence of L 50 on the considered variable complicates the development of scale-aware parameterizations for models with grid spacing in the “terra incognita”. In simulations using a larger domain with heterogeneous surface properties, the development of internal mesoscale patterns significantly affects variance distributions inside analyzed subdomains. Sizes of boundary layer structures also strongly depend on the LES grid spacing and, in case of heterogeneous surface and topography, on location of the subdomain inside a larger computational domain.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Evaluation of Turbulence and Dispersion in Multiscale Atmospheric Simulations over Complex Urban Terrain during the Joint Urban 2003 Field Campaign

Abstract This paper evaluates the representation of turbulence and its effect on transport and dispersion within multiscale and microscale-only simulations in an urban environment. These simulations, run using the Weather Research and Forecasting Model with the addition of an immersed boundary method, predict transport and mixing during a controlled tracer release from the Joint Urban 2003 field campaign in Oklahoma City, Oklahoma. This work extends the results of a recent study through analysis of turbulence kinetic energy and turbulence spectra and their role in accurately simulating wind speed, direction, and tracer concentration. The significance and role of surface heat fluxes and use of the cell perturbation method in the numerical simulation setup are also examined. Our previous study detailed the model development necessary for our multiscale simulations, examined model skill at predicting wind speeds and tracer concentrations, and demonstrated that dynamic downscaling from mesoscale to microscale through a sequence of nested simulations can improve predictions of transport and dispersion relative to a microscale-only simulation forced by idealized meteorology. Here, predictions are compared with observations to assess qualitative agreement and statistical model skill at predicting wind speed, wind direction, tracer concentration, and turbulent kinetic energy at locations throughout the city. We also investigate the scale distribution of turbulence and the associated impact on model skill, particularly for predictions of transport and dispersion. Our results show that downscaled large-scale turbulence, which is unique to the multiscale simulations, significantly improves predictions of tracer concentrations in this complex urban environment. Significance Statement Simulations of atmospheric transport and mixing in urban environments have many applications, including pollution modeling for urban planning or informing emergency response following a hazardous release. These applications include phenomena with spatial scales spanning from millimeters to kilometers. Most simulations resolve flow only within the urban area of interest, omitting larger scales of turbulence and regional influences. This study examines a method that resolves both the small and large-scale flow features. We evaluate simulation accuracy by comparing predictions with observations from an experiment involving the release of a tracer gas in Oklahoma City, Oklahoma, with emphasis on correctly modeling turbulent fluctuations. Our results demonstrate the importance of resolving large-scale flow features when predicting transport and dispersion in urban environments.

42 ENGINEERING↗
DOE OSTI2024

Refining HPCToolkit for application performance analysis at exascale

As part of the US Department of Energy’s Exascale Computing Project (ECP), Rice University has been refining its HPCToolkit performance tools to better support measurement and analysis of applications executing on exascale supercomputers. To efficiently collect performance measurements of GPU-accelerated applications, HPCToolkit employs novel non-blocking data structures to communicate performance measurements between tool threads and application threads. To attribute performance information in detail to source lines, loop nests, and inlined call chains, HPCToolkit performs parallel analysis of large CPU and GPU binaries involved in the execution of an exascale application to rapidly recover mappings between machine instructions and source code. To analyze terabytes of performance measurements gathered during executions at exascale, HPCToolkit employs distributed-memory parallelism, multithreading, sparse data structures, and out-of-core streaming analysis algorithms. To support interactive exploration of profiles up to terabytes in size, HPCToolkit’s hpcviewer graphical user interface uses out-of-core methods to visualize performance data. The result of these efforts is that HPCToolkit now supports collection, analysis, and presentation of profiles and traces of GPU-accelerated applications at exascale. These improvements have enabled HPCToolkit to efficiently measure, analyze and explore terabytes of performance data for executions using as many as 64K MPI ranks and 64K GPU tiles on ORNL’s Frontier supercomputer. HPCToolkit’s support for measurement and analysis of GPU-accelerated applications has been employed to study a collection of open-science applications developed as part of ECP. This paper reports on these experiences, which provided insight into opportunities for tuning applications, strengths and weaknesses of HPCToolkit itself, as well as unexpected behaviors in executions at exascale.

Adhianto, Laksono↗
DOE OSTI2024

Emergence of two distinct phase transitions in monolayer CoSe2 on graphene

Abstract Dimensional modifications play a crucial role in various applications, especially in the context of device miniaturization, giving rise to novel quantum phenomena. The many-body dynamics induced by dimensional modifications, including electron-electron, electron-phonon, electron-magnon and electron-plasmon coupling, are known to significantly affect the atomic and electronic properties of the materials. By reducing the dimensionality of orthorhombic CoSe 2 and forming heterostructure with bilayer graphene using molecular beam epitaxy, we unveil the emergence of two types of phase transitions through angle-resolved photoemission spectroscopy and scanning tunneling microscopy measurements. We disclose that the 2 × 1 superstructure is associated with charge density wave induced by Fermi surface nesting, characterized by a transition temperature of 340 K. Additionally, another phase transition at temperature of 160 K based on temperature dependent gap evolution are observed with renormalized electronic structure induced by electron-boson coupling. These discoveries of the electronic and atomic modifications, influenced by electron-electron and electron-boson interactions, underscore that many-body physics play significant roles in understanding low-dimensional properties of non-van der Waals Co-chalcogenides and related heterostructures. Graphical Abstract

36 MATERIALS SCIENCE↗
DOE OSTI2023

Incubation temperature alters stripe formation and head colouration in American alligator hatchlings and is unaffected by estradiol-induced sex reversal

ABSTRACT Considerations of the impact climate change has on reptiles are typically focused on habitat change or loss, range shifts and skewed sex ratios in species with temperature-dependent sex determination. Here, we show that incubation temperature alters stripe number and head colouration of hatchling American alligators (Alligator mississippiensis). Animals incubated at higher temperatures (33.5°C) had, on average, one more stripe than those at lower temperatures (29.5°C), and also had significantly lighter heads. These patterns were not affected by estradiol-induced sex reversal, suggesting independence from hatchling sex. Therefore, increases in nest temperatures as a result of climate change have the potential to alter pigmentation patterning, which may have implications for offspring fitness.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2026

Dynamic Transmission Line Switching Amid Wildfire-Prone Weather Under Decision-Dependent Uncertainty

During dry and windy seasons, environmental conditions significantly increase the risk of wildfires, exposing power grids to disruptions caused by transmission line failures. Wildfire propagation exacerbates grid vulnerability, potentially leading to prolonged power outages. To address this challenge, we propose a multistage optimization model that dynamically adjusts transmission grid topology in response to wildfire propagation, aiming to develop an optimal response policy. By accounting for decision-dependent uncertainty, where line survival probabilities depend on usage, we employ distributionally robust optimization to model uncertainty in line survival distributions. We adapt the stochastic nested decomposition algorithm and derive a deterministic upper bound for its finite convergence. To enhance computational efficiency, we exploit the Lagrangian dual problem structure for a faster generation of Lagrangian cuts. Using realistic data from the California transmission grid, we demonstrate the superior performance of dynamic response policies against two-stage alternatives through a comprehensive case study. In addition, after solving the multistage formulation, we construct easy-to-implement policies that significantly reduce computational burden while maintaining good performance in real-time deployment. History: Accepted by Russell Bent, Area Editor for Network Optimization: Algorithms and Applications. Funding: This work was supported by the U.S. Department of Energy, Office of Electricity [Grant DE-AC02-05CH11231]. The work of R. Jiang was supported in part by the U.S. National Science Foundation, Division of Electrical, Communications and Cyber Systems [Grant ECCS-1845980] and the U.S. Air Force Office of Scientific Research [Grant FA9550-23-1-0323]. Supplemental Material: The software that supports the findings of this study is available within the paper and its Supplemental Information ( https://pubsonline.informs.org/doi/suppl/10.1287/ijoc.2025.1210 ) as well as from the IJOC GitHub software repository ( https://github.com/INFORMSJoC/2025.1210 ). The complete IJOC Software and Data Repository is available at https://informsjoc.github.io/ .

Estrada-Garcia, Juan-Alberto↗
DOE OSTI2020

Gaussian Process Optimization of Sensitivity-Based Similarity Metrics between New Nuclear Applications and New/Existing Benchmarks

Confidently designing safe, new nuclear criticality experiments requires expert judgement, which could take years of experience. Sensitivity/uncertainty (S/U) analysis can be utilized by less experienced individuals to conservatively estimate uncertainties in important parameters, such as k eff , in newly proposed nuclear applications. This type of analysis relies on matching new nuclear applications with existing benchmark experiments. The Whisper-1.1 software package included in MCNP6.2 ®1 contains more than 1,100 International Criticality Safety Benchmark Evaluation Project (ICSBEP) benchmarks. These benchmarks however rarely match new nuclear applications. The number of benchmarks available to match a given set of materials or geometric configurations varies significantly. Furthermore, recently performed benchmark experiments may not have had enough time to be properly documented and published. Benchmarks are vital for determining the accuracy of nuclear data and can assist nuclear physics and evaluators in improving nuclear data libraries. Exhaustively exploring the parameter space using simulations with software such as MCNP is too computationally expensive. In this work, Gaussian process optimization was implemented to reduce the number of simulations needed for optimization over multiple parameters. This optimization scheme was designed to selectively generate new benchmarks with high sensitivity-based similarity metrics to user-defined nuclear applications. Two benchmark models of spherically nested shells containing plutonium, uranium, tantalum, and water were used in the optimization to match an application containing plutonium plates, stacked in a tantalum reflector, surrounded by water.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Deep-learning-aided forward optical coherence tomography endoscope for percutaneous nephrostomy guidance

Percutaneous renal access is the critical initial step in many medical settings. In order to obtain the best surgical outcome with minimum patient morbidity, an improved method for access to the renal calyx is needed. In our study, we built a forward-view optical coherence tomography (OCT) endoscopic system for percutaneous nephrostomy (PCN) guidance. Porcine kidneys were imaged in our experiment to demonstrate the feasibility of the imaging system. Three tissue types of porcine kidneys (renal cortex, medulla, and calyx) can be clearly distinguished due to the morphological and tissue differences from the OCT endoscopic images. To further improve the guidance efficacy and reduce the learning burden of the clinical doctors, a deep-learning-based computer aided diagnosis platform was developed to automatically classify the OCT images by the renal tissue types. Convolutional neural networks (CNN) were developed with labeled OCT images based on the ResNet34, MobileNetv2 and ResNet50 architectures. Nested cross-validation and testing was used to benchmark the classification performance with uncertainty quantification over 10 kidneys, which demonstrated robust performance over substantial biological variability among kidneys. ResNet50-based CNN models achieved an average classification accuracy of 82.6%±3.0%. The classification precisions were 79%±4% for cortex, 85%±6% for medulla, and 91%±5% for calyx and the classification recalls were 68%±11% for cortex, 91%±4% for medulla, and 89%±3% for calyx. Interpretation of the CNN predictions showed the discriminative characteristics in the OCT images of the three renal tissue types. The results validated the technical feasibility of using this novel imaging platform to automatically recognize the images of renal tissue structures ahead of the PCN needle in PCN surgery.

Wang, Chen↗
DOE OSTI2023

Integrated photonic molecule Brillouin laser with a high-power sub-100-mHz fundamental linewidth

Photonic integrated lasers with an ultra-low fundamental linewidth and a high output power are important for precision atomic and quantum applications, high-capacity communications, and fiber sensing, yet wafer-scale solutions have remained elusive. Here we report an integrated stimulated Brillouin laser (SBL), based on a photonic molecule coupled resonator design, that achieves a sub-100-mHz fundamental linewidth with greater than 10-mW output power in the C band, fabricated on a 200-mm silicon nitride (Si 3 N 4 ) CMOS-foundry compatible wafer-scale platform. The photonic molecule design is used to suppress the second-order Stokes (S2) emission, allowing the primary lasing mode to increase with the pump power without phase noise feedback from higher Stokes orders. The nested waveguide resonators have a 184 million intrinsic and 92 million loaded Q, over an order of magnitude improvement over prior photonic molecules, enabling precision resonance splitting of 198 MHz at the S2 frequency. We demonstrate S2-suppressed single-mode SBL with a minimum fundamental linewidth of 71±18 mHz, corresponding to a 23±6-mHz 2 /Hz white-frequency-noise floor, over an order of magnitude lower than prior integrated SBLs, with an ∼11-mW output power and 2.3-mW threshold power. The frequency noise reaches the resonator-intrinsic thermo-refractive noise from 2-kHz to 1-MHz offset. The laser phase noise reaches −155 dBc/Hz at 10-MHz offset. The performance of this chip-scale SBL shows promise not only to improve the reliability and reduce size and cost but also to enable new precision experiments that require the high-speed manipulation, control, and interrogation of atoms and qubits. Realization in the silicon nitride ultra-low loss platform is adaptable to a wide range of wavelengths from the visible to infrared and enables integration with other components for systems-on-chip solutions for a wide range of precision scientific and engineering applications including quantum sensing, gravitometers, atom interferometers, precision metrology, optical atomic clocks, and ultra-low noise microwave generation.

Liu, Kaikai (ORCID:0000000213609646)↗
DOE OSTI2026

Semi-Analytical Hierarchical Bayesian Inference of Nonlinear Model Structure in Stochastic Dynamics: Applied to Compartmental Models of Infectious Diseases

A Bayesian computational framework for parsimonious inference in stochastic nonlinear dynamical systems is presented. This framework enables the concurrent estimation of system states, time-varying parameters, time-invariant parameters, and the optimal sparsity structure of the model parameters. Because differential equation-based models are often simplified mechanistic or phenomenological representations, robust inference from noisy measurement data requires explicit treatment of model error and uncertainty. Model error and time-varying parameters can be represented as random processes, enabling inference while making minimal assumptions about the underlying sources of discrepancy and variability. Adopting stochastic differential equation representations affords the model significant flexibility, but can also render it susceptible to overfitting during statistical inversion, where the inferred model may track noise rather than the underlying signal. To alleviate the effects of overfitting and to enable the discovery of the optimal sparse representation of the time-invariant parameters, a Bayesian sparse learning algorithm is embedded within the framework. This sparse learning framework adopts an approximate hierarchical Bayesian setting defined by a series of semi-analytical expressions. The model structure inference framework is validated using a stochastic compartmental model for tracking and forecasting active cases of an infectious disease. Compartmental models describe population-level infectious disease dynamics through interactions among population fractions grouped by disease state. Mathematically, such models consist of a system of coupled ordinary differential equations. This example adopts an expressive compartmental model that includes multiple possible interactions between disease states, motivated by early uncertainty surrounding COVID-19 reinfection dynamics and their implications for long-term epidemic forecasting. The sparse learning exercise permits the inference of a priori unknown epidemiological dynamics from simulated public health data, discovering the nested compartmental model that optimizes the trade-off between average data-fit and model complexity. It is shown that inducing sparsity among the model parameters eliminates redundant interactions between compartments, equivalently revealing the optimal coupling structure between differential equations.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Overview of the MOSAiC expedition: Snow and sea ice

Year-round observations of the physical snow and ice properties and processes that govern the ice pack evolution and its interaction with the atmosphere and the ocean were conducted during the Multidisciplinary drifting Observatory for the Study of Arctic Climate (MOSAiC) expedition of the research vessel Polarstern in the Arctic Ocean from October 2019 to September 2020. This work was embedded into the interdisciplinary design of the 5 MOSAiC teams, studying the atmosphere, the sea ice, the ocean, the ecosystem, and biogeochemical processes. The overall aim of the snow and sea ice observations during MOSAiC was to characterize the physical properties of the snow and ice cover comprehensively in the central Arctic over an entire annual cycle. This objective was achieved by detailed observations of physical properties and of energy and mass balance of snow and ice. By studying snow and sea ice dynamics over nested spatial scales from centimeters to tens of kilometers, the variability across scales can be considered. On-ice observations of in situ and remote sensing properties of the different surface types over all seasons will help to improve numerical process and climate models and to establish and validate novel satellite remote sensing methods; the linkages to accompanying airborne measurements, satellite observations, and results of numerical models are discussed. We found large spatial variabilities of snow metamorphism and thermal regimes impacting sea ice growth. We conclude that the highly variable snow cover needs to be considered in more detail (in observations, remote sensing, and models) to better understand snow-related feedback processes. The ice pack revealed rapid transformations and motions along the drift in all seasons. The number of coupled ice–ocean interface processes observed in detail are expected to guide upcoming research with respect to the changing Arctic sea ice.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Evaluation of simulations of near-surface variables using the regional climate model CCLM for the MOSAiC winter period

The ship-based experiment MOSAiC 2019/2020 was carried out during a full year in the Arctic and yielded an excellent data set to test the parameterizations of ocean/sea-ice/atmosphere interaction processes in regional climate models (RCMs). In the present paper, near-surface data during MOSAiC are used for the verification of the RCM COnsortium for Small-scale MOdel–Climate Limited area Mode (COSMO-CLM or CCLM). CCLM is used in a forecast mode (nested in ERA5) for the whole Arctic with 15 km resolution and is run with different configurations of sea ice data. These include the standard sea ice concentration taken from passive microwave data with around 6 km resolution, sea ice concentration from Moderate Resolution Imaging Spectroradiometer (MODIS) thermal infrared data and MODIS sea ice lead fraction data for the winter period. CCLM simulations show a good agreement with the measurements. Relatively large negative biases for temperature occur for November and December, which are likely associated with a too large ice thickness used by CCLM. The consideration of sea ice leads in the sub-grid parameterization in CCLM yields improved results for the near-surface temperature. ERA5 data show a large warm bias of about 2.5°C and an underestimation of the temperature variability.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Storm and Annual Time Scale Hydrological Data for the Russian River Watershed 1996-2022

These files contains observed and simulated hydrological data (discharge, precipitation) that have been aggregated to the storm-event scale and to the annual time scale for the Russian River Watershed (RRW), California. Observed data are obtained from 12 different USGS hydrological stations located throughout the watershed, and discharge is simulated using the GR5H hourly hydrological model. Precipitation data are obtained from the NASA NLDAS NOAH Community Land Surface Model. All simulated and observed discharge are aggregated to the storm and annual time scales. We investigated changes in watershed hydrological conditions in the Russian River Watershed, a Mediterranean, drought prone, wildfire-adapted ecosystem, following eleven wildfires that occurred from 2017-2020. We ask two research questions: 1) How do interacting wildfire events, drought, and atmospheric rivers impact hydrological conditions of the watershed and in particular streamflow?, and 2) What percentage of new and/or cumulative wildfire disturbance is required to initiate hydrological change? We hypothesize that sub-watersheds of the RRW Mediterranean ecoregion have not burned beyond an intrinsic, and still unknown, threshold required to initiate change. Using a series of paired burned/unburned catchments nested within the larger watershed, we examined temporal and spatial patterns of pre-and-post wildfire water hydrology using a rainfall runoff hydrological model compared with data.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Groundwater level elevation and temperature data, Oct 2018-Dec 2021, Slate River Floodplain, Crested Butte, CO

This data package includes a time series of water level and temperature measurements from October 2018 to December 2021 in groundwater and surface water from the Slate River floodplain outside Crested Butte, CO, a focus field site for the SLAC Floodplain Hydro-Biogeochemistry SFA. The data was generated as part of the work targeting the overarching research question for the SLAC SFA: How do ubiquitous subsurface interfaces mediate molecular-scale biogeochemical processes and groundwater quality in floodplains and watersheds? The data were recorded by pressure transducers were installed in four types of piezometers: 1) a network of gravel bed ("GB") drive-point piezometers with a 6-inch screen interval installed all at ~330 cm below ground surface. 2) a network of piezometers screened across the water table, used to measure evapotranspiration ("ET") using the White method. Each of these piezometers is screened along almost its entire length.3) a suite nested piezometers used the measure the vertical hydraulic gradient ("VHG") across the fines-cobble interface. Each of these piezometers uses a 6-inch screen length.4) a group of piezometers scattered across the boundaries of the floodplain, used to monitor boundary conditions ("BC") flowing into and out of the floodplain. With the exception of "SR-BD-WT", each of these piezometers is screened along its entire length. Within the data package, "FLMD.csv" describes file-level metadata and "dd.csv" defines column headers and universal terms across the dataset. The data package includes 11 "*data.csv" files, one for each piezometer type for each calendar year. Because piezometers have been added over time, not every sensor has data dating back to Oct 2018. Each "*data.csv" file has a corresponding "*_InstallationMethods.csv" file that describes the location, elevation, screen depth, sediment type and sensor metadata for each piezometer and pressure transducer.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Field and Model Data Associated with the Manuscript “Drivers of Streamflow Intermittency in Humid Regions: 2. Evaluating Controls on Flow Persistence in an Urbanized Catchment”

This package contains field data, modeling files, and scripts supporting the investigation of the drivers of streamflow intermittency in an urbanized catchment. It includes the field data collected from electrical resistivity tomography (ERT) surveys, distributed temperature sensing (DTS), continuous self-potential (SP) monitoring, groundwater and stilling well. In addition, it contains the data and results of the coupled water- and electrical-flow model developed using the COMSOL Multiphysics and Advanced Terrestrial Simulator (ATS), as well as software files and Jupyter notebooks used to process the data and generate figures in the manuscript submitted for peer review. The data archive is organized in the following directories: 1) Climate Includes hourly precipitation and daily evapotranspiration time series (2024 – 2025) provided as CSV files, alongside a text file detailing dataset units. 2) Coupled_model Field_Application subfolder contains the ATS XML input scripts, data files, output data for the SP site. It also contains the Jupyter notebook (Plot_final_calib.ipynb) to visualize the results of the modeled SP, stream-groundwater exchange and moisture content. The flow model simulation is executed using the ATS XML scripts and the included Python script (generate_data_set.py) to convert ATS output to COMSOL-ready input. COMSOL Multiphysics template (.m can only be used with COMSOL with MATLAB) is executed using the ATS output data to simulate the potential field. 3) Discharge Includes the electrical conductivity (EC) time series (provided as CSV files) from salt slug injections. It also includes the Jupyter notebook (Discharge_process.ipynyb) used to estimate discharge. All discharge measurements collated into rating_curve_processed.csv 4) DTS Contains collated DTS data including raw Stokes and anti-Stokes measurement (provided as .h5 file). It also includes DTS processing.ipynb, a Jupyter notebook for calibrating the DTS data using dts_calibration Python package. cooler_calibration.csv is the DTS calibration CSV used in the calibration sequence. 5) ERT Contains raw resistivity data (provided as CSV files), spatial location of each of the electrodes (provided as CSV files), and files used for the resistivity inversion. 6) Slug_test Includes the slug test data at all the groundwater wells provided as CSV files, as well as the Jupyter notebook (Slug_test.ipynb) for calculating hydraulic conductivity. 7) SP Contains the SP data collected in field at the SP sites (provided as CSV files). 8) Well_data Contains two subfolders: 1) Raw, which provides unprocessed pressure, electrical conductivity and temperature timeseries downloaded from the loggers in all the groundwater and stilling wells, and 2) Processed, which contains sorted, QA/QC timeseries data for each well. The data archive also contains data_process.ipynb, a Jupyter notebook used for field data analysis and generating figures (plotting well, SP, climate, and discharge data, as well as calculating head gradient at sites with nested groundwater wells). Note: Code files (.ipynb, .py, .xml) can be opened in any standard code editor, .exo file can be viewed using Paraview, .h5 files can be opened using HDFView software and h5py Python package, and .resipy file can be opened with the open-source ResIPy software.

ATS↗
DOE OSTI2026

Field and Model Data Associated with the Manuscript “Drivers of Streamflow Intermittency in Humid Regions: 1. Evaluating Above- and Below-ground Controls of Flow Persistence in a Forested Catchment”

This package contains field data, modeling files, and scripts supporting the investigation of the drivers of streamflow intermittency in a forested catchment. It includes the field data collected from electrical resistivity tomography (ERT) surveys, ground penetrating radar (GPR), continuous self-potential (SP) monitoring, electromagnetic (EM) imaging, groundwater and stilling well. In addition, it contains the data and results of the coupled water- and electrical-flow model developed using the COMSOL Multiphysics and Advanced Terrestrial Simulator (ATS), as well as software files and Jupyter notebooks used to process the data and generate figures in the manuscript submitted for peer review. The data archive is organized in the following directories: 1) Climate Includes hourly precipitation and daily evapotranspiration time series (2024 – 2025) provided as CSV files, alongside a text file detailing dataset units. 2) Coupled_model Contains two subfolders: Synthetic and Field_Application subfolder. Synthetic subfolder contains the ATS XML input script (can be opened using any code editor) for the four synthetic hydrological cases tested (Connected and gaining, Connected and losing, Disconnected and losing, and dry stream). It also includes other experimental cases to test the influence of precipitation and concentration gradient. For each synthetic case, the flow model simulation is executed using the ATS XML scripts and the included Python script (generate_data_set.py) to convert ATS output to COMSOL-ready input. COMSOL Multiphysics template (.mph can be opened with the commercial software COMSOL and requires a license) is executed using the ATS output data to simulate the potential field. It also includes the Synthetic_model_plot.ipynb (can be opened using any code editor) to visualize the SP result and generate manuscript figures. The data subfolder contains mesh files to run both the ATS (.exo and .stl files can be viewed using Paraview; .h5 files can be opened using HDFView software and h5py Python package) and COMSOL models. Field_Application subfolder contains two subfolders: ES_MDA_inversion and Final_Model. ES_MDA_inversion contains the Python script (.py can be opened using any code editor) and SP observation data used to run the Ensemble Smoother with Multiple Data Assimilation (ES-MDA) inversion sequence to get the optimal model parameters. The Final_model subfolder contains the ATS XML input scripts, data files, output data for the two SP sites. The same workflow steps outlined for the Synthetic subfolder apply here. It also contains the Jupyter notebook (Plot_final_calib.ipynb) to visualize the results of the modeled SP, stream-groundwater exchange and moisture content. 3) Discharge Includes the electrical conductivity (EC) time series (provided as CSV files) from salt slug injections. It also includes the Jupyter notebook (Discharge_process.ipynyb) used to estimate discharge. All discharge measurements collated into rating_curve_processed.csv 4) EM Contains the CSV file of the EM data from the DUALEM-42, including spatial coordinates (x, y, z), apparent conductivity, and in-phase measurements at 2 m coil separations for horizontal coplanar (HCP) and perpendicular (PRP) geometries. 5) ERT Contains raw resistivity data (provided as CSV files), spatial location of each of the electrodes (provided as CSV files), and files used for the resistivity inversion (.resipy can be opened with the open-source ResIPy software). 6) GPR Includes GPR field datasets collected at 100 MHz and 250 MHz antenna frequencies, along with the processing/interpretation project file (GPR_process.gpz can be viewed using EKKO_Project 6, a commercial software by Sensors & Software that requires a license). 7) Slug_test Includes the slug test data at all the groundwater wells provided as CSV files, as well as the Jupyter notebook (Slug_test.ipynb) for calculating hydraulic conductivity. 8) SP Contains the SP data collected in field at the two SP sites (one in the perennial reach and the other in the intermittent reach), provided as DAT files. 9) Well_data Contains two subfolders: 1) Raw, which provides unprocessed pressure, electrical conductivity and temperature timeseries downloaded from the loggers in all the groundwater and stilling wells, and 2) Processed, which contains sorted, QA/QC timeseries data for each well. The data archive also contains data_process.ipynb, a Jupyter notebook used for field data analysis and generating figures (plotting well, SP, climate, and discharge data, as well as calculating head gradient at sites with nested groundwater wells). It also includes DTW.ipynb, a Jupyter notebook containing the code for the dynamic time warping (DTW) with sliding window to evaluate SP signal synchronicity.

ATS↗