Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “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 685 records · Page 38
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Strain-Modified Raman Responses in Monolayer MoS 2 Nanobubbles Resolved at 5 nm

The formation of nanoscale bubbles is an unavoidable consequence during the transfer of two-dimensional materials onto target substrates, driven by van der Waals interactions at the interface. While often viewed as imperfections, these nanoscale bubbles have garnered considerable scientific interest due to the substantial in-plane strain gradients they induce, which in turn give rise to a variety of intriguing optoelectronic effects, particularly in semiconducting transition metal dichalcogenides. Determining and analyzing the strain distribution within nanobubbles at the nanoscale is crucial for advancing our understanding of these underlying strain-induced effects. Here, we present a high-resolution scanning tunneling microscopy-based tip-enhanced Raman spectroscopic investigation of localized nanoscale strain distribution within the nanobubbles formed between monolayer MoS 2 and Au interface. By employing cryogenic temperature (78 K), we successfully differentiate the nanoscale Raman signatures between nanobubble edge and pristine MoS 2 . We verify a maximum tensile strain of ∽1.15–1.34% at the nanobubble edge, which gradually diminishes toward the center, yielding a cross-sectional strain profile consistent with a doughnut-shaped distribution. Furthermore, we report to achieve ∽5 nm spatial resolution in probing such edge-localized strain within the nanobubble. In addition, comparative average strain analysis of such MoS 2 nanobubbles is conducted via geometric mechanistic analysis such as membrane and nonlinear plate theories, providing key insight into the geometric nature near the bubble edge. Our work provide fundamental information about strain-induced nanoscale chemical understanding of 2D materials on the nanometer scale, paving the way for practical applications of nanobubbles in strain-engineered optoelectronic devices.

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

iButton and Tinytag snow/ground interface temperature measurements at Teller 27 and Kougarok 64 from 2022-2023, Seward Peninsula, Alaska

Snow/ground interface temperature measurements were collected at the NGEE Arctic Teller Road Site at mile marker 27 (TL_MM27) and at the Kougarok Road Site at mile marker 64 (KG_MM64) on the Seward Peninsula, Alaska. Data were collected between October 1, 2022 to September 18, 2023 using iButton Link DS1921G-F5# Thermochron miniature temperature sensors (https://www.ibuttonlink.com/products/ds1921g) and Tinytag TGP-4017 internal sensors (https://www.micronmeters.com/product/tgp-4017-internal-sensor-40-to-85-c-40-f-to-185-f) deployed across the Kougarok and Teller sites. These sensors are a cost-efficient way to collect snowpack temperatures at a higher spatial resolution than what is normally achieved. iButton data were collected every 4 hours, while Tinytag data were collected every 30 minutes. In total, data were collected from 196 iButtons and 26 Tinytags. This dataset contains four *.csv files of near-ground surface temperatures at various locations throughout each study site and two *.kml files of sensor locations. Data were collected throughout the snow cover season so that snowpack characteristics could be derived using the temperature data. Sensors were placed both inside and outside of vegetation to better capture the spatial variability of snow properties across each domain.The Next-Generation Ecosystem Experiments: Arctic (NGEE Arctic) was a 15-year research effort (2012-2027) to reduce uncertainty in Earth System Models by developing a predictive understanding of carbon-rich Arctic ecosystems and feedbacks to climate. NGEE Arctic was supported by the Department of Energy’s Office of Biological and Environmental Research.The NGEE Arctic project had two field research sites: 1) located within the Arctic polygonal tundra coastal region on the Barrow Environmental Observatory (BEO) and the North Slope near Utqiagvik (Barrow), Alaska and 2) multiple areas on the discontinuous permafrost region of the Seward Peninsula north of Nome, Alaska.Through observations, experiments, and synthesis with existing datasets, NGEE Arctic provided an enhanced knowledge base for multi-scale modeling and contributed to improved process representation at global pan-Arctic scales within the Department of Energy’s Earth system Model (the Energy Exascale Earth System Model, or E3SM), and specifically within the E3SM Land Model component (ELM).

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Estimating snow cover from high-resolution satellite imagery by thresholding blue wavelengths: Supporting Data

The extent and duration of snow cover is predicted to be altered as the climate changes. Developing high-resolution estimates of snow cover change is crucial for estimating changes in snow cover and the effects of these changes on watershed and ecosystems processes. Remote sensing tools have been a common method for rapidly mapping snow covered area (SCA) across a landscape. The most common remote sensing method for estimating SCA uses satellite-based calculations of the normalized difference snow index (NDSI), which relies on spectral measurements in the shortwave-infrared wavelengths (SWIR). NDSI is effective at catchment- to regional-scale estimates of SCA, but due to spatial resolution limitations of SWIR measurements, NDSI cannot be used to assess fine-scale SCA. In this work, we develop a new algorithm, called the Blue Snow Threshold (BST) algorithm, that maps high-resolution SCA by calculating a threshold on the blue wavelengths from high-resolution satellite imagery. This data package includes Orthorectified IKONOS-2 imagery (IkonosTestImage.tif) from August 14, 2004 at 1.00 meters Ground Sample Distance for Cook Inlet, Alaska (59.966414 , -152.982975). The Blue Snow Threshold algorithm (BST.py) was then used to produce a snow cover estimate (IkonosTestImage_BST.tif) for this study area. Additional imagery metadata is included in the ImageInfo.txt file. See Thaler et al., 2023 (https://doi.org/10.1016/j.rse.2022.113403) for more information about the BST algorithm.The Next-Generation Ecosystem Experiments: Arctic (NGEE Arctic), was a 15-year research effort (2012-2027) to reduce uncertainty in Earth System Models by developing a predictive understanding of carbon-rich Arctic ecosystems and feedbacks to climate. NGEE Arctic was supported by the Department of Energy's Office of Biological and Environmental Research.The NGEE Arctic project had two field research sites: 1) located within the Arctic polygonal tundra coastal region on the Barrow Environmental Observatory (BEO) and the North Slope near Utqiagvik (Barrow), Alaska and 2) multiple areas on the discontinuous permafrost region of the Seward Peninsula north of Nome, Alaska.Through observations, experiments, and synthesis with existing datasets, NGEE Arctic provided an enhanced knowledge base for multi-scale modeling and contributed to improved process representation at global pan-Arctic scales within the Department of Energy's Earth system Model (the Energy Exascale Earth System Model, or E3SM), and specifically within the E3SM Land Model component (ELM).

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Emergent interface vibrational structure of oxide superlattices

As the length scales of materials decrease, the heterogeneities associated with interfaces become almost as important as the surrounding materials. This has led to extensive studies of emergent electronic and magnetic interface properties in superlattices. However, the interfacial vibrations that affect the phonon-mediated properties, such as thermal conductivity, are measured using macroscopic techniques that lack spatial resolution. Although it is accepted that intrinsic phonons change near boundaries, the physical mechanisms and length scales through which interfacial effects influence materials remain unclear. Here we demonstrate the localized vibrational response of interfaces in strontium titanate–calcium titanate superlattices by combining advanced scanning transmission electron microscopy imaging and spectroscopy, density functional theory calculations and ultrafast optical spectroscopy. Structurally diffuse interfaces that bridge the bounding materials are observed and this local structure creates phonon modes that determine the global response of the superlattice once the spacing of the interfaces approaches the phonon spatial extent. Our results provide direct visualization of the progression of the local atomic structure and interface vibrations as they come to determine the vibrational response of an entire superlattice. Direct observation of such local atomic and vibrational phenomena demonstrates that their spatial extent needs to be quantified to understand macroscopic behaviour. Tailoring interfaces, and knowing their local vibrational response, provides a means of pursuing designer solids with emergent infrared and thermal responses.

36 MATERIALS SCIENCE↗
DOE OSTI2022

The short ionizing photon mean free path at z = 6 in Cosmic Dawn III, a new fully coupled radiation-hydrodynamical simulation of the Epoch of Reionization

ABSTRACT Recent determinations of the mean free path of ionizing photons (mfp) in the intergalactic medium (IGM) at z = 6 are lower than many theoretical predictions. In order to gain insight, we investigate the evolution of the mfp in our new massive fully coupled radiation-hydrodynamics cosmological simulation of reionization: Cosmic Dawn III (CoDa III). CoDa III’s scale ($\rm 94^3 \, cMpc^3$) and resolution ($\rm 8192^3$ grid) make it particularly suitable to study the IGM during reionization. The simulation was performed with ramses-cudaton on Summit, and used 13 1072 processors coupled to 24 576 GPUs, making it the largest reionization simulation, and largest ever ramses simulation. A superior agreement with global constraints on reionization is obtained in CoDa III over Cosmic Dawn II (CoDa II), especially for the evolution of the neutral hydrogen fraction and the cosmic photoionization rate, thanks to an improved calibration, later end of reionization (z = 5.6), and higher spatial resolution. Analysing the mfp, we find that CoDa III reproduces the most recent observations very well, from z = 6 to z = 4.6. We show that the distribution of the mfp in CoDa III is bimodal, with short (neutral) and long (ionized) mfp modes, due to the patchiness of reionization and the coexistence of neutral versus ionized regions during reionization. The neutral mode peaks at sub-kpc to kpc scales of mfp, while the ionized mode peak evolves from 0.1 Mpc h−1 at z = 7 to ∼10 Mpc h−1 at z = 5.2. Computing the mfp as the average of the ionized mode provides the best match to the recent observational determinations. The distribution reduces to a single neutral (ionized) mode at z > 13 (z < 5).

Lewis, Joseph S. W. (ORCID:0000000179178474)↗
DOE OSTI2022

Continental United States may lose 1.8 petagrams of soil organic carbon under climate change by 2100

Abstract Aims High‐resolution information on soils’ vulnerability to climate‐induced soil organic carbon (SOC) loss can enable environmental scientists, land managers, and policy makers to develop targeted mitigation strategies. This study aims to estimate baseline and decadal changes in continental US surface SOC stocks under future emission scenarios. Location Continental United States. Time period 2014–2100. Methods We used recent SOC field observations ( n = 6,213 sites), environmental factors ( n = 32), and an ensemble machine learning (ML) approach to estimate baseline SOC stocks in surface soils across the continental United States at 100‐m spatial resolution, and decadal changes under the projected climate scenarios of Coupled Model Intercomparison Project Phase Six (CMIP6) earth system models (ESMs). Results Baseline SOC projections from ML approaches captured more than 50% of variability in SOC observations, whereas ESMs represented only 6–16% of observed SOC variability. ML estimates showed a mean total loss of 1.8 Pg C from US surface soils under the high‐emission scenario by 2100, whereas ESMs showed no significant change in SOC stocks with wide variation among ESMs. Both ML and ESM predictions agree on the direction of SOC change (net emissions or sequestration) across 46–51% of continental US land area. These differences are attributable to the high‐resolution site‐specific data used in the ML models compared to the relatively coarse grid represented in CMIP6 ESMs. Main conclusions Our high‐resolution estimates of baseline SOC stocks, identification of key environmental controllers, and projection of SOC changes from US land cover types under future climate scenarios suggest the need for high‐resolution simulations of SOC in ESMs to represent the heterogeneity of SOC. We found that the SOC change is sensitive to key soil related factors (e.g. soil drainage and soil order) that have not been historically considered as input parameters in ESMs, because currently more than 95% variability in the SOC of CMIP6 ESMs is controlled by net primary productivity, temperature, and precipitation. Using additional environmental factors to estimate the baseline SOC stocks and predict the future trajectory of SOC change can provide more accurate results.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Implementation of an impurity diagnostic suite on the Pegasus-III experiment

A suite of diagnostics used to assess impurity content and dynamics has been updated, upgraded, and installed on the Pegasus-III Experiment. Typical plasma parameters during local helicity injection start-up are τ shot ~10 ms, n e ~ 1 × 10 19 m –3 , and T e ~ 50 eV. The deployed diagnostics are compatible with this modest temperature and density regime and provide species identification, source localization, and estimation of radiation losses. Impurity species are determined by recording time-evolving, single line-of-sight spectra at 1.25 kfps using a SPRED (Survey, Poor Resolution, Extended Domain) vacuum ultraviolet spectrometer. SPRED is equipped with 450 g/mm grating, giving a spectral resolution of 0.33 nm and a spectral range from ~10 to 110 nm, useful to identify light impurity species in this temperature and density range. An absolutely calibrated spectrometer that collects light from the plasma at R tan = 15.9 cm and Δt ≥ 2 ms is used as a visible survey spectrometer and for continuum measurements. The radiated power from the plasma is estimated with a photodiode-based diagnostic. Two 16-channel absolute extreme ultraviolet diode arrays are placed behind pinhole apertures, resulting in 32 lines of sight at Z = 0, with a spatial resolution of 2–3 cm and a time response of 60 kHz. A photometrically calibrated collinear D α /near infrared filtered photodiode-based system measures the D α emission and around 1040 nm. All these instruments have been designed to suppress electromagnetic interference from megawatt-class switching power supplies.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Insights from application of a hierarchical spatio-temporal model to an intensive urban black carbon monitoring dataset

Existing regulatory pollutant monitoring networks rely on a small number of centrally located measurement sites that are purposefully sited away from major emission sources. While informative of general air quality trends regionally, these networks often do not fully capture the local variability of air pollution exposure within a community. Recent technological advancements have reduced the cost of sensors, allowing air quality monitoring campaigns with high spatial resolution. The 100×100 black carbon (BC) monitoring network deployed 100 low-cost BC sensors across the 15 km 2 West Oakland, CA community for 100 days in the summer of 2017, producing a nearly continuous site-specific time series of BC concentrations which we aggregated to one-hour averages. Leveraging this dataset, we employed a hierarchical spatio-temporal model to accurately predict local spatio-temporal concentration patterns throughout West Oakland, at locations without monitors (average cross-validated hourly temporal R 2 =0.60). Using our model, we identified spatially varying temporal pollution patterns associated with small-scale geographic features and proximity to local sources. In a sub-sampling analysis, here we demonstrated that fine scale predictions of nearly comparable accuracy can be obtained with our modeling approach by using ~30% of the 100×100 BC network supplemented by a shorter-term high-density campaign.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Estimating snow cover from high-resolution satellite imagery by thresholding blue wavelengths

We report the extent and duration of snow cover, a critical component of the hydrologic cycle and the global climate system, is expected to shift dramatically under climate change. Therefore, developing high-resolution assessments of snow cover change is crucial for estimating the impact of changing snow cover on watershed and ecosystems processes in cold regions. Remote sensing tools provide a powerful method for mapping snow-covered area (SCA) across a landscape. The most common method for estimating SCA utilizes the normalized difference snow index (NDSI), which relies on spectral measurements in the shortwave-infrared wavelengths (SWIR). NDSI can effectively estimate catchment- to regional-scale SCA, but it cannot be used to assess fine-scale SCA because of current limitations on the spatial resolution of satellite-derived SWIR measurements. Here, we map SCA using a threshold of blue wavelengths and high-resolution satellite imagery. The thresholding method, which we call the Blue Snow Threshold algorithm (BST), has previously been used with digital camera imagery. We refine and automate the algorithm for use with cloud-free high-resolution satellite imagery and find that the BST can be used to assess fine-scale SCA. For validation, we compared BST-derived estimates of SCA to a) airborne lidar surveys, b) Landsat fractional SCA, and c) snow disappearance dates from Snow Telemetry (SNOTEL) stations. When compared to airborne lidar surveys of SCA, the BST predicted SCA had a range of F-scores between 0.81 and 0.94 in four study areas in California and Colorado. We also found general agreement between SCA and snow disappearance at multiple SNOTEL sites across the western United States. Given the relatively recent availability of high-resolution satellite imagery with spectral measurements in the visible wavelengths but lacking in SWIR, the BST offers a reliable and easy-to-apply tool for examining fine-scale snow-related processes.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Development of a River Dynamical Core for E3SM to simulate compound flooding on Exascale-class heterogeneous supercomputers

Flooding events pose significant risk to human life, property, and infrastructure. Physically-consistent quantification of altered flood risks in global models requires hyper-resolution (~1 km) or fine flood simulations using two-dimensional (2D) physics schemes, both of which are unavailable in the current generation Earth System Models. Here, in this work, we have developed the River Dynamical Core (RDycore), which is an open-source, 2D shallow water equation (SWE) library for the U.S. Department of Energy's Energy Exascale Earth System Model (E3SM). RDycore uses PETSc and libCEED libraries that allows it to run efficiently on CPUs and GPUs, as well as select a time-integration algorithm at runtime without requiring any code modifications. RDycore achieves spatial error convergence rates for problems with analytical and manufactured solutions similar to those reported previously in the literature, or consistent with the implemented first-order spatial discretization scheme. RDycore's accuracy in predicting flooding for a well-studied dam break problem is comparable to existing SWE models. For a problem with 471 million grid cells, RDycore achieves a speedup of 6.6x and 7.6x on GPUs compared to CPUs when using 320 compute nodes on DOE's Perlmutter and Frontier supercomputers, respectively. The one-way coupling of the RDycore library within E3SM is demonstrated by performing multiple 5-day flooding simulations during Hurricane Harvey driven by five precipitation datasets. The E3SM--RDycore simulations at 30 m spatial resolution accurately simulate maximum water height during the hurricane when benchmarked against a previously published study and achieve a speedup of 15x (Perlmutter) and 21x (Frontier) on GPUs relative to CPUs. The work presented here is the foundational step in providing hardware and algorithmic portability framework for simulating kilometer-scale river dynamics within E3SM.

Flood Simulation↗
DOE OSTI2023

Imaging the electron charge density in monolayer MoS 2 at the Ångstrom scale

Four-dimensional scanning transmission electron microscopy (4D-STEM) has recently gained widespread attention for its ability to image atomic electric fields with sub-Ångstrom spatial resolution. These electric field maps represent the integrated effect of the nucleus, core electrons and valence electrons, and separating their contributions is non-trivial. In this paper, we utilized simultaneously acquired 4D-STEM center of mass (CoM) images and annular dark field (ADF) images to determine the projected electron charge density in monolayer MoS 2 . We evaluate the contributions of both the core electrons and the valence electrons to the derived electron charge density; however, due to blurring by the probe shape, the valence electron contribution forms a nearly featureless background while most of the spatial modulation comes from the core electrons. Our findings highlight the importance of probe shape in interpreting charge densities derived from 4D-STEM and the need for smaller electron probes.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Bulk and Short-Circuit Anion Diffusion in Epitaxial Fe 2 O 3 Films Quantified Using Buried Isotopic Tracer Layers

Self-diffusion is a fundamental physical process that, in solid materials, is intimately correlated with both microstructure and functional properties. In this work, a universal approach is presented to precisely characterize self-diffusion in ionic solids by isotopically enriching anions and/or cations at specific locations within an epitaxial film stack, and characterize their redistribution at high spatial resolution with atom probe tomography. Nanoscale anion diffusivity is quantified in epitaxial α-Fe 2 O 3 thin films deposited by molecular beam epitaxy with a thin (10 nm) buried layer highly enriched in 18 O. The isotopic sensitivity of the atom probe allows precise measurement of 18 O distribution across the sharp interfaces between this layer and the surrounding Fe 2 O 3 after annealing. Short-circuit anion diffusion through 1D and 2D structural defects in Fe 2 O 3 are also directly visualized in 3D. This versatile approach to study precisely tailored thin film samples at high spatial and mass fidelity will facilitate a deeper understanding of atomic-scale diffusion phenomena.

36 MATERIALS SCIENCE↗
DOE OSTI2020

X-ray and gamma-ray tomographic imaging of fuel relocation inside sodium fast reactor test assemblies during severe accidents

The present work reports on x-ray and gamma-ray high-spatial resolution computerized tomography measurements of the Pin Bundle Metallic Fuel Relocation (PBR) assemblies tested in the Metallic Uranium Safety Experiment (MUSE) facility at Argonne National Lab (ANL). The aim of the study was to characterize fuel relocation structures that develop during severe core accidents pertaining to SFR assemblies; these include but are not limited to advanced core disruption recreated in the PBR-1 assembly, and cladding breach recreated in the PBR-2 assembly. We report the x-ray tomography measurements were able to resolve small quantities of relocation fuel; with increased presence of relocation fuel, the x-ray 1measurements spatially mapped the material but could not resolve the inner regions of these. The gamma-tomography measurements showed improved results, resolving the relocation structures in great detail. The upper plenum of the PBR-1 assembly where the molten uranium was initially inserted presented high structural damage, reflected by the partial and complete disintegration of the central rods. Relocation fuel filled the subchannels, adhering to surviving cladding walls and the assembly casing. In the lower portion of the measured section, the tomogram degrades due to photon starvation effects hinting at the increased amount of relocation fuel potentially plugging the assembly; flow blockage in this section was difficult to determine due to the tomogram’s degradation from photon starvation. Small fragments were observed further down the assembly, dislodged from the initial insertion of the molten material. This section was used as an unperturbed assembly reference, with a calculated blockage of less than 1% from the present fragments. The PBR-2 assembly was characterized by columnar relocation structures propagating through the subchannels. Three of the relocation structures were captured in the measured section, with evidence of cross migration on to adjacent subchannels. The measured section in this assembly captures the leading edge of two structures. The calculated flow blockage was 16% in the planes where the three relocation structures are present, but this quickly decreases to approximately 5% past the leading edge of two of the structures.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2022

LAI, EVI, NDVI, and kNDVI in 23 pantropical forests affected by 21 cyclones

Statement of purpose: Cyclones alter the function and composition of tropical forests, making effects of intensifying cyclones on carbon-rich forests a critical topic of study. Here, we quantified cyclone-induced damage and recovery of 21 cyclone disturbances affecting 23 pantropical forest sites between 1988-2017 utilizing leaf area index (LAI), enhanced vegetation index (EVI), normalized difference vegetation index (NDVI), and transformed NDVI (kNDVI) values from Google Earth Engine. Field observations collected in a meta-analysis (Bomfim et al., 2022, in review) were used to ground-truth and test effects of soil resource availability and disturbance factors on damage and recovery. This meta-analysis also served as the basis to begin vegetation index extraction, utilizing unique site and date combinations, from tropical forests effect by cyclone disturbances. We began collecting NDVI (5km resolution) from the NOAA Climate Data Record (CDR) of AVHRR Normalized Difference Vegetation Index (NDVI), Version 5 data product (Vermote, 2019) for all case studies included, 42. Next, we began extracting Landsat data from Landsat 4, 5, and 8, courtesy of the U.S. Geological Survey, in search of higher resolution data. We selected a 3 by 3 Landsat pixel area, leading to a 90m resolution data extraction. The specific imagery used includes Landsat 4 USGS Landsat 4 TM Collection 1 Tier 1 TOA (top of atmosphere) Reflectance, Landsat 5 USGS Landsat 5 TM (thematic mapper) Collection 1 Tier 1 TOA Reflectance, and Landsat 8 USGS Landsat 8 Collection 1 Tier 1 TOA Reflectance. Within Google Earth Engine, we selected the date and location (latitude and longitude), calculated NDVI, kNDVI, and EVI utilizing Landsat bands (see metadata_NGEE-tropics_cyclones), and extracted post- and pre-cyclone values for each case study to calculate cyclone-induced change in the vegetative index. Due to limited spatial resolution of Landsat remote sensing data, MODIS products were investigated next. First, the MOD13Q1.006 Terra Vegetation Indices 16-Day Global 250m product was used to extract 250m EVI and NDVI (Didan, 2015) and then the MCD15A3H.006 MODIS Leaf Area Index/FPAR 4-Day Global 500m product product was used to extract LAI 500m (Myneni et al., 2015). Pre- and post-cyclone values, change in the vegetative index, and standard deviation for all values are included in the main csv (see case_study_data.csv) for all vegetative indices collected, including LAI 500m, EVI 250m, NDVI 250m, NDVI 90m, kNDVI 90m, EVI 90m, and NDVI 5km. Lastly, recovery values were calculated utilizing a standardization method (see metadata_NGEE-tropics_cyclones) and recovery values for MODIS (see MODIS_recovery.csv) and Landsat (Landsat_recovery.csv) data are included.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Improving Subsurface Stress Characterization for Carbon Dioxide Storage Projects by Incorporating Machine Learning Techniques

The overall objective of this project is to develop a framework for reliable characterization and prediction of the state of stress in the overburden and underburden (including the basement) in CO 2 storage reservoirs using machine learning and integrated geomechanics and geophysical methods. Specifically, we propose to develop workflow encompassing of technologies and/or methods to predict stress and pressure changes due to CO 2 injection in an active tertiary recovery site and their impacts on subtle fault activation, fractures and occurrence of microseismic events and compare responses to field observations. In this project, we anticipate using dataset from the Farnsworth field Unit (FWU) which is operated by Purdure Petroleum. A novel elastic-waveform VSP inversion technique will be used to estimate high-resolution spatial and temporal changes of elastic moduli in CO 2 storage reservoirs, which will be combined with velocity-stress relationship derived from laboratory tests to obtain subsurface pressure and stress. Clustered microseismic data will be jointly inverted for improved focal mechanisms. Least-squares reverse-time migration of microseismic waveform data will be performed to directly image fracture/fault zones. Additionally, a deep neural network machine learning technique with convolutional and recurrent layers will be used for learning the spectro-temporal structures in microseismic waveforms. The results of this geotechnical data analysis will be integrated to develop a high-resolution 3D mechanical earth model extending from the overburden sealing formations to the underburden including the basement. Mechanical properties will be derived through integration of mechanical logs, tests, available results from chemo-mechanical laboratory tests, and elastic inversion of seismic data using a combination of Bayesian and stochastic methods as well as machine learning technique. Failure features (faults/fractures) will be represented and/or modeled based on seismic and core data analysis. A transient hydrodynamic-geomechanical model will be developed through coupling with the calibrated FWU reservoir simulation model. The full physics coupled model will be used to train a reduced order proxy model using machine learning algorithm for estimating stress which will then be used with appropriate constitutive relationships and forward seismological models to simulate pressure changes and induced microseismicity. An advanced optimization framework will be developed to perform a history match to minimize error between field observations and simulated. The history matched proxy model will be verified against the full-physics equivalent. The field observations that will be used in the coupled model calibration process include pressure/stress inverted from VSP, moment magnitude from microseismic analysis, real time downhole pressure measurements, production and injection data. Parameter sensitivity and uncertainty analysis will be performed to characterize the impact of model parameter uncertainty on stress estimates. The proposed project will have significant impact on future field implementation of the proposed technology. Because the project field site is an ongoing CO 2 EOR development, the value of the new technology will be demonstrated in an operational context and evaluated as a viable risk mitigation strategy. Cost/benefit will be evaluated together with the various commercial incentives for CO 2 sequestration available to oil and gas operators. The extensive available dataset and ongoing data acquisition under the SWP Phase III work plan provides flexibility for investigation of multiple approaches and reduces technical risk.

58 GEOSCIENCES↗
DOE OSTI2020

Characterization of Extreme Wave Conditions for Wave Energy Converter Design and Project Risk Assessment

Best practices and international standards for determining n-year return period extreme wave (sea states) conditions allow wave energy converter designers and project developers the option to apply simple univariate or more complex bivariate extreme value analysis methods. The present study compares extreme sea state estimates derived from univariate and bivariate methods and investigates the performance of spectral wave models for predicting extreme sea states at buoy locations within several regional wave climates along the US East and West Coasts. Two common third-generation spectral wave models are evaluated, a WAVEWATCH III® model with a grid resolution of 4 arc-minutes (6–7 km), and a Simulating WAves Nearshore model, with a coastal resolution of 200–300 m. Both models are used to generate multi-year hindcasts, from which extreme sea state statistics used for wave conditions characterization can be derived and compared to those based on in-situ observations at National Data Buoy Center stations. Comparison of results using different univariate and bivariate methods from the same data source indicates reasonable agreement on average. Discrepancies are predominantly random. Large discrepancies are common and increase with return period. There is a systematic underbias for extreme significant wave heights derived from model hindcasts compared to those derived from buoy measurements. This underbias is dependent on model spatial resolution. However, simple linear corrections can effectively compensate for this bias. A similar approach is not possible for correcting model-derived environmental contours, but other methods, e.g., machine learning, should be explored.

16 TIDAL AND WAVE POWER↗
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↗