Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “state dependent variable”

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 829 records · Page 46

Water Balance in the Amazon Basin from a Land Surface Model Ensemble

Despite recent advances in land surfacemodeling and remote sensing, estimates of the global water budget are still fairly uncertain. This study aims to evaluate the water budget of the Amazon basin based on several state-ofthe- art land surface model (LSM) outputs. Water budget variables (terrestrial water storage TWS, evapotranspiration ET, surface runoff R, and base flow B) are evaluated at the basin scale using both remote sensing and in situ data. Meteorological forcings at a 3-hourly time step and 18 spatial resolution were used to run 14 LSMs. Precipitation datasets that have been rescaled to matchmonthly Global Precipitation Climatology Project (GPCP) andGlobal Precipitation Climatology Centre (GPCC) datasets and the daily Hydrologie du Bassin de l'Amazone (HYBAM) dataset were used to perform three experiments. The Hydrological Modeling and Analysis Platform (HyMAP) river routing scheme was forced with R and B and simulated discharges are compared against observations at 165 gauges. Simulated ET and TWS are compared against FLUXNET and MOD16A2 evapotranspiration datasets andGravity Recovery and ClimateExperiment (GRACE)TWSestimates in two subcatchments of main tributaries (Madeira and Negro Rivers).At the basin scale, simulated ET ranges from 2.39 to 3.26 mm day(exp -1) and a low spatial correlation between ET and precipitation indicates that evapotranspiration does not depend on water availability over most of the basin. Results also show that other simulated water budget components vary significantly as a function of both the LSM and precipitation dataset, but simulated TWS generally agrees with GRACE estimates at the basin scale. The best water budget simulations resulted from experiments using HYBAM, mostly explained by a denser rainfall gauge network and the rescaling at a finer temporal scale.

GPCP↗

Singular perturbation analysis of AOTV-related trajectory optimization problems

The problem of real time guidance and optimal control of Aeroassisted Orbit Transfer Vehicles (AOTV's) was addressed using singular perturbation theory as an underlying method of analysis. Trajectories were optimized with the objective of minimum energy expenditure in the atmospheric phase of the maneuver. Two major problem areas were addressed: optimal reentry, and synergetic plane change with aeroglide. For the reentry problem, several reduced order models were analyzed with the objective of optimal changes in heading with minimum energy loss. It was demonstrated that a further model order reduction to a single state model is possible through the application of singular perturbation theory. The optimal solution for the reduced problem defines an optimal altitude profile dependent on the current energy level of the vehicle. A separate boundary layer analysis is used to account for altitude and flight path angle dynamics, and to obtain lift and bank angle control solutions. By considering alternative approximations to solve the boundary layer problem, three guidance laws were derived, each having an analytic feedback form. The guidance laws were evaluated using a Maneuvering Reentry Research Vehicle model and all three laws were found to be near optimal. For the problem of synergetic plane change with aeroglide, a difficult terminal boundary layer control problem arises which to date is found to be analytically intractable. Thus a predictive/corrective solution was developed to satisfy the terminal constraints on altitude and flight path angle. A composite guidance solution was obtained by combining the optimal reentry solution with the predictive/corrective guidance method. Numerical comparisons with the corresponding optimal trajectory solutions show that the resulting performance is very close to optimal. An attempt was made to obtain numerically optimized trajectories for the case where heating rate is constrained. A first order state variable inequality constraint was imposed on the full order AOTV point mass equations of motion, using a simple aerodynamic heating rate model.

Calise, Anthony J.↗

A Framework for Deep Learning Emulation of Numerical Models With a Case Study in Satellite Remote Sensing

Numerical models based on physics represent the state of the art in Earth system modeling and comprise our best tools for generating insights and predictions. Despite rapid growth in computational power, the perceived need for higher model resolutions overwhelms the latest generation computers, reducing the ability of modelers to generate simulations for understanding parameter sensitivities and characterizing variability and uncertainty. Thus, surrogate models are often developed to capture the essential attributes of the full-blown numerical models. Recent successes of machine learning methods, especially deep learning (DL), across many disciplines offer the possibility that complex nonlinear connectionist representations may be able to capture the underlying complex structures and nonlinear processes in Earth systems. A difficult test for DL-based emulation, which refers to function approximation of numerical models, is to understand whether they can be comparable to traditional forms of surrogate models in terms of computational efficiency while simultaneously reproducing model results in a credible manner. A DL emulation that passes this test may be expected to perform even better than simple models with respect to capturing complex processes and spatiotemporal dependencies. Here, we examine, with a case study in satellite-based remote sensing, the hypothesis that DL approaches can credibly represent the simulations from a surrogate model with comparable computational efficiency. Our results are encouraging in that the DL emulation reproduces the results with acceptable accuracy and often even faster performance. We discuss the broader implications of our results in light of the pace of improvements in high-performance implementations of DL and the growing desire for higher resolution simulations in the Earth sciences.

Bayesian Deep Learning↗

SPRUCE Soil Volumetric Water Content in Experimental Plots, Marcell Experimental Forest, Minnesota, 2017-2025

This dataset contains soil volumetric water content (VWC) measurements from the Spruce and Peatland Responses Under Changing Environments (SPRUCE) Experiment. Measurements were made inside SPRUCE experimental plots in the S1 Bog at the US Forest Service Marcell Experimental Forest in northern Minnesota, USA from 2018-2025 (2018-06-14 to 2025-12-31). Observations were made with METER 10HS soil moisture sensors. To reduce bulk density related variability, 10HS sensors were placed inside mesh tubes filled with peat at a standard bulk density. The observations under standard bulk density represent relative differences in water content between hummock and hollow positions, and in response to experimental treatments. Productivity of peatlands, and their keystone species sphagnum, are highly dependent on water availability, which is affected by lateral inputs, precipitation, ground water depth and evapotranspiration. Knowledge of near surface and sphagnum water content is useful to understand peatland function. This dataset contains 8 data files in comma-separate values (.csv) format. Additional metadata are provided: one data dictionary and a file-level metadata file in comma-separate values (*.csv) format and a user guide in PDF (*.pdf) format.

Warren, Jeffrey [ORNL] (ORCID:0000000206804697)↗

Simulating Isotropic Rayleigh Wave Horizontal to Vertical Ratios at Rock Valley

Elastic anisotropy in the subsurface can manifest from pre-existing structure or stress-induced cracking – both phenomena are of interest at the Rock Valley site in Nevada. The azimuthal dependence of Rayleigh wave ellipticity, also known as H/V, can reveal elastic anisotropy, but the effect of isotropic variability in H/V must first be accounted for. To that end, we simulate synthetic seismograms from a ring of sources in a fully 3-D isotropic velocity model that captures the geology at Rock Valley. We measure H/V from these waveforms and compare them with H/V predictions obtained using 1D models extracted at specific locations from the 3D model. We find that an average of the azimuthal H/V measurements from the 3-D simulations and 1-D predictions agree, but that the azimuthal variability of H/V is a large fraction of the mean value. These results estimate the noise character in azimuthal H/V so it can be factored out of future anisotropic simulations.

58 GEOSCIENCES↗

Pipe Parity Analysis of Seawater Desalination in the United States: Exploring Costs, Energy, and Reliability via Case Studies and Scenarios of Emerging Technology

As climate change-induced variables exacerbate water scarcities, the use of seawater reverse osmosis (SWRO) membrane desalination technology for treating seawater could potentially provide a long-term drought-proof source of drinking water. This study carried out a technoeconomic assessment (TEA) for three SWRO desalination plants in the U.S. and one in Israel to set baselines for the cost and energy consumption for seawater desalination. In addition, a breakeven curve for implementing SWRO desalination was estimated in relation to the cost of water conservation measures to meet drought-induced reduction of the traditional water supply. The results show that the cost of SWRO water production scales with the plant capacity while energy intensity is not dramatically different across the plants. The higher cost in some U.S. plants is due to high capital investment, including land acquisition and permitting. Variations in plant capacity utilization have the greatest impact on the levelized cost of water (LCOW) over the plant service life, suggesting the importance of reducing fouling and maintenance-related downtime. Scenario analysis of fixed labor cost reduction through process automation indicates investments in automation and sensing technology could result in long-term savings. Breakeven analysis shows the decision to adopt SWRO is highly dependent on the local cost associated with water conservation to meet water supply reduction. Moreover, a small reduction in SWRO cost can influence a shift toward the adoption of SWRO over water conservation measures. Incorporation of future water demand, water conservation potential, and water stress data around the nation indicates SWRO desalination could be an important contributor to the future municipal drinking water portfolio in the U.S.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗

Qubit Assignment Using Time Reversal

As quantum computers with large numbers of qubits become increasingly available, experiments executed on a given device may not utilize all available qubits. In this case, the outcome of executing a quantum program will depend on the ability to efficiently select a subset of high-performing physical qubits. For any given quantum program and device there are many ways to assign physical qubits for execution of the program, and assignments will differ in performance due to the variability in quality across qubits and entangling operations on a single device. Evaluating the performance of each assignment using fidelity estimation introduces significant experimental overhead and will be infeasible for many applications, while relying on standard device benchmarks provides incomplete information about the performance of any specific program. Furthermore, the number of possible assignments grows combinatorially in the number of qubits on the device and in the program, motivating the use of heuristic optimization techniques. We demonstrate a practical solution to the problem of qubit assignment by using simulated annealing with a cost function based on the Loschmidt echo, a diagnostic that measures the reversibility of a quantum process. We provide theoretical justification for this choice of cost function by demonstrating that the optimal qubit assignment coincides with the optimal qubit assignment based on state fidelity in the weak error limit, and we provide experimental justification using diagnostics performed on Google’s superconducting qubit devices. We then establish the performance of simulated annealing for qubit assignment using classical simulations of noisy devices as well as optimization experiments performed on a quantum processor. Our results demonstrate that the use of Loschmidt echoes and simulated annealing provides a scalable and flexible approach to optimizing qubit assignment on near-term hardware.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Stakeholder Support for Wildlife Conservation Funding Policies

State wildlife management agencies in the United States have depended on a “user-pay” funding model for conservation efforts that relies on revenue from hunting license sales and a federal excise tax on firearms, ammunition, and archery equipment. Declines in hunting participation, however, jeopardize sustainability of the current funding model. Ensuring support among stakeholders for wildlife management and conservation may require expanding sources of funding and incorporating the perspectives and values of a diversifying constituency into decision making processes. We used a web-based survey of wildlife-associated recreationists in Michigan, USA to evaluate support for a range of conservation funding policies. Respondents self-identified primarily as hunters (n = 2,558) or wildlife watchers (n = 942). We used binary logistic regression to evaluate support for four conservation funding policy options: state sales tax, lottery proceeds, extractive industry revenue, and a user-based tax on outdoor gear (i.e., “backpack tax”). Determinants of support varied by type of policy and stakeholder characteristics. We found no statistically significant differences between hunters and wildlife watchers in their support for conservation funding policies when accounting for other variables such as wildlife value orientations, engagement in stewardship behaviors, age, and gender. The industry-based policy achieved the greatest level of approval, while the backpack tax had the lowest. Respondents were mixed in their support of the sales tax and lottery proceeds options. Cluster analysis revealed three homogenous groups related to conservation funding policies: “strong support,” “mixed/opposed,” and “anti-backpack tax.” Clusters differed in their support for conservation funding policies and on psychological and demographic variables. The “strong support” and “anti-backpack tax” groups differed in their levels of stewardship engagement, knowledge of conservation funding mechanisms, and support for the backpack tax option. The “mixed-opposed” group tended to be older, less educated, and less likely to be a member of a conservation organization. Results suggest support for conservation funding differs by policy type and social and psychological characteristics of stakeholders. Based on differences in policy support revealed in this study, we suggest a multi-tiered approach to funding conservation and building on support among wildlife stakeholders to mitigate the looming funding crisis for state wildlife agencies.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗

Optimal short-range trajectories for helicopters

An optimal flight path algorithm using a simplified altitude state model and a priori climb cruise descent flight profile was developed and applied to determine minimum fuel and minimum cost trajectories for a helicopter flying a fixed range trajectory. In addition, a method was developed for obtaining a performance model in simplified form which is based on standard flight manual data and which is applicable to the computation of optimal trajectories. The entire performance optimization algorithm is simple enough that on line trajectory optimization is feasible with a relatively small computer. The helicopter model used is the Silorsky S-61N. The results show that for this vehicle the optimal flight path and optimal cruise altitude can represent a 10% fuel saving on a minimum fuel trajectory. The optimal trajectories show considerable variability because of helicopter weight, ambient winds, and the relative cost trade off between time and fuel. In general, reasonable variations from the optimal velocities and cruise altitudes do not significantly degrade the optimal cost. For fuel optimal trajectories, the optimum cruise altitude varies from the maximum (12,000 ft) to the minimum (0 ft) depending on helicopter weight.

Slater, G. L.↗

Assessment of impact damage of composite rocket motor cases

This contract reviewed the available literature on mechanisms of low velocity impact damage in filament wound rocket motor cases, MDE methods to quantify damage, critical coupon level test methods, manufacturing and material process variables and empirical and analytical modeling off impact damage. The critical design properties for rocket motor cases are biaxial hoop and axial tensile strength. Low velocity impact damage is insidious because it can create serious nonvisible damage at very low impact velocities. In thick rocket motor cases the prevalent low velocity impact damage is fiber fracture and matrix cracking adjacent to the front face. In contrast, low velocity loading of thin wall cylinders induces flexure, depending on span length and the flexure induces delamination and tensile cracking on the back face wall opposed to impact occurs due to flexural stresses imposed by impact loading. Important NDE methods for rocket motor cases are non-contacting methods that allow inspection from one side. Among these are vibrothermography, and pulse-echo methods based on acoustic-ultrasonic methods. High resolution techniques such as x-ray computed tomography appear to have merit for accurate geometrical characterization of local damage to support development of analytical models of micromechanics. The challenge of coupon level testing is to reproduce the biaxial stress state that the full scale article experiences, and to determine how to scale the composite structure to model full sized behavior. Biaxial tensile testing has been performed by uniaxially tensile loading internally pressurized cylinders. This is experimentally difficult due to gripping problems and pressure containment. Much prior work focused on uniaxial tensile testing of model filament wound cylinders. Interpretation of the results of some studies is complicated by the fact that the fabrication process did not duplicate full scale manufacturing. It is difficult to scale results from testing subscale cylinders since there are significant differences in out time of the resins relative to full scale cylinder fabrication, differences in hoop fiber tensioning and unsatisfactory coupon configurations. It appears that development of a new test method for subscale cylinders is merited. Damage tolerance may be improved by material optimization that uses fiber treatments and matrix modifications to control the fiber matrix interface bonding. It is difficult to develop process optimization in subscale cylinders without also modeling the longer out times resins experience in full scale testing. A major breakthrough in characterizing the effect of impact damage on residual strength, and understanding how to scale results of subscale evaluations, will be a sound micromechanical model that described progressive failure of the composite. Such models will utilize a three dimensional stress analysis due to the complex nature of low velocity impact stresses in thick composites. When these models are coupled with non-contact NDE methods that geometrically characterize the damage and acoustic methods that characterize the effective local elastic properties, accurate assessment of residual strength from impact damage may be possible. Directions for further development are suggested.

Paris, Henry G.↗

Are You In or Out: the Influence of Sea Ice Drift on Sea Ice Survivability

The state of Arctic sea ice influences aspects of many global systems including ecosystems, economies, geopolitics, and climate. And, it is in severe decline. Numerous studies have considered the factors that influence Arctic sea ice decline using monthly mean or gridded data. However, these methods struggle to cleanly separate the influences on sea ice melt because the sea ice cover is always moving; thus, different sea ice floes and sea ice regimes constitute the monthly average sea ice state within the Eulerian frame. This study uses a LaGrangian sea ice parcel tracking satellite database to investigate the factors that influence sea ice parcel survivability. The concept of sea ice survivability is defined as the likelihood that a sea ice parcel will last through the summer melt season. Survivability can be stratified in many ways to analyze how it changes as a function of sea ice regime, region, and sea ice characteristics. How is the survivability of first year (FY) and multi-year (MY) ice classes influenced by sea ice motion? This is the central question addressed in this presentation. Specifically, this analysis determines sea ice survivability as a function of distance traveled by sea ice parcels and compares the survivability statistics for sea ice parcels that stay in their starting region versus those that leave their starting region. The preliminary results indicate that the influence of sea ice movement on survivability depends strongly on which region the sea ice resides in at the beginning of the growth season. Lastly, we find substantial inter-annual variability in the regional distribution of sea ice survivability that seems to be strongly linked to sea ice parcel drift between regions of the Arctic Ocean. Moreover, we use these results to discuss the potential contribution of changes in sea ice motion on the observed slowdown in Arctic sea ice extent loss.

Patrick C Taylor↗

Synthesizing land use and demographic change in Southeast Asia’s smaller urbanized areas from 2000–2015

The majority of the human population now reside in urban areas today. The United Nations estimates that nearly half of all urban dwellers currently live in cities smaller than 500 000 persons and the majority of future urban growth will take place in Asia and Africa, likely in these smaller urban areas, not mega cities. Thus, understanding the factors that influence urban demographic trajectories in small urban areas is critical to address sustainable and equitable policy initiatives related to food security, changing climate hazard exposure, and economic opportunities. Here we focus on Southeast Asia—a region historically characterized by lower urban population proportions, yet with a rapidly shifting dynamic demographic—to examine correlates of demographic change among smaller cities. We combine two open-source satellite-informed datasets: GHS urban center database (2015) and age-sex gridded data from WorldPop to calculate socio-demographic characteristics to model drivers of change in annualized urban population growth from 2000–2015 for 505 urbanized places. We find a general pattern of decreasing dependency ratios as city-size increases for most urban areas in Southeast Asia. Higher rates of growth and more variation is observed for smaller cities—those with fewer than 300 000 persons, the lowest population limit for UN data on urbanization. When examining covariates of urban population growth, we find significant statistical associations of population change in smaller urbanized areas with climatic, economic, and land cover/land use variables, but with country-specific variations. Characterizing a continuum of urban population development in the context of changing environmental, economic and climate conditions has been an important sustainable development and equity issue for decades, but newer analysis of city-level drivers allows for systematic inquiry thus moving beyond total population counts for policy-relevant insight.

Southeast Asia synthesis↗

High-resolution modeling of indoor radon exposure with uncertainty quantification in Utah

Indoor radon accounts for 37% of population-level exposure to ionizing radiation in the United States. However, radon metrics are typically reported at coarse spatial scales, potentially obscuring meaningful local variation. We developed a high-resolution modeling framework to estimate indoor radon concentrations across Utah while explicitly quantifying predictive uncertainty. A total of 19,497 residential radon measurements collected between 2006 and 2017 were combined with environmental and housing characteristics and analyzed using a geospatial neural network that accommodates spatial dependence and nonlinear associations. Predictions were generated on a uniform hexagonal grid at 0.73 km2 resolution (H3 level 8). Out-of-sample predictions aggregated to the H3 level 8 grid showed good agreement with observed concentrations (Pearson r=0.64), while household-level predictions exhibited more moderate agreement (r=0.45). The model produced well-calibrated uncertainty estimates, with 24.1% of held-out observations exceeding the predicted 75th-percentile threshold. Maps of predicted radon concentrations and the probability of exceeding the U.S. EPA action level of 148 Bq/m3 (4 pCi/L) revealed substantial fine-scale spatial heterogeneity that was not apparent in conventional coarse-resolution summaries, with greater local variability observed in densely monitored urban counties than in sparsely sampled regions. High-resolution radon models that explicitly quantify uncertainty provide a useful framework for characterizing the spatial distribution of indoor radon and identifying areas of elevated exceedance risk. These findings highlight the value of fine-scale monitoring data and uncertainty-aware modeling approaches for radon exposure assessment, environmental risk characterization, and radon-related health research.

Wu, Yunhan [ORNL] (ORCID:0000000178842994)↗

Nonlinear development and secondary instability of Gortler vortices in hypersonic flows

In a hypersonic boundary layer over a wall of variable curvature, the region most susceptible to Goertler vortices is the temperature adjustment layer over which the basic state temperature decreases monotonically to its free stream value. Except for a special wall curvature distribution, the evolution of Goertler vortices trapped in the temperature adjustment layer will in general be strongly affected by the boundary layer growth through the O(M sup 3/2) curvature of the basic state, where M is the free stream Mach number. Only when the local wavenumber becomes as large as of order M sup 3/8, do nonparallel effects become negligible in the determination of stability properties. In the latter case, Goertler vortices will be trapped in a thin layer of O(epsilon sup 1/2) thickness which is embedded in the temperature adjustment layer; here epsilon is the inverse of the local wavenumber. A weakly nonlinear theory is presented in which the initial nonlinear development of Goertler vortices in the neighborhood of the neutral position is studied and two coupled evolution equations are derived. From these, it can be determined whether the vortices are decaying or growing depending on the sign of a constant which is related to wall curvature and the basic state temperature.

Fu, Yibin B.↗

Unlocking plant-microbial interactions in deep Mollisols in the Midwestern US: Linking depth gradients in roots, microbial activity, and soil carbon in agroecosystems

Deep-rooted plants may build soil carbon (C) stocks, but most research has focused on shallow soils, leaving gaps in our understanding of how shifts in the balance between decomposition and C inputs drive soil C accumulation with depth. Thus, our objectives were to: (1) link depth gradients in root biomass with microbial activity and soil C stocks down to 1 m, and (2) examine the potential of simple C inputs to prime soil C across depths. To this end, we dug 5 quantitative soil pits in Argiudolls under mature perennial miscanthus plots in the SoyFACE Farm (Champaign-Urbana, IL). We added 13 C labeled glucose to our soils to determine the fate of simple C inputs with depth. We found that fine root biomass, total soil C, mineral-associated organic C (MAOC), particulate organic C (POC), and microbial activity (as measured by potential enzyme activity) declined with depth. POC declined more rapidly than MAOC, resulting in an increase in the ratio of MAOC-to-POC. Root biomass, enzyme activity (either acid phosphatase or n-acetyl-glucosaminadase) activity, and microbial respiration explained 74% and 38% of the variability in soil total C and MAOC, respectively, while POC was dependent on root biomass and microbial respiration (47%). Although the incorporation of simple 13 C inputs into MAOC was similar across depths, these inputs led to greater net MAOC losses in shallow soils than in deeper soils between 50 and 100 cm. The divergent impact of simple C inputs across depths may suggest that MAOC in shallow soils is more susceptible to priming losses, while C inputs into deep soils may instead be more persistent. Collectively, our results suggest that depth gradients in soil C stocks represents a balance between inputs, decomposition, and microbial necromass production and that increases in root C inputs by deep-rooted plants may have the potential to build stable MAOC.

60 APPLIED LIFE SCIENCES↗

Modelling the high-latitude ionosphere

Results of an ionospheric model program are presented which demonstrate the extreme variability of the steady state, daytime, ionospheric F region electron density and ion composition due to both neutral atmospheric changes with solar cycle, season and magnetic activity, and to the effects of ionospheric drifts caused by perpendicular electric fields. Consideration is given to the time history of the ionospheric plasma as it undergoes convective motion due to the combined effects of corotation forces and electromagnetic forces which results from the mapping of the magnetospheric cross tail electric field to the rotating ionosphere. A simple model of the convection pattern is described. The model calculates the net effect of the tendency for the plasma to corotate about the geographic pole and the E sub Bar times B sub Bar velocity induced by a perpendicular electric field mapped to a circle centered about a point 5 deg antisunward of the geomagnetic pole and oriented such that the equipotentials are parallel to the noon midnight meridian. This convection pattern shows the generally accepted features of high latitude convection, but because of the offset between the geographic and geomagnetic poles a marked universal time dependence in these features is predicted.

Raitt, W. J.↗

Disentangling the Impacts of Microtopography and Shrub Distribution on Snow Depth in a Subarctic Watershed: Toward a Predictive Understanding of Snow Spatial Variability

Snow plays a critical role in carbon cycling, vegetation dynamics, and permafrost hydrology at high latitudes by influencing surface energy exchange. Predicting snow distribution patterns is essential for understanding the evolution of Arctic ecosystems, yet scaling process-level knowledge to landscape predictions remains challenging. Here, we analyze snow depth (2019 and 2022), terrain elevation, and vegetation height from a watershed on the Seward Peninsula, Alaska, to examine how topography and shrubs shape snow redistribution across spatial scales. We find that snow depth is strongly coupled to terrain at scales below ∼60 m but becomes increasingly decoupled at larger scales. The topographic model of snow depth variation, which transforms terrain data to align with these scale-dependent snow patterns, is well correlated with local snow depth variations (linear fit R 2 > 0.5 for 85% of 100-m patches). A machine learning reconstruction of shrub canopy snow trapping reveals a simple exponential relationship between canopy structure and snow accumulation ( R 2 = 0.59), highlighting the combined influence of topography and vegetation on snow distribution. Together, these empirical relationships capture much of the observed snow variability in the watershed ( R 2 = 0.49, root mean square error (RMSE) = 30 cm), though systematic limitations persist in areas of strong scour and at coarser scales where wind-terrain interactions are more complex. These findings provide a framework for more efficient snow depth prediction and offer insights to improve snow-vegetation feedback representation in Earth System Models.

54 ENVIRONMENTAL SCIENCES↗

Ti-48Al-2Cr-2Nb Evaluated Under Fretting Conditions

Material parameters govern many of the design decisions in any engineering task. When two materials are in contact and microscopically small, relative motions (either vibratory or creeping) occur, and fretting fatigue can result. Fretting fatigue is a material response influenced by the materials in contact as well as by such variables as loading and vibratory conditions. Fretting produces fresh, clean interacting surfaces and induces adhesion, galling, and wear in the contact zone. Time, money, and materials are unnecessarily wasted when galling and wear result in excessive fretting fatigue that leads to poorly performing, unreliable mechanical systems. Fretting fatigue is a complex problem of significant interest to aircraft engine manufacturers. It can occur in a variety of engine components. Numerous approaches, depending on the component and the operating conditions, have been taken to address the fretting problems. The components of interest in this investigation were the low-pressure turbine blades and disks. The blades in this case were titanium aluminide, Ti-48Al-2Cr- 2Nb, and the disk was a nickel-base superalloy, Inconel 718 (IN 718). A concern for these airfoils is the fretting in fitted interfaces at the dovetail where the blade and disk are connected. Careful design can reduce fretting in most cases, but not completely eliminate it, because the airfoils frequently have a skewed (angled) blade-disk dovetail attachment, which leads to a complex stress state. Furthermore, the local stress state becomes more complex when the influence of the metal-metal contact and the edge of contact are considered.

Miyoshi, Kazuhisa↗