SEARCH · Search NASA

Results for “Time Constant”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 379 records · Page 21

All sources NASA NTRS DOE OSTI

DOE OSTI2026

Fast solvers for tokamak fluid models with PETSc

Multigrid (MG) is widely recognized as a highly effective solver for the model problem, the Laplacian, but textbook MG fails on most problems of interest. MG methods have been applied to complex, real-world applications with careful consideration of the physical model and discretization. In this work we develop the first step in applying MG methods to science and engineering relevant magnetohydrodynamics (MHD) tokamak models in the M3D-C1 (https://m3dc1.pppl.gov) fusion energy science code. The semi-implicit time integrator in M3D-C1 is composed of many linear solves. The implicit advance of the momentum equation is the most challenging and is the focus of this work. The current production solver in M3D-C1 is a block Jacobi (BJ) preconditioner within a Krylov solver, where blocks group degrees of freedom on planes of constant toroidal coordinate. BJ convergence degrades as the number of planes increases due to the spectral properties of the matrix preconditioned with BJ. The partially magnetic field-aligned, regular toroidal grid structure in M3D-C1 is amenable to semi-coarsening geometric MG in the toroidal direction. This paper develops such a solver and demonstrates competitive performance on a runaway electron model of a SPARC (https://cfs.energy/technology/sparc) disruption, and superior robustness on a stellarator model on which the BJ solver fails to converge.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2024

Nonlinear optical properties of lanthanum titanate glasses prepared by levitation melting

Lanthanum titanate glasses (17 La 2 O 3 -83 TiO 2 ) are fabricated via aerodynamic levitation melting. Here, nonlinear refraction and absorption were measured with the Z scan technique using ps optical pulses at 532 nm wavelength with peak intensities in the range of 0.25–2.5 GW/cm 2 . The two-photon absorption coefficient (β) and the effective nonlinear refraction coefficient (⁠$n^{eff}_{2}$) are found to be, respectively, 2.19 cm/GW and 152 × 10 -16 cm 2 /W (85 × 10 -13 esu). The Raman gain is measured to be 106 × 10 -11 cm/W. The nonlinearity strength is found to be nearly 60 times larger in lanthanum titanate glasses relative to silica.

36 MATERIALS SCIENCE↗

DOE OSTI2026

FLASH-therapy suitable single-pulse proton generation using TiH 2 under nanosecond laser irradiation

We investigated proton emission from titanium hydride (TiH 2 ) targets irradiated in vacuum by a 6-ns, 1064-nm Nd:YAG laser. Time-of-flight measurements with a Faraday cup and an electrostatic ion analyzer resolved distinct proton peaks at high pulse energies, with yields on the order of 10 9 protons per shot. Systematic scans over pulse energy and up to 1000 repeated shots showed that, while the peak amplitude gradually decreased, the integrated proton number remained nearly constant. This behavior is consistent with surface-induced broadening of the plasma expansion while bulk hydrogen is replenished by diffusion and repeated ablation of fresh TiH 2 layers. Using the measured proton numbers and pulse widths, a simple direct-plasma-injection-style scaling indicates peak currents of ∼100 mA and sub-microsecond pulse durations. The pulse structure and yield satisfy key ultra-high-dose-rate criteria and, together with sustained output over 1000 shots, support TiH 2 -based laser ion sources as practical candidates for FLASH-therapy (ultra-high-dose-rate)-oriented studies and injector development using direct plasma injection.

43 PARTICLE ACCELERATORS↗

DOE OSTI2024

Asymmetric drift in MaNGA: mass and radially dependent stratification rates in galaxy discs

ABSTRACT We measure the age–velocity relationship from the lag between ionized gas and stellar tangential speeds in ∼500 nearby disc galaxies from MaNGA in Sloan Digital Sky Survey IV (SDSS-IV). Selected galaxies are kinematically axisymmetric. Velocity lags are asymmetric drift, seen in the Milky Way’s (MW) solar neighbourhood and other Local Group galaxies; their amplitude correlates with stellar population age. The trend is qualitatively consistent in rate ($\dot{\sigma }$) with a simple power-law model where σ ∝ tb that explains the dynamical phase-space stratification in the solar neighbourhood. The model is generalized based on disc dynamical times to other radii and other galaxies. We find in-plane radial stratification parameters σ0,r (dispersion of the youngest populations) in the range of 10–40 km s−1 and 0.2 < br < 0.5 for MaNGA galaxies. Overall, brincreases with galaxy mass, decreases with radius for galaxies above 10.4 dex (M⊙) in stellar mass, but is ∼constant with radius at lower mass. The measurement scatter indicates the stratification model is too simple to capture the complexity seen in the data, unsurprising given the many possible astrophysical processes that may lead to stellar population dynamical stratification. None the less, the data show dynamical stratification is broadly present in the galaxy population, with systematic trends in mass and density. The amplitude of the asymmetric drift signal is larger for the MaNGA sample than the MW, and better represented in the mean by what is observed in the discs of M31 and M33. Either typical discs have higher surface-density or, more likely, are dynamically hotter (hence thicker) than the MW.

Bershady, Matthew A. (ORCID:0000000231314374)↗

DOE OSTI2024

Phylometagenomics of cycad coralloid roots reveals shared symbiotic signals

Cycads are known to host symbiotic cyanobacteria, includingNostocalesspecies, as well as other sympatric bacterial taxa within their specialized coralloid roots. Yet, it is unknown if these bacteria share a phylogenetic origin and/or common genomic functions that allow them to engage in facultative symbiosis with cycad roots. To address this, we obtained metagenomic sequences from 39 coralloid roots sampled from diverse cycad species and origins in Australia and Mexico. Culture-independent shotgun metagenomic sequencing was used to validate sub-community co-cultures as an efficient approach for functional and taxonomic analysis. Our metanalysis shows a host-independent microbiome core consisting of seven bacterial orders with high species diversity within the identified taxa. Moreover, we recovered 43 cyanobacterial metagenome-assembled genomes, and in addition toNostocspp., symbiotic cyanobacteria of the genusAulosirawere identified for the first time. Using this robust dataset, we used phylometagenomic analysis to reveal three monophyletic cyanobiont clades, two host-generalist and one cycad-specific that includesAulosiraspp. Although the symbiotic clades have independently arisen, they are enriched in certain functional genes, such as those related to secondary metabolism. Furthermore, the taxonomic composition of associated sympatric bacterial taxa remained constant. Our research quadruples the number of cycad cyanobiont genomes and provides a robust framework to decipher cyanobacterial symbioses, with the potential of improving our understanding of symbiotic communities. This study lays a solid foundation to harness cyanobionts for agriculture and bioprospection, and assist in conservation of critically endangered cycads.

Genetics & Heredity↗

DOE OSTI2025

In-Situ Laser Diagnostics to Understand Plasma-Surface Interactions in Titanium Thin Film Deposition

During the current project period the following tasks were completed: (1) We demonstrated the feasibility of a new way of measuring CH3 radicals in a plasma environment using a single, femtosecond (fs) pulsed laser for the first time, in a CH4 + Ar flow. A single fs pulsed Ti:Sa laser was split into two beams and with beam conversion using crystals we generatedtwo beams of 275 nm and 205 nm. The 275 nm was used to predissociate the CH3 radicals to CH2+H, followed by 205 nm Two Photon Laser Induced Fluorescence (TALIF) detection of H atoms, with a 11 ns delay from 275 nm predissociation pulse. (2) We demonstrated preliminary CH3 concentration measurements with spatial (< 100 µm) and time resolution (< 10 ns) by image processing of H and Kr TALIF signal images. (3) The methodology to quantify CH3 concentrations was developed. (4) Preliminary measurements of laser dissociation cross sections for CH4 and CH3 to (CH3+H), and (CH2+H), respectively, were conducted. (5) 275 nm laser dissociation of CH4 is a multiphoton process, and the CH4 dissociation cross section is a more sensitive function of the laser pulse energy than CH4 concentration. In the coming months, we will make more precise measurements of the dissociation cross sections, 𝑘275, 𝐶𝐻4 and 𝑘275, 𝐶𝐻3, by varying the laser energy while keeping the flow composition and conditions constant. This will enable more precise quantification of CH3 radical concentrations. More detailed measurements of CH3 radicals will be conducted in the following months of 2025 (Sept-Dec) in a Titanium Isopropoxide (TIP) + H2 + Ar flow meant for pure phase Ti thin film deposition by Chemical Vapor Deposition (CVD). More precise measurements of 275 nm laser dissociation cross sections for CH4 (to CH3 + H) and CH3 (to CH2 + H) will be conducted in the future. Detailed 275 nm laser predissociation cross sections for CH4 and CH3 dissociation to (CH3+H), and (CH2+H), respectively, as a function of gas heating, plasma power, and laser pulse energy will be evaluated. This will enable quantification of CH3 concentrations in the precursor flow over a substrate with time and spatial resolution. Comparison of the concentrations with 2D numerical simulations can lead to understanding of plasma surface interactions during pure metallic Ti thin film deposition using environmentally safer non-halogen TIP precursor. Flow conditions for pure Ti thin film deposition using non-halogen precursor (TIP) will be identified in the following months.

Uddi, Mruthunjaya [Advanced Cooling Technologies, ↗

DOE OSTI2025

Revisiting the Role of Mobile Ion Size for the Activation Barrier of Charge Transport in Single Ion Conducting Polymers

The development of polymer electrolytes for energy storage and conversion technologies requires a fundamental understanding of the material parameters controlling the energy barriers for ion transport. In glassy polymers, these activation barriers are usually extracted by using Arrhenius procedures. However, our recent studies of single ion conducting polymers reveal that this traditional Arrhenius description provides anomalously small prefactors, an issue that is widely disregarded in the literature. Here, to address it, an alternative approach was introduced for extracting the effective activation barrier E by imposing a physically valid limit (∼10 –13 s) for the ion hopping attempt time. Two consequences of this recent approach are that in polymer electrolytes (i) E is significantly (∼30–40%) lower than the barrier estimated using traditional Arrhenius fits and (ii) E displays significant temperature dependence even in the glassy state. Under these circumstances, we are revisiting the role of ion size, glass transition temperature T g , and dielectric and elastic constants of the polymer matrix on the effective activation barriers for ion transport as extracted using the recently proposed approach and highlight the differences from the picture that previously emerged based on simplistic apparent Arrhenius analysis. To this end, we investigate cation transport in three families of single ion conducting trifluoromethane sulfonimide-based polyanions with varied T g . Our results indicate that E decreases with mobile cation size and is highly sensitive to changes in the dielectric permittivity of the matrix, even for large cations. These insights call for revisions of many earlier results based on apparent Arrhenius fits as the proposed approach can provide more accurate guidance for the design of polymer electrolytes with enhanced ionic conductivity.

energy barrier↗

DOE OSTI2025

Hydroxy-substituted electron deficient Pd porphyrin cofactors illuminate ultrafast proton transfer reactions

Coordinated electron and proton movement drive bioenergetic functions. Relative to electron transfer reactions, tracking proton transport over fast-to-ultrafast time scales is challenging. Optical resolution of proton transfer dynamics can take advantage of chromophoric photoacids that not only trigger proton migration upon photoexcitation, but produce distinct spectroscopic changes associated with protonation/deprotonation. Here, in this work, we report the design of a hydroxy-substituted electron deficient Pd porphyrin, PPd(C 6 F 5 ) 3 OH; upon photoexcitation, the acidity constant of this weak acid (pKa = 6.49) dramatically drops (pKa* = 0.97). Electronic excitation of PPd(C 6 F 5 ) 3 OH triggers an ultrafast proton transfer reaction (PPd(C 6 F 5 ) 3 OH +:B + hυ → 1[PPd(C 6 F 5 ) 3 OH]* +:B → 1 [PPd(C 6 F 5 ) 3 O]−* + HB + ; τ PT = 342 fs) to a H-bonded base (B) in solution. Both PPd(C 6 F 5 ) 3 OH and its conjugate anion PPd(C 6 F 5 ) 3 O − exhibit distinct vis-NIR spectral features for their respective ground and excited states. Because the electroni-cally excited triplet lifetimes of these species exceed tens of microseconds, the PPd(C 6 F 5 ) 3 OH photoacid defines an ideal cofactor to probe light-triggered proton release and track long-range proton migration in protein environments.

Palladium↗

DOE OSTI2025

Xerogel-Derived Ni Electrocatalysts for the Hydrogen Evolution Reaction in Alkaline Media

Anion exchange membrane water electrolyzers (AEMWEs) represent a promising technology for hydrogen production. The big advantage of the technology is that it allows for the use of platinum group metal-free (PGM-free) electrocatalysts at both electrodes, including catalysts for the hydrogen evolution reaction (HER) at the cathode. In addition to fulfilling the cost requirement, PGM-free HER catalysts need to meet the activity and durability targets of the AEMWEs. Here, in this work, we developed several carbon-supported, xerogel-derived nickel (Ni) HER electrocatalysts and evaluated the effect of various synthesis conditions, such as the type of carbon support, Ni-to-carbon ratio, and heat-treatment temperature and time, on their performance. Scanning transmission electron microscopy combined with energy-dispersive X-ray spectroscopy (STEM-EDS), X-ray diffraction spectroscopy (XRD), and X-ray photoelectron spectroscopy (XPS) revealed the formation of Ni nanoparticles with an oxygen-rich layer on the outside. Durability of the best-performing catalyst was assessed via a constant-current hold at 10 mA cm –2 over 100 h. This catalyst was found to be more active and durable than the reference PGM-free material, a commercial Ni catalyst supported on a Vulcan XC-72. The catalyst was also tested in the cathode of a fully PGM-free AEMWE, allowing to reach 1.90 V (1.84 V HFR-free) at 1 A cm –2 at 80 °C.

08 HYDROGEN↗

DOE OSTI2025

Catalytic Resonance Theory: Experimental and Kinetic Interpretation of Programmable Catalysis

The complexity of programmable catalysts that change in the physical or electronic state on the time scale of a catalytic turnover results in uncertainty in the kinetic outcome of common measurements with dynamic catalyst experiments. To provide an interpretation of utility to experimentalists, a general programmable catalytic reaction exhibiting high turnover efficiency was simulated to understand the kinetic behavior of catalysts with varying frequency, binding energy, oscillation amplitude, and minimum oscillation binding energy at different temperatures and reactant gas pressures. The time-averaged catalytic rate for varying temperature in the form of an Arrhenius plot exhibited three distinct kinetic regimes, with the intermediate temperatures or applied frequencies exhibiting zero slope, indicative of barrier-free kinetic control. Transitions in experimentally measurable kinetic regimes with both temperature and applied catalyst oscillation frequency were associated with transitions in the sensitivity of reaction rate constants and degrees of rate control. Furthermore, these distinct kinetic macroscopic features specific to programmable catalysts were identified for observation in experimental dynamic kinetic measurements.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2025

Metabolic Modification of Sphingobium lignivorans SYK-6 for Lignin Valorization Through the Discovery of an Unusual Transcriptional Repressor of Lignin-Derived Dimer Catabolism

Sphingobium lignivorans SYK-6 catabolizes guaiacylglycerol-..beta..-guaiacyl ether (GGE, a ..beta..-O-4-type dimer) and 1,2-diguaiacylpropane-1,3-diol (DGPD, a ..beta..-1-type dimer) derived from lignin. Recently, SLG_35860 containing TetR- and MarR-type transcriptional regulator motifs was suggested to be involved in the regulation of GGE and DGPD catabolism. Here we investigated the role of SLG_35860 in the transcriptional regulation of GGE and DGPD catabolism genes. SLG_35860 designated ligS repressed 11 genes involved in GGE and DGPD catabolism. LigS binds directly to specific sequences in the promoter region of each gene. The MarR domain was shown to be involved in these bindings; however, GGE, DGPD, and their metabolites did not function as effectors of LigS. We discovered unidentified compound(s) in the black liquor of oxygen-soda anthraquinone pulping of Japanese cedar that SYK-6 cannot metabolize and that acted as effector(s). Therefore, LigS constantly represses the transcription of the GGE and DGPD catabolism genes to low levels. Based on these findings, we examined the productivity of a polymer building block, 2-pyrone-4,6-dicarboxylic acid (PDC), from GGE, DGPD, and a GGE metabolite using an engineered ligS mutant. The rates of PDC production from each compound by this strain were 1.5-6.0 times higher than those of a PDC-producing strain carrying ligS.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2025

Novel Technology of Non-Contact Real-Time Radiation Damage Sensor for High Power Targets

High-power proton beams planned for forthcoming long-baseline neutrino experiments will subject solid targets to unprecedented radiation damage, threatening reliability and increasing costs. To address this challenge, we are building a real-time, non-contact sensor that tracks damage by measuring broadband laser reflectivity changes from the target surface. A low-power super-continuum Class 3B laser illuminates the sample inside a vacuum test chamber while a high-resolution fiber-coupled spectrometer captures S- and P-polarized light; spectral shifts reveal defect-driven variations in optical constants. My internship goal is to design, build, and commission this prototype by implementing a laser-safety interlock and light shield, integrating remote-operation cameras, and fabricating modular 3-D-printed mounts that enable tool-free swaps without disturbing alignment. The laser, spectrometer, and vacuum chamber have been delivered; interlock hardware, cameras, and mounts are in final assembly, and leak testing of the chamber is underway. Upcoming work will focus on initial calibration of the sensor and executing first beam-irradiation studies.

Pumarino, Rafael↗

DOE OSTI2026

Characterization of Spurious Charge in SENSEI Skipper-CCDs

Skipper Charge-Coupled Devices (Skipper-CCDs) are a leading technology in the search for sub-GeV dark matter and coherent elastic neutrino-nucleus scattering. A key background for rare-event searches with these detectors arises from "spurious charge" -- single-electron events generated when charges are transferred through the active region to the serial register, and across the serial register to the readout stage. We present a characterization of spurious charge in both the active region and the serial register of SENSEI Skipper-CCDs, and show that, in a well-shielded low-background environment, the dominant contribution originates in the serial register during Skipper readout, when horizontal clocks are held at constant voltage between pixel transfers. Motivated by this finding, we develop a "tri-level" clocking scheme in which the held-low phase is raised to an intermediate voltage during readout to suppress trap-mediated charge generation. Using the SENSEI detector near the MINOS cavern, we measure a serial-register single-electron density of $(2.9 \pm 0.1) \times 10^{-5}$ electrons/pixel/image under standard SENSEI readout conditions, reduced to $(4.0 \pm 0.4) \times 10^{-6}$ electrons/pixel/image with tri-level clocking -- a factor of $\sim$7 improvement. This technique offers a promising path to lower backgrounds in current and future Skipper-CCD experiments.

Wu, Yikai [YITP, Stony Brook; SUNY, Stony Brook]↗

DOE OSTI2025

Randomized Adiabatic Quantum Linear Solver Algorithm with Optimal Complexity Scaling and Detailed Running Costs

Solving linear systems of equations is a fundamental problem with a wide variety of applications across many fields of science, and there is increasing effort to develop quantum linear solver algorithms. Subaşı et al. [Phys. Rev. Lett. 122, 060504 (2019)] proposed a randomized algorithm inspired by adiabatic quantum computing, based on a sequence of random Hamiltonian simulation steps, with suboptimal scaling in the condition number 𝜅 of the linear system and the target error 𝜖. Here we go beyond these results in several ways. Firstly, using filtering [Lin and Tong, Quantum 4, 361 (2020)] and Poissonization techniques [Cunningham and Roland, ArXiv:2406.03972 (2024)], the algorithm complexity is improved to the optimal scaling 𝑂⁡(𝜅⁢log (1/𝜖))—an exponential improvement in 𝜖, and a shaving of a log 𝜅 scaling factor in 𝜅. Secondly, the algorithm is further modified to achieve constant factor improvements, which are vital as we progress towards hardware implementations on fault-tolerant devices. We introduce a cheaper randomized walk operator method replacing Hamiltonian simulation—which also removes the need for potentially challenging classical precomputations; randomized routines are sampled over optimized random variables; circuit constructions are improved. We obtain a closed formula rigorously upper bounding the expected number of times one needs to apply a block-encoding of the linear system matrix to output a quantum state encoding the solution to the linear system. The upper bound is 837⁢𝜅 at 𝜖 = 10 −10 for Hermitian matrices.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2025

A Dynamic Pricing Method to Manage the Impact of EV Charging on the Grid Using RL

This work addresses the challenge of managing electrical vehicle (EV) charging loads on distribution feeders with the increase in deployment of fast charging stations. To mitigate the adverse impacts on feeder health, a novel dynamic grid-informed pricing approach is proposed. This approach leverages reinforcement learning (RL) to determine hourly charging prices based on real-time grid conditions. A synthetic environment was developed to train the reinforcement learning agent. A model of an IEEE 34-bus distribution feeder with EV charging stations has been developed in OpenDSS utilizing Caldera for realistic EV charging profiles. Test cases demonstrate that the dynamic pricing strategy achieves higher energy delivery to the EV end user at a lower cost compared to constant pricing methods, while lowering voltage deviations and congestion. This approach offers more granular price adjustments, responding dynamically to feeder conditions and potentially improving grid stability and efficiency. The communication architecture to implement this dynamic pricing method is described. This research contributes to the development of smart grid-informed charging solutions that can reduce the cost of charging to the end user while also helping the grid.

EV charging, dynamic pricing, grid-informed chargi↗

DOE OSTI2025

A dynamic pricing method to manage the impact of EV charging on the grid using RL

24 - POWER TRANSMISSION AND DISTRIBUTION↗

DOE OSTI2024

Incorporation of rare earth elements $\text{Sc}$, $\text{Y}$ and $\text{La}$ into gibbsite

Rare earth elements (REE) are increasingly critical resources in technological applications, but the current understanding of their separation by uptake into minerals remains limited. We examined the adsorption and solubilization of scandium (Sc), yttrium (Y), and lanthanum (La) into aluminum hydroxide (gibbsite, Al(OH) 3 ). Based on spectroscopic investigations such as inductively coupled plasma-optical emission spectroscopy (ICP-OES) and time of flight-secondary ion mass spectrometry (ToF-SIMS), only Sc exhibited favorable incorporation and adsorption (0.66 and 0.06 atom%, respectively). Analysis of solid state 45 Sc magic angle spinning-nuclear magnetic resonance spectroscopy (MAS-NMR) indicated that Sc exhibits an octahedral coordination consistent with isolated substitutions for Al in the gibbsite structure, leading to a increase in lattice constants detectable by X-ray diffraction. Density functional theory (DFT) calculations reinforced this interpretation by reproducing the detected structural distortion as well as demonstrating the relatively favorable energetic basis for Sc incorporation into gibbsite by substitution.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2025

Longitudinal Patterns in Nitrate Revealed Through Drone‐Based Measurements in an Agriculturally Influenced Midwestern U.S. River

Nitrate concentrations in streams and rivers in the Midwestern United States are often elevated, reflecting the predominance of agriculture in the surrounding landscape. Recent advances in technology, including surface water drones and more precise sensors, provide opportunities to investigate nitrate dynamics with high spatial and temporal resolution. We deployed an aquatic drone, the AquaBOT, in a sixth-order, agriculturally influenced river to examine longitudinal patterns in water quality. Our goal was to measure the spatial and temporal heterogeneity in nitrate and nitrate removal processes and determine the influence of tributary inputs on main stem chemistry. We navigated the drone along a 12-km reach of the Des Moines River (Iowa, USA) nine times between June 2021 and August 2022. Across the deployments, mean nitrate concentration was positively related to discharge and was nearly two orders of magnitude higher in spring than summer. We observed contrasting patterns in main stem nitrate, which decreased downstream during some runs (e.g., 3.1–2.7 mg N L −1 in June 2021), demonstrating net nitrate uptake along the reach, and remained constant on other dates. Similarly, tributaries to the Des Moines had a varied influence on riverine nitrate. Tributaries either increased or decreased main stem nitrate concentrations depending on the tributary and the date. Nitrate removal rates were spatially and temporally variable but showed some consistency at the subreach (2 km) scale, with two subreaches often showing elevated rates of nitrate removal across dates. In conclusion, our study reveals nuanced heterogeneity in nitrate dynamics of the Des Moines River despite the homogeneity of agricultural land cover in the watershed.

54 ENVIRONMENTAL SCIENCES↗