Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Multiscale systems”

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

Simulation Models for Exploring Magnetic Reconnection

Simulations have played a critical role in the advancement of our knowledge of magnetic reconnection. However, due to the inherently multiscale nature of reconnection, it is impossible to simulate all physics at all scales. For this reason, a wide range of simulation methods have been crafted to study particular aspects and consequences of magnetic reconnection. This article reviews many of these methods, laying out critical assumptions, numerical techniques, and giving examples of scientific results. Plasma models described include magnetohydrodynamics (MHD), Hall MHD, Hybrid, kinetic particle-in-cell (PIC), kinetic Vlasov, Fluid models with embedded PIC, Fluid models with direct feedback from energetic populations, and the Rice Convection Model (RCM).

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Microstructural impacts on ionic conductivity of oxide solid electrolytes from a combined atomistic-mesoscale approach

Abstract Although multiple oxide-based solid electrolyte materials with intrinsically high ionic conductivities have emerged, practical processing and synthesis routes introduce grain boundaries and other interfaces that can perturb primary conduction channels. To directly probe these effects, we demonstrate an efficient and general mesoscopic computational method capable of predicting effective ionic conductivity through a complex polycrystalline oxide-based solid electrolyte microstructure without relying on simplified equivalent circuit description. We parameterize the framework for Li 7- x La 3 Zr 2 O 12 (LLZO) garnet solid electrolyte by combining synthetic microstructures from phase-field simulations with diffusivities from molecular dynamics simulations of ordered and disordered systems. Systematically designed simulations reveal an interdependence between atomistic and mesoscopic microstructural impacts on the effective ionic conductivity of polycrystalline LLZO, quantified by newly defined metrics that characterize the complex ionic transport mechanism. Our results provide fundamental understanding of the physical origins of the reported variability in ionic conductivities based on an extensive analysis of literature data, while simultaneously outlining practical design guidance for achieving desired ionic transport properties based on conditions for which sensitivity to microstructural features is highest. Additional implications of our results are discussed, including a possible connection between ion conduction behavior and dendrite formation.

25 ENERGY STORAGE↗
DOE OSTI2024

The Role of Radiation-Induced Non-Equilibrium Plutonium Oxidation States in Solution

Plutonium plays a key role in global actinide research and nuclear fuel cycle technologies, and yet, our fundamental understanding of its inherent radiation-induced chemical behavior is limited. These radiation-induced processes cannot simply be switched off, as they are as fundamentally inherent to plutonium as the impact of relativistic effects on its f-electrons. In less chemically complex actinide systems, such as aqueous solutions of neptunium and americium, both steady-state and transient radiolysis products play a significant role in the redox cycling of their oxidation states, which ultimately impacts their chemistry and transport. However, plutonium's multiple, coexisting, and chemically active oxidation states, which comprise of bare ions and dioxo cations, provide additional redox pathways that complicate the aforementioned radiation-induced competition kinetics. Mechanistically understanding the response of these multiple oxidation states, both equilibrium and non-equilibrium, to intense ionizing radiation fields is essential for predicting the behavior of plutonium under a host of conditions that support technological innovation in global nuclear energy and non-proliferation efforts. Here, advances in radiation-induced plutonium solution chemistry will be presented, including new results from time-resolved electron pulse and steady-state gamma and alpha irradiations.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗
DOE OSTI2026

In situ catalyst activation and regeneration enable energy-efficient high-current CO 2 reduction to ethanol-rich C 2+ mixtures

Electrochemical conversion of dissolved CO 2 in bicarbonate electrolytes, i.e., bicarbonate electrolysis, offers distinct advantages over gas diffusion electrode systems by enabling direct utilization of the CO 2 capture electrolyte while bypassing the energy-intensive CO 2 release step. However, bicarbonate electrolysis faces challenges such as CO 2 mass-transfer limitation, local pH-driven CO 2 depletion, and high cathodic potentials. The higher potential often causes catalyst surface reorganization, leading to a gradual loss of active sites and variations in selectivity during CO 2 reduction. Here, we report a directed, in situ activation and regeneration method that allows precatalysts to equilibrate under dynamic (pulsed) electrolysis conditions. We demonstrate in situ activation of a scalable Cu 2 O/Cu mesh that, under short-width (t = 4 s) pulsed electrolysis, provides stable mixed oxidation states of Cu, favoring the formation of an ethanol-rich crude mixture. The pulsed electrolysis waveform, consisting of six distinct segments, is tuned to form Cu + oxides, which are then reduced to generate local alkaline conditions favoring C–C coupling. This synergistic effect results in FEs of 73% for C2+ products and 39% for ethanol at an applied current density of −150 mA cm −2 and a cathodic potential of −1.45 V (vs. RHE). The overall half-cell energy efficiency is ∼30% for C 2+ products. The in situ Raman experiments confirm the role of pCO 2 R in dynamically regenerating Cu+-containing surface species during pulsed operation, thereby steering selectivity towards C 2+ products. A comprehensive multiscale, multiphysics model is developed to investigate the dynamic behavior of copper surface species (Cu, Cu + , and Cu 2+ ) and local microenvironmental conditions during the pCO 2 R. The results reveal that the coexistence of different copper oxidation states, especially the Cu+ intermediate, is critical in steering selectivity towards multicarbon (C 2+ ) products. The dynamic modulation of surface redox states via tailored pulsing strategies favors C–C coupling pathways by inducing localized alkaline conditions and stabilizing reactive intermediates. This work establishes a predictive modeling platform that links pulse waveform design with mechanistic insights into catalyst state evolution and product selectivity. Overall, this study provides valuable insights into the synergistic effect of in situ activation of pre-catalysts and pulsed electrolysis for higher selectivity towards C 2+ products.

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

Multi-scale model based design of membrane reactor/separator processes for intensified hydrogen production through the water gas shift reaction

This paper aims to first quantify the impact of various diffusion models (Maxwell-Stefan, Wilke, Dusty-Gas) on the predictions of a multi-scale membrane reactor/separator mathematical model, and to then demonstrate this model's use for the design and process intensification of membrane reactor/separator systems for hydrogen production. This multi-scale model captures velocity, temperature and species' concentration profiles along the catalyst pellet's radial direction, and along the reactor's axial direction, by solving the momentum, energy, and species transport equations, accounting for convection, conduction, reaction, and diffusion mechanisms. In the first part of work, the effect of pellet-scale design parameters (mean pore diameter, volumetric porosity, tortuosity factor, etc.) and various species' flux models on the model predictions is studied. In the second part, the study focuses on the comparison, in terms of their process intensification characteristics, of various hydrogen production processes. These include a conventional high-temperature shift reactor (HTSR)/low-temperature shift reactor (LTSR) sequence, a novel HTSR/membrane separator (MS)/LTSR/MS sequence, and a process that involves low-temperature shift membrane reactors-LTSMR in a series.

08 HYDROGEN↗
DOE OSTI2024

Solution of the Schrödinger equation for quasi-one-dimensional materials using helical waves

We formulate and implement a spectral method for solving the Schrödinger equation, as it applies to quasi-one-dimensional materials and structures. This allows for computation of the electronic structure of important technological materials such as nanotubes (of arbitrary chirality), nanowires, nanoribbons, chiral nanoassemblies, nanosprings and nanocoils, in an accurate, efficient and systematic manner. Our work is motivated by the observation that one of the most successful methods for carrying out electronic structure calculations of bulk/crystalline systems — the plane-wave method — is a spectral method based on eigenfunction expansion. Our scheme avoids computationally onerous approximations involving periodic supercells often employed in conventional plane-wave calculations of quasi-one-dimensional materials, and also overcomes several limitations of other discretization strategies, e.g., those based on finite differences and atomic orbitals. The basis functions in our method — called helical waves (or twisted waves) — are eigenfunctions of the Laplacian with symmetry adapted boundary conditions, and are expressible in terms of plane waves and Bessel functions in helical coordinates. We describe the setup of fast transforms to carry out discretization of the governing equations using our basis set, and the use of matrix-free iterative diagonalization to obtain the electronic eigenstates. Miscellaneous computational details, including the choice of eigensolvers, use of a preconditioning scheme, evaluation of oscillatory radial integrals and the imposition of a kinetic energy cutoff are discussed. We have implemented these strategies into a computational package called HelicES (Helical Electronic Structure). We demonstrate the utility of our method in carrying out systematic electronic structure calculations of various quasi-one-dimensional materials through numerous examples involving nanotubes, nanoribbons and nanowires. We also explore the convergence properties of our method, and assess its accuracy and computational efficiency by comparison against reference finite difference, transfer matrix method and plane-wave results. We anticipate that our method will find applications in computational nanomechanics and multiscale modeling, for carrying out transport calculations of interest to the field of semiconductor devices, and for the discovery of novel chiral phases of matter that are of relevance to the burgeoning quantum hardware industry.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

A Stochastic Reduced-Order Model for Statistical Microstructure Descriptors Evolution

Integrated computational materials engineering (ICME) models have been a crucial building block for modern materials development, relieving heavy reliance on experiments and significantly accelerating the materials design process. However, ICME models are also computationally expensive, particularly with respect to time integration for dynamics, which hinders the ability to study statistical ensembles and thermodynamic properties of large systems for long time scales. To alleviate the computational bottleneck, we propose to model the evolution of statistical microstructure descriptors as a continuous-time stochastic process using a non-linear Langevin equation, where the probability density function (PDF) of the statistical microstructure descriptors, which are also the quantities of interests (QoIs), is modeled by the Fokker–Planck equation. In this work, we discuss how to calibrate the drift and diffusion terms of the Fokker–Planck equation from the theoretical and computational perspectives. The calibrated Fokker–Planck equation can be used as a stochastic reduced-order model to simulate the microstructure evolution of statistical microstructure descriptors PDF. Considering statistical microstructure descriptors in the microstructure evolution as QoIs, we demonstrate our proposed methodology in three integrated computational materials engineering (ICME) models: kinetic Monte Carlo, phase field, and molecular dynamics simulations.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Revealing 3D Morphological and Chemical Evolution Mechanisms of Metals in Molten Salt by Multimodal Microscopy

Growing interest in molten salts as effective high-temperature heat-transfer fluids for sustainable energy systems drives a critical need to fundamentally understand the interactions between metals and molten salts. This work utilizes the multimodal microscopy methods of synchrotron X-ray nanotomography and electron microscopy to investigate the 3D morphological and chemical evolution of two-model systems, pure nickel metal and Ni-20Cr binary alloy, in a representative molten salt (KCl-MgCl2 50–50 mol %, 800 °C). In both systems, unexpected shell-like structures formed because of the presence of more noble tungsten, suggesting a potential route of using Ni–W alloys for enhanced molten-salt corrosion resistance. The binary alloy Ni-20Cr developed a bicontinuous porous structure, reassembling functional porous metals manufactured by dealloying. This work elucidates better mechanistic understanding of corrosion in molten salts, which can contribute to the design of more reliable alloys for molten salt applications including next-generation nuclear and solar power plants and opens the possibility of using molten salts to fabricate functional porous materials.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Interfacial Electromechanics Predicts Phase Behavior of 2D Hybrid Halide Perovskites

Quasi-two-dimensional (2D) mixed-cation hybrid halide perovskites (A' 2 A N–1 M N X 3N+1 ; A' = large organic molecule with cationic group, A = [Cs + , CH 3 NH 3 + , HC(NH 2 ) 2 + ], M = [Pb, Sn, Ge], X = [I – , Br – , Cl – ]) have rapidly emerged as candidates to improve the structural stability and device lifetime of 3D perovskite semiconductor devices under operating conditions. The addition of the large A' cation to the traditional AMX 3 structure introduces several synthetic degrees of freedom and breaks M–X bonds, giving rise to peculiar critical phase behavior in the phase space of these complex materials. In this work, we propose a thermodynamic model parametrized by first-principles calculations to generate the phase diagram of 2D and quasi-2D perovskites (q-2DPKs) based on the mechanics and electrostatics of the interface between the A' cations and the metal halide octahedral network. Focusing on the most commonly studied methylammonium lead iodide system where A' is n-butylammonium (BA; CH 3 (CH 2 ) 3 NH 3 + ), we find that the apparent difficulty in synthesizing phase-pure samples with a stoichiometric index N > 5 can be attributed to the energetic competition between repulsion of opposing interfacial dipole layers and mechanical relaxation induced by interfacial stress. Our model shows quantitative agreement with experimental observations of the maximum phase-pure stoichiometric index (N crit ) and explains the nonmonotonic evolution of the lattice parameters with increasing stoichiometric index (N). This model is generalizable to the entire family of q-2DPKs and can guide the design of photovoltaic and optical materials that combine the structural stability of the q-2DPKs while retaining the charge carrier properties of their 3D counterparts.

14 SOLAR ENERGY↗
DOE OSTI2026

Detecting and Characterizing Fracture Zones Using a Convolutional Neural Network

This project directly supports the Geothermal Technologies Office (GTO) objectives outlined in the Multi-Year Program Plan (MYPP) by advancing two key research areas: “Exploration and Characterization” and “Data, Modeling, and Analysis.” This project has successfully demonstrated a pre-drilling ability to image and characterize the distribution and connectivity of subsurface faults and fractures, key parameters for identifying permeable pathways that enable geothermal fluids to circulate and produce energy. Specifically, we developed and implemented innovative machine learning methodologies to enhance geothermal exploration. Large-scale faults were detected using a Convolutional Neural Network (CNN), while small-scale fractures were characterized using a novel Double-Beam Neural Network (DBNN). These tools have proven both technically effective and cost-efficient by reducing reliance on expensive exploratory drilling. Through collaboration with our geothermal industry partner, this research has significantly advanced techniques for identifying hidden geothermal systems and extending the productive lifespan of existing geothermal fields. We applied our methods to two geothermal fields—Soda Lake (Nevada) and Lightning Dock (New Mexico)—to identify shallow steam-charged fracture zones and characterize deep faults at depths of 1.5-2 km. The steam zone identified at the Soda Lake geothermal field showed excellent agreement with prior drilling data, validating the effectiveness of our approaches. In addition, the analysis revealed three new prospective drilling targets for further development and verification. The outcomes of this project improve our scientific understanding of geothermal reservoir behavior, enhance exploration efficiency, extend the economic life of existing geothermal plants. Ultimately, these advancements contribute to GTO’s goal of achieving more sustainable, affordable, and data-driven geothermal energy development across the United States.

15 GEOTHERMAL ENERGY↗
DOE OSTI2021

The Colorado East River Community Observatory Data Collection

Abstract The U.S. Department of Energy's (DOE) Colorado East River Community Observatory (ER) in the Upper Colorado River Basin was established in 2015 as a representative mountainous, snow‐dominated watershed to study hydrobiogeochemical responses to hydrological perturbations in headwater systems. The ER is characterized by steep elevation, geologic, hydrologic and vegetation gradients along floodplain, montane, subalpine, and alpine life zones, which makes it an ideal location for researchers to understand how different mountain subsystems contribute to overall watershed behaviour. The ER has both long‐term and spatially‐extensive observations and experimental campaigns carried out by the Watershed Function Scientific Focus Area (SFA), led by Lawrence Berkeley National Laboratory, and researchers from over 30 organizations who conduct cross‐disciplinary process‐based investigations and modelling of watershed behaviour. The heterogeneous data generated at the ER include hydrological, genomic, biogeochemical, climate, vegetation, geological, and remote sensing data, which combined with model inputs and outputs comprise a collection of datasets and value‐added products within a mountainous watershed that span multiple spatiotemporal scales, compartments, and life zones. Within 5 years of collection, these datasets have revealed insights into numerous aspects of watershed function such as factors influencing snow accumulation and melt timing, water balance partitioning, and impacts of floodplain biogeochemistry and hillslope ecohydrology on riverine geochemical exports. Data generated by the SFA are managed and curated through its Data Management Framework. The SFA has an open data policy, and over 70 ER datasets are publicly available through relevant data repositories. A public interactive map of data collection sites run by the SFA is available to inform the broader community about SFA field activities. Here, we describe the ER and the SFA measurement network, present the public data collection generated by the SFA and partner institutions, and highlight the value of collecting multidisciplinary multiscale measurements in representative catchment observatories.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Machine-learning-assisted deciphering of microstructural effects on ionic transport in composite materials: A case study of Li 7 La 3 Zr 2 O 12 -LiCoO 2

The effective diffusivity of ionic species in multiphase materials is critical for the design and function of composite materials for electrochemical energy storage. In practice, effective diffusivity depends sensitively not only on the intrinsic diffusivities of constituting materials but also on their topological arrangement; nevertheless, these coupled contributions are oversimplified in most analytical models. Here, we combine atomistically informed mesoscale modeling and machine learning (ML) analysis to unravel how such features affect effective diffusivity in two-phase composites. Using the Li 7 La 3 Zr 2 O 12 -LiCoO 2 composite solid-state battery cathode as a model system, we compute effective diffusivity for 600 distinct dense polycrystalline microstructures with different topological configurations of grains, grain boundaries, and heterointerfaces. We verify that in addition to atomic-scale variabilities, microstructural feature diversity can significantly impact effective transport properties. Across the ensemble of test microstructures, this often results in bimodal distributions of effective diffusivity that encompass two qualitatively distinct operating mechanisms, which we identify via flux analysis. An ML approach reveals that the most critical determining factors for effective diffusivity are the connectivity of bulk phases and their heterointerfaces. The role of ionic mobility at the heterointerfaces is also discussed. These insights highlight the combined importance of microstructure and interface engineering in tuning the transport properties of ionic species in composite materials. In conclusion, our framework can also be extended for understanding generic microstructure-property relationships in other complex multiphase materials.

25 ENERGY STORAGE↗
DOE OSTI2024

D2NO: Efficient handling of heterogeneous input function spaces with distributed deep neural operators

Neural operators have been applied in various scientific fields, such as solving parametric partial differential equations, dynamical systems with control, and inverse problems. However, challenges arise when dealing with input functions that exhibit heterogeneous properties, requiring multiple sensors to handle functions with minimal regularity. To address this issue, discretization-invariant neural operators have been used, allowing the sampling of diverse input functions with different sensor locations. However, existing frameworks still require an equal number of sensors for all functions. We propose a novel distributed approach to further relax the discretization requirements and solve the heterogeneous dataset challenges. Our method involves partitioning the input function space and processing individual input functions using independent and separate neural networks. A centralized neural network is used to handle shared information across all output functions. This distributed methodology reduces the number of gradient descent back-propagation steps, improving efficiency while maintaining accuracy. Here, we demonstrate that the corresponding neural network is a universal approximator of continuous nonlinear operators and present three numerical examples to validate its performance.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

The Observed and Projected Changes of Global Monsoons: Current Status and Future Perspectives

The global monsoon system, encompassing the Asian-Australian, African, and American monsoons, sustains two-thirds of the world’s population by regulating water resources and agriculture. Monsoon anomalies pose severe risks, including floods and droughts. Recent research associated with the implementation of the Global Monsoons Model Intercomparison Project under the umbrella of CMIP6 has advanced our understanding of its historical variability and driving mechanisms. Observational data reveal a 20th-century shift: increased rainfall pre-1950s, followed by aridification and partial recovery post-1980s, driven by both internal variability (e.g., Atlantic Multidecadal Oscillation) and external forcings (greenhouse gases, aerosols), while ENSO drives interannual variability through ocean-atmosphere interactions. Future projections under greenhouse forcing suggest long-term monsoon intensification, though regional disparities and model uncertainties persist. Models indicate robust trends but struggle to quantify extremes, where thermodynamic effects (warming-induced moisture rise) uniformly boost heavy rainfall, while dynamical shifts (circulation changes) create spatial heterogeneity. Volcanic eruptions and proposed solar radiation modification (SRM) further complicate predictions: tropical eruptions suppress monsoons, whereas high-latitude events alter cross-equatorial flows, highlighting unresolved feedbacks. The emergent constraint approach is booming in terms of correcting future projections and reducing uncertainty with respect to the global monsoons. Critical challenges remain. Model biases and sparse 20th-century observational data hinder accurate attribution. The interplay between natural variability and anthropogenic forcings, along with nonlinear extreme precipitation risks under warming, demands deeper mechanistic insights. Additionally, SRM’s regional impacts and hemispheric monsoon interactions require systematic evaluation. Addressing these gaps necessitates enhanced observational networks, refined climate models, and interdisciplinary efforts to disentangle multiscale drivers, ultimately improving resilience strategies for monsoon-dependent regions.

climate extreme events↗
DOE OSTI2025

Tethys Water Demand Data

U.S. water demand varies sharply by sector and region as land use, population, weather patterns, and economic activity co-evolve. High-resolution water demand data is required to capture these dynamics, support integrated energy-water-land modeling, and local-to-regional water scarcity assessments. This dataset contains gridded (1/8 degree), monthly, multi-sector water demand dataset for the contiguous United States (CONUS) covering 1980-2100 across eight future scenarios of human-Earth system change. The dataset covers irrigation, thermoelectric, municipal (public-supply and domestic), livestock, manufacturing, and mining demands, separately for withdrawals and consumption, and includes per-cell renewable vs. non-renewable water source attributions. The dataset is validated against the latest USGS 2010-2020 water-use data for the three largest water demand sectors (Domestic, Electricity, and Irrigation), with correlations ranging from 0.73-0.95 at the HUC6 scale. The two datasets largely agree on an aggregate basis with per-sector bias falling within +/-7%, but they disagree on the spatial allocation of water with individual HUC6 basins having normalized RMSE from 68-171% and median absolute percent difference from 37-86%. This dataset advances prior global products by combining state-resolved sectoral demands from GCAM-USA, future power-plant siting from the CERF model, and scenario-consistent high-resolution climate and population forcing data across the eight scenarios.

GCAM-USA↗
DOE OSTI2025

IM3 Projected U.S. Western Interconnection Grid Stress Dataset

This dataset provides projected grid stress and reliability results (including all model inputs and outputs from GO WEST and TEP) for Integrated Multisector, Multiscale Modeling (IM3) Phase 2 simulations across eight different scenarios for the U.S. Western Interconnection through 2055. The scenarios include combinations of two Shared Socioeconomic Pathways (SSP3 and SSP5) with four high-resolution climate projections specific to the United States from a set of Thermodynamic Global Warming (TGW) simulations. These climate projections include "hotter" and "cooler" variants for two Representative Concentration Pathways (RCP4.5 and RCP8.5). The resulting eight simulations are: rcp45cooler_ssp3 rcp45cooler_ssp5 rcp45hotter_ssp3 rcp45hotter_ssp5 rcp85cooler_ssp3 rcp85cooler_ssp5 rcp85hotter_ssp3 rcp85hotter_ssp5 GO WEST is an open-source power grid modeling framework for the U.S. Western Interconnection, which allows users to tailor the model depending on their research study and science questions. It covers 28 balancing authorities (BAs) and 12 states in U.S. Western Interconnection. GO WEST allows users to select different number of nodes and come up with a simplified network by utilizing 10,000 nodal topology of the U.S. Western Interconnection (ACTIVSg10k). Users can select different number of nodes, mathematical formulations (linear programming vs. mixed-integer linear programming), transmission line limit scaling factors, and hurdle rate scaling factors. GO WEST offers a unit commitment and economic dispatch (UC/ED) module to simulate grid operations on an hourly scale. In this sense, users can calibrate and validate their model versions by comparing model outputs to historical datasets. TEP is an open-source transmission capacity expansion model, built on the GO WEST framework. It utilizes linear programming to optimize transmission capacity addition investment on existing lines within the GO WEST framework. The TEP model only increases the thermal capacity of existing transmission lines and does not add new lines to the system, which leaves the topology preserved. In order to use TEP model, users need to create scenarios with the GO WEST framework. Please refer to README file for a detailed description of the dataset including individual files and references.

Capacity Expansion Model↗
DOE OSTI2024

A network approach for multiscale catchment classification using traits

Abstract. The classification of river catchments into groups with similar biophysical characteristics is useful to understand and predict their hydrological behavior. The increasing availability of remote sensing and other large-scale geospatial datasets has enabled the use of advanced data-driven approaches to classify catchments using traits such as topography, geology, climate, land cover, land use, and human influence. Unsupervised clustering algorithms based on the Euclidean distance are commonly used for trait-based classification but are not suitable for highly dimensional data. In this study we present a new network-based method for multi-scale catchment classification, which can be applied to large datasets and used to determine the traits associated with different catchment groups. In this framework, two networks are analyzed in parallel: the first being where the nodes are traits and the second being where the nodes are catchments. In both cases, edges represent pairwise similarity, and a network cluster detection algorithm is used for the classification. The trait network is used to investigate redundancy in the trait data and to condense this information into a small number of interpretable categories. The catchments network is used to classify the catchments into clusters and to identify representative catchments for the different groups using the degree centrality metric. We apply this method to classify 9067 river catchments across the contiguous United States at both regional and continental scales using 274 non-categorical traits. At the continental scale, we identify 25 interpretable trait categories and 34 catchment clusters of sizes greater than 50. We find that catchments with similar trait categories are typically located in the same region, with different spatial patterns emerging among clusters dominated by natural and anthropogenic traits. We also find that the catchment clusters exhibit distinct hydrological behavior based on an analysis of streamflow indices. This network approach provides several advantages over traditional means of classification, including better separation of clusters, the use of alternate similarity metrics that are more suitable for highly dimensional data, and reducing redundancy in the trait information. The paired catchment–trait networks enable analysis of hydrological behavior using the dominant trait categories for each catchment cluster. The approach can be used at multiple spatial scales since the network topologies adjust automatically to reflect the trait patterns at the scale of investigation. Finally, the representative catchments identified as hub nodes in the network can be used to guide transferable observational and modeling strategies. The method is broadly applicable beyond hydrology for classification of other complex systems that utilize different types of trait datasets.

54 ENVIRONMENTAL SCIENCES↗