SEARCH · Engineering Papers

Results for “response surface methods”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 487 records · Page 27

All sources NASA NTRS DOE OSTI

DOE OSTI2025

Groundwater flowpath characteristics drive variability in per- and polyfluoroalkyl substances (PFAS) loading across a stream-wetland system

Groundwater dependent ecosystems in areas with industrial and military land use are at risk of direct exposure to a wide range of contaminants, including PFAS chemicals. Glaciated terrain often has mixed high and low permeability sediments coupled with groundwater flow-through lake features. These hydrogeologic attributes create highly complex ‘source to seep’ dynamics that make spatiotemporal contaminant transport patterns difficult to predict. We investigated one such system in detail using a suite of heat-tracing and chemical methods. Numerous (n=57) preferential groundwater discharge zones (vertical flux rates ranging 0.2 to 3.2 m/d) were identified across the upper Quashnet River stream-wetland system in Mashpee, MA, USA, adjacent to an Air Force Base with several known PFAS source areas. Surface-water and groundwater samples were collected and analyzed (n=145) for precursors and terminal PFAS compounds between March and September 2022. Samples were collected at identified seeps along the Quashnet River (n=59), from wells upgradient from the stream-wetland system (n= 44), from contributing flow-through kettle lakes (n=8), and at multiple locations along the Quashnet River (n=34). Samples from seeps and wells had measured PFAS concentrations ranging from non-detect to approximate 3,500 ng/L (mean= 1,650 ng/L), and a range of deuterium excess values (3.2 to 15.9 per mil) indicative of varying degrees of groundwater-lake interaction prior to emergence at the discharge zones. Groundwater-lake interaction along flowpaths that sourced the sampled seeps was farther supported by significant correlations (p < 0.01) between deuterium excess and %PFAS precursors, and between %PFAS precursors and multiple terminal PFAS compounds (e.g., PFPeS, PFBS, PFHxS). However, some sampled seeps contributing groundwater to the stream-wetland system had much higher total PFAS concentrations (>1000 ng/L) than the upgradient kettle lakes, despite showing lake (evaporative) isotopic signatures, indicating the potential for groundwater flowpath convergence at wetland discharge zones and the influence of lakebed PFAS precursor reactions. PFAS compounds and water isotopic composition at sampled multilevel groundwater wells, rivers, lakes, and seeps suggest that a complex mixture of source groundwater and flowpath characteristics are responsible for diverse observed PFAS mixtures at preferential discharge zones across the stream-wetland system. Further, total PFAS loading patterns to the Quashnet River via groundwater discharge remained remarkably similar from winter to summer to fall conditions, despite a regional dry period in late summer 2022 with the upper river channel completely drying. This work addresses gaps in the existing PFAS literature by demonstrating the importance of subsurface fate and transport on PFAS compound concentrations in controlling contaminant mass loading in preferential groundwater discharge zones and presents a transferrable field toolkit for efficient characterization of spatially preferential PFAS transport dynamics.

Contaminant transport↗

DOE OSTI2022

Reinforcement Learning for feedback-enabled cyber resilience

The rapid growth in the number of devices and their connectivity has enlarged the attack surface and made cyber systems more vulnerable. As attackers become increasingly sophisticated and resourceful, mere reliance on traditional cyber protection, such as intrusion detection, firewalls, and encryption, is insufficient to secure the cyber systems. Cyber resilience provides a new security paradigm that complements inadequate protection with resilience mechanisms. A Cyber-Resilient Mechanism (CRM) adapts to the known or zero-day threats and uncertainties in real-time and strategically responds to them to maintain the critical functions of the cyber systems in the event of successful attacks. Feedback architectures play a pivotal role in enabling the online sensing, reasoning, and actuation process of the CRM. Reinforcement Learning (RL) is an important gathering of algorithms that epitomize the feedback architectures for cyber resilience. It allows the CRM to provide dynamic and sequential responses to attacks with limited or without prior knowledge of the environment and the attacker. In this work, we review the literature on RL for cyber resilience and discuss the cyber-resilient defenses against three major types of vulnerabilities, i.e., posture-related, information-related, and human-related vulnerabilities. Here we introduce moving target defense, defensive cyber deception, and assistive human security technologies as three application domains of CRMs to elaborate on their designs. The RL algorithms also have vulnerabilities themselves. We explain the major vulnerabilities of RL and present develop several attack models where the attacker target the information exchanged between the environment and the agent: the rewards, the state observations, and the action commands. We show that the attacker can trick the RL agent into learning a nefarious policy with minimum attacking effort. The paper introduces several defense methods to secure the RL-enabled systems from these attacks. However, there is still a lack of works that focuses on the defensive mechanisms for RL-enabled systems. Last but not least, we discuss the future challenges of RL for cyber security and resilience and emerging applications of RL-based CRMs.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2020

Toward Understanding the Simulated Phase Partitioning of Arctic Single-Layer Mixed-Phase Clouds in E3SM

21 Mixed-phase clouds are widespread in the Arctic and have substantial impacts on the 22 surface energy budget, sea ice extent, and regional climate. Significant changes are found in the 23 modeled phase partitioning of Arctic mixed-phase clouds in the newly developed U.S. 24 Department of Energy (DOE) Energy Exascale Earth System Model (E3SM) Atmosphere Model 25 version 1 (EAMv1) compared to its predecessor (CAM5/EAMv0). The mixed-phase clouds 26 simulated by EAMv1 is overly dominated by super-cool liquid which is often underestimated by 27 CAM5. Compared to CAM5, EAMv1 has made several important changes in its cloud physical 28 parameterizations, including using CLUBB to unify shallow convection, PBL turbulence, and 29 cloud macrophysics, as well as updating the MG two-moment cloud microphysics scheme and 30 using the CNT heterogeneous ice nucleation scheme to replace the Meyers scheme for mixed- 31 phase clouds. To understand which of the newly-implemented physical parameterizations is 32 responsible for the change of modeled mixed-phase cloud phase partitioning in EAMv1, 33 sensitivity tests are conducted by utilizing the short-term hindcast approach and results are 34 compared to the observations from the ARM Mixed-Phase Arctic Cloud Experiment (M-PACE) 35 field campaign. 36 It is found that although each scheme contributes to the decrease of cloud ice water 37 content (IWC) and increase of cloud liquid water content (LWC), only the combined effect of 38 CLUBB and CNT can explain the significant underestimation of IWC. By examining the budgets 39 of cloud microphysical processes, CLUBB is found to play a dominant role due to the fact that 40 CLUBB does not diagnoses cloud ice phase water in its assumed joint PDF method and thus may 41 not provide sufficient cloud ice to initialize ice phase microphysical processes such as the growth 42 of ice crystals through vapor deposition. Although a more effective heterogeneous ice nucleation 43 may alleviate the underestimation of IWC, the ice phase process in CLUBB requires urgent 44 attention in the future model development.

Zhang, Meng↗

DOE OSTI2020

Evaluating the Effectiveness of an Ultrasonic Acoustic Deterrent in Reducing Bat Fatalities at Wind Energy Facilities

This project was designed to use thermal video cameras and fatality monitoring to evaluate the effectiveness of an ultrasonic acoustic deterrent (UAD) on bat activity and mortality, respectively. Our goals were to redesign the UAD device and installation infrastructure, determine the placement on wind turbines to optimize safety, compatibility and functionality, and to compare the mortality among the following conditions: Control (deterrents off and turbines feathered up to the manufacturer’s cut-in speed of 3.5 m/s), Deterrent (deterrents on and turbines feathered up to the manufacturer’s cut-in speed of 3.5 m/s), Curtailment (deterrents off and turbines feathered up to 5/ m/s), and combination (deterrents on and turbines feathered up to 5 m/s). The project was divided into a Feasibility Study and Comparative Study. The objectives for the Feasibility Study were to develop an installation strategy, redesign a previous iteration of a UAD to improve performance and weatherization, and test the effectiveness of the deterrents on bat activity. The Feasibility Study was intended to work out potential issues using a relatively small number of devices prior to manufacturing and installing numerous devices for a larger-scale comparative study. The division of the project into these separate studies was based on previous experience and the challenges of assessing the capabilities of an untested UAD. During the initial development and manufacturing of the UAD, NRG Systems decided to use a piezoelectric transducer rather than an electrostatic transducer, which was used by a previous device (i.e. Deaton UAD). NRG Systems conducted lab testing (i.e., IP67 or Ingress Protection) to ensure no water or dust ingress. In addition, shock/drop trials and variations in temperature exposure were conducted as part of the reliability testing. NRG Systems also developed a communications system to allow for continuous performance monitoring of the UADs. For the Feasibility Study, we installed 6 UADs on each of 2 Gamesa G90 2-MW wind turbines (Turbine 14 and 8) at the South Chestnut Wind Energy Facility, Pennsylvania and monitored activity under control and treatment (i.e. Deterrent) conditions using thermal video monitoring. There are two major sources of variation in bat activity (beyond the anticipated treatment effect): 1) environment around the turbines might inherently favor more activity at one than the other; 2) weather conditions on any given night or within season difference (e.g., migration later during the study period) might favor more activity on some nights than on other nights. Because we could only monitor two turbines on any night, we sought to control the potential influences of these two sources by alternating the turbine on which deterrents were activated each night. If there were no loss of data due to technical failures, this design would result in an equal number of deterrent and control nights at each turbine through the monitoring period, balancing the effects of both sources of variation. We compared the time bats spent and the number of events that occurred in overlapping cameras FOV as an indicator of risk, since 80% of the overlapping FOV of the cameras was in the RSA. Equipment failures, majority due to lightning, resulted in only 17 nights with useable data, with unbalanced treatment assignment within turbines and uneven distribution of treatment assignment throughout the observational period. This resulted in a confounding of treatment assignment and seasonal change. Deterrent treatment was measured at Turbine 14 on only 2 of the first 8 usable nights (spanning the period from 8/18-9/17), whereas 6 times on Turbine 8. From 9/18-927, deterrent was on at Turbine 14 on 6 of the remaining 10 nights, and 4 on Turbine 8. We recorded a total of 1,057 bats and observed a reduction in number of events and duration of events at the UAD-activated turbine when it was Turbine 14. When Turbine 8 had the UAD activated, we observed no difference in number or duration of events. Variation between turbines is not unusual and can cause issues when study designs have no true replication, i.e., multiple turbines per treatment. Within-turbine differences suggested a trend for reduced activity when UADs were activated, particularly for Turbine 14. These results may be caused by the overall higher bat activity at Turbine 8 and confounding of treatment assignment and seasonal trends. We mapped 58 bat events in 3-dimensional space (3D), 30 and 28 events during control and treatment conditions, respectively. We observed bats crossing the rotor plane under both control and treatment conditions and observed a total of 40 confirmed or near-collisions (i.e., target close to blade but no visual confirmation of a strike) out of a total of 1,491 medium and high confidence bat observations (880 control, 611 at treatment). Twice as many collisions/possible collisions were observed during control conditions. Given the challenges with the equipment and potential confounding of the data (i.e., different activity levels at the two wind turbines), we were unable to determine whether this initial turbine placement and orientation was optimal. Given no new information on how best to install the devices, we elected to use the same placement and orientation for the comparative study. For the Comparative Study, the objectives were to investigate the relative mortality rates among 4 treatments. We searched the area within 90 m of each turbine daily to recover the highest number of fresh fatalities possible. We were unable to detect a clear reduction in mortality from deterrents alone for any individual species. Surprisingly, mortality rate of the eastern red bat (Lasiurus borealis) was estimated to be 1.3–4.2 times as much when turbines were operating normally and UADs were on than when UADs were off. Reduction in mortality of all bat species combined due to curtailment of turbines was estimated to be between 0%–38%. This effect was nullified when, in addition to curtailment, UADs were on, with 95% confidence interval ranging from a 45% reduction to a 36% increase in mortality. This was likely due to the large proportion of eastern red bats in the total carcass population. Mortality of all low-frequency echolocating bats combined (i.e. hoary bat [L. cinereus], big brown bat [Eptesicus fuscus], silver-haired bat [Lasionycteris noctivagans]) relative to control was lower when curtailed (95% CI: 0%–74%), but the addition of UADs had no detectable effect (95%CI: 13%–79%). The combined treatment reduced mortality in silver-haired bats relative to control by 11%–99%, compared to curtailment (81% reduction–67% increase) or deterrent (82% reduction–67% increase) alone. Because silver-haired bats comprised a large proportion of low-frequency calling bats found during this study, a similar effect was seen for that group. The higher mortality observed for eastern red bats at UAD compared to control could have been caused by several factors, such as the effective range of the UAD, particularly at higher frequencies, behavior, positioning of the devices on the nacelle, or a combination of these. We used 3D thermal videography to compare control and UAD bat behavior from two turbines using a total of 203 3D bat-tracks across 34 nights. We recorded a similar number of bat-tracks between treatment groups, with 51% and 49% for control and UAD, respectively. Due to potential differences in bat behavior around spinning vs stationary turbine blades, we examined UAD effectiveness separately for non-operating turbines (feathered below cut-in speed of 3.5 m/s) and operating (normal operation above wind speed of 3.5 m/s). We found a higher proportion of bat-tracks at operating turbines (82%) compared to non-operating turbines (18%), although this does not account for overall time turbines were operating versus not. At non-operating turbines the UAD appears to be effective at reducing the amount of time, flight length, and number of passes through the rotor plane, compared to control. In addition, we found bats approached turbines similarly between control and treatment turbines, with 61% of control and 63% of UAD bat-tracks originating leeward of the hub. In contrast, at operating turbines, we saw little change in bat behavior in response to UADs. For example, we found an increase in the average duration of bat-tracks between non-operating and operating turbines for UAD but at control turbines we found average duration decreased once turbines became operational. Both control and UAD had a high proportion of bat-tracks that crossed the rotor plane (i.e., collision risk) originate from the windward side when turbines were operational 65% to 92%, respectively. Given that the UAD devices closest to the blades were orientated parallel to the blades, its possible bats were not exposed to the signal until they were close to the turbine blades, as suggested by the slightly higher mean duration within 5 meters of the blades for UAD turbines. Future research should consider concentrating UAD intensity on the areas of risk (i.e. blades) with enough buffer to allow bats to react to the sound before entering the rotor-swept area (RSA). In addition, investigating the potential of installing UAD units windward of the turbine blades (e.g. a hub-mounted UAD), particularly since even under control conditions, a relatively high proportion (65%) of crosses through the blade plane originated windward. 3D thermal videography provided valuable information on future testing strategies (e.g. device placement, UAD orientation) to improve UAD effectiveness when bats are at risk (i.e., operating wind turbines). Across the entire project we experienced issues with the operation and communication with the UADs. Most of the issues occurred during the Feasibility Study and were resolved prior to the Comparability Study. Additional challenges surfaced during the Comparability Study but were remedied immediately and are thought to have little impact on the results. We had logistical constraints at the project that limited our ability use traditional methods in our camera calibration. Several calibrations showed inaccurate scales, which may have been related to inadequate spatial coverage of “points” in the camera calibration volume. Because we were using actual video recordings of bats at a wind turbine, the behavior of bats could have concentrated “points” in specific areas of the turbine (i.e. leeward of nacelle), and limited “points” in other areas, resulting in camera calibration issues. We have plans to address these inconsistencies and improving the software and related methodologies by early 2020.

17 WIND ENERGY↗

DOE OSTI2024

Regional climate change: consensus, discrepancies, and ways forward

Climate change has emerged across many regions. Some observed regional climate changes, such as amplified Arctic warming and land-sea warming contrasts have been predicted by climate models. However, many other observed regional changes, such as changes in tropical sea surface temperature and monsoon rainfall are not well simulated by climate model ensembles even when taking into account natural internal variability and structural uncertainties in the response of models to anthropogenic radiative forcing. This suggests climate model predictions may not fully reflect what our future will look like. The discrepancies between models and observations are not well understood due to several real and apparent puzzles and limitations such as the “signal-to-noise paradox” and real-world record-shattering extremes falling outside of the possible range predicted by models. Addressing these discrepancies, puzzles and limitations is essential, because understanding and reliably predicting regional climate change is necessary in order to communicate effectively about the underlying drivers of change, provide reliable information to stakeholders, enable societies to adapt, and increase resilience and reduce vulnerability. The challenges of achieving this are greater in the Global South, especially because of the lack of observational data over long time periods and a lack of scientific focus on Global South climate change. To address discrepancies between observations and models, it is important to prioritize resources for understanding regional climate predictions and analyzing where and why models and observations disagree via testing hypotheses of drivers of biases using observations and models. Gaps in understanding can be discovered and filled by exploiting new tools, such as artificial intelligence/machine learning, high-resolution models, new modeling experiments in the model hierarchy, better quantification of forcing, and new observations. Conscious efforts are needed toward creating opportunities that allow regional experts, particularly those from the Global South, to take the lead in regional climate research. This includes co-learning in technical aspects of analyzing simulations and in the physics and dynamics of regional climate change. Finally, improved methods of regional climate communication are needed, which account for the underlying uncertainties, in order to provide reliable and actionable information to stakeholders and the media.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2022

Detection of surface water temperature variations of Mongolian lakes benefiting from the spatially and temporally gap-filled MODIS data

Lakes provide critical water resources for human activities and ecosystems, particularly in the Mongolian Plateau (MP), which is characterized by a dry climate and a harsh environment. As a region that is sensitive to anthropogenic warming, tracking lake surface water temperature (LSWT) changes in Mongolian lakes is crucial for understanding the consequences of a warming climate on lake ecosystems. However, the long-term monitoring of LSWT is restricted by the spatiotemporal gaps in the raw imagery of remote sensing-based land surface temperature (LST), e.g., the commonly used Moderate Resolution Imaging Spectroradiometer (MODIS) LST products. This study applied an improved gap-filling method by utilizing the discrete cosine transform-based penalized least squares (DCT-PLS) strategy in the spatial domain combined with the linear interpolation (LI) algorithm in the temporal domain. The method was applied to fill gaps in the LSWT imagery of 12 representative lakes across MP. The randomly sampled high-quality MODIS LSWT values in the spatial and temporal domains were excavated as false data gaps and considered “virtual true” validation datasets. The spatial validation results showed that the estimated LSWT for all the lake cases were comparable with the “virtual true” LSWT values, with the average values of the coefficient of determination, mean absolute error, mean square error, and root mean square error being 0.98, 0.38 °C, 0.45 °C, and 0.59 °C, respectively. Meanwhile, the error of nighttime LSWT results was relatively lower than that of daytime LSWT. For temporal interpolation validation, the LI algorithm exhibited relatively better performance and could more objectively indicate the variation in LSWT. Benefiting from the spatially and temporally well-constrained data, we analyzed the interannual and intra-annual change characteristics of the LSWTs of the 12 lakes. The long-term variations of annual and seasonal mean LSWTs in the 12 selected lakes exhibited no evident trends in 2000–2020, while presented apparent interannual fluctuations. The slight changes in the average LSWTs of the 12 selected lakes were in excellent synchronization with the surrounding LST derived from the reanalysis datasets, confirming the widely reported phenomenon of “global warming hiatus” that occurred in the early 21st century. This study improves the understanding of the LSWT variations in Mongolian lakes in response to global climate change. It has the potential to provide an effective approach for monitoring LSWT changes in other large-scale studies.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2023

Solution combustion synthesis of iron-based alumina catalysts for dehydrogenation of fossil fuels

Microwave-assisted thermocatalytic dehydrogenation of fossil fuels has been proposed as a method for hydrogen production with no CO2 emissions. The byproduct is carbon, which has a commercial value. Iron-based alumina nanocomposites are promising catalysts for this process as they both absorb microwave radiation and catalyze the dehydrogenation. In addition, they may have magnetic properties, which could be used for their separation from the carbon byproduct. Solution combustion synthesis (SCS) is an attractive technique for the fabrication of iron-based alumina nanocomposites because it can produce high specific surface area oxides in a facile manner. However, the effects of the heating mode, the fuel, and the iron/aluminum ratio on the combustion synthesis behavior and the product properties have not been studied yet. In the present work, the iron-based alumina nanocomposites were obtained by SCS using a hotplate, a muffle furnace, and a microwave oven. The initial mixtures were aqueous solutions of iron nitrate, aluminum nitrate, and an organic fuel. The concentrations of the two nitrates were varied in proportion to Fe:Al molar ratios of 1:2, 1:1, and 2:1. Two fuels were tested: citric acid and glycine. The combustion was less vigorous when citric acid was used as the fuel, which is explained by the lower exothermicity of the reaction. X-ray diffraction analysis of the products confirmed the formation of complex oxides of Fe and Al, specifically it detected hercynite and magnetite phases. Brunauer-Emmett-Teller surface area analysis has shown that the powders obtained using citric acid have specific surface areas as high as 276 m2/g, significantly higher than those obtained using glycine. Muffle furnace and hotplate heating led to comparable specific surface areas, while microwave heating resulted in significantly lower specific surface areas. The magnetic properties of the products increase with increasing the Fe/Al molar ratio from 1:2 (no response to a magnetic field) to 1:1 (weak response) and 2:1 (strong response).

solution combustion synthesis, nanoscale oxides, f↗

DOE OSTI2021

Evaluating fine-resolution, regional outputs of a variable resolution global climate model

Climate models have been used to study water resources and regional hydrologic responses to climate change, but climate model outputs must be downscaled to provide relevant regional data. However, the accuracy of this regional data is limited by uncertainties across and within downscaling methods, uncertainty across global outputs, and discontinuities at downscaled boundaries. A new alternative to traditional downscaling is a variable resolution model that incorporates fine-resolution regions directly into a coarse-resolution, global climate simulation in order to capture contiguous dynamics across resolution boundaries. In this study, we used the Variable-Resolution Community Earth System Model (VR-CESM) to generate one-eighth degree (14 km) fine-resolution outputs for the western U.S. and eastern China from 1970-2006. We focus our evaluation on precipitaiton, temperature, snow pack, solar radiation, and wind. We compare the model outputs with remote-sensing-based precipitation data, and both reanalysis and gridded weather station data for precipitation and temperature. VR-CESM generally has a cold bias in winter and a warm bias in summer in the western U.S., which compensate each other to reduce the annual bias. In eastern China, however, the sign of temperature biases are more consistent throughout the year with cold biases in the higher mountains and warm biases throughout most of the rest of the region. Precipitation biases are dependent upon reference data, and show slight overestimation in high mountain regions in both the U.S. and China with respect to gridded weather station data. Simulated snow cover in the western U.S. is reasonable compared to remote sensing data, but snow cover and snow water equivalent have larger biases when compared to reanalysis data. In eastern China there are widespread snow cover biases compared to remote sensing data. VR-CESM underestimates downward shortwave radiation to a greater degree in summer than in winter, and underestimates surface layer windspeed over mountains to a greater degree than in other areas. Comparison between VR-CESM and a coarser simulation (1-degree Beijing Climate Center model) shows reduced precipitation biases in the mountainous regions with finer resolution, indicating the value of variable-resolution modeling for reigonal studies.

Di Vittorio, Alan↗

DOE OSTI2023

Reversible Switch in Charge Storage Enabled by Selective Ion Transport in Solid Electrolyte Interphase

Solid-electrolyte interphases (SEIs) in advanced rechargeable batteries ensure reversible electrode reactions at extreme potentials beyond the thermodynamic stability limits of electrolytes by insulating electrons while allowing working ions to transport. Such selective ion transport occurs naturally in biological cell membranes as a ubiquitous prerequisite of many life processes and a foundation of biodiversity. In addition, cell membranes can selectively open and close the ion channels in response to external stimuli (e.g., electrical, chemical, mechanical, thermal), giving rise to “gating” mechanisms that help manage intracellular reactions. We wondered whether the chemistry and structure of SEIs can mimic cell membranes, such that ion gating can be replicated. That is, can SEIs realize a reversible switching between two electrochemical behaviors, i.e., the ion intercalation chemistry of batteries and the ion adsorption of capacitors? Herein, we report such SEIs that result in thermally activated selective ion transport. The function of open/close gate switches is governed by the chemical and structural dynamics of SEIs under different thermal conditions, with precise behaviors as conducting and insulating interphases that enable battery and capacitive processes within a finite temperature window. Such an ion gating function is synergistically contributed by Arrhenius-activated ion transport and SEI dissolution/regrowth. Following the understanding of this new mechanism, we then develop an electrochemical method to heal the SEI layer in situ. As a result, the knowledge acquired in this work reveals the possibility of hitherto unknown biomimetic properties of SEIs, which will guide us to leverage such complexities to design better SEIs for future battery chemistries.

25 ENERGY STORAGE↗

DOE OSTI2024

WEC-Sim Modeling of Laminar Scientific's Patented Seesaw Wave Energy Converter: Preprint

Laminar Scientific's patented seesaw wave energy converter was modeled in WEC-Sim to predict performance. The device operates by utilizing ocean surface waves to rotate a truss in pitch about a pivot. The pivot is located at the top of two pylons, which are embedded in the seafloor. The seesaw has a float on either end, and the buoyancy forces from each float cause the system to rise or fall with passing waves. Device performance relies upon seesaw length and ocean wavelength creating an antiphase effect. The seesaw truss has an adjustable length intended to achieve this effect. The operation method enforces a narrow band of wavelengths which induce the largest rotational motion from the device. The hydrodynamic analysis of the device was performed using Capytaine, and the results confirmed that the device operates best in a narrow frequency band. Four float-to-float spacing cases and three pylon radii were examined. The hydrodynamic results indicate a match between the model and the physical expectations for the device, and that varying the pylon radii by 0.1-m increments for three instances creates minimal changes in hydrodynamic properties. Power matrices for three float spacing cases of the device were simulated with Joint North Sea Wave Project spectra waves and optimal power take-off damping in WEC-Sim. The maximum average power production for the 15-m spacing case was 14.1 kW with a 5.0-s peak wave period and 4-m significant wave heights. Plots of capture-width ratios indicated that the device performance was linear and confirmed that the device is optimal in a narrow frequency band. The maximum percentage of the available wave power produced by the 15-m device was approximately 16%. Simulations of the device in regular waves were used to produce plots of average power compared to a ratio of float spacing to wavelength. These plots indicate that the power production is maximized at a ratio of 0.5, and further confirm that the device has a narrow frequency response. The device was simulated at an example field location, where the device produced an annual average power rating of 1.6 kW given an average omnidirectional wave climate of 10.3 kW m-1 and an optimal, linearized power take-off model. While the maximum predicted device performance is reliant upon a narrow band of wave frequencies, the conducted analysis provides an opportunity to improve device design prior to prototyping and testing. Modifying the design to respond to a broader frequency range would improve device performance.

numerical modeling↗

DOE OSTI2025

Second-generation downscaled earth system model data using generative machine learning

The second-generation Sup3rCC dataset provides high-resolution meteorological data generated through the downscaling of multiple earth system models (ESMs) from the Coupled Model Intercomparison Project Phase 6 (CMIP6). This downscaling is performed through application of a generative machine learning approach called Super-Resolution for Renewable Resource Data (sup3r). This dataset builds on the first-generation Sup3rCC data by applying improved bias correction methods and adding downscaled precipitation to the output variables. As with the first Sup3rCC version, the data still include temperature, wind speed and direction at multiple heights, pressure, three components of downwelling solar radiation, and relative humidity—all at 4-kilometer (km) hourly resolution over the contiguous United States. This is a 25x spatial enhancement and 24x temporal enhancement of the source 100-km daily-average ESM data. This extension of the Sup3rCC dataset includes data from six ESMs from two shared socioeconomic pathways (SSPs) totaling 400 years of data with multiple future projections of changing meteorological conditions. The scenario selection was based on a structured evaluation of historical ESM skill and comprehensive representation of possible trajectories of future climate change in temperature, humidity, precipitation, solar irradiance, and near-surface wind speeds. The inclusion of multiple future projections is intended to enable users to assess key drivers of un 36 certainty and variability. All data are double-bias corrected, resulting in a product that can be used out-of-the-box for energy system analysis with minimal historical bias. The potential applications of Sup3rCC data extend to various topics in renewable energy resource assessment, energy systems modeling, and grid resilience studies. High-resolution future meteorological projections are critical for evaluating the effects of changing meteorological conditions on renewable energy generation, energy demand, and for optimizing energy storage and grid infrastructure. The 4-km hourly resolution of the downscaled data enables understanding of spatial and temporal variability at the scales necessary for energy system operational planning. In addition, the dataset can support risk assessments by providing detailed information on possible future extreme weather events and long-term meteorological variability at scales relevant to energy infrastructure. By offering an enhanced representation of possible future meteorological conditions, the second-generation Sup3rCC dataset enables more precise modeling of energy resilience and adaptation strategies in response to changing meteorological conditions.

24 POWER TRANSMISSION AND DISTRIBUTION↗

DOE OSTI2020

Geophysical monitoring using active seismic techniques at the Citronelle Alabama CO 2 storage demonstration site

Between August 2012 and September 2014, about 114,000 metric tonnes of CO 2 was captured from the coal-fired Plant Barry Power Station at Bucks Alabama and injected into the Paluxy Formation above the oil pool in the southeast unit of the Citronelle Oilfield. Various monitoring methods were deployed at land surface and in project wells to measure system performance, comply with permit requirements and test new and innovative monitoring tools. The monitoring program relied heavily on active seismic methods for subsurface imaging of geologic structure and time-lapse seismic techniques to track the CO 2 migration in the injection interval. Both conventional geophone/hydrophone and fiber-optic based Distributed Acoustic Sensing (DAS) arrays were deployed and tested, allowing a side by side comparison of the equipment and techniques. Geophysical imaging of the subsurface was successful using DAS in the offset vertical seismic profile (OVSP) survey configuration. A high resolution OVSP image of the subsurface was obtained in 2014 with DAS, which exceeded project expectations in comparison to a lower resolution image obtained in 2012 using a conventional 80-level geophone array. A time-lapse image of the redistribution of CO 2 after injection ended in September 2014 was obtained with two DAS OVSP surveys from June 2014 and December 2015, thus successfully demonstrating its proof-of-concept. Unfortunately, a pre-injection baseline survey with DAS, which was in its initial stage of technology development in 2012, did not have sufficient quality for use, making it difficult to interpret the acquired DAS time-lapse difference. Additional research in this area has since demonstrated the utility of time-lapse DAS OVSP. DAS data were also acquired during a cross-well seismic survey conducted in 2014. Unfortunately, the DAS technique was not success in the cross-well survey configuration because the system noise level was too high in the crosswell frequency output range (100–1200 Hz) of the piezoelectric source (increasing by a factor of ten compared to VSP frequency band). Additionally, the cross-well geometry causes sub-horizontal (broadside) incidence on the vertical DAS fiber cable, which is known to be problematic. Current research is focused on improving the DAS cable response to broadside acoustic energy. Time-lapse seismic surveys using commercially available conventional arrays were also acquired. In contrast to the DAS acquired data, the cross-well seismic results obtained with the conventional array was highly successful and clearly showed the CO 2 remained in zone at the end of injection. Time-lapse differencing of the OSVP surveys acquired with the conventional arrays proved to be inconclusive. Finally, changes in wellbore conditions between surveys and unavoidable changes in equipment (the array used for the baseline survey was retired) affected data quality, making it difficult to interpret the OVSP results.

58 GEOSCIENCES↗

DOE OSTI2025

Design and Production for Maximum Structural Efficiency With Respect to Fiber Orientation With Increased Understanding of Hybrid Fiber Flow Behavior

Discontinuous fiber-reinforced thermoplastic composites have gained considerable attention in automotive, aerospace, and other industries, due to their high-rate of production combined with their ability to attain complex and intricate shapes. Among other high-rate thermoplastic manufacturing processes, injection-molding is one of the most common manufacturing methods due to fast production and high surface finishing of complex geometries. Fiber orientation in discontinuous fiber composites plays a pivotal role in determining the mechanical, electrical, and thermomechanical properties, underscoring the necessity to comprehend fiber orientation in injection molded parts. Among different fiber types, glass and carbon fibers are most common in the composite industries. The recent trend of hybrid composites comprising both glass fiber (GF) and carbon fiber (CF) is also gaining importance in the automotive industry. Hybrid fiber options allow designers to optimize the balance between glass and carbon fibers by leveraging the high durability and low cost of GF while the strength and lightweight properties of CF. Consequently, comparing the fiber oriented distribution (FOD) of injection molded composites containing GF, CF, and a hybrid of GF/CF is critical to investigating the local mechanical properties of intricate structures for high-end applications. In Phase I of this project, FOD in injection-molded panels with respect to distance from the gate was analyzed using X-ray computed tomography (X-CT) for GF, CF, and hybrid CF/GF (CGF) reinforced nylon 66. To understand the reason behind the FOD with different fiber types, computational fluid dynamics (CFD) and rheology were performed. Samples were extracted at three locations: near the gate, center, and opposite end. Thickness of the layers of typical skin-shell-core type FOD varies with fiber type and location. GF achieved flow direction alignment (in shell) earlier than viscous CF and CGF near the gate, whereas CF showed the highest flow-direction alignment at the center due to shear induced orientation. At the opposite end, GF experienced more backflow than others indicating faster mold filling owing to its lower viscosity. Hybrid CGF exhibited GF-dominated center and CF-dominated end region. The numerical model used to obtain FOD and rheological predictions for the CF and GF composites served to corroborate the trends observed in the experimental trials. The FOD responses across fiber types and location were reflected in their longitudinal and transverse properties. Only GF showed higher longitudinal modulus over transverse modulus near the gate attributed to rapid alignment, whereas CF and CGF exhibited opposite trend. However, fountain flow enhanced the longitudinal modulus over transverse modulus with the distance for all, particularly for CF. This study offers insights into mold filling behavior of different fibers which are critical in optimizing injection molding conditions for tailored final properties.

36 MATERIALS SCIENCE↗

DOE OSTI2024

Global Carbon Budget 1800-2018 Loss of Additional Sink Capacity and Present Transient Difference

A dataset of two global carbon fluxes: the loss of additional sink capacity (LASC) and present versus transient difference (PTD). The loss of additional sink capacity is the loss of indirect anthropogenic sink capacity in land ecosystems caused by the loss of ecosystems that can increase carbon in response to changing environments, primarily forests. The present transient difference refers to the direct anthropogenic flux in land ecosystems (from land cover and land use change) and the difference in that flux when calculated assuming "present day" (actually more like 1980) ecosystem carbon stocks or when calculated assuming transient carbon stocks (i.e. caused by environmental change). In this dataset, the LASC and PTD fluxes were calculated using the TRENDY ensemble of Dynamic Global Vegetation Models (DGVMs; Sitch et al., 2024) following the methods described therein and in Obermeier et al. (2021). For further details see Obermeier et al. (2021), Friedlingstein et al. (2019), Sitch et al. (2024), and Walker et al. (2025). The dataset includes ensemble means and standard deviations for annual fluxes, fluxes cumulated annually over the whole time period, and cumulated since 1959. These data were originally collected to compare different approaches (DGVMs versus bookkeeping models) for estimating land use and land cover change (LULCC) emissions (Obermeier et al. 2021). More recently, and similarly, these data were used to harmonize the estimates of the direct and indirect anthropogenic fluxes in the global carbon cycle and to recalculate the global carbon budget (Walker et al., 2025).

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2022

Kimberlina 1.2 CCUS Geophysical Models and Synthetic Data Sets

This synthetic multi-scale and multi-physics data set was produced in collaboration with teams at the Lawrence Berkeley National Laboratory, National Energy Technology Laboratory, Los Alamos National Laboratory, and Colorado School of Mines through the Science-informed Machine Learning for Accelerating Real-Time Decisions in Subsurface Applications (SMART) Initiative. Data are associated with the following publication: Alumbaugh, D., Gasperikova, E., Crandall, D., Commer, M., Feng, S., Harbert, W., Li, Y., Lin, Y., and Samarasinghe, S., “The Kimberlina Synthetic Geophysical Model and Data Set for CO2 Monitoring Investigations”, The Geoscience Data Journal, 2023, DOI: 10.1002/gdj3.191. The dataset uses the Kimberlina 1.2 CO2 reservoir flow model simulations based on a hypothetical CO2 storage site in California (Birkholzer et al., 2011; Wainwright et al., 2013). Geophysical properties models (P- and S-wave seismic velocities, saturated density, and electrical resistivity) were produced with an approach similar to that of Yang et al. (2019) and Gasperikova et al. (2022) for 100 Kimberlina 1.2 reservoir models. Links to individual resources are provided below: [CO2 Saturation Models](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-co2-saturation-models); Resistivity Models – [part 1](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-resistivity-models-part-1), [part 2](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-resistivity-models-part-2), and [part 3](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-resistivity-models-part-3); [Vp Velocity Models](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-vp-velocity-models); [Vs Velocity Models](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-vs-velocity-models); [Density Models](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-density-models). The 3D distributions of geophysical properties for the 33 time stamps of the SIM001 model were used to generate synthetic seismic, gravity, and electromagnetic (EM) responses for 33 times between zero and 200 years. Synthetic surface seismic data were generated using 2D and 3D finite-difference codes that simulate the acoustic wave equation (Moczo et al., 2007). 2D data were simulated for six point-pressure sources along a 2D line with 10 m receiver spacing and a time spacing of 0.0005 s. 3D simulations were completed for 25 surface pressure sources using a source separation of 1 km in both the x and y directions and a time spacing of 0.001 s. Links to individual resources are provided below: [2D velocity models](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-2d-velocity-models) and [2D surface seismic data](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-2d-surface-seismic-data). [3D velocity models](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-velocity-models), and 3D seismic data [year0](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year0), [year1](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year1), [year2](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year2), [year5](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year5), [year10](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year10), [year15](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year15), [year20](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year20), [year25](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year25), [year30](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year30), [year35](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year35), [year40](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year40), [year45](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year45), [year49](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year49), [year50](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year50), [year51](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year51), [year52](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year52), [year55](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year55), [year60](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year60), [year65](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year65), [year70](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year70), [year75](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year75), [year80](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year80), [year85](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year85), [year90](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year90), [year95](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year95), [year100](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year100), [year110](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year110), [year120](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year120), [year130](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year130), [year140](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year140), [year150](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year150), [year175](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year175), [year200](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-seismic-data-year200). The Python scripts to read these models and data are provided [here](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-python-scripts). EM simulations used a borehole-to-surface survey configuration, with the source located near the reservoir level and receivers on the surface using the code developed by Commer and Newman (2008). Pseudo-2D data for the source at [2500 m](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-pseudo-2d-csem-data-tz2500m) and [3025 m](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-pseudo-2d-csem-data-tz3025m), used a 2D inline receiver configuration to simulate a response over 3D resistivity models. The [3D data](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-3d-csem-data) contain electric fields generated by borehole sources at monitoring well locations and measured over a surface receiver grid. Vector gravity data, both on the surface and in boreholes, were simulated using a modeling code developed by Rim and Li (2015). The simulation scenarios were parallel to those used for the EM: [pseudo-2D data](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-gravity-data) were calculated along the same lines and within the same boreholes, and [3D data](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-gravity-data) were simulated over 3D models on the surface and in three monitoring wells. A series of [synthetic well logs](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-well-logs) of CO2 saturation, acoustic velocity, density, and induction resistivity in the injection well and three monitoring wells are also provided at 0, 1, 2, 5, 10, 15, and 20 years after the initiation of injection. These were constructed by combining the low-frequency trend of the geophysical models with the high-frequency variations of actual well logs collected in the Kimberlina 1 well that was drilled at the proposed site. Measurements of permeability and pore connectivity were made on cores of Vedder Sandstone, which forms the primary reservoir unit: [CT micro scans](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-ct-micro-scans-of-vedder-formation) and [Industrial CT Images](https://edx.netl.doe.gov/dataset/kimberlina-1-2-ccus-geophysical-models-and-synthetic-data-sets-industrial-ct-images-vedder-formation). These measurements provide the range of scales in the otherwise synthetic data set to be as close to a real-world situation as possible. References: Birkholzer, J.T., Zhou, Q., Cortis, A. and Finsterle, S., 2011. A sensitivity study on regional pressure buildup from large-scale CO2 storage projects. Energy Procedia, 4, 4371-4378. Commer, M., and Newman, G.A., 2008. New advances in three-dimensional controlled-source electromagnetic inversion, Geophysical Journal International, 172, 513-535. Gasperikova, E., Appriou, D., Bonneville, A., Feng, Z., Huang, L., Gao, K., Yang, X., Daley, T., 2022, Sensitivity of geophysical techniques for monitoring secondary CO2 storage plumes, Int. J. Greenh. Gas Control, Volume 114, 103585, ISSN 1750-5836, https://doi.org/10.1016/j.ijggc.2022.103585. Moczo, P., J.O. Robertsson and L. Eisner, 2007, The finite-difference time-domain method for modeling of seismic wave propagation: Advances in geophysics, 48, 421-516. Rim, H., and Y. Li, 2015, Advantages of borehole vector gravity in density imaging, Geophysics, 80, G1-G13. Wainwright, H. M.; Finsterle, S.; Zhou, Q.; Birkholzer, J. T., 2013. Modeling the Performance of Large-Scale CO2 Storage Systems: A Comparison of Different Sensitivity Analysis Methods. International Journal of Greenhouse Gas Control, 17, 189205. https://doi.org/10.1016/j.ijggc.2013.05.007, DOI: 10.18141/1603331. Yang, X., Buscheck, T.A., Mansoor, K., Wang, Z., Gao, K., Huang, L., Appriou, D., and Carroll, S.A., 2019. Assessment of geophysical monitoring methods for detection of brine and CO2 leakage in drinking water aquifers, International Journal of Greenhouse Gas Control, 90, 102803, https://doi.org/10.1016/j.ijggc.2019.102803.

CCUS↗

DOE OSTI2023

Patterns and rates of soil movement and shallow failures across several small watersheds on the Seward Peninsula, Alaska

Abstract. Thawing permafrost can alter topography, ecosystems, and sediment and carbon fluxes, but predicting landscape evolution of permafrost-influenced watersheds in response to warming and/or hydrological changes remains an unsolved challenge. Sediment flux and slope instability in sloping saturated soils have been commonly predicted from topographic metrics (e.g., slope, drainage area). In addition to topographic factors, cohesion imparted by soil and vegetation and melting ground ice may also control spatial trends in slope stability, but the distribution of ground ice is poorly constrained and hard to predict. To address whether slope stability and surface displacements follow topographic-based predictions, we document recent drivers of permafrost sediment flux present on a landscape in western Alaska that include creep, solifluction, gullying, and catastrophic hillslope failures ranging in size from a few meters to tens of meters, and we find evidence of rapid and substantial landscape change on an annual timescale. We quantify the timing and rate of surface movements using a multi-pronged, multi-scalar dataset including aerial surveys, interannual GPS surveys, synthetic aperture radar interferometry (InSAR), and climate data. Despite clear visual evidence of downslope soil transport of solifluction lobes, we find that the interannual downslope surface displacement of these features does not outpace downslope displacement of soil in locations where lobes are absent (downslope movement means: 7 cm yr−1 for lobes over 2 years vs. 10 cm yr−1 in landscape positions without lobes over 1 year). Annual displacements do not appear related to slope, drainage area, or modeled total solar radiation but are likely related to soil thickness, and volumetric sediment fluxes are high compared to temperate landscapes of comparable bedrock lithology. Time series of InSAR displacements show accelerated movement in late summer, associated with intense rainfall and/or deep thaw. While mapped slope failures do cluster at slope–area thresholds, a simple slope stability model driven with hydraulic conductivities representative of throughflow in mineral and organic soil drastically overpredicts the occurrence of slope failures. This mismatch implies permafrost hillslopes have unaccounted-for cohesion and/or throughflow pathways, perhaps modulated by vegetation, which stabilize slopes against high rainfall. Our results highlight the breadth and complexity of soil transport processes in Arctic landscapes and demonstrate the utility of using a range of synergistic data collection methods to observe multiple scales of landscape change, which can aid in predicting periglacial landscape evolution.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2021

Relaxation behavior in low-frequency complex conductivity of sands caused by bacterial growth and biofilm formation by Shewanella oneidensis under a high-salinity condition

Complex electrical conductivity is increasingly used to monitor subsurface processes associated with microbial activities because microbial cells mostly have surface charges and thus electrical double layers. Although highly saline environments are frequently encountered in coastal and marine sediments, there are limited data available on the complex conductivity associated with microbial activities under a high-salinity condition. Therefore, we have developed the spectral responses of complex conductivity of sand associated with bacterial growth and biofilm formation under a highly saline condition of approximately 1% salinity and approximately 2 S/m pore water conductivity with an emphasis on relaxation behavior. A column test is performed, in which the model bacteria Shewanella oneidensis MR-1 are stimulated for cell growth and biofilm formation in a sand pack, whereas the complex conductivity is monitored from 0.01 Hz to 10 kHz. The test results indicate that the real conductivity increases in the early stage due to the microbial metabolites and the increased surface conduction with cell growth but soon begin to decrease because of the reduction of charge passages due to bioclogging. However, the imaginary conductivity significantly increases with time, and clear bell-shaped relaxation behaviors are observed with the peak frequency of 0.1–1 Hz, associated with the double-layer polarization of cells and electrically conductive pili and biofilms. The Cole-Cole relaxation model appears to capture such relaxation behaviors well, and the modeling results indicate gradual increases in normalized chargeability and decreases in relaxation time during bacterial growth and biofilm formation in the highly saline condition. Comparison with previous literature confirms that the high-salinity condition further increases the normalized chargeability, whereas it suppresses the phase shift and thus the imaginary conductivity. Our results suggest that the complex conductivity can effectively capture microbial biomass formation in sands under a highly saline condition.

58 GEOSCIENCES↗

DOE OSTI2026

Post-fire time series of sensor and geochemistry sample data from surface water, groundwater, precipitation, soil, and vegetation across Oak Creek watershed, Washington

This dataset supports a broader study examining wildfire impacts on hydrologic connectivity across 5 sites within the Oak Creek watershed and the resulting biogeochemical impacts. Stream sites were selected using the Advanced Terrestrial Simulator (ATS) hydrologic model to identify locations with varying groundwater contributions and hydrologic responses across different burn severity scenarios. The Retreat Fire burned from July 23 to August 2 in 2024, affecting the five study sites at varying burn severities. Each site is equipped with YSI EXO2 sondes logging sub-hourly throughout the year, and grab samples are collected approximately every six weeks. YSI sondes are used to measure temporally resolved proxies for groundwater inputs (specific conductivity) and organic matter (fluorescent dissolved organic matter; fDOM) along with basic water quality and depth. Grab samples of surface water, groundwater, and precipitation are analyzed for water stable isotopes and conductivity to understand endmembers for hydrologic mixing Grab samples of surface water, groundwater, soil water, and litter/vegetation/soil leachates are analyzed for organic matter composition measured by Fourier-Transform Ion Cyclotron Resonance Mass Spectrometry (FTICR-MS) to understand organic matter dynamics. Game camera photos are provided in a separate data package available at https://data.ess-dive.lbl.gov/datasets/doi:10.15485/3018598. Future versions of this dataset will include time series data from YSI EXO2 sondes (fDOM, dissolved oxygen, temperature, depth, specific conductance, turbidity, pH), BaroTROLL sensors (air temperature and barometric pressure), rain gauges (precipitation), and data from the soil and vegetation samples. Because this study is ongoing, this data package will be updated regularly to include newly collected data and the additional data types. For details on how to navigate data packages generated by this project, see https://data.ess-dive.lbl.gov/portals/PNNLRiverCorridorSFA/About. In addition to a readme, this data package also includes a file-level metadata (FLMD) file that describes each file and a data dictionary (DD) that describes all column/row headers and variable definitions. This dataset is comprised of (1) a folder of field photos; (2) a folder of surface water sample data; (3) a folder of raw Fourier transform ion cyclotron resonance mass spectrometry (FTICR-MS) data; (4) a data checks report; (5) file-level metadata; (6) data dictionary; (7) field metadata; (8) readme; (9) international generic sample number (IGSN) mapping file; and (10) field protocols. The sample data subfolder contains (1) dissolved organic carbon (DOC, measured as non-purgeable organic carbon, NPOC) data and averages; (2) total dissolved nitrogen data and averages; (3) stable water isotopes and averages; (4) methods codes; (5) FTICR-MS methods; and (15) a subfolder of 9.4 Tesla (9.4T) FTICR-MS data. This folder contains the processed data and three subfolders, one containing the .xml files, one containing the water CoreMS output files, and the other containing instructions and scripts for processing the files in CoreMS (https://github.com/EMSL-Computing/CoreMS). All files are .csv, .pdf, .R, .xml, .d, .html, .Rmd, .py, .cal, .json, .jpg, .jpeg, .png, .mov, or .mp4.

Biogeochemistry↗