Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “spatial heterogeneity”

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 91 records · Page 5
All sourcesNASA NTRSDOE OSTI
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↗
NASA NTRS2005

Isotopic Investigations of Nebular and Parent Body Processes with a High Sensitivity Ion Microprobe

NASA supported the development of the CAMECA ims 1270 ion microprobe at UCLA for applications in cosmochemistry. The primary investigations centered on measuring the microscopic distributions of key isotopic abundances in primitive meteoritic materials as a means of constraining the nature of important thermal and chemical processes in the solar nebula and the timescales associated with those processes. Our prior work on oxygen isotope anomalies in a wide variety of meteoritic materials had led us to a view of a spatially heterogeneous nebula, and in particular, a restricted region for CAI formation that is characterized by O-16-rich gas. Because of its production of CAIs in the energetic local environment near the protosun, the existence of a natural transport mechanism via bipolar outflows, and a general astrophysical plausibility, we were attracted to the fluctuating X-wind model which had been put forward by Frank Shu, Typhoon Lee, and colleagues. With our collaborators, we undertook a series of investigations to test the viability of this hypothesis; this work led directly to the discovery of live Be in CAIs and a clear demonstration of the existence of 160-rich condensates, which necessarily implies an O-16-rich gaseous reservoir in the nebula. Both of these observations fit well within the context of X-wind type models, i.e. formation of CAIs (or condensation of their precursors) in the reconnection ring sunward of the inner edge of the accretion disk, however much work remains to be done to test whether the physical parameters of the model can quantitatively predict not only the thermal histories of CAIs but also their radioactivity. The issue of spatial heterogeneity in the nebula, central to the X-wind model, is also at the heart of any chronology based on short-lived radioisotopes. In this work, we followed up on strong hints for presence of exireme:j: (53 day) short-lived Be-7, and have prepared a manuscript (in revision). We also measured A1-Mg systematics by a combined approach of high-precision multiple-collector SIMS analyses, traditional analyses on the UCLA ims 1270, and high-spatial resolution analyses using a NanoSIMS instrument. The data help to deconvolve effects due to partial resetting of the A1-Mg system by multiple thermal events. Finally, we initiated investigations related to nebular heterogeneity with a new initiative of in situ high-precision sulfur isotope analyses of sulfides from a wide variety of components of chondrites. The ultimate goal of all this work is to help develop a better understanding of the relationships between CAIs and chondrules, the astrophysical environments in which they formed, and the timescales of nebular processes. As detailed in Table 1, for the project period, 14 manuscripts were published and 17 abstracts were presented describing the work.

McKeegan, Kevin D.↗
NASA NTRSDate not supplied

Radiative Flux Measurements from ARISE: A Comparison with CERES Top-of-Atmosphere Radiative Fluxes

Uncertainty in top-of-atmosphere (TOA) radiation fluxes observations are larger in the Arctic than in other regions. These uncertainties are due to the low sun angles and the highly reflective, anisotropic, and heterogeneous surface conditions. Quantifying, attributing, and reducing Arctic TOA radiative flux uncertainty enables a better understanding of the rapidly changing Arctic. To advance this goal, we compare the Cloud and Earth’s Radiant Energy System (CERES) TOA radiative fluxes with Arctic Radiation-IceBridge Sea and Ice Experiment (ARISE) campaign measurements collected in September 2014. We compare CERES TOA and aircraft radiative flux measurements using two complementary approaches: grid box average fluxes and instantaneously matched footprints. The grid box mean flux comparison indicates an agreement between CERES and aircraft measurements within 2 uncertainty (calibration and inversion) in the longwave for all five grid boxes and for four-of-five grid boxes in the shortwave; shortwave and longwave mean differences are -7.9 and +2.3 Wm‑2, respectively. The comparison of 36 instantaneously matched footprints with aircraft measurements reveals mean differences of -10.5 and 0.4 Wm‑2 in the shortwave and longwave, respectively. To further explore the persistent negative difference in the shortwave, we further quantify the effects of temporal and spatial sampling differences, angular distribution models, and scene identification to CERES-aircraft differences. Our analysis indicates that sampling differences (including scene evolution) account for an additional 1.8 and 1.7% uncertainty in the shortwave and longwave, respectively and indicates no bias. After accounting for this sampling uncertainty, all CERES-aircraft grid box mean fluxes agree within 2 uncertainty. Scene identification errors due to sea ice concentration data set differences exhibit no bias in the shortwave flux difference and indicate the possibility of substantial differences in the CERES fluxes in specific cases with large spatial heterogeneity. Considering the instantaneously matched footprints, we find that the angular distribution models account may account for up to ‑7.3 Wm-2 of the persistent CERES-aircraft shortwave flux difference due to systematic differences in the anisotropy for sea ice partly cloudy scenes. Additional analysis using a special scan model with the CERES FM2 instrument suggests a significant view zenith angle dependence of the CERES fluxes for sea ice partly cloudy scenes where shortwave fluxes systematically decrease with increasing view zenith angle; no dependence is found for other scene types. We conclude that (1) spatial heterogeneity and scene temporal evolution substantially limit our ability to use aircraft measurements to place strong constraints on CERES TOA fluxes and (2) that the representation of anisotropy in sea ice partly cloudy scenes is a significant factor contributing to the persistent negative CERES-aircraft shortwave flux difference in this comparison and require additional data to analysis fully quantify the potential bias.

Patrick C Taylor↗
NASA NTRSDate not supplied

A comparison between CERES TOA Radiative Fluxes and Airborne Radiative Flux Measurements from ARISE

Uncertainty in top-of-atmosphere (TOA) radiation fluxes observations are larger in the Arctic than in other regions. These uncertainties are due to the low sun angles and the highly reflective, anisotropic, and heterogeneous surface conditions. Quantifying, attributing, and reducing Arctic TOA radiative flux uncertainty enables a better understanding of the rapidly changing Arctic. To advance this goal, we compare the Cloud and Earth’s Radiant Energy System (CERES) TOA radiative fluxes with Arctic Radiation-IceBridge Sea and Ice Experiment (ARISE) campaign measurements collected in September 2014. We compare CERES TOA and aircraft radiative flux measurements using two complementary approaches: grid box average fluxes and instantaneously matched footprints. The grid box mean flux comparison indicates an agreement between CERES and aircraft measurements within 2σ uncertainty (calibration and inversion) in the longwave for all five grid boxes and for four-of-five grid boxes in the shortwave; shortwave and longwave mean differences are -7.9 and +2.3 Wm 2, respectively. The comparison of 36 40 instantaneously matched footprints with aircraft measurements reveals mean differences of -120.25 and -1.00.4 Wm 2 in the shortwave and longwave, respectively. To further explore the persistent negative difference in the shortwave, Wwe further quantify the effects of temporal and spatial sampling differences, angular distribution models, and scene identification to CERES-aircraft differences. Our analysis indicates that sampling differences (including scene evolution) account for an additional 1.8 and 1.7% uncertainty in the shortwave and longwave, respectively and , but indicates no bias. After accounting for this sampling uncertainty, all CERES-aircraft grid box mean fluxes agree within 2σ uncertainty. Scene identification errors due to sea ice concentration data set differences exhibit no bias in the shortwave flux difference and indicate the possibility of substantial differences in the CERES fluxes in specific cases with large spatial heterogeneity. Considering the instantaneously matched footprints, we find that the angular distribution models account may account for up to 7.3 Wm-2 of the persistent CERES-aircraft shortwave flux difference due to systematic differences in the anisotropy for sea ice partly cloudy scenes. Additional analysis using a special programmable scan model with the CERES FM2 instrument suggests a significant view zenith angle dependence of the CERES fluxes for sea ice partly cloudy scenes where shortwave fluxes systematically decrease with increasing view zenith angle; no dependence is found for other scene types. We conclude that (1) spatial heterogeneity and scene temporal evolution substantially limit our ability to use aircraft measurements to place strong constraints on CERES TOA fluxes and (2) that the representation of anisotropy in sea ice partly cloudy scenes is likely a significant factor contributing to the persistent negative CERES-aircraft shortwave flux difference in this comparison and require additional data to analysis fully quantify the potential bias.

Patrick C. Taylor↗
DOE OSTI2025

Future Intensity‐Duration‐Frequency Curves of Extreme Precipitation in the Midwest United States From Convection‐Permitting Modeling

Abstract During the last four decades, global warming has statistically significant intensified extreme precipitation events in the Midwestern United States (defined here as the region covering Illinois, Indiana, Ohio, and Kentucky), leading to increased risks to human life, property, and infrastructure. To enable climate change adaptation and resilience across various economic and social sectors in this region, updated information about future climate changes, specifically at finer spatial scales, is essential. Leveraging a new 150‐year dynamical downscaling data set at convection‐permitting resolution, this study introduces a framework to construct the projected future intensity‐duration‐frequency (IDF) curves of heavy precipitation, which are prominent tools for infrastructure design and water resources management. This framework generates IDF curves at both sub‐daily and multi‐day duration utilizing hourly in situ observations as well as quantile‐based statistical techniques in bias‐correction and return levels selection. The assumption of non‐stationarity in the distribution parameter fitting process is also implemented in this workflow. Compared to historical IDF curves for 1980–2022, future projected IDF curves for 2058–2100 under Representative Concentration Pathway (RCP) 4.5 and RCP 8.5 scenarios indicate an average intensity increase of approximately 15% and 25%, respectively, across 74 stations, considering both annual and seasonal timescales. Future projections suggest that extreme precipitation events may become more severe across six investigated return periods, with longer return periods showing a greater increase. The frequency of future extreme precipitation events in the Midwest region is also projected to double. Furthermore, current results reveal spatial heterogeneity of future trends across stations owing to the high‐resolution input data set. Plain Language Summary This study investigates the evolving nature of extreme precipitation events in the Midwestern United States under a changing climate. By leveraging a high‐resolution dynamical downscaling data set, we construct projected intensity‐duration‐frequency (IDF) curves for future extreme rainfall events. These curves serve as vital tools for infrastructure planning and water resource management. Our analysis reveals a significant increase in both the intensity and frequency of extreme precipitation events in the region. Future projected IDF curves for the late century indicate an average intensity increase of approximately 15%–25% compared to historical values. Moreover, the frequency of such events is expected to double. Spatial heterogeneity in future trends is observed across different stations within the Midwest, highlighting the importance of high‐resolution modeling in capturing localized climate variability. These findings underscore the urgent need for climate adaptation strategies to mitigate the increasing risks associated with extreme precipitation events in the region. Key Points This study introduces a workflow to construct future intensity‐duration‐frequency (IDF) curves over the Midwest United States using a new convection‐permitting modeling data set The current IDF construction workflow reproduces well the historical observed IDF 30 curves in summer months with median relative errors of 2.4% among 74 stations and 6 investigated durations The projected IDF curves show diverse future trends of extreme precipitation across stations, with intensity increases of approximately 15% and 25% under RCP4.5 and RCP8.5 climate scenarios, respectively, and a doubling of frequency on average

Nguyen, Trung↗
DOE OSTI2022

Parameterizing Convective Organization Effects With a Moisture-PDF Approach in Climate Models: Concept and a Regional Case Simulation

We propose a parameterization scheme of convective organization effects based on a moisture-distribution approach, which can reflect aggregation of convective cells within a model grid as well as the interaction between convection and spatial heterogeneity in free-troposphere moisture. With this concept, convective cells in an aggregated state are surrounded by air that is moister than the grid-mean condition, which provides a “shielding” effect favorable for deeper convective updrafts. Such effects are represented quantitatively via utilizing a schematic diagram about the geometric interpretation of convective cluster, dry area, and their contact area, in which both the convective cluster size and dry area fraction increase during the development of convection. Our sensitivity analyses indicate that the new scheme performs well in capturing the mean precipitation features. Excluding the convective organization effects in the model leads to a considerable reduction in the simulated precipitation magnitude. Notably, there exists an inverse relationship between background mean moisture condition and sub-grid moisture variability for a given rain rate, suggesting that with a high spatial heterogeneity in free-troposphere moisture, convection can still develop and maintain its strength under relatively dry background conditions, consistent with the convection-permitting model simulation. Overall, our new parameterization of convective organization effects can successfully reproduce the relationship between precipitation and sub-grid moisture variability, which is a missing element in traditional convection parameterization schemes but important for the simulations of precipitation variability at various scales in climate models.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Challenges in quantifying unparameterized spatial uncertainties in deep geologic repositories for nuclear waste

Spatially heterogeneous uncertainties are prevalent in geophysical modeling applications, such as probabilistic post-closure performance assessment (PA) of deep geologic repositories for nuclear waste. Such uncertainties are often highly influential to model outputs, so it is desirable to identify the most important mechanisms by which they influence model predictions. However, these uncertainties are often unparameterized in the sense that there is no set of parameters that can be specified to yield a particular realization of the uncertainty. Additionally, the uncertainty is not intrinsically endowed with a parameterization that captures a realization’s mechanistic influence on model outputs. Therefore, in this work we present a novel methodology to develop and assess a set of proxy variables that aim to represent this influence. We show how they can be computed, downselected, and incorporated into the construction of statistical surrogate models mapping model inputs to outputs. We present our methodology in the context of a motivating application problem in deep geologic repository PA and discuss the challenges in capturing the effects of these spatial heterogeneities in uncertainty analyses.

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

Nucleation kinetics model for primary crystallization in Al–Y–Fe metallic glass

The high density of aluminum nanocrystals (>10 21 m -3 ) that develop during the primary crystallization in Al-based metallic glasses indicates a high nucleation rate (~10 18 m -3 s -1 ). Several studies have been advanced to account for the primary crystallization behavior, but none have been developed to completely describe the reaction kinetics. Recently, structural analysis by fluctuation electron microscopy has demonstrated the presence of the Al-like medium range order (MRO) regions as a spatial heterogeneity in as-spun Al 88 Y 7 Fe 5 metallic glass that is representative for the class of Al-based amorphous alloys that develop Al nanocrystals during primary crystallization. From the structural characterization, an MRO seeded nucleation configuration is established, whereby the Al nanocrystals are catalyzed by the MRO core to decrease the nucleation barrier. The MRO seeded nucleation model and the kinetic data from the delay time (τ) measurement provide a full accounting of the evolution of the Al nanocrystal density (N v ) during the primary crystallization under isothermal annealing treatments. Moreover, the calculated values of the steady state nucleation rates (J ss ) predicted by the nucleation model agree with the experimental results. Moreover, the model satisfies constraints on the structural, thermodynamic, and kinetic parameters, such as the critical nucleus size, the interface energy, and the volume-free energy driving force that are essential for a fully self-consistent nucleation kinetics analysis. In conclusion, the nucleation kinetics model can be applied more broadly to materials that are characterized by the presence of spatial heterogeneities.

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

Gas Production from Hot Water Circulation through Hydraulic Fractures in Methane Hydrate-Bearing Sediments: THC-Coupled Simulation of Production Mechanisms

Methane hydrates, widely found in permafrost and deep marine sediments, have great potential as a future energy 6 source. Conventional production schemes perform poorly for challenging hydrate reservoirs with low permeability. We propose an 7 efficient production scheme by combining hydraulic fracturing from horizontal wells and hot water circulation through fractures. A 8 fully coupled thermo-hydro-chemical (THC) model is developed to simulate the key physical processes during gas production from 9 a hydrate reservoir representative of typical geological settings in Shenhu, South China Sea. We found that the gas production 10 process has two distinct stages divided by thermal breakthrough: a relatively short prebreakthrough stage and a postbreakthrough 11 stage yielding stable gas production. Heat advection along and near the hydraulic fracture dominates the prebreakthrough stage, 12 whereas conduction-driven thermal recovery in the volume around fractures dominates the postbreakthrough stage. We identified 13 that the steady-state injection temperature has a strong effect on the performance of the proposed scheme while the fluid mass 14 circulation rate has a moderate impact beyond a threshold. The proposed scheme proves to be efficient and robust over a range of 15 reservoir conditions with respect to initial hydrate saturation and intrinsic permeability, including their spatial heterogeneities, 16 thereby offering a promising solution for challenging reservoir conditions.

gas hydrate, fracture flow, gas production, numeri↗
DOE OSTI2023

Sensitivities of subgrid-scale physics schemes, meteorological forcing, and topographic radiation in atmosphere-through-bedrock integrated process models: a case study in the Upper Colorado River basin

Abstract. Mountain hydrology is controlled by interacting processes extending from the atmosphere through the bedrock. Integrated process models (IPMs), one of the main tools needed to interpret observations and refine conceptual models of the mountainous water cycle, require meteorological forcing that simulates the atmospheric process to predict hydroclimate then subsequently impacts surface–subsurface hydrology. Complex terrain and extreme spatial heterogeneity in mountainous environments drive uncertainty in several key considerations in IPM configurations and require further quantification and sensitivity analyses. Here, we present an IPM using the Weather Research and Forecasting (WRF) model which forces an integrated hydrologic model, ParFlow-CLM, implemented over a domain centered over the East River watershed (ERW), located in the Upper Colorado River basin (UCRB). The ERW is a heavily instrumented 300 km2 region in the headwaters of the UCRB near Crested Butte, CO, with a growing atmosphere-through-bedrock observation network. Through a series of experiments in the water year 2019 (WY19), we use four meteorological forcings derived from commonly used reanalysis datasets, three subgrid-scale physics scheme configurations in WRF, and two terrain shading options within WRF to test the relative importance of these experimental design choices for key hydrometeorological metrics including precipitation and snowpack, as well as evapotranspiration, groundwater storage, and discharge simulated by the ParFlow-CLM. Our hypothesis is that uncertainty from synoptic-scale forcings produces a much larger spread in surface–subsurface hydrologic fields than subgrid-scale physics scheme choice. Results reveal that the WRF subgrid-scale physics configuration leads to larger spatiotemporal variance in simulated hydrometeorological conditions, whereas variance across meteorological forcing with common subgrid-scale physics configurations is more spatiotemporally constrained. Despite reasonably simulating precipitation, a delay in simulated discharge peak is due to a systematic cold bias across WRF simulations, suggesting the need for bias correction. Discharge shows greater variance in response to the WRF simulations across subgrid-scale physics schemes (26 %) rather than meteorological forcing (6 %). The topographic radiation option has minor effects on the watershed-average hydrometeorological processes but adds profound spatial heterogeneity to local energy budgets (±30 W m−2 in shortwave radiation and 1 K air temperature differences in late summer). This is the first presentation of sensitivity analyses that provide support to help guide the scientific community to develop observational constraints on atmosphere-through-bedrock processes and their interactions.

54 ENVIRONMENTAL SCIENCES↗
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 OSTI2024

Polysulfide-incompatible additive suppresses spatial reaction heterogeneity of Li-S batteries

Rational electrolyte engineering for practical pouch cells remains elusive because the correlation between the cathode/solid-electrolyte interphase layer and cell-level reaction behavior is poorly understood. Here, by combining multiscale characterization and computational modeling, we show that—counter to the conventional perception of polysulfide-incompatible additives—the spontaneous reaction of sparingly solvated polysulfides with Lewis acid additives (LAAs) can induce in situ formation of a homogeneous interphase on thick and tortuous S cathode. Multiscale synchrotron X-ray characterization consistently affirms that such interface design could effectively eliminate the notorious problems of polysulfide shuttle and lithium corrosion and, more importantly, provide an interconnected “ion transport highway” to alleviate the uneven ion transport within the tortuous S cathode. Hence, this design dramatically reduces the reaction heterogeneity of lithium-sulfur (Li-S) pouch cells under lean electrolyte conditions. Further, this work resolves controversy around the role of polysulfide-incompatible additives in high-energy Li-S pouch cells and highlights the importance of suppressing reaction heterogeneity for practical batteries.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Electric-Field-Induced Spatially Dynamic Heterogeneity of Solvent Motion and Cation Transference in Electrolytes

While electric fields primarily result in migration of charged species in electrolytic solutions, the solutions are dynamically heterogeneous. Solvent molecules within the solvation shells of the cation will be dragged by the field while free solvent molecules will not. Here we combine electrophoretic NMR measurements of ion and solvent velocities under applied electric fields with molecular dynamics simulations to interrogate different solvation motifs in a model liquid electrolyte. Measured values of the cation transference number $(t^0_+)$ agree quantitatively with simulation-based predictions over a range of electrolyte concentrations. Solvent-cation interactions strongly influence the concentration-dependent behavior of $t^0_+$. We identify a critical concentration at which most of the solvent molecules lie within solvation shells of the cations. The dynamic heterogeneity of solvent molecules is minimized at this concentration where $t^0_+$ is approximately equal to zero.

42 ENGINEERING↗
DOE OSTI2021

Nutrient availability and senescence spatially structure the dynamics of a foundation species

Significance Landscapes are spatially heterogeneous, and the measurement of spatial pattern is dependent upon observation scale. Understanding plant populations requires assessing their extrinsic interactions with the environment as well as intrinsic biological processes and has been difficult because of the inability to track both plant abundance and health on appropriate scales. We introduce remote sensing observations that assess the abundance and health of giant kelp, an important ecosystem-structuring species, over regional and local scales. We find that both extrinsic nutrient availability and intrinsic senescence processes regulate population dynamics but on regional and local scales, respectively. This suggests that future satellite missions will be able to assess plant abundance and health and their interactions with the environment on local to global scales.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Integrating NDVI-Based Within-Wetland Vegetation Classification in a Land Surface Model Improves Methane Emission Estimations

Earth system models (ESMs) are a common tool for estimating local and global greenhouse gas emissions under current and projected future conditions. Efforts are underway to expand the representation of wetlands in the Energy Exascale Earth System Model (E3SM) Land Model (ELM) by resolving the simultaneous contributions to greenhouse gas fluxes from multiple, different, sub-grid-scale patch-types, representing different eco-hydrological patches within a wetland. However, for this effort to be effective, it should be coupled with the detection and mapping of within-wetland eco-hydrological patches in real-world wetlands, providing models with corresponding information about vegetation cover. In this short communication, we describe the application of a recently developed NDVI-based method for within-wetland vegetation classification on a coastal wetland in Louisiana and the use of the resulting yearly vegetation cover as input for ELM simulations. Processed Harmonized Landsat and Sentinel-2 (HLS) datasets were used to drive the sub-grid composition of simulated wetland vegetation each year, thus tracking the spatial heterogeneity of wetlands at sufficient spatial and temporal resolutions and providing necessary input for improving the estimation of methane emissions from wetlands. Our results show that including NDVI-based classification in an ELM reduced the uncertainty in predicted methane flux by decreasing the model’s RMSE when compared to Eddy Covariance measurements, while a minimal bias was introduced due to the resampling technique involved in processing HLS data. Our study shows promising results in integrating the remote sensing-based classification of within-wetland vegetation cover into earth system models, while improving their performances toward more accurate predictions of important greenhouse gas emissions.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

A hybrid data–model approach to map soil thickness in mountain hillslopes

Abstract. Soil thickness plays a central role in the interactions between vegetation, soils, and topography, where it controls the retention and release of water, carbon, nitrogen, and metals. However, mapping soil thickness, here defined as the mobile regolith layer, at high spatial resolution remains challenging. Here, we develop a hybrid model that combines a process-based model and empirical relationships to estimate the spatial heterogeneity of soil thickness with fine spatial resolution (0.5 m). We apply this model to two aspects of hillslopes (southwest- and northeast-facing, respectively) in the East River watershed in Colorado. Two independent measurement methods – auger and cone penetrometer – are used to sample soil thickness at 78 locations to calibrate the local value of unconstrained parameters within the hybrid model. Sensitivity analysis using the hybrid model reveals that the diffusion coefficient used in hillslope diffusion modeling has the largest sensitivity among all input parameters. In addition, our results from both sampling and modeling show that, in general, the northeast-facing hillslope has a deeper soil layer than the southwest-facing hillslope. By comparing the soil thickness estimated between a machine-learning approach and this hybrid model, the hybrid model provides higher accuracy and requires less sampling data. Modeling results further reveal that the southwest-facing hillslope has a slightly faster surface soil erosion rate and soil production rate than the northeast-facing hillslope, which suggests that the relatively less dense vegetation cover and drier surface soils on the southwest-facing slopes influence soil properties. With seven parameters in total for calibration, this hybrid model can provide a realistic soil thickness map with a relatively small amount of sampling dataset comparing to machine-learning approach. Integrating process-based modeling and statistical analysis not only provides a thorough understanding of the fundamental mechanisms for soil thickness prediction but also integrates the strengths of both statistical approaches and process-based modeling approaches.

58 GEOSCIENCES↗
DOE OSTI2024

Graph network heterogeneity predicts interplant wake losses

Wind plants generate large-scale wakes, which can affect the performance of neighboring installations. Such wakes are challenging to estimate due to the inherent complexity in modeling wake interactions between large quantities of turbines at various distances. Weighted directed graph networks can inform complex models by linking turbine pairs into chains of upstream and downstream neighbors for a given wind direction. A novel interpretation of the graph network adjacency matrix is proposed where each element of the matrix represents the cumulative impact of upstream turbines on an individual. In this study, wake losses were estimated with an engineering wake model across a range of inflow conditions for nine parametric variations of a system containing two neighboring wind plants. The parametric nature of the study isolates turbine spacing within the plant, separation distance between plants, and wind direction as the main drivers of wake losses. Spatial heterogeneity is computed from the weighted average adjacency matrix of each plant arrangement. The proposed method is orders of magnitude faster than wake modeling and does not require detailed turbine information or atmospheric conditions. Furthermore, the weighted average adjacency matrix provides insight on the spatial organization of wake losses at various scales. Plant heterogeneity is correlated with wake losses within and among plants. Framing wind plant wake interaction in terms of graph network spatial heterogeneity provides an efficient approach for predicting wake losses within and among neighboring wind plants with applications to other complex systems where wake interactions are key factor.

17 WIND ENERGY↗
DOE OSTI2024

Effects of spatial variability in vegetation phenology, climate, landcover, biodiversity, topography, and soil property on soil respiration across a coastal ecosystem

Coastal terrestrial-aquatic interfaces (TAIs) are crucial contributors to global biogeochemical cycles and carbon exchange. A systematic evaluation of the interaction between coastal catchment properties and carbon dioxide (CO2) emission by soil respiration is significant for assessing carbon dynamics and predicting the future trajectory of atmospheric CO2 concentrations in coastal TAIs. The soil CO2 efflux in these transition zones is however poorly understood due to the high spatiotemporal dynamics of TAIs, as various sub-ecosystems in this region are compressed and expanded by complex influences of tides, changes in river levels, climate, and land use. We focus on the Chesapeake Bay region to (i) investigate the spatial heterogeneity of the coastal ecosystem and identify spatial zones with similar environmental characteristics based on the spatial data layers, including vegetation index (kNDVI), climate, landcover, diversity, topography, soil property, and relative tidal elevation; (ii) understand the primary driving factors affecting soil respiration within sub-ecosystems of the coastal ecosystem. Specifically, we employed hierarchical clustering analysis to identify spatial regions with distinct environmental characteristics, followed by the determination of main driving factors using Random Forest regression and SHapley Additive exPlanations. Maximum and minimum temperature are the main drivers common to all sub-ecosystems, while each region also has additional unique major drivers that differentiate them from one another. Precipitation exerts an influence on vegetated lands, while soil pH value holds importance specifically in forested lands. In croplands characterized by high clay content and low sand content, the significant role is attributed to bulk density. Wetlands demonstrate the importance of both elevation and sand content, with clay content being more relevant in non-inundated wetlands than in inundated wetlands. The topographic wetness index significantly contributes to the mixed vegetation areas, including shrub, grass, pasture, and forest. Additionally, our research reveals that dense vegetation land covers and urban/developed areas exhibit distinct soil property drivers. Overall, there is no one-size-fits-all approach to modeling carbon fluxes in coastal TAIs, and our study highlights the importance of further research and monitoring practices to improve our understanding of carbon dynamics and promote the sustainable management of coastal TAIs.

54 ENVIRONMENTAL SCIENCES↗