Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “High spatial resolution”

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,711 records · Page 95
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

A multiphysics coupling framework for exascale simulation of fracture evolution in subsurface energy applications

Predicting the evolution of fractured media is challenging due to coupled thermal, hydrological, chemical and mechanical processes that occur over a broad range of spatial scales, from the microscopic pore scale to field scale. We present a software framework and scientific workflow that couples the pore scale flow and reactive transport simulator Chombo-Crunch with the field scale geomechanics solver in GEOS to simulate fracture evolution in subsurface fluid-rock systems. This new multiphysics coupling capability comprises several novel features. An HDF5 data schema for coupling fracture positions between the two codes is employed and leverages the coarse resolution of the GEOS mechanics solver which limits the size of data coupled, and is, thus, not taxed by data resulting from the high resolution pore scale Chombo-Crunch solver. The coupling framework requires tracking of both before and after coarse nodal positions in GEOS as well as the resolved embedded boundary in Chombo-Crunch. We accomplished this by developing an approach to geometry generation that tracks the fracture interface between the two different methodologies. The GEOS quadrilateral mesh is converted to triangles which are organized into bins and an accessible tree structure; the nodes are then mapped to the Chombo representation using a continuous signed distance function that determines locations inside, on and outside of the fracture boundary. The GEOS positions are retained in memory on the Chombo-Crunch side of the coupling. The time stepping cadence for coupled multiphysics processes of flow, transport, reactions and mechanics is stable and demonstrates temporal reach to experimental time scales. The approach is validated by demonstration of 9 days of simulated time of a core flood experiment with fracture aperture evolution due to invasion of carbonated brine in wellbore-cement and sandstone. We also demonstrate usage of exascale computing resources by simulating a high resolution version of the validation problem on OLCF Frontier.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2026

Revisiting the Supernova Engines in the 3C 397 and W49B Supernova Remnants

The nature of the supernova remnants (SNRs) 3C 397 and W49B has long been a subject of debate, with prior studies offering conflicting interpretations between thermonuclear and core-collapse scenarios. To help settle this debate, we present a systematic, spatially resolved, spectroscopic analysis of both remnants using XMM-Newton. By applying multicomponent thermal models, we derive key physical properties including elemental abundances, ejecta temperatures, ambient densities, and explosion energetics. We compare the inferred metal abundance ratios to a wide range of core-collapse and thermonuclear nucleosynthesis models, including new models whose explosion energies differ from the canonical value of 10 51 erg. We find that the observed Fe/Si and Ca/Si ratios in both SNRs are best matched by certain thermonuclear models. However, no model fully reproduces the complete set of observed abundance patterns. In 3C 397, high Fe enrichment and spatial abundance variations suggest interaction with a dense progenitor environment, and W49B’s composition is overall consistent with a thermonuclear origin; however, both require a low-energy (∼10 50 erg) supernova explosion. We additionally map the Fe Kα line centroid energies and find a spread, with W49B falling within the core-collapse region—highlighting both environmental complexity and the limitations of this diagnostic for supernova classification. Our results highlight the need for caution in relying on any single diagnostic or nucleosynthesis model for supernova typing, underscore the need for improved nucleosynthesis models, and motivate future high-resolution, high-throughput observations.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2019

A 25 micron-thin microscope for imaging upconverting nanoparticles with NIR-I and NIR-II illumination

Rationale: Intraoperative visualization in small surgical cavities and hard-to-access areas are essential requirements for modern, minimally invasive surgeries and demand significant miniaturization. However, current optical imagers require multiple hard-to-miniaturize components including lenses, filters and optical fibers. These components restrict both the form-factor and maneuverability of these imagers, and imagers largely remain stand-alone devices with centimeter-scale dimensions. Methods: We have engineered INSITE (Immunotargeted Nanoparticle Single-Chip Imaging Technology), which integrates the unique optical properties of lanthanide-based alloyed upconverting nanoparticles (aUCNPs) with the time-resolved imaging of a 25-micron thin CMOS-based (complementary metal oxide semiconductor) imager. We have synthesized core/shell aUCNPs of different compositions and imaged their visible emission with INSITE under either NIR-I and NIR-II photoexcitation. We characterized aUCNP imaging with INSITE across both varying aUCNP composition and 980 nm and 1550 nm excitation wavelengths. To demonstrate clinical experimental validity, we also conducted an intratumoral injection into LNCaP prostate tumors in a male nude mouse that was subsequently excised and imaged with INSITE. Results: Under the low illumination fluences compatible with live animal imaging, we measure aUCNP radiative lifetimes of 600 μs - 1.3 ms, which provides strong signal for time-resolved INSITE imaging. Core/shell NaEr 0.6 Yb 0.4 F 4 aUCNPs show the highest INSITE signal when illuminated at either 980 nm or 1550 nm, with signal from NIR-I excitation about an order of magnitude brighter than from NIR-II excitation. The 55 μm spatial resolution achievable with this approach is demonstrated through imaging of aUCNPs in PDMS (polydimethylsiloxane) micro-wells, showing resolution of micrometer-scale targets with single-pixel precision. INSITE imaging of intratumoral NaEr 0.8 Yb 0.2 F 4 aUCNPs shows a signal-to-background ratio of 9, limited only by photodiode dark current and electronic noise. Conclusion: This work demonstrates INSITE imaging of aUCNPs in tumors, achieving an imaging platform that is thinned to just a 25 μm-thin, planar form-factor, with both NIR-I and NIR-II excitation. Based on a highly paralleled array structure INSITE is scalable, enabling direct coupling with a wide array of surgical and robotic tools for seamless integration with tissue actuation, resection or ablation.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2025

Non-intrusive reduced-order modeling for dynamical systems with spatially localized features

This work presents a non-intrusive reduced-order modeling framework for dynamical systems with spatially localized features characterized by slow singular value decay. The proposed approach builds upon two existing methodologies for reduced and full-order non-intrusive modeling, namely Operator Inference (OpInf) and sparse Full-Order Model (sFOM) inference. We decompose the domain into two complementary subdomains that exhibit fast and slow singular value decay. The dynamics of the subdomain exhibiting slow singular value decay are learned with sFOM while the dynamics with intrinsically low dimensionality on the complementary subdomain are learned with OpInf. The resulting, coupled OpInf-sFOM formulation leverages the computational efficiency of OpInf and the high resolution of sFOM, and thus enables fast non-intrusive predictions for conditions beyond those sampled in the training data set. A novel regularization technique with a closed-form solution based on the Gershgorin disk theorem is introduced to promote stable sFOM and OpInf models. We also provide a data-driven indicator for subdomain selection and ensure solution smoothness over the interface via a post-processing interpolation step. We evaluate the efficiency of the approach in terms of offline and online speedup through a quantitative, parametric computational cost analysis. We demonstrate the coupled OpInf-sFOM formulation for two test cases: a one-dimensional Burgers’ model for which accurate predictions beyond the span of the training snapshots are presented, and a two-dimensional parametric model for the Pine Island Glacier ice thickness dynamics, for which the OpInf-sFOM model achieves an average prediction error on the order of 1% with an online speedup factor of approximately 8$\times$ compared to the numerical simulation.

42 ENGINEERING↗
DOE OSTI2023

Comparing the cooling effectiveness of operationalisable urban surface combination scenarios for summer heat mitigation

Extreme summer heat in cities exacerbates the vulnerability of urban communities to heatwaves. Vegetative and reflective urban surfaces can help reduce urban heat. This study investigated the impacts of urban trees, green roofs and cool roofs on heat mitigation during average and extreme summer conditions in temperate oceanic Melbourne, Australia. We simulated the city climate using ‘The Air Pollution Model’ (TAPM) at a 1 km spatial resolution over 10 years, which according to our review of the literature, was the most prolonged period for simulation in Melbourne. During a widespread heatwave event, some of the tested scenarios with combined surface parameters could reduce the extreme values of the energy budget components- sensible heat, latent heat, and storage heat fluxes up to seasonal averages compared to the existing situation for Melbourne (control). The scenario with the highest (reasonable maximum) ground-level vegetation, green roofs, and cool roofs could reduce air temperatures up to 2.4 °C. The simulations suggest that a combined strategy with vegetative and high-albedo surfaces will deliver higher effectiveness with maximum cooling benefits and cost-effectiveness than individual strategies in cities. These results suggest the importance of collaborative strategic planning of urban surfaces to make cities healthier, sustainable, and liveable.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

A conserved epitope III on hepatitis C virus E2 protein has alternate conformations facilitating cell binding or virus neutralization

Epitope III, a highly conserved amino acid motif of 524 APTYSW 529 on the hepatitis C virus (HCV) E2 glycoprotein, resides in the critical loop that binds to the host receptor CD81, thus making it one of the most important antibody targets for blocking HCV infections. Here, we have determined the X-ray crystal structure of epitope III at a 2.0-Å resolution when it was captured by a site-specific neutralizing antibody, monoclonal antibody 1H8 (mAb1H8). The snapshot of this complex revealed that epitope III has a relatively rigid structure when confined in the binding grooves of mAb1H8, which confers the residue specificity at both ends of the epitope. Such a high shape complementarity is reminiscent of the “lock and key” mode of action, which is reinforced by the incompatibility of an antibody binding with an epitope bearing specific mutations. By subtly positioning the side chains on the three residues of Tyr 527 , Ser 528 , and Trp 529 while preserving the spatial rigidity of the rest, epitope III in this cocrystal complex adopts a unique conformation that is different from previously described E2 structures. With further analyses of molecular docking and phage display–based peptide interactions, we recognized that it is the arrangements of two separate sets of residues within epitope III that create these discrete conformations for the epitope to interact selectively with either mAb1H8 or CD81. These observations thus raise the possibility that local epitope III conformational dynamics, in conjunction with sequence variations, may act as a regulatory mechanism to coordinate “mAb1H8-like” antibody-mediated immune defenses with CD81-initiated HCV infections.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

Photoionization of the I 4d and valence orbitals of methyl iodide

Here, the photoabsorption and photoionization dynamics of the I 4d and valence orbitals in methyl iodide have been studied both experimentally and theoretically. Synchrotron radiation has been employed to measure the total ion yield in the vicinity of the I 4d ionization thresholds. The observed structure, due to excitations into Rydberg or valence states, has been assigned using transition energies and relative intensities computed with time-dependent density functional theory within the Tamm–Dancoff approximation. Photoelectron spectra, recorded with plane polarized radiation in two polarization geometries, have allowed the effect of autoionization on the valence electron angular distributions to be investigated. The spectra obtained at photon energies of 50.62 and 52.34 eV, coinciding respectively with the I 4d 5/2 → σ* and 4d 3/2 → σ* transitions, reveal, in addition to valence shell photoelectron bands, features not associated with simple photoionization of the parent molecule. High resolution photoelectron spectra of the I 4d main-lines display structure resulting from spin–orbit coupling and molecular field splitting. The binding energies of the five states contributing to the (I 4d) -1 ionization have been determined. The iodine (in CH 3 I) N 45 VV Auger spectrum has been measured and the observed structure has been assigned using the core hole binding energies derived in the present work together with established ionization energies of the doubly charged ion. The experimentally determined Auger electron angular distributions have been discussed in relation to the theoretical angular distribution parameter characterizing the spatial alignment of molecular axes in the (I 4d) -1 state.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2024

Assessing the Atmospheric Response to Subgrid Surface Heterogeneity in the Single-Column Community Earth System Model, Version 2 (CESM2)

Land-atmosphere interactions are central to the evolution of the atmospheric boundary layer and the subsequent formation of clouds and precipitation. Existing global climate models represent these connections with bulk approximations on coarse spatial scales, but observations suggest that small-scale variations in surface characteristics and co-located turbulent and momentum fluxes can significantly impact the atmosphere. Recent model development efforts have attempted to capture this phenomenon by coupling existing representations of subgrid-scale (SGS) heterogeneity between land and atmosphere models. Such approaches are in their infancy and it is not yet clear if they can produce a realistic atmospheric response to surface heterogeneity. Here, we implement a parameterization to capture the effects of SGS heterogeneity in the Community Earth System Model (CESM2), and compare single-column simulations against high-resolution Weather Research and Forecasting (WRF) large-eddy simulations (LESs), which we use as a proxy for observations. The CESM2 experiments increase the temperature and humidity variances in the lowest atmospheric levels, but the response is weaker than in WRF-LES. In part, this is attributed to an underestimate of surface heterogeneity in the land model due to a lack of SGS meteorology, a separation between deep and shallow convection schemes in the atmosphere, and a lack of explicitly represented mesoscale secondary circulations. These results highlight the complex processes involved in capturing the effects of SGS heterogeneity and suggest the need for parameterizations that communicate their influence not only at the surface but also vertically.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Novel Proppant Logging Technique for Infill Drilling of Unconventional Shale Wells

Summary During the development of an unconventional play, wells are drilled and completed in batches, and depending on the development plans, current and expected energy market trends, as well as other developmental considerations, new wells are drilled and hydraulically fractured later near existing producing laterals. This creates challenges in terms of optimizing resource recovery and reducing interwell communication. A novel approach is proposed that utilizes systematic composite sampling and analysis of drilling mud returns to look for and quantitatively identify sand particles. The workflow involves cleaning, drying, and segregation of samples into sizes of interest to us (size distribution of pumped proppant in offset parent wells). These samples are imaged at a very high resolution and analyzed for grains using characteristic optical imaging properties to classify proppant sand particles using computer vision algorithms. Further analysis, such as elemental compositional analysis, is used to validate the results from the imaging workflow. We present a case study from the Permian Basin, where a new child well was used as a test case to prove this technology at the Hydraulic Fracturing Test Site (HFTS-2) in Delaware Basin. We introduce new proppant parameters that help identify sustained proppant zones vs. localized propped fractures. We have used additional diagnostics and data collected at the test site to validate observations from the proppant log and have successfully interpreted significantly propped vs. unpropped zones. A key finding from this test has been the significant proppant transport distances observed away from parent wells. Observable proppant was found at a lateral distance of approximately 425 m for one set of parent wells and more than 915 m for another set of parent wells. While a major limitation of this technique is the sampling rate, given adequate sampling, the proposed technology represents a systematic and one-of-a-kind interpretation of spatial proppant distribution while drilling infill wells. It provides us with unique opportunities to better understand the current state of the reservoir being targeted, including zones that are likely highly drained relative to others, and how the planned hydraulic fracturing of child wells can be improved.

Energy & Fuels↗
DOE OSTI2025

Decoding crops one cell at a time: from cell atlases to single-cell genetics

Understanding the mechanisms underlying key agricultural traits remains a central challenge in crop research, but recent advances in technologies are providing powerful tools to address this issue. Among these, single-cell and spatial transcriptomics have revealed tissue heterogeneity and spatial organization, offering unique insights into cellular gene expression dynamics and the coordinated activity of multiple cell types. These approaches help uncover how specific cell types contribute to agricultural traits and refine candidate loci lists through integration with trait-associated loci. Additionally, single-cell and spatial transcriptomics have the potential to serve as cell-level readout platforms integrating cellular perturbations, enabling high-throughput discovery of causal relationships between genotype and gene expression at the cellular level in plants. Successful implementation will accelerate the identification of key genetic variants for crop improvement. Furthermore we review lessons learned from application of single-cell screening in mammalian cells, highlight major technical and biological barriers to its use in plants, and outline potential strategies to overcome these challenges. Together, the widespread application and integration of single-cell and spatial transcriptomics with other technologies enable not only the descriptive cataloging of cell states but also the causal interrogation of sequence functions and regulatory networks at cell type resolution, ultimately advancing gene function studies and accelerating crop improvement.

Cellular heterogeneity↗
DOE OSTI2019

Compressed video camera with a moving patterned disk

A high-speed video system is disclosed that includes a moving image absorbing disk at an image plane. The disk has a pattern that passes and blocks image data. The disk is located between an event and an image sensor, or reflects an image to the image sensor. The disk is rotated at a speed that matches the desired reconstructed image frame rate. The image sensor frame data is processed using image reconstruction techniques, such as the D-AMP or TWIST algorithm, to recover a time sequence of reconstructed images. Additional images can be reconstructed for each image sensor frame if some spatial resolution is sacrificed. For continuous video, the disk speed is adjusted to the sensor frame rate. For burst mode, a single sensor image is acquired and a short image sequence is reconstructed. This image capture system works with a variety of radiations, including infrared, light, UV and X-rays.

42 ENGINEERING↗
DOE OSTI2025

Self-Supervised T-GCN for Detection of Disturbance and Propagation in Power Grid

Urban power systems increasingly rely on dense sensing to monitor grid reliability, yet disturbance labels are scarce and events are rare. We present a self-supervised spatio-temporal method that detects, localizes, and characterizes grid frequency disturbances across urban areas using only unlabeled data. Our approach trains a tiny Temporal Graph Convolutional Network (T-GCN) to forecast per-site frequency residuals (deviation from 60 Hz). The sensor graph is constructed directly from signals using pre-event Pearson correlation with a cross-correlation lag penalty without geocoding. At inference, node-level anomalies are the model's forecast errors; region-level alarms arise from connected components of high-score nodes. We estimate disturbance propagation by computing per-node arrival times (first persistent exceedance), then fit a planar or time-of-arrival model to obtain direction, speed, and an epicenter proxy. With only three real events collected at decisecond resolution across U.S. cities, we evaluate the T-GCN and report time-to-detect, footprint size, and propagation consistency. We further show that short-window embeddings from the T-GCN's hidden states enable few-shot event-vs-background recognition via a simple prototypical classifier. Despite minimal data and no labels, our system yields fast, spatially coherent detection and interpretable propagation maps, offering a practical, lightweight pathway to city-scale grid resilience analytics.

Niu, Haoran [ORNL] (ORCID:0000000155228297)↗
DOE OSTI2019

Optimization of SiV Defect Yield in Diamond Substrates

Color centers (defect complexes such as SiV) in diamond have shown potential in fields ranging from metrology, cybersecurity to quantum computation. Demonstrations in these fields have pushed the envelope of state-of-the-art operations – for example, single photon sources (SPS) making use of SiV centers in diamond for quantum key distribution have demonstrated all the requirements for SPS operation including: (1) stable operation with second correlation function <<1, (2) electrically driven single photon emission and (3) compatibility with frequenc y conversion to telecommunication frequencies. To-date, however, all these demonstrations have been on lab-scale one-off devices. The key question behind how to deterministically fabricate these devices, namely activation yield has been overlooked. For context, Si based semiconductor devices are hugely successful because we have a high activation yield for implanted dopants. This is not yet true for diamond color centers. As currently understood, the color center yield is dominated by a lack of vacancies in the immediate area of the implantation. We propose to optimize the activation yield of color center using a combination of (1) focused single ion implantation with in-situ detection to count the number of implanted Si ions and (2) localized point defect (vacancy) creation using a focused Li ion beam to improve the yield. These experiments build on the unique capabilities of the SNL nanoImplanter (nI) to produce focused ion beam with spatial resolution of < 10 nm of both Si and Li ions. This work will also leverage our world-leading single ion implantation and detection capabilities.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Measurement of the Neutron Elastic Electromagnetic Form Factor Ratio at Large Momentum Transfer

Exploring nucleon structure is vital both for understanding its origin and existence as well as for the advancement of the sciences. It helps us answer key questions such as how quark and gluon dynamics create 99% of the nucleon mass. Electron- nucleon scattering has been widely used for precision studies of the nucleon and nuclear structure since the Nobel Prize winning investigations by Robert Hofstadter and collaborators in the 1950s. These studies provide information about the spatial charge and current densities of the nucleon in terms of the electromagnetic form factors. The form factors are functions of four momentum transfer squared (Q2). Extending the electromagnetic form factor measurements to higher Q2 plays a critical role in furthering the understanding of nucleon structure. This motivated the Super BigBite Spectrometer (SBS) program at Jefferson Lab. The open nature of the spectrometers and the direct line of sight from the target to the tracking detector locations in experimental setups such as SBS creates high levels of background at the detectors. This necessitates the use of tracking detectors with high rate capability and good position resolution. Gas Electron Multiplier (GEM) detectors are an excellent choice for tracking detectors in such experiments. Understanding the performance of the GEM detectors is important not just for SBS experiments but also for future high-luminosity experiments. This thesis reports the exploratory results from the measurement of the neutron elastic electromagnetic form factor ratio (Gn E/Gn M) at high momentum transfer. A longitudinally polarized electron beam was scattered off a polarized 3He target, used as an effective polarized neutron target. In this experiment, the polarized 3He target achieved a world record polarization weighted luminosity at a beam current of 45 µA. Double spin asymmetry of the scattered neutron events is used to extract the neutron form factor ratio. Measurements were taken at Q2 = 3.0, 6.8, 9.8 (GeV/c)2. The lowest Q2 measurement is in good agreement with the existing world data, and the higher-Q2 measurements extend the Q2 reach well beyond the existing world data and are expected to remain unmatched for a long time.

Gamage, Vimukthi Haththotuwa [Univ. of Virginia, C↗
DOE OSTI2019

Spatial Heterodyne Raman Spectrometer (SHRS) for In Situ Chemical Sensing Using Sapphire and Silica Optical Fiber Raman Probes

A spatial heterodyne Raman spectrometer (SHRS), constructed using a modular optical cage and lens tube system, is described for use with a commercial silica and a custom single-crystal (SC) sapphire fiber Raman probe. The utility of these fiber-coupled SHRS chemical sensors is demonstrated using 532 nm laser excitation for acquiring Raman measurements of solid (sulfur) and liquid (cyclohexane) Raman standards as well as real-world, plastic-bonded explosives (PBX) comprising 1,3,5- triamino- 2,4,6- trinitrobenzene (TATB) and octahydro-1,3,5,7-tetranitro-1,3,5,7-tetrazocine (HMX) energetic materials. The SHRS is a fixed grating-based dispersive interferometer equipped with an array detector. Each Raman spectrum was extracted from its corresponding fringe image (i.e., interferogram) using a Fourier transform method. Raman measurements were acquired with the SHRS Littrow wavelength set at the laser excitation wavelength over a spectral range of ∼1750 cm −1 with a spectral resolution of ∼8 cm −1 for sapphire and ∼10 cm −1 for silica fiber probes. The large aperture of the SHRS allows much larger fiber diameters to be used without degrading spectral resolution as demonstrated with the larger sapphire collection fiber diameter (330 μm) compared to the silica fiber (100 μm). Unlike the dual silica fiber Raman probe, the dual sapphire fiber Raman probe did not include filtering at the fiber probe tip nearest the sample. Even so, SC sapphire fiber probe measurements produced less background than silica fibers allowing Raman measurements as close as ∼85 cm −1 to the excitation laser. Despite the short lengths of sapphire fiber used to construct the sapphire probe, well-defined, sharp sapphire Raman bands at 420, 580, and 750 cm −1 were observed in the SHRS spectra of cyclohexane and the highly fluorescent HMX-based PBX. SHRS measurements of the latter produced low background interference in the extracted Raman spectrum because the broad band fluorescence (i.e., a direct current, or DC, component) does not contribute to the interferogram intensity (i.e., the alternating current, or AC, component). SHRS spectral resolution, throughput, and signal-to-noise ratio are also discussed along with the merits of using sapphire Raman bands as internal performance references and as internal wavelength calibration standards in Raman measurements.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2023

Mapping Glacier Basal Sliding Applying Machine Learning

During the RESOLVE project (“High-resolution imaging in subsurface geophysics: development of a multi-instrument platform for interdisciplinary research”), continuous surface displacement and seismic array observations were obtained on Glacier d’Argentière in the French Alps for 35 days in May 2018. The data set is used to perform a detailed study of targeted processes within the highly dynamic cryospheric environment. In particular, the physical processes controlling glacial basal motion are poorly understood and remain challenging to observe directly. Especially in the Alpine region for temperate based glaciers where the ice rapidly responds to changing climatic conditions and thus, processes are strongly intermittent in time and heterogeneous in space. Spatially dense seismic and Global Positioning System (GPS) measurements are analyzed applying machine learning to gain insight into the processes controlling glacial motions of Glacier d’Argentière. Using multiple bandpass-filtered copies of the continuous seismic waveforms, we compute energy-based features, develop a matched field beamforming catalog and include meteorological observations. Features describing the data are analyzed with a gradient boosting decision tree model to directly estimate the GPS displacements from the seismic noise. We posit that features of the seismic noise provide direct access to the dominant parameters that drive displacement on the highly variable and unsteady surface of the glacier. The machine learning model infers daily fluctuations and longer term trends. The results show on-ice displacement rates are strongly modulated by activity at the base of the glacier. The techniques presented provide a new approach to study glacial basal sliding and discover its full complexity.

58 GEOSCIENCES↗
DOE OSTI2022

Photosynthesis phenology, as defined by solar-induced chlorophyll fluorescence, is overestimated by vegetation indices in the extratropical Northern Hemisphere

Vegetation phenology is highly sensitive to climate change, although the data and methods used to estimate key phenological states can influence this sensitivity. Because of its direct relation to leaf photosynthetic carbon uptake, remotely sensed solar-induced chlorophyll fluorescence (SIF) can provide new insight assessing changes in vegetation phenology. In this work, we investigated the potential of using a SIF time series product named contiguous SIF (CSIF) to estimate spring, summer, and autumn phenology in the extratropical Northern Hemisphere (>30°N) and compared the results with those based on Moderate Resolution Imaging Spectroradiometer (MODIS) Normalized Difference Vegetation Index (NDVI) for the period 2001–2018. Overall, we found similar spatial patterns in phenological states. However, specific dates of key phenological events differed when using CSIF vs. MODIS NDVI data. NDVI data indicated that the growing season started earlier (by 10.1 days on average) and ended later (11.5 days on average) relative to CSIF data. This implies that actual periods for photosynthetic activity are shorter (by 21.6 days on average) than those estimated from vegetation indices more directly related to changes in canopy structure. These large differences between results from NDVI and that from CSIF suggest that vegetation indices such as NDVI seem to overestimate the period for active photosynthesis over the extratropical Northern Hemisphere. Furthermore, while phenology of the early growing season is dominated by temperature for both NDVI and CSIF data, phenology of the late growing season is mainly controlled by temperature for NDVI but by precipitation for CSIF. Our findings were further confirmed by other SIF (GOME-2 SIF) and vegetation index (MODIS EVI) datasets. Phenology modes in Earth system modelling are often parameterized using leaf unfolding and senescence from either station or satellite observations. Our results imply that canopy structure-based parameterization schemes may have overestimated photosynthesis active period, and thus productivity responses. We conclude that SIF data offers a novel and unique approach for assessing phenological change - one that is more directly tied to the carbon cycle and how it is being influenced by climate change.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

CERF: IM3 Projected Western US Power Plant Locations

Overview The Capacity Expansion Regional Feasibility (CERF) model is an open-source geospatial python package that provides new power plant locations at a 1km resolution. The model ingests U.S. state or regional-scale electricity system capacity expansion plans, such as those produced by the Global Change Analysis Model (GCAM-USA), and identifies feasible, site-specific locations for individual new power plants (renewable and non-renewable). CERF combines high-resolution geospatial suitability analyses with an economic algorithm that selects individual plant siting locations based on grid interconnection costs and the locational marginal value of new generation. The model incorporates a wide range of dynamic constraints and opportunities, such as protected lands, population density, existing infrastructure, and water availability. This dataset provides CERF power plant siting results for IM3 Phase 2 simulations across eight different scenarios for the Western US through 2055. The scenarios include combinations of two Shared Socioeconomic Pathways (SSP3 and SSP5) with four high-resolution climate projections specific to the United States (see, https://tgw-data.msdlive.org/). These climate projections include "hotter" and "cooler" variants for two Representative Concentration Pathways (RCP4.5 and RCP8.5). The resulting eight simulations are: rcp45cooler_ssp3 rcp45cooler_ssp5 rcp45hotter_ssp3 rcp45hotter_ssp5 rcp85cooler_ssp3 rcp85cooler_ssp5 rcp85hotter_ssp3 rcp85hotter_ssp5 CERF siting results in this dataset correspond to capacity expansion plans in the GCAM-USA IM3 Phase 2 simulation data and are available for each of the above scenarios. Data Details Temporal Range: 2015-2055 in 5-year timesteps. Note that 2015 is the experiment base year and 2020 and beyond represent model simulation years. Spatial Range: Plant locations are provided for the eleven states in the Western US including Arizona, California, Colorado, Idaho, Montana, New Mexico, Nevada, Oregon, Utah, Washington, and Wyoming. Spatial Resolution: 1 km-squared, provided in x and y coordinates Geospatial Projection: Albers Equal Area Conic (ESRI:102003) File Type: csv The dataset contains subdirectories for each of the eight scenarios described in the overview. Each scenario folder contains two subfolders with the following information: 1. Power Plant Data This directory contains a single .csv file of power plant locations for both pre-existing (non-CERF sited plants in operation in 2015) and new (CERF-sited) power plants across the temporal range along with additional CERF model output parameters for CERF-sited plants. Plant with a siting year earlier than 2020 correspond to facilities that are operational leading into the first timestep CERF simulation. For a more detailed description of CERF model output parameters, see the CERF model documentation. Note that the cerf_plant_id parameter is unique within each scenario file but not across scenario files. Parameter Descriptions scenario - Name of scenario cerf_plant_id - Unique siting identifier cerf_sited - If True, indicates that plant was sited by CERF model. If False, indicates pre-existing facility region_name - Name of region (state) tech_id - Technology ID tech_name - Full generation technology name inclusive of cooling type (if applicable) and additional characteristics tech_simple - Simplified generation technology type unit_size_mw - Power plant unit size (MW) xcoord - X coordinate in the default CRS (meters) ycoord - Y coordinate in the default CRS (meters) index - Index position in the flattend 2D array buffer_in_km - Exclusion buffer around site (km) sited_year - Year of siting retirement_year - Year of retirement lmp_zone - Locational marginal price (LMP) zone ID locational_marginal_price_usd_per_mwh - Locational marginal price ($/MWh) generation_mwh_per_year - Generation output (MWh/yr) operating_cost_usd_per_year - Cost of plant operations ($/yr) net_operational_value - Net operational value based on LMP and and operating costs ($/yr) interconnection_cost - Cost of interconnection for transmission & gas pipeline (if applicable) net_locational_cost -- Difference of interconnection cost and operating value ($/yr) capacity_factor_fraction - Capacity factor (fraction) carbon_capture_rate_fraction - Carbon capture rate (fraction) fuel_co2_content_tons_per_btu - Fuel CO2 content (tons/Btu) fuel_price_usd_per_mmbtu - Fuel price ($/MMBtu) fuel_price_esc_rate_fraction - Fuel price escalation rate (fraction) heat_rate_btu_per_kWh - Heat rate (Btu/kWh) lifetime_yrs - Technology lifetime for annuity (years) operational_life_yrs - Operational lifetime for retirement (years) variable_om_usd_per_mwh - Variable operation and maintenance costs of yearly capacity use ($/MWh) variable_om_esc_rate_fraction - Variable operation and maintenance costs escalation rate (fraction) carbon_tax_usd_per_ton - Carbon tax ($/ton) carbon_tax_esc_rate_fraction - Carbon tax escalation rate (fraction) 2. Storage Data This directory contains information on new and pre-existing energy storage facilities operational in each timestep along with various storage operational parameters. The 2015 timestep provides pre-existing energy storage data and corresponds with facilities that are operational leading into the first model simulation timestep. Note that coordinates in the storage files correspond to the interconnection point on the grid (substation location), not individual energy storage locations. Energy storage is added in a cumulative process at each given interconnection point. That is, each individual file provides the total operational storage capacity interconnected to the specified substation for the given timestep, inclusive of previously installed storage at that location and new storage installed in that timestep at that location. Parameters scenario - Name of scenario timestep - Simulation timestep name - Unique storage identifier s_typ - Type of energy storage technology (battery or pumped storage hydro) s_node - Node ID of interconnecting substation xcoord - X coordinate in the default CRS (meters) ycoord - Y coordinate in the default CRS (meters) charge_rate - Maximum charge rate (power capacity) of storage system (MW) discharge_rate - Maximum discharge rate (power capacity) of storage system (MW) duration - Duration of storage system (hours) max_SoC - Allowed maximum state of charge (energy capacity) of storage system (MWh) min_SoC -Allowed minimum state of charge (energy capacity) of storage system (MWh) charge_eff - Efficiency of charge (fraction between 0 and 1) discharge_eff - Efficiency of discharge (fraction between 0 and 1) Acknowledgment IM3 is a multi-institutional effort led by Pacific Northwest National Laboratory and supported by the U.S. Department of Energy's Office of Science as part of research in MultiSector Dynamics, Earth and Environmental Systems Modeling Program.

CERF↗