Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “data gap analysis”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

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

At least 1,099 records · Page 61
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Technology Development for Dry Storage of Aluminum-Clad Spent Nuclear Fuel - 20490

A candidate disposition pathway for the > 13 MTHM of aluminum-clad spent nuclear fuel (ASNF), owned and managed by the U.S. Department of Energy, is the drying and placement of the SNF into sealed-canister dry storage, with the ASNF-in- canisters 'road-ready' for transportation to and final direct disposal in a repository waste package. Technical information gaps in fuel drying, and fuel dry storage behavior, have previously challenged the declaration of technology readiness for drying and placement of this fuel into the DOE Standard Canister design for > 50 years of safe dry storage. The principal technical information gaps included: i) characterization and thermal dehydration behavior of aluminum (oxy)hydroxide films attendant on the cladding due to film formation during reactor operation and during post-discharge up to long-term wet storage histories; and ii) G-value data to enable estimation of the radiolytic gas generation from the cladding with its (oxy)hydroxide films. Thus, the oxide films on the ASNF challenged the safety of a sealed storage canister with thermal and radiolytic decomposition of the waters on the films that can lead to corrosion, pressurization, and flammability issues. These gaps in the technical information base have largely been closed. This paper discusses the investigations at the Idaho National Laboratory (INL) and the Savannah River National laboratory (SRNL), and outlines the pending technology development work for input to an engineering design to enable a road-ready dry storage system for ASNF. The ASNF inventory considered for road-ready dry storage is stored at the Savannah River Site (SRS) and at the Idaho Nuclear Technology and Engineering Center (INTEC) at the INL. The ASNF inventory in the SRS L Basin is from foreign and domestic research reactors (FRR and DRR), and is diverse in terms of design, irradiation, and post-reactor-discharge storage conditions; these factors yield a range of characteristics of cladding oxide films on ASNF. Mixed aluminum (oxy)hydroxide (boehmite and bayerite/gibbsite) films, non-uniform in thicknesses up to a maximum local thickness bounded by 25 μm, were observed on ASNF materials removed from wet storage in the L Basin and in non-sealed dry storage at the SRS following reactor service and a long-term (up to 40+ years) interim storage history. The ASNF inventory at INTEC, in both wet and vented dry storage, is predominantly from the Advanced Test Reactor (ATR), but it also includes DRR and FRR fuel. To address a profound behavior of these films, radiolysis testing of aluminum specimens with mixed type boehmite/bayerite oxide films was performed to develop basic data on G-values for production of radiolytic hydrogen under dry storage conditions with nominal relative humidity, temperature, and cover gas. Modeling and simulation of canister internal environments with postulated inventories of oxide films provides estimation of the evolution of the conditions of the canister loaded with ASNF. Simultaneous Thermal Analysis (TGA/DSC) of hydrated oxide powders, and laboratory-scale aluminum specimens with a bayerite film (∼10 μm) using TGA methods, inform drying time/temperature conditions to be used for the ASNF. A demonstration project is recommended for Verification and Validation of the drying and storage of the ASNF. Remaining major tasks leading up to the hot demonstration include scale-up radiolysis testing and scale-up drying testing. Engineering design with the information from the technology program will establish the safety basis and enable long-term (> 50 years) dry storage compatibility with ASNF in the DOE Standard Canister pending its transportation to and disposal in a repository. This full capability would show the ASNF-in-canister storage system to be road-ready. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
NASA NTRS2001

Soft X-Ray Properties of the Binary Millisecond Pulsar J0437-4715

We obtained a light curve for the 5.75 ms pulsar J0437-4715 in the 65-120 A range with 0.5 ms time resolution using the Deep Survey instrument on the EUVE satellite. The single-peaked profile has a pulsed fraction of 0.27 +/- 0.05, similar to the ROSAT data in the overlapping energy band. A combined analysis of the EUVE and ROSAT data is consistent with a power-law spectrum of energy index alpha = 1.2 - 1.5, intervening column density N(sub H) = (5 - 8) x 10(exp 19)/sq cm, and luminosity 5.0 x 10(exp 30) ergs/s in the 0.1 - 2.4 keV band. We also use a bright EUVE/ROSAT source only 4.2 min. from the pulsar, the Seyfert galaxy RX J0437.4-4711 (= EUVE J0437-471 = IES 0435-472), to obtain an independent upper limit on the intervening absorption to the pulsar, N(sub H) less than 1.2 x 10(exp 20)/sq cm. Although a blackbody spectrum fails to fit the ROSAT data, two-component spectral fits to the combined EUVE/ROSAT data are used to limit the temperatures and surface areas of thermal emission that might make partial contributions to the flux. A hot polar cap of radius 50 - 600 m and temperature (1.0 - 3.3) x 10(exp 6) K could be present. Alternatively, a larger region with T = (4 - 12) x 10(exp 5) K and area less than 200 sq km, might contribute most of the EUVE and soft X-ray flux, but only if a hotter component were present as well. Any of these temperatures would require some mechanism(s) of surface reheating to be operating in this old pulsar, the most plausible being the impact of accelerated electrons and positrons onto the polar caps. The kinematically corrected spin-down power of PSR J0437-4715 is only 4 x 10(exp 33) ergs/s, which is an order of magnitude less than that of the lowest-luminosity gamma-ray pulsars Geminga and PSR B1055-52. The absence of high-energy gamma-rays from PSR J0437-4715 might signify an inefficient or dead outer gap accelerator, which in turn accounts for the lack of a more luminous reheated surface such as those intermediate-age gamma-ray pulsars may have.

Halpern, Jules P.↗
NASA NTRS2021

The Scientific Legacy of NASA’s Operation IceBridge

The National Aeronautics and Space Administration (NASA)’s Operation IceBridge (OIB) was a 13-year (2009–2021) airborne mission to survey land and sea ice across the Arctic, Antarctic, and Alaska. Here, we review OIB’s goals, instruments, campaigns, key scientific results, and implications for future investigations of the cryosphere. OIB’s primary goal was to use airborne laser altimetry to bridge the gap in fine-resolution elevation measurements of ice from space between the conclusion of NASA’s Ice, Cloud, and land Elevation Satellite (ICESat; 2003–2009) and its follow-on, ICESat-2 (launched 2018). Additional scientific requirements were intended to contextualize observed elevation changes using a multisensor suite of radar sounders, gravimeters, magnetometers, and cameras. Using 15 different aircraft, OIB conducted 968 science flights, of which 42% were repeat surveys of land ice, 42% were surveys of previously unmapped terrain across the Greenland and Antarctic ice sheets, Arctic ice caps, and Alaskan glaciers, and 16% were surveys of sea ice. The combination of an expansive instrument suite and breadth of surveys enabled numerous fundamental advances in our understanding of the Earth’s cryosphere. For land ice, OIB dramatically improved knowledge of interannual outlet-glacier variability, ice-sheet, and outlet-glacier thicknesses, snowfall rates on ice sheets, fjord and sub-ice-shelf bathymetry, and ice-sheet hydrology. Unanticipated discoveries included a reliable method for constraining the thickness within difficult-to-sound incised troughs beneath ice sheets, the extent of the firn aquifer within the Greenland Ice Sheet, the vulnerability of many Greenland and Antarctic outlet glaciers to ocean-driven melting at their grounding zones, and the dominance of surface-melt-driven mass loss of Alaskan glaciers. For sea ice, OIB significantly advanced our understanding of spatiotemporal variability in sea ice freeboard and its snow cover, especially through combined analysis of fine-resolution altimetry, visible imagery, and snow radar measurements of the overlying snow thickness. Such analyses led to the unanticipated discovery of an interdecadal decrease in snow thickness on Arctic sea ice and numerous opportunities to validate sea ice freeboards from satellite radar altimetry. While many of its data sets have yet to be fully explored, OIB’s scientific legacy has already demonstrated the value of sustained investment in reliable airborne platforms, airborne instrument development, interagency and international collaboration, and open and rapid data access to advance our understanding of Earth’s remote polar regions and their role in the Earth system.

Joseph A. MacGregor↗
DOE OSTI2023

Composite Analysis for Low-Level Waste Disposal in the Hanford Site Central Plateau (FY 2020) (Rev. 2)

This document provides an updated Hanford Site composite analysis (CA). A Hanford Site CA was prepared and issued in 1998 (PNNL-11800) with an addendum provided in 2002 (PNNL-11800 Addendum 1). The CA was approved in 2002 (Frei, 2002) and has been maintained (DOE/RL-2000-29 and subsequent revisions) to support low-level waste disposal performance assessments (PAs) and disposal authorizations for facilities at the Hanford Site, including the following: Continued operation of the Environmental Restoration Disposal Facility (ERDF) and the 200 East and 200 West Low-Level Burial Grounds; Construction of the Integrated Disposal Facility (IDF); Forthcoming closure of tank residual waste systems such as Waste Management Area (WMA) C. The CA maintenance program resulted in a determination in 2015 (DOE/RL-2015-66) that the Hanford Site CA needed an update for the following reasons: While the initial Hanford Site CA has been maintained since 2001, the accumulation of basis changes reported in the annual summary reports over the succeeding 14 years merit evaluation in an updated analysis; The U.S. Department of Energy (DOE) Headquarters requested in a memorandum in 2015 (Gilbertson and Marcinowski, 2015) that “as soon as the relevant PAs are complete, the CA will be revised to account for all of the new information.” This updated Hanford Site CA accounts for the following new information: 1. Inclusion of a detailed Hanford Site baseline disposition that projects remedial activities through site closure. There have been significant changes through decision making in the Comprehensive Environmental Response, Compensation, and Liability Act of 1980 (CERCLA) process that were not available when the original CA was produced; 2. Inclusion of an updated inventory basis, new modeling capabilities, and new decisions reached in the associated record of decision (ROD) that was provided by issuance of DOE/EIS-0391 in fiscal year (FY) 2013. Development of a Hanford Site groundwater model from the baseline provided in a technical transfer of models for the Final Tank Closure and Waste Management Environmental Impact Statement (DOE/EIS-0391) commenced in FY 2014 and resulted in the plateau-to-river (P2R) groundwater flow model that is used in this updated Hanford Site CA; 3. Accounting for pump and treat systems, which were not evaluated in the initial CA. Pump and treat systems have had significant impact on groundwater flow system behavior, contaminant transport, and contaminant removal from Hanford Site groundwater. This process is accounted for in the use of the P2R model for this updated Hanford Site CA; 4. Inclusion of water-level data collected since the initial CA was completed. Data collected as the unconfined aquifer water levels continue to recede since the cessation of large liquid discharges in the late 1990s have led to marked improvement in understanding of the flow system for future conditions, particularly regarding northward flow potential in the critical Gable Gap area. The calibration in the latest version of the P2R model accounts for this information; 5. Inclusion of results from the revised ERDF PA completed in 2013 (WCH-5209) that account for updated inventory and expansion of the ERDF facility to about twice the size that was evaluated in the original CA; 6. Use of updated geoframeworks to provide the structural basis for numerical fate and transport models in the groundwater pathway of this updated CA. The geologic basis for groundwater models has continued to improve with additional data collection and interpretation with the creation and maintenance of the Hanford South Geoframework and the Central Plateau Vadose Zone Geoframework tools; 7. Incorporation of updated tank residual inventory estimates. Tank residual inventory estimates have improved with the incorporation of tank retrieval inventory data for those tanks that have completed retrieval. The CA inventory data package includes this updated information; 8. Incorporation of WMA C PA results. Two additional PAs for tank farm closure decisions are in preparation during the period required to prepare an updated Hanford Site CA: WMA C (FY 2016) and WMA A-AX (in preparation). The WMA C PA results are incorporated into the updated Hanford Site CA, and its grouted residuals model is used as the basis for a release model to account for the other tank farm systems modeling in this updated CA; 9. An update to the IDF PA was submitted in FY 2017 and has been reviewed and approved. The results of the IDF PA are incorporated into this updated CA; 10. Updated risk assessment scenarios. The risk assessment scenarios currently in use for Hanford Site CERCLA and Resource Conservation and Recovery Act of 1976 (RCRA) analyses differ from those evaluated in the initial Hanford Site CA. The representative person exposure scenario evaluated in this updated Hanford Site CA is consistent with recent PAs and CERCLA and RCRA analyses as the Hanford Site. This updated Hanford Site CA provides the following: A comparison of the updated Hanford Site CA all-pathway dose results with the performance measures during the compliance period, which is assumed to begin with site closure in calendar 2070 with the last scheduled disposal action and continue for 1,000 years postclosure (to calendar year 3069); A comparison of Hanford Site CA all-pathway dose results with the performance measures during the postcompliance period to address potential peaks beyond the compliance period. This is accomplished by evaluating dose in the period for 9,000 years following the compliance period (i.e., from calendar year 3070 to 12070).

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2023

Shared micromobility as a first- and last-mile transit solution? Spatiotemporal insights from a novel dataset

The first- and last-mile (FM/LM) problem is a major deterrent to public transit use. With the rise of shared micromobility options such as shared e-scooters in recent years, there is a growing interest in understanding their potential to serve as a last-mile transit solution. However, empirical data regarding the integrated use of shared micromobility and public transit have been limited so far. As a result, much is unknown regarding the spatiotemporal patterns and characteristics of shared micromobility trips serving as an FM/LM connection to transit. Here, this paper addresses these knowledge gaps by leveraging a novel dataset (i.e., the Spin post-ride survey dataset) that records thousands of transit-connecting shared e-scooter trips in Washington DC. Specifically, we used the dataset to reveal the spatiotemporal patterns of transit-connecting shared e-scooter trips in Washington DC, resulting in some major policy insights regarding the integral use of shared e-scooters and public transit. We further leveraged the dataset to validate if and to what extent a commonly applied buffer-zone approach can infer FM/LM micromobility trips accurately. Statistical tests showed that the actual FM/LM Spin e-scooter trips differ from inferred FM/LM Spin e-scooter trips in both spatial and temporal dimensions. This indicates that the common practice of inferring FM/LM micromobility trips with a buffer-zone approach can lead to inaccurate estimates of transit-connecting micromobility trips.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗
DOE OSTI2020

A High-Throughput Solver for Marginalized Graph Kernels on GPU

Here, we present the design and optimization of a solver for efficient and high-throughput computation of the marginalized graph kernel on General Purpose GPUs. The graph kernel is computed using the conjugate gradient method to solve a generalized Laplacian of the tensor product between a pair of graphs. To cope with the large gap between the instruction throughput and the memory bandwidth of the GPUs, our solver forms the graph tensor product on-the-fly without storing it in memory. This is achieved by using threads in a warp cooperatively to stream the adjacency and edge label matrices of individual graphs by small square matrix blocks called tiles, which are then staged in registers and the shared memory for later reuse. Warps across a thread block can further share tiles via the shared memory to increase data reuse. We exploit the sparsity of the graphs hierarchically by storing only non-empty tiles using a coordinate format and nonzero elements within each tile using bitmaps. We propose a new partition-based reordering algorithm for aggregating nonzero elements of the graphs into fewer but denser tiles to further exploit sparsity. We carry out extensive theoretical analyses on the graph tensor product primitives for tiles of various density and evaluate their performance on synthetic and real-world datasets. Our solver delivers three to four orders of magnitude speedup over existing CPU-based solvers such as GraKeL and GraphKernels. The capability of the solver enables kernel-based learning tasks at unprecedented scales.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

The Influence of Soil Properties on Sea-Breeze Circulations in the Southeast U.S.

Sea-breeze circulations (SBCs) are common weather phenomena at and near coastal regions. They form because of a thermal gradient between the land surface at the coast and the sea surface. In a mid-day regime, a “thermal low” generated at the warm coast will lead to rising air motion, creating a wind shift coming from the sea near the surface displacing the coastal air. A “return flow” moving back towards the sea is generated by upper-level divergence because of the rising motion from the thermal low. SBCs propagate and serve as a method of urban pollutant dispersion in the Los Angeles region of California and are constrained to the coast due to the topography of surrounding mountains serving as a boundary for further inland propagation. Within the northeast U.S. SBCs are seen in the warm season but tend to remain coastally bound due to Coriolis distortion over long distances. Within the southeast U.S. (SEUS), the paradigmatic example of SBCs occurs over the Florida peninsula, where thunderstorms form on a nearly daily occurrence due to the convergence of SBCs from the east and west sides of the peninsula. However, there are further examples of sea-breezes in the SEUS that warrant study. Within the region bordering the SEUS and the Mid-Atlantic, just east of the southern Appalachian Mountains, warm-season SBCs form at the coast of Georgia and the Carolinas. Relatively flat topography ~150-200km inland allows for mostly unimpeded inland SBC propagation. Through visual analysis, Viner et al. catalogued several SBCs that propagated as far inland as the Central Savannah River Area surrounding Augusta, Georgia. Wermter et al. found that while the land-sea thermal gradient at the coast can influence coastal SBC genesis, the inland-coastal thermal gradient over the land is the primary influencer on the speed and depth of inland propagation of SBCs in this region. Additionally, soil moisture itself is a known correlative factor to sea-breeze formation, as it influences the soil temperature and the thermal gradient needed for SBC genesis and inland penetration. Physick determined that higher latent heat fluxes associated with wetter soil dampen the land-sea thermal gradient and suppress the formation of a SBC. Physick determined that higher latent heat fluxes associated with wetter soil dampen the land-sea thermal gradient and suppress the formation of a SBC. Conversely, drier soil enhances the thermal gradient and promotes SBC formation. However, while there is an inverse relationship between soil moisture and SBCs, higher soil moisture can actually promote more convective rainfall following a SBC if it does not significantly impact the thermal gradient. While the relationship between soil moisture and SBC formation has been conceptually explored and modeled numerically, there is a research gap in observed connections. The Soil Moisture Active Passive (SMAP) satellite mission has been operational since 2015 and has been used to create high-resolution re-analytical Level 4 (L4) datasets of soil moisture and soil temperature at different soil depths: the surface (0-5cm) and rootzone (0-1m). The surface soil temperature effectively acts as the “skin temperature” of the surface at these levels, and a spatial map of the land-sea as well as the coastal-inland thermal gradients can be represented. SMAP data are also assimilated in some atmospheric models such at the High Resolution Rapid Refresh (HRRR) mesoscale model. We propose leveraging the use of SMAP products to fill in spatial gaps left by weather and mesonet stations within the SEUS region, as well as assessing the effectiveness of utilizing SMAP products towards SBC forecasting in both deterministic and machine learning (ML) models.

58 GEOSCIENCES↗
DOE OSTI2026

Final DOE-ASR Report for the Project “Using LASSO to bridge the gap between model and observations and to learn about atmospheric convection”

Atmospheric convection spans a wide range of spatial and temporal scales and involves complex interactions with the surrounding dynamic and thermodynamic environment, particularly over tropical continental regions. These processes remain a major source of uncertainty in weather and climate models, including persistent biases in the diurnal cycle of convective precipitation that directly affect estimates of climate sensitivity. Addressing these challenges requires the combined use of high-resolution observations and cloud-resolving modeling frameworks. In this context, the DOE Atmospheric Radiation Measurement (ARM) program’s Large-Eddy Simulation ARM Symbiotic Simulation and Observation (LASSO) activity provides a powerful platform that pairs comprehensive observations with numerical simulations to enable process-level understanding of atmospheric convection. Within this context, this Research and Development Partnership Pilot (RDPP) project was designed to initiate and expand DOE ARM/ASR research capacity at minority-serving institutions, while advancing scientific understanding of convective processes over the Amazon rainforest. Consistent with the RDPP mission, the project emphasized partnership development, training, and workforce capacity building alongside exploratory research activities. On the scientific side, the project produced two peer-reviewed journal articles, and one manuscript currently under review (see list in section 3.1). Together, these studies combine long-term ARM observations and cloud-resolving and convection-permitting modeling to investigate the environmental controls on the shallow-to-deep convective transition during the Amazon wet season. The results demonstrate the central role of early-day moisture preconditioning and large-scale dynamical forcing in regulating isolated deep convection, provide mechanistic insight into convective evolution, and establish physically informed modeling frameworks for future sensitivity experiments. These scientific outcomes are described in sections 2.1 to 2.3 and were disseminated in 8 conference presentations (see section 3.2) and 5 invited talks (see section 3.3), reflecting broad engagement with our community. Equally important, the project achieved its RDPP capacity-building objectives (see section 2.4). A sustained research partnership was established among the University of Maryland, Baltimore County (UMBC), Morgan State University (MSU), and Howard University (HU), and extended to include collaboration with Pacific Northwest National Laboratory (PNNL). The project organized multiple multi-day training events focused on ARM data, LASSO simulations, and quantitative analysis methods, directly engaging students, postdoctoral researchers, and faculty across institutions. These activities broadened participation in ASR research and led to independent adoption of LASSO workflows by students beyond the immediate project team. Finally, the project successfully positioned the participating institutions to pursue future DOE research. Preliminary scientific results, coupled with strengthened partnerships and technical capacity, enabled the submission of follow-on proposals to DOE ASR funding opportunities. In this way, the project fulfilled the RDPP goal of seeding durable research capacity and laying the foundation for larger-scale, sustained engagement with DOE ARM and ASR programs.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

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

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

atmospheric warming↗
NASA NTRSDate not supplied

Improved Benchmarking of Cohesive Elements in Abaqus Standard for Predicting Disbond and Delamination in Composite Structures

Traditional approaches for aircraft certification require the assumption of an initial flaw condition, either represented as barely visible impact damage (BVID) or through inclusion of a Teflon insert to serve as surrogate damage. Based on the initial composite damage state, the structure must be shown to demonstrate structural durability and damage tolerance (DaDT) according to the following criteria: a. Damage displays no detrimental growth under cyclic loading b. The structure is able to sustain design limit load (DLL) Currently, the only available manner for validating structural performance is through test. Since damage can occur over a wide variety of areas within a structure, this approach has proven to be increasingly expensive and time consuming for composite airframes and acreage structure within the design-test-certification building block. A further complicating factor is the requirement to accurately capture the most critical damage morphologies as a starting condition. To understand the severity of the damage, it is either required to experimentally determine the most critical areas at tremendous expense or rely on legacy data of similar structural testing, which limits design space expansion. A preferred solution is to use advanced analysis to provide improved understanding of load margins for critical locations based on a wide variety of potential starting damage conditions. The standard industry approach for DaDT certification adheres to the use of the traditional virtual crack closure technique (VCCT) method. VCCT is generally a preferred method because it conforms to the current certification principles of damage from a known flaw, and when used correctly, can be effective at predicting delamination propagation under static and cyclic loading. The VCCT method requires the inclusion of an initial flaw in the finite element (FE) model requiring a-priori knowledge of the flaw location. This in turn requires a plethora of analysis cases to be examined to cover a reasonable span of potential damage states. Additionally, the VCCT approach requires node-to-node connectivity rendering it incompatible with the best practices and approaches for using continuum damage mechanics (CDM) based progressive damage and failure analysis (PDFA) tools within a typical FE solver. Alternatives to VCCT have emerged in the form of cohesive elements which utilize the cohesive zone model (CZM). Unlike VCCT which models linear elastic fracture mechanics, cohesive elements couples continuum and fracture based responses through the use of bilinear traction separation laws. These laws are defined based on a penalty stiffness, a cohesive strength, and a strain energy release rate. The approach can be mesh regularized with native cohesive elements within many FE solvers such as Abaqus and LS-DYNA. In Phase I of the NASA Advanced Composites Consortium (ACC) post-buckled stiffened panel with BVID, Strength and Life [1], the performance of cohesive elements were benchmarked in comparison to VCCT and LEFM solutions and showed good agreement using Abaqus explicit [2]. To realize savings on current and future programs, it is still necessary to close technical gaps related to the use of cohesive elements with Abaqus Standard. Within a program environment, standard finite element analysis is the preferred analytical capability for quasi-static loading as it eliminates uncertainty due to oscillatory behavior commonly seen with explicit analysis. This oscillatory behavior creates difficulties in writing margins of safety based on the analysis. The use of negative tangent stiffness material models complicates convergence which typically requires the use of numerical controls such as viscous damping to overcome. To date, there has not been a comprehensive study on how to establish best practices for cohesive element convergence for predictive capability within the Abaqus implicit solver. In pursuit of these goals, under the NASA ACC program, several numerical benchmark problems were proposed including pure mode I (double cantilevered beam – DCB), pure mode II (end notch flexure – ENF), and symmetric/unsymmetric evolving mixed mode (single leg bend – SLB). This paper focuses on the use of cohesive elements to model the delamination through the use of CZM. Specifically, finite element models for the DCB, ENF, symmetric SLB, and unsymmetric SLB, are developed and various solution controls for convergence are studied to develop a best practice. Once the best practice has been developed, the predictive capability of the objective CZM model is used to analyze the hat pull-off strength of a standard hat stiffened configuration under various loading conditions.

Abaqus↗
DOE OSTI2026

Integrative Modeling and Analysis of Fungal Central Carbon Metabolism

Over a thousand fungal genomes have been sequenced, yet manually curated genome-scale metabolic models (GEMs) are available for only a limited number of species. Moreover, these models have often been developed independently, leading to inconsistencies in namespaces, compartment definitions, and pathway representations that hinder comparative analysis, the systematic reuse of prior curation efforts, and the integration of consolidated metabolic knowledge. Here, we present the Consolidated Fungal Core Metabolism Model (CFCMM), constructed by integrating thirteen published fungal models spanning Ascomycota, Mucoromycota, and both Crabtree-positive and Crabtree-negative yeasts. We harmonized metabolites and reactions into a non-redundant shared ModelSEED ontological space, standardized compartmentalization, and refined gene–protein–reaction (GPR) rules. Using pathway-level visualization and systematic gap detection, we further improved the integrated network through literature-guided curation to correct stoichiometry, stereospecificity, and pathway architecture. Orthologous protein family reconstruction and functional annotation workflows were used to validate and inform GPR associations, with particular emphasis on ambiguous enzyme superfamilies and membrane-associated components. Using the resulting CFCMM, we built high-quality central carbon core models for each fungus and performed flux balance analysis to quantify ATP-yield variation under aerobic and anaerobic conditions, explicitly evaluating scenarios driven by differences in electron transport chain (ETC) composition. Simulations reproduced the expected fermentative yield of approximately 2 mmol ATP per mmol glucose under anaerobic conditions and separated the thirteen fungi into two bioenergetic groups under aerobic respiration based on Complex I status, with predicted yields of approximately 30 versus 22 mmol ATP per mmol glucose. Forcing flux through the alternative oxidase bypass further reduced ATP yields to approximately 12 and 4 mmol ATP per mmol glucose in Complex I-containing and Complex I-lacking fungi, respectively. Collectively, this work provides a manually curated, ModelSEED-consistent, and extensible fungal core metabolic template, deployed in DOE KBase as a resource for automated reconstruction of central carbon core models from any sequenced fungal genome. In addition, the CFCMM provides modular components for developing GEMs with more accurate energy predictions and enables robust comparative analyses of fungal bioenergetics and core metabolic diversity

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2026

Data for Clumping Index Estimation With 30°-tilted Cameras in Row Crops: Evaluation of Methods and Segment Size Effects

The clumping index (CI) quantifies the spatial distribution of foliage elements and is essential for accurately estimating the plant area index (PAI), canopy radiative transfer, and photosynthesis. Traditionally, the finite-length averaging method (LX), the gap size distribution method (CC), and a combined approach of CC and LX (CLX) have been applied to instruments like TRAC and digital hemispherical photography to estimate CI. However, a comprehensive evaluation of these methods in row crops remains limited, especially regarding the influence of segment size on CI. Meanwhile, digital cameras offer a cost-effective and user-friendly solution for canopy measurements in row crops, yet their application in this context remains underexplored. In this study, we employed a new approach using a 30°-tilted digital camera to estimate CI in corn and soybean fields, applying the LX, CC, and CLX methods. We systematically assessed the performance of these three methods by combining field measurements in real-world fields with simulations using the LESS 3D radiative transfer model. Our results showed that CLX applied to the whole image and 45° segment offered accurate estimation of CI (bias within ±0.1, RMSE < 0.2) and PAI (bias within ±0.4, RMSE < 1) in real-world fields and LESS simulations. The accuracy of the LX method was highly sensitive to segment size, with the best performance observed at the 15° segment (PAI bias within ±0.4). In contrast, the CC method remained stable across different segment sizes, and its performance was generally comparable to that of LX, except at the 15° segment. Across view zenith angles, CI derived from CC generally showed a continuous increase, while those from LX and CLX followed a rising trend at small zenith angles but began to decline at 68°, likely due to an increasing proportion of no-gap segments. Seasonally, LX tended to show decreasing CI during early growth stages but increased as the canopy matured, whereas CC and CLX showed gradually increasing CI before plateauing at peak PAI. The 30°-tilted camera effectively captured CI variations across different angles and growth stages, making it a practical and robust instrument for row crop canopy structure analysis. Applying these CI methods to digital cameras offers a low-cost and accessible CI estimation alternative, improving canopy structure monitoring accuracy in row crops.

Modeling↗
NASA NTRS2025

Components Refurbishment and Chemical Analysis Facility, SWMU #041 - Per- and Polyfluoroalkyl Substances Site Assessment Report Kennedy Space Center, Florida

This PFAS Site Assessment (SA) Report presents the activities and results associated with PFAS investigations at the Components Refurbishment and Chemical Analysis (CRCA) facility located at Kennedy Space Center (KSC), Florida. In 2022, CRCA was identified as an Area of Potential Concern because the facility stores several potential PFAS-containing chemicals. A groundwater sample collected from an onsite monitoring well detected PFOA and PFOS at concentrations greater than State of Florida provisional Groundwater Cleanup Target Levels (pGCTLs). PFAS SA activities were conducted between March 2022 and February 2024. During this timeframe, a total of 16 direct-push technology (DPT) locations and 82 discrete samples were collected from these locations, along with 90 monitoring well samples. The analytical data screening process focused on State of Florida pGCTLs and United States Environmental Protection Agency (USEPA) Regional Screening Level (RSLs) from November 2023 to evaluate the data. Analytical results identified PFOS, PFOA, and PFBA at concentrations exceeding their respective RSLs at each depth interval (shallow, intermediate, and deep). PFOS and PFOA had the largest footprint of RSL exceedances in each depth interval, but pGCTL exceedances were only observed at limited locations in the shallow and intermediate intervals. PFBA had the highest detections of any PFAS compound, with concentrations exceeding 180,000 nanograms per liter (ng/L), which is 100-times the RSL of 1,800 ng/L in the shallow and intermediate intervals and 10-times the RSL in the deep interval. The maximum PFBA detection was 841,000 ng/L at shallow monitoring well, MW0006. The PFAS SA also included samples collected from the onsite hydraulic containment system (HCS) which was installed to control and treat the onsite chlorinated volatile organic compound (CVOC) plume. Influent and effluent aqueous samples were collected monthly from the system. PFAS concentrations were relatively the same for both influent and effluent samples, indicating that while the HCS has been effective for CVOC treatment, it does not provide any additional treatment for PFAS compounds. Since the HCS has achieved its objectives, the system was shut down in December 2024. PFAS data gaps still exist, to include surface water and soil, which were not sampled during this SA. Soil sampling near the Chemical Process Area is recommended. Surface water and soil sampling at select stormwater outfalls is also recommended. Additionally, further groundwater sampling is recommended (DPT and monitoring well) in all depth intervals to delineate the extent of PFAS impacts at CRCA.

K Alex Murphy↗
DOE OSTI2021

Demystifying asynchronous I/O Interference in HPC applications

With increasing complexity of HPC workflows, data management services need to perform expensive I/O operations asynchronously in the background, aiming to overlap the I/O with the application runtime. However, this may cause interference due to competition for resources: CPU, memory/network bandwidth. The advent of multi-core architectures has exacerbated this problem, as many I/O operations are issued concurrently, thereby competing not only with the application but also among themselves. Furthermore, the interference patterns can dynamically change as a response to variations in application behavior and I/O subsystems (e.g. multiple users sharing a parallel file system). Without a thorough understanding, I/O operations may perform suboptimally, potentially even worse than in the blocking case. To fill this gap, here we investigate the causes and consequences of interference due to asynchronous I/O on HPC systems. Specifically, we focus on multi-core CPUs and memory bandwidth, isolating the interference due to each resource. Then, we perform an in-depth study to explain the interplay and contention in a variety of resource sharing scenarios such as varying priority and number of background I/O threads and different I/O strategies: sendfile, read/write, mmap/write underlining trade-offs. The insights from this study are important both to enable guided optimizations of existing background I/O, as well as to open new opportunities to design advanced asynchronous I/O strategies.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Connect the Dots: In Situ 4-D Seismic Monitoring of CO 2 Storage With Spatio-Temporal CNNs

4-D seismic imaging has been widely used in CO 2 sequestration projects to monitor the fluid flow in the volumetric subsurface region that is not sampled by wells. Ideally, real-time monitoring and near-future forecasting would provide site operators with great insights to understand the dynamics of the subsurface reservoir and assess any potential risks. However, due to obstacles such as high deployment cost, availability of acquisition equipment, exclusion zones around surface structures, only very sparse seismic imaging data can be obtained during monitoring. That leads to an unavoidable and growing knowledge gap over time. The operator needs to understand the fluid flow throughout the project lifetime and the seismic data are only available at a limited number of times. This is insufficient for understanding reservoir behavior. To overcome those challenges, we have developed spatio-temporal neural-network-based models that can produce high-fidelity interpolated or extrapolated images effectively and efficiently. Specifically, our models are built on an autoencoder, and incorporate the long short-term memory (LSTM) structure with a new loss function regularized by optical flow. We validate the performance of our models using real 4-D post-stack seismic imaging data acquired at the Sleipner CO 2 sequestration field. We employ two different strategies in evaluating our models. Numerically, we compare our models with different baseline approaches using classic pixel-based metrics. We also conduct a blind survey and collect a total of 20 responses from domain experts to evaluate the quality of data generated by our models. Finally, via both numerical and expert evaluation, we conclude that our models can produce high-quality 2-D/3-D seismic imaging data at a reasonable cost, offering the possibility of real-time monitoring or even near-future forecasting of the CO 2 storage reservoir.

4-D seismic imaging↗
NASA NTRS2011

Launch Vehicle Demonstrator Using Shuttle Assets

The Marshall Space Flight Center Advanced Concepts Office (ACO) has the leading role for NASA s preliminary conceptual launch vehicle design and performance analysis. Over the past several years the ACO Earth-to-Orbit Team has evaluated thousands of launch vehicle concept variations for a multitude of studies including agency-wide efforts such as the Exploration Systems Architecture Study (ESAS), Constellation, Heavy Lift Launch Vehicle (HLLV), Heavy Lift Propulsion Technology (HLPT), Human Exploration Framework Team (HEFT), and Space Launch System (SLS). NASA plans to continue human space exploration and space station utilization. Launch vehicles used for heavy lift cargo and crew will be needed. One of the current leading concepts for future heavy lift capability is an inline one and a half stage concept using solid rocket boosters (SRB) and based on current Shuttle technology and elements. Potentially, the quickest and most cost-effective path towards an operational vehicle of this configuration is to make use of a demonstrator vehicle fabricated from existing shuttle assets and relying upon the existing STS launch infrastructure. Such a demonstrator would yield valuable proof-of-concept data and would provide a working test platform allowing for validated systems integration. Using shuttle hardware such as existing RS-25D engines and partial MPS, propellant tanks derived from the External Tank (ET) design and tooling, and four-segment SRB s could reduce the associated upfront development costs and schedule when compared to a concept that would rely on new propulsion technology and engine designs. There are potentially several other additional benefits to this demonstrator concept. Since a concept of this type would be based on man-rated flight proven hardware components, this demonstrator has the potential to evolve into the first iteration of heavy lift crew or cargo and serve as a baseline for block upgrades. This vehicle could also serve as a demonstration and test platform for the Orion Program. Critical spacecraft systems, re-entry and recovery systems, and launch abort systems of Orion could also be demonstrated in early test flights of the launch vehicle demo. Furthermore, an early demonstrator of this type would provide a stop-gap for retaining critical human capital and infrastructure while affording the current emerging generation of young engineers opportunity to work with and capture lessons learned from existing STS program offices and personnel, who were integral in the design and development of the Space Shuttle before these resources are no longer available. The objective of this study is to define candidate launch vehicle demonstration concepts that are based on Space Shuttle assets and determine their performance capabilities and how these demonstration vehicles could evolve to a heavy lift capability to low earth orbit.

Threet, Grady E., Jr.↗
DOE OSTI2020

Climate-Related Trends of Within-Storm Intensities Using Dimensionless Temporal-Storm Distributions

Huff curves are probabilistic time distributions of rainfall expressed as dimensionless cumulative percentages of storm depth and duration. Previous studies have documented development factors, spatial robustness, and the utility of Huff curves in practical applications. However, the effects of trending rainfall on Huff curve intensity patterns have not yet been studied. As such, the goal of this paper is to fill this gap by studying Huff curve patterns in a watershed with demonstrated increasing trends of temperature and precipitation, with the intention that it can be generalized to other areas in the US and the world. To achieve this goal, the high temporal resolution precipitation data collected from a high spatial density, 72-year precipitation-gauge network on the 4.25-km 2 North Appalachian Experimental Watershed in east-central Ohio were used. Seasonal storm pattern trends from 1939 to 2010 were investigated using dimensionless depth (with the frequency of 50%, $d_{50}$) and the curve variability $(V = d_{80} - d_{20})$ at three dimensionless within-storms time periods (three verticals). The Spearman rank correlation procedure (correlation coefficient, ρ and significance probability, p) was used to statistically determine trends over time using 8 periods of 4-season sets of Huff curves over the 72 years. Two averages of ρ and p were computed: (1) by averaging the individual ρ and p obtained from 10 gauges (AvgI); and (2) by grouped averaging of all individual gauge values of $d_{50}$ and V and then computing ρ and p (AvgG). The test results of individual gauges showed that 4 cases for $d_{50}$ and 23 cases for V were significant for all seasons and verticals (total of 120 cases for each variable). The test results of AvgI for $d_{50}$ and V and AvgG for $d_{50}$ showed no significant trends in all seasons and verticals. Only the AvgG for V led to a significant trend for V in spring and fall at different times within storm patterns. The data do not provide sufficient evidence at the p=0.05 significance level to reject the null hypothesis of unchanging position of the dimensionless depth of the 50% Huff curves for individual or averages for all seasons and verticals. Also, there is insufficient evidence to reject the null hypothesis of unchanging variability, V, using the average ρ of individual gauges (AvgI for V) for all seasons and verticals. AvgG results showed a significant trend in V; however, this analysis may be affected by the nonindependence of storm data. The results suggest that it is likely that there is little if any effect of trending climate over approximately 70 years on Huff curve patterns. The results of this study add to the robustness characteristics of Huff curves and to their potential use in hydrological practice as design storms, as the foundation of stochastic storm generation, and for storm disaggregation, and they deserve further investigation. These different forms of inputs to watershed models have the potential to improve runoff estimation. The results suggest that, if verified in other studies, they have applicability to provide useful stationary precipitation patterns across the US and other areas of the world in areas of nonstationary climate. Also, individual rain gauge data may not be representative of trends even over small areas, and seasonal differences were noticeable as found in other studies. Recommendations are provided.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS2017

Ground-Based Observing Campaign of Briz-M Debris

In 2015, NASA's Orbital Debris Program Office (ODPO) completed the installation of the Meter Class Autonomous Telescope (MCAT) on Ascension Island. MCAT is a 1.3m optical telescope designed with a fast tracking capability for observing orbital debris at all orbital regimes (Low-Erath orbits to Geosyncronous (GEO) orbits) from a low latitude site. This new asset is dedicated year-round for debris observations, and its location fills a geographical gap in the Ground-based Electro Optical Space Surveillance (GEODSS) network. A commercial off the shelf (COTS) research grade 0.4m telescope (named the Benbrook telescope) will also be installed on Ascension at the end of 2016. This smaller version is controlled by the same master software, designed by Euclid Research, and can be tasked to work independently or in concert with MCAT. Like MCAT, it has a the same suite of filters, a similar field of view, and a fast-tracking Astelco mount, and is also capable of tracking debris at all orbital regimes. These assets are well suited for targeted campagins or surveys of debris. Since 2013, NASA's ODPO has also had extensive access to the 3.8m infrared UKIRT telescope, located on Mauna Kea. At nearly 14,000-ft, this site affords excellent conditions for collecting both photometery and spectroscopy at near-IR (0.9 - 2.5 micrometers SWIR) and thermal-IR (8 - 25 micrometers; LWIR) regimes, ideal for investigating material properties as well as thermal characteristics and sizes of debris. For the purposes of understanding orbital debris, taking data in both survey mode as well as targeting individual objects for more in-depth characterizations are desired. With the recent break-ups of Briz-M rocket bodies, we have collected a suite of data in the optical, near-infrared, and mid-infrared of in-tact objects as well as those classified as debris. A break-up at GEO of a Briz-M rocket occurred in January, 2016, well timed for the first remote observing survey-campaign with MCAT. Access to the 3.8m UKIRT telescope has also allowed for investigating this break-up in the near-infrared at wavelengths where debris is often much more reflective, allowing for the potential detection of a smaller population of these debris. In addition, a suite of near-IR reflectance spectroscopy (0.8-2.5 micrometers) and thermal-IR (8-15 micrometers) of individual in-tact and debris Briz-M objects has been collected. Analysis of the survey data will be discussed.

Lederer, S. M.↗