Search NASA⌕ Search

SEARCH · Search NASA

Results for “heterogeneous spatial distributions”

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 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Quantifying the chemical, electrochemical heterogeneity and spatial distribution of (poly) sulfide species using Operando SANS

In-situ and operando measurements are done to gain a better understanding of the precipitation mechanisms during charge and discharge in lithium sulfur batteries. A carbon felt networked with microfibers consisting of pores, 2 nm and smaller, is used as a freestanding sulfur host. Three different methods of sulfur infiltration are explored in order to determine the best one to fill most of the pores. It was identified to be a melt infiltration method in a vacuum oven, where the pores including the ultra-micropores, are successfully filled. In-situ electrochemical impedance spectroscopy measurements show a solid product formation occurring at the sulfur cathode, both during the high voltage plateau and at the end of discharge. In a 3-electrode EIS measurement, a similar solid product formation on the Li counter electrode due to its reaction with polysulfides is also observed. Operando small angle neutron scattering measurements show the solid product formation, in the carbon, both near the beginning and at the end of discharge, confirming the precipitation data via contrast changes as a function of charge and discharge. It is shown that Li 2 S precipitates in the pores at the beginning and the end of discharge, whereas S8 precipitates on the surface of the carbon felt. Finally, this lithium sulfur system shows both the quasi-solid-state and the solid-liquid-solid reactions in a typical ethereal electrolyte solution with a two-plateau discharge profile.

25 ENERGY STORAGE↗
DOE OSTI2023

Derivation of heterogeneous material distributions and their sensitivity to HM-coupled two-phase flow models exemplified with the LASGIT experiment

Abstract Advective gas transport in bentonite, a possible buffer material in repositories for radioactive materials, is difficult to simulate in numerical continuum models, partly due to the complicated microstructure of bentonite. To generate reliable models of repositories nevertheless, spatially distributed heterogeneous material properties can be used to allow localization of gas flow. In this study, a pore-size-dependent stochastic approach of the gas entry pressure is derived from Mercury Intrusion Porosimetry, which is used to replicate measurements from the LASGIT experiment. In addition, three benchmark tests are simulated to investigate the dependence of heterogeneous distributions of material properties on the mesh discretization, the temporal dependence, and the coupling between the processes influenced by the heterogeneous parameters. The numerical modeling results of the LASGIT experiment show that the onset of gas flow into the system and the subsequent increase in pressure and stress can be well reproduced using heterogeneous distributions. Compared to a model with homogeneous material properties, heterogeneous distributions may allow the generation of dilatancy-controlled microfractures—an important feature with regard to the advective gas flow in bentonites. However, it can be observed that the heterogeneous distributions in LASGIT are less significant, as technical gaps or differences in material types could have a greater impact.

Environmental Sciences & Ecology↗
DOE OSTI2022

Detection and attribution of long-term and fine-scale changes in spring phenology over urban areas: A case study in New York State

Spring phenology plays an essential role in climate change, terrestrial ecosystem, and public health. Field-based monitoring and understanding of changes in spring phenology for long periods and in large regions are challenging due to the limited in-site observations. Space-based remotely sensed observations offer great potentials for monitoring decadal spring phenology changes from regional to global scales. However, the coarse-scale remotely sensed observations are insufficient to capture fine-scale spring phenology dynamics, especially in urban areas, and this makes it challenging for understanding the combined effects of climate change and urbanization on spring phenology. We derived the start of phenology season (SOS) in New York State using 30 m Landsat observations from 1990 to 2015 to understand the impact of the environment and urbanization on SOS. The results show that SOS for different years reveals heterogeneous spatial distribution. Most regions of New York State have been experiencing significant spring phenology changes in form of earlier onset of vegetation greening, ranging from 0.2 to 0.6 day/year during 1990 to 2015, and this trend varies slightly with latitudes and urbanization levels. Further, spatial correlation analysis shows that the increase in temperature and urbanization could both promote the advancement of SOS. However, the effect of urbanization (partial correlation coefficient (R) ranges from −0.289 to −0.542) on SOS is greater than the effect of temperature (R ranges from 0.006 to −0.192). The study generates a high spatio-temporal resolution spring phenology dataset for ecological, environmental and public health studies, especially in urban areas, and reveals the importance of better accounting for the urbanization effects when quantifying the SOS dynamics in phenology models.

Landsat↗
DOE OSTI2020

Oxophilicity Drives Oxygen Transfer at a Palladium–Silver Interface for Increased CO Oxidation Activity

A single-layer AgO x phase grown on Ag(111) efficiently transfers oxygen to Pd domains at room temperature, rendering the Pd-decorated surface highly reactive toward CO oxidation. Here, oxygen transfer from AgO x to Pd and the surface reactivity toward CO were investigated as a function of the Pd coverage using X-ray photoelectron spectroscopy, surface infrared spectroscopy of adsorbed CO, temperature-programmed reaction spectroscopy, and density functional theory (DFT) calculations. Our results show that all of the oxygen from the AgO x layer (~0.375 monolayer) migrates to the surface of Pd during formation of a nearly complete Pd bilayer at 300 K and that the oxygen coverages generated on Pd increase as the Pd cluster size decreases, reaching values that exceed the oxygen concentration in the AgO x layer by as much as a factor of 2. Experimental measurements and DFT calculations show that preferential binding of oxygen on the edges of the Pd clusters enhances the oxygen coverage on Pd clusters of decreasing size and produces a heterogeneous spatial distribution of oxygen. CO adsorbs in high coverages at 100 K by binding on both the terraces and O-rich edges of the Pd clusters. During subsequent heating, oxidation of the adsorbed CO consumes nearly all of the oxygen that transferred from AgO x to the Pd domains; in contrast, the pure AgO x layer exhibits limited reactivity toward CO adsorbed at 100 K. These results demonstrate that differences in oxophilicity drive facile oxygen transfer from Ag to the edges of Pd nanoclusters and thereby give rise to an efficient pathway for CO oxidation on bimetallic PdAg surfaces. The cooperation between the Pd and Ag domains results in near-interfacial chemistry that may be broadly important in catalysis by bimetallic alloys.

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

Delayed Neutron Temporal Signatures for Uranium Enrichment Measurement NA-241 SGTech (Final Report)

Nondestructive determination of uranium enrichment is a core capability for nuclear material accounting and control (NMAC) and safeguards verification measurements; however, traditional gamma spectroscopy-based techniques for enrichment measurement rely on significant assumptions of material composition and geometry, precluding their use in scenarios where a heterogeneous spatial distribution of enrichments is encountered. As an alternative, we are developing a technique to use delayed neutron temporal signatures for the measurement of uranium enrichment. Each uranium isotope has unique delayed neutron group yields, resulting in a unique delayed neutron decay time profile which can be analyzed to determine enrichment without the need for calibration sources. As part of this effort, we performed a series of measurement campaigns in which we used an active well coincidence counter (AWCC) retrofitted with commercial D-D and D-T generators to evaluate the operational characteristics of this method in response to a set of uranium enrichment and mass standards, as well as representative diversion scenarios in which either “concealed” enriched uranium is shielded by depleted uranium or declared enriched uranium is “hollowed out” and replaced with a central region of depleted uranium. A standard operating procedure and best practices were compiled to facilitate the use of delayed neutron-based enrichment measurements for international safeguards inspections.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2026

Permafrost Thaw, Uneven Subsidence and Projected Drying of Ice-wedge Polygon Tundra: Modeling Archive

This dataset is a model archive of the paper Permafrost Thaw, Uneven Subsidence and Projected Drying of Ice-wedge Polygon Tundra (in prep) to support a modeling study investigating how projected increases in Arctic temperature and precipitation will jointly influence hydrologic conditions in ice-rich tundra landscapes. With this dataset, this study is to address the research question: Will Arctic tundra landscapes become wetter or drier with increasing precipitation and temperature in the future when thaw-induced ground subsidence and associated microtopographic evolution are represented? The simulations focus on ice-wedge polygon tundra, a widespread form of ice-rich permafrost terrain that is highly sensitive to thaw-driven landscape change. This dataset contains model input and output data for four study watersheds in Alaska: Anaktuvuk, Utqiagvik (formerly Barrow), Brooks Foothills, and Prudhoe Bay. Simulations were performed using the Advanced Terrestrial Simulator (ATS, v1.5), a physics-rich integrated surface–subsurface hydrologic model. For each watershed, ten modeling cases were performed representing two landscape evolution conditions (with subsidence and without subsidence) combined with five climate forcing scenarios derived from Shared Socioeconomic Pathways (SSP5, SSP5 with precipitation trend, SSP2, SSP2 with precipitation trend, and SSP2 with double precipitation trend). Particularly, for each watershed under the forcing SSP2 with precipitation trend, there are two additional simulations considering spatially heterogeneous subsidence distributions: one assumes randomly distributed scaling and the other includes elevation dependent distribution scaling. These simulations span 1980–2099 and include spin-up runs (1980–2009) followed by transient projections (2010–2099). To facilitate reproducibility of simulations, all datasets are organized by watershed. For each study watershed, the dataset contains: (1) Pre-partitioned mesh files for 32-core modeling (.par.32.XX), located in EACH_WATERSHED/mesh/basin; and also a non-partitioned mesh file (.exo) located in EACH_WATERSHED/mesh; (2) Climate forcings corresponding to the five SSP scenarios (.h5), located in EACH_WATERSHED/data; (3) Final states (.h5) from column spin-up modeling used to initialize historical watershed-scale spin-up runs from 1980 to 2009, located in EACH_WATERSHED/PreSpinupHistorical; (4) Final states (.h5) of historical watershed-scale spin-up runs from 1980 to 2009 used to initialize projection runs, located in EACH_WATERSHED/Spinup_daymetERA5; (5) ATS modeling input files (.xml), located in EACH_WATERSHED/EACH_SIMULATION_SCENARIO/inputfiles; (6) ATS modeling output files (.dat), located in in EACH_WATERSHED/EACH_SIMULATION_SCENARIO/combined_obs; (7) For the Brooks Foothills watershed, additional spatial model outputs are provided (.h5) for selected years (2033 and 2093) used to generate spatial figures in this study, located in Brooksfoothills/EACH_SIMULATION_SCENARIO/results-WITH/WITHOUT_SUBSIDENCE-year2033/2093. All data files with suffix .h5 can be accessible through Python h5py, and all data files with suffix of .dat can be imported by Python pandas. Mesh file with .exo can be visualized through Paraview or read by Python netCDF. The Next-Generation Ecosystem Experiments in the Arctic (NGEE Arctic) project is a research effort to reduce uncertainty in the Department of Energy’s Energy Exascale Earth System Model (E3SM) by developing a predictive understanding of Arctic tundra ecosystems underlain by permafrost and to quantify feedbacks from the Arctic tundra to the Earth system. NGEE Arctic is supported by the Department of Energy's Office of Biological and Environmental Research. Over Phases 1–3, observations made by the NGEE Arctic team across a gradient of permafrost landscapes in Arctic Alaska improved the representation of tundra processes in the land surface component of E3SM (the E3SM Land Model, ELM). Model improvements emphasized unique aspects of permafrost environments and explored reductions in model complexity while retaining predictive power. The Arctic-informed ELM developed by NGEE Arctic has been used to make novel predictions on processes ranging from permafrost thaw to soil biogeochemical cycling to Earth system feedbacks associated with the unique characteristics of tundra plants. In Phase 4, the NGEE Arctic team is evaluating our new predictive understanding under novel conditions across the Arctic domain. In collaboration with partners at long-term pan-Arctic research sites we are examining whether an Arctic-informed ELM can faithfully simulate interactions among surface and subsurface processes at site, regional, and pan-Arctic scales. In turn, we are using variety of tools to dynamically extend and evaluate ELM inference, with an emphasis on data synthesis and pan-Arctic model evaluation, reintegration of code with an evolving E3SM, scaling across heterogeneous Arctic landscapes, and the appropriate representation of the impacts of increasingly frequent Arctic disturbances.

EARTH SCIENCE > ATMOSPHERE > PRECIPITATION↗
DOE OSTI2021

Competition within low-density bacterial populations as an unexpected factor regulating carbon decomposition in bulk soil

Bacterial decomposition of organic matter in soils is generally believed to be mainly controlled by the access bacteria have to organic substrate. The influence of bacterial traits on this control has, however, received little attention. Using the substrate-dependent Monod growth model, we develop a bioreactive transport model to screen the interactive impacts of spatial dispersion and bacterial traits on mineralization. Bacterial traits primarily involved in the bacterial response to the substrate concentration, such as the maximum specific uptake rate and efficiency, the adaptation time of the uptake rate and the initial population density, are considered. We compare the model results with two sets of previously performed cm-scale soil-core experiments in which the mineralization of the pesticide 2,4-D was measured under well-controlled initial distributions and transport conditions. Bacterial dispersion away from the initial substrate location induced a significant increase in 2,4-D mineralization. Additionally, it reveals an increase of specific uptake rates at lower bacterial densities, more than compensating the decrease of specific uptake rates caused by substrate dilution. This regulation of bacterial activities by density, caused by the local depletion of substrate by competing bacteria, becomes dominant for bacteria with an efficient uptake of substrate at low substrate concentrations (a common feature of oligotrophs). Such oligotrophs, commonly found in soils, compete with each other for substrate even at remarkably low population densities. The ratio-dependent Contois growth model, which includes a density regulation in the expression of the uptake efficiency, is more accurate and convenient to calibrate than the substrate-dependent Monod model, at least under these conditions. In view of their strong interactions, bioreactive and transport processes cannot be handled independently but should be integrated, in particular when reactive processes of interest are carried out by oligotrophs.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

Cloud Atlas: Unraveling the Vertical Cloud Structure with the Time-series Spectrophotometry of an Unusually Red Brown Dwarf

Rotational modulations of emission spectra in brown dwarf and exoplanet atmospheres show that clouds are often distributed non-uniformly in these ultracool atmospheres. The spatial heterogeneity in cloud distribution demonstrates the impact of atmospheric dynamics on cloud formation and evolution. In this study, we update the Hubble Space Telescope (HST) time-series data analysis of the previously reported rotational modulations of WISEP J004701+680352—an unusually red late-L brown dwarf with a spectrum similar to that of the directly imaged planet HR 8799e. We construct a self-consistent spatially heterogeneous cloud model to explain the HST and the Spitzer time-series observations, as well as the time-averaged spectra of WISE 0047. Additionally, in the heterogeneous cloud model, a cloud thickness variation of around one pressure scale height explains the wavelength dependence in the HST near-IR spectral variability. By including disequilibrium CO/CH 4 chemistry, our models also reproduce the redder $J-{K}_{{\rm{s}}}$ color of WISE 0047 compared to that of field brown dwarfs. We discuss the impact of vertical cloud structure on atmospheric profile and estimate the minimum eddy diffusivity coefficient for other objects with redder colors. Our data analysis and forward modeling results demonstrate that time-series spectrophotometry with a broad wavelength coverage is a powerful tool for constraining heterogeneous atmospheric structure.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Characterization of Temperature Heterogeneity in Utility-Scale Power Plant Boilers by Spatially Distributed Ultrasonic Measurements

In extreme environments, even hardened insertion sensors fail quickly. For such environments, we have developed an ultrasonic (US) method for measuring the spatial distribution of temperatures in solid materials and, specifically, across containments of extreme processes. By deploying US sensors in multiple locations, spatial heterogeneity in temperatures inside harsh environments may be characterized. In this note, we update our progress on testing the developed approach in an industrial setting and its performance in describing the temperature distribution inside a 500 MW (electrical) coalfired utility boiler of an electrical power generation plant. We use waveguides (WG) structured to contain echogenic features that produce a train of echoes in response to an excitation pulse. The time of flight between echoes encodes the information on the temperature distribution in the corresponding segment of the WG, which we reconstruct using its parametrization. Five waveguides were welded to the boiler’s heat exchange surface (waterwall) and different locations at the same boiler elevation and produced reliable US waveforms 18 months after the installation. For several weeks during the latest trial, we performed US measurements in five locations spanning half of the boiler’s 14-meter width. The measurements were acquired during regular operation, including load cycling to adjust for demand and intermittent contributions of renewable power sources, and during boiler shutdown for emergency maintenance and the subsequent restart. The interpretation of US waveforms by signal processing resulted in an accurate estimation of temperature distribution along the waveguides. We captured daily cyclical variations in demand-following load. Measurements performed in multiple locations revealed unexpected temperature variations across the boiler. We conclude that an array of developed sensors can provide responsive and spatial temperature measurements while maintaining functionality despite prolonged exposure to an extreme environment.

Walton, Kenneth↗
DOE OSTI2023

DUET - Distribution of Understory using Elliptical Transport: A mechanistic model of leaf litter and herbaceous spatial distribution based on tree canopy structure

Heterogeneity in surface fuels produced by overstory trees and understory vegetation is a major driver of fire behavior and ecosystem dynamics. Previous attempts at predicting tree leaf and needle litter accumulation over time have been constrained in scope to probabilistic models that consider a limited number of key factors influencing tree litter dispersal patterns and decomposition processes. Here we present a mechanistic model for estimating variation in surface fuels called the Distribution of Understory using Elliptical Transport (DUET). DUET uses a pre-generated voxelated canopy array and distributes the leaf litter from the trees in an elliptical shape based on tree species characteristics, wind data, and location-specific features. DUET then calculates grass growth between trees and around the litter patterns and sets moisture levels based on relative humidity. In this study, we use DUET to simulate spatial distributions of leaf litter and illustrate model sensitivity to several parameters, providing inputs for testing fuel spatial heterogeneity impacts on fire behavior using process-based models.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Dust-planetary boundary layer interactions amplified by entrainment and advections

Mineral dust contributes to more than half of the global aerosol loading. However, the radiative impacts of dust aerosols on planetary boundary layer (PBL) structure have not been explored sufficiently. During a typical dust storm event over Tarim Basin, dust aerosols exhibit a well-mixed distribution during the daytime in spite of a shallow layer of dust particles accumulated at higher altitudes. By contrast, nocturnal dust plumes are located near the surface due to stable stratification. We demonstrate that these differentiated vertical distributions determine the spatial heterogeneity of dust loading, radiative fluxes and PBL height variations. Dust aerosols cause daytime PBL suppression and nighttime PBL promotion through modulating surface and atmospheric radiative budgets. Specifically, dust-induced cooling effect within PBL directly inhibits the daytime PBL development. PBL suppression effect is then amplified by entrainment processes resulting in excessively low PBL height, especially for dust particles below but near the PBL top. Dust plumes weaken updrafts from PBL and downdrafts of the free atmosphere, which further reduce the entrainment mixing through attenuating horizontal and vertical advection, and eventually amplify PBL suppression. At night, near-surface dust aerosols stimulate a warm and unstable lower atmosphere, generate warm advection heating and promote the PBL development. Our study highlights the importance of specifying entrainment parameters and nighttime advection activities in quantifying the dust-PBL interactions.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Geophysical methods reveal the soil architecture and subsurface stratigraphic heterogeneities across land-lake interfaces along Lake Erie

Abstract Purpose The land-lake interface is a unique zone where terrestrial and aquatic ecosystems meet, forming part of the Earth’s most geochemically and biologically active zones. The unique characteristics of this interface are yet to be properly understood due to the inherently high spatiotemporal variability of subsurface properties, which are difficult to capture with the traditional soil sampling methods. Geophysical methods offer non-invasive techniques to capture variabilities in soil properties at a high resolution across various spatiotemporal scales. Methods We combined electromagnetic induction (EMI), electrical resistivity tomography (ERT), and ground penetrating radar (GPR) with data from soil cores and in situ sensors to investigate hydrostratigraphic heterogeneities across land-lake interfaces along the western basin of Lake Erie. Results EMI revealed high spatial heterogeneities in ECa distribution across the land-lake interfaces, with higher values in the wetland and transition zones compared to the upland zone. Soil ECa maps matched soil maps from a public database with the hydric soil units delineated as high conductivity zones (ECa > 40 mS/m). ERT and GPR showed vertical variation in soil properties with clear stratigraphic boundaries, and correlation of ERT profiles with lithologs from piezometers revealed the stratigraphic units of silt–clay and till sequence down to 3.5 m depth which are consistent with the surficial geology of the study area. Conclusions These results validate the use of multiple geophysical methods for extrapolating soil properties and mapping stratigraphic structures at land-lake interfaces, thereby providing the missing information required to improve the earth system model (ESM) of coastal interfaces.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

A hierarchical stochastic modeling approach for representing point bar geometries and petrophysical property variations

The flow of fluids in point bars is affected by the existence of heterogeneities like shale drapes that are found on the surfaces of inclined heterolithic stratifications. In fact, these shale drapes can act as fluid flow baffles; therefore, developing a framework for modeling point bars and their associated heterogeneities is vital. In this study, a stochastic process-based modeling approach is presented for capturing the main point bar heterogeneities: accretion surfaces (i.e., the aerial heterogeneity) and inclined heterolithic stratifications (i.e., the vertical heterogeneity). The former was modeled using a sine-generation function and the latter, with a sigmoidal function after which they were combined into a 3D point bar model. To ensure proper modeling of petrophysical properties, we developed a more representative gridding scheme which generates curvilinear grids representative of the point bar geometry. This grid was then transformed into a rectilinear grid to allow for geostatistical simulation after which all petrophysical properties were mapped back into the original curvilinear grid. An essential element of this modeling approach is the stochastic representation of shale drapes at the interface between successive accretion surfaces. The workflow was tested using a real field dataset for the Cranfield field, Mississippi. The constructed model was then subjected to a flow simulation study mimicking a CO 2 storage scenario. Various sensitivities were simulated to evaluate the effect of heterogeneities on CO 2 flow within the point-bar. Results demonstrate the importance of representing point-bar related heterogeneity and the spatial distribution of shale drapes on CO 2 plume migration and storage.

58 GEOSCIENCES↗
DOE OSTI2021

Distributed Solar Adoption in Orlando: A Household-Level Model for Distribution Resource Planning

Potential for rooftop solar in Florida is massive (47% of retail sales, 3rd overall nationally), yet adoption lags (12th nationally). A 2018 Florida Public Service Commission ruling authorizing solar third-party ownership (leasing) has substantially increased attention on distributed solar in the state. The city of Orlando has committed to a 100% clean-energy target by 2050 and deployment of solar and storage are expected to contribute significantly to reaching the goal. Deployment of customer-adopted solar, unlike utility-procured solar, is uncertain, but known to be spatially correlated with demographic factors and existing adoption. We develop a new method to adapt NREL's dGen model in order to represent building-level agents in adoption forecasts for the Orlando Utility Commission (OUC) service territory. Using the agent-based model we develop projections of solar adoption, subject to scenarios varying future solar costs and valuation, and aggregate adoption predictions by OUC distribution feeder. We find substantial spatial heterogeneity in the projected level of adoption by OUC distribution feeder. For instance, 25% of all projected adoption through 2050 would be concentrated on just 5% of feeders and 88% of projected adoption on 50% of feeders. Because of the uncertainty in adoption, bottoms-up solar adoption forecasting methods at the household-level are integral to long-term resource planning by anticipating system needs as customers increasingly adopt distributed solar, storage, electric vehicles, and other distributed energy resources.

14 SOLAR ENERGY↗
DOE OSTI2026

Comparison of automated chemical-guided segmentation and human annotation of soil organic matter in X-ray microcomputed tomography imaging in contrasted soil types

Soil organic matter (OM) formation and persistence is strongly influenced by the spatial distribution of organic substrates and microscale soil heterogeneity by dictating OM accessibility to microorganisms. However, traditional size and/or density fractionation techniques disrupt aggregate architecture, eliminating spatial information needed to fully understand intra-aggregate OM distribution. To quantify three-dimensional OM spatial distribution and automate segmentation in X-ray microcomputed tomography (µCT) imaging without human annotation bias, we developed an iodine gas vapor (I2) based staining workflow that eliminates labor-intensive manual annotation while maintaining segmentation accuracy, using aggregates from four taxonomically diverse soils (Xerofluvent, Haploxeroll Sphagnofibrist, Palehumult) with an 8-fold range of soil organic carbon. Human annotation of 10 µCT slices by the experienced and inexperienced annotators resulted in variations up to 3% in the Dice similarity coefficient (DSC), reflecting a degree of inherent subjectivity of manual labeling. Such inconsistencies are expected to compound as the number of manually annotated slices increases. Dual-energy µCT imaging at 33.1 keV (below the iodine (I) K-edge) and 33.2 keV (above the I K-edge) was used to resolve aggregate microstructure following I2 staining. The automated image subtraction pipeline identified OM regions by the I Kedge induced brightness increases, achieving DSC values of 0.58–0.83 relative to an experienced annotator. Sensitivity analyses revealed that the reconstruction alpha value—optimized via the open-source tool TomocuPy—and the 3D registration slice count were the primary determinants of accuracy, providing a novel benchmark for dual-energy soil imaging. The pipeline without GPU acceleration achieved 9.6 to 43.2 times faster than manual annotation. Using GPU-accelerated image post-processing and affine transformation matrices, the pipeline successfully segmented OM elements for large-scale datasets (3232×3232 pixel, 2048 slices) within ~5200 s from raw file acquisition to segmented output. The high-throughput approach enables the quantification of OM spatial distribution across diverse and heterogeneous soil.

Soil microbial biomass↗
DOE OSTI2021

Sub-10 nm Probing of Ferroelectricity in Heterogeneous Materials by Machine Learning Enabled Contact Kelvin Probe Force Microscopy

Reducing the dimensions of ferroelectric materials down to the nanoscale has strong implications on the ferroelectric polarization pattern and on the ability to switch the polarization. As the size of ferroelectric domains shrinks to the nanometer scale, the heterogeneity of the polarization pattern becomes increasingly pronounced, enabling a large variety of possible polar textures in nanocrystalline and nanocomposite materials. Critical to the understanding of fundamental physics of such materials and hence their applications in electronic nanodevices is the ability to investigate their ferroelectric polarization at the nanoscale in a nondestructive way. We show that contact Kelvin probe force microscopy (cKPFM) combined with a k-means response clustering algorithm enables to measure the ferroelectric response at a mapping resolution of 8 nm. In a BaTiO 3 thin film on silicon composed of tetragonal and hexagonal nanocrystals, we determine a nanoscale lateral distribution of discrete ferroelectric response clusters, fully consistent with the nanostructure determined by transmission electron microscopy. Moreover, we apply this data clustering method to the cKPFM responses measured at different temperatures, which allows us to follow the corresponding change in the polarization pattern as the Curie temperature is approached and across the phase transition. This work opens up perspectives for mapping complex ferroelectric polarization textures such as curled/swirled polar textures that can be stabilized in epitaxial heterostructures and more generally for mapping the polar domain distribution of any spatially highly heterogeneous ferroelectric materials.

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

Microscale spatial distribution and soil organic matter persistence in top and subsoil

The spatial distribution of organic substrates and microscale soil heterogeneity significantly influence organic matter (OM) persistence as constraints on OM accessibility to microorganisms. However, it is unclear how changes in OM spatial heterogeneity driven by factors such as soil depth affect the relative importance of substrate spatial distribution on OM persistence. This work evaluated the decomposition and persistence of 13 C and 15 N labeled water-extractable OM inputs over 50 days as either hotspot (i.e., pelleted in 1–2 mm-size pieces) or distributed (i.e., added as OM < 0.07 μm suspended in water) forms in topsoil (0–0.2 m) and subsoil (0.8–0.9 m) samples of an Andisol. We observed greater persistence of added C in the subsoil with distributed OM inputs relative to hotspot OM, indicated by a 17% reduction in cumulative mineralization of the added C and a 10% higher conversion to mineral-associated OM. A lower substrate availability potentially reduced mineralization due to OM dispersion throughout the soil. NanoSIMS (nanoscale secondary ion mass spectrometry) analysis identified organo-mineral associations on cross-sectioned aggregate interiors in the subsoil. On the other hand, in the topsoil, we did not observe significant differences in the persistence of OM, suggesting that the large amounts of particulate OM already present in the soil outweighed the influence of added OM spatial distribution. Here, we demonstrated under laboratory conditions that the spatial distribution of fresh OM input alone significantly affected the decomposition and persistence of OM inputs in the subsoil. On the other hand, spatial distribution seems to play a lower role in topsoils rich in particulate OM. The divergence in the influence of OM spatial distribution between the top and subsoil is likely driven by differences in soil mineralogy and OM composition.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

Joint impedance and facies inversion of time-lapse seismic data for improving monitoring of CO 2 incidentally stored from CO 2 EOR

Time-lapse seismic monitoring is an effective and proven technology for mapping the distribution of CO 2 in a subsurface reservoir. When injected CO 2 displaces other reservoir fluids, porous-medium properties are changed and thus the seismic impedance changes, causing time-lapse seismic amplitude differences in the injection zones. The analysis and interpretation of images created from these amplitude differences can provide information about reservoir architecture and the CO 2 migration within the reservoir. Incorporating seismic inversion and rock physics into the interpretation of time-lapse seismic data can considerably improve the modeling and monitoring to detect and assess the location of CO 2 over time. The joint inversion method presented in this paper has an integral representation of the geology in the inversion algorithm using elastic facies, which provides information about the spatial distribution of the geologic heterogeneities controlling the movement of fluids in the reservoir. The method was successfully applied to time-lapse seismic data from a mature oil field undergoing CO 2 enhanced oil recovery. The estimated seismic acoustic impedances and facies reflect the characteristics of individual geologic facies and fluid conditions of the reservoir subject to CO 2 injection. Finally, the probabilities estimated by the joint impedance and facies inversion for the reservoir's litho-fluid facies can be used for forecasting CO 2 saturation and pressure changes within the target reservoir.

4D seismic↗