Search NASA⌕ Search

SEARCH · Search NASA

Results for “constrained estimation”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 361 records · Page 20
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Observational constraints on wet and dry deposition of black carbon to land surfaces (Final Report)

This project focused on providing observational constraints on dry deposition of black carbon and aerosol particles, and the extent to which wet and dry deposition impacts the lifetime of aerosols in the atmosphere. The proposed work tackled the following questions: To what extent dry versus wet deposition controls the removal of BC from the atmosphere, and thus atmospheric BC concentrations? Is BC dry deposition controlled by the same processes as non-refractory aerosol dry deposition? We successfully made eddy covariance flux measurements at the DOE Southern Great Plains site of size-resolved aerosol number fluxes with an ultra high sensitivity aerosol spectrometer (UHSAS) and of black carbon mass and number fluxes using a single particle soot photometer (SP2). We also collected precipitation samples for off-line measurement of rBC using a wet deposition collector. Using these measurements, we identified that wet deposition was the dominant control for black carbon lifetime, but that dry deposition was poorly constrained in models. We developed a new parameterization for particle dry deposition, and applied it to global models to demonstrate the importance of dry deposition in air pollution and radiative effect estimates.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

PV Module Design for Recycling Guidelines

The global growth of clean energy technology deployment will be inexorably followed by a parallel growth of end-of-life (EOL) products that bring both challenges and opportunities. Cumulatively, by 2050, estimates project 78 million tonnes of raw materials embodied in the mass of EOL photovoltaic (PV) modules. Owing partly to concern that the projected growth of clean energy technologies could become constrained by availability of raw materials, despite ongoing dematerialization efforts, significant attention under the umbrella of circular economy has been brought to recycling these technologies at EOL. Yet PV has not been designed with recycling at EOL in mind, and it presents challenges to returning embodied raw materials back to use in new products through recycling. This study aims to inform future designs to improve recyclability through synthesis of prior published works augmented by novel recommendations that result in a set of general design for recycling (DfR) guidelines, with a subset specific to PV modules. We further discuss how established trends in design of PV modules could affect recyclability. If adopted today, application of these DfR guidelines could help to mitigate tomorrow's resource scarcity, lower the barriers and cost for PV recycling, and enable a circular economy during the energy transition.

14 SOLAR ENERGY↗
DOE OSTI2023

Constraining neutrino oscillation and interaction parameters with the NOvA Near Detector and Far Detector data using Markov Chain Monte Carlo

This thesis reports a constraint of the neutrino oscillation parameters $\Delta m^{2}_{32}$, $\sin^2 \theta_{23}$, and $\delta_{CP}$ using the NuMI Off-Axis $\nu$ Appearance (NOvA) experiment's Near Detector (ND) data and Far Detector (FD) fake data set simultaneously. This thesis also reports a constraint on NOvA's systematic uncertainty model solely with its Near Detector data. The Hamiltonian Monte Carlo algorithm is used to estimate Bayesian Credible Intervals for the oscillation and interaction parameters. The $1\sigma$ Credible Intervals for $\sin^2 \theta_{23}$ are $(0.44, 0.512)$ $\cup$ $(0.536, 0.56)$, for $\Delta m^{2}_{32}$ $(2.41 \times 10^{-3}$ eV$^2,\ 2.52 \times 10^{-3}$ eV$^2)$, and for $\delta_{CP}$ $(0.74\pi,\ 1.1\pi)$ $\cup$ $(1.38\pi,\ 1.58\pi)$. The statistical power of the ND data constrains NOvA's interaction parameters, while the FD fake data constrains the oscillation parameters. This is the first analysis within NOvA to constrain the ND and FD prediction sim ultaneously, and to investigate the neutrino interaction modeling in the context of constraining the oscillation parameters. To constrain the ND data requires a sophisticated understanding of the neutrino interaction modeling and its uncertainties. The interested reader is advised to focus on Chapters 4 and 6, which discuss the ND selection, uncertainties, and ND-only fits to data. The reader interested in oscillation parameter constraints will find this in Chapter 7.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

Primordial black hole mass functions as a probe of cosmic origin

Here, we discuss a novel window to probe the origin of our universe via the mass functions of primordial black holes (PBHs). The mass functions of PBHs are simply estimated using the conventional Press-Schechter formalism for two paradigms of cosmic origin, including inflationary ΛCDM and bounce cosmology. The standard inflationary ΛCDM model cannot generate an appreciable number of massive PBHs; however, non-trivial inflation models with blue-tilted power spectra at small scales and matter bounce cosmology provide formation mechanisms for heavy PBHs, which in turn, may seed the observed supermassive black holes (SMBHs). By fitting the SMBH mass functions at high redshift (z ~ 6) derived from Sloan Digital Sky Survey (SDSS) and Canada-France High-z Quasar Survey (CFHQS) quasars, for two paradigms of cosmic origin, we derive constraints on the PBH density fraction $f$ PBH at z ~ 6 and the characteristic mass M $\ast$ , with the prior assumption that all SMBHs stem from PBHs. We demonstrate that this newly proposed procedure, relying on astronomical measurements that utilize deep-field surveys of SMBHs at high redshift, can be used to constrain models of cosmic origin. Additionally, although not the main focus of this paper, we evolve the mass function from z ~ 6 to z ~ 0 through an assumption of 3 × 10 8 -year Eddington’s accretion, and give a rough estimation of $f$ PBH at z ~ 0.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Concurrent multi-peak Bragg coherent x-ray diffraction imaging of 3D nanocrystal lattice displacement via global optimization

Abstract In this paper we demonstrated a method to reconstruct vector-valued lattice distortion fields within nanoscale crystals by optimization of a forward model of multi-reflection Bragg coherent diffraction imaging (MR-BCDI) data. The method flexibly accounts for geometric factors that arise when making BCDI measurements, is amenable to efficient inversion with modern optimization toolkits, and allows for globally constraining a single image reconstruction to multiple Bragg peak measurements. This is enabled by a forward model that emulates the multiple Bragg peaks of a MR-BCDI experiment from a single estimate of the 3D crystal sample. We present this forward model, we implement it within the stochastic gradient descent optimization framework, and we demonstrate it with simulated and experimental data of nanocrystals with inhomogeneous internal lattice displacement. We find that utilizing a global optimization approach to MR-BCDI affords a reliable path to convergence of data which is otherwise challenging to reconstruct.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Constraining primordial non-Gaussianity from the large scale structure two-point and three-point correlation functions

Surveys of cosmological large-scale structure (LSS) are sensitive to the presence of local primordial non-Gaussianity (PNG), and may be used to constrain models of inflation. Local PNG, characterized by f NL ⁠, the amplitude of the quadratic correction to the potential of a Gaussian random field, is traditionally measured from LSS two-point and three-point clustering via the power spectrum and bi-spectrum. We propose a framework to measure f NL using the configuration space two-point correlation function (2pcf) monopole and three-point correlation function (3pcf) monopole of survey tracers. Our model estimates the effect of the scale-dependent bias induced by the presence of PNG on the 2pcf and 3pcf from the clustering of simulated dark matter haloes. We describe how this effect may be scaled to an arbitrary tracer of the cosmological matter density. The 2pcf and 3pcf of this tracer are measured to constrain the value of f NL ⁠. In LSS surveys, the effect of imaging systematics on two-point statistics is often degenerate with the PNG signal. Our proposed model employs three-point statistics primarily to break this degeneracy. Using simulations of luminous red galaxies observed by the Dark Energy Spectroscopic Instrument (DESI), we demonstrate the accuracy and constraining power of our method. Our forecast indicates the ability to constrain f NL to a precision of σf NL ≈ 22 with one year of DESI survey data, as well as the ability to constrain the imaging systematic weights in situ.

early Universe↗
DOE OSTI2022

Load Margin Constrained Moving Target Defense against False Data Injection Attacks

Cyber physical security of power systems with high penetration of renewable generation has attracted attention from researchers. One critical issue is that cyber-physical attacks, disguised as uncertain renewable generation, can target conventional power system state estimation (SE). Moving target defense (MTD) is a promising defense strategy to detect stealthy false data injection (FDI) attacks against SE. However, all existing studies myopically perturb the reactance of transmission lines equipped with distributed flexible AC transmission system (D-FACTS) devices without adequately considering the system voltage stability. Exacerbated by the renewable generation uncertainty, existing MTD may cause voltage instability when the power grid is under stress. To address this issue, we propose a novel MTD framework that explicitly considers system voltage stability by using continuation power flow. We utilize the sensitivity matrix of power injection to line impedance, on which an optimization problem for maximizing load margin is formulated. This framework is validated on the IEEE 14-bus system and the IEEE 118-bus system, in which net load redistribution attacks are launched by sophisticated attackers. Steady-state simulations and dynamic simulations on PSS/E show the effectiveness of the proposed framework in circumventing the voltage instability while maintaining the detection effectiveness of MTD. The impact of the proposed method on attack detection effectiveness is also revealed.

Zhang, Hang↗
DOE OSTI2025

In-Situ Synchrotron X-Ray Diffraction of Ultrasonic Microstructural Refinement During Solidification in a Commercial Al–Si–Mg Alloy

This study reports the first use of in-situ synchrotron X-ray diffraction (SXRD) to study the effects of ultrasonic melt processing (USMP) on phase and grain size evolution during solidification in a commercial Al–Si–Mg casting alloy. USMP is a technique that, when applied to aluminum as it solidifies, can be used to refine the local microstructure of large-scale castings. Analysis of the in-situ SXRD data to estimate the average grain size of primary α-Al grains during USMP demonstrates that USMP slows the growth rate of α-Al grains and reduces grain size by 36 pct. Furthermore, there is also evidence that USMP causes the primary α-Al grains to move relative to the X-ray beam; such motion increases the probability of primary α-Al grains colliding and fragmenting. This movement becomes constrained at the onset of the Al–Si binary eutectic, suggesting that USMP ceases to effectively refine the microstructure once the Al–Si binary eutectic begins to form. Complementary laboratory-scale X-ray diffraction (XRD) data were used to correlate the lattice parameters of the α-Al and Si (D-A4) phases with temperature to estimate cooling rate during solidification. Thus, this study can guide the design of novel castings with spatially distributed fine-grained regions produced using local ultrasonic processing.

Aluminum Alloys↗
DOE OSTI2023

An ecologically-constrained deep learning model for tropical leaf phenology monitoring using PlanetScope satellites

In tropical forests, leaf phenology signals leaf-on/off status and exhibits considerable variability across scales from a single tree-crown to the entire forest ecosystem. Such phenology signals importantly regulate large-scale biogeochemical cycles and regional climate. PlanetScope CubeSats data with a 3-m resolution and near-daily global coverage provide an unprecedented opportunity to monitor both fine- and ecosystem-scale phenology variability along large environmental gradients. However, a scalable method that accurately characterizes leaf phenology from PlanetScope with biophysically meaningful metrics remains lacking. Here we developed an index-guided, ecologically constrained autoencoder (IG-ECAE) method to automatically derive a deciduousness metric (i.e. percentage of tree canopies with leaf-off status within an image pixel) from PlanetScope. The IG-ECAE first estimated the reflectance spectra of leafy/leafless canopies based on their spectral indices characteristics, then used the derived reflectance spectra to guide an autoencoder deep learning method with additional ecological constraints to refine the reflectance spectra, and finally used linear spectral unmixing to estimate the relative abundance of leafless canopies (or deciduousness) per PlanetScope image pixel. We tested the IG-ECAE method at 16 tropical forest sites spanning multiple continents and a large precipitation gradient (1470-2819 mm year -1 ). Among these sites, we evaluated the PlanetScope-derived deciduousness with corresponding measures derived from WorldView-2 (n=9 sites) and local phenocams (n=9 sites). Our results show that PlanetScope-derived deciduousness agrees: 1) with that derived from WorldView-2 at the patch level (90m×90m) with r 2 =0.89 across all sites; and 2) with that derived from phenocams to quantify ecosystem-scale seasonality with r 2 ranging from 0.62-0.96. These results demonstrate the effectiveness and scalability of IG-ECAE in characterizing the wide variability in deciduousness across scales from pixels to forest ecosystems, and from a single date to the full annual cycle, indicating the potential for using high-resolution satellites to track the large-scale phenological patterns and response of tropical forests to climate change.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Analytic Neural Network Gaussian Process Enabled Chance-Constrained Voltage Regulation for Active Distribution Systems with PVs, Batteries and EVs

This paper proposes an analytic neural network Gaussian process (NNGP)-based chance-constrained real-time voltage regulation method for active distribution systems with photovoltaics (PVs), batteries, and electric vehicles (EVs). NNGP can utilize historical measurement data to achieve real-time probabilistic node voltage estimation through Bayesian inference. Then, NNGP is fully analytically embedded into the optimal power flow model to perform voltage regulation and adapt to various topological changes. The uncertainties of voltage estimations are easily considered via the chance constraint, and it has been shown that the adoption of this chance constraint can significantly improve the reliability of voltage regulation under various scenarios. The comparison results with other methods, carried out on a real 759-node distribution system located in western Colorado, U.S., show that the proposed method can achieve accurate voltage estimation across different topologies and reliably perform voltage regulation considering PVs, batteries, and EVs.

active distribution systems↗
DOE OSTI2025

Partnership Center for High-Fidelity Boundary Plasma Simulation (Final Report)

Within the Partnership Center for High-Fidelity Boundary Plasma Simulation (HBPS), work at UT-Austin was aimed at improved verification, validation, and uncertainty quantification (VVUQ) for edge plasma simulations and on performing gyrokinetics simulations of pedestal instabilities and turbulence in order to expand foundational understanding of pedestal transport. Regarding VVUQ, the accomplishments can be summarized as follows. First, it was shown that the Moment Preserving Constrained Resampling technique, when applied periodically in particle-in-cell simulations in the XGC code, can dramatically improve the accuracy of the simulation at essentially equivalent computational cost. Second, a technique for estimating model correlations, which are required to solve the model selection and sample allocation problem in multifidelity UQ techniques, without sampling the highest fidelity, most computationally expensive model, was developed and demonstrated. Third, previously developed methods for estimating statistical and discretization errors were applied to numerical methods relevant to edge plasma simulations, namely in particle-in-cell-based approaches, and shown to work. Finally, benchmark studies for comparing gyrokinetic codes were developed and performed, leading to reasonable agreement between four commonly used codes. Regarding physics studies, gyrokinetic simulations to investigate microtearing modes in the DIII-D pedestal were performed using the GENE code.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Simulation of Physics-Based 0-10Hz Strong Motion Using High Performance Computing Supporting Refinements to Regional Ground Motion Models for the Central Eastern US

In collaboration with the U.S. Nuclear Regulatory Commission (NRC) the LLNL has developed a computationally efficient simulation platform designed to perform physics-based ground motion simulations for crustal earthquakes in the Stable Continental Regions of Central and Eastern US (CEUS), using high-performance computing. The main objective of the earthquake simulations was to use synthetic ground motion to provide constrains to refinements of existing ergodic Ground Motion Models (GMMs), for large magnitude earthquakes and near-fault distances, for which these models are less reliable. Physics-based broadband (0-10Hz) ground motion simulations were used to estimate the near-fault ground motion amplitudes and within event and between-event variabilities associated with fault rupture characteristics. In our simulations we used a 3D regional velocity model that was based on Saikia’s 1D velocity model (1994). In simulations performed during the first stage of this project the Saikia’s velocity model demonstrated better performance in modelling high frequency regional wave propagation for the CEUS region recorded during the Mw5.0 November 7, 2016, Cushing Oklahoma (Taylor et al., 2017), and Mw5.8 September 3, 2016, Pawnee Oklahoma earthquakes. The proposed regional 3D model includes random perturbations to the 1D background model using the stochastic scheme of Pitarka and Mellors (2021). In addition, validation analysis of the rupture generator and regional wave propagation models, using comparisons with different GMMs for Mw6.5 and Mw7.0 scenario earthquakes in the CEUS region resulted in a very good match between the simulated and empirical ground motion models. For the purposes of seismic hazard assessment at the existing and planned nuclear power plants, NRC is interested in studies aimed at improving the current ground motion models (GMM) for both Stable Continental Regions (SCR) in the Central and Eastern US and Active Crustal Regions (ACR) in the Western US. Due to lack of recorded data, these improvements require synthetic data for short fault distances and large magnitude earthquakes for which the existing recorded data is not enough to uniquely constrain the GMMs. The need for simulations and strong motion data is especially critical for the CEUS region where we do not have recorded data from potentially large damaging earthquakes with moment magnitudes 6.0 and higher. In this the project, we focused on 10Hz simulations of Mw7.0 scenario earthquakes with strike slip and thrust faulting mechanisms. We used more than 50 Mw7.0 earthquake rupture scenarios to investigate the ground motion uncertainty due to unknown earthquake rupture parameters, in particular, the slip distribution, rupture velocity, and faulting mechanism, and their implication on ground motion amplification due to forward rupture directivity effects.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2020

Where did the tumor start? An inverse solver with sparse localization for tumor growth models

In this work, we present a numerical scheme for solving an inverse problem for parameter estimation in tumor growth models for glioblastomas, a form of aggressive primary brain tumor. The growth model is a reaction–diffusion partial differential equation (PDE) for the tumor concentration. We use a PDE-constrained optimization formulation for the inverse problem. The unknown parameters are the reaction coefficient (proliferation), the diffusion coefficient (infiltration), and the initial condition field for the tumor PDE. Segmentation of magnetic resonance imaging (MRI) scans drive the inverse problem where segmented tumor regions serve as partial observations of the tumor concentration. Like most cases in clinical practice, we use data from a single time snapshot. Moreover, the precise time relative to the initiation of the tumor is unknown, which poses an additional difficulty for inversion. We perform a frozen-coefficient spectral analysis and show that the inverse problem is severely ill-posed. We introduce a biophysically motivated regularization on the structure and magnitude of the tumor initial condition. In particular, we assume that the tumor starts at a few locations (enforced with a sparsity constraint on the initial condition of the tumor) and that the initial condition magnitude in the maximum norm is equal to one. We solve the resulting optimization problem using an inexact quasi-Newton method combined with a compressive sampling algorithm for the sparsity constraint. Our implementation uses PETSc and AccFFT libraries. We conduct numerical experiments on synthetic and clinical images to highlight the improved performance of our solver over a previously existing solver that uses standard two-norm regularization for the calibration parameters. The existing solver is unable to localize the initial condition. Our new solver can localize the initial condition and recover infiltration and proliferation. In clinical datasets (for which the ground truth is unknown), our solver results in qualitatively different solutions compared to the two-norm regularized solver.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Simple self-consistent method for excited states in density functional theory to characterize defect-derived behavior in wide-band-gap-based microelectronic materials

This final report summarizes the results of the Laboratory Direct Research and Development (LDRD) Project Number 229740. Wide band gap semiconductors such as gallium nitride (GaN) have features highly desirable for multiple mission electronic applications. Realization of their potential requires atomic-scale understanding of electronic behavior. The principal experimental tools for electronically probing defects in GaN are chemically undifferentiating and lack a practical theoretical counterpart needed to identify and characterize specific defects. This project investigated whether a simple idea for modeling defect excited states and their associated photoluminescence (PL) energies is viable, as a path to accelerate the understanding of defect behavior and gain valuable insights into engineering new electronic materials and devices. The research implemented a non-self-consistent total-energy evaluation of a Koopmans-type estimation of an excited electronic state energy in density functional theory (DFT) calculations, and proceeded to design, implement, and assess a self-consistent method for computing excited states based upon an OCcupation-Constrained-DFT (occ-DFT). The occ-DFT was verified in test calculations of defect excited states and validated against well-characterized PL data for 3d transition metal defects in GaN. The method proved stable and robust in computing excited states and gave accurate predictions compared to experimental PL data. The combined ground state/excited-state capability proved capable of chemically differentiating defect species in GaN. In application to 3d dopants in GaN, we reinterpreted extensive experimental literature, proposed new defects as prospective candidates for use in quantum information applications, and outlined design strategies to create and exploit these potentially useful functional defects in GaN.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Large reductions in Permian Basin methane intensity shown in multi-year comparison of aerially-visible methane emissions

Spanning the US states of Texas and New Mexico, the Permian Basin has been a hotspot of methane emissions from oil and natural gas activity 1–4, although studies disagree over the magnitude of these emissions. The most comprehensive measurement campaigns published were conducted in 2019 1–5 .There have been large changes in the energy industry since then, including in the prices of oil and gas, both state and federal regulatory environments, investor and activist pressure over methane emissions, and the adoption of new technologies and policies by energy operators. Understanding how any or all of these might influence methane emissions is important for policy makers, oil and gas operators, and other stakeholders. We characterize the time evolution of Permian Basin methane emissions using a series of comprehensive aerial surveys conducted every year from 2020-2023 and compare them to the 2019 results cited above. To maintain comparability, all the data sets are from surveys using Insight M point source methane sensing technology. The scope of these surveys expanded over time: from 33-46% of wells, oil production, and gas production in 2020 to 60-65% in 2021, to 84% of wells and over 90% of both oil and gas production in 2023. These surveys by Insight M also include hundreds of gas processing plants and compressor stations as well as 1000s of km of gathering and transmission pipelines. Considering only the aerially detected portion of emissions (typically the majority of the total in such surveys 4), we find reductions of more than 70% in methane emissions intensity compared to the 2019 New Mexico-only Insight M survey, with variation depending on the year 3,4. Notably, although sources below 100 kg/hr contributed less than 10% of aerially measured emissions the 2019 New Mexico survey 3, these smaller sources constitute a larger proportion of total aerially measured emissions (although not the majority) in 2020-2023. Production facilities and gathering pipelines are responsible for the larges shares of total emissions, followed by compressor stations and gas processing plants. Permian methane emissions were also measured in a comprehensive 2019 Permian-wide survey by the Carbon Mapper team 2. That analysis led to a lower total emissions estimate at the time 4. These new Insight M-based emission rates are still roughly 30-70% lower than the aerially measured portion of the 2019 Carbon Mapper-based estimates 4. Further work is needed to harmonize these surveys in space and time to create the most intercomparable numbers possible 5. Additional analysis is needed to compare our findings to the more spatially constrained 2020, 2021, and 2023 Carbon Mapper surveys in the Permian 4,6. The evidence is strong from these two survey teams that emissions intensity has declined significantly since 2019. Reasons for this trend are currently unclear but point to possible success of emissions control programs. Future work investigating frequency, source, and operator-specific intensities could provide insights into the causes of this promising trend.

methane, oil and gas, data science, remote sensing↗
DOE OSTI2024

Deriving cloud droplet number concentration from surface-based remote sensors with an emphasis on lidar measurements

Abstract. Given the importance of constraining cloud droplet number concentrations (Nd) in low-level clouds, we explore two methods for retrieving Nd from surface-based remote sensing that emphasize the information content in lidar measurements. Because Nd is the zeroth moment of the droplet size distribution (DSD), and all remote sensing approaches respond to DSD moments that are at least 2 orders of magnitude greater than the zeroth moment, deriving Nd from remote sensing measurements has significant uncertainty. At minimum, such algorithms require the extrapolation of information from two other measurements that respond to different moments of the DSD. Lidar, for instance, is sensitive to the second moment (cross-sectional area) of the DSD, while other measures from microwave sensors respond to higher-order moments. We develop methods using a simple lidar forward model that demonstrates that the depth to the maximum in lidar-attenuated backscatter (Rmax⁡) is strongly sensitive to Nd when some measure of the liquid water content vertical profile is given or assumed. Knowledge of Rmax⁡ to within 5 m can constrain Nd to within several tens of percent. However, operational lidar networks provide vertical resolutions of > 15 m, making a direct calculation of Nd from Rmax⁡ very uncertain. Therefore, we develop a Bayesian optimal estimation algorithm that brings additional information to the inversion such as lidar-derived extinction and radar reflectivity near the cloud top. This statistical approach provides reasonable characterizations of Nd and effective radius (re) to within approximately a factor of 2 and 30 %, respectively. By comparing surface-derived cloud properties with MODIS satellite and aircraft data collected during the MARCUS and CAPRICORN II campaigns, we demonstrate the utility of the methodology.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Constraining effects of aerosol-cloud interaction by accounting for coupling between cloud and land surface

Aerosol-cloud interactions (ACIs) are vital for regulating Earth’s climate by influencing energy and water cycles. Yet, effects of ACI bear large uncertainties, evidenced by systematic discrepancies between observed and modeled estimates. This study quantifies a major bias in ACI determinations, stemming from conventional surface or space measurements that fail to capture aerosol at the cloud level unless the cloud is coupled with land surface. We introduce an advanced approach to determine radiative forcing of ACI by accounting for cloud-surface coupling. By integrating field observations, satellite data, and model simulations, this approach reveals a drastic alteration in aerosol vertical transport and ACI effects caused by cloud coupling. In coupled regimes, aerosols enhance cloud droplet number concentration across the boundary layer more homogeneously than in decoupled conditions, under which aerosols from the free atmosphere predominantly affect cloud properties, leading to marked cooling effects. Our findings spotlight cloud-surface coupling as a key factor for ACI quantification, hinting at potential underassessments in traditional estimates.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Reparameterizing Litter Decomposition Using a Simplified Monte Carlo Method Improves Litter Decay Simulated by a Microbial Model and Alters Bioenergy Soil Carbon Estimates

Litter decomposition determines soil organic matter (SOM) formation and plant-available nutrient cycles. Therefore, accurate model representation of litter decomposition is critical to improving soil carbon (C) projections of bioenergy feedstocks. Soil C models that simulate microbial physiology (i.e., microbial models) are new to bioenergy agriculture, and their parameterization is often based on small datasets or manual calibration to reach benchmarks. Here, we reparameterized litter decomposition in a microbial soil C model (CORPSE - Carbon, Organisms, Rhizosphere, and Protection in the Soil Environment) using the continental-scale Long-term Inter-site Decomposition Experiment Team (LIDET) dataset which documents decomposition across a range of litter qualities over a decade. Here, we conducted a simplified Monte Carlo simulation that constrained parameter values to reduce computational costs. The LIDET-derived parameters improved modeled C and nitrogen (N) remaining, decomposition rates, and litter mean residence times as compared to Baseline parameters. We applied the LIDET litter decomposition parameters to a microbial bioenergy model (Fixation and Uptake of Nitrogen – Bioenergy Carbon, Rhizosphere, Organisms, and Protection) to examine soil C estimates generated by Baseline and LIDET parameters. LIDET parameters increased estimated soil C in bioenergy feedstocks, with even greater increases under elevated plant inputs (i.e., by increasing residue, N fertilization). This was due to the integrated effects of plant litter quantity, quality, and agricultural practices (tillage, fertilization). Collectively, we developed a simple framework for using large-scale datasets to inform the parameterization of microbial models that impacts projections of soil C for bioenergy feedstocks.

54 ENVIRONMENTAL SCIENCES↗