Search NASASearch

SEARCH · Search NASA

Results for “12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE WASTES FROM NUCLEAR FACILITIES”

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.

3,403 records · Page 72

Second LDEF Post-Retrieval Symposium Abstracts

These abstracts from the symposium represent the data analysis of the 57 experiments flown on the LDEF. The experiments include materials, coatings, thermal systems, power and propulsion, science, (cosmic ray, interstellar gas, heavy ions, micrometeoroids, etc.), electronics, optics, and life science.

Arlene S. Levine

Ejecta Generation and Redistribution on 433 Eros: Modeling Ejecta Launch Conditions

The NEAR-Shoemaker mission to asteroid 433 Eros presents an unprecedented opportunity to gain fundamental new knowledge about the processes governing regolith formation and redistribution on small bodies. NEAR-Shoemaker’s high-resolution imaging of the surface of Eros makes the asteroid a valuable and heretofore unparalleled laboratory for the detailed study of impact ejecta reaccretion and regolith redistribution on low-gravity (of order 10 -3 g) objects. Regolith is produced on asteroids by impact cratering, and the existence of regolith on the smallest solar system bodies supports the view that some of the ejecta from impact events on such objects may be retained. Impact craters and retained ejecta on low-gravity objects like Eros represent valuable natural laboratories for evaluating various models of impact cratering processes, since they may present crater structures or ejecta features that either do not form or are hidden on higher-gravity bodies like the Moon. Further, quantifying the extent to which impact processes generate and redistribute regoliths on small body surfaces (excavation depths, retained fraction, turnover timescales, etc.) is pivotal to the issue of how to relate meteoritical samples to their asteroidal parent bodies when surficial processes ( i.e., “space weathering”) may disguise or cover up underlying material and confound the ability of remote sensing techniques to provide reliable mineralogical assays of the parent objects. The rich variety of data on Eros’ regolith properties and distribution returned by NEAR-Shoemaker now require detailed analysis in order to take full advantage of the clues these observations offer for elucidating details of the impact cratering process on small bodies. Complicating simple interpretations of crater and ejecta morphology are dynamical effects on ejecta emplacement resulting from Eros’ irregular shape, rapid (5.27 hr) rotation, and low gravity. Figure 1 shows the very different ejecta deposit morphology that can result if the effects of rotation alone are neglected. Considering the additional complicating factors of Eros’ irregular shape and complex gravitational field, simple calculations of the extent and thickness of ejecta blankets and the spatial distribution of ejecta blocks from basic crater scaling laws or numerical hydrocodes alone do not suffice. In order to fully interpret the suite of NEAR-Shoemaker observations of regolith features across the surface of Eros and to evaluate various impact models for specific craters on the asteroid, detailed dynamical modeling of the deposition of crater ejecta from those craters is required . Here, I describe some modifications and improvements to the dynamical model being used for these studies.

D D Durda

A Review of Pump Package Assembly Anomalies and Demonstration to Gain Motor Controller Replacement Capability

The International Space Station (ISS) generates heat within the US modules from systems used to operate the station, maintain crew health, and conduct science experiments in the orbiting laboratory. The Internal Thermal Control System (ITCS) collects this heat in water loops and transfers it to the external ammonia system for rejection to space. The ITCS uses six Pump Package Assemblies (PPA) as the motive force to circulate the water within the Lab, Node 2, and Node 3 (Node 1 is parasitic to the Lab and Node 3 ITCS loops). In the life of ISS, 3 PPAs have had to be replaced and 2 PPAs are currently showing signs of degradation. The first PPA to fail in 2003 was due to early coolant chemistry issues; the remaining pump failures and current pump degradations are determined to be related to the Pump/Fan Motor Controller (PFMC), with their individual signatures all being unique to this point. The PFMC of one of the failed PPAs was returned and testing confirmed its sub-components to be the root cause of the PPA failure. Due to the relative size and weight of the standalone PFMC in comparison to the integrated assembly, it is desired, from a launch and on-orbit stowage perspective, to replace the PFMC component of the degraded PPA instead of the entire PPA. The PPA was not designed for this capability, so on the ground demonstrations are planned using crew tools to show that the PFMC can be replaced on-orbit while maintaining full functionality. This paper will focus on the signatures of the two most recent PPA failures and the two PPAs showing signs of degradation, in addition to presenting the findings from the PFMC demonstration effort.

Aaron Rodriguez

A Review of Pump Package Assembly Anomalies and Demonstration to Gain Motor Controller Replacement Capability

The International Space Station (ISS) generates heat within the US modules from systems used to operate the station, maintain crew health, and conduct science experiments in the orbiting laboratory. The Internal Thermal Control System (ITCS) collects this heat in water loops and transfers it to the external ammonia system for rejection to space. The ITCS uses six Pump Package Assemblies (PPA) as the motive force to circulate the water within the Lab, Node 2, and Node 3 (Node 1 is parasitic to the Lab and Node 3 ITCS loops). In the life of ISS, 3 PPAs have had to be replaced and 2 PPAs are currently showing signs of degradation. The first PPA to fail in 2003 was due to early coolant chemistry issues; the remaining pump failures and current pump degradations are determined to be related to the Pump/Fan Motor Controller (PFMC), with their individual signatures all being unique to this point. The PFMC of one of the failed PPAs was returned and testing confirmed its sub-components to be the root cause of the PPA failure. Due to the relative size and weight of the standalone PFMC in comparison to the integrated assembly, it is desired, from a launch and on-orbit stowage perspective, to replace the PFMC component of the degraded PPA instead of the entire PPA. The PPA was not designed for this capability, so on the ground demonstrations are planned using crew tools to show that the PFMC can be replaced on-orbit while maintaining full functionality. This paper will focus on the signatures of the two most recent PPA failures and the two PPAs showing signs of degradation, in addition to presenting the findings from the PFMC demonstration effort.

Aaron Rodriguez

Design Trade-Offs in Composite Fuel Cell Membranes: Effects of Reinforcement and Chemical Additives

Perfluorosulfonic acid (PFSA) membranes are critical components in proton exchange membrane fuel cells, where performance depends on balancing ionic conductivity, mechanical durability, and chemical stability. This study characterizes a composite membrane (NC700) featuring PFSA-impregnated expanded polytetrafluoroethylene (ePTFE) reinforcement and cerium-based radical scavengers, benchmarked against unreinforced NR211. Complementary techniques, including electron microscopy, X-ray scattering, infrared spectroscopy, thermogravimetric analysis, and dynamic mechanical analysis, identify the structural and compositional strategies employed in NC700. Water sorption isotherms reveal lower water uptake for NC700 across all conditions, attributed to reinforcement and cerium incorporation. Reinforcement reduces in-plane swelling from 11% to 2.1% at 90% RH, confirming strong swelling anisotropy, while maintaining mechanical properties at elevated temperatures. While the ionic conductivity of NC700 is approximately 10% lower than that of NR211, the reduced thickness yields a 40% decrease in calculated area-specific resistance, suggesting the composite architecture can favorably shift the conductivity-stability trade-off. The composite structure also reduces gas permeability, indicating potential for improved separator function alongside favorable transport properties. Systematic deconvolution of reinforcement and additive contributions shows that conductivity losses from cerium incorporation are largely offset by gains from the lower equivalent-weight polymer, providing quantitative relationships that may guide composite membrane design for fuel cells and other electrochemical applications.

25 ENERGY STORAGE

Improving the Performance of NEML2 with Modern Graph Compilation Backends

NEML2 vectorizes constitutive-model evaluation for large-scale multiphysics simulation, using PyTorch as its tensor backend so that a batch of material-point updates runs on CPU or GPU through a single implementation. In the two prior reports in this series it was a C++-native library, deployed through TorchScript tracing and just-in-time (JIT) compilation; it has since been rewritten from the ground up into a Python-native library deployed through Ahead-of-Time Inductor (AOTInductor), a modern PyTorch graph-compilation backend. The rewrite is driven by a persistent tension, not a language preference: NEML2 composes constitutive models at runtime from a registry of small, independently-authored pieces, and that flexibility is difficult to reconcile with the compile-time knowledge an efficient GPU kernel needs. This report documents the rewrite and the investment that accompanied it: the AOTInductor export pipeline that turns a Python-authored model into a portable, Python-free compiled artifact loadable from pure C++; the eager and compiled runtimes and the new implicit solver layer built on them; a head-to-head benchmark of legacy JIT against AOTInductor; the physics-model catalog and its worked examples; the developer tooling; and the corresponding overhaul of MOOSE’s NEML2 integration that lets MOOSE consume it. A central objective is to examine whether modern PyTorch graph-compilation backends are effective for MOOSE GPU integration. The benchmark answers directly: AOTInductor outperforms legacy JIT on every GPU scenario measured, by 1.0–4.5×. Modern graph-compilation backends are effective for MOOSE GPU integration, and AOTInductor specifically – not compilation in the abstract – is why.

Hu, Gary (Tianchen) [Argonne National Laboratory (

Predicting Operational Performance of xEMU Boot at Lunar South Pole Temperatures using Thermal Desktop ®

The spacesuit boots that will be used on Artemis lunar south pole surface missions will be exposed to extremely cold temperatures (down to ~50 K). To assess the performance of the government’s Exploration Extravehicular Mobility Unit (xEMU) lunar boot in these permanently shadowed regions, testing was performed at the Jet Propulsion Lab (JPL) in the Cryogenic Ice Transfer, Acquisition Development, and Excavation Laboratory (CITADEL) thermal vacuum (TVAC) chamber. This paper documents the data analysis, thermal boot model correlation, and operational predictions conducted using data from the xEMU CITADEL TVAC test. Expected thermal conductances within the boot and between the boot and environment were calculated from test data, which was then used as an initial guess for conductances within a Thermal Desktop (TD) model. Correlation of the TD model using the internal SOLVER feature was performed across 10 different test points which varied external temperature, internal boot ventilation flowrate, and contact pressure. Operational performance at the lunar south pole was then predicted using results from the correlated model. While the predictions provide evidence for acceptable performance of the boots at the 100K environment test point, there is still substantial uncertainty in performance, especially at the 48K test point. This uncertainty is due in part to testing limitations such as contacting the foot to a hard metal plate rather than granular regolith, and model limitations such as the lack of a realistic foot model. These limitations and their impacts are addressed in detail in this paper. The results of this test series and model correlation underscore the importance of additional improved testing and modeling for characterizing the expected thermal resistance between the outside of the boot and the lunar surface.

Spacesuit

Thermal Design

An essential phase of the Explorer XIII development was a thermal study and design to provide an acceptable temperature environment for the electronics components and external surfaces of the satellite during ascent and in orbit. It was necessary to establish by preflight analysis and tests that tolerable temperatures could be maintained in three regimes of flight. In the first regime (from launch until release of the heat shield) the satellite was subjected to radiative and conductive heat from the shield. In the second regime, after release of the shield at 350,000 feet, the satellite was heated by free molecular flow. In the third regime - the orbit phase - it was necessary that temperatures within limits for a 1-year lifetime be established. In addition to these three regimes, a study was made to investigate the effect of elevated rocket-motor temperatures during launch and after burnout. This chapter will deal with some pertinent preflight estimates and correlation of these estimates with flight data.

Explorer XIII Satellite

Deep Learning and Photogrammetric Reconstruction for Automated Crack Detection and Dimensional Measurement in Mining Operations

Surface crack detection and dimensional measurement at active mining sites present significant safety and operational challenges. Manual inspection methods are labor-intensive, spatially incomplete, and expose personnel to hazardous environments, while existing automated approaches have been developed primarily for concrete civil infrastructure and have not been validated on the complex, variable surfaces characteristic of mining environments. This dissertation presents an automated pipeline that integrates deep learning semantic segmentation with Structure-from-Motion photogrammetry to detect surface cracks and measure their aperture, length, and vertical displacement from standard RGB imagery acquired during routine Uncrewed Aerial Vehicle (UAV) survey operations, without requiring additional sensor hardware or manual measurement. The pipeline combines a U-Net architecture with an EfficientNet-B0 encoder, pretrained on the SDNET2018 concrete crack dataset and fine-tuned on a mining-specific dataset spanning laboratory concrete specimens, coal refuse impoundment embankments, and post-blast limestone quarry benches. Photogrammetric reconstruction is performed using COLMAP Structure-from-Motion and Multi-View Stereo, with crack segmentation masks projected into the reconstructed point cloud to enable three-dimensional vertical displacement measurement through local plane fitting and bimodal surface detection. The pipeline was validated across 36 controlled laboratory specimens at three imaging distances and four vertical displacement levels, achieving aperture measurement RMSE of 0.047 cm and R² of 0.954, and vertical displacement RMSE of 0.140 cm and R² of 0.966, against independent caliper measurements. Field application at a coal refuse impoundment in southwestern Pennsylvania detected 71 crack components across the embankment crest, with a dominant longitudinal crack exhibiting aperture values reaching 28 cm and a 95th percentile vertical displacement of 35.53 cm, consistent in magnitude and spatial distribution with simultaneously acquired LiDAR-derived estimates. Application across four post-blast limestone quarry bench datasets in California successfully characterized blast-induced fracture networks at ground sampling distances ranging from 0.59 to 1.23 cm/pixel, with detected crack geometries physically consistent with observable surface conditions at each site. The results demonstrate that deep learning-based crack detection and photogrammetric measurement can be integrated into routine UAV inspection workflows at mining sites, providing repeatable, scalable, and quantitative crack characterization across surface types, crack scales, and displacement magnitudes not previously addressed in the literature. The pipeline requires no dedicated surveying equipment beyond the UAV platforms already deployed at mine sites for survey and monitoring purposes, supporting practical adoption within existing operational workflows.

Crack detection, Dimensional Measurement

Beyond Melting: Amorphous Bonding for Joining and Consolidation

Crystallization may be the hidden constraint in thermoplastic composite manufacturing. It requires tightly controlled cooling, induces residual stresses through shrinkage, and introduces path-dependent behavior that complicates predictive modeling yet remains essential for structural performance. This work asks: can bonding be achieved without relying on melt-driven crystallization? To address this, thin (5–20 μm) polyetherimide (PEI) interlayers are pre-healed to slow-cooled polyaryletherketone (PAEK) in two contexts. The first, Thermabond®, is sub-melt joining of low melt-PAEK laminates. Results show that bond quality is governed primarily by processing (i.e., adequate healing and film handling) rather than modest changes in interlayer thickness. This concept is then extended to laminate-scale manufacturing through an architecture known as OATMEAL (Out-of-autoclave Amorphous/semicrystalline Thermoplastic Material for Energy-efficient Aerospace-grade Laminates). PEI is healed to carbon fiber reinforced polyetheretherketone (PEEK) at the prepreg and excess PEI is then ablated from the surface. Crystallinity is developed off-line during prepreg fabrication, while subsequent consolidation occurs below the melt temperature to preserve it. Cross-ply warpage experiments show that, contrary to intuition, repeated amorphous interfaces reduce global curvature by lowering the effective stress lock-in temperature and eliminating crystallization shrinkage from the lamina response. Correspondingly, laminate behavior is accurately predicted using classical laminate theory (CLT) with a single effective stress-free temperature, whereas conventional CF/PEEK requires accounting for crystallization-driven effects. By decoupling interfacial healing from crystallization, OATMEAL enables sub-melt consolidation, reduces energy consumption by up to 75%, and increases manufacturing throughput by fivefold. These results demonstrate that amorphous bonding is not only a joining strategy, but a pathway to more predictable and scalable thermoplastic composite manufacturing.

solidification

Weak localization and universal conductance fluctuations in large-area twisted bilayer graphene

We study diffusive magnetotransport in highly 𝑝-doped large-area twisted bilayer graphene in 1∘, 7∘, 9∘, and 20∘ samples. All samples exhibit weak localization, from which we extract the phase coherence length and intervalley scattering lengths, and from that determine that dephasing is caused by electron-electron scattering and intervalley scattering is caused by point defects. We observe signatures of universal conductance fluctuations in the 9∘ sample, which has high mobility and is near the van Hove singularity. Further improvements in sample quality and applications to large-area moiré materials will open new avenues to observe quantum interference effects.

Talkington, Spenser [University of Pennsylvania]

A Predictive Bubble Point Pressure Model for Porous Liquid Acquisition Device Screens

This article presents a simplified model for porous screen channel liquid acquisition devices based on a maximum bubble point pressure method from Adamson and Gast (1997). To validate the model, three 304 stainless steel (325 × 2300, 450 × 2750, and 510 × 3600) mesh samples were tested in methanol, acetone, isopropyl alcohol, and water. Screen pores are estimated based on analysis from scanning electron microscopy, historical data, and current test data. Results show that the bubble point pressure is proportional to the surface tension of the fluid only when accounting for nonzero contact angles. The previous assumption that bubble point pressure scales inversely with effective pore diameter is shown to be invalid, as the second finest 450 × 2750 produced the highest bubble point of the three screens. The simplified bubble point model can be used to make predictions for any pure fluid when pore diameters are based on bubble point tests and not SEM analysis.

Liquid Acquisition Device

Aragonite saturation horizon variability along North Pacific seamounts and implications for deep-sea coral reefs

The 2014 discovery of living deep-sea coral reefs along the Northwest Hawaiian Islands (NWHI) and lower Emperor Seamount Chain (ESC), despite the North Pacific’s shallow aragonite saturation horizon (ASH) and high CaCO3 dissolution rates, underscores the need to understand the local seawater chemistry that allows these reefs to persist. We investigated seawater carbonate chemistry using discrete samples along the NWHI and ESC from two cruises ~1 year apart (08/26/21 – 09/26/21, 09/09/22 – 10/24/22). Across the two cruises, ASH depth difference ranged from 15 to 77 m. Since the Pacific ASH shoals by 1–2 m yr?¹, this long-term trend cannot explain the magnitude of ASH change observed between cruises. Temperature-salinity plots similarly indicate no major shifts in intermediate water masses that could account for these changes. Instead, ASH depth variability was primarily governed by localized biogeochemical processes, namely changes in intermediate water respiration and CaCO3 dissolution. Indicators for dissolution (TA*) and respiration (AOU) suggest changes in ASH depth were driven by changes in dissolution at the northern- and southern-most sites, whereas respiration exerted stronger control at central sites. Combining 2021 and 2022 data with data from 2014 – 2019 (excluding 2018) revealed high interannual ASH variability, by as much as >200 m at one site. Deep-sea coral reefs across the NWHI and ESC currently reside close to the ASH depth and likely experience seasonal and interannual shifts between under- and supersaturation. As ocean acidification progresses however, persistent exposure to undersaturated seawater could further threaten these critical deep-sea ecosystems.

coral reefs

Host galaxies of ultra-strong Mg ii absorbers at z ∼ 0.7

ABSTRACT We report spectroscopic identification of the host galaxies of 18 ultra-strong Mg ii systems (USMg ii) at 0.6 ≤ z ≤ 0.8. We created the largest sample by merging these with 20 host galaxies from our previous survey within 0.4 ≤ z ≤ 0.6. Using this sample, we confirm that the measured impact parameters ($\rm 6.3\leqslant D[kpc] \leqslant 120$ with a median of 19 kpc) are much larger than expected, and the USMg ii host galaxies do not follow the canonical $\rm {\it W}_{2796}-{\it D}$ anticorrelation. We show that the presence and significance of this anticorrelation may depend on the sample selection. The $\rm {\it W}_{2796}-{\it D}$ anticorrelation seen for the general Mg ii absorbers show a mild evolution at low $\rm W_{2796}$ end over the redshift range 0.4 ≤ z ≤ 1.5 with an increase of the impact parameters. Compared to the host galaxies of normal Mg ii absorbers, USMg ii host galaxies are brighter and more massive for a given impact parameter. While the USMg ii systems preferentially pick star-forming galaxies, they exhibit slightly lower ongoing star-forming rates compared to main sequence galaxies with the same stellar mass, suggesting a transition from star-forming to quiescent states. For a limiting magnitude of mr < 23.6, at least 29 per cent of the USMg ii host galaxies are isolated, and the width of the Mg ii absorption in these cases may originate from gas flows (infall/outflow) in isolated haloes of massive star forming but not starbursting galaxies. We associate more than one galaxy with the absorber in $\ge 21~{{\ \rm per\ cent}}$ cases, where interactions may cause wide velocity spread.

Astronomy & Astrophysics

Comparison of Cryogenic Flow Boiling in Liquid Nitrogen and Liquid Hydrogen Chilldown Experiments

This paper presents a comparison between experimental results from recent liquid hydrogen (LH 2 ) transfer line chilldown experiments at high Reynolds (Re) numbers versus liquid nitrogen (LN 2 ) experiments conducted at low Re numbers. Parasitic heat leak, inner wall temperatures, inner wall heat fluxes, and heat transfer coefficients are computed to compare between the two systems. Analysis of temperature traces and flow visualization indicates that the chilldown process evolves much more rapidly at higher Re numbers due to a quick transition from vapor flow to annular liquid flow and near immediate liquid contact along the pipe walls. The lower kinematic viscosity and surface tension of LH 2 , along with reduced parasitic heat leak and higher Re numbers relative to the LN 2 experiments, causes chilldown to proceed almost immediately into the nucleate boiling regime, in comparison to low Re flows where >75% of the chilldown is spent in vapor film boiling.

Cryogenic Heat Transfer Coefficient

The Flight Performance of the Galileo Orbiter USO

Results are presented from an analysis of radio metric data received by the DSN stations from the Galileo spacecraft using an Ultrastable Oscillator (USO) as a signal source. These results allow the health and performance of the Galileo USO to be evaluated, and are used to calibrate this Radio Science instrument and the data acquired for Radio Science experiments such as the Red-shift Observation, Solar Conjunction, and Jovian occultations. Estimates for the USO-referenced spacecraft-transmitted frequency and frequency stability were made for 82 data acquisition passes conducted between launch (October 1989) and November 1991. Analyses of the spacecraft-transmitted frequencies show that the USO is behaving as expected. The USO was powered off and then back on in August 1991 with no adverse effect on its performance. The frequency stabilities measured by Allan deviation are consistent with expected values due to thermal wideband noise and the USO itself at the appropriate time intervals. The Galileo USO appears to be healthy and functioning normally in a reasonable manner.

D D Morabito

Catalytic disproportionation on carbon superstructures enables long-life, high-loading Li–S batteries

Electrocatalysis has been widely explored as an effective strategy to accelerate polysulfide (PS) conversion and suppress the shuttle effect in lithium–sulfur (Li–S) batteries. However, the underlying mechanisms remain elusive, and electrocatalytic reactions are inactive during cell resting. In this work, we reveal and quantitatively analyze a previously unrecognized sulfur reduction route (SRR) driven by catalytic disproportionation at the carbon cathode surface—fundamentally distinct from conventional electrocatalysis. Unlike conventional stepwise pathways, this SRR enables high-order polysulfides (Sₓ²⁻, x = 5–8) to directly convert into S₈ and Li₂S₂, bypassing low-order intermediates. This sulfur-reduction shortcut is systematically elucidated through high-performance liquid chromatography, revealing the intrinsic catalytic contribution of carbon frameworks and the dynamic evolution of PS species. We demonstrate that carbon superstructures (CSS-0.5), assembled from nanosheet subunits with abundant N/O functionalities and interconnected charge-migration channels, synergistically promote this catalytic process. Benefiting from these features, CSS-0.5 delivers superior electrochemical performance under practical conditions, enabling high sulfur loading (6.0 mg cm⁻²) pouch cells with 80.5% capacity retention over 210 cycles. This study provides the first quantitative evidence of electrocatalytic disproportionation in Li–S batteries, offering mechanistic insights and design principles for advanced sulfur cathodes.

25 ENERGY STORAGE

The Significance of Water Vapor Isotopes in Improving Weather Prediction

Water vapor isotopes carry the integrated history of evaporation, condensation, mixing, and transport. Although previous studies have shown potential to improve forecasts under controlled conditions, real-world applications have been limited by systematic biases in both models and satellite retrievals arising from sparse measurements in the free troposphere. Here we assimilate mid‑tropospheric δD retrievals (peak sensitivity ~4.2 km) from the Infrared Atmospheric Sounding Interferometer into the Isotope‑incorporated Global Spectral Model and evaluate the added value beyond co‑assimilated temperature and specific humidity with identical spatial and temporal coverage. Assimilating δD improves 0–120 h forecasts of wind, temperature, specific humidity, and geopotential height, with the largest gains in the midlatitudes; heavy‑precipitation skill also increases for thresholds >3 mm per 6 h. Demonstrated in a coarse‑resolution configuration with limited observations, the results indicate that isotopic information strengthens transport tracking and hydrological constraints, motivating evaluation in operational high‑resolution forecasting systems.

Hydrology