Search NASA⌕ Search

SEARCH · Search NASA

Results for “pose estimation”

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 289 records · Page 16
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

An automated and portable method for selecting an optimal GPU frequency

Power consumption poses a significant challenge in current and emerging graphics processing unit (GPU) enabled high-performance computing systems. In modern GPUs, dynamic voltage frequency scaling (DVFS) appears to be a reliable control to regulate power consumption and performance. However, the DVFS design space is large - hence, brute-force approaches are infeasible to select the optimal frequency. Furthermore, no single frequency can be universally optimal for applications with varying computational intensities. Thus, the application's complexity and the availability of a wide range of frequency settings are a challenge in selecting the optimal frequency configuration for a given GPU workload. To that end, this paper proposes a systematic approach that consists of three steps. The feature characterization study identifies the fine-grain GPU utilization metrics that influence the power consumption and execution time of a given workload. To understand the performance, power, and energy consumption behaviors of a workload across GPU's DVFS design space, we derived analytical power and performance models using the identified fine-grain features. Here, it is shown that the same set of GPU utilization metrics can estimate both the power consumption and execution time while being agnostic of changes to frequency and input sizes. Applying a power control with the single objective of reducing power may cause performance degradation, leading to more energy consumption. A multi-objective approach is proposed to select the optimal GPU DVFS configuration for a workload that reduces power consumption with negligible degradation in performance. The evaluation was conducted using SPEC ACCEL benchmarks and three real applications - NAMD LAMMPS, and LSTM on NVIDIA GV100, GA100, and AMD MI210 GPUs. On average, real applications showed 29.6% energy savings with a performance loss of 5.2% on GA100 and 22.6% energy savings with a performance loss of 4.7% on GV100. Moreover, the proposed models are portable to real applications, GPU architectures, and vendors, and require metric collection at only the default frequency rather than all supported DVFS configurations. Additionally, we conducted a comparison between our models and the GPU assembly instructions (PTX)-based static models. The results revealed a significant reduction in the average error rates, with a decrease from 19.7% to 3.1% for power models and from 29.4% to 5.2% for performance models.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Underwater unexploded ordnance discrimination based on intrinsic target polarizabilities – A case study

Seabed unexploded ordnance that resulted partly from the high failure rate among munitions from more than 80 years ago and from decades of military training and testing of weapons systems poses an increasing concern all around the world. Although existing magnetic systems can detect clusters of debris, they are not able to tell whether a munition is still intact requiring special removal (e.g. in situ detonation) or is harmless scrap metal. The marine environment poses unique challenges, and transferring knowledge and approaches from land to a marine environment has not been easy and straightforward. On land, the background soil conductivity is much lower than the conductivity of the unexploded ordnance and the electromagnetic response of a target is essentially the same as that in free space. For those frequencies required for target characterization in the marine environment, the seawater response must be accounted for and removed from the measurements. The system developed for this study uses fields from three orthogonal transmitters to illuminate the target and four three-component receivers to measure the signal arranged in a configuration that inherently cancels the system's response due to the enclosing seawater, the sea–bottom interface and the air–sea interface for shallow deployments. The system was tested as a cued system on land and underwater in San Francisco Bay – it was mounted on a simple platform on top of a support structure that extended 1 m below and allowed the diver to place metal objects to a specific location even in low-visibility conditions. The measurements were stable and repeatable. Furthermore, target responses estimated from marine measurements matched those from land acquisition, confirming that the seawater and air–sea interface responses were removed successfully. Thirty-six channels of normalized induction responses were used for the classification, which was done by estimating the target principal dipole polarizabilities. Our results demonstrated that the system can resolve the intrinsic polarizabilities of the target, with clear distinctions between those of symmetric intact unexploded ordnance and irregular scrap metal. The prototype system was able to classify an object based on its size, shape and metal content and correctly estimate its location and orientation.

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗
DOE OSTI2024

Carbon storage cost modeling for the offshore Gulf of Mexico

Groundbreaking for geologic carbon storage (GCS) projects in the offshore Gulf of Mexico is imminent, and there is great interest in utilizing this region for GCS projects. Offshore saline reservoirs provide a significant and accessible resource for GCS. However, conducting GCS in the offshore environment will pose distinct challenges pertaining to site selection, operations, infrastructure use, and monitoring compared to operating onshore that ultimately affect technoeconomic assessment of offshore GCS projects. Carbon storage and transport costs are critical to project developers looking to deploy carbon storage in the offshore environment. We present CO2_S_COM_Offshore, a model developed by the National Energy Technology Laboratory (NETL) as a screening-level offshore saline GCS cost modeling tool. Based on NETL’s widely used CO2_S_COM cost model for onshore saline CS, CO2_S_COM_Offshore enables technoeconomic analysis of GCS in offshore areas. This model comprehensively incorporates multiple facets of offshore GCS projects, from regional evaluation and site selection to permitting, transport, operations, monitoring, site closure, and decommissioning. In general, the model can explore the cost implications for potential offshore GCS project(s) by enabling the user to change several project operational and financial attribute configurations. Key inputs include offshore storage formation options, CO2 injection rate and duration, infrastructure types, monitoring intensity, project financing, and post-injection site care duration. Supporting cost algorithms within CO2_S_COM_Offshore were compiled utilizing S&P Global’ s QUE$TORTM cost estimation software alongside a variety of open-source scientific literature. In addition to reviewing key model components, we discuss several sensitivity analyses, input variabilities, and results on analysis of break-even CO2 price required by a project based on different regulation/policy and operational scenarios for the offshore Gulf of Mexico. These results indicate the value of modeling offshore GCS specifically, and the potential of offshore GCS within a decarbonization value chain. Presented at the 41st USAEE/IAEE North American Conference, 3-6 November 2024, Baton Rouge, LA, United States.

Mark-Moser, Mackenzie K.↗
DOE OSTI2024

Dynamic Cooling on Contemporary Quantum Computers

We study the problem of dynamic cooling whereby a target qubit is cooled at the expense of heating up N − 1 further identical qubits by means of a global unitary operation. A standard back-of-the-envelope high-temperature estimate establishes that the target qubit temperature can be dynamically cooled by at most a factor of 1 / N . Here we provide the exact expression for the minimum temperature to which the target qubit can be cooled and reveal that there is a crossover from the high initial temperature regime, where the scaling is 1 / N , to a low initial temperature regime, where a much faster scaling of 1 / N occurs. This slow, 1 / N scaling, which was relevant for early high-temperature NMR quantum computers, is the reason dynamic cooling was dismissed as ineffectual around 20 years ago; the fact that current low-temperature quantum computers fall in the fast, 1 / N scaling regime, reinstates the appeal of dynamic cooling today. We further show that the associated work cost of cooling is exponentially more advantageous in the low-temperature regime. We discuss the implementation of dynamic cooling in terms of quantum circuits and examine the effects of hardware noise. We successfully demonstrate dynamic cooling in a three-qubit system on a real quantum processor. Since the circuit size grows quickly with N , scaling dynamic cooling to larger systems on noisy devices poses a challenge. We therefore propose a suboptimal cooling algorithm, whereby relinquishing a small amount of cooling capability results in a drastically reduced circuit complexity, greatly facilitating the implementation of dynamic cooling on near-future quantum computers. Published by the American Physical Society2024

Physics↗
DOE OSTI2025

A derecho climatology (2004–2021) in the United States based on machine learning identification of bow echoes

Due to their persistent widespread severe winds, derechos pose significant threats to human safety and property, with impacts comparable to many tornadoes and hurricanes. Yet, automated detection of derechos remains challenging due to the absence of spatiotemporally continuous observations and the complex criteria employed to define the phenomenon. This study presents an objective derecho detection approach capable of automatically identifying derechos through both observations and model results. The approach is grounded in a physically based definition of derechos and integrates three algorithms: (1) the Python Flexible Object Tracker (PyFLEXTRKR) algorithm to track mesoscale convective systems (MCSs), (2) a semantic segmentation convolutional neural network to identify bow echoes, and (3) a comprehensive classification algorithm to detect derechos within MCS life cycles and distinguish derecho-producing from non-derecho-producing MCSs. Using this approach, we developed a novel high-resolution (4 km and hourly) observational dataset of derechos and accompanying derecho-producing MCSs over the United States east of the Rocky Mountains from 2004 to 2021. The dataset consists of two subsets based on different gust speed data sources and is analyzed to document the climatology of derechos in the United States. On average, 12–15 derechos are identified per year, aligning with previous estimations (∼6–21 events annually). The spatial distribution and seasonal variation patterns are consistent with prior studies, showing peak occurrences in the Great Plains and the Midwest during the warm season. Additionally, during the study period, derechos account for approximately 3.1 % of measured damaging gusts (≥25.93 m s−1) over the eastern United States. The dataset is publicly available at https://doi.org/10.5281/zenodo.14835362 (Li et al., 2025).

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Assessment of oil and gas fields in California as potential CO 2 storage sites

California's total annual greenhouse gas (GHG) emissions (425.3 MtCO 2 e) in 2018 were about 6.4% of the US total (6,677 MtCO 2 e) and around 1% of global emissions. About 39% of 2018 GHG emissions in California were from the industrial and electrical sectors. Many of these emissions were from large stationary point sources and were suitable for carbon capture retrofit with subsequent storage of the captured carbon dioxide (CO 2 ) in geological formations. Previous studies of California found suitable geology and CO 2 storage resource. This study refines and furthers prior work using a three-stage screening process of oil fields, gas fields, and underground natural gas storage (UGS) sites by combining criteria from previous studies while excluding sites that pose technical risk or are located in regions with surface restrictions including sensitive habitats and dense populations. In the first stage, 129 CO 2 storage sites in California were identified using qualification criteria based upon formation properties including geological conditions and pore pressure. The second stage identified sensitive sites by applying conservative screens including seismic activity, faulting, population density, restricted lands, and sensitive habitats. During the third stage, 61 CO 2 potential storage sites were identified by subtraction of stage 2 areas from stage 1. The potential storage volume in the third stage ranged from 1.0 to 2.0 GtCO 2 . Finally, we applied a scoring system with seven parameters to rank the 61 potential sites based on subsurface technical criteria. The scored sites are classified as high priority, medium priority, and sites for future study. Prospective CO 2 storage sites with high and moderate priority were selected and linked to CO 2 sources. There are 14 prospective sites (above 20 MtCO 2 storage resource per site) with a total storage resource of 1024 MtCO 2 distributed in Northern and Southern California. Of these sites, there are 9 potential CO 2 -EOR sites and 1 depleted oil field with a total estimated CO 2 storage volume of ~800 MtCO 2 in the Southern San Joaquin and Ventura Basin. These 10 prospective sites with a storage resource greater than 20 MtCO 2 could potentially deliver more than 20 years of storage with an average injection rate of 40 MtCO 2 /year. The remaining 4 highly prospective sites are in Northern California. Additionally, study results suggest that saline formations should be re-evaluated in concert with storage in oil, gas, and natural gas storage reservoirs.

58 GEOSCIENCES↗
DOE OSTI2025

Evaluation and calibration of MERRA-2 and CAMS reanalysis for PM 2.5 in a semi-urbanized area in the south of the Amazon

Air pollution has significant implications for the climate and poses irreversible risks to human health. The Amazon region of Brazil is severely affected by biomass burning (BB) emissions, yet air quality monitoring remains highly inadequate. Given the scarcity of surface-based observations, reanalysis models have become essential tools for assessing air pollution. Although MERRA-2 and CAMS PM 2.5 products are widely utilized, their validation and comprehensive evaluation for the Amazon Basin remain limited. Here, this study assesses the performance of these products in a semi-urbanized region in the southern Amazon. The calibrated time series was employed to analyze PM 2.5 concentrations from 2003 to 2023. Our results showed satisfactory performance of both products for the 24-h averages of PM 2.5 , with linear correlations above 0.76. However, it was found that both products overestimate surface concentrations. MERRA-2 performed better, with approximately 30% lower bias than CAMS. Time series analysis showed that the study area is strongly impacted by emissions BB in the dry period, mainly in August and September. Furthermore, our findings indicate a positive trend in increasing PM 2.5 concentrations, with a notable rise observed since 2014. The average PM 2.5 levels frequently exceed the daily air quality guidelines established by the WHO in 2021. It has been estimated that the population of this region is exposed to concentrations above 15 μg.m -3 , on average, more than 30 days per year. Our results contribute to the evaluation of MERRA-2 and CAMS products for Amazon and provide a corrected estimate for surface PM 2.5 . Recent concerns about air quality and the implementation of new surface monitoring networks may improve the evaluation of reanalysis products. In the short term, the need for this information makes our assessments indispensable.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Universal dwell time optimization for deterministic optics fabrication

Computer-Controlled Optical Surfacing (CCOS) has been greatly developed and widely used for precision optical fabrication in the past three decades. It relies on robust dwell time solutions to determine how long the polishing tools must dwell at certain points over the surfaces to achieve the expected forms. However, as dwell time calculations are modeled as ill-posed deconvolution, it is always non-trivial to reach a reliable solution that 1) is non-negative, since CCOS systems are not capable of adding materials, 2) minimizes the residual in the clear aperture 3) minimizes the total dwell time to guarantee the stability and efficiency of CCOS processes, 4) can be flexibly adapted to different tool paths, 5) the parameter tuning of the algorithm is simple, and 6) the computational cost is reasonable. In this study, we propose a novel Universal Dwell time Optimization (UDO) model that universally satisfies these criteria. First, the matrix-based discretization of the convolutional polishing model is employed so that dwell time can be flexibly calculated for arbitrary dwell points. Second, UDO simplifies the inverse deconvolution as a forward scalar optimization for the first time, which drastically increases the solution stability and the computational efficiency. Finally, the dwell time solution is improved by a robust iterative refinement and a total dwell time reduction scheme. The superiority and general applicability of the proposed algorithm are verified on the simulations of different CCOS processes. A real application of UDO in improving a synchrotron X-ray mirror using Ion Beam Figuring (IBF) is then demonstrated. The simulation indicates that the estimated residual in the 92.3 mm × 15.7 mm CA can be reduced from 6.32 nm Root Mean Square (RMS) to 0.20 nm RMS in 3.37 min. After one IBF process, the measured residual in the CA converges to 0.19 nm RMS, which coincides with the simulation.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Synthetic Evaluation of Infrasonic Multipole Waveform Inversion

Acoustic source inversions estimate the mass flow rate of volcanic explosions or yield of chemical explosions and provide insight into potential source directionality. However, the limitations of applying these methods to complex sources and their ability to resolve a stable solution have not been investigated in detail. We perform synthetic infrasound waveform inversions that use 3-D Green’s functions for a variety of idealized and realistic deployment scenarios using both a flat plane and Yasur volcano, Vanuatu as examples. We investigate the ability of various scenarios to retrieve the input source functions and relative amplitudes for monopole and multipole (monopole and dipole) inversions. Infrasound waveform inversions appear to be a robust method to quantify mass flow rates from simple sources (monopole) using deployments of infrasound sensors placed around a source, but care should be taken when analyzing and interpreting results from more complex acoustic sources (multipole) that have significant directional components. In the examples we consider the solution is stable for monopole inversions with a signal-to-noise ratio greater than five and the dipole component is small. For most scenarios investigated, the vertical dipole component of the multipole explosion source is poorly constrained and can impact the ability to recover the other source term components. Because multipole inversions are ill-posed for many deployments, a low residual does not necessarily mean the proper source vector has been recovered. Synthetic studies are able to help investigate the limitations and place bounds on information that may be missing using monopole and multipole inversions for potentially directional sources.

58 GEOSCIENCES↗
DOE OSTI2021

Modeling 233 Pa Generation in Thorium-fueled Reactors for Safeguards

Thorium has been considered as a possible alternative to uranium for nuclear fuel for many decades. It is three to four times more abundant in the earth than uranium and produces significantly less long-lived transuranic nuclear waste. Some claim thorium poses fewer proliferation concerns than other fuel types largely due to 232 U buildup (and associated high energy gamma-emitting decay products) in the irradiated thorium fuel. However, to fully explore potential proliferation concerns, generation and subsequent decay of 233 Pa produced in the reactor core still must be studied. With its half-life of 27 days, 233 Pa decays to 233 U, which is an International Atomic Energy Agency (IAEA) defined special fissionable material that can be used for nuclear weapons production. With more research being dedicated to thorium-fueled reactors, and several of these reactor designs possessing online fuel processing (allowing for on-site protactinium separation), it is important to understand this potential proliferation pathway. In particular, it is theoretically possible to extract protactinium from the irradiated fuel salt before it decays into 233 U. This hypothetical potential diversion can become an even greater proliferation concern if the extracted protactinium is purified through a second separation of protactinium approximately ten days later to remove the short half-life decay products of 232 Pa and 234 Pa, thus resulting in a higher concentration of the 233 Pa isotope, which decays into weapons usable 233 U with hardly any 232 U or 234 U in it. To estimate the concern of this potential proliferation challenge of thorium, different nuclear material accountancy techniques were reviewed for their viability to quantify 233 Pa if extracted from used thorium fuel. Characteristics of interest included technology maturity, cost, precision, and time taken to acquire results. Some technologies, like hybrid K-edge densitometry and passive gamma spectroscopy, appear to be viable techniques based on current literature. Due to the limited scope of this project, only passive gamma spectroscopy was further investigated. Three different reactor types (PWR, CANDU, MSR) were modeled with mixed thorium-uranium oxide fuels that were burned until the fuel was spent. The protactinium in the used fuel was extracted at the time of shutdown and the change in isotopic content of the protactinium quantified. Gamma spectroscopy simulations were performed for the protactinium isotopes and their decay products at various decay times. Given the simplicity of the models and large assumptions made (e.g. no background, no shielding, no self-attenuation), the initial results indicate that though 233 Pa is detectible for all the reactor types modeled at all decay times (0 to 300 days), more work should be done with higher fidelity models.

07 ISOTOPE AND RADIATION SOURCES↗
DOE OSTI2025

Modeling of hepatitis B virus infection spread in primary human hepatocytes

ABSTRACT Chronic hepatitis B virus (HBV) infection poses a significant global health threat, causing severe liver diseases including cirrhosis and hepatocellular carcinoma. We characterized HBV DNA kinetics in primary human hepatocytes (PHHs) over 32 days post-inoculation (p.i.) and modified ourin-vivoagent-based modeling (ABM) to gain insights into the HBV lifecycle and spreadin vitro. Parallel PHH cultures were mock-treated or treated with HBV entry inhibitor Myr-preS1 (6.25 µg/mL) was initiated 24 h p.i. In untreated PHH, three viral DNA kinetic patterns were identified: (i) an initial decline, followed by (ii) rapid amplification and (iii) slower amplification/accumulation. In the presence of Myr-preS1, viral DNA and infected cell numbers in phase 3 were effectively blocked, with minimal to no increase. This suggests that phase 2 represents viral amplification in initially infected cells, while phase 3 corresponds to viral spread to naïve cells. The ABM reproduced well the HBV kinetic patterns observed and predicted that the viral eclipse phase lasts between 18 and 38 h. After the eclipse phase, the viral production rate increased over time, starting with a slow production cycle of 1 virion per day, which gradually accelerated to 1 virion per hour after 3 days. Approximately 4 days later, virion production reached a steady state production rate of 4 virions/h. The estimated median efficacy of Myr-preS1 in blocking HBV spread was 91% (range: 90–92%). The HBV kinetics and the predicted estimates of the HBV eclipse phase duration and HBV production cycles in PHH are similar to those predicted in uPA/SCID mice with human livers. IMPORTANCE While primary human hepatocytes (PHHs) are the most physiologically relevant culture system for studying HBV infectionin vitro, a comprehensive understanding of HBV infection kinetics and spread in PHH is lacking. In this study, we characterize HBV viral kinetics and modify ourin vivoagent-based modeling (ABM) to provide quantitative insights into the HBV production cycle and viral spread in PHH. The ABM provides an estimate of the HBV eclipse phase duration, HBV production cycles, and Myr-preS1 efficacy in blocking HBV spread in PHH. The results resemble those predicted in uPA/SCID mice with human livers, demonstrating that estimated HBV infection kinetic parameters in PHHin vitromirror those observed in thein vivoHBV infection chimeric mouse model.

Virology↗
DOE OSTI2025

Co-Firing Switchgrass and Waste Coal in A Power Plant: A Techno-Economic and Life Cycle Evaluation for The Ohio River Valley (SWITCH) (Final Technical Report for Ohio State/FE0032204)

Abandoned coal mine lands (AMLs) represent one of the most persistent environmental challenges in the United States. Prior to the enactment of the Surface Mining Control and Reclamation Act (SMCRA) in 1977, coal mining operations were not legally required to reclaim disturbed lands, leaving behind approximately 500,000 AML sites nationwide. These sites pose severe environmental and health risks, including acid mine drainage, soil and water contamination, and spontaneous combustion of waste coal piles. Millions of Americans live within one mile of these AMLs, underscoring the urgency of remediation. Traditional reclamation practices, such as planting cool-season grasses, often fail to fully restore ecological function or leverage the economic potential of these lands. This project addressed these challenges by developing integrated strategies for resource recovery, land reclamation, and sustainable energy production. This project evaluated an integrated strategy to convert this liability into an opportunity by recovering waste coal and co-firing it with switchgrass (Panicum virgatum L.) cultivated on reclaimed or marginal AML areas in existing coal-fired power plants. Switchgrass not only provides a renewable feedstock but also aids in land reclamation and carbon sequestration. 1) Remote Sensing and Machine Learning for Waste Coal Identification Using Sentinel-2 satellite imagery and supervised classification, we applied four machine learning models to detect historical waste coal piles. Random Forest achieved the highest accuracy (precision: 86%, recall: 77%). Time-series analysis revealed gradual vegetation recovery since 1986, indicating natural reclamation processes in historical sites, while active mining areas showed ongoing disturbance. This workflow enables scalable monitoring and prioritization of reclamation efforts. 2) UAS-Based Stockpile Volume Estimation To quantify recoverable waste coal, we evaluated Unmanned Aerial Systems (UAS) equipped with Light Detection and Ranging (LiDAR) and multispectral sensors. Structure-from-Motion (SfM) photogrammetry combined with interpolated Digital Terrain Models (DTMs) achieved strong agreement with LiDAR reference volumes (Root Mean Square Error (RMSE) ≈147 m 3 , Mean Absolute Percentage Error (MAPE) ≈2%). Sensitivity analysis confirmed that spatial resolution significantly influences accuracy, emphasizing the need for high-resolution data for precise volume estimation. This approach offers a scalable, cost-effective, and accurate alternative to conventional ground-based surveys. 3) Switchgrass Cultivation for Bioenergy and Water Quality Improvement We assessed the hydrological and water quality impacts of converting AMLs to switchgrass production areas using the Soil and Water Assessment Tool (SWAT). Results showed that converting 10% of the watershed area into the switchgrass production zone reduced streamflow by 3.1%, total suspended solids by 18.1%, total nitrogen by 7.6%, and total phosphorus by 6.2%, while achieving biomass yields of 8.6–9.2 metric tons per hectare. These findings highlight switchgrass as a dual-benefit strategy for land reclamation and bioenergy feedstock production. 4) Integrated Co-Firing and CCS for Carbon-Negative Power Generation We modeled co-firing scenarios using the Power Plant Flexible Model (PPFM) to evaluate plant efficiency, greenhouse gas (GHG) emissions, and levelized cost of electricity (LCOE). Without carbon capture and storage (CCS), increasing switchgrass co-firing ratios reduced LCOE from $\$$150/MWh at 0% biomass to $\$$110/MWh at full substitution. Under CCS, costs remained higher (~$\$$250/MWh at 0% biomass) but decreased to $\$$200/MWh at 100% biomass, while enabling net-zero or carbon-negative electricity due to switchgrass sequestration benefits. Although CCS introduced efficiency penalties, pairing it with biomass co-firing offset these impacts and maximized climate benefits. Overall, optimizing co-firing ratios between 60-100%, supported by reliable logistics and storage strategies, emerged as a practical pathway to balance affordability, sustainability, and net-zero or negative GHG emissions while promoting productive reuse of AMLs.

01 COAL, LIGNITE, AND PEAT↗
DOE OSTI2020

Characterization of Partially Observed Epidemics - Application to COVID-19

This report documents a statistical method for the "real-time" characterization of partially observed epidemics. Observations consist of daily counts of symptomatic patients, diagnosed with the disease. Characterization, in this context, refers to estimation of epidemiological parameters that can be used to provide short-term forecasts of the ongoing epidemic, as well as to provide gross information for the time-dependent infection rate. The characterization problem is formulated as a Bayesian inverse problem, and is predicated on a model for the distribution of the incubation period. The model parameters are estimated as distributions using a Markov Chain Monte Carlo (MCMC) method, thus quantifying the uncertainty in the estimates. The method is applied to the COVID-19 pandemic of 2020, using data at the country, provincial (e.g., states) and regional (e.g. county) levels. The epidemiological model includes a stochastic component due to uncertainties in the incubation period. This model-form uncertainty is accommodated by a pseudo-marginal Metropolis-Hastings MCMC sampler, which produces posterior distributions that reflect this uncertainty. We approximate the discrepancy between the data and the epidemiological model using Gaussian and negative binomial error models; the latter was motivated by the over-dispersed count data. For small daily counts we find the performance of the calibrated models to be similar for the two error models. For large daily counts the negative-binomial approximation is numerically unstable unlike the Gaussian error model. Application of the model at the country level (for the United States, Germany, Italy, etc.) generally provided accurate forecasts, as the data consisted of large counts which suppressed the day-to-day variations in the observations. Further, the bulk of the data is sourced over the duration before the relaxation of the curbs on population mixing, and is not confounded by any discernible country-wide second wave of infections. At the state-level, where reporting was poor or which evinced few infections (e.g., New Mexico), the variance in the data posed some, though not insurmountable, difficulties, and forecasts were able to capture the data with large uncertainty bounds. The method was found to be sufficiently sensitive to discern the flattening of the infection and epidemic curve due to shelter-in-place orders after around 90% quantile for the incubation distribution (about 10 days for COVID-19). The proposed model was also used at a regional level to compare the forecasts for the central and north-west regions of New Mexico. Modeling the data for these regions illustrated different disease spread dynamics captured by the model. While in the central region the daily counts peaked in the late April, in the north-west region the ramp-up continued for approximately three more weeks.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2025

Machine learning and process-based modeling of spatiotemporal changes in active layer thickness across Alaska

Permafrost degradation poses a growing threat to infrastructure stability and ecosystem resilience in the rapidly warming Arctic. We investigated the spatiotemporal dynamics of active layer thickness (ALT) across Alaska by integrating field observations, environmental datasets, a physically based Stefan model, and machine learning (ML) techniques. Using weather projections from the Coupled Model Intercomparison Project Phase 6 under two Shared Socioeconomic Pathways (SSP 2-4.5 and SSP 5-8.5), we assessed ALT sensitivity to projected future weather conditions. The random forest (RF) model outperformed the Stefan approach in predicting ALT on the training dataset (R² = 0.84 vs. 0.53) but demonstrated lower generalizability on the test dataset (R² = 0.24 vs. 0.54). The root mean square error (RMSE) for the RF model for training and testing ranged from 14 to 22 cm, compared to 17 and 18 cm for the Stefan model. Variable importance analysis revealed that mean annual temperature and slope angle were the strongest predictors of ALT, accounting for 19% and 18% of the variance, respectively, followed by sediment transport index (14%) and stream power index (11%). Comparative analysis of baseline ALT predictions showed the Stefan model tended to project a thicker active layer (mean ± SD: 65 ± 16 cm), compared to the RF model (mean ± SD: 59 ± 8.8) cm). Both models indicated a latitudinal gradient in ALT, with shallower depths at higher latitudes. Projected ALT increases by 2100 were estimated at 3.3 ± 2.2 cm under SSP 2-4.5 and 5.9 ± 4.0 cm under SSP 5-8.5 for the ML model, whereas the Stefan model projected substantially larger increases of 13 ± 2.6 cm (SSP 2-4.5) and 28 ± 4.4 cm (SSP5-8.5). Spatial analysis showed the greatest ALT increases in northern Alaska, with relatively smaller changes in southern regions. These findings highlight the complex, multifactorial nature of ALT dynamics and the value of hybrid modeling approaches. As rising temperatures accelerate permafrost thaw, changes in ALT can disrupt ecosystems, damage infrastructures, and enhance the release of stored soil carbon, highlighting the urgent need for improved predictive capabilities to inform adaptation strategies in the Arctic.

Climate sciences↗
DOE OSTI2024

Modeling the impact of extreme summer drought on conventional and renewable generation capacity: Methods and a case study on the Eastern U.S. power system

Across recent years, there has been a growing prevalence of extreme weather events throughout the United States, posing significant challenges to the reliable and resilient operation of power systems. Specifically, summer droughts threaten to severely reduce available generation capacity to meet regional electricity demand, potentially leading to power outages. This underscores the importance of accurate resource adequacy (RA) assessment to ensure the reliable operation of the nation’s energy infrastructure. Accurately evaluating the usable capacity of regional generation fleets is a challenging undertaking due to the intricate interactions between power systems and hydro-climatic systems. Here, this paper proposes a systematic and analytical framework to evaluate the impacts of extreme summer drought events on the available capacity of various generating technologies, incorporating both meteorological and hydrologic factors. The framework provides detailed plant-level capacity derating models for hydroelectric, thermoelectric, and renewable power plants, facilitating evaluations with high temporal and spatial resolution. The application of the proposed impact assessment framework to the 2025 generation fleet of the real-world power system within the PJM and SERC regions of the United States yields insightful results. By analyzing the daily usable capacity of 6,055 at-risk generators across the study region, it shows that the summer capacity deration is most significant for hydroelectric and once-through thermal power plants, followed by recirculating thermal power plants and combustion turbines. In the event of the recurrence of the 2007 southeastern summer drought event in the near future, the generation fleet could experience a substantial reduction in available capacity, estimated at approximately 8.5 GW, compared to typical summer conditions. The sensitivity analysis reveals that the usable capacity of the generation fleet would suffer an even more significant decrease under conditions of increasingly severe summer droughts. The proposed approach and the findings of this study provide valuable methodologies and insights, empowering stakeholders to bolster the resilience of power systems against the potentially devastating effects of future extreme drought events.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2022

Integrating Deep Learning and Hydrodynamic Modeling to Improve the Great Lakes Forecast

The Laurentian Great Lakes, one of the world’s largest surface freshwater systems, pose a modeling challenge in seasonal forecast and climate projection. While physics-based hydrodynamic modeling is a fundamental approach, improving the forecast accuracy remains critical. In recent years, machine learning (ML) has quickly emerged in geoscience applications, but its application to the Great Lakes hydrodynamic prediction is still in its early stages. This work is the first one to explore a deep learning approach to predicting spatiotemporal distributions of the lake surface temperature (LST) in the Great Lakes. Our study shows that the Long Short-Term Memory (LSTM) neural network, trained with the limited data from hypothetical monitoring networks, can provide consistent and robust performance. The LSTM prediction captured the LST spatiotemporal variabilities across the five Great Lakes well, suggesting an effective and efficient way for monitoring network design in assisting the ML-based forecast. Furthermore, we employed an explainable artificial intelligence (XAI) technique named SHapley Additive exPlanations (SHAP) to uncover how the features impact the LSTM prediction. Our XAI analysis shows air temperature is the most influential feature for predicting LST in the trained LSTM. The relatively large bias in the LSTM prediction during the spring and fall was associated with substantial heterogeneity of air temperature during the two seasons. In contrast, the physics-based hydrodynamic model performed better in spring and fall yet exhibited relatively large biases during the summer stratification period. Finally, we developed a statistical integration of the hydrodynamic modeling and deep learning results based on the Best Linear Unbiased Estimator (BLUE). The integration further enhanced prediction accuracy, suggesting its potential for next-generation Great Lakes forecast systems.

Xue, Pengfei (ORCID:000000025702421X)↗
DOE OSTI2023

Exploring the scaling limitations of the variational quantum eigensolver with the bond dissociation of hydride diatomic molecules

Abstract Materials simulations involving strongly correlated electrons pose fundamental challenges to state‐of‐the‐art electronic structure methods but are hypothesized to be the ideal use case for quantum computing algorithms. To date, no quantum computer has simulated a molecule of a size and complexity relevant to real‐world applications, despite the fact that the variational quantum eigensolver (VQE) algorithm can predict chemically accurate total energies. Nevertheless, because of the many applications of moderately sized, strongly correlated systems, such as molecular catalysts, the successful use of the VQE stands as an important waypoint in the advancement toward useful chemical modeling on near‐term quantum processors. In this paper, we take a significant step in this direction. We lay out the steps, write, and run parallel code for an (emulated) quantum computer to compute the bond dissociation curves of the TiH, LiH, NaH, and KH diatomic hydride molecules using the VQE. TiH was chosen as a relatively simple chemical system that incorporates d orbitals and strong electron correlation. Because current VQE implementations on existing quantum hardware are limited by qubit error rates, the number of qubits available, and the allowable gate depth, recent studies using it have focused on chemical systems involving s and p block elements. Through VQE + UCCSD calculations of TiH, we evaluate the near‐term feasibility of modeling a molecule with d‐orbitals on real quantum hardware. We demonstrate that the inclusion of d‐orbitals and the use of the UCCSD ansatz, which are both necessary to capture the correct TiH physics, dramatically increase the cost of this problem. We estimate the approximate error rates necessary to model TiH on current quantum computing hardware using VQE + UCCSD and show them to likely be prohibitive until significant improvements in hardware and error correction algorithms are available.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

tih_vqe [SWR-23-32]

This software supports the paper, "Exploring the scaling limitations of the variational quantum eigensolver with the bond dissociation of hydride diatomic molecules," published in the International Journal of Quantum Chemistry, whose abstract is as follows: Materials simulations involving strongly correlated electrons pose fundamental challenges to state-of-the-art electronic structure methods but are hypothesized to be the ideal use case for quantum computing. To date, no quantum computer has simulated a molecule of a size and complexity relevant to real-world applications, despite the fact that the variational quantum eigensolver (VQE) algorithm can predict chemically accurate total energies. Nevertheless, because of the many applications of moderately-sized, strongly correlated systems, such as molecular catalysts, the successful use of the VQE stands as an important waypoint in the advancement toward useful chemical modeling on near-term quantum processors. In this paper, we take a significant step in this direction. We lay out the steps, write, and run parallel code for an (emulated) quantum computer to compute the bond dissociation curves of the TiH, LiH, NaH, and KH diatomic hydride molecules using VQE. TiH was chosen as a relatively simple chemical system that incorporates d orbitals and strong electron correlation. Because current VQE implementations on existing quantum hardware are limited by qubit error rates, the number of qubits available, and the allowable gate depth, recent studies have focused on chemical systems involving s and p block elements. Through VQE + UCCSD calculations of TiH, we evaluate the near-term feasibility of modeling a molecule with d-orbitals on real quantum hardware. We demonstrate that the inclusion of d-orbitals and the use of the UCCSD ansatz, which are both necessary to capture the correct TiH physics, dramatically increase the cost of this problem. We estimate the approximate error rates necessary to model TiH on current quantum computing hardware using VQE+UCCSD and show them to likely be prohibitive until significant improvements in hardware and error correction algorithms are available.

Graf, Peter↗