Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Observing System Simulation”

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 541 records · Page 30
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Power spectrum and form factor in random diagonal matrices and integrable billiards

Triggered by a controversy surrounding a universal behavior of the power spectrum in quantum systems exhibiting regular classical dynamics, we focus on a model of random diagonal matrices (RDM), often associated with the Poisson spectral universality class, and examine how the power spectrum and the form factor get affected by two-sided truncations of RDM spectra. Having developed a nonperturbative description of both statistics, we perform their detailed asymptotic analysis to demonstrate explicitly how a traditional assumption (lying at the heart of the controversy) – that the power spectrum is merely determined by the spectral form factor – breaks down for truncated spectra. This observation has important consequences as we further argue that bounded quantum systems with integrable classical dynamics are described by heavily truncated rather than complete RDM spectra. High-precision numerical simulations of semicircular and irrational rectangular billiards lend independent support to these conclusions.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Origin of Mechanical Enhancement in Polymer Nanoparticle (NP) Composites with Ultrahigh NP Loading

Polymer nanoparticle composites (PNCs) with ultrahigh loading of nanoparticles (NPs) (>50%) have been shown to exhibit markedly improved strength, stiffness, and toughness simultaneously compared to the neat systems of their components. Recent experimental studies on the effect of polymer fill fraction in these highly loaded PNCs reveal that even at low polymer fill fractions, hardness and modulus increase significantly. In this work, we aim to understand the origin of these performance enhancements by examining the dynamics of both polymer and NPs under tensile deformation. We perform molecular dynamics simulations of coarse-grained, glassy polymer in random-close-packed NP packings with a varying polymer fill fraction. We characterize the mechanical properties of the PNC systems, compare the NP rearrangement behavior, and study the polymer segmental and chain-level dynamics during deformation below the polymer glass transition. Finally, our simulation results confirm the experimentally observed increase in modulus at low polymer fill fractions, and we provide evidence that the source of mechanical enhancement is the polymer bridging effect.

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

Development and Experimental Validation of a Heat Transfer Model for Spilled Molten Salt Pools

A spill of radionuclide-bearing molten salt is one of the major postulated events that needs to be analyzed for liquid fluorine salt-cooled high-temperature reactor (FHR) or molten salt reactor licensing purposes. In this postulated event, radioactive source term materials (RSTMs) in the molten salt are discharged from the reactor vessel to the reactor building. The release of RSTMs from the spilled salt pool to the gas space in the reactor building is expected to be controlled by the cooling behavior of the spilled salt, including the growth and shrinkage of the solid crust on the surface of the spilled salt pool. This paper presents a simulation model for spilled salt pool heat transfer and validation efforts. The validation data come from two molten salt spill tests that were performed recently: the PELE2 test by the Rapid Experimental Laboratory of Kairos Power LLC (KP) and the Argonne salt cooling test conducted by Argonne National Laboratory. The former was a large-scale test involving kilograms of molten spilled FLiNaK salt, and the latter was a relatively smaller-scale test targeting various processes associated with a salt spill event. Both tests generated valuable data sets that can be used to assess salt cooling and validate evaluation models. This paper provides a new one-dimensional model that can simulate the cooling process of a spilled salt pool as well as the thermal responses of heat structures, such as the stainless steel liner and the concrete below the salt. The model has been implemented as part of KP-SAM code, which is a branch of the systems code SAM specific to KP FHR. In conclusion, the simulation results of the model are compared with the data of the PELE2 and Argonne tests, and reasonable agreements are observed between the model and test data.

heat transfer model↗
DOE OSTI2022

Scaling Superconducting Quantum Computers with Chiplet Architectures

Fixed-frequency transmon quantum computers (QCs) have advanced in coherence times, addressability, and gate fidelities. Unfortunately, these devices are restricted by the number of on-chip qubits, capping processing power and slowing progress toward fault-tolerance. Although emerging transmon devices feature over 100 qubits, building QCs large enough for meaningful demonstrations of quantum advantage requires overcoming many design challenges. For example, today’s transmon qubits suffer from significant variation due to limited precision in fabrication. As a result, barring significant improvements in current fabrication techniques, scaling QCs by building ever larger individual chips with more qubits is hampered by device variation. Severe device variation that degrades QC performance is referred to as a defect. Here, we focus on a specific defect known as a frequency collision. When transmon frequencies collide, their difference falls within a range that limits two-qubit gate fidelity. Frequency collisions occur with greater probability on larger QCs, causing collision-free yields to decline as the number of on-chip qubits increases. As a solution, we propose exploiting the higher yields associated with smaller QCs by integrating quantum chiplets within quantum multi-chip modules (MCMs). Yield, gate performance, and application-based analysis show the feasibility of QC scaling through modularity. Our results demonstrate that chiplet architectures, relative to monolithic designs, benefit from average yield improvements ranging from 9.6 – 92.6 × for ≲5 qubit machines. In addition, our simulations explore the design space of chiplet systems and discover configurations that demonstrate average two-qubit gate infidelity reductions that are at best 0.815 × their monolithic counterpart. Lastly, we observe that carefully-selected modular systems achieve fidelity improvements on a range of benchmark circuits.

quantum architecture↗
DOE OSTI2025

Challenges in Continuous In-Field Critical Current Testing of High-Temperature Superconducting Tapes: Thermal and Mechanical Perspectives

High-temperature superconductors (HTS) are essential for ultra-high-field applications requiring exceptional current-carrying capacity under extreme conditions. However, systematic characterization of critical current in long-length conductors remains challenging due to complex thermal, electromag netic, and mechanical interactions during continuous testing. This study reports the development of a continuous in-field magnetization testing system for position-dependent critical current measurement in HTS tapes at 20 K under 7.5 T fields applied normal to the tape plane, enabling identification of performance-limiting regions that could compromise magnet stability. Here, the system addresses two fundamental challenges inherent to cryogenic reel to-reel testing. First, thermal management requires continuous cooling of a moving conductor to 20 K, achieved through liquid nitrogen precooling combined with a 100 W@20 K Gifford McMahon cryocooler. Second, screening currents in high fields generate Lorentz forces that induce twisting, bowing, and potential delamination. To mitigate these risks, we propose mechanical reinforcement and active current density suppression strategies. Numerical simulations using the stream function formulation reveal four primary failure modes: frictional heating at guide interfaces, unstable equilibria causing deformation, transverse current-induced stresses at guide transitions, and unsupported forces in vertical spans. Our mitigation strategies include PTFE coated guides to minimize friction, spring-loaded stabilization mechanisms to maintain tape alignment, controlled pre-heating using the liquid nitrogen thermal jacket to suppress critical current at stress points, and optimized guide positioning to minimize force accumulation. The experimental system is nearing completion, with testing planned to commence within two months. Preliminary validation at 65 K under 0.5 T demonstrates strong correlation between simulation-predicted mechanical instabilities and observed critical current variations during conductor tran sitions through the measurement region. These findings establish a robust foundation for quality assurance protocols essential to next-generation superconducting magnet applications.

Chen, Siwei [Princeton Plasma Physics Laboratory (↗
DOE OSTI2023

Cyber Physical Subsystems under Adverse Conditions

The United States power grid is built upon complex cyber-physical systems (CPS), including many distributed energy resources (DERs) and microgrid structures. Due to cross-disciplinary domains that govern the interplay between cyber and physical technologies that sustain DERs, it is difficult to confidently predict how these systems will react under adverse conditions such as faults and cyberattacks. This project utilized power hardware in the loop (PHiL) to gain system behavior knowledge and robust understanding through experimentation and observational science. It also produced a large volume of datasets to feed into other efforts to advance high-fidelity model simulations and to support the design of novel control approaches.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2023

Surface impacts and associated mechanisms of a moisture intrusion into the Arctic observed in mid-April 2020 during MOSAiC

Distinct events of warm and moist air intrusions (WAIs) from mid-latitudes have pronounced impacts on the Arctic climate system. We present a detailed analysis of a record-breaking WAI observed during the MOSAiC expedition in mid-April 2020. By combining Eulerian with Lagrangian frameworks and using simulations across different scales, we investigate aspects of air mass transformations via cloud processes and quantify related surface impacts. The WAI is characterized by two distinct pathways, Siberian and Atlantic. A moist static energy transport across the Arctic Circle above the climatological 90th percentile is found. Observations at research vessel Polarstern show a transition from radiatively clear to cloudy state with significant precipitation and a positive surface energy balance (SEB), i.e., surface warming. WAI air parcels reach Polarstern first near the tropopause, and only 1–2 days later at lower altitudes. In the 5 days prior to the event, latent heat release during cloud formation triggers maximum diabatic heating rates in excess of 20 K d -1 . For some poleward drifting air parcels, this facilitates strong ascent by up to 9 km. Based on model experiments, we explore the role of two key cloud-determining factors. First, we test the role moisture availability by reducing lateral moisture inflow during the WAI by 30%. This does not significantly affect the liquid water path, and therefore the SEB, in the central Arctic. The cause are counteracting mechanisms of cloud formation and precipitation along the trajectory. Second, we test the impact of increasing Cloud Condensation Nuclei concentrations from 10 to 1,000 cm -3 (pristine Arctic to highly polluted), which enhances cloud water content. Resulting stronger longwave cooling at cloud top makes entrainment more efficient and deepens the atmospheric boundary layer. Finally, we show the strongly positive effect of the WAI on the SEB. This is mainly driven by turbulent heat fluxes over the ocean, but radiation over sea ice. The WAI also contributes a large fraction to precipitation in the Arctic, reaching 30% of total precipitation in a 9-day period at the MOSAiC site. However, measured precipitation varies substantially between different platforms. Therefore, estimates of total precipitation are subject to considerable observational uncertainty.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

SDSS-IV MaNGA: Cannibalism Caught in the Act—On the Frequency of Occurrence of Multiple Cores in Brightest Cluster Galaxies

Although it is generally accepted that massive galaxies form in a two-phased fashion, beginning with a rapid mass buildup through intense starburst activities followed by primarily dry mergers that mainly deposit stellar mass at outskirts, the late time stellar mass growth of brightest cluster galaxies (BCGs), the most massive galaxies in the universe, is still not well understood. Several independent measurements have indicated a slower mass growth rate than predictions from theoretical models. We attempt to resolve the discrepancy by measuring the frequency of BCGs with multiple cores, which serve as a proxy of the merger rates in the central region and facilitate a more direct comparison with theoretical predictions. Using 79 BCGs at z = 0.06–0.15 with integral field spectroscopic data from the Mapping Nearby Galaxies at Apache Point Observatory (MaNGA) project, we obtain a multiple-core fraction of 0.11 ± 0.04 at z ≈ 0.1 within an 18 kpc radius from the center, which is comparable to the value of 0.08 ± 0.04 derived from mock observations of 218 simulated BCGs from the cosmological hydrodynamical simulation IllustrisTNG. We find that most cores that appear close to the BCGs from imaging data turn out to be physically associated systems. Anchoring on the similarity in the multiple-core frequency between the MaNGA and IllustrisTNG, we discuss the mass growth rate of BCGs over the past 4.5 Gyr.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Model Form Error Correction for a Black-Box Thermal Battery Heat Transfer Simulation

Thermal batteries are crucial for supplying power to high-consequence engineering applications such as rockets. Computational simulations have been developed to predict thermal battery behavior, but these simulations often suffer from modeling errors, including model form uncertainty. Addressing this uncertainty can be achieved by quantifying either the model discrepancy in the output or the model form error (MFE) in the governing equation. MFE is particularly valuable as it can be better extrapolated beyond observed outputs, which is essential for predictions involving changes in external system loading, system configuration and geometry, or output quantities. This paper employs a state estimation approach to estimate MFE using experimental data and then utilizes machine learning (ML) to model its relationship with state variables. A nonintrusive technique is used to estimate MFE in a black-box thermal battery heat transfer simulation. The trained machine learning model for MFE is then applied to correct simulation predictions under extrapolated initial conditions and battery configurations. In conclusion, the methodology's performance is evaluated using additional experimental data, demonstrating its effectiveness in improving prediction accuracy.

Batteries↗
DOE OSTI2020

Superfluid condensate fraction and pairing wave function of the unitary Fermi gas

The unitary Fermi gas is a many-body system of two-component fermions with zero-range interactions tuned to infinite scattering length. Despite much activity and interest in unitary Fermi gas and its universal properties, there have been great difficulties in performing accurate calculations of the superfluid condensate fraction and pairing wave function. In this paper, we present auxiliary-field lattice Monte Carlo simulations using a lattice interaction which accelerates the approach to the continuum limit, thereby allowing for robust calculations of these difficult observables. As a benchmark test, we compute the ground-state energy of 33 spin-up and 33 spin-down particles. As a fraction of the free Fermi gas energy $E_{\text{FG}}$, we find $E_0/E_{\text{FG}}$ = 0.369(2), 0.372(2), using two different definitions of the finite-system energy ratio, in agreement with the latest theoretical and experimental results. We then determine the condensate fraction by measuring off-diagonal long-range order in the two-body density matrix. We find that the fraction of condensed pairs is α = 0.43(2). Further, we also extract the pairing wave function and find the pair correlation length to be $ζ_pk_F$ = 1.8(3)ℏ, where $k_F$ is the Fermi momentum. Provided that the simulations can be performed without severe sign oscillations, the methods we present here can be applied to superfluid neutron matter as well as more exotic $\textit{P}$-wave and $\textit{D}$-wave superfluids.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Kinetics of Carbon Condensation in Detonation of High Explosives: First-Order Phase Transition Theory Perspective

We report the kinetics of carbon condensation, or carbon clustering, in detonation of carbon-rich high explosives is modeled by solving a system of rate equations for concentrations of carbon particles. Unlike previous efforts, the rate equations account not only for the aggregation of particles but also for their fragmentation in a thermodynamically consistent manner. Numerical simulations are performed, yielding the distribution of particle concentrations as a function of time. In addition to that, analytical expressions are obtained for all the distinct steps and regimes of the condensation kinetics, which facilitates the analysis of the numerical results and allows one to study the sensitivity of the kinetic behavior to the variation of system parameters. The latter is important because the numerical values of many parameters are not reliably known at present. The theory of the kinetics of first-order phase transitions is found adequate to describe the general kinetic trends of carbon condensation, as described by the rate equations. Such physical phenomena and processes as the coagulation, nucleation, growth, and Ostwald ripening are observed, and their dependence on various system parameters is studied and reported. It is believed that the present work will become useful when analyzing the present and future results for the kinetics of carbon condensation, obtained from experiments or atomistic simulations.

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

Role of Dynamic Polarization Interactions in the Electrical Double Layer at Calcite (104) Interfaces with Aqueous Solutions

Reactions at mineral interfaces with aqueous solutions control many geochemical and biogeochemical processes in the Earth’s critical zone. At the molecular level, insights into important properties such as the structure of the electrical double layer (EDL) at specific mineral interfaces continue to improve due to the increasing fidelity of laboratory instrumentation and computational approaches. However, molecular simulation approaches suffer from limited reach into relevant scales of time, length, and system complexity. To span this gap, a novel hybrid approach that couples first principle plane-wave density functional theory (DFT) with classical DFT (cDFT) is demonstrated and applied to calcite (104) interfaces with various electrolytes. In this approach, a region of interest described using DFT interacts with the surrounding medium described using cDFT to arrive at a self-consistent ground state. Benchmarking against experimental observations and entirely first principle DFT simulations demonstrates that this hybrid model efficiently encompasses the key short-range and collective interactions in the EDL. Simulations of calcite (104)/solution interfaces reveal the key static and dynamic polarization interactions that give rise to structuring of ions and water. Ion hydration interactions have the strongest effect on the depth of the first minimum in the density distribution of counterions at the surface, and the position and width of the first density peak is largely determined by the strength of ion-correlation forces. Finer details of ion distributions are controlled by mutual polarization of the calcite surface and interfacial electrolyte. Finally, this new ability to efficiently and rigorously predict EDL structure at mineral surfaces in contact with complex solutions paves the way to accurately modeling sorption, nucleation, dissolution, and growth in realistic systems.

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

Hybrid Basis and Multi-Center Grid Method for Strong-Field Processes

We present a time-dependent framework that combines a hybrid basis, consisting of Gaussian-type orbitals (GTOs) and finite-element discrete-variable representation (FEDVR) functions, with a multicenter grid to simulate strong-field and attosecond dynamics in atoms and molecules. The method incorporates the construction of the orthonormal hybrid basis, the evaluation of electronic integrals, a unitary time-propagation scheme, and the extraction of optical and photoelectron observables. Its accuracy and robustness are benchmarked on one-electron systems such as atomic hydrogen and the dihydrogen cation (H$^+_2$) through comparisons with essentially-exact reference results for bound-state energies, high-harmonic generation spectra, photoionization cross sections, and photoelectron momentum distributions. This work establishes the groundwork for its integration with quantum-chemistry methods, which is already operational but will be detailed in future work, thereby enabling ab initio simulations of correlated polyatomic systems in intense ultrafast laser fields.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2022

Internal variability and forcing influence model–satellite differences in the rate of tropical tropospheric warming

Climate-model simulations exhibit approximately two times more tropical tropospheric warming than satellite observations since 1979. The causes of this difference are not fully understood and are poorly quantified. Here, we apply machine learning to relate the patterns of surface-temperature change to the forced and unforced components of tropical tropospheric warming. This approach allows us to disentangle the forced and unforced change in the model-simulated temperature of the midtroposphere (TMT). In applying the climate-model-trained machine-learning framework to observations, we estimate that external forcing has produced a tropical TMT trend of 0.25 ± 0.08 K⋅decade −1 between 1979 and 2014, but internal variability has offset this warming by 0.07 ± 0.07 K⋅decade −1 . Using the Community Earth System Model version 2 (CESM2) large ensemble, we also find that a discontinuity in the variability of prescribed biomass-burning aerosol emissions artificially enhances simulated tropical TMT change by 0.04 K⋅decade −1 . The magnitude of this aerosol-forcing bias will vary across climate models, but since the latest generation of climate models all use the same emissions dataset, the bias may systematically enhance climate-model trends over the satellite era. Our results indicate that internal variability and forcing uncertainties largely explain differences in satellite-versus-model warming and are important considerations when evaluating climate models.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Influence of the Dielectric Constant on the Ionic Current Rectification of Bipolar Nanopores

In this paper, we investigate how the dielectric constant, ϵ, of an electrolyte solvent influences the current rectification characteristics of bipolar nanopores. It is well recognized that bipolar nanopores with two oppositely charged regions rectify current when exposed to an alternating electric potential difference. Here, we consider dilute electrolytes with NaCl only and with a mixture of NaCl and charged nanoparticles. These systems are studied using two levels of description, all-atom explicit water molecular dynamics (MD) simulations and coarse-grained implicit solvent MD simulations. The charge density and electric potential profiles and current-voltage relationship predicted by the implicit solvent simulations with ϵ = 11.3 show good agreement with the predictions from the explicit water simulations. Under nonequilibrium conditions, the predictions of the implicit solvent simulations with a dielectric constant closer to the one of bulk water are significantly different from the predictions obtained with the explicit water model. Further, these findings are closely aligned with experimental data on the dielectric constant of water when confined to nanometric spaces, which suggests that ϵ decreases significantly compared to its value in the bulk. Moreover, the largest electric current rectification is observed in systems containing nanoparticles when ϵ = 78.8. Using enhanced sampling, we have shown that this larger rectification arises from the presence of a significantly deeper minimum in the free energy of the system with a larger ϵ, and when a negative voltage bias is applied. Since implicit solvent models and mean-field continuum theories are often used to design Janus membranes based on bipolar nanopores, this work highlights the importance of properly accounting for the effects of confinement on the dielectric constant of the electrolyte solvent. The results presented here indicate that the dielectric constant in implicit solvent simulations may be used as an adjustable parameter to approximately account for the effects of nanometric confinement on aqueous electrolyte solvents.

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

Future projections of wind patterns in California with the variable-resolution CESM: a clustering analysis approach

Wind energy production is expected to be affected by shifts in wind patterns that will accompany climate change. However, many questions remain on the magnitude and character of this impact, especially on regional scales. In this study, clustering is used to group and analyze large-scale wind patterns in California using model simulations from the variable-resolution Community Earth System Model (VR-CESM). Specifically, simulations have been produced that cover historical (1980–2000), mid-century (2030–2050), and end-of-century (2080–2100) time periods. Once clustered, observed changes to wind patterns can be analyzed in terms of both the change in frequency of those clusters and changes to winds within-clusters. Statistically significant capacity factors changes have been found at all five wind plant sites. Further, decomposition of the capacity factor changes into frequency changes and within-cluster changes enables a better understanding of their drivers. A further examination of the synoptic-scale fields associated with each cluster then provides a better understanding of how changes to large-scale meteorological fields are important for driving changes in localized wind speeds.

17 WIND ENERGY↗
DOE OSTI2020

Role of the Terrestrial Biosphere in Atmospheric Chemistry and Climate

The terrestrial biosphere–atmosphere interface provides a key chemical, biological, and physical lower boundary for the atmosphere. The presence of vegetation itself modifies the physical boundary, or the biogeophysical aspects of the system, by controlling important climate drivers such as soil moisture, light environment, and temperature. The leaf surface area of the terrestrial biosphere provides additional surface area for emissions, and it can be up to 55% of the total Earth’s surface area during the boreal summer. Vegetation also influences the biogeochemical aspects of the system by emitting a broad suite of reactive trace gases such as biogenic volatile organic compound (BVOC) emissions and climate-relevant primary biological aerosol particles (PBAP). Many of these emissions are a function of meteorological and climatological conditions at the surface, including temperature, light environment, soil moisture, and winds. Once emitted, they can be processed in the troposphere through a suite of chemical reactions. BVOC can contribute to the formation of ozone and secondary organic aerosols (SOA), and PBAP can rupture to form smaller particles with climatic relevance. These emissions and subsequent aerosol products can influence atmospheric processes that affect the surface climate, such as the attenuation of radiation, the formation of greenhouse gases such as ozone that can feedback to surface air temperature, and the alteration of clouds and subsequent precipitation. These atmospheric changes can then feedback to the land surface and emissions themselves, creating positive or negative feedback loops that can dampen or amplify the emission response. For the dominant BVOC isoprene, the feedback response to temperature can be positive or negative depending on ambient temperatures that drive isoprene emissions. The feedback response to soil moisture and precipitation can be positive, negative, or uncoupled depending on the moisture content of the soil and the total atmospheric aerosol loading. For light, the isoprene response can be positive or negative depending on the role of diffuse light. Overall, these feedbacks highlight the dynamical response of the biosphere to changing atmospheric conditions across a range of time scales, from minutes for trace gases and aerosols, to months for phenological changes, to years for land cover and land use change. Here, the dynamic aspect of this system requires us to understand, simulate, and predict the complex feedbacks between the biosphere and atmosphere and understand their role in the simulation and understanding of climate and global change. From the observational perspective, these feedbacks are challenging to identify in observations, and predictive modeling tools provide a crucial link for understanding how these feedbacks will change under warming climate scenarios.

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

The Effect of Physically Based Ice Radiative Processes on Greenland Ice Sheet Albedo and Surface Mass Balance in E3SM

Abstract A significant portion of surface melt on the Greenland Ice Sheet (GrIS) is due to dark ice regions in the ablation zone, where solar absorption is influenced by the physical properties of the ice, light absorbing constituents (LACs), and the overlying crustal surface or melt ponds. Earth system models (ESMs) typically prescribe the albedo of ice surfaces as a constant value in the visible and near‐infrared spectral regions. This work advances ESM ice radiative transfer modeling by (a) incorporating a physically based radiative transfer model (SNow, ICe and Aerosol Radiation model Adding‐Doubling Version 4; SNICAR‐ADv4) into the Energy Exascale Earth System Model (E3SM), (b) determining spatially and temporally varying bare ice physical properties over the GrIS ablation zone from satellite observations to inform SNICAR‐ADv4, and (c) assessing the impacts on simulated GrIS albedo and surface mass balance associated with modeling of more realistic bare ice albedo. GrIS‐wide bare ice albedo in E3SMv2 is overestimated by ∼4% in the visible and ∼7% in the near‐infrared wavelengths compared to the Moderate Resolution Imaging Spectroradiometer. Our bare ice physical property retrieval method found that LACs, ice crustal surfaces, and melt ponds reduce visible albedo by 30% in the bare ice region of the GrIS ablation zone. The realistic bare ice albedo reduces surface mass balance by ∼145 Gt, or 0.4 mm of sea‐level equivalent between 2000 and 2021 compared to the default E3SM. This work highlights the importance of simulating bare ice albedo accurately and realistically to improve our ability to quantify changes in the GrIS surface mass and radiative energy budgets.

58 GEOSCIENCES↗