Search NASA⌕ Search

SEARCH · Search NASA

Results for “Heterogeneous Computing”

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 703 records · Page 39
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Using the ATLAS experiment software on heterogeneous resources

With the large dataset expected from 2030 onwards by the HL-LHC at CERN, the ATLAS experiment is reaching the limits of the current data processing model in terms of traditional CPU resources based on x86_64 architectures and an extensive program for software upgrades towards the HL-LHC has been set up. The ARM CPU architecture is becoming a competitive and energy efficient alternative. Accelerators like GPUs are available in any recent HPC. In the past years ATLAS has successfully ported its full data processing and simulation software framework Athena to ARM and has invested significant effort in porting parts of the reconstruction and simulation algorithms to GPUs. We report on the successful usage of the ATLAS experiment offline and online software framework Athena on ARM and GPUs through the PanDA workflow management system at various WLCG sites. Furthermore we report on performance optimizations of the builds for ARM CPUs and the GPU integration efforts. We will discuss performance comparisons of different ARM and x86_64 architectures on WLCG resources and Cloud compute providers like GCP and AWS using ATLAS productions workflows as used in the Hep-Score23 benchmark suite.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Adsorption, Orientation, and Speciation of Amino Acids at Air–Aqueous Interfaces for the Direct Air Capture of CO 2

Amino acids make up a promising family of molecules capable of direct air capture (DAC) of CO 2 from the atmosphere. Under alkaline conditions, CO 2 reacts with the anionic form of an amino acid to produce carbamates and deactivated zwitterionic amino acids. The presence of the various species of amino acids and reactive intermediates can have a significant effect on DAC chemistry, the role of which is poorly understood. In this study, all-atom molecular dynamics (MD) based computational simulations and vibrational sum frequency generation (vSFG) spectroscopy studies were conducted to understand the role of competitive interactions at the air–aqueous interface in the context of DAC. We find that the presence of potassium bicarbonate ions, in combination with the anionic and zwitterionic forms of amino acids, induces concentration and charge gradients at the interface, generating a layered molecular arrangement that changes under pre- and post-DAC conditions. In parallel, an enhancement in the surface activity of both anionic and zwitterionic forms of amino acids is observed, which is attributed to enhanced interfacial stability and favorable intermolecular interactions between the adsorbed amino acids in their anionic and zwitterionic forms. The collective influence of these competitive interactions, along with the resulting interfacial heterogeneity, may in turn affect subsequent capture reactions and associated rates. Finally, these effects underscore the need to consider dynamic changes in interfacial chemical makeup to enhance DAC efficiency and to develop successful negative emission and storage technologies.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRSDate not supplied

Spatial and Temporal Deconfliction of Trajectories in the Presence of Uncertainties

Demonstration of conflict-free movement for multi-agent teams in challenging scenarios is crucial in developing trust and trustworthiness in an autonomous transport system. Tolerance verification queries are explored as a mechanism to enforce spatial and temporal deconfliction for a cooperating team of Unmanned Aerial Systems (UAS) with prescribed heterogeneous path-following performance. Obstacles in the environment are modelled as set of polyhedra, whereas each vehicle’s trajectory is represented as a sequence of polynomial curves with C2 continuity, expressed in a Bernstein basis. Each vehicle is modelled as a point mass and a safety distance, informed by the geometry of the UAS and the worst-case path-following error. This defines a tube around the trajectories where each UAS is most likely to fly through. In addition, obstacles in the environment have an associated safety buffer around them to account for the uncertainty in their location and geometric description. The tolerance verification queries explored in this paper combine the safety distance information from each UAS and environmental hazard to compute trajectories that are contained within the safe configuration space. Tolerance verification is also compared with other proximity queries to determine the suitability of each method along the different steps of the trajectory generation algorithm. This paper analyzes the fitness and performance of three proximity queries – collision, tolerance verification, and distance computations – between polyhedral and polynomial curves to ensure deconfliction between obstacles and vehicles, but also between polynomial curves to guarantee safe separation among cooperating UAS.

trajectory generation↗
DOE OSTI2024

Active doping controls the mode of failure in dense colloidal gels

Mechanical properties of disordered materials are governed by their underlying free energy landscape. In contrast to external fields, embedding a small fraction of active particles within a disordered material generates nonequilibrium internal fields, which can help to circumvent kinetic barriers and modulate the free energy landscape. In this work, we investigate through computer simulations how the activity of active particles alters the mechanical response of deeply annealed polydisperse colloidal gels. We show that the “swim force” generated by the embedded active particles is responsible for determining the mode of mechanical failure, i.e., brittle vs. ductile. We find, and theoretically justify, that at a critical swim force the mechanical properties of the gel decrease abruptly, signaling a change in the mode of mechanical failure. The weakening of the elastic modulus above the critical swim force results from the change in gel porosity and distribution of attractive forces among gel particles, while below the critical swim force, the ductility enhancement is caused by an increase of gel structural disorder. Above the critical swim force, the gel develops a pronounced heterogeneous structure characterized by multiple pore spaces, and the mechanical response is controlled by dynamical heterogeneities. We contrast these results with those of a simulated monodisperse gel that exhibits a nonmonotonic trend of ductility modulation with increasing swim force, revealing a complex interplay between the gel energy landscape and embedded activity.

Zhou, Tingtao (ORCID:000000021766719X)↗
DOE OSTI2024

Arbitrary Low-Dimensional Film Transfer Enabled by GeO 2 Release Layer

Low-dimensional materials show great promise for enhanced computing and sensing performance in mission-relevant environments. However, integrating low-dimensional materials into conventional electronics remains a challenge. Here, we demonstrate a novel transfer method by which low-dimensional materials and their heterostructures can be transferred onto any arbitrary substrate. Our method relies on a water soluble GeO 2 substrate from which lowdimensional materials are transferred without significant perturbation. We apply the method to transfer a working electronic device based on a low-dimensional material. Process developments are achieved to enable the fabrication and transfer of a working electronic device, including the growth of high-k dielectric on GeO 2 by atomic layer deposition and inserting an indium diffusion barrier into the device gate stack. This work supports Sandia’s heterogeneous integration strategy to broaden the implementation of low-dimensional films and their devices.

36 MATERIALS SCIENCE↗
NASA NTRS1993

Development of deterministic transport methods for low energy neutrons for shielding in space

Transport of low energy neutrons associated with the galactic cosmic ray cascade is analyzed in this dissertation. A benchmark quality analytical algorithm is demonstrated for use with BRYNTRN, a computer program written by the High Energy Physics Division of NASA Langley Research Center, which is used to design and analyze shielding against the radiation created by the cascade. BRYNTRN uses numerical methods to solve the integral transport equations for baryons with the straight-ahead approximation, and numerical and empirical methods to generate the interaction probabilities. The straight-ahead approximation is adequate for charged particles, but not for neutrons. As NASA Langley improves BRYNTRN to include low energy neutrons, a benchmark quality solution is needed for comparison. The neutron transport algorithm demonstrated in this dissertation uses the closed-form Green's function solution to the galactic cosmic ray cascade transport equations to generate a source of neutrons. A basis function expansion for finite heterogeneous and semi-infinite homogeneous slabs with multiple energy groups and isotropic scattering is used to generate neutron fluxes resulting from the cascade. This method, called the FN method, is used to solve the neutral particle linear Boltzmann transport equation. As a demonstration of the algorithm coded in the programs MGSLAB and MGSEMI, neutron and ion fluxes are shown for a beam of fluorine ions at 1000 MeV per nucleon incident on semi-infinite and finite aluminum slabs. Also, to demonstrate that the shielding effectiveness against the radiation from the galactic cosmic ray cascade is not directly proportional to shield thickness, a graph of transmitted total neutron scalar flux versus slab thickness is shown. A simple model based on the nuclear liquid drop assumption is used to generate cross sections for the galactic cosmic ray cascade. The ENDF/B V database is used to generate the total and scattering cross sections for neutrons in aluminum. As an external verification, the results from MGSLAB and MGSEMI were compared to ANISN/PC, a routinely used neutron transport code, showing excellent agreement. In an application to an aluminum shield, the FN method seems to generate reasonable results.

Ganapol, Barry↗
DOE OSTI2025

LABQ3: Bayesian method for quantification of mineral compositions and nano-scale elemental mapping of 3D synchrotron XCT data

Quantitative analysis of mineral compositions is essential in understanding geochemical, mineralogical and environmental processes. Fine-resolution 3D imaging is widely done using synchrotron X-ray computed tomography (XCT), but existing analyses are limited to visualization and segmentation. This paper presents a new method, Linear Attenuation Bayesian Quantitative 3D-mapper (LABQ3), based on the linearity of X-ray attenuation with respect to elemental concentrations. To address the random variability in attenuation measurements, LABQ3 employs Bayesian decision theory to minimize classification error, using reference attenuation distributions from scans of pure mineral standards. To demonstrate LABQ3 and test its performance, we studied precipitated carbonate samples. XCT scans were done at multiple energies using the transmission X-ray microscope (TXM) at beamline 32-ID-C of the Advanced Photon Source at Argonne National Laboratory. The reconstructed 3D images have a voxel size of 20 nm. Analyses revealed rich nano-scale compositional heterogeneity within individual particles. A mixture of calcium and cadmium produced an overall stoichiometric composition of (Ca 0.78 ,Cd 0.22 )CO 3 , with some voxels containing nearly pure CdCO 3 . The addition of zinc led to an overall stoichiometric composition of 33% Ca, 28% Cd, 39% Zn, with a nearly pure CaCO 3 core and compositional zonation through the rim. These compositional gradients are related to temporal sequences of carbonate mineral formation where Cd precipitated at the beginning in (Ca,Cd)CO 3 , while Cd and Zn precipitated at the end in (Ca, Cd,Zn)CO 3 . Results differ from bulk analyses using Inductively Coupled Plasma-Mass Spectrometry (ICP-MS), showing that LABQ3 provides particle-specific insights. LABQ3 distinguishes itself by quantifying chemical compositions along a continuum, making it different from XCT analyses based on segmentation. LABQ3 allows simultaneous acquisition of morphology and chemical composition in 3D, facilitating the interpretation of chemical gradients of trace elements, quantification of solid solution compositions, inferences about temporal sequences of mineral precipitation, and addressing other concerns about solid-phase chemistry.

58 GEOSCIENCES↗
NASA NTRS1976

A satellite-borne ion mass spectrometer for the energy range 0 to 16 keV

The Ion Composition Experiment (ICE) on GEOS represents the first comprehensive attempt to measure the positive ion composition at high altitudes in the magnetosphere. Due to the heterogeneous nature of the magnetospheric plasma a novel mass spectrometer has been developed to cover the mass per charge range from H-1(+) to beyond Ba-138(+) and the energy per charge range from 0 to 16 keV/e. The ICE consists primarily of a cylindrical electrostatic analyzer followed by a curved analyzer incorporating crossed magnetic and electric fields. This combination has limited angular and energy focusing properties, but it maintains a mass resolution of about 4 over a wide range in energy and mass, sufficient for the objectives of measuring plasmas of both solar and terrestrial origin. High sensitivity and low background should allow measurements of rarer ion constituents down to flux levels of 0.01 ions/sq cm sec ster eV. A sophisticated electronics combined with powerful ground computer and telecommand systems allow for very efficient scanning of the mass-energy space.

Balsiger, H.↗
DOE OSTI2026

Regional inertia dynamics of U.S. interconnections: An event-based measurement approach

Power grid inertia plays a vital role in frequency stability following large disturbances, yet its distribution across the U.S. grid is highly uneven. While interconnection-wide inertia benchmarks are useful, they can mask regional variability driven by resource mix, network coupling, and geographic separation. This paper extends event-driven inertia estimation to the regional scale using field measurements from the Frequency Monitoring Network (FNET/GridEye). Starting from balancing authority and independent system operator footprints, candidate regions are refined using a composite coherency score that combines frequency-trajectory shape similarity, timing spread, and lead/lag behavior to ensure dynamic consistency. A filtered sliding difference method (FSDM) is then used to construct regional frequency trajectories, detect disturbance onset, and compute robust regional rate-of-change of frequency (RoCoF). Regional, local, and interconnection inertia are estimated by combining RoCoF with event power imbalance, and additional indicators (regional-to-system inertia ratio and inertial-support arrival time) quantify regional-to-interconnection coupling and relative regional contributions. The method is demonstrated on eleven regions across the Eastern Interconnection (EI) and the Western Electricity Coordinating Council (WECC), with the Electric Reliability Council of Texas (ERCOT) used for validation. In ERCOT, estimates compared against energy management system (EMS) values achieve a mean absolute percentage error of 17.94%. WECC exhibits consistently shorter inertial-support arrival times (0.15–0.3 s) than EI (0.7–1.1 s), highlighting contrasting coupling and disturbance-propagation behavior. Overall, the results reveal pronounced spatial heterogeneity in inertia and coupling, underscoring the value of regional monitoring for both operational decision-making and long-term system planning.

Disturbance events↗
NASA NTRS1972

Studies in interactive communication. I - The effects of four communication modes on the behavior of teams during cooperative problem-solving.

Two-man teams solved credible, 'real-world' problems for which computer assistance has been or could be useful. Conversations were carried on in one of four modes of communication: (1) typewriting, (2) handwriting, (3) voice, and (4) natural, unrestricted communication. Two groups of subjects (experienced and inexperienced typists) were tested in the typewriting mode. Performance was assessed on three classes of dependent measures: time to solution, behavioral measures of activity, and linguistic measures. Significant and meaningful differences among the communication modes were found in each of the three classes of dependent variable. This paper is concerned mainly with the results of the activity analyses. Behavior was recorded in 15 different categories. The analyses of variance yielded 34 statistically significant terms of which 27 were judged to be practically significant as well. When the data were transformed to eliminate heterogeneity, the analyses of variance yielded 35 statistically significant terms of which 26 were judged to be practically significant.

Chapanis, A.↗
DOE OSTI2024

Development of Whole System Digital Twins for Advanced Reactors: Leveraging Graph Neural Networks and SAM Simulations

Here, in this work, we introduce a novel method to develop whole system digital twins (DTs) for advanced nuclear reactors. This method treats a complex reactor system as a heterogeneous graph: with the system components as different types of graph nodes and their physical interconnections as edges. Based on the heterogeneous graph, a graph neural network combining graph convolution and temporal node attention is developed as the DT, facilitating a comprehensive understanding of the system's dynamic behavior. By utilizing the System Analysis Module (SAM) code for simulating various operational transients, we develop a graph-based database that trains the DT. This DT is characterized by two primary functions: It can infer the entire system's status using sparse node information, and it can predict the progress of transients based on current and historical system information. Our approach is validated through case studies on the Experimental Breeder Reactor II (EBR-II) system and a generic Fluoride-salt-cooled High-temperature Reactor (gFHR), demonstrating the DT's accuracy in forecasting operational transients. The DT's rapid computation capabilities enhance its potential for supporting advanced reactor operations, offering benefits in intelligent simulation, autonomous control, and anomaly detection, paving the way for improved safety analysis and intelligent component health management for advanced reactor systems and reducing their operations and maintenance cost.

EBR-II↗
DOE OSTI2024

Upsampling Monte Carlo Reactor Simulation Tallies in Depleted Sodium-Cooled Fast Reactor Assemblies Using a Convolutional Neural Network

The computational demand of neutron Monte Carlo transport simulations can increase rapidly with the spatial and energy resolution of tallied physical quantities. Convolutional neural networks have been used to increase the resolution of Monte Carlo simulations of light water reactor assemblies while preserving accuracy with negligible additional computational cost. Here, we show that a convolutional neural network can also be used to upsample tally results from Monte Carlo simulations of sodium-cooled fast reactor assemblies, thereby extending the applicability beyond thermal systems. The convolutional neural network model is trained using neutron flux tallies from 300 procedurally generated nuclear reactor assemblies simulated using OpenMC. Validation and test datasets included 16 simulations of procedurally generated assemblies, and a realistic simulation of a European sodium-cooled fast reactor assembly was included in the test dataset. We show the residuals between the high-resolution flux tallies predicted by the neural network and high-resolution Monte Carlo tallies on relative and absolute bases. The network can upsample tallies from simulations of fast reactor assemblies with diverse and heterogeneous materials and geometries by a factor of two in each spatial and energy dimension. The network’s predictions are within the statistical uncertainty of the Monte Carlo tallies in almost all cases. This includes test assemblies for which burnup values and geometric parameters were well outside the ranges of those in assemblies used to train the network.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
NASA NTRS1997

A regional ultrastructural analysis of the cellular and synaptic architecture in the chinchilla cristae ampullares

The chinchilla crista ampullaris was studied in 10 samples, each containing 32 consecutive ultrathin sections of the entire neuroepithelium. Dissector methods were used to estimate the incidence of various synaptic features, and results were confirmed in completely reconstructed hair cells. There are large regional variations in cellular and synaptic architecture. Type I and type II hair cells are shorter, broader, and less densely packed in the central zone than in the intermediate and peripheral zones. Complex calyx endings are most common centrally. On average, there are 15-20 ribbon synapses and 25-30 calyceal invaginations in each type I hair cell. Synapses and invaginations are most numerous centrally. Central type II hair cells receive considerably fewer afferent boutons than do peripheral type II hair cells, but have similar numbers of ribbon synapses. The numbers are similar because central type II hair cells make more synapses with the outer faces of calyx endings and with individual afferent boutons. Most afferent boutons get one ribbon synapse. Boutons without ribbon synapses were only found peripherally, and boutons getting multiple synapses were most frequent centrally. Throughout the neuroepithelium, there is an average of three to four efferent boutons on each type II hair cell and calyx ending. Reciprocal synapses are rare. Most synaptic ribbons in type I hair cells are spherules; those in type II hair cells can be spherical or elongated and are particularly heterogeneous centrally. Consistent with the proposal that the crista is concentrically organized, the intermediate and peripheral zones are each similar in their cellular and synaptic architecture near the base and near the planum. An especially differentiated subzone may exist in the middle of the central zone.

Non-NASA Center↗
NASA NTRS1998

Geology of Lunar Landing Sites and Origin of Basin Ejecta from a Clementine Perspective

The goals of this research were to examine Clementine multispectral data covering the Apollo landing sites in order to: (1) provide ground truth for the remotely sensed observations, (2) extend our understanding of the Apollo landing sites to the surrounding regions using the empirically calibrated Clementine data, and (3) investigate the composition and distribution of impact-basin ejecta using constraints based upon the remotely sensed data and the Apollo samples. Our initial efforts (in collaboration with P. Lucey and coworkers) to use the Apollo soil compositions to "calibrate" information derived from the remotely sensed data resulted in two extremely useful algorithms for computing estimates of the concentrations of FeO and TiO2 from the UV-VIS 5-band data. In this effort, we used the average surface soil compositions from 37 individual Apollo and 3 Luna sample stations that could be resolved using the Clementine data. We followed this work with a detailed investigation of the Apollo 17 landing site, where the sampling traverses were extensive and the spectral and compositional contrast between different soils covers a wide range. We have begun to investigate the nature and composition of basin ejecta by comparing the thick deposits on the rim of Imbrium in the vicinity of the Apollo 15 site and those occurring southeast of the Serenitatis basin, in the Apollo 17 region. We continue this work under NAG5-6784, "Composition, Lithology, and Heterogeneity of the lunar crust using remote sensing of impact-basin uplift structures and ejecta as probes. The main results of our work are given in the following brief summaries of major tasks. Detailed accounts of these results are given in the attached papers, manuscripts, and extended abstracts.

Jolliff, Bradley L.↗
DOE OSTI2026

Solvent Stabilization of Protic Oxonium/Ammonium Intermediates in Cation Radical Cyclization Reactions Investigated via Computational Approaches

Anodic intramolecular cyclization reactions have substantial synthetic utility for formation of cyclic carbon–carbon or carbon–heteroatom bonds. For cases of intramolecular trapping of a cation radical by a protic nucleophile, the cyclization step coincides with a substantial increase in substrate acidity and thus may exhibit particularly pronounced solvent effects. In this computational work, we employ both quantum chemical (QM) and quantum mechanics/molecular mechanics (QM/MM) methods to compute solvent effects on free energy profiles for cyclization and deprotonation reaction steps for cation radical intermediates of substrates representative for anodic intramolecular cyclizations. We find substantial solvent contribution to the thermodynamic driving force for cation radical cyclization; for example, methanol and tetrahydrofuran solvents provide ∼30–35 kJ/mol driving force to form cyclic oxonium cation radicals and ∼15–25 kJ/mol driving force to form cyclic ammonium cation radicals, compared to baseline reactions in dichloromethane solvent. Given that these solvent shifts are on par with the innate cyclization reaction thermodynamics, the choice of solvent plays a crucial role in promoting/driving the cation radical cyclization step. Methanol is particularly effective at facilitating rapid deprotonation of the cyclic cation radical intermediate, which may lead to the full electrochemical process (e.g., second electron transfer) proceeding heterogeneously at the anode.

Khan, Shahriar [ORNL] (ORCID:0000000289138430)↗
NASA NTRS1994

Methodologies and Tools for Tuning Parallel Programs: Facts and Fantasies

The need for computing power has forced a migration from serial computation on a single processor to parallel processing on multiprocessors. However, without effective means to monitor (and analyze) program execution, tuning the performance of parallel programs becomes exponentially difficult as program complexity and machine size increase. The recent introduction of performance tuning tools from various supercomputer vendors (Intel's ParAide, TMC's PRISM, CRI's Apprentice, and Convex's CXtrace) seems to indicate the maturity of performance tool technologies and vendors'/customers' recognition of their importance. However, a few important questions remain: What kind of performance bottlenecks can these tools detect (or correct)? How time consuming is the performance tuning process? What are some important technical issues that remain to be tackled in this area? This workshop reviews the fundamental concepts involved in analyzing and improving the performance of parallel and heterogeneous message-passing programs. Several alternative strategies will be contrasted, and for each we will describe how currently available tuning tools (e.g. AIMS, ParAide, PRISM, Apprentice, CXtrace, ATExpert, Pablo, IPS-2) can be used to facilitate the process. We will characterize the effectiveness of the tools and methodologies based on actual user experiences at NASA Ames Research Center. Finally, we will discuss their limitations and outline recent approaches taken by vendors and the research community to address them.

Yan, Jerry C.↗
NASA NTRS1994

Performance Evaluation Methodologies and Tools for Massively Parallel Programs

The need for computing power has forced a migration from serial computation on a single processor to parallel processing on multiprocessors. However, without effective means to monitor (and analyze) program execution, tuning the performance of parallel programs becomes exponentially difficult as program complexity and machine size increase. The recent introduction of performance tuning tools from various supercomputer vendors (Intel's ParAide, TMC's PRISM, CSI'S Apprentice, and Convex's CXtrace) seems to indicate the maturity of performance tool technologies and vendors'/customers' recognition of their importance. However, a few important questions remain: What kind of performance bottlenecks can these tools detect (or correct)? How time consuming is the performance tuning process? What are some important technical issues that remain to be tackled in this area? This workshop reviews the fundamental concepts involved in analyzing and improving the performance of parallel and heterogeneous message-passing programs. Several alternative strategies will be contrasted, and for each we will describe how currently available tuning tools (e.g., AIMS, ParAide, PRISM, Apprentice, CXtrace, ATExpert, Pablo, IPS-2)) can be used to facilitate the process. We will characterize the effectiveness of the tools and methodologies based on actual user experiences at NASA Ames Research Center. Finally, we will discuss their limitations and outline recent approaches taken by vendors and the research community to address them.

Yan, Jerry C.↗
NASA NTRS2014

Assimilation of MODIS Snow Cover Through the Data Assimilation Research Testbed and the Community Land Model Version 4

To improve snowpack estimates in Community Land Model version 4 (CLM4), the Moderate Resolution Imaging Spectroradiometer (MODIS) snow cover fraction (SCF) was assimilated into the Community Land Model version 4 (CLM4) via the Data Assimilation Research Testbed (DART). The interface between CLM4 and DART is a flexible, extensible approach to land surface data assimilation. This data assimilation system has a large ensemble (80-member) atmospheric forcing that facilitates ensemble-based land data assimilation. We use 40 randomly chosen forcing members to drive 40 CLM members as a compromise between computational cost and the data assimilation performance. The localization distance, a parameter in DART, was tuned to optimize the data assimilation performance at the global scale. Snow water equivalent (SWE) and snow depth are adjusted via the ensemble adjustment Kalman filter, particularly in regions with large SCF variability. The root-mean-square error of the forecast SCF against MODIS SCF is largely reduced. In DJF (December-January-February), the discrepancy between MODIS and CLM4 is broadly ameliorated in the lower-middle latitudes (2345N). Only minimal modifications are made in the higher-middle (4566N) and high latitudes, part of which is due to the agreement between model and observation when snow cover is nearly 100. In some regions it also reveals that CLM4-modeled snow cover lacks heterogeneous features compared to MODIS. In MAM (March-April-May), adjustments to snowmove poleward mainly due to the northward movement of the snowline (i.e., where largest SCF uncertainty is and SCF assimilation has the greatest impact). The effectiveness of data assimilation also varies with vegetation types, with mixed performance over forest regions and consistently good performance over grass, which can partly be explained by the linearity of the relationship between SCF and SWE in the model ensembles. The updated snow depth was compared to the Canadian Meteorological Center (CMC) data. Differences between CMC and CLM4 are generally reduced in densely monitored regions.

data assimilation↗