Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “distributed data structures”

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 1,729 records · Page 96

SU-E-T-318: Dosimetric Evaluation of ArcCHECK and 3DVH System Using Customized Polymer Gel Dosimeter

Purpose: ArcCHECK and 3DVH system (Sun Nuclear) can reconstruct the three-dimensional (3D) dose distribution and provide the DVH analysis in a patient. The aim of this study was to evaluate dosimetric accuracy of this system using customized polymer gel dosimeter, and also Gafchromic EBT3 films. Methods: Polyacrylamide-based gel contained magnesium chloride as a sensitizer (iPAGAT) was used in this study. Volumetric-modulated arc therapy (VMAT) plan was performed for the C-shape structure by the Eclipse treatment planning system (Varian) and used to irradiate the ArcCHECK by the Novalis Tx linear accelerator (Varian/BrainLAB). The cubic phantom filled with iPAGAT and EBT3 films placed in three orthogonal planes (axial, sagittal, and coronal) inserted into the I’mRT Phantom (IBA Dosimetry) simulated a patient were irradiated with the same VMAT plan. The measurement-guided 3D dose distribution was reconstructed using 3DVH software from the measured data of the ArcCHECK. The 3D dose distribution in iPAGAT was read out by Signa 1.5 T MRI system (GE), and 2D dose distribution on EBT3 was read out by color scanner (Epson). The comparison of all the dose distributions was performed with dose profiles and gamma index analysis in orthogonal planes using in-house developed software. Results: A good agreement was observed by overlaying the dose profiles of 3DVH, EBT3, and iPAGAT. The mean pass rates by gamma index analysis with 3%/3 mm criteria in orthogonal planes were 94.3% (3DVH vs EBT3), 91.1% (3DVH vs iPAGAT), and 96.4% (iPAGAT vs EBT3), respectively. Conclusion: 3D dose distribution reconstructed by ArcCHECK and 3DVH system was estimated accurately in a patient. However, slightly differences were observed between 3DVH and iPAGAT because of MRI noise, therefore further study is required to improve the accuracy of MRI based polymer gel dosimetry for the DVH analysis.

60 APPLIED LIFE SCIENCES↗

Study of the Effective Refractive Index Profile in Self-Assembling Nanostructured ITO Films

The optical characteristics and structural features of self-assembling nanostructured indium–tin oxide (ITO) films deposited onto substrates heated above 400°C are studied. Estimations of the distribution profiles of the material and effective refractive index in the films under study are obtained by computer simulation of a medium having spectral dependences of light transmission and reflection coefficients closest to those observed in the experiment. These results are in good agreement with scanning electron microscopy data for these films and make it possible to not only confirm the gradient character of such coatings but also restore some of its features. The further overgrowth of voids in self-assembling nanostructured films by means of the deposition of an additional material amount by magnetron sputtering at room temperature leads to a variation in the effective refractive index profile in them. Thus, it is possible to form transparent conductive films with an efficient refractive index profile adjusted for specific tasks. Adjustment of the refractive index profile in self-assembled ITO films acquires a special value for designing coatings with specified properties by virtue of a limited set of transparent conductive materials.

36 MATERIALS SCIENCE↗

Charge Carrier Transport and Deep Levels Recharge in Avalanche S-Diodes Based on GaAs

Carrier transport and deep-level recharging in semiconductor avalanche S-diode structures have been investigated. Gallium-arsenide n{sup +}–π–ν–n structures with the diffusion distribution of deep iron acceptors have been studied. It has been found by solving the continuity and Poisson equations with the use of a commercial software that the electron injection affects the avalanche breakdown voltage and the spacecharge region broadens due to capture of avalanche holes on negative iron ions in the π-region. It is demonstrated by comparing the results of numerical calculation with the experimental data that the S-shaped I–V characteristic of the diffusion avalanche S-diodes cannot be explained within the previously proposed mechanism of capture of avalanche holes on the deep iron levels.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Mapping All the Mass in the Universe (with Weak Gravitational Lensing) - Oral Presentation [Slides]

Recent discoveries have shown that most of the universe is made of an exotic dark matter that behaves much differently than the normal matter that we experience in everyday life. As we cannot detect dark matter directly, we must infer its location in the universe by indirect effects, such as the distortion of light from distant galaxies as it travels through large clouds of dark matter. While the degree of distortion should be proportional to the amount of dark matter present, we don't know the original shape of the galaxy so the distortion is difficult to quantify. If we had a model of how dark matter is linked to galaxies, and thus be able predict the amount of distortion that should occur, we could apply the model to galaxy surveys to map out the dark matter in our universe. In this research, I attach a spherical 'halo' of dark matter to each galaxy in a simulated universe to approximate its known complex dark matter structure. I then predict how the halos distort the light from distant galaxies generated behind the halos. As the data is simulated, the true distortion of the light is known which is compared to the halo-predicted distortion. I find that, on average, the model under-predicts the degree of distortion on all scales and fails to capture distortions from large-scale dark matter structure. These issues are likely due to missing features in the model, as the halo model is a greatly simplified version of the actual distribution of dark matter. Potential improvements to the model for future work are discussed.

79 ASTRONOMY AND ASTROPHYSICS↗

Final Scientific/Technical Report

A variety of calcifying organisms produce a transient or metastable amorphous calcium carbonate (ACC) precursor phase that is assembled and subsequently transformed into a crystalline biomineral, typically calcite or aragonite. The complex shapes, hierarchical structures, and unique physical properties of the biominerals that result from this calcification pathway have stimulated interest in adapting these concepts for the design and creation of bio-inspired functional materials in the laboratory. ACC also forms as a reactive precursor in diverse inorganic systems and is likely to play a much broader role in calcium carbonate formation. Knowledge of the structure, composition, and behavior of this metastable phase is critical for establishing a structural and mechanistic framework for calcium carbonate formation and its role in biogeochemical processes, including carbon cycling. Minor additives, such as magnesium, phosphorus, and organic macromolecules, are known to play important roles in controlling ACC stability, transformation kinetics, and selection of final crystalline polymorph. Molecular water also occurs in many types of ACC and is thought to play a structural role in its stability and transformation behavior. One of the major challenges that remain unresolved is identification of the structural basis for the role of these minor additives and molecular water. The absence of long-range order in ACC, and other amorphous phases, has posed a challenge for study by techniques commonly used for crystalline solids. Preliminary studies in our group show that the combination of two techniques, synchrotron X-ray-based pair distribution function (PDF) analysis and nuclear magnetic resonance (NMR) spectroscopy can provide entirely new insight to structural properties of synthetic ACC over length scales that are most relevant for understanding its transformation properties. Building on preliminary experiments, we propose a systematic study of synthesis, structure, and transformation behavior in abiotic systems. The work will specifically address the influence of phosphate as a minor additive. PDF analysis will utilize total X-ray scattering data collected at synchrotron facilities optimized for this method and will provide direct characterization of the short- and intermediate-range structure of ACC synthesized under controlled conditions. Parallel computational work using reverse Monte Carlo methods will allow structural models to be constructed for a more complete analysis of PDF results. NMR spectroscopy, using a variety of single- and double-resonance techniques, will provide information on H and CO3 components, including dynamical properties, and their relationship to stabilizing additives. PDF and NMR results will be complemented by parallel studies using X-ray absorption and FT-IR spectroscopy to allow direct comparison to previous studies. These techniques will be used to follow the transformation of different ACC samples, with specific additives, to crystalline phases under controlled condition (e.g., relative humidity). This work will provide the structural and mechanistic basis for understanding ACC stability, its transformation behavior, and the factors that govern polymorph selection. This new insight will directly benefit researchers in diverse fields, as well as adding to the framework of knowledge for understanding and controlling calcium carbonate formation in natural and engineered systems.

58 GEOSCIENCES↗

Data-Driven Optimization for the Design and Control of Large-Scale Systems (LDRD Final Report)

Engineering decisions are often formulated as optimization problems such as the optimal design or control of physical systems. In these applications, the resulting optimization problems are constrained by large-scale simulations involving systems of partial differential equations (PDEs), ordinary differential equations (ODEs), and differential algebraic equations (DAEs). In addition, critical components of these systems are fraught with uncertainty, including unverifiable modeling assumptions, unknown boundary and initial conditions, and uncertain coefficients. Typically, these components are estimated using noisy and incomplete data from a variety of sources (e.g., physical experiments). The lack of knowledge of the true underlying probabilistic characterization of model inputs motivates the need for optimal solutions that are robust to this uncertainty. In this report, we introduce a framework for handling "distributional" uncertainties in the context of simulation-based optimization. This includes a novel measure discretization technique that will lead to an adaptive optimization algorithm tailored to exploit the structures inherent to simulation- based optimization.

42 ENGINEERING↗

ARM-IRL: Adaptive Resilience Metric Quantification Using Inverse Reinforcement Learning

The resilience of safety-critical systems is gaining importance due to the rise in cyber and physical threats, especially within critical infrastructure. Traditional static resilience metrics may not capture dynamic system states, leading to inaccurate assessments and ineffective responses to cyber threats. This work aims to develop a data-driven, adaptive method for resilience metric learning. We propose a data-driven approach using inverse reinforcement learning (IRL) to learn a single, adaptive resilience metric. The method infers a reward function from expert control actions. Unlike previous approaches using static weights or fuzzy logic, this work applies adversarial inverse reinforcement learning (AIRL), training a generator and discriminator in parallel to learn the reward structure and derive an optimal policy. The proposed approach is evaluated on multiple scenarios: optimal communication network rerouting, power distribution network reconfiguration, and cyber–physical restoration of critical loads using the IEEE 123-bus system. The adaptive, learned resilience metric enables faster critical load restoration in comparison to conventional RL approaches.

97 MATHEMATICS AND COMPUTING↗

Microstructural, Optical, and Magnetic Properties of Vanadium-Substituted Nickel Spinel Nanoferrites

The current study investigates the impact of vanadium substitution on the structural, magnetic, and optical properties of NiFe{sub 2−x}V{sub x}O{sub 4} (x ≤ 0.3) nanoparticles (NPs) produced by the cost-effective sol-gel route. The as-prepared spinel ceramic powders were examined by X-ray diffraction (XRD), UV-visible diffuse reflectance spectroscopy (DRS), Fourier transform infrared (FT-IR) spectroscopy, scanning electron microscopy (SEM), and vibrating sample magnetometry (VSM). The functional groups, spinel phase, and crystal structure were confirmed by XRD and FT-IR, respectively. The crystallites size decreased from 45.24 to 36.56 nm as the doping process increases. The plots of Tauc were drawn to determine optical band gap magnitudes of 1.291, 1.302, and 1.312 eV for x = 0.0, 0.2, and 0.3, respectively. The estimated saturation magnetization is maximum for pristine NiFe{sub 2}O{sub 4} NPs and decreases to minimum for NiFe{sub 1.7}V{sub 0.3}O{sub 4} NPs. The σ-H hysteresis loops have finite coercivity (between 125 and 169 Oe) and retentivity (between 9.36 and 14.04 emu/g) values. The calculated σ{sub r}/σ{sub s} ratios are lower than 0.500, assigning the uniaxial anisotropy for NiFe{sub 2−x}V{sub x}O{sub 4}. The effective anisotropy constants (K{sub eff}) are in the range of 0.824 × 10{sup 5} and 1.303 × 10{sup 5} Erg/g. The magnetocrystalline anisotropy field (H{sub a}) values are around 5.0 kOe. The characteristics of hysteresis (σ-H) curves and the order of magnetic data reveal the soft ferrimagnetic feature of as-prepared nanoparticle samples. From Mossbauer analysis, the variations in hyperfine magnetic field, quadrupole splitting, line width, and isomer shift have been evaluated. The distribution of cations showed that the octahedral B sites are occupied by all the ions of V{sup 3+}. Mossbauer spectra are composed of four Zeeman sextets and one doublet.

77 NANOSCIENCE AND NANOTECHNOLOGY↗

Joint Resource Modeling and Assessment for Hybrid Distributed Solar and Wind Systems

The inherent variability and uncertainty in distributed energy resources can presents myriad challenges to the planning and operations of power systems. These risks are poised to become larger as the penetration of renewable energy sources rises in the power generation mix. Hybrid solar-wind energy systems are able to mitigate some of these risks by their complementary resource availability. Surface solar and wind fields are coupled and correlated in both space and time. Appropriately estimating the hybrid solar wind energy system requires simulating the spatio-temporal structure of these fields that can be produced for each time horizon. We introduce a novel joint spatio-temporal stochastic differential equation (SPDE) approach that captures the spatio-temporal dynamics of solar and wind fields and their joint dependency over a domain for each time step. In the case study on Colorado, we consider nonstationary three-level hierarchical spatio temporal models for both hourly solar irradiance data and wind speed data in Colorado. Dependence between the solar irradiance data and wind speed data is captured by a shared spatio-temporal random effect. Our approach performs well in terms of the prediction score criterion.

joint modeling↗

Seismicity and structure of Akutan and Makushin Volcanoes, Alaska, using joint body and surface wave tomography

Joint inversions of seismic data recover models that simultaneously fit multiple constraints while playing upon the strengths of each data type. Here, we jointly invert 14 years of local earthquake body wave arrival times from the Alaska Volcano Observatory catalog and Rayleigh wave dispersion curves based upon ambient noise measurements for local V p , V s , and hypocentral locations at Akutan and Makushin Volcanoes using a new joint inversion algorithm.The velocity structure and relocated seismicity of both volcanoes are significantly more complex than many other volcanoes studied using similar techniques. Seismicity is distributed among several areas beneath or beyond the flanks of both volcanoes, illuminating a variety of volcanic and tectonic features. The velocity structures of the two volcanoes are exemplified by the presence of narrow high-V p features in the near surface, indicating likely current or remnant pathways of magma to the surface. A single broad low-V p region beneath each volcano is slightly offset from each summit and centered at approximately 7 km depth, indicating a potential magma chamber, where magma is stored over longer time periods. Differing recovery capabilities of the Vp and Vs datasets indicate that the results of these types of joint inversions must be interpreted carefully.

58 GEOSCIENCES↗

Analysis of thermal and mechanical properties with inventory level of the molten salt storage tank in central receiver concentrating solar power plants

Molten salt thermal energy storage (TES) tanks ensure steady power output of concentrating solar power (CSP) plants; however, recent tank failures have highlighted the need for further analysis. Current studies primarily focus on analyzing the molten salt flow, heat transfer, and thermal efficiency. Additionally, research on the latest tank structures is limited and lacks newest experimental validation. This study measures temperature and molten salt inventory levels in the high-temperature tank at a 50 MW central receiver CSP plant, connected to the power grid in 2019. A multi-physics model was developed to evaluate thermal and mechanical properties of TES tanks by combining computational fluid dynamics and finite element modeling using real plant data. Heat loss, temperature, displacement, and stress distribution of the tank at different inventory levels were investigated. Results show that ambient air velocity near the tank roof reaches 2.14 m/s, much higher than 0.2 m/s near the wall. The temperatures of inventory fluid and tank are close, varying slightly at different levels due to thermal conduction and radiation. Because the heat loss strongly depends on temperature, the total tank loss remains nearly constant across inventory levels. Larger temperature gradients and thermal stresses are primarily localized along the tank floor edge and the air-salt interface. Notably, the maximum thermal stress at the tank edge is three times higher than that at the interface. The magnitude of total stress changes by less than 5 MPa with and without thermal load, indicating that high temperatures exert only a minor impact on tank stress. In contrast, thermal load significantly affects tank deformation, particularly at the roof edge, where values exceed 150 mm. Despite the large variation in molten salt levels, tank wall temperatures and displacements present a minor change, suggesting a weak correlation with inventory levels. In conclusion, the findings obtained in this study provide important insights on the TES tank that could be used to optimize tank design and operation strategies.

14 SOLAR ENERGY↗

Residual service life of carbonated structures based on site non-destructive tests

Assessing the residual service life of carbonated reinforced concrete structures is typically conducted by measuring cover depths and carbonation of drilled cores, assuming the ''square root law'' for the progress of the carbonation front. Given the large spatial variability of carbonation depths, the required number of cores becomes prohibitive. A model is presented, based on numerous experimental data, by which the carbonation rate can be estimated by means of non-destructive site measurements of the air-permeability kT. A Weibull distribution of the carbonation rate K{sub c} is fitted to the results, which is a function of kT. Hence, for each kT value measured on site coupled with the cover depth x{sub d}, expected, optimistic and pessimistic values of the corrosion initiation time can be obtained without the need of drilling cores. The method is validated on 111 paired results of kT and x{sub d}, obtained on the emblematic Tokyo's Museum of Western Art.

36 MATERIALS SCIENCE↗

Dynamically Disordered Lattice in a Layered Pb-I-SCN Perovskite Thin Film Probed by Two-Dimensional Infrared Spectroscopy

The dynamically flexible lattices in lead-halide perovskites may play important roles in extending carrier recombination lifetime in 3D perovskite solar-cell absorbers and in exciton self-trapping in 2D perovskite white-light phosphors. Here, two-dimensional infrared (2D IR) spectroscopy was applied to study a recently reported Pb-I-SCN layered perovskite. The Pb-ISCN perovskite was spin-coated on a SiO 2 surface as a thin film, with a thickness of ~100 nm, where the S 12 CN – anions were isotopically diluted with the ratio of S 12 CN:S 13 CN=5:95 or 1:99 to avoid vibrational coupling and excitation transfer between adjacent SCN – anions. The minor S 12 CN – component was the principal vibrational probe that reported on the structural evolution through 2D IR spectroscopy. Spectral diffusion was observed with a time constant of 4.1 ± 0.3 ps. Spectral diffusion arises from small structural changes that result in sampling of frequencies within the distribution of frequencies that comprise the inhomogeneously broadened absorption spectrum of the in the Pb-I-SCN layered perovskite. These transitions among discrete local structures are distinct from oscillatory phonon motions of the lattice. To accurately evaluate the structural dynamics through measurement of spectral diffusion, the vibrational coupling between adjacent SCN – anions had to be carefully treated. Although the inorganic layers of typical 2D perovskites are structurally isolated from each other, the 2D IR data demonstrated that the layers of the Pb-I-SCN perovskite are vibrationally coupled. When both S 12 CN – and S 13 CN – were pumped simultaneously, cross-peaks between S 12 CN and S 13 CN vibrations and an oscillating 2D band shape of the S 12 CN – vibration were observed. Both observables demonstrate vibrational coupling between the closest SCN – anions, which reside in different inorganic layers. The thin films and the isotopic dilution produced exceedingly small vibrational echo signal fields; measurements were made possible using the near-Brewster’s angle reflection pump-probe geometry.

14 SOLAR ENERGY↗

X-ray Mapping of Nanoparticle Superlattice Thin Films

We combine grazing-incidence and transmission small-angle X-ray diffraction with electron microscopy studies to characterize the structure of nanoparticle films with long-range order. Transmission diffraction is used to collect in-plane diffraction data from single grains and locally aligned nanoparticle superlattice films. Systematic mapping of samples can be achieved by translating the sample in front of the X-ray beam with a spot size selected to be on the order of superlattice grain features. This allows a statistical determination of superlattice grain size and size distribution over much larger areas than typically accessible with electron microscopy. Transmission X-ray measurements enables spatial mapping of the grain size, orientation, uniformity, strain, or crystal projections and polymorphs. Furthermore, we expand this methodology to binary nanoparticle superlattice and nanorod superlattice films. Our study provides a framework for characterization of nanoparticle superlattices over large areas which complements or expands microstructure information from real-space imaging.

GISAXS↗

Neutron diffraction, Mössbauer effect and electron paramagnetic resonance studies on multiferroic Pb(Fe{sub 2/3}W{sub 1/3})O{sub 3}

Multiferroic Pb(Fe{sub 2/3}W{sub 1/3})O{sub 3} ceramics were synthesized via a modified two-stage Columbite method. Single phase formation was confirmed from the analysis of x-ray and neutron diffraction patterns recorded at room temperature. Structural analysis of the diffraction data reveals cubic phase (space group Pm-3m) for the title compound. Magnetic structure of the title compound at room temperature exhibits G-type antiferromagnetic structure. The Mössbauer spectroscopy and Electron Paramagnetic Resonance (EPR) studies were carried out at 300 K. The isomer shift and quadrupole splitting of the Mössbauer spectra confirms the trivalent state of iron (Fe{sup 3+}). The Mössbauer spectra also suggest that the iron and tungsten are randomly distributed at the octahedral, B site. EPR spectra show a single broad line associated with Fe{sup 3+} ions. Both spectra clearly exhibit weak ferromagnetic behaviour of Pb(Fe{sub 2/3}W{sub 1/3})O{sub 3} ceramic at 300 K. Considering neutron diffraction, Mössbauer and EPR results together, it may be stated here that Pb(Fe{sub 2/3}W{sub 1/3})O{sub 3} exhibits antiferromagnetic behavior along with weak ferromagnetism at room temperature.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

SMALL-SCALE PRESSURE-BALANCED STRUCTURES DRIVEN BY MIRROR-MODE WAVES IN THE SOLAR WIND

Recently, small-scale pressure-balanced structures (PBSs) have been studied with regard to their dependence on the direction of the local mean magnetic field B{sub 0} . The present work continues these studies by investigating the compressive wave mode forming small PBSs, here for B{sub 0} quasi-perpendicular to the x-axis of Geocentric Solar Ecliptic coordinates (GSE-x). All the data used were measured by WIND in the quiet solar wind. From the distribution of PBSs on the plane determined by the temporal scale and angle θ{sub xB} between the GSE-x and B{sub 0} , we notice that at θ{sub xB} = 115° the PBSs appear at temporal scales ranging from 700 s to 60 s. In the corresponding temporal segment, the correlations between the plasma thermal pressure P{sub th} and the magnetic pressure P{sub B}, as well as that between the proton density N{sub p} and the magnetic field strength B, are investigated. In addition, we use the proton velocity distribution functions to calculate the proton temperatures T and T{sub ∥}. Minimum Variance Analysis is applied to find the magnetic field minimum variance vector B{sub N} . We also study the time variation of the cross-helicity σ{sub c} and the compressibility C{sub p} and compare these with values from numerical predictions for the mirror mode. In this way, we finally identify a short segment that has T > T{sub ∥}, proton β ≅ 1, both pairs of P{sub th}-P{sub B} and N{sub p}-B showing anti-correlation, and σ{sub c} ≈ 0 with C{sub p} > 0. Although the examination of σ{sub c} and C{sub p} is not conclusive, it provides helpful additional information for the wave mode identification. Additionally, B{sub N} is found to be highly oblique to B{sub 0} . Thus, this work suggests that a candidate mechanism for forming small-scale PBSs in the quiet solar wind is due to mirror-mode waves.

79 ASTRONOMY AND ASTROPHYSICS↗

New scaling and nuclear structure aspects in heavy-ion fusion reactions

Three new behaviors have been found in comparisons of fusion cross sections for different collision systems. root (1) Replacing the energy E with a scaling one, E scal = (E-V g )/($\sqrt{2}$W g ), is successful for washing out the Coulomb interaction in the spectra of fusion cross sections, where V g and W g are barrier height and width of the single-Gaussian barrier distribution model. (2) In a representation of σE vs the scaling energy, E scal , all data sets display in parallel. Here, the ratio for sigma E from any two fusion systems over the whole range is a constant value. That behavior is also studied in another representation, in which the data sets display as parallel horizontal lines for any heavy-ion fusion system. (3) The constant ratio value is the ratio of parameter products, $R^2_gW_g$, of the two systems; where R g is the barrier radius obtained in the single-Gaussian barrier distribution model. Moreover, when comparing neighboring collision systems at the same E scal , the ratio of sigma is near a constant value within a few percent over the whole range. Thus a quantitative comparison for the fusion enhancement for neighboring systems is developed. The present finding could be beneficial for predicting unmeasured fusion cross sections.

Jiang, C. L. [Argonne National Laboratory (ANL), A↗

Hydrologic Model Files (PRMS) for Historical Conditions in the East River Watershed, Colorado between 1987-2019

The data package contains model input and output file for the East River watershed (85 km2) located in the headwaters of the Colorado River. The model is the U.S. Geological Survey (USGS) Precipitation-Modeling Runoff System (PRMS) for the historical, or baseline, condition (1987-2019) to assess the efficiency of summer monsoon rains to generate streamflow. The finite difference grid resolution is 100 m with elevations resampled from the USGS National Elevation Dataset. LANDFIRE raster data is used to derive parameters of dominant cover type, summer and winter cover density, canopy interception characteristics for snow and rain and transmission coefficient for shortwave radiation. Climate forcing uses minimum and maximum daily temperature lapse rates defined by the two Snow Telemetry (SNOTEL) stations located in the vicinity of the East River adjusted for aspect. Schofield SNOTEL snowfall is spatially distributed using LiDAR derived snow depth observations from the Airborne Snow Observatory (ASO) flown 4 April 2016. Snow depths are converted to snow water equivalent based on ground surveys and density modeling. Rainfall is spatially distributed using the monthly Parameter-elevation Relationships on Independent Slopes Model 30-year monthly averages. Simulated solar radiation is calibrated to match weather station observations. Model verification of snow water equivalent (SWE) accumulation and ablation relies on the 2018 and 2019 ASO maps and weather station snow depth. The data package contains model input files in text format, model documentation in pfd and a readme.txt file to help users understand the file structure to run the program.

54 ENVIRONMENTAL SCIENCES↗