An enterprise control assessment case study of the energy–water nexus for the ISO New England system
Not Available
SEARCH · Search NASA
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.
Not Available
As countries accelerate their energy transitions, understanding how renewable energy (RE) systems structurally integrate into national economies is essential. This study presents a longitudinal economic input-output (EIO) analysis of the renewable energy sector in South Korea from 2016 to 2022. We develop a novel EIO-based framework that disaggregates the RE sector both by energy source (thermal, hydro, nuclear and renewable) and by industrial function (manufacturing, generation, and services), allowing for a detailed assessment of production dynamics, value-added creation, and import dependency. By quantifying backward and forward linkages and induced economic effects, the analysis reveals persistent structural vulnerabilities in renewable manufacturing and increasing sectoral interdependencies. Results reveal that while the renewable energy sector's production and value-added shares have increased, critical segments remain highly import-dependent, particularly in equipment manufacturing. The analysis highlights systemic gaps in domestic supply chain resilience and offers sector-specific insights for reducing vulnerability and enhancing energy security. Although applied to South Korea as a case study, the proposed framework is designed to be transferable to other national contexts where renewable energy planning requires economic structural insights. The findings offer policy-relevant guidance for enhancing domestic energy resilience and aligning industrial strategy with long-term decarbonization goals.
Coastal environments are dynamic interfaces that mediate carbon and nutrient exchanges between terrestrial landscapes and open waters, and understanding the biogeochemical factors controlling these exchanges, particularly iron (Fe) redox transformations, is crucial for predicting coastal ecosystem functions. Here, we investigated the mechanisms controlling Fe speciation changes across upland-to-shoreline gradients in freshwater and estuarine soils using Fe K-edge X-ray absorption spectroscopy, solid and porewater composition analysis, and 16S rRNA sequencing analysis. We show that Fe transformations depend primarily on inundation patterns. In unsaturated uplands, Fe occurs as Fe(III) oxyhydroxides, mainly goethite (9–35 %), Fe(II,III)-phyllosilicates (39–89 %), and Fe(III)-organic species (0–61 %). Soils influenced by estuarine waters exhibit porewater sulfide concentrations reaching up to 221 μM, Fe- and S-cycling bacteria, and up to 81 % pyrite (FeS 2 ), indicating that sulfur-driven redox dynamics control Fe transformations. In lacustrine wetlands, Fe(III) reduction is indicated by porewater Fe(II) concentrations increasing to 1.0–2.1 mM, and ~10–15 % of Fe as Fe(II,III)-(hydr)oxides (green rust), vivianite (Fe 3 (PO 4 ) 2 ·8H 2 O), and/or adsorbed Fe(II) species. EXAFS data also indicate reduction of structural Fe(III) to Fe(II) in phyllosilicates. The presence of Fe- and S-cycling bacteria, as well as sulfide (0–10 μM), suggests that Fe-cycling is microbially driven and potentially coupled with cryptic S-cycling. Fe(II) oxidation was indicated above/near the water table by the presence of Fe(III) oxyhydroxides (ferrihydrite, lepidocrocite). Furthermore, negligible Fe(III) or sulfate reduction was observed at some water-saturated sites located at the upland-wetland transition, likely due to oxic (sub-)surface water inputs. Overall, our results highlight the importance of considering both Fe-speciation and hydro-biogeochemical dynamics when predicting Fe-cycling at coastal interfaces.
We apply thermoporoelasticity and a sequentially coupling technique for modeling thermally-driven coupled Thermo-Hydro-Mechanical (THM) processes in tight claystone. A THM benchmark case with a corresponding analytic solution for thermoporoelasticity under a constant heat loading verifies the model. Thereafter, two in situ heating experiments are simulated for model validation: a smaller-scale heating experiment (TED experiment) and a larger-scale experiment (ALC experiment) in Callovo-Oxfordian (COx) claystone at the Meuse/Haute-Marne underground research laboratory in France. The model exhibits good performance to match the observed temperature and pore pressure evolution for the smaller-scale TED experiment. For the larger-scale ALC experiment, general trends of thermal-pressurization are captured in the modeling, but pressure is underestimated at some monitoring points during cool-down. This indicates that the THM response in the field may be affected by the variability of rock's properties or irreversible or time-dependent mechanical processes that are not included in the current thermoporoelastic model. The main contributions of this work are as follows: (1) we verify and validate the numerical simulator, TOUGH-FLAC, to be a valuable coupled THM modeling tool; (2) prove that the laboratory determined material parameters can be used as reference values for upscaling experiments. However, to better identify and quantify THM processes with modeling of in situ tests, more emphasize should be dedicated to obtaining high-quality mechanical deformation data.
We conduct coupled thermo-hydro-mechanical modeling of a KBS-3V repository design in crystalline rocks, using data and conditions from the Forsmark in Olkiluoto repository sites in Sweden and Finland. The study focuses on repository performance related to the impact of thermal and hydraulic evolution on the potential for thermal–mechanical damage to underground repository excavations. For the designs and conditions considered at the Forsmark and Olkiluoto repository sites, the simulations show a peak temperature well under the adopted performance target of a 100°C maximum temperature, whereas there is still a high potential for thermal–mechanical damage to the KBS-3V waste deposition holes. The thermal–mechanical damage is much more likely if rock permeability is so low that it delays saturation and swelling of bentonite-clay-based backfill beyond the time for the thermal–mechanical peak, which occurs 50 to 100 years after nuclear waste deposition. We also found that sidewalls of the KBS-3V emplacement tunnels are vulnerable to tensile fracturing due to the combined effect of thermal stressing and backfill swelling. The study highlights a strong interaction between bentonite-based backfill and host rock through capillary suction along with induced rock desaturation. A careful design and selection of the bentonite-clay-based backfill materials for KBS-3V tunnels and deposition holes can facilitate a timely saturation and backfill swelling that in turn can minimize thermal–mechanical damage.
In the manufacturing sector, water has been often considered too cheap to conserve. Such thinking relies on water valuations that limit the value of water to the price paid. Using such simple methods, the share of water cost to total manufacturing cost is significantly small, <3%. As a result, conserving water and enabling technology uptake is difficult to justify economically and slow to advance, hindering progress toward sustainable water use. However, the value of water to a manufacturer is far greater than the price paid. Valuations such as the true cost of water consider the additional in-plant treatment and energy costs and have been gaining greater traction in the manufacturing sector. However, true cost alone still undervalues water by not accounting for economic and social costs related to scarcity and environmental externalities. This paper makes the case and presents a framework for valuing manufacturing water beyond the price paid and the true cost. The proposed fuller valuation of manufacturing water takes into account the internal and opportunity costs associated with the realization of water risks. The paper follows with a review of a wide range of water valuation metrics, both at the specific industry level and regional/economy-wide level. The use of various valuation metrics incorporating the relationship between the change in value with change in water use, such as marginal value of water, shadow price, and elasticity at the specific industry level, has been limited in the U.S. manufacturing sector. Further, a limited number of studies exist on data-intensive subjective evaluation techniques such as computable general equilibrium modeling and input-output modeling for regional water valuation. After reviewing water value metrics, several recent case studies from manufacturers from the literature are presented to illustrate both the promise and challenges of a fuller valuation of water as proposed here. Some large multinational corporations have moved toward assessing the value of water via supply chain sustainability initiatives, environmental profit and loss accounting, estimating risk-adjusted values of water, hydro-economic modeling, natural capital asset valuation, and developing value chain indices. This paper provides policymakers and technology developers a framework for monetizing water value beyond its true cost and current metrics. If adopted, such fuller water valuations can help make the business case for the development and deployment of cost-effective water-conserving technologies, thereby improving the sustainability of the manufacturing sector with respect to water.
The large production volumes of commodity polyolefins (specifically, polyethylene, polypropylene, polystyrene, and poly(vinyl chloride)), in conjunction with their low unit values and multitude of short-term uses, have resulted in a significant and pressing waste management challenge. Only a small fraction of these polyolefins is currently mechanically recycled, with the rest being incinerated, accumulating in landfills, or leaking into the natural environment. Since polyolefins are energy-rich materials, there is considerable interest in recouping some of their chemical value while simultaneously motivating more responsible end-of-life management. An emerging strategy is catalytic depolymerization, in which a portion of the C–C bonds in the polyolefin backbone is broken with the assistance of a catalyst and, in some cases, additional small molecule reagents. When the products are small molecules or materials with higher value in their own right, or as chemical feedstocks, the process is called upcycling. This review summarizes recent progress for four major catalytic upcycling strategies: hydrogenolysis, (hydro)cracking, tandem processes involving metathesis, and selective oxidation. Key considerations include macromolecular reaction mechanisms relative to small molecule mechanisms, catalyst design for macromolecular transformations, and the effect of process conditions on product selectivity. Metrics for describing polyolefin upcycling are critically evaluated, and an outlook for future advances is described.
Uncontrolled sand production impedes continuous gas production from a hydrate reservoir as observed in the past field-scale gas production tests. Sand mobilization is strongly linked with sediment deformation and high pressure gradient. Sand production occurs when the mobilized sands reach the well. Throughout gas production from a hydrate reservoir, because of hydrate dissociation, deformation and pressure gradient change both in time and space. Therefore, it is necessary to consider the entire process to identify where sand is likely mobilized and how much mobilized sand could reach the well. Here, this study utilizes a coupled thermal, hydrological, chemical, and mechanical (thermo-hydro-chemo-mechanical) sand migration model to simulate a year long gas production from the hydrate reservoir at the Kuparuk 7–11–12 site in the Prodhoe Bay Unit, Alaska. It is found that sand would mainly come from the lower portion of the production zone where faster hydrate dissociation occurs. The relatively faster hydrate dissociation coupled with the fact that hydrate-bearing sediments and hydrate-free sediments have similar stiffness but different strengths causes complex stress transfer, which results in excessive shear deformation when the upper portion of the production zone starts to dissociate. This is evident not only near the well but also away from the well, leading to a large amount of sand mobilization. Another focus of the modeling study is evaluation of the effect of sand migration on gas production. The comparison of two extreme cases–with sandscreen and without sandscreen–suggests that migrated and settled solids around a sand control device can prevent pressure drop across the near-wellbore zone and lead to reduction of the gas production rate.
Hydraulic fracturing of shale reservoirs resulted in significant opportunity for increased oil and gas production in the United States. Rock-fluid interactions can cause mineral dissolution and precipitation reactions that lead to permeability changes in the shale matrix, which ultimately may affect transport pathways and hydrocarbon production. Understanding the distribution of secondary precipitates, such as barite and Fe(III) (hydro)oxides, and cation leaching at the rock-fluid interface is an important step to further investigate how these geochemical processes can change permeability and transport pathways. In this study, thin sections of the fracture-matrix interface were made from reacted Marcellus shale cores. The thin sections were characterized using synchrotron X-ray fluorescence imaging and synchrotron X-ray absorption spectroscopy. Fe species with different oxidation states were identified in the maps, together with barite and Ca distribution. The results show that ferrihydrite, as newly formed Fe(III)-bearing precipitates, aligned well with the border of the Ca (e.g., calcite) leaching region in the reaction front. Some Fe-containing clay also dissolved, but the dissolution region for the clay was not as deep as the calcite. Further, the reaction front is about three times deeper in the direction parallel to the shale bedding than that perpendicular to the bedding. The Ca leaching region can be an index for reaction front detection for Marcellus shale. Reactive transport modeling was conducted and the predicted Ca leaching boarder align well with ferrihydrite precipitation, consistent with the experimental observation. The carbonate mineral dissolution can be crucial to promote fluid access into the shale matrix. Together with our previous study on the shale reactive surface, this follow-up study showed similar Ca leaching region and Fe(III) precipitates distribution in the reaction front regardless of barite precipitation on the surface, indicating that the barite coatings on the surface may not pose a significant impact on reactive transport at the shale-fluid interface.
This paper presents a coupled techno-economic and life cycle assessment of “blue” hydrogen to be produced at a hydrogen facility through steam methane reforming (SMR) equipped with carbon capture and storage (CCS). Blue hydrogen was modeled in ChemCAD, while an integrated asset model represented the carbon capture and storage chain. An unabated carbon dioxide (CO 2 ) configuration release 11.99 kgCO 2 -eq/ kgH 2 . Capturing ≥95% of the CO 2 stream lowers the carbon footprint to 6.59 kgCO 2 -eq/kgH2 but raises the levelized cost of hydrogen (LCOH) from $\$$1.82/kgH 2 (no CO 2 capture) to $\$$3.22/kgH 2 ; the U.S. 45Q tax credit reduces it to $\$$2.59/kgH 2 . Incorporating CCS reduces the levelized net present value from $\$$0.87/kgH 2 to $\$$0.74/kgH 2 , owing to additional capture, transport, and storage costs. Supplying SMR with low-carbon electricity, especially nuclear, wind, or hydro, delivers the lowest carbon footprint relative to geothermal or grid mixes. Sensitivity analysis identifies that hydrogen sales price, internal rate of return, and CCS cost as the strongest economic levers, while electricity demand dominates residual lifecycle emissions. The results underscore a clear trade-off; substantial CO 2 reductions are achievable, but only with higher production costs, making supportive policy instruments, access to clean power, and robust hydrogen markets essential for large-scale deployment of blue hydrogen.
A new donor–acceptor biradical complex, Tp Cum,Me Zn(SQ-VD) (Tp Cum,Me Zn + = zinc(II) hydro-tris(3-cumenyl-5-methylpyrazolyl)borate complex cation; SQ = orthosemiquinone; VD = oxoverdazyl), which is a ground-state analogue of a charge-separated excited state, has been synthesized and structurally characterized. Here, the magnetic exchange interaction between the S = 1/2 SQ and the S = 1/2 VD within the SQ-VD biradical ligand is observed to be ferromagnetic, with $J_{SQ-VD}$ = +77 cm –1 (H = -2$J_{SQ-VD}\bar{S}_{SQ}·\bar{S}_{VD}$) determined from an analysis of the variable-temperature magnetic susceptibility data. The pairwise biradical exchange interaction in Tp Cum,Me Zn(SQ-VD) can be compared with that of the related donor–acceptor biradical complex Tp Cum,Me Zn(SQ-NN) (NN = nitronyl nitroxide, S = 1/2), where $J_{SQ-NN}$ ≅ +550 cm –1 . This represents a dramatic reduction in the biradical exchange by a factor of ~7, despite the isolobal nature of the VD and NN acceptor radical SOMOs. Computations assessing the magnitude of the exchange were performed using a broken-symmetry density functional theory (DFT) approach. These computations are in good agreement with those computed at the CASSCF NEVPT2 level, which also reveals an S = 1 triplet ground state as observed in the magnetic susceptibility measurements. A combination of electronic absorption spectroscopy and CASSCF computations has been used to elucidate the electronic origin of the large difference in the magnitude of the biradical exchange coupling between Tp Cum,Me Zn(SQ-VD) and Tp Cum,Me Zn(SQ-NN). A Valence Bond Configuration Interaction (VBCI) model was previously employed to highlight the importance of mixing an SQ SOMO → NNLUMO charge transfer configuration into the electronic ground state to facilitate the stabilization of the high-spin triplet (S = 1) ground state in Tp Cum,Me Zn(SQ-NN). Here, CASSCF computations confirm the importance of mixing the pendant radical (e.g., VD, NN) LUMO (VD LUMO and NN LUMO ) with the SOMO of the SQ radical (SQ SOMO ) for stabilizing the triplet, in addition to spin polarization and charge transfer contributions to the exchange. An important electronic structure difference between Tp Cum,Me Zn(SQ-VD) and Tp Cum,Me Zn(SQ-NN), which leads to their different exchange couplings, is the reduced admixture of excited states that promote ferromagnetic exchange into the Tp Cum,Me Zn(SQ-VD) ground state, and the intrinsically weaker mixing between the VD LUMO and the SQ SOMO compared to that observed for Tp Cum,Me Zn(SQ-NN), where this orbital mixing is significant. The results of this comparative study contribute to a greater understanding of biradical exchange interactions, which are important to our understanding of excited-state singlet–triplet energy gaps, electron delocalization, and the generation of electron spin polarization in both the ground and excited states of (bpy)Pt(CAT-radical) complexes.
Selective oxidation of hydrocarbons over metal oxide catalysts frequently invokes Mars-van Krevelen (MvK) mechanisms involving lattice oxygens. At typical reaction temperatures, oxidative dehydrogenation reactions predominate and compete with total oxidation; production of more useful oxygenates is challenging. However, because regeneration of lattice oxygens has been previously hypothesized to involve transient formation of surface hydroperoxides or peroxides, MvK cycles should be capable of forming products derived from both lattice oxygen and surface (hydro)peroxide species, under certain conditions. In this report, we study the low temperature reaction of cyclohexene with O 2 over CeO 2 -supported VO x catalysts. Below 150°C and for higher VO x loading, cyclohexene converts to benzene and trans-1,2-cyclohexanediol with no appreciable CO 2 formation. Control experiments show similar behavior in other supported vanadia catalysts. Conversely, using vaporized H 2 O 2 as oxidant exclusively gives 1,2-cyclohexanedione, suggesting that MvK reoxidation cycles do not directly involve H 2 O 2 . Furthermore, these findings contribute to the current understanding of selective oxidation mechanisms, and demonstrate routes to the formation of industrially useful oxygenates with O 2 or vaporized H 2 O 2 over metal oxide catalysts.
In the search for energy storage materials, metal octahydrotriborates, M(B3H8)n, n=1,2, are promising candidates but their synthesis suffers from residual solvents which tend to interact and greatly alter their decomposition mechanism. Therefore, we studied the thermal conversion of unsolvated Mg(B3H8)2 to BH4 -: as synthesized, and in the presence of MgH2. The conversion of our unsolvated Mg(B3H8)2 starts at ~100°C and yields ~22 wt% of BH4 - along with the formation of (closo-hydro)borates and volatile boranes. This loss of boron (B) is a sign of poor cyclability of the system. However, the addition of MgH2 to unsolvated Mg(B3H8)2 drastically increases the thermal conversion to 85-88wt% of BH4 - while simultaneously decreasing the amounts of B-losses. Our results strongly indicate that the presence of activated MgH2 substantially decreases the formation of (closohydro) borates and provides the necessary H2 for the B3H8-to-BH4 conversion. This is the first report of a metal octahydrotriborate system to selectively convert to BH4 - under moderate conditions of temperature (200°C) in less than 1h, making the MgB3H8-MgH2 system very promising for energy storage applications.
Direct regeneration of spent Li-ion batteries based on the hydrothermal relithiation of cathode materials is a promising next-generation recycling technology. In order to demonstrate the feasibility of this approach at a large scale, we systematically design and optimize the process parameters to minimize both energy and raw material costs. Specifically, the effects of regenerative processing parameters on the composition, structure, and electrochemical performance of the regenerated cathode materials are investigated via systematic characterization and testing. From this analysis, it was found that the raw material costs can be substantially reduced by either replacing the typically employed 4 M LiOH solution by a cost-effective mixture of 0.1 M LiOH and 3.9 M KOH or recycling of the concentrated 4 M LiOH for continuous relithiation processes. Finally, life cycle analysis suggests that this strategy results in reduced energy consumption and greenhouse gas emissions, leading to an increased potential revenue, particularly when compared with hydro- and pyrometallurgical recycling methods.
The commercialization of biorefineries has been met with significant obstacles due to the technical difficulties of handling solid biomass feedstocks, which are highly variable in physical properties and chemical compositions. An understanding of the principles of unit operations along the supply chain is key to the success of the biofuel industry. This study applies a dynamic life-cycle analysis (DLCA) methodology by developing quantitative relationships between the inputs and outputs for each unit operation based on scientific understanding of the causal effects. Using the DLCA, we assessed system sustainability of drop-in fuel production from fast pyrolysis of pine residues followed by hydro-processing. Life-cycle greenhouse gas (GHG) emissions were calculated for 4641 runs involving two key feedstock parameters: moisture content after field drying and particle size of the feed to pyrolyzer. Moisture contents varied between 25% and 35%, while particle sizes varied between 0.5 and 5 mm in these runs. The life-cycle GHG emissions of these cases vary from 22 to 40 g/MJ. Close examination of each unit operation's contribution to the total GHG emissions reveals that low moisture content after field drying produces savings in energy consumption during feedstock transportation and on-site drying. However, although small particle size leads to overall higher fuel yield, it also requires a significant amount of energy for feedstock size reduction and fuel production. The energy penalty outweighs the benefit of increased fuel yields, especially when fine particles (<1 mm) are used for pyrolysis; thus, small particle size leads to increased GHG emissions overall. The results highlight the trade-offs between the energy demand for preprocessing and the conversion yields, which can be addressed with DLCA.
Heterostructures of three-dimensional (3D) halide perovskites are unstable because of facile anion interdiffusion at halide interfaces. Two-dimensional (2D) Ruddlesden–Popper halide perovskites (RPPs) show suppressed and anisotropic ion diffu-sion that could enable stable RPP heterostructures, yet the direct and general growth of lateral RPP heterostructures remains challenging. Here we show that halide miscibility in RPPs decreases with perovskite layer thickness (n), enabling the for-mation of sharp halide lateral heterostructures from n = 1 and 2 RP lead iodide microplates via anion exchange with hydro-gen bromide vapor. In contrast, RPPs with n = 3 form more diffuse lateral heterojunctions, more similar to those in 3D perov-skites. The anion exchange behaviors are further modulated by the spacer and A-site cations in the RPP structures. These new insights, and kinetic studies of the exchange reactions, enable the preparation of lateral heterostructures from various n = 2 RPPs that are more stable against anion interdiffusion and degradation for potential device applications.
Here, transmetallation of [VCl 3 (THF) 3 ] and [TlTp tBu,Me ] afforded [(Tp tBu,Me )VCl 2 ] (1, Tp tBu,Me = hydro-tris(3-tert-butyl-5-methylpyrazol-1-yl)borate), which was reduced with KC 8 to form a $C_{3v}$ symmetric V II complex, [(Tp tBu,Me )VCl] (2). Complex 1 has a high-spin ($\textit{S}$ = 1) ground state and displays rhombic high-frequency and -field electron paramagnetic resonance (HFEPR) spectra, while complex 2 has an $\textit{S}$ = 3/2 4 A 2 ground state observable by conventional EPR spectroscopy. Complex 1 reacts with NaN 3 to form the V V nitride-azide complex [(Tp tBu,Me )V≡N(N 3 )] (3). A likely V III azide intermediate en route to 3, [(Tp tBu,Me )VCl(N 3 )] (4), was isolated by reacting 1 with N 3 SiMe 3 . Complex 4 is thermally stable but reacts with NaN3 to form 3, implying a bis-azide intermediate, [(Tp tBu,Me )V(N 3 ) 2 ] (A), leading to 3. Reduction of 3 with KC 8 furnishes a trinuclear and mixed-valent nitride, [{(Tp tBu,Me )V} 2 ($μ_{4-}$VN 4 )] (5), conforming to a Robin–Day class I description. Complex 5 features a central vanadium ion supported only by bridging nitride ligands. Contrary to 1, complex 2 reacts with NaN 3 to produce an azide-bridged dimer, [{(Tp tBu,Me )V} 2 (1,3-$μ_2$-N 3 ) 2 ] (6), with two antiferromagnetically coupled high-spin V II ions. Complex 5 could be independently produced along with [($κ_2$-Tp tBu,Me ) 2 V] upon photolysis of 6 in arene solvents. The putative {V IV ≡N} intermediate, [(Tp tBu,Me )V≡N] (B), was intercepted by photolyzing 6 in a coordinating solvent, such as tetrahydrofuran (THF), yielding [(Tp tBu,Me )V≡N(THF)] (B-THF). In arene solvents, B-THF expels THF to afford 5 and [($κ_2$-Tp tBu,Me ) 2 V]. A more stable adduct (B-OPPh 3 ) was prepared by reacting B-THF with OPPh 3 . These adducts of B are the first neutral and mononuclear V IV nitride complexes to be isolated.
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.