SEARCH · Engineering Papers

Results for “Partially explicit”

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.

SourceResearch connectionFrom yearThrough year

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

At least 73 records · Page 4

All sources NASA NTRS DOE OSTI

DOE OSTI2026

Gauged soft recursion: on-shell construction of Goldstone-gauge amplitudes

We present a new on-shell recursion relation for scattering amplitudes involving Nambu-Goldstone bosons with a gauged unbroken symmetry. A central challenge is that gauge interactions break Adler’s zero condition for charged scalars, invalidating the standard soft recursion. To overcome this, we introduce a “gauged soft recursion” that leverages the soft theorems of the gauge bosons themselves, combined with a novel decomposition of amplitudes into gauge-invariant components where Adler’s zero is partially restored. The formalism, which also incorporates internal gauge bosons via angular momentum constraints, enables the systematic construction of tree-level amplitudes with arbitrary numbers of Goldstone bosons and gauge bosons in both Abelian and non-Abelian theories, as we demonstrate with explicit examples.

Chiral Lagrangian↗

DOE OSTI2021

Multilevel Graph Partitioning for Three-Dimensional Discrete Fracture Network Flow Simulations

We present a topology-based method for mesh-partitioning in three-dimensional discrete fracture network (DFN) simulations that takes advantage of the intrinsic multi-level nature of a DFN. DFN models are used to simulate flow and transport through low-permeability fractured media in the subsurface by explicitly representing fractures as discrete entities. The governing equations for flow and transport are numerically integrated on computational meshes generated on the interconnected fracture networks. Modern high-fidelity DFN simulations require high-performance computing on multiple processors where performance and scalability depends partially on obtaining a high-quality partition of the mesh to balance work-loads and minimize communication across all processors. The discrete structure of a DFN naturally lends itself to various graph representations, which can be thought of as coarse-scale representations of the computational mesh. Using this concept, we develop two applications of the multilevel graph partitioning algorithm to partition the mesh of a DFN. In the first, we project a partition of the graph based on the DFN topology onto the mesh of the DFN and in the second, this DFN-based projection is used as the initial condition for further partitioning refinement of the mesh. We compare the performance of these methods with standard multi-level graph partitioning using graph-based metrics (cut, imbalance, partitioning time), computational-based metrics (FLOPS, iterations, solver time), and total run time. The DFN-based and the mesh-based partitioning methods are comparable in terms of the graph-based metrics, but the time required to obtain the partition is several orders of magnitude faster using the DFN-based partitions. The computation-based metrics show comparable performance between both methods so, in combination, the DFN-based partitions are several orders of magnitude faster than the mesh-based partition. Furthermore, the method which uses the DFN-partition solution as the initial condition of the mesh partition provided cut and imbalance values that were close to the mesh-based partition but in a fraction of the time. In turn, this hybrid method outperformed both of the other methods in terms of the total run time.

58 GEOSCIENCES↗

DOE OSTI2023

A general non-Fourier Stefan problem formulation that accounts for memory effects

The Stefan problem is the classical model of a melting phase change. In heterogeneous systems, such phase changes can exhibit non-Fourier (anomalous) behaviors, where the advance of the melt interface does not follow the expected time scaling. These situations can be modeled by replacing the derivatives, in the governing partial differential equations, with fractional order derivatives. In particular, replacing the time derivatives leads to non-Fourier models that account for memory effects in the system. In this work, by using appropriate time convolution integrals, a general thermodynamic balance statement for melting phase problems, explicitly accounting for memory effects, is developed. From this balance, a general model formulation applicable to problems involving melting over a temperature range (i.e., a mushy region) is derived. A key component in this model is the representation of memory effects through the use of fractional derivative based constitutive models of the enthalpy and heat flux. Further, on shrinking the mushy region to a single isotherm, a general sharp interface melting model is obtained. Here, in contrast to the classic Stefan problem, the fractional derivatives induce a natural regularization, such that the constitutive models for enthalpy and heat flux are continuous at the melt interface; a result confirmed through numerical simulation. To further support the theoretical findings, a physical example of a non-Fourier Stefan problem is presented. Overall the development and results in this paper underscore the importance of explicitly relating the development of fractional calculus models to the appropriate thermodynamic balance statements.

42 ENGINEERING↗

DOE OSTI2022

Phase Field Dislocation Dynamics (PFDD) version 2.x

This disclosure is for version 2.x of a mesoscale model called Phase Field Dislocation Dynamics (PFDD). PFDD is used for investigating deformation in nanoscale (grain sizes of ~300 nm and less) materials, such as metals and alloys. This approach models the motion and interaction of individual defects, namely dislocations, in the material using scalar-valued phase field variables, also called order parameters. The system is evolved through energy minimization thus the model calculates the total energy density in terms of the phase field variables. The energy minimization is completed using the Ginzburg-Landau equation, and is implemented with explicit time integration. The total system energy can be comprised of several terms, including the strain energy (which describes dislocation-dislocation interactions), the energy due to an applied stress (dislocation interactions with the applied stress), and a core/lattice (perfect dislocations) or generalized stacking fault (partial dislocations) energy (described the dislocation core structure). The latter term in particular may vary based on the crystal structure being modeled and is typically informed using lower length scale (e.g., atomistic) approaches, although no such (atomistic) calculations are completed within the PFDD algorithm. This basic formulation was previously reviewed by Los Alamos National Laboratory and released under license number C17113. This previously reviewed version we will henceforth refer to as PFDD v1.0. PFDD v1.0 consisted of 2 codes (one parallel and one serial) plus input files, all written in the C language. This new disclosure is addressing the next versions of the PFDD, versions 2.x. There have been several enhancements of PFDD v1.0, which are described here and included in the attached code, which we will refer to as PFDD v2.0. There are also several new features described here that are either planned or already in process and are expected to be subsequent releases, i.e., v2.1, v2.2, ...v2.x.

Hunter, Abigail↗

DOE OSTI2024

Low and moderate x gluon contribution to exclusive Compton scattering processes

We revisit the high energy semi-classical description of the exclusive processes DVCS, TCS, and Double DVCS by explicitly keeping track of the Feynman x dependence in both the hard and the hadronic matrix elements. This is achieved by a modification of the standard shock wave approximation to derive the effective Feynman rules, which leads to a generic expression on which we then perform a partial twist expansion to get rid of quantities suppressed by the proper physical scales. We obtain a compact factorized master formula that can be used to investigate the Bjorken limit at leading twist. In particular, we recover the full one-loop result in the collinear limit for pure gluon exchange with the target. Finally, we discuss the subtleties in taking the simultaneous collinear and small x limit.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2021

Robust and Spectrally Selective Aerogels for Solar Receivers

The development of next generation concentrated solar thermal (CST) plants requires solar receivers to reach high temperatures (> 600°C) while minimizing costs. Efficiently generating such high temperatures using a line-focusing collector and an air-stable receiver has the potential to change the current CSP paradigm because it would: (1) result in large improvements in system-level efficiency of line-focusing systems by enabling the use of high-efficiency supercritical CO 2 power cycles; (2) decrease the LCOE of line-focusing systems to meet SETO 2030 goals; (3) maintain performance over the lifetime of the CSP plant by avoiding issues related to breakdown of vacuum. The use of aerogels as transparent insulating materials (TIMs) in solar thermal receivers has the potential to enable such improvements. Specifically, aerogels can be placed in front of a black (or partially selective) high-temperature absorber to allow sunlight to transmit to the absorber but block heat from escaping. These materials can improve the efficiency and simplify the complexity of thermal transport systems by enabling operation at moderate or no vacuum levels. These potential benefits have been explicitly identified by SETO as ‘high impact’. The use of mesoporous silica as a TIM has already led to large improvements in solar collection efficiency. Nevertheless, these materials inherently lose structural integrity and the ability to suppress radiation losses at high temperatures. Here, we address these challenges by leveraging 1-cycle atomic layer deposition of aluminum oxide onto a silica aerogel to form a thermally stable interfacial layer that is transparent to sunlight but broadly absorbs in the mid-infrared. This interfacial absorption results in a two-fold reduction in measured heat losses from an absorber at 700°C and solar thermal efficiency of 81% at 700°C under 60 Suns, based on measurements of heat loss and solar transmittance. Extended-duration heat treatment reveals that the multicomponent aerogel is stable at high temperatures relative to a silica aerogel. Furthermore, we experimentally demonstrate the concept of using infrared plasmon resonances to selectively enhance the thermal absorption coefficient of TIMs in order to strongly suppress thermal radiative losses at high temperatures. We term this mechanism plasmon-enhanced greenhouse selectivity (PEGS). Unlike silica aerogels, where much of the IR absorption is lost at high temperatures, local surface plasmon resonances (LSPRs) in doped oxide nanoparticles overlap a significant portion of the blackbody spectrum and are maintained at high temperatures. Overall, this work paves the way for next-generation solar thermal energy using thermally robust transparent insulating materials.

14 SOLAR ENERGY↗

DOE OSTI2022

On the energy landscape of symmetric quantum signal processing

Symmetric quantum signal processing provides a parameterized representation of a real polynomial, which can be translated into an efficient quantum circuit for performing a wide range of computational tasks on quantum computers. For a given polynomial f , the parameters (called phase factors) can be obtained by solving an optimization problem. However, the cost function is non-convex, and has a very complex energy landscape with numerous global and local minima. It is therefore surprising that the solution can be robustly obtained in practice, starting from a fixed initial guess Φ 0 that contains no information of the input polynomial. To investigate this phenomenon, we first explicitly characterize all the global minima of the cost function. We then prove that one particular global minimum (called the maximal solution) belongs to a neighborhood of Φ 0 , on which the cost function is strongly convex under the condition ‖ f ‖ ∞ = O ( d − 1 ) with d = d e g ( f ) . Our result provides a partial explanation of the aforementioned success of optimization algorithms.

Wang, Jiasu↗

DOE OSTI2023

QuadConv: Quadrature-based convolutions with applications to non-uniform PDE data compression

We present a new convolution layer for deep learning architectures which we call QuadConv — an approximation to continuous convolution via quadrature. Our operator is developed explicitly for use on non-uniform, mesh-based data, and accomplishes this by learning a continuous kernel that can be sampled at arbitrary locations. Moreover, the construction of our operator admits an efficient implementation which we detail and construct. As an experimental validation of our operator, we consider the task of compressing partial differential equation (PDE) simulation data from fixed meshes. Here, we show that QuadConv can match the performance of standard discrete convolutions on uniform grid data by comparing a QuadConv autoencoder (QCAE) to a standard convolutional autoencoder (CAE). Further, we show that the QCAE can maintain this accuracy even on non-uniform data. In both cases, QuadConv also outperforms alternative unstructured convolution methods such as graph convolution.

Compression↗

DOE OSTI2023

Mechanistic Insights into Aldehyde Production from Electrochemical CO 2 Reduction on CuAg Alloy via Operando X-ray Measurements

CO 2 electrolysis converts the greenhouse gas CO 2 into valuable fuels and chemicals, such as carbon monoxide, ethylene, ethanol, etc. Currently, Cu is the only known monometallic catalyst capable of producing multicarbon products from electrochemical CO 2 reduction reaction (eCO2RR), while the poor selectivity limits its further use. It has been found that introducing Ag atoms into the Cu lattice can modulate product preference. However, the synergistic effects between Cu and Ag, and thus, the catalytic performance, are strongly influenced by catalyst morphology, electrolyzer configuration, reaction conditions, etc. Operando measurements can provide explicit information on the catalyst dynamic variation during the reaction, but their operation and analysis are challenging. Herein, we prepared CuAg multiphase alloy catalysts by magnetron sputtering, which allowed for investigating the intrinsic interaction between Cu and Ag. eCO2RR performance exhibited an improved selectivity toward carbonyls at the expense of hydrogen and hydrocarbons. The partially alloyed Cu and Ag phases were confirmed by operando X-ray diffraction. By means of combining operando X-ray measurements and density functional theory (DFT) calculations, the preferred carbonyl production is attributed to the reduced electron density and compressive strain of Cu due to Ag incorporation, which leads to a deeper d-band center and therefore weakened intermediate adsorption and oxophilicity. In conclusion, this work provides evidence of the intrinsic structural and electronic interaction between Cu and Ag during eCO2RR. The obtained information will facilitate the design of bi/multi-phase metallic or alloy electrocatalysts.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2025

Generalized entropy of gravitational fluctuations

The corrections to holographic entanglement entropy from bulk quantum fields in a classical gravitational background are now well understood. They lead, in particular, to unitary Page curves for evaporating black holes. However, the correct treatment of quantum fluctuations of the metric, including graviton excitations, is a longstanding problem. We provide a gauge-invariant prescription for the generalized entropy of gravitons in anti-de Sitter space in terms of areas and bulk entanglement entropy, generalizing the quantum extremal surface prescription to accommodate fluctuations in the semiclassical spacetime geometry. This task requires a careful treatment of the area operator on the graviton Hilbert space and the definition of a “quantum extremal gauge” in which the extremal surface is unperturbed. It also requires us to determine the correct vacuum modular Hamiltonian for the graviton field, which we fix by requiring that it doesn’t contain a boundary term in extremal gauge. We check our prescription with an explicit computation of the vacuum-subtracted generalized entropy of states containing a graviton in an AdS-Rindler background. Our results exactly match vacuum-subtracted von Neumann entropies for stress-tensor excited states in holographic conformal field theory with d > 2 dimensions. We also use covariant phase space techniques to give a partial proof of our prescription when the entanglement wedge for the background spacetime has a bifurcate Killing horizon. Along the way, we identify a class of perturbative graviton states that have parametrically larger generalized entropy, in the small G N expansion, than any low-energy excitations of an ordinary quantum field.

1/N expansion↗

DOE OSTI2025

The role of intermediate ΔΔ states in nucleon–nucleon scattering in the large-Nc and unitary limits, and ΔΔ and ΩΩ scattering

We explore potential explanations for why using large-Nc (Nc is the number of colors) scaling to determine the relative size of few-nucleon low-energy operators agrees with experiment even when dynamical Δ’s are not explicitly included. Given that the large-Nc analysis is predicated on the nucleons and Δ’s being degenerate, this is a curious result. We show that for purely S-wave interactions the relationships dictated by large-Nc scaling are unaffected whether the Δ is included or not. In the case of higher partial waves that do not mix with S-waves, the impact of the Δ is perturbative, which makes the agreement with naive (Δ-less) large-Nc ordering unsurprising. For higher partial waves that mix with S-waves, the nucleon and Δ would need to decouple to get agreement with naive large-Nc ordering. We find all NN, ΔN, and ΔΔ low energy coefficients for leading-order baryon–baryon scattering in Δ-full pionless effective field theory in terms of the two independent parameters dictated by the SU(2F) spin-flavor symmetry that arises in the Nc → ∞ limit. Because of recent lattice quantum chromodynamics results and experimental interest, we extend our analysis to the three-flavor case to study ΩΩ scattering. We show that in the unitary limit (where scattering lengths become infinite) one of the two SU(2F) parameters is driven to zero, resulting in enhanced symmetries, which agree with those found in spin-1/2 entanglement studies.

Richardson, Thomas R↗

DOE OSTI2022

Physics constrained learning for data-driven inverse modeling from sparse observations

Deep neural networks (DNN) have been used to model nonlinear relations between physical quantities. Those DNNs are embedded in physical systems described by partial differential equations (PDE) and trained by minimizing a loss function that measures the discrepancy between predictions and observations in some chosen norm. This loss function often includes the PDE constraints as a penalty term when only sparse observations are available. As a result, the PDE is only satisfied approximately by the solution. However, the penalty term typically slows down the convergence of the optimizer for stiff problems. We present a new approach that trains the embedded DNNs while numerically satisfying the PDE constraints. We develop an algorithm that enables differentiating both explicit and implicit numerical solvers in reverse-mode automatic differentiation. This allows the gradients of the DNNs and the PDE solvers to be computed in a unified framework. We demonstrate that our approach enjoys faster convergence and better stability in relatively stiff problems compared to the penalty method. Furthermore, our approach allows for the potential to solve and accelerate a wide range of data-driven inverse modeling, where the physical constraints are described by PDEs and need to be satisfied accurately.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2021

Dynamic Linking of Upstream Energy and Freight Demands for Bio and Fossil Energy Pathways in the Global Change Analysis Model

Comprehensive study of the environmental impacts associated with demand for an energy resource or carrier in any one sector requires a full consideration of the direct and indirect impacts on the rest of the regional and global energy system. It is important to consider the energy that is consumed in producing primary energy resources and energy consumed in transporting these energy resources internationally and domestically in order to produce a complete picture of these impacts. Biofuels are especially complex since they have feedbacks not just to the energy system but also to regional and global crop markets. Different modeling and analysis strategies have been applied to consider these impacts, which occur “upstream” of final energy consumption, with some success. Traditional life cycle models allow for rich technological detail in linking upstream impacts of energy demand with downstream final consumption; however, they cannot by themselves represent dynamic economic feedbacks across multiple sectors and regions over time. Computable general equilibrium (CGE) modeling does account for economic feedback among all sectors and regions; however, CGE modeling does not lend itself to physical representation of technological detail. Partial equilibrium (PE) modeling is a heterogeneous category describing economic models that focus on a subset of the economy, and they differ in coverage and ease in incorporating more sectors and economic links. In this study, we present a strategy for dynamically including the direct and indirect impacts of energy demand with physical and technology detail by explicitly adding these upstream energy and transportation links to GCAM, a PE model of global energy, land use, and emissions. We incorporate the following inter-sectoral linkages: energy inputs to crop production, energy inputs to fossil resource production, and freight transport requirements of energy and agricultural commodities. We assess the implications of explicitly including these links by measuring the global impacts of increased corn ethanol demand in the United States with and without these links included. Although the net global impact of the upstream links on energy and emissions are relatively modest in the scenarios we studied, the inclusion of these links illustrates interesting trade-offs in energy and transportation demand among fossil fuel and agriculture sectors. These sectoral interactions suggest that this level of modeling detail could be important in evaluating future analytical questions.

Sampedro Martinez de Estivariz, Jon↗

DOE OSTI2021

Ab Initio Study of Atomic and Electronic Structure of Promising Ba4XMn3O12 (X = Nb, Ce, Pr) Oxides for Solar Thermochemical Hydrogen Production

The two-step metal oxide water-splitting cycle is one of the most viable approach for Solar Thermochemical Hydrogen (STCH) production. Challenges exist in finding suitable oxides that can satisfy thermodynamics of the STCH redox cycle under viable range of temperatures and partial pressures. Recently, BXM family members, i.e., Ba4NbMn3O12 (BNM), Ba4CeMn3O12 (BCM), Ba4PrMn3O12 (BPM) are discovered to exhibit promising STCH performance. However, the magnetic degrees of freedom of their experimental crystal structures is not characterized, and only hypothetical models exist for their electronic structure. We performed Monte-Carlo sampling for the magnetic spin arrangement of Mn atoms to provide explicit atomic structure models for their ground state 12R polytype and determined the most stable spin configuration for each member after structure relaxation using density functional theory (DFT) calculations. We also elucidate the most stable charge configuration (Mn3+/Mn4+ arrangement) for BNM among the existing hypothesis in literature. We investigate these strongly correlated oxides using different levels of theory (DFT+U, Hybrid-DFT) to develop better understanding of their electronic structure. These calculations allowed us to develop a model based on molecular orbital interactions that help explain the observed changes in their experimental x-ray spectroscopy measurements. Thereby, help link differences in electronic structure of BXM family members to their varying water-splitting behavior. These structure models will further facilitate future in-depth defect studies to model STCH redox cycle of these oxides.

atomic and magnetic structure↗

DOE OSTI2025

Analytic reconstruction with massive particles: one-loop amplitudes for $0\to \overline{q} qt\overline{t}H$

We present an analytic reconstruction of one-loop amplitudes for the process $0\to \overline{q} qt\overline{t}H$. Our calculation is a novel use of analytic reconstruction, retaining explicit covariance in the massive spin states through the massive spinor-helicity formalism. The analytic reconstruction relies on embedding the massive five-point kinematics in a fully massless eight-point phase space while still building a minimal ansatz directly in the five-point phase space. In order to obtain compact analytic expressions it is necessary to identify suitable partial fraction decompositions and extract common numerator factors, which we achieve through careful inspection of limits in which pairs of denominators vanish. We find that the resulting amplitudes are more numerically efficient than ones computed using automatic methods but that the gains are not as significant as in the massless case, at least at present. The method opens the door to applications at two-loop order, where numerical efficiency and improvements in the reconstruction methodology are more crucial, especially with regards to the number of free parameters in the ansatz.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2021

Non-intrusive data-driven model reduction for differential–algebraic equations derived from lifting transformations

In this paper we present a non-intrusive data-driven approach for model reduction of nonlinear systems. The approach considers the particular case of nonlinear partial differential equations (PDEs) that form systems of partial differential–algebraic equations (PDAEs) when lifted to polynomial form. Such systems arise, for example, when the governing equations include Arrhenius reaction terms (e.g., in reacting flow models) and thermodynamic terms (e.g., the Helmholtz free energy terms in a phase-field solidification model). Using the known structured form of the lifted algebraic equations, the approach computes the reduced operators for the algebraic equations explicitly, using straightforward linear algebra operations on the basis matrices. The reduced operators for the differential equations are inferred from lifted snapshot data using operator inference, which solves a linear least squares regression problem. The approach is illustrated for the nonlinear model of solidification of a pure material. The lifting transformations reformulate the solidification PDEs as a system of PDAEs that have cubic structure. The operators of the lifted system for this solidification example have affine dependence on key process parameters, permitting us to learn a parametric reduced model with operator inference. Numerical experiments show the effectiveness of the resulting reduced models in capturing key aspects of the solidification dynamics.

42 ENGINEERING↗

DOE OSTI2025

A higher-order finite-element implementation of the nonlinear Fokker–Planck collision operator for charged particle collisions in a low density plasma

Collisions between particles in a low density plasma are described by the Fokker–Planck collision operator. In applications, this nonlinear integro-differential operator is often approximated by linearised or ad-hoc model operators due to computational cost and complexity. In this work, we present an implementation of the nonlinear Fokker–Planck collision operator written in terms of Rosenbluth potentials in the Rosenbluth–MacDonald–Judd (RMJ) form. The Rosenbluth potentials may be obtained either by direct integration or by solving partial differential equations (PDEs) similar to Poisson's equation: we optimise for performance and scalability by using sparse matrices to solve the relevant PDEs. We represent the distribution function using a tensor-product continuous-Galerkin finite-element representation and we derive and describe the implementation of the weak form of the collision operator. We present tests demonstrating a successful implementation using an explicit time integrator and we comment on the speed and accuracy of the operator. Finally, we speculate on the potential for applications in the current and next generation of kinetic plasma models.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2026

Forest aboveground biomass estimation through integration of sentinel-2 and PALSAR-2 time series: assessing models trained on GEDI and field inventory benchmarks

Accurate and spatially explicit forest Aboveground Biomass (AGB) mapping through remote sensing is critical for quantifying terrestrial carbon stocks and informing effective forest management strategies. However, AGB estimation in dense forests with complex terrain remains challenging due to satellite sensor signal saturation problem (saturation issue occurs in high biomass forests), structural complexity, and limited ground truth for calibration. This study presents a novel framework that integrates multi-temporal Sentinel-2 optical imagery, ALOS PALSAR-2 Synthetic Aperture Radar (SAR) data, and topographic variables with explainable Machine Learning to map AGB across mountainous forests within subtropical and temperate oceanic climate zones of Mexico. We evaluate the effects of temporal granularity and sensor synergy by comparing multiple temporal inputs and sensor configurations (Sentinel-2, PALSAR-2, and their fusion), and assess model performance using two reference datasets: NASA GEDI LiDAR-derived biomass and Mexico’s National Forest and Soil Inventory (INFyS). Our results showed that models trained on INFyS consistently outperformed those trained on GEDI, highlighting limitations in GEDI’s reliability in biomass estimates within this study region. Furthermore, the integration of Sentinel-2 and PALSAR-2 provided improved predictions compared to single-sensor models, particularly when combined with temporally explicit yearly statistics. The best-performing model, which was trained on INFyS data, and considered both Sentinel-2 and PALSAR-2 yearly statistics, as well as topographic variables, achieved an R2 of 0.64, RMSE of 51.10 Mg/ha, and relative RMSE (rRMSE) of 58.69%. Explainable ML analysis identified Sentinel-2 spectral indices and topographic features as key predictors, while PALSAR-2 metrics provided complementary information, partially mitigating saturation effects in high-biomass areas. Specifically, integrating both sensors substantially improved AGB estimation in high biomass forest (≥200 Mg/ha), yielding 98% gains over optical-only model, with resulting estimates exceeding GEDI L4B by 29% and ESA-CCI-BIOMASS by 174%. Terrain-stratified analysis indicated close agreement with GEDI in low-slope areas, with increasing divergence as slope steepness increased, while estimates remained consistently higher than ESA-CCI-BIOMASS across all slope classes. The proposed approach advances multi-sensor fusion and temporal feature engineering for AGB mapping using open-access satellite datasets, providing a scalable and reproducible framework for annual biomass monitoring in topographically complex mountainous forests. The resulting 25 m resolution biomass product has the potential to provide spatially detailed information for forest monitoring and may support applications in carbon accounting and forest management.

54 ENVIRONMENTAL SCIENCES↗