SEARCH · Search NASA

Results for “forcing”

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 145 records · Page 8

All sources NASA NTRS DOE OSTI

DOE OSTI2026

Investigation of the Effect of Framework Flexibility on Adsorption in SIFSIX-3-Cu using a Machine-Learned Force Field

Metal-organic frameworks (MOFs) are a promising class of adsorbents. The performance of MOF sorbents relies on high selectivity and low regeneration energy. This work focuses on the use of machine learned force fields (MLFFs) to model adsorption in a flexible MOF, SIFSIX-3-Cu. A DeePMD-based MLFF was trained to reproduce DFT (PBE+D3) energies, forces, and stresses, using an iterative sampling scheme combining sampling based on molecular dynamics, Monte Carlo, and geometry optimization to capture both attractive and repulsive regions of the potential energy surface. Flexibility of the MOF was explicitly included in this model. Hybrid Monte Carlo/molecular dynamics (MC/MD) simulations using the MLFF predicted adsorption isotherms in good agreement with experimental data for a range of pressures (40 Pa – 104 Pa) in contrast to rigid models, which overpredict CO2 adsorption at low pressures. The improvement was the result of a description of the variability of fluorine-fluorine diagonal distances at adsorption sites. This detailed description of flexibility afforded by the MLFF resulted in more accurate predictions adsorption isotherms when compared to the experimentally measured values. These results underscore the importance of including framework flexibility when modeling adsorption phenomena in MOFs, particularly for low pressure applications and provide a robust procedure for training MLFF models for MOFs.

Atomistic Simulation↗

DOE OSTI2026

Mapping Interfacial Solution Forces that Drive Fluorescent Nanoparticle Incorporation into Crystals

Interfacial solution structure governs processes ranging from catalytic and electrochemical reactions to particle aggregation and crystallization. Here, we design bio-inspired nanocomposites by mapping interfacial forces and controlling nanoparticle incorporation during crystal growth. Our analytical model contains measured kinetic barriers for surface approaches and equilibrium binding constants with the crystal surface. We validate this model using fluorescent silica nanoparticles and calcite. Our results show that surface chemistry dictates incorporation pathways: methoxy groups are least favorable, hydroxyls reduce kinetic barriers by interacting favorably with hydration layers, carboxylates bind strongly but must overcome a kinetic barrier, and amines outperform all kinetically and thermodynamically. This framework resolves the interplay between interfacial solution structures, particle forces and dynamics, and growth kinetics, enabling predicative control of nanocomposite formation and multiplexed mixed-particle systems.

Zhang, Mingyi [Pacific Northwest National Laborato↗

DOE OSTI2025

A new method for diagnosing effective radiative forcing from aerosol–cloud interactions in climate models

Aerosol–cloud interactions (ACIs) are a leading source of uncertainty in estimates of the historical effective radiative forcing (ERF). One reason for this uncertainty is the difficulty in estimating the ERF from aerosol–cloud interactions (ERFaci) in climate models, which typically requires multiple calls to the radiation code. Most commonly used methods also cannot disentangle the contributions from different processes to ERFaci. Here, we develop a new, computationally efficient method for estimating the shortwave (SW) ERFaci from liquid clouds using histograms of monthly averaged cloud fraction partitioned by cloud droplet effective radius (r e ) and liquid water path (LWP). Multiplying the histograms with SW cloud radiative kernels gives the total SW ERFaci from liquid clouds, which can be decomposed into contributions from the Twomey effect, LWP adjustments, and cloud fraction (CF) adjustments. We test the method with data from five CMIP6-era models, using the Moderate Resolution Imaging Spectroradiometer (MODIS) satellite instrument simulator to generate the histograms. Our method gives similar total SW ERFaci estimates to other established methods in regions of prevalent liquid cloud and indicates that the Twomey effect, LWP adjustments, and CF adjustments have contributed -0.34 ± 0.23, -0.22 ± 0.13, and -0.09 ± 0.11 W m -2 , respectively, to the effective radiative forcing of the climate since 1850 in the ensemble mean (95 % confidence). These results demonstrate that widespread adoption of a MODIS r e –LWP joint histogram diagnostic would allow the SW ERFaci and its components to be quickly and accurately diagnosed from climate model outputs, a crucial step for reducing uncertainty in the historical ERF.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2025

Development of a forced advection sampling technique (FAST) for quantification of methane emissions from orphaned wells

Abstract. Orphaned wells, meaning unplugged and non-producing wells lacking responsible owners, pose a significant and undersampled environmental challenge due to their vast number and unknown associated emissions. We propose, develop and test an alternative method for estimating emissions from orphaned wells using a forced advection sampling technique (FAST) that can overcome many of the limitations in current methods (cost, accuracy, safety). In contrast to existing ambient Gaussian plume methods, our approach uses a fan-generated flow to force advection between the emission source and a point methane (CH4) sensor. The fan flow field is characterized using a colocated sonic anemometer to measure the 3D wind profile generated by the fan. Using time-series measurements of CH4 concentration and wind, a simple estimate of the CH4 emission rate of the source can be inferred. The method was calibrated using outdoor controlled-release experiments and then tested on four orphaned wells in Lufkin, TX, and Osage County, OK. Our results suggest that the FAST method can provide a low-cost, portable, fast and safe alternative to existing methods with reasonable estimates of orphaned well emissions over a range of leak rates below 40 g h−1 and within certain geometric and atmospheric constraints.

Dubey, Mohit L↗

DOE OSTI2025

Legacy Survey of Space and Time Data Preview 1: dia_object_forced_source dataset type

The Legacy Survey of Space and Time Data Preview 1 (DP1) is the first release of data from the NSF-DOE Vera C. Rubin Observatory. It consists of raw and calibrated single-epoch images, co-adds, difference images, detection catalogs, and other derived data products. DP1 is based on 1792 science-grade optical/near-infrared exposures acquired over 48 distinct nights by the Rubin Commissioning Camera, LSSTComCam, on the Simonyi Survey Telescope at the Summit Facility on Cerro Pachón, Chile during the first on-sky commissioning campaign in late 2024. DP1 covers a total of approximately 15 sq. deg. over seven roughly equally-sized non-contiguous fields, each independently observed in six broad photometric bands, ugrizy, spanning a range of stellar densities and latitudes and overlapping with external reference datasets. This dataset is a subset of the full data release consisting of the dia_object_forced_source dataset type. These are forced measurements in visit and difference images, at the coordinates of all DIA objects. This release contains 597 datasets of this type.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2025

Legacy Survey of Space and Time Data Preview 1: object_forced_source dataset type

The Legacy Survey of Space and Time Data Preview 1 (DP1) is the first release of data from the NSF-DOE Vera C. Rubin Observatory. It consists of raw and calibrated single-epoch images, co-adds, difference images, detection catalogs, and other derived data products. DP1 is based on 1792 science-grade optical/near-infrared exposures acquired over 48 distinct nights by the Rubin Commissioning Camera, LSSTComCam, on the Simonyi Survey Telescope at the Summit Facility on Cerro Pachón, Chile during the first on-sky commissioning campaign in late 2024. DP1 covers a total of approximately 15 sq. deg. over seven roughly equally-sized non-contiguous fields, each independently observed in six broad photometric bands, ugrizy, spanning a range of stellar densities and latitudes and overlapping with external reference datasets. This dataset is a subset of the full data release consisting of the object_forced_source dataset type. These are forced measurements in visit and difference images, at the coordinates of all objects. This release contains 636 datasets of this type.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2026

Legacy Survey of Space and Time Data Preview 2: dia_object_forced_source dataset type

We present Rubin Data Preview 2 (DP2), the second data preview from the NDF-DOE Vera C. Rubin Observatory. Data Preview 2 (DP2) comprises coadds, detection catalogs, and ancillary data products; and when fully released will also include single-epoch images and difference images. DP2 is derived from observations acquired by the LSST Science Camera (LSSTCam) on the Simonyi Survey Telescope at the Summit Facility on Cerro Pachón, Chile, primarily during the on-sky commissioning campaign between 2025-04-16 and 2025-09-21, supplemented by observations taken between 2025-10-25 and 2026-01-06 that overlap the commissioning footprint. The DP2 footprint comprises the Science Validation wide-area survey, five Deep Drilling Fields, and a number of targeted small-field regions, including Trifid and Lagoon, Prawn, M49, and New Horizons, all observed as part of the Rubin First Look campaign. Each field was imaged in up to six broad photometric bands, ugrizy, and coadded to produce deep imaging covering an estimated 3,000 deg2. The addition of single-visit-only areas expands the total DP2 footprint to an estimated 15,000 deg2, with coverage in at least one filter. The median per-visit PSF FWHM across the wide-area survey ranges from 1.17 arcsec in the z band to 1.26 arcsec in g and r bands. The deepest field, reaches estimated coadded 5σ depths of u=26 mag, g=26.8 mag, r=26.3 mag, i=26.1 mag, z=25.3 mag, y=23.9 mag. Based on a roughly five-month primary observing baseline and covering only part of the eventual LSST footprint, DP2's area, depth, and multiband coverage nonetheless support a broad range of early science investigations ahead of LSST Data Release This dataset is a subset of the full data release consisting of the dia_object_forced_source dataset type. These are forced measurements in visit and difference images, at the coordinates of all DIA objects. This release contains 183,169 datasets of this type.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2026

Legacy Survey of Space and Time Data Preview 2: object_forced_source dataset type

We present Rubin Data Preview 2 (DP2), the second data preview from the NDF-DOE Vera C. Rubin Observatory. Data Preview 2 (DP2) comprises coadds, detection catalogs, and ancillary data products; and when fully released will also include single-epoch images and difference images. DP2 is derived from observations acquired by the LSST Science Camera (LSSTCam) on the Simonyi Survey Telescope at the Summit Facility on Cerro Pachón, Chile, primarily during the on-sky commissioning campaign between 2025-04-16 and 2025-09-21, supplemented by observations taken between 2025-10-25 and 2026-01-06 that overlap the commissioning footprint. The DP2 footprint comprises the Science Validation wide-area survey, five Deep Drilling Fields, and a number of targeted small-field regions, including Trifid and Lagoon, Prawn, M49, and New Horizons, all observed as part of the Rubin First Look campaign. Each field was imaged in up to six broad photometric bands, ugrizy, and coadded to produce deep imaging covering an estimated 3,000 deg2. The addition of single-visit-only areas expands the total DP2 footprint to an estimated 15,000 deg2, with coverage in at least one filter. The median per-visit PSF FWHM across the wide-area survey ranges from 1.17 arcsec in the z band to 1.26 arcsec in g and r bands. The deepest field, reaches estimated coadded 5σ depths of u=26 mag, g=26.8 mag, r=26.3 mag, i=26.1 mag, z=25.3 mag, y=23.9 mag. Based on a roughly five-month primary observing baseline and covering only part of the eventual LSST footprint, DP2's area, depth, and multiband coverage nonetheless support a broad range of early science investigations ahead of LSST Data Release This dataset is a subset of the full data release consisting of the object_forced_source dataset type. These are forced measurements in visit and difference images, at the coordinates of all objects. This release contains 192,147 datasets of this type.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2026

Force-triggered, Bio-based, Sealants for Prefabricated Building Components: Towards Improved Efficiency, Performance and Sustainability

The prefabricated building construction industry has made extensive progress in expediting the manufacture of prefabricated components at offsite plants. However, this progress has not translated to the assembly of the prefabricated components at the construction site. Case in point, sealing the joints between components to prevent air leaks requires the manual application of tape, caulk, or spray foam at the jobsite, and performance is highly dependent on the skills of the installer. To reduce assembly time and improve assembly quality of prefabricated components, we developed a sealant that can be installed at the plant and triggered at the jobsite. Additionally, we used this opportunity to lower the use of fossil fuel derived feedstocks and introduced bio-based alternatives to decrease the embodied carbon of the new sealant. We evaluated a series of force-triggered, bio-based, high strength, and fast curing sealants, consisting of a one-part heterogeneous system. These sealants are derived from formulations with ≥80% of biogenic carbon, consisting of a cardanol derived diepoxy that is microencapsulated in a polymer shell and embedded in a cardanol derived amine curing agent. The microcapsule shell allows separation of the reactive species in the one part sealant allowing shelf stability to an otherwise fast curing system. When the microcapsules are broken and activated by force, the highly reactive species mix and cure, exhibiting peel strengths up to 143 ppi. The shelf stability of these sealants and the on-demand triggering of the curing reactions enable installation on prefabricated components and storage prior to delivery and assembly at a jobsite, which could result in consistent sealant application with lower installation time than tapes and caulks at the construction site.

Cortes Guzman, Karen [ORNL] (ORCID:000000028793468↗

DOE OSTI2024

Unveiling the nanoscale architectures and dynamics of protein assembly with in situ atomic force microscopy

Proteins play a vital role in different biological processes by forming complexes through precise folding with exclusive inter- and intra-molecular interactions. Understanding the structural and regulatory mechanisms underlying protein complex formation provides insights into biophysical processes. Furthermore, the principle of protein assembly gives guidelines for new biomimetic materials with potential applications in medicine, energy, and nanotechnology. Atomic force microscopy (AFM) is a powerful tool for investigating protein assembly and interactions across spatial scales (single molecules to cells) and temporal scales (milliseconds to days). It has significantly contributed to understanding nanoscale architectures, inter- and intra-molecular interactions, and regulatory elements that determine protein structures, assemblies, and functions. This review describes recent advancements in elucidating protein assemblies with in situ AFM. We discuss the structures, diffusions, interactions, and assembly dynamics of proteins captured by conventional and high-speed AFM in near-native environments and recent AFM developments in the multimodal high-resolution imaging, bimodal imaging, live cell imaging, and machine-learning-enhanced data analysis. These approaches show the significance of broadening the horizons of AFM and enable unprecedented explorations of protein assembly for biomaterial design and biomedical research.

36 MATERIALS SCIENCE↗

DOE OSTI2025

Global warming potential estimates of mass timber constructions beyond the first life: A dynamic radiative forcing modeling approach

In the face of a warming planet, steps must be taken to reduce the greenhouse gas emissions (GHG) associated with our building industry, which is a significant contributor to global emissions. Large, prefabricated wood elements such as mass timber panels (MTP) have great potential to achieve these reductions as they help displace high-embodied‑carbon materials like concrete and steel. Furthermore, storing the wood's biogenic carbon in buildings benefits the climate because it delays the eventual release of the carbon into the atmosphere. While these climate impacts have been assessed for the construction phase of mass timber buildings, relatively few life cycle assessment (LCA) studies have evaluated the climate impacts for the buildings' end-of-life (EOL) phase. This research estimates the climate impacts of four EOL scenarios for MTP: reusing as MTP, recycling into particleboard, incinerating, and landfilling. Using dynamic radiative forcing modeling and factoring in temporal GHG emissions and biogenic carbon storage, the global warming potential impacts are calculated for construction, deconstruction, and EOL processing of hybrid mass timber buildings in the U.S. Pacific Northwest for 160 years (GWP 160 ). The 160-year temporal scale used in this paper is an arbitrary scale, with the first 80 years being the assumed life of the building, followed by a series of reuse, recycle, or disposal scenarios over the second 80 years of that temporal scale. Of the four EOL scenarios considered in this paper, the ‘reuse’ scenario has the lowest net GWP 160 impact (calculated by summing the GWP 160 and carbon storage benefits, i.e., GWP$^{bioCS}_{160}$ ), emerging as a climate-preferred scenario, followed by ‘landfill’, ‘incinerate’, and ‘recycle’ scenarios. The lower net GWP 160 impact associated with the reuse scenario is due to the low fossil carbon emissions during EOL processing, as well as the biogenic carbon storage benefits. The results of this study also highlight the importance of efficient reuse and recycling strategies for wood in MTP.

42 ENGINEERING↗

DOE OSTI2024

Atomic force microscopy monitoring of subcritical slow crack growth in soda lime silicate glass

Prediction of brittle fracture in silicate glasses requires a fundamental understanding of slow crack growth rates, far below the critical stress intensity (K IC ) of the material. Here, in this study, atomic force microscopy (AFM) was used to characterize slow crack propagation in soda lime silicate glasses. Using large-sample AFM systems with specialized tips and scanners, fractures less than 20 nm wide were measured as a function of relative humidity (RH). The results were found to be sensitive to the type of AFM tip used and to the dryness of the surface. Continuous scanning and image capture over several days allowed the determination of the crack tip velocity. The dependance of crack tip velocity on RH was measured by enclosing the sample and AFM head in a humidity-controlled chamber. The crack tip velocity was measured as low 4.0×10 -12 m/s at 15 % RH and as high as 3.7×10 -10 m/s at 40 % RH. Overall crack growth was smooth during the AFM scan. Additionally, recession of the crack tip was noted as the sample dried, indicating crack tip closure may be occurring. Further studies on alternative compositions will provide details on subcritical fracture processes in a range of different environmental conditions.

Atomic force microscopy↗

DOE OSTI2025

High-fidelity forced convection simulations of the University of Wisconsin–Madison air-cooled reactor cavity cooling system

Among the various types of advanced nuclear technologies that exist, high-temperature gas-cooled reactors (HTGRs) stand out for their co-generation capabilities and exceptional passive safety systems. HTGRs rely on a reactor cavity cooling system (RCCS) to remove decay heat from the reactor pressure vessel (RPV) during the course of long-term transients. Several studies have focused on experimental and numerical modeling of the RCCS; however, the literature reveals a notable gap when it comes to high-fidelity simulations such as direct numerical simulations (DNS) or large eddy simulation (LES) models of the RCCS. An LES model was developed for forced convection inside the RCCS so as to obtain a high-fidelity representation of the RCCS at a scale never before attempted. Furthermore, the LES model was compared against three different Reynolds-averaged Navier–Stokes (RANS) models (Realizable k − ϵ, k − ω SST, and Reynolds stress model [RSM]). This comparison verified the robustness and accuracy of the turbulence modeling approach of the RANS models for the RCCS. Finally, based on the comparison between the LES and the RANS models, the k − ω SST RANS model of the experimental facility was developed and compared against the available experimental data.

22 - GENERAL STUDIES OF NUCLEAR REACTORS↗

DOE OSTI2024

Characterization of Polyamide Thin Films by Atomic Force Microscopy

This study directly compares the mechanical behavior of novel molecular layer deposition (MLD) and analogous interfacial polymerization (IP) polyamide thin films in environments relevant to reverse osmosis (RO) membrane operation. The elastic modulus of the films was determined using atomic force microscopy (AFM) in dry, hydrated, and chlorinated states. Surface roughness characteristics were also obtained given their potential influence on AFM modulus measurements. The much smoother MLD films demonstrated a statistically higher modulus in all states as compared to their IP counterparts. The MLD films maintained a modulus ~3X and ~5X greater than that of IP films after hydration and chlorination, respectively. Such differences in behavior may be due to the higher density and correspondingly lower void content of the MLD films. Results from this study provide a rationale for future development of MLD for fabrication of polyamide films for incorporation in RO membranes.

atomic force microscopy↗

DOE OSTI2026

A Transferable Force Field for Simulating Adsorption in Metal–Organic Frameworks with Open Metal Sites Based on the 12–6–4 Lennard-Jones Potential

Metal−organic frameworks (MOFs) that contain coordinatively unsaturated open metal sites (OMSs) provide strong host− guest interactions, making them promising sorbents for low-concentration gas adsorption applications such as direct air capture and atmospheric water harvesting. However, accurately modeling host−guest interactions involving OMSs remains challenging for classical force fields (FFs) based on the 12−6 Lennard−Jones (LJ) potential, as the polarization effect of the guest molecule induced by the positively charged OMS is not considered. Here, we introduce an FF based on the 12−6−4 LJ potential, which incorporates charge−induced dipole interactions and is parametrized against a diverse set of host−guest potential energy surfaces (PESs) obtained from density functional theory (DFT). The resulting FF, trained on a generic trimetallic cluster, performs well in both host−guest binding energetics and gas adsorption isotherms across different OMS-containing MOFs, including MOF-74 series and Cu-BTC. These results highlight the excellent transferability of our approach and its potential to enhance the accuracy and robustness of high-throughput MOF discovery workflows, particularly for gas adsorption and separation in large and diverse MOF databases.

36 MATERIALS SCIENCE↗

DOE OSTI2024

ReaxFF Reactive Force Field for Exploring Electronically Switchable Polarization in Zn 1– x Mg x O Ferroelectric Semiconductors

Cation misfit in traditional ferroelectric crystals offers a new material platform that can drive electronic components toward structural miniaturization and high-density integration, enabling deviation from von-Neumann architectures. Here, we explore ferroelectricity in Zn 1–x Mg x O, a nontraditional ferroelectric material with tunable properties. Using data from density-functional theory calculations, we have developed a ReaxFF reactive force field to explore the ferroelectric properties of Zn 1–x Mg x O and reveal the hysteresis behavior. We discover that ferroelectric switching can be observed at a critical thickness of 10 nm with a residual polarization of ~100 μC/cm 2 . Our analysis indicates that an increase in Mg-substitution correlates with a decrease in the coercive field. We also observe a strong temperature dependence of the coercive field in Zn 1–x Mg x O, with values decreasing as the temperature increases. Additionally, we find that the distribution of Mg atoms significantly impacts the coercive field, with a clustered distribution leading to a substantial increase. In particular, a decrease in coercive field values is observed when Mg atoms are randomly distributed, compared to uniform distribution. Furthermore, leveraging tunable hysteresis behavior offered by varied percentages and distribution of Mg-substitution provides valuable insights into the design of next-generation functional devices and will inspire further investigations.

36 MATERIALS SCIENCE↗

DOE OSTI2025

Kelvin Probe Force Microscopy Imaging of Plasticity in Hydrogenated Perovskite Nickelate Multilevel Neuromorphic Devices

Ion drift in nanoscale electronically inhomogeneous semiconductors is among the most important mechanisms being studied for designing neuromorphic computing hardware. However, nondestructive imaging of the ion drift in operando devices directly responsible for multiresistance states and synaptic memory represents a formidable challenge. Here, we present Kelvin probe force microscopy imaging of hydrogen-doped perovskite nickelate device channels subject to high-speed electric field pulses to directly visualize proton distribution by monitoring surface potential changes spatially, which is also supported with finite element-based electric field distribution studies. First-principles calculations provide mechanistic insights into the origin of surface potential changes as a function of hydrogen donor doping that serves as the contrast mechanism. We demonstrate 128 (7-bit) nonvolatile conductance levels in such devices relevant to in-memory computing applications. The synaptic plasticity measurements are implemented in spiking neural networks and show promising results for classification (SciKit Learn’s Iris and Wine data sets) and control (OpenAI’s CartPole-v1 and BipedalWalker-v3) simulation tasks.

Kelvin probe force microscopy↗

DOE OSTI2024

Impacts of Topography-Based Subgrid Scheme and Downscaling of Atmospheric Forcing on Modeling Land Surface Processes in the Conterminous US

The effects of small-scale topography-induced land surface heterogeneity are not well represented in current Earth System Models (ESMs). In this study, a new topography-based subgrid structure referred to as topographic units (TGU) designed to better capture subgrid topographic effects, and methods to downscale atmospheric forcing to the land TGUs have been implemented in the Energy Exascale Earth System Model (E3SM) Land Model (ELM). Effects of the subgrid scheme and downscaling methods on ELM simulated land surface processes are evaluated over the conterminous United States (CONUS). For this purpose, ELM simulations are performed using two configurations without (NoD ELM) and with (D ELM) downscaling, both using TGUs derived for the 0.5-degree grids and the same land surface parameters. Simulations using the two ELM configurations are compared over the CONUS domain, regional levels, and at observational sites (e.g., SNOTEL). The CONUS-level results suggest that D ELM simulates more snowfall and snow water equivalent (SWE), higher runoff, and less ET during spring and summer. Regional-level results suggest more pronounced impacts of downscaling over regions dominated by higher elevation TGUs and regions with maximum precipitation occurring during cool seasons. Results at the SNOTEL sites suggest that D ELM has superior capability of reproducing the observed SWE at 83% of the sites, with more pronounced performance over topographically heterogeneous TGUs with their maximum precipitation occurring during cool seasons. The results highlight the importance of improving representation of small-scale surface heterogeneity in ESMs and motivate future research to understand their effects on land-atmosphere interactions, streamflow, and water resources management over mountainous regions.

54 ENVIRONMENTAL SCIENCES↗