Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “structural alignment”

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 721 records · Page 40
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Effect of Solvent Composition on Non-DLVO Forces and Oriented Attachment of Zinc Oxide Nanoparticles

Oriented attachment (OA) occurs when nanoparticles in solution align their crystallographic axes prior to colliding and subsequently fuse into single crystals. Traditional colloidal theories such as DLVO provide a framework for evaluating OA but fail to capture key particle interactions due to the atomistic details of both the crystal structure and the interfacial solution structure. Using zinc oxide as a model system, we investigated the effect of the solvent on short-ranged and long-ranged particle interactions and the resulting OA mechanism. In situ TEM imaging showed that ZnO nanocrystals in toluene undergo long-range attraction comparable to 1kT at separations of 10 nm and 3kT near particle contact. These observations were rationalized by considering non-DLVO interactions, namely dipole-dipole forces and torques between the polar ZnO nanocrystals. Langevin dynamics simulations showed stronger interactions in toluene compared to methanol solvents, consistent with the experimental results. Concurrently, we performed atomic force microscopy measurements using ZnO-coated probes for the short-ranged interaction. Our data provided valuable insights into another type of non-DLVO interaction, namely the repulsive solvation force. Specifically, the solvation force was stronger in water compared to ethanol and methanol, due to the stronger hydrogen bonding and denser packing of water molecules at the interface. In conclusion, our results highlight the importance of non-DLVO forces in a general framework for understanding and predicting particle aggregation and attachment.

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

Recyclable Graphene Sheets as a Growth Template for Crystalline ZnO Nanowires

Recent advances in nanoscience have opened ways of recycling substrates for nanomaterial growth. Novel materials, such as atomically thin materials, are highly desirable for the recycling substrates. In this work, we report recycling of monolayer graphene as a growth template for synthesis of single crystalline ZnO nanowires. Selective nucleation of ZnO nanowires on graphene was elucidated by scanning electron microscopy and density functional theory calculation. Growth and subsequent separation of ZnO nanowires was repeated up to seven times on the same monolayer graphene film. Raman analyses were also performed to investigate the quality of graphene structure along the recycling processes. The chemical robustness of graphene enables the repetitive ZnO nanowire growth without noticeable degradation of the graphene quality. This work presents a route for graphene as a multifunctional growth template for diverse nanomaterials such as nanocrystals, aligned nanowires, other two-dimensional materials, and semiconductor thin films.

36 MATERIALS SCIENCE↗
DOE OSTI2026

2025 EIC-France Workshop: Physics Highlights and Perspectives

This document presents a synthesis of the theory contributions and discussions from the 2nd EIC-France Workshop, held at IJCLab (Orsay) on 1-3 December 2025. The workshop brought together members of the French hadron-physics community to review recent theoretical developments relevant to the future Electron-Ion Collider (EIC) and to coordinate national efforts in preparation for its early physics program. The report first summarizes the collider's initial running conditions and luminosity performance, as outlined in the EIC Early Science Matrix. It then provides concise overviews of the theoretical presentations on inclusive, semi-inclusive, exclusive, heavy-flavor, and small-x physics. Based on these discussions, two measurements emerged as especially well suited for early EIC operation and strongly aligned with areas of established French expertise: inclusive diffraction and inclusive quarkonium production. These channels offer clean signatures, robust theoretical interpretability, and direct sensitivity to fundamental QCD phenomena such as gluon saturation, heavy-quark dynamics, and the small-x structure of hadrons and nuclei. In addition, the workshop identified longer-term physics opportunities that will benefit from the full capabilities of the EIC after its ramp-up phase. These include accessing the three-dimensional structure of the pion through the Sullivan process and a broader program of exclusive three-body final states, both of which represent high-impact avenues for exploring hadronic structure and non-perturbative QCD. Together, the elements summarized in this report provide a coherent overview of the strategic priorities and scientific ambitions shaping the French community's contribution to the EIC physics program.

FOS: Physical sciences↗
DOE OSTI2023

Neutron Skin Measurement of 208Pb and 48Ca Using Parity Violating Electron Scattering

Parity-violating electron scattering experiments (PVES) provides a clean probe of neutron densities that is model independent and free from strong interaction uncertainties in interpretation. The PREX-2 and CREX experiments were run in 2019 and 2020 at Jefferson laboratory measured the nucleon skin thickness, the difference between the r.m.s. neutron radius Rn and the r.m.s. proton radius Rp, of 208Pb and 48Ca via parity violating electroweak asymmetry in the elastic scattering of longitudinally polarized electrons. PREX-2 experiment was performed with 950 MeV electrons scattered at a 5¿ angle with Q2 = 0.00616 ± 0.00004 (GeV/c)2, while CREX used 2182 MeV electrons at the same angle with Q2 = 0.0297 ± 0.0002 (GeV/c)2. For PREX-2 the measured asymmetry was AP V = 550 ± 16 (stat.) ± (syst.) ppb, which corresponds to Rskin = 0.278 ± 0.078 (exp.) ± 0.012 (theo.) fm. The CREX asymmetry was AP V = 2668 ± 106 (stat.) ± 40 (syst.) ppb, which corresponds to Rskin = 0.121 ± 0.026 (exp.) ±0.024 (model) fm. One of the crucial systematic uncertainty that PREX-2 and CREX were sensitive to was the non-parity violating asymmetries that resulted from the helicity correlated false asymmetries in the polarized electron beam. There was a lot of work put towards understanding and suppressing the asymmetries arising from these effects. The parity violating asymmetry measurement required a very precise determination of the electron beam polarization. To accurately determine the beam polarization, a Compton polarimeter was used during both PREX-2 and CREX. A careful alignment of the laser to the Fabry-Perot cavity, data analysis and systematic control was employed to get a precise beam polarization result for the experiments. The PREX-2 measurement has broad implications for increasing our knowledge about neutron star structure and the equation of state of nuclear matter. The combined PREX and CREX results will have implications for future energy density functional calculations and the theory of nuclear structure.

Premathilake, Sachinthani↗
DOE OSTI2021

Tightly-coupled camera/LiDAR integration for point cloud generation from GNSS/INS-assisted UAV mapping systems

Unmanned aerial vehicles (UAVs) equipped with integrated global navigation satellite systems/inertial navigation systems (GNSS/INS) together with cameras and/or LiDAR sensors are being widely used for topographic mapping in a variety of applications such as precision agriculture, coastal monitoring, and archaeological documentation. Integration of image-based and LiDAR point clouds can provide a comprehensive 3D model of the area of interest. For such integration, ensuring a good alignment between data from the different sources is critical. Although many works have been conducted on this topic, there is still a need for a rigorous integration approach that minimizes the discrepancy between camera and LiDAR data caused by inaccurate system calibration parameters and/or trajectory artifacts. This study proposes an automated tightly-coupled camera/LiDAR integration workflow for GNSS/INS-assisted UAV systems. The proposed strategy is conducted in three main steps. First, an image-based point cloud is generated using a LiDAR/GNSS/INS-assisted structure from motion (SfM) strategy. Then, feature correspondences between image-based and LiDAR point clouds are automatically identified. Finally, an integrated-bundle adjustment procedure including image points, LiDAR raw measurements, and GNSS/INS information is conducted to minimize the discrepancy between point clouds from different sensors while estimating system calibration parameters and refining the trajectory information. The proposed SfM strategy and integration framework are evaluated using five datasets. The SfM results show that using LiDAR data can facilitate feature matching and further increase the number of reconstructed 3D points. The experimental results also illustrate that the developed automated camera/LiDAR integration strategy is capable of accurately estimating system calibration parameters to achieve good alignment among camera/LiDAR data from single/multiple systems. Finally, an absolute accuracy in the range of 3–5 cm is achieved for the image/LiDAR point clouds after the integration process.

42 ENGINEERING↗
DOE OSTI2022

In Situ Studies of the Formation of MoP Catalysts and Their Structure under Reaction Conditions for Higher Alcohol Synthesis: The Role of Promoters and Mesoporous Supports

Three formulations of a molybdenum phosphide (MoP) catalyst system were characterized for the higher alcohol synthesis (HAS) reaction using in situ X-ray absorption spectroscopy (XAS) and X-ray diffraction (XRD), monitoring the chemical phase evolution during activation and under reaction conditions. Here, the in situ study herein provides important insights into the effect of the support and of K that lead to high performance in HAS for K-promoted MoP supported on carbon, as evidenced by previous studies (ethanol selectivity: 29%; conversion: 5%). During the activation process, XAS shows that the carbon-supported samples reduce and reach a highly crystalline state at a lower temperature than the SBA-15-supported sample, indicating a substantial difference in catalyst activation. After activation, the samples are introduced to relevant reaction conditions resulting in spectra fairly similar to one another. XRD results corroborate the difference in the degree of crystallinity of these samples, in alignment with the XAS, and reveal the formation of crystalline potassium pyrophosphate (K 4 P 2 O 7 ) during the activation period of the K-promoted samples. This K 4 P 2 O 7 phase remains present under reaction conditions. Taken together, these results provide insights into the roles played by the carbon support and K promotion, connecting activity to electronic and crystal structure.

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

Exploring microstructures and anisotropies of serpentinites

Serpentine minerals have received a lot of attention because of their unique crystal structures, their wide occurrence in orogenic belts and their potential role in contributing seismic anisotropy in subducting slabs. Several studies have investigated crystal preferred orientation (CPO) in high temperature antigorite serpentinites from Japan, the Alps, Spain, Cuba and Tibet, documenting significant crystal alignment. However, only a limited number of lower grade serpentines have been explored to date. Mainly because of submicroscopic microstructural heterogeneities CPO cannot be measured with conventional methods such as optical microscopy and EBSD. In this study 15 serpentinites from different tectonic settings in California, the Central Alps and Northern Spain have been investigated, mainly with high energy synchrotron X-ray diffraction, to quantify bulk crystal alignment. We find that CPO is strong on sheared surfaces of fractured blocks and secondary veins but the bulk of most serpentinite samples, except high-grade recrystallized antigorite serpentinite, show only weak crystal alignment. Correspondingly calculated seismic anisotropy based on CPO is not very significant. This is supported by very heterogeneous microstructures as documented with SEM and TEM analyses.

58 GEOSCIENCES↗
DOE OSTI2025

Improved damage tolerance of SiC-based nuclear fuel cladding with novel multi-layered SiC coating design at 1200 °C

Continuous SiC fibre reinforced SiC matrix composites (SiC f -SiC m ) with monolithic SiC outer coatings are considered as a damage-tolerant cladding design for loss of coolant accident (LOCA) conditions in light water reactors. However, monolithic SiC coatings are brittle and prone to catastrophic failure. In this study, a SiC f -SiC m cladding with a novel multi-layer SiC outer coating (11 sub-layers, ∼260 µm in total thickness) was investigated under C-ring compression at room temperature and 1200 °C in argon environment. Real-time synchrotron X-ray computed tomography (XCT) was employed to capture crack initiation and propagation processes. Compared to conventional monolithic SiC outer coatings, the multi-layer coating structure facilitated crack deflection and bifurcation enhancing its damage tolerance at both temperatures. Despite pre-existing surface cracks, claddings exhibited stable mechanical-performance at both temperatures. These initial cracks did not critically affect the failure processes as they were not aligned with the maximum stress direction. Furthermore, the microstructure, distribution of residual stresses, and local properties of individual components in the material were thoroughly characterized, and compared with open literature on conventional claddings with monolithic outer coatings. These results provide new insights into the failure mechanisms of multi-layer SiC coatings and offer guidance for the future design of accident-tolerant nuclear fuel claddings.

36 - MATERIALS SCIENCE↗
DOE OSTI2022

Perovskite superlattices with efficient carrier dynamics

Compared with their three-dimensional (3D) counterparts, low-dimensional metal halide perovskites (2D and quasi-2D; B 2 A n-1 M n X 3n+1 , such as B = R-NH 3 + , A = HC(NH 2 ) 2 + , Cs + ; M = Pb 2+ , Sn 2+ ; X = Cl - , Br - , I - ) with periodic inorganic-organic structures have shown promising stability and hysteresis-free electrical performance. However, their unique multiple-quantum-well structure limits the device efficiencies because of the grain boundaries and randomly oriented quantum wells in polycrystals. In single crystals, the carrier transport through the thickness direction is hindered by the layered insulating organic spacers. Furthermore, the strong quantum confinement from the organic spacers limits the generation and transport of free carriers. Also, lead-free metal halide perovskites have been developed but their device performance is limited by their low crystallinity and structural instability 11 . Here, in this paper, we report a low-dimensional metal halide perovskite BA 2 MA n-1 Sn n I 3n+1 (BA, butylammonium; MA, methylammonium; n = 1, 3, 5) superlattice by chemical epitaxy. The inorganic slabs are aligned vertical to the substrate and interconnected in a criss-cross 2D network parallel to the substrate, leading to efficient carrier transport in three dimensions. A lattice-mismatched substrate compresses the organic spacers, which weakens the quantum confinement. The performance of a superlattice solar cell has been certified under the quasi-steady state, showing a stable 12.36% photoelectric conversion efficiency. Moreover, an intraband exciton relaxation process may have yielded an unusually high open-circuit voltage (V OC ).

36 MATERIALS SCIENCE↗
DOE OSTI2025

LatticeAnalytics: Strut-Level Visualization and Inspection of Additively Manufactured Lattice Structures

Additive manufacturing (AM) is revolutionizing the production of custom components with complex internal geometries, essential for high-performance applications in diverse fields such as medicine and defense. These AM parts optimize strength while minimizing weight by utilizing internal lattice structures consisting of large quantities of small interconnected struts. However, the complexity of these structures, combined with the challenges of using X-ray Computed Tomography (XCT) data, makes validation of part reliability difficult. This ultimately inhibits the development of novel parts for our collaborating material scientists. Here, we introduce LatticeAnalytics, a novel framework specifically designed for visual inspection of defects in these lattice structures. Our framework offers an end-to-end solution that includes the data management of XCT scans, enables remote access for geographically dispersed teams through a web-based dashboard, and incorporates novel visualizations. Our analysis is facilitated by a coarse alignment between the lattice’s nominal model, a spatial graph, and the XCT data. We employ a simple VR-based approach for fast and rough alignment, followed by an offline registration and identification of the struts. With the nodes and struts aligned and identified in the volume, our framework allows querying of subvolumes containing a single strut at multiple resolutions. This avoids computation over the entire lattice and also allow for easy parallelization of down-stream computations, such as strut-specific metrics. To depict a fast overview of the strut quality, we introduce two innovative visual encodings, crucial for our collaborators’ research in creating novel AM parts: the Contour View and the Roughness Map, which depict critical geometrical and surface features of individual struts in standardized two 2D views. We evaluated the integrated system through expert interviews. The feedback confirms the framework’s practicality and its effectiveness in enhancing current inspection workflows. It solves major bottlenecks for our collaborators, ultimately helping them create novel parts with advanced properties.

Miao, Haichao [Lawrence Livermore National Laborat↗
DOE OSTI2024

Additively Manufactured, Lightweight, Low-Cost Composite Vessels for Compressed Natural Gas Fuel Storage

This project will develop a process to combine AM via direct ink writing (DIW) technology for CFC printing and to use design optimization tools pioneered at LLNL, with advances in resin/composite formulation enabled by chemical and nano-material modification to produce lightweight low-cost CNG tanks. Our approach will yield sub-scale prototype composite pressure tanks equivalent to Type-5 CNG vessel designs that demonstrate a potential cost-benefit advantage. Central to our vision is using agile AM and design based on computationally informed DIW of both short and continuous CF, further coupled with high-performance thermoset polymer matrixes modified by emergent nanomaterials. Our single-stage, multi-material AM technology, combined with a decreased volume fraction of CF and an increased proportion of economically advantaged short fiber, all together drive the reduction in manufacturing time and overall cost. Importantly, reductions in continuous fiber and overall fiber volume fraction will be achieved without detriment to the mechanical strength of the composite vessel. This will be achieved by employing a single process using multi materials grading involving a thermoset resin “ink” modified with aligned nanoplatelets to leverage the efficient tortuous-path gas barrier effect, printed as an inner flexible gas barrier as the initial stage in our manufacturing process before compositionally grading the AM feedstock in real-time to transition to a rigid, structural CF-filled resin. The proposed hybrid construction is projected to achieve pressure ratings at a service range of 2,900–3,600 psi with a 3× burst safety factor comparable to conventional filament-wound composite tanks with an estimated 30–50% reduction in total manufacturing cost.

03 NATURAL GAS↗
DOE OSTI2022

Chiral-like doublet band structure and octupole correlations in 104 Ag

The nature of the yrast negative-parity band and its chiral-like partner band in 104 Ag is investigated experimentally and theoretically. Lifetimes of states in the negative-parity yrast band and positive-parity band based on the 4424-keV level are measured using Doppler shift attenuation technique. Lifetimes of three more states have been determined along with the upper limit for the lifetime of the highest observed yrast states. Further, lifetimes known from earlier studies are determined with better precision. The level scheme of 104 Ag has also been extended with the addition of new enhanced E1 transitions linking the positive-parity band based on the 4424-keV levels and the yrast negative-parity and its partner band. B(E1) and/or B(E1)/B(M1) values for the transitions from the positive-parity band to the yrast and its partner band have been determined for the first time; these suggest strong octupole correlation between the positive-parity and the negative-parity bands. Here, calculations based on the triaxial projected shell model (TPSM) and covariant density functional theory (CDFT) have been performed to unravel the intrinsic structures of the partner band and the excited positive-parity band. TPSM calculations predict that doublet bands have significant angular momentum contributions along the three principle axes, suggesting that bands could have chiral symmetry breaking origin. The CDFT calculations predict a π(g 9/2 ) -1 Ⓧν(h 11/2 )(g 7/2 ,d 5/2 ) 2 aligned quasiparticle configuration for the negative-parity doublet bands with deformation parameters β≈0.20 and γ≈5°. The partner band could be interpreted as a chiral vibration mode built on top of the yrast band. The excited positive-parity band is predicted to have aligned four quasiparticle configurations, namely, π(g 9/2 ) -1 Ⓧν(h 11/2 ) 2 (g 7/2 ,d 5/2 ) 1 . Further, these calculations predict significant octupole softness in 104 Ag which could be the reason for enhanced E1 transitions between the four quasiparticle positive-parity bands and the doublet negative-parity bands.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2026

Hydration-Controlled Proton Transport in Respiratory Complex I

Proton pumping by respiratory complex I is one essential element for generating the proton motive force that drives ATP synthesis in mitochondria. Although it is understood that electrons from NADH reduce ubiquinone at the peripheral arm and that four protons are transferred in the membrane domain, the mechanism by which this redox reaction initiates proton translocation remains unclear. A lateral pathway linking the quinone binding site to the membrane domain via ND1, ND3, and ND4L subunits has been proposed as a possible initial path of an excess proton. However, experimental structures indicate that the hydration connectivity between D66 ND3 and E34 ND4L is comparatively weaker than in neighboring segments, suggesting a potential regulatory point for proton transfer. Using multiscale reactive molecular dynamics (MS-RMD) and a water wire connectivity metric, we directly simulate proton transport through this region as coupled to the hydration by water molecules. Our results reveal that proton transfer is thermodynamically feasible when transient hydration aligns with the presence of an excess proton, revealing the strong coupling between hydration and proton transfer (PT) in this region of Complex I. These findings support a model where proton injection enhances local hydration, dynamically opening the pathway for proton transfer and regulating the onset of proton pumping in Complex I.

Chemical structure↗
DOE OSTI2025

Characterization of thermally heat-treated polyacrylonitrile carbon fibers

This study investigates the graphitization process of polyacrylonitrile (PAN) carbon fibers by subjecting commercial fibers to thermal heat treatment at temperatures ranging from 1400 to 2100 °C in 100 °C increments, using either argon or nitrogen gas atmospheres. Changes in crystallinity, surface morphology, and lattice parameters were analyzed for two commercial carbon fibers using X-ray diffraction, scanning electron microscopy, and Raman spectroscopy. Results indicated minimal changes in surface morphology with increasing heat-treatment temperature; however, crystallinity significantly increased. Crystallinity changes were more strongly dependent on temperature rather than gas atmosphere or fiber type. At intermediate heat-treatment temperatures (1600–1800 °C), fibers treated in argon showed a slight preference for graphitization. The highest level of graphitization was measured at 2100 °C. Crystallite size increased as the intensity ratio of the D1 to G Raman peaks increased, reaching a peak around ~1800 °C, after which the ratio started to decrease. This behavior aligns with Ferrari's three-stage model of carbon crystallization and is consistent with both the Marie-Mering degree of graphitization and Brubaker's Integrated Absolute Differential models, all of which describe the transformation from an amorphous to a more graphitic structure. At the higher heat-treatment temperatures, the changes between atmospheres and fiber types were measured to converge to similar levels of graphitization. In conclusion, this study evaluates the progressive change in commercial grade carbon fibers when heat-treated.

Characterization↗
DOE OSTI2024

Elucidating the Role of Carbon Conductive Additive in the Processing and Electrochemical Behavior of Surface-Modified Si Anodes

In the present work, we explore the impact of conducting carbon on the processing of silicon (Si)-rich electrodes, their resulting cycling performance, and parasitic side reactions. We employed three different carbon additives, Super C45 (carbon black), Super C65 (carbon black), and TIMCAL C-NERGY KS 6L (graphite), with varying sizes and morphologies to determine how these parameters influence the structure and the electrochemical behavior of the cast electrodes. Raman mapping indicates that the improved performance achieved using C45 can be explained by the homogeneous distribution of Si and carbon, enabling the formation of continuous electrical pathways throughout the electrode. Further, this finding aligns with the results of ζ potential measurements, which indicate that C45 can maintain more stable dispersions in N-methyl-2-pyrrolidone (NMP) compared with the other carbon additives. However, electrodes with C45 as the conducting agent exhibited increased parasitic side reactions, as evidenced by the leakage current obtained from voltage hold experiments. These side reactions can have adverse effects on the calendar life of the electrode. Furthermore, the Raman maps of the electrodes reveal heterogeneities in Si and carbon distribution with lower carbon content, especially with only 2% carbon. This drastically impacts the cyclability of the electrodes, with portions of the low-carbon electrodes that remain unable to undergo cycling even at slow rates (0.2 Ah g –1 ). Additionally, we found that slower formation rates are favored in the context of parasitic side reactions, although slower formation rates can potentially create thicker passivating films that could hinder Li + -ion transport. Therefore, we emphasize the need for an intricate balance in materials selection and cycling protocol for optimal electrode processing and electrochemical behavior.

25 ENERGY STORAGE↗
DOE OSTI2022

Observation of an Orientational Glass in a Superlattice of Elliptically-Faceted CdSe Nanocrystals

Extensive prior work has shown that colloidal inorganic nanocrystals coated with organic ligand shells can behave as artificial atoms and, as such, form superlattices with different crystal structures and packing densities. Although ordered superlattices present a high degree of long-range positional order, the relative crystallographic orientation of the inorganic nanocrystals with respect to each other tends to be random. Recent works have shown that superlattices can achieve orientational alignment through combinations of nanocrystal faceting and ligand modification, as well as selective metal particle attachment to particular facets. These studies have focused on the assembly of high-symmetry nanocrystals, such as cubes and cuboctahedra. In this work, we study the assembly of elliptically faceted CdSe/CdS core/shell nanocrystals with one distinctive crystallographic orientation along the major elliptical axis. We show that the nanocrystals form an unexpectedly well-ordered translational superlattice, with a degree of order comparable to that achieved with higher-symmetry nanocrystals. Additionally, we show that, due to the particles' faceted shape, the superlattice is characterized by an orientational glass phase in which only certain orientations are possible due to entropically frustrated crystallization. In this phase, the nanocrystals do not exhibit a local orientational ordering but rather have distinct orientations that emerge at different locations within the same domain. The distinct orientations are a result of a facet-to-facet lock-in mechanism that occurs during the self-assembly process. These facet-to-facet alignments force the nanocrystals to tilt on different lattice planes forming different projections that we termed apparent polydispersity. Our experimental realization of an orientational glass phase for multifaceted semiconducting nanocrystals can be used to investigate how this phase is formed and how it can be utilized for potential optical, electrical, and thermal transport applications.

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

Disentangling the Impacts of Microtopography and Shrub Distribution on Snow Depth in a Subarctic Watershed: Toward a Predictive Understanding of Snow Spatial Variability

Snow plays a critical role in carbon cycling, vegetation dynamics, and permafrost hydrology at high latitudes by influencing surface energy exchange. Predicting snow distribution patterns is essential for understanding the evolution of Arctic ecosystems, yet scaling process-level knowledge to landscape predictions remains challenging. Here, we analyze snow depth (2019 and 2022), terrain elevation, and vegetation height from a watershed on the Seward Peninsula, Alaska, to examine how topography and shrubs shape snow redistribution across spatial scales. We find that snow depth is strongly coupled to terrain at scales below ∼60 m but becomes increasingly decoupled at larger scales. The topographic model of snow depth variation, which transforms terrain data to align with these scale-dependent snow patterns, is well correlated with local snow depth variations (linear fit R 2 > 0.5 for 85% of 100-m patches). A machine learning reconstruction of shrub canopy snow trapping reveals a simple exponential relationship between canopy structure and snow accumulation ( R 2 = 0.59), highlighting the combined influence of topography and vegetation on snow distribution. Together, these empirical relationships capture much of the observed snow variability in the watershed ( R 2 = 0.49, root mean square error (RMSE) = 30 cm), though systematic limitations persist in areas of strong scour and at coarser scales where wind-terrain interactions are more complex. These findings provide a framework for more efficient snow depth prediction and offer insights to improve snow-vegetation feedback representation in Earth System Models.

54 ENVIRONMENTAL SCIENCES↗