Search NASA⌕ Search

SEARCH · Search NASA

Results for “Connection length”

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 289 records · Page 16

Partial connectomes of labeled dopaminergic circuits reveal non-synaptic communication and axonal remodeling after exposure to cocaine

Dopaminergic (DA) neurons exert profound influences on behavior including addiction. However, how DA axons communicate with target neurons and how those communications change with drug exposure remains poorly understood. We leverage cell type-specific labeling with large volume serial electron microscopy to detail DA connections in the nucleus accumbens (NAc) of the mouse (Mus musculus) before and after exposure to cocaine. We find that individual DA axons contain different varicosity types based on their vesicle contents. Spatially ordering along individual axons further suggests that varicosity types are non-randomly organized. DA axon varicosities rarely make specific synapses (<2%, 6/410), but instead are more likely to form spinule-like structures (15%, 61/410) with neighboring neurons. Days after a brief exposure to cocaine, DA axons were extensively branched relative to controls, formed blind-ended ‘bulbs’ filled with mitochondria, and were surrounded by elaborated glia. Finally, mitochondrial lengths increased by ~2.2 times relative to control only in DA axons and NAc spiny dendrites after cocaine exposure. We conclude that DA axonal transmission is unlikely to be mediated via classical synapses in the NAc and that the major locus of anatomical plasticity of DA circuits after exposure to cocaine are large-scale axonal re-arrangements with correlated changes in mitochondria.

59 BASIC BIOLOGICAL SCIENCES↗

Neural-network quantum states for the nuclear many-body problem

A long-standing goal of nuclear theory is to explain how the structure and dynamics of atomic nuclei and neutron-star matter emerge from the underlying interactions among protons and neutrons. Achieving this goal requires solving the nuclear quantum many-body problem with high accuracy across a wide range of length scales and density regimes. In this review, we discuss how artificial neural network representations of the nuclear many-body wave function have significantly extended the capabilities of continuum quantum Monte Carlo methods. In particular, neural network quantum states enable calculations of larger systems than were previously accessible and provide a flexible framework for capturing phenomena that challenge conventional approaches, including the emergence of nuclear clusters and superfluid phases in dense matter. We highlight recent applications to finite nuclei, infinite nuclear and neutron matter, and dynamical processes relevant to lepton-nucleus and nucleus-nucleus scattering. We also discuss conceptual and methodological connections with condensed matter physics, emphasizing developments in neural network quantum states that bridge strongly correlated systems across disciplines. Together, these developments demonstrate how neural-network methods open new avenues toward unified and accurate descriptions of nuclear structure, matter, and reactions.

Lovato, Alessandro [Argonne; TIFPA-INFN, Trento; V↗

Vibrational spectroscopy analysis of silica and silicate glass networks

Vibrational spectroscopy has been widely used to investigate various structural aspects of the glass network, and there are a plethora of papers reporting subtle but consistent changes in infrared and Raman spectral features of glass upon alterations of glass compositions, thermal histories, mechanical stresses, or surface treatments. However, interpretations of such spectral features are still obscured due to the lack of well-established physical principles accurately describing vibrational modes of the non-crystalline glass network. Due to the non-equilibrium nature of the glass network, three-dimensionally connected without any long-range orders, vibrational spectral features of glass cannot be interpreted using the analogy with those of isolated molecular moieties or crystalline counterparts. Here, this feature article explains why such comparisons are outdated and describes the recent advances made from theoretical calculations of vibrational spectral features of amorphous networks or comparisons of computational results with experimental data. For the interpretation of vibrational spectral features of silica and silicate glasses, the following empirical relationships are suggested: (i) the intensity-weighted peak position of the Si-O-Si stretch mode negatively correlates with the weighted average of the Si-O bond length distribution, and (ii) the broad band of the Si-O-Si bending mode negatively correlates with the Si-O-Si bond angle distribution. Selected examples of vibrational spectroscopic imaging of surface defects are discussed to deliberate the implication of these findings in the structure-property relationship of silica and silicate glass materials. Unanswered questions and continuing research challenges are identified.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Computing diffraction patterns of microstructures from phase-field simulations

The diffraction pattern of a material contains information not only on the crystal structures of its constituting phases, but also on its mesoscale spatial distributions of phases, grains, and ferroelastic, ferroelectric, and ferromagnetic domains. While diffraction patterns from experiments such as X-ray diffraction are presented in the reciprocal or Fourier space, mesoscale microstructure models such as the phase-field method naturally produce real-space images of spatial distribution of chemical composition, structural, and ferroic domains. Although one could rather readily compute the Fourier amplitudes of chemical and structural domain distributions generated by mesoscale simulations, they only contain information about the length scale and alignment of the real-space chemical and structure domains. Therefore, a direct comparison between diffraction experiments and mesoscale microstructure simulations is not possible. Here, we develop a theoretical approach to directly compute the crystal diffraction patterns of microstructures predicted by phase-field simulations. In particular, we consider five representative examples of microstructure patterns involving purely compositional domains, a single pair of tetragonal twin structures, multiple twin variants in a hexagonal system, ferroelectric polar vortices, and polycrystalline grains. The results are compared with previous experimental observations as well as X-ray diffraction experiments performed in the present study. Furthermore, the theoretical framework allows one to directly connect material microstructures and diffraction patterns predicted from phase-field simulations and the corresponding diffraction patterns from experiments, and thus providing guidance to experimental diffraction characterization and interpretation of microstructures.

36 MATERIALS SCIENCE↗

Comparing three generations of D-Wave quantum annealers for minor embedded combinatorial optimization problems

Abstract Quantum annealing (QA) is a novel type of analog computation that aims to use quantum mechanical fluctuations to search for optimal solutions of Ising problems. QA in the transverse Ising model, implemented on D-Wave quantum processing units, are available as cloud computing resources. In this study we report concise benchmarks across three generations of D-Wave quantum annealers, consisting of four different devices, for the NP-hard discrete combinatorial optimization problems unweighted maximum clique and unweighted maximum cut on random graphs. The Ising, or equivalently quadratic unconstrained binary optimization, formulation of these problems do not require auxiliary variables for order reduction, and their overall structure and weights are not highly variable, which makes these problems simple test cases to understand the sampling capability of current D-Wave quantum annealers. All-to-all minor embeddings of size 52, with relatively uniform chain lengths, are used for a direct comparison across the Chimera, Pegasus, and Zephyr device topologies. A grid-search over annealing times and the minor embedding chain strengths is performed in order to determine the level of reasonable performance for each device and problem type. Experiment metrics that are reported are approximation ratios for non-broken chain samples, chain break proportions, and time-to-solution for the maximum clique problem instances. How fairly the quantum annealers sample optimal maximum cliques, for instances which contain multiple maximum cliques, is quantified using entropy of the measured ground state distributions. The newest generation of quantum annealing hardware, which has a Zephyr hardware connectivity, performed the best overall with respect to approximation ratios and chain break frequencies.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Deep Neural Network Algorithm for CMC Microstructure Characterization and Variability Quantification

Microstructure characterization and variability quantification are crucial for understanding ceramic matrix composites (CMCs) mechanical behavior and deformation mechanisms across length scales. Traditionally, analyses of the micrographs obtained from microscopy are labor-intensive. However, with the vast improvement in computer vision (CV) and deep learning (DL), an automated algorithm can be designed to extract essential microstructure variability from micrographs which can then be used to construct a statistically representative volume element (SRVE). The DL-based algorithm spans the taxonomy of microstructure analyses, including semantic segmentation of microstructure constituents, secondary phases, matrix/fiber interface, and defects, and quantifying the microstructure variability in terms of probability distributions. In this work, C/SiNC and SiC/SiNC CMCs microstructures are semantically segmented through a deep convolutional neural network, followed by variability quantification through the implementation of a fully connected regression layer, hence forming a deep regression network. The deep regression network operates in a feedforward regime, in which the neuron output signal traverses through the network in a unidirectional manner. The weight tensor associated with each layer is updated through a backpropagation stochastic gradient descent approach. The input gray-scale image obtained through in-house scanning electron microscope and confocal microscope micrographs is augmented through affine transformations to increase the training set size, which is then processed through four strided convolutional layers. This compresses the image resolution by half at each layer while increasing the image depth by applying different filters (image encoding). The class activation maps (CAMs) corresponding to the applied filters highlight the key architectural features and assist with the semantic segmentation of the microstructure.

Hamza, Mohamed H.↗

Fencing affects movement patterns of two large carnivores in Southern Africa

The erection of peripheral fences around protected areas has been shown to adversely affect the movement behavior, genetic connectivity and long-term viability of animal populations. However, fencing is also used extensively to reduce human encroachment, limit poaching activities, and mitigate human–wildlife conflicts. Limited information on the movement behavior of carnivores is available along the periphery of protected areas, particularly in semi-arid environments. To quantify the potential effects of the 824 km perimeter fence surrounding the Etosha National Park (Etosha, Namibia) on the movement and space use of carnivores, we collected GPS movement data from 36 lions (Panthera leo) and seven spotted hyenas (Crocuta crocuta) over a 14-year period. For these two apex carnivores, we measured the impact of the Etosha fence on their movements (step lengths and path straightness), calculated the likelihood of fence crossings, and evaluated the spatial–temporal characteristics of these crossings. Our results indicated that lions and hyenas moved faster (longer step lengths) when closer to the fence. Lions showed no variation in path straightness with respect to the fence but hyenas had straighter movement paths when closer to the fence. When moving within the vicinity of the fence, lions had a 9% likelihood of crossing the fence, while hyenas were much more likely (18%) to do so, which suggests that the Etosha fence was more permeable for hyenas than for lions. Fence crossings predominantly occurred at night for both species. Lions were more likely to cross the fence during the cold dry season (May–August) whereas hyenas crossed more often during the warm wet season (December–April). Female and adult lions had a higher probability of crossing the fence than males and subadults. These findings could be useful in the development of appropriate strategies to mitigate human–carnivore conflicts and thus promote the conservation of carnivores in multiple-use areas along the periphery of the park. Further research is needed to better elucidate the impact of perimeter fences on a broader suite of carnivores and other large mammals, including species-specific variation in the response to fences, as well as the economic benefits of fencing provided to local communities.

59 BASIC BIOLOGICAL SCIENCES↗

First-principles comparative study of Cr migration in O3 and O/P hybrid-phased NaCrO 2

In layered Na transition-metal (TM) oxides, TM migration usually occurs at highly charged states and severely deteriorates the capacity and reversibility. Meanwhile, the formation of hybrid phases with the intergrowth of octahedral (O-type) and prismatic (P-type) Na layers also takes place at highly charged states. These hybrid phases are often more stable than simple O3 or P3 stackings. However, little is known about the mechanism and impact of TM migration in these hybrid phases. In this work, a comparative first-principles study is performed to understand the connections between structural changes and Cr migration in layered O3 and hybrid-phased NaCrO 2 . After Cr migration, the hybrid-phased NaCrO 2 suffers from greater layer shrinkage than the O3 phase. Three factors affect the Cr migration energy E mig : the Na concentration, local 3D configurations, and 2D in-plane geometries. Low Na concentrations and certain local 3D configurations facilitate the Cr migration. The Cr migration barriers in both O3 and hybrid-phased NaCrO 2 are positively correlated with the Cr E mig . The Cr migration in 17 doped O3 and hybrid-phased NaCrO 2 is surveyed. In these doped NaCrO 2 , a more uniform distribution of the Cr–O bond lengths usually suggests suppressed Cr migration. Optimal dopants for suppressing Cr migration are identified by considering E mig for both Cr and the dopant. Finally, our comparative study on Cr migration in O3 and hybrid-phased NaCrO 2 reveals the significant role of hybrid-phased structures in the development of layered cathode materials.

36 MATERIALS SCIENCE↗

Evaluating an Earth system model from a water manager perspective

Abstract. The large spatial scale of global Earth system models (ESMs) is often cited as an obstacle to using the output by water resource managers in localized decisions. Recent advances in computing have improved the fidelity of hydrological responses in ESMs through increased connectivity between model components. However, the models are seldom evaluated for their ability to reproduce metrics that are important for and resonate with practitioners or that allow practitioners to situate higher-resolution model outputs within a cascade of uncertainty stemming from different models and scenarios. We draw on the combined experience of the author team and water manager workshop participants to identify salient water management metrics and evaluate whether they are credibly reproduced over the conterminous USA by the Community Earth System Model v2 (CESM2) Large Ensemble. We find that, while the exact values may not match the observations, aspects such as interannual variability can be reproduced by CESM2 for the mean wet day precipitation and length of dry spells. CESM2 also captures the proportion of total annual precipitation that derives from the heaviest rain days in watersheds that are not snow-dominated. Aggregating the 7 d mean daily runoff to two-digit Hydrological Unit Code (HUC2) watersheds also shows that rain-dominated regions capture the timing and interannual variability of annual maximum and minimum flows. We conclude that there is potential for far greater use of large-ensemble ESMs, such as CESM2, in long-range water management decisions to supplement high-resolution regional projections.

Tye, Mari R. (ORCID:0000000324911020)↗

Layer-by-layer growth of graphene oxide multilayers using robust interlayer linking chemistry. 1. Zr-bisphosphates

Graphene oxide (GO) holds substantial interest because of its utility in applications ranging from chemical sensing to catalysis and energy storage. We report on the layer-by-layer formation of graphene oxide phosphate (P-GO) multilayers on silica and oxidized silicon supports, with the layers connected via Zr-bisphosphate (ZP) linking chemistry. Layers were grown either directly by ZP linkages between P-GO sheets, or with the use of 1,4-phenylene bisphosphate (BP1) or [1,1′-biphenyl]4,4′-diyl bisphosphate (BP2). The layers have been characterized using optical null ellipsometry, X-ray photoelectron spectroscopy (XPS), UV-visible absorption spectroscopy and scanning electron microscopy (SEM). There is a linear dependence of assembly thickness on the number of layers deposited and the integrated area of the optical absorbance bands also increases linearly with number of layers deposited. XPS data provides information on interlayer linking stoichiometry. SEM images provide insight into the morphology of the adlayers, suggesting the structure and length of the interlayer linking moieties used in the multilayer assembly play a significant role in the organization of the resulting system.

Mukherjee, Neelanjana [Michigan State University, ↗

Universal scaling of electron transmission for nearly ballistic and quantum dragon nanodevices

Here, we predict two different universal scaling regimes for the quantum transmission of metallic nanodevices following the addition of a small amount of uncorrelated disorder. A nanodevice is connected to two thin semi-infinite uniform leads, and the Non-Equilibrium Green’s Function (NEGF) methodology yields the electron transmission $\mathscr{T}(E)$ as a function of the injected electron energy $E$. Ballistic nanodevices have no disorder and have $\mathscr{T}(E)=1$for all $E$ that allow electron propagation in the leads. Quantum dragon nanodevices can have extremely strong properly correlated disorder, and still have $\mathscr{T}(E)=1$for all $E$. Additional uncorrelated site disorder leads to Fano resonances in $\mathscr{T}(E)$. Averaging over the uncorrelated disorder we predict using perturbation theory two universal scaling regimes for $\mathscr{T}_{ave}(E)$. The functional form of both universal scaling regimes depend on the device length and width, energy, and variance of the uncorrelated disorder. The second scaling regime, valid for small but somewhat larger uncorrelated disorder than the first scaling regime, also has the form dependent on the density of states of the system. These two scaling regimes are demonstrated to be valid via large scale computer calculations.

Ballistic transport↗

Tradeoffs and Synergies in Tropical Forest Root Traits and Dynamics for Nutrient and Water Acquisition: Field and Modeling Advances

Vegetation processes are fundamentally limited by nutrient and water availability, the uptake of which is mediated by plant roots in terrestrial ecosystems. While tropical forests play a central role in global water, carbon, and nutrient cycling, we know very little about tradeoffs and synergies in root traits that respond to resource scarcity. Tropical trees face a unique set of resource limitations, with rock-derived nutrients and moisture seasonality governing many ecosystem functions, and nutrient versus water availability often separated spatially and temporally. Root traits that characterize biomass, depth distributions, production and phenology, morphology, physiology, chemistry, and symbiotic relationships can be predictive of plants’ capacities to access and acquire nutrients and water, with links to aboveground processes like transpiration, wood productivity, and leaf phenology. In this review, we identify an emerging trend in the literature that tropical fine root biomass and production in surface soils are greatest in infertile or sufficiently moist soils. We also identify interesting paradoxes in tropical forest root responses to changing resources that merit further exploration. For example, specific root length, which typically increases under resource scarcity to expand the volume of soil explored, instead can increase with greater base cation availability, both across natural tropical forest gradients and in fertilization experiments. Also, nutrient additions, rather than reducing mycorrhizal colonization of fine roots as might be expected, increased colonization rates under scenarios of water scarcity in some forests. Efforts to include fine root traits and functions in vegetation models have grown more sophisticated over time, yet there is a disconnect between the emphasis in models characterizing nutrient and water uptake rates and carbon costs versus the emphasis in field experiments on measuring root biomass, production, and morphology in response to changes in resource availability. Closer integration of field and modeling efforts could connect mechanistic investigation of fine-root dynamics to ecosystem-scale understanding of nutrient and water cycling, allowing us to better predict tropical forest-climate feedbacks.

54 ENVIRONMENTAL SCIENCES↗

Chemical timescale effects on detonation convergence

Numerical simulations of detonation-containing flows have emerged as crucial tools for designing next- generation power and propulsion devices. As these tools mature, it is important for the combustion community to properly understand and isolate grid resolution effects when simulating detonations. To this end, the objective of this work is to provide a comprehensive analysis of the numerical convergence of unsteady detonation simulations, with focus on isolating the impacts of chemical timescale modifications on convergence characteristics in the context of operator splitting. With the aid of an AMReX-based adaptive mesh refinement flow solver-which enables resolutions up to ($\mathcal{O}$ (1000) cells-per-induction length-the convergence analysis is conducted using two kinetics configurations: (1) the simplified three-step Arrhenius-based model mechanism of Short and Quirk (1997), where chemical timescales in the detonation are modified by adjusting activation energies in the initiation and branching reactions, and (2) a detailed hydrogen- air mechanism, where the chemical timescales are adjusted by varying the ambient pressure. The convergence of unsteady self-sustained detonations in one-dimensional channels is then analyzed with reference to steady-state theoretical baseline solutions using these mechanisms. The goal of the analysis is to provide a detailed comparison of the effects of grid resolution on both macroscopic (peak pressures and wave speeds) and microscopic (wave structure) quantities of interest, drawing connections between the deviations from steady-state baselines and minimum chemical timescales. In particular, chemical timescale reductions were found to have minimal impact on the convergence of macroscopic properties. Furthermore, analyses of microscopic convergence trends, particularly in the reaction front location, revealed a key insight: maintaining the induction time while eliminating prohibitive chemical timescales through mechanism simplifications and combustion modeling can significantly enhance detonation convergence properties. Ultimately, this work uncovers resolution-dependent unsteady detonation convergence regimes and highlights the important role played by not only the chemical timescales, but also the ratio between the chemical timescale and induction time on the numerical convergence of the detonation wave structure.

Adaptive mesh refinement↗

Bridging Atomic Solvation Environment with Electrochemical Properties for the Bis(trifluoromethylsulfonyl)imide-Based Divalent Cation Electrolytes for the Next-Generation Energy Storage Systems

A deep molecular-level understanding of the multivalent electrolyte and its correlation with the electrochemical properties is crucial for designing optimized electrolytes for next-generation rechargeable batteries. Comprehensive knowledge of the atomic level of the solvation structure and its connection with electrochemical stability and ion transport properties is especially critical. However, the interaction of these three components coupled with clear atomistic insights is lacking in the literature. Here, our current contribution evaluates representative electrolytes with the bis(trifluoromethanesulfonyl)imide (TFSI) anions for multivalent cations of Mg, Ca, and Zn, at different ionic conditions with and without a cosolvated environment in ether-based solvent. Two critical problems are investigated: first, resolving the solvation structures in the electrolyte solutions as a function of concentrations through pair distribution function analysis and the corresponding electrochemical transport properties; second, unmasking the quantitative correlation of the atomistic environment with both electrochemical kinetics and cation dependence. We discovered that the magnesium- and calcium-based electrolytes display versatile coordination lengths but poor average anodic stability due to ion pairing with TFSI - . On the contrary, the zinc-based electrolytes show the shortest solvent coordination lengths, shielding the Zn cation from rigid solvent interactions and resulting in the highest anodic stabilities. Calcium-based electrolytes exhibit the longest and most concentration-independent coordination lengths. This work provides valuable insights into the molecular structural and electrochemical features of diverse multivalent electrolyte systems with cations in various solvation environments, emphasizing the importance of the solvation structure and construction in designing high-performance electrolytes.

cation coordination↗

Designing Peptide Fossils That Model the Evolution of the Bacterial Ferredoxin Fold

Electron transfer coupled to redox chemistry is at the heart of metabolism. The proteins responsible for moving electrons (protein electron carriers) must have emerged at the origin of life. The small iron–sulfur-binding bacterial ferredoxins were likely among these first proteins. Embedded within the ferredoxin sequence and structure is a symmetry that points to an ancient gene duplication event. Little is understood about the nature of ferredoxins prior to this duplication event or what environmental factors may have driven the selection for more complex forms. The deep-time molecular history of ferredoxins goes back billions of years and cannot be reconstructed by phylogenetic analyses based on amino acid sequences. Here, we use structure-guided protein design to model a fossil half-ferredoxin stage in the evolution of this fold, the semidoxins, and their symmetric full-length counterparts, the symdoxins. Semidoxin designs homodimerize, exhibiting structural, thermodynamic, and electrochemical behaviors in most cases identical to cognate symdoxins. However, the semi- and symdoxin fossil stages behave differently when incorporated into an in vivo electron transfer complementation assay. Both can support bacterial growth dependent on protein expression. Growth rates of bacteria expressing the semidoxins are much more sensitive to oxygen than those of bacteria expressing symdoxins. Motivated by the in vivo functionality of designed semidoxins, we identified putative naturally occurring semidoxins in extant anaerobic microorganisms. This is consistent with the observed in vivo oxygen sensitivity of the semidoxin designs. One natural semidoxin is shown to be folded and redox active. However, it exists as a mixture of monomers and dimers, suggesting a potential connection between semidoxins and even simpler single iron–sulfur cluster-binding peptides.

59 BASIC BIOLOGICAL SCIENCES↗

Understanding interfacial segregation in polymer blend films with random and mixed side chain bottlebrush copolymer additives

Bottlebrush polymers are complex macromolecules with tunable physical properties dependent on the chemistry and architecture of both the side chains and the backbone. Prior work has demonstrated that bottlebrush polymer additives can be used to control the interfacial properties of blends with linear polymers but has not specifically addressed the effects of bottlebrush side chain microstructures. Here, using a combination of experiments and self-consistent field theory (SCFT) simulations, we investigated the effects of side chain microstructures by comparing the segregation of bottlebrush additives having random copolymer side chains with bottlebrush additives having a mixture of two different homopolymer side chain chemistries. Specifically, we synthesized bottlebrush polymers with either poly(styrene-ran-methyl methacrylate) side chains or with a mixture of polystyrene (PS) and poly(methyl methacrylate) (PMMA) side chains. Furthermore, the bottlebrush additives were matched in terms of PS and PMMA compositions, and they were blended with linear PS or PMMA chains that ranged in length from shorter to longer than the bottlebrush side chains. Experiments revealed similar behaviors of the two types of bottlebrushes, with a slight preference for mixed side-chain bottlebrushes at the film surface. SCFT simulations were qualitatively consistent with experimental observations, predicting only slight differences in the segregation of bottlebrush additives driven by side chain microstructures. Specifically, these slight differences were driven by the chemistries of the bottlebrush polymer joints and side chain end-groups, which were entropically repelled and attracted to interfaces, respectively. Using SCFT, we also demonstrated that the interfacial behaviors were dominated by entropic effects with high molecular weight linear polymers, leading to enrichment of bottlebrush near interfaces. Surprisingly, the SCFT simulations showed that the chemistry of the joints connecting the bottlebrush backbones and side chains played a more significant role compared with the side chain end groups in affecting differences in surface excess of bottlebrushes with random and mixed side chains. This work provides new insights into the effects of side chain microstructure on segregation of bottlebrush polymer additives.

36 MATERIALS SCIENCE↗

RivGraph: Automatic extraction and analysis of river and delta channel network topology

River networks sustain life and landscapes by carrying and distributing water, sediment, and nutrients throughout ecosystems and communities. At the largest scale, river networks drain continents through tree-like tributary networks. At typically smaller scales, river deltas and braided rivers form loopy, complex distributary river networks via avulsions and bifurcations.In order to model flows through these networks or analyze network structure, the topology, or connectivity, of the network must be resolved. Additionally, morphologic properties of each river channel as well as the direction of flow through the channel inform how fluxes travel through the network’s channels. Riv Graphis a Python package that automates the extraction and characterization of river channel networks from a user-provided binary image, or mask, of a channel network (Fig. 1). Masks may be derived from (typically remotely-sensed) imagery, simulations, or even hand-drawn. RivGraph will create explicit representations of the channel network by resolving river centerlines as links, and junctions as nodes. Flow directions are solved for each link of the network without using auxiliary data, e.g., a digital elevation model (DEM). Morphologic properties are computed as well, including link lengths, widths, sinuosities, branching angles,and braiding indices. If provided,RivGraph will preserve georeferencing information of the mask and will export results as ESRI shapefiles, GeoJSONs, and GeoTIFFs for easy import into GIS software.RivGraph can also return extracted networks as networkx objects for convenient interfacing with the full-featured networkx package (Hagberg et al., 2008). Finally, RivGraph offers a suite of topologic metrics that were specifically designed for river channel network analysis (Tejedor et al., 2015b).

54 ENVIRONMENTAL SCIENCES↗