Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Matrix factorization”

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

Machine learning to identify geologic factors associated with production in geothermal fields: a case-study using 3D geologic data, Brady geothermal field, Nevada

Abstract In this paper, we present an analysis using unsupervised machine learning (ML) to identify the key geologic factors that contribute to the geothermal production in Brady geothermal field. Brady is a hydrothermal system in northwestern Nevada that supports both electricity production and direct use of hydrothermal fluids. Transmissive fluid-flow pathways are relatively rare in the subsurface, but are critical components of hydrothermal systems like Brady and many other types of fluid-flow systems in fractured rock. Here, we analyze geologic data with ML methods to unravel the local geologic controls on these pathways. The ML method, non-negative matrix factorization with k -means clustering (NMF k ), is applied to a library of 14 3D geologic characteristics hypothesized to control hydrothermal circulation in the Brady geothermal field. Our results indicate that macro-scale faults and a local step-over in the fault system preferentially occur along production wells when compared to injection wells and non-productive wells. We infer that these are the key geologic characteristics that control the through-going hydrothermal transmission pathways at Brady. Our results demonstrate: (1) the specific geologic controls on the Brady hydrothermal system and (2) the efficacy of pairing ML techniques with 3D geologic characterization to enhance the understanding of subsurface processes.

58 GEOSCIENCES↗
DOE OSTI2024

Investigating New Particle Formation and Growth Over an Urban Location in the Eastern Mediterranean

Abstract This study investigates the new particle formation (NPF) events at an urban location in the Eastern Mediterranean. Particle size distribution, particulate chemical composition, and gaseous pollutants were monitored in Rehovot, Israel (31°53″N 34°48″E) during two campaigns: from April 29 to 3 May 2021 (Campaign 1) and from May 3 to 11 May 2023 (Campaign 2), coinciding with an intensive bonfire burning festival. The organic aerosols (OA) source apportionment identified two major factors—Hydrocarbon‐like OA and Biomass‐burning OA—as well as two secondary factors—MO‐OOA (more oxidized‐oxygenated OA) and LO‐OOA (low oxidized oxygenated OA). NPF events were frequently observed during the day (mostly well‐defined nucleation events) and at night (burst of ultrafine mode particles without any discernible growth). A condensation sink value of (9.4 ± 4.0) × 10 −3 s −1 during Campaign 1 and (14.2 ± 6.0) × 10 −3 s −1 during Campaign 2 was obtained. The daytime events were associated with enhanced sulfuric acid proxy concentrations of (2–12) × 10 6 molecules cm −3 , suggesting the role of gas‐phase photochemistry in promoting NPF. A novel approach of hybrid positive matrix factorization analysis was used to deconvolve the chemical species responsible for the observed events. The results suggest the involvement of multiple components, including ammonium sulfate and MO‐OOA, in the nucleation; Nitrate, HOA and LO‐OOA participate in the subsequent particle growth for the daytime events. Nighttime events involve only semi‐volatile species (LO‐OOA, HOA and nitrate) along with ammonium sulfate.

Meteorology & Atmospheric Sciences↗
DOE OSTI2022

Submicrometer spectromicroscopy of UO 2 aged under high humidity conditions

We report the oxidation of uranium dioxide is a complicated process, depending on factors including humidity, temperature, and microstructure. To further determine the characteristics of this process, UO 2 particles were allowed to age and agglomerate under 98% relative humidity at room temperature for 378 days. A focused ion beam (FIB) section of this agglomeration was then measured at the O K-edge, U N 5 -edge, and C K-edge using the scanning transmission x-ray microscope (STXM) at the Advanced Light Source. O K-edge and U N 5 -edge x-ray absorption measurements allowed for the elemental and chemical species mapping of the agglomerates and indicated the formation of schoepite at the submicrometer scale in specific locations. Non-negative matrix factorization was employed to elucidate the main components at the O K-edge, which were uranyl (schoepite) formed primarily at the interface of the sample with controlled atmosphere, a UO 2 -like bulk component present in the majority of the sample, and an oxygen species present at the surface of the FIB section, which is likely adsorbed water. STXM spectromicroscopy measurements at the U N 5 -edge measurements also confirmed the location of oxidized uranium. This analysis is a valuable insight into the formation of schoepite on UO 2 and shows the sensitivity to and utility of STXM spectromicroscopy for uranium speciation.

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

Compression of tokamak boundary plasma simulation data using a maximum volume algorithm for matrix skeleton decomposition

This report demonstrates satisfactory data compression of SOLPS-ITER simulation output ranging from 2D fields, 1D profiles, and 0D scalar variables with a novel matrix decomposition approach. The singular value decomposition (SVD) scales poorly for large matrix sizes and is unsuited to the application on high dimensional data common to fusion plasma physics simulation. In this work, we employ the columns-submatrix-rows (CUR) matrix factorization technique in order to compute a low-rank approximation up to two orders of magnitude faster than the SVD, but within a nominal L2-norm relative error of ε = 10 –2 . In addition, the CUR approach maintains the original format of the data, in its extracted columns and rows, allowing for interpretable data storage at the original resolution of the simulation. We utilize an iterative algorithm to compute the CUR decomposition of simulation output by maximizing the volume, or linearly independent information content, of a low-rank submatrix contained within the data. Experiments over $\textit{n} × \textit{n}$ randomized test matrices with embedded rank-deficient features show that this maximum volume implementation of CUR matrix approximation has reduced asymptotic computational complexity on the order of n compared to the SVD, which scales approximately as $n^3$. These results show that the CUR technique can be used to effectively select time step snapshots (columns) of over 140 SOLPS-ITER output variables and the associated discretized coordinate timeseries (rows) allowing for reconstruction of the complete simulation dynamics.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS2011

Graph Theory Roots of Spatial Operators for Kinematics and Dynamics

Spatial operators have been used to analyze the dynamics of robotic multibody systems and to develop novel computational dynamics algorithms. Mass matrix factorization, inversion, diagonalization, and linearization are among several new insights obtained using such operators. While initially developed for serial rigid body manipulators, the spatial operators and the related mathematical analysis have been shown to extend very broadly including to tree and closed topology systems, to systems with flexible joints, links, etc. This work uses concepts from graph theory to explore the mathematical foundations of spatial operators. The goal is to study and characterize the properties of the spatial operators at an abstract level so that they can be applied to a broader range of dynamics problems. The rich mathematical properties of the kinematics and dynamics of robotic multibody systems has been an area of strong research interest for several decades. These properties are important to understand the inherent physical behavior of systems, for stability and control analysis, for the development of computational algorithms, and for model development of faithful models. Recurring patterns in spatial operators leads one to ask the more abstract question about the properties and characteristics of spatial operators that make them so broadly applicable. The idea is to step back from the specific application systems, and understand more deeply the generic requirements and properties of spatial operators, so that the insights and techniques are readily available across different kinematics and dynamics problems. In this work, techniques from graph theory were used to explore the abstract basis for the spatial operators. The close relationship between the mathematical properties of adjacency matrices for graphs and those of spatial operators and their kernels were established. The connections hold across very basic requirements on the system topology, the nature of the component bodies, the indexing schemes, etc. The relationship of the underlying structure is intimately connected with efficient, recursive computational algorithms. The results provide the foundational groundwork for a much broader look at the key problems in kinematics and dynamics. The properties of general graphs and trees of nodes and edge were examined, as well as the properties of adjacency matrices that are used to describe graph connectivity. The nilpotency property of such matrices for directed trees was reviewed, and the adjacency matrices were generalized to the notion of block weighted adjacency matrices that support block matrix elements. This leads us to the development of the notion of Spatial Kernel Operator SKO kernels. These kernels provide the basis for the development of SKO resolvent operators.

Jain, Abhinandan↗
DOE OSTI2022

Latent Mechanisms of Polarization Switching from In Situ Electron Microscopy Observations

In situ scanning transmission electron microscopy enables observation of the domain dynamics in ferroelectric materials as a function of externally applied bias and temperature. The resultant data sets contain a wealth of information on polarization switching and phase transition mechanisms. However, identification of these mechanisms from observational data sets has remained a problem due to a large variety of possible configurations, many of which are degenerate. Here, an approach based on a combination of deep learning-based semantic segmentation, rotationally invariant variational autoencoder (VAE), and non-negative matrix factorization to enable learning of a latent space representation of the data with multiple real-space rotationally equivalent variants mapped to the same latent space descriptors is introduced. By varying the size of training sub-images in the VAE, the degree of complexity in the structural descriptors is tuned from simple domain wall detection to the identification of switching pathways. Importantly, this yields a powerful tool for the exploration of the dynamic data in mesoscopic electron, scanning probe, optical, and chemical imaging. Moreover, this work adds to the growing body of knowledge of incorporating physical constraints into the machine and deep-learning methods to improve learned descriptors of physical phenomena.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Developing Reaction Chemistry Models from Reactive Molecular Dynamics: TATB

Reactive Molecular Dynamic (RMD) are used to simulate the cook-off chemistry of TATB at a variety of fixed density and fixed temperature conditions. The chemical transformations are monitored using a Coordination Geometry Analysis (CGA) approach which tracks which atom types are bonded to each specific atom. This particularly identifies oxidation state changes that occur during the transformations. Correlations between these different chemical changes are identified using a Non-negative Matrix Factorization (NMF) approach. These identify reduced order chemistry models for the TATB system which contains six components whose concentration profiles are a function of both the temperature and density/pressure. The time histories of these transformations appear to show exponential growth/decay properties that could be fit with Arrhenius rates. These components should form the basis of deflagration rate models for these materials which could then be used in mesoscale simulations to analyze accidental initiation, shock-to-detonation and detonation propagation.

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

Compact representations of structured BFGS matrices

For general large-scale optimization problems compact representations exist in which recursive quasi-Newton update formulas are represented as compact matrix factorizations. For problems in which the objective function contains additional structure, recent structured quasi-Newton methods exploit available second-derivative information and approximate unavailable second derivatives. Here, this article develops the compact representations of two structured Broyden-Fletcher-Goldfarb-Shanno update formulas. The compact representations enable efficient limited memory and initialization strategies. Two limited memory line search algorithms are described for which extensive numerical results demonstrate the efficacy of the algorithms, including comparisons to IPOPT on large machine learning problems, and to L-BFGS on a real world large scale ptychographic imaging application.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

DOME: Directional medical embedding vectors from Electronic Health Records

Motivation: The increasing availability of Electronic Health Record (EHR) systems has created enormous potential for translational research. Recent developments in representation learning techniques have led to effective large-scale representations of EHR concepts along with knowledge graphs that empower downstream EHR studies. However, most existing methods require training with patient-level data, limiting their abilities to expand the training with multi-institutional EHR data. On the other hand, scalable approaches that only require summary-level data do not incorporate temporal dependencies between concepts. Methods: We introduce a DirectiOnal Medical Embedding (DOME) algorithm to encode temporally directional relationships between medical concepts, using summary-level EHR data. Specifically, DOME first aggregates patient-level EHR data into an asymmetric co-occurrence matrix. Then it computes two Positive Pointwise Mutual Information (PPMI) matrices to correspondingly encode the pairwise prior and posterior dependencies between medical concepts. Following that, a joint matrix factorization is performed on the two PPMI matrices, which results in three vectors for each concept: a semantic embedding and two directional context embeddings. They collectively provide a comprehensive depiction of the temporal relationship between EHR concepts. Results: We highlight the advantages and translational potential of DOME through three sets of validation studies. First, DOME consistently improves existing direction-agnostic embedding vectors for disease risk prediction in several diseases, for example achieving a relative gain of 5.5% in the area under the receiver operating characteristic (AUROC) for lung cancer. Second, DOME excels in directional drug-disease relationship inference by successfully differentiating between drug side effects and indications, correspondingly achieving relative AUROC gain over the state-of-the-art methods by 10.8% and 6.6%. Finally, DOME effectively constructs directional knowledge graphs, which distinguish disease risk factors from comorbidities, thereby revealing disease progression trajectories. The source codes are provided at https://github.com/celehs/Directional-EHRembedding.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2025

Speciation mapping of the oxidation layer on aged uranium dioxide using scanning transmission x-ray spectromicroscopy

In this study, UO 2 was aged in humid air and prepared as a thin section using a focused ion beam (FIB) instrument. The specimen was measured using synchrotron radiation spectromicroscopy techniques at the Beamline 11.0.2 STXM end station of the Advanced Light Source (ALS). Non-negative matrix factorization (NMF) methods were used to identify and map three component x-ray absorption near-edge structure (XANES) spectra in the oxygen K-edge data, revealing a surface layer of U 4 O 9 with a thickness of 206 ± 21 nm, and the bulk of the sample remaining as UO 2 . Uranium N 4,5 -edge XANES spectromicroscopy supports these results. The diffusion-controlled parabolic rate constant for UO 2 oxidation to U 4 O 9 was calculated from the observed layer thickness and compared to literature values. Complementary transmission electron microscopy (TEM) was used to image the sample and identify the phases present in various regions, confirming the STXM results.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2022

Analysis of Interpretable Data Representations for 4D-STEM Using Unsupervised Learning

Abstract Understanding the structure of materials is crucial for engineering devices and materials with enhanced performance. Four-dimensional scanning transmission electron microscopy (4D-STEM) is capable of mapping nanometer-scale local crystallographic structure over micron-scale field of views. However, 4D-STEM datasets can contain tens of thousands of images from a wide variety of material structures, making it difficult to automate detection and classification of structures. Traditional automated analysis pipelines for 4D-STEM focus on supervised approaches, which require prior knowledge of the material structure and cannot describe anomalous or deviant structures. In this article, a pipeline for engineering 4D-STEM feature representations for unsupervised clustering using non-negative matrix factorization (NMF) is introduced. Each feature is evaluated using NMF and results are presented for both simulated and experimental data. It is shown that some data representations more reliably identify overlapping grains. Additionally, real space refinement is applied to identify spatially distinct sample regions, allowing for size and shape analysis to be performed. This work lays the foundation for improved analysis of nanoscale structural features in materials that deviate from expected crystallographic arrangement using 4D-STEM.

Bruefach, Alexandra (ORCID:0000000209323477)↗
DOE OSTI2025

Sources and Composition of Organic Aerosols in the Central Arctic during Spring and Summer

Organic aerosols (OA) are a major component of Arctic aerosol mass and influence the region’s radiation budget, yet their sources and physicochemical properties remain largely unknown. We investigate OA sources and climate-relevant characteristics over the central Arctic Ocean in spring and summer, by applying positive matrix factorization to aerosol mass spectrometry data from two ship-based expeditions (2018 and 2020), complemented by total and interstitial aerosol measurements during fog periods. Six distinct OA factors were identified: haze-related OA, Arctic oxygenated OA, two mixed-OA types resembling biomass burning and primary marine OA linked to warm-air intrusions, marine OA, and hydrocarbon-related OA. Seasonal transitions strongly shaped OA composition. Following polar sunrise, highly oxygenated OA, likely formed secondarily from photo-oxidized volatile organic compounds, became dominant. After the polar vortex collapsed in May, episodic spikes in marine OA from the marginal ice zone appeared, alongside a reduced influence from Eurasian anthropogenic sources. These transitions influenced OA oxidation state and related properties including volatility, acidity, and hygroscopicity, highlighting the role of the Arctic spring atmosphere as an active photochemical reactor. Overall, OA was highly oxidized, with particles activated in fog during summer showing even greater oxidation, suggesting that central Arctic OA can be highly cloud active.

Arctic haze↗
DOE OSTI2022

Gamma-Ray-Induced Formation of Uranyl Peroxide Cage Clusters

Aqueous solutions of lithium uranyl triperoxide, Li 4 [UO 2 (O 2 ) 3 ] (LiUT), were irradiated with gamma rays at room temperature and found to form the uranyl peroxide cage cluster, Li 24 [(UO 2 )(O 2 )(OH)] 24 (Li–U 24 ). Raman spectroscopy and 18 O labeling were used to identify the Raman-active vibrations of LiUT. With these assignments, the concentration of LiUT was tracked as a function of radiation dose. A discrepancy between monomer removal and cluster formation suggests that the reaction proceeds by the assembly of an intermediate. Non-negative matrix factorization was used to separate Raman spectra into components and resulted in the identification of a unique intermediate species. Much of the conversion appears to be driven by water radiolysis products, particularly the hydroxyl radical. Furthermore, this differs from the 18 O-labeled copper-catalyzed formation of U 24 , which progresses at a steady rate with no observation of intermediates. Li–U 24 in solution decomposes at high radiation doses resulting in a solid insoluble product similar to Na-compreignacite, Na 2 (UO 2 ) 6 O 4 (OH) 6 ·7H 2 O, which contains uranyl oxyhydroxy sheets.

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

Novel Application of Machine Learning Techniques for Rapid Source Apportionment of Aerosol Mass Spectrometer Datasets

In this work, we apply machine learning approaches sparse multinomial logistic regression to classify aerosol mass spectrometer (AMS) unit mass resolution (UMR) data followed by an ensemble regression technique for source apportionment of organic aerosols (OA). The classifier was trained on 60 well characterized laboratory and positive matrix factorization (PMF) deconvolved reference spectra to identify eight OA types. These include four laboratory-derived secondary organic aerosol (SOA) spectra, which include isoprene photooxidation SOA, isoprene epoxydiols (IEPOX) SOA, a monoterpene SOA type that includes a-pinene and ß-pinene SOA, and aromatic SOA from oxidation of naphthalene and m-xylene precursors, as well as PMF deconvolved spectra for three primary organic aerosol (POA) types, namely, hydrocarbon-like organic aerosol (HOA), biomass burning organic aerosol (BBOA), and cooking OA (COA), and a more oxidized oxygenated OA type (MO-OOA). A 5-fold cross-validation strategy, repeated 10 times, was used to assess the classifier’s performance. The classifier had high classification accuracy for COA, aromatic SOA, and isoprene SOA spectra but incorrectly classified ~9% by number of MO-OOA spectra as BBOA, 12% of BBOA spectra as HOA (and vice versa), and 18% of IEPOX-SOA spectra as aromatic SOA. Next, an ensemble regression model was trained on an artificially generated dataset consisting of mixtures of different OA types to assess its ability to predict fractional mass abundances from classification probabilities of various OA species obtained from the multinomial logistic regression classifier trained on the reference spectra. Ultimately, the proposed approach was applied for source apportionment of aircraft-based AMS measurements of OA UMR spectra during the HI-SCALE field campaign. On two representative days (May 6th and 18th, 2016), the algorithm determined that ~50-60% of OA by mass was MO-OOA, which represented a highly aged organic aerosol mixture from different sources. On both days, BBOA was determined to contribute less than 10% to OA by mass. However, on May 18th, the aromatic SOA fraction was higher compared to that on May 6th. The proposed approach is capable of rapidly analyzing AMS data in real time, making it suitable for applications where rapid source apportionment of AMS OA spectra is desirable.

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

Chemical Robotics Enabled Exploration of Stability in Multicomponent Lead Halide Perovskites via Machine Learning

Metal halide perovskites have attracted immense interest as a promising material for a variety of optoelectronic and sensing applications. However, issues regarding long-term stability have emerged as the key bottleneck for commercialization. Here, we develop an automated experimental workflow based on combinatorial synthesis and rapid throughput characterization to explore long-term stability of these materials in ambient conditions. We apply it to four model perovskite systems: MA x FA y Cs 1–x–y PbBr 3 , MA x FA y Cs 1–x–y PbI 3 , Cs x FA y MA 1–x–y Pb(Br x+y I 1–x–y ) 3 , and Cs x MA y FA 1–x–y Pb(I x+y Br 1–x–y ) 3 . Non-negative matrix factorization and Gaussian process regression are used to interpolate the photoluminescent behavior of the phase diagram. This interpolative regression analysis helps to distinguish mixtures that form solid solutions from those that segregate into multiple materials, pointing out the most stable regions of the phase diagram. Here, we find stability dependence on composition to be nonuniform within the composition space, suggesting the presence of potential preferential compositional regions. This proposed workflow is universal and can be applied to other solution-processable materials.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Application of a Chemical Index to Aerosol Mass Spectrometry: Delta Plots and Functional Group Distributions

A better understanding of the chemical properties of organic aerosol (OA) particles will improve our ability to characterize their sources and predict their lifetime. The high-resolution time-of-flight aerosol mass spectrometer (HR-ToF-AMS) is widely used to measure OA in real time using thermal vaporization followed by electron ionization (EI). EI creates fragment ions that can be assigned to functional groups using delta analysis, a method of classifying mass spectra according to the presence of different chemically related ion series. In this study, we demonstrate the application of delta analysis to characterize molecular structures using a new visualization method. We also use delta analysis to quantify the functional group distribution with an average absolute error of ∼5–6% for individual standard molecules, comparable to the error observed for OA mixtures from biomass and coal combustion fit with Fourier transform infrared spectroscopy. Finally, we apply delta functional group analysis to AMS positive matrix factorization (PMF) factors across seven different field campaigns and find a similar composition across the more oxidized factors with about 55% acid and 26% alcohol groups. The analysis method described here can be applied to any HR-ToF-AMS data set to provide quantitative relative functional group distributions for OA mixtures.

aerosol↗
DOE OSTI2022

Node Distortion as a Tunable Mechanism for Negative Thermal Expansion in Metal–Organic Frameworks

Chemically functionalized series of metal–organic frameworks (MOFs), with subtle differences in local structure but divergent properties, provide a valuable opportunity to explore how local chemistry can be coupled to long-range structure and functionality. Using in situ synchrotron X-ray total scattering, with powder diffraction and pair distribution function (PDF) analysis, we investigate the temperature dependence of the local- and long-range structure of MOFs based on NU-1000, in which Zr 6 O 8 nodes are coordinated by different capping ligands (H 2 O/OH, Cl – ions, formate, acetylacetonate, and hexafluoroacetylacetonate). We show that the local distortion of the Zr 6 nodes depends on the lability of the ligand and contributes to a negative thermal expansion (NTE) of the extended framework. Using multivariate data analyses, involving non-negative matrix factorization (NMF), we demonstrate a new mechanism for NTE: progressive increase in the population of a smaller, distorted node state with increasing temperature leads to global contraction of the framework. The transformation between discrete node states is noncooperative and not ordered within the lattice, i.e., a solid solution of regular and distorted nodes. Density functional theory calculations show that removal of ligands from the node can lead to distortions consistent with the Zr···Zr distances observed in the experiment PDF data. Control of the node distortion imparted by the nonlinker ligand in turn controls the NTE behavior. Furthermore, these results reveal a mechanism to control the dynamic structure of MOFs based on local chemistry.

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

Characterizing Drought Behavior in the Colorado River Basin Using Unsupervised Machine Learning

Drought is a pressing issue for the Colorado River Basin (CRB) due to the social and economic value of water resources in the region and the significant uncertainty of future drought under climate change. Here, we use climate simulations from various Earth System Models (ESMs) to force the Variable Infiltration Capacity hydrologic model and project multiple drought indicators for the sub-watersheds within the CRB. We apply an unsupervised machine learning (ML) based on Non-Negative Matrix Factorization using K-means clustering (NMFk) to synthesize the simulated historical, future, and change in drought indicators. The unsupervised ML approach can identify sub-watersheds where key changes to drought indicator behavior occur, including shifts in snowpack, snowmelt timing, precipitation, and evapotranspiration. While changes in future precipitation vary across ESMs, the results indicate that the Upper CRB will experience increasing evaporative demand and surface-water scarcity, with some locations experiencing a shift from a radiation-limited to a water-limited evaporation regime in the summer. Large shifts in peak runoff are observed in snowmelt-dominant sub-watersheds, with complete disappearance of the snowmelt signal for some sub-watersheds. The work demonstrates the utility of the NMFk algorithm to efficiently identify behavioral changes of drought indicators across space and time and to quickly analyze and interpret hydro climate model results.

54 ENVIRONMENTAL SCIENCES↗