Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “variability analysis”

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,063 records · Page 59

Techno-Economic Analysis for Co-Processing Fast Pyrolysis Liquid in Fossil Refineries

Recent work at the National Renewable Energy Laboratory (NREL) and throughout the bioenergy community has highlighted incentives associated with co-processing bio-intermediates in existing refineries to reduce capital costs associated with renewable fuels and chemicals production and reduce overall risk for emerging biomass conversion technologies. This presentation summarizes the techno-economic analysis results from the NREL-Petrobras collaboration that focused on refinery integration of fast pyrolysis oil through co-processing in the fluid catalytic cracking (FCC) process. The NREL-Petrobras work highlights the economic opportunity for refiners to engage in co-processing for low-carbon products and introduces the risks for refiners based on variability in crude and fossil-product markets. The NREL-Petrobras results also serve as a basis to inform policy design that reduces economic risk for refinery co-processing and repurposing opportunities. The final topics in the presentation highlight experimental capabilities and emerging analysis approaches at NREL and partner laboratories that support the development and commercial deployment of refinery utilization strategies.

bio-fuels↗

Application of FARM to an IES scenario within the FORCE ecosystem

The FARM (Feasible Actuator Range Modifier) software module is a component of the RAVEN-based FORCE framework for analysis of Integrated Energy Systems (IES). FARM supports the HERON software module in the evaluation of the optimal dispatch by evaluating feasible set-points for the different IES unit components. Set-points are required to satisfy limits on both production variables (i.e., the variables to be optimized such as the electrical power, the hydrogen production rate, etc.) and process variables tied to the service life of equipment (e.g., steam flowrate, vessel pressure, turbine firing temperature, etc.). This problem is addressed by adopting a two-stage approach. First, the HERON power dispatcher determines set-points that meet the constraints on the former variables (e.g., power levels and power ramp rate limits). These constraints are called explicit constraints. Then, FARM adjusts these set-points to ensure the respect of the limits on the latter variables given the knowledge of the system physics acquired through machine learning algorithms. These constraints are called implicit constraints. From this standpoint, FARM constitutes a bridge between the HERON power dispatcher that adopts a simplified description of the IES unit (low-resolution physics) and the HYBRID high-fidelity models (high-resolution physics). The original version of the FARM software module (FARM-Alpha) was released by Argonne National Laboratory in January 2021. In the latest version of the code released in July 2022 (FARM-Delta), the Reference Governor (RG) algorithm was upgraded to a Multi-Input Multi-Output version from its original Single-Input Single-Output form. The RG algorithm acts to enforce constraints. With this improvement an IES unit is now treated as a single dynamic system from the standpoint of control. The crosstalk among components in an IES unit is now fully considered thereby ensuring a true optimization is obtained for those units that have multiple set-points. In this report, the capabilities of FARM-Delta operating within the FORCE ecosystem are demonstrated for an IES test case. The specific configuration of IES unit for this case was selected by the IES team with consultation from the Advanced Reactor IES Expert Group. A full TEA analysis that invoked HERON, HYBRID, FARM, and RAVEN was performed and serves to demonstrate how the latest modification to FARM algorithms (i.e., state variable selection, state-space matrices derivation, set-point verification) can shape setpoints that might otherwise compromise the health of equipment through accelerated wear and tear. In this specific test case, it was demonstrated that these algorithms ensure a more efficient utilization of steam resources to be shared by two different subsystems, namely Balance of Plant (BOP) and High-Temperature Steam Electrolysis (HTSE). Finally, some code improvements that can further enhance the user-friendliness are suggested.

97 MATHEMATICS AND COMPUTING↗

Reliability-Informed Life Cycle Warranty Cost and Life Cycle Analysis of Newly Manufactured and Remanufactured Units

Standard life cycle techniques such as life cycle warranty cost (LCWC) analysis and life cycle analysis (LCA) are used to respectively quantify the relative economical and environmental advantages of remanufactured goods while simultaneously identifying avenues for improvement. In this paper, we contribute to the literature on life cycle studies by incorporating reliability into LCWC analysis and LCA with the goal of improving long-term/multiple life cycle decision making. We develop a branched power-law model to incorporate the physical degradation mechanisms leading to reduced reuse rates of system parts over multiple life cycles. We then follow a standard LCA protocol to quantify the difference between a new unit and its remanufactured version in terms of environmental impact items such as abiotic depletion potential, global warming potential, and energy consumption. We then devise four practical warranty policies that vary in the choice of replacement and/or provision for extended warranty. All possible replacement scenarios for multiple life cycles are explored for each policy and a mathematically rigorous framework is provided, where the reliability information is used to calculate probabilistic LCWC and life cycle impact items. This reliability-informed LCWC analysis and LCA framework enables design engineers to compare design options and warranty policies by quantifying both economical and environmental impacts to aid in decision making. Although the framework is presented in a general form applicable to any engineered system, we demonstrate the utility of this framework by using a case study of an infinitely variable transmission used in agricultural equipment.

42 ENGINEERING↗

Advanced Method Optimization with Categorical and Constrained Continuous Parameters

Traditional approaches to analytical method optimization (e.g., univariate and “guess-and-check”) can be time-consuming, costly, and often fail to identify true optima within the parameter space. Previous work defined and implemented a generalized technique for method optimization for continuous method parameters, but a knowledge gap remains for the incorporation of categorical variables into these advanced method optimization schemes. This work presents and validates a generalized optimization approach that incorporates both continuous and categorical variables while also utilizing a multivariate, multiobjective optimization scheme with Karush–Kuhn–Tucker conditions to bound the optimization space to solutions within the physical limitations of the parameter space. Method optimization from a case study using GC–MS for the analysis of 11 analytical standards with objectives to minimize peak width and maximize peak height resulted in a 3 orders of magnitude improvement in the average peak height and a 2 orders of magnitude improvement in the average peak width compared to the least optimal (but reasonable) instrumental parameters utilized in this study. This approach to optimization allows for a customizable method optimization in which users can include both continuous and categorical variables to achieve objectives specific to their analytical goals. This approach significantly reduces the labor and cost associated with traditional method development approaches and can be applied in a variety of scientific fields across a range of laboratory techniques (e.g., instrument method development, sample preparation, and extraction techniques).

Amorphous materials↗

Dimuon production in the MicroBooNE detector

This analysis studies the production of the rare signal with two muons (dimuon) in the finalstate at µB using the Booster Neutrino Beam (BNB) as source of neutrino interactions. Themain background source for the dimuon signal are pions. Differentiating muons and pions hasbeen challenging at µB because the tracks left by these particles are almost identical. Thisanalysis faces this problem using a supervised classifier algorithm called boosted decisiontree (BDT). A Monte Carlo simulation (MC) enhanced dimuon signal sample was producedto train a BDT capable of recognizing this signal. A set of variables that show somediscrimination characteristics was used to perform the training. A MC neutrino interactionsimulation sample was used to estimate the background. Estimations of dimuon productionin real data samples were made for different cases. Additionally, an analysis of the systematicuncertainties, fake data studies, and BDT performance checks were obtained.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Simulation and analysis of slip flow of water at hydrophobic silica surfaces of nanometer slit pores

Water flow in nanopores is controlled fundamentally by pore wettability and interfacial water features. In the research herein reported, Molecular dynamics (MD) was used to simulate water flow in nanometer slit pores of silica surfaces with variable wetting characteristics. Compared to Couette flow, the simulated slit-pores of the talc (001) and octadecyltrichlorosilane (OTS) monolayer surfaces exhibited water flow enhancement and a positive slip length at the interface, whereas water flow velocity at the quartz (001) and amorphous silica surfaces was zero, adhering to the no slip condition. More intensive interactions between the interfacial water and hydrophilic silica surfaces were confirmed by MD simulation and sum-frequency vibrational spectroscopy (SFVS) results. A water exclusion zone of approximately 3 Å was found at the hydrophobic silica surfaces by quantitative MD interfacial water analysis and validated by SFVS. This exclusion zone dimension was similar in magnitude to the slip length (4–5 Å) as determined from simulated water flow in hydrophobic silica slit pores. The simulated slip lengths at the hydrophobic surfaces (0.2–1 nm) were significantly less than the experimental slip lengths (>40 nm) reported in the literature. In the case of the experimental results, this difference appears to be due to nanobubbles expected to be present at hydrophobic surfaces.

42 ENGINEERING↗

Performance and Commissioning of the BigBite Timing Hodoscope for Nucleon Form Factor Measurements at Jefferson Lab

The BigBite Timing Hodoscope detector is the primary subject of this thesis. The Super BigBite Spectrometer is a Jefferson Lab Hall A Collaboration project that has and will continue to measure nucleon electromagnetic form factors. This spectrometer includes the Timing Hodoscope which provides high resolution particle timing data for scattered electrons in the electron arm of BigBite. The Timing Hodoscope utilizes 90, 25 × 25 × 600 mm3 scintillator bars stacked on top of each other to form a single detector plane, and these bars are connected to 180 photo-multiplier tubes via light guides. Particles collide with the scintillating material creating a shower of optical photons and these particle events in the bars are collected to generate signals that are readout by the data acquisition (DAQ) electronics. NINO ASIC amplifier-discriminator cards output signals from the photo-multiplier tubes into analogue and logic signals, which are sent to analogue-to-digital (ADC) and time-to-digital (TDC) converter data acquisition readout modules. This data is then used for analysis of the detector. The focus of this thesis is the construction, commissioning, calibration, and performance of the BigBite Timing Hodoscope before and during the first of five nucleon electromagnetic form factor experiments at Jefferson Lab Hall A. Before the neutron magnetic form factor, G n M, experiment, cosmic ray data was collected during commissioning to confirm proper operation of the Timing Hodoscope electronics by observing the ADC and TDC data. Commissioning studies for charge normalization, gain matching, and other ADC and TDC detector data variables were performed before moving the detector into Hall A. Following installation in Hall A, several calibration studies were implemented to fine-tune the detector in preparation for use in the experiment. The calibration studies included analysis of timing cuts, TDC alignment, the time-walk effect, time difference offsets, and scintillator velocity corrections. Once the Timing Hodoscope was well-calibrated, data-taking during the experiment commenced and the beam-on-target data was used to characterize the Timing Hodoscope performance during the G n M experiment run-time. The performance analysis included studies observing energy deposit, cluster size, rates, accidentals, pile-up, tracking efficiency, position resolution, and time resolution. After application of physics cuts to ensure a data set comprised of particle tracks corresponding to elastic electrons, which is the main data of interest for measurement of G n M, the Timing Hodoscope is shown on average across all kinematic settings to have a >98% tracking efficiency, a position resolution of 4-6 cm in the non-dispersive plane and 1.5-2 cm in the dispersive plane, and a time resolution of 500-750 ps. These performance results are compared to a GEANT4 based performance simulation of the BigBite Timing Hodoscope for reference, showing to what degree the measured performance values match those taken from the simulation.

Marinaro, Ralph↗

Fault Detection via Occupation Kernel Principal Component Analysis

Reliable operation of automatic systems is heavily dependent on the ability to detect faults in the underlying dynamics. While traditional model-based methods have been widely used for fault detection, data-driven approaches have garnered increasing attention due to their ease of deployment and minimal need for expert knowledge. In this letter, we present a novel principal component analysis (PCA) method that uses occupation kernels. Occupation kernels result in feature maps that are tailored to the measured data, have inherent noise-robustness due to the use of integration, and can utilize irregularly sampled system trajectories of variable lengths for PCA. The occupation kernel PCA method is used to develop a reconstruction error approach to fault detection and its efficacy is validated using numerical simulations.

42 ENGINEERING↗

Seasonal forecasting skill for the High Mountain Asia region in the Goddard Earth Observing System

Seasonal variability of the global hydrologic cycle directly impacts human activities, including hazard assessment and mitigation, agricultural decisions, and water resources management. This is particularly true across the High Mountain Asia (HMA) region, where availability of water resources can change depending on local seasonality of the hydrologic cycle. Forecasting the atmospheric states and surface conditions, including hydrometeorologically relevant variables, at subseasonal-to-seasonal (S2S) lead times of weeks to months is an area of active research and development. NASA's Goddard Earth Observing System (GEOS) S2S prediction system has been developed with this research goal in mind. Here, we benchmark the forecast skill of GEOS-S2S (version 2) hydrometeorological forecasts at 1–3-month lead times in the HMA region, including a portion of the Indian subcontinent, during the retrospective forecast period, 1981–2016. To assess forecast skill, we evaluate 2 m air temperature, total precipitation, fractional snow cover, snow water equivalent, surface soil moisture, and terrestrial water storage forecasts against the Modern-Era Retrospective analysis for Research and Applications, Version 2 (MERRA-2) and independent reanalysis data, satellite observations, and data fusion products. Anomaly correlation is highest when the forecasts are evaluated against MERRA-2 and particularly in variables with long memory in the climate system, likely due to the similar initial conditions and model architecture used in GEOS-S2S and MERRA-2. When compared to MERRA-2, results for the 1-month forecast skill range from an anomaly correlation of R anom =0.18 for precipitation to R anom =0.62 for soil moisture. Anomaly correlations are consistently lower when forecasts are evaluated against independent observations; results for the 1-month forecast skill range from R anom =0.13 for snow water equivalent to R anom =0.24 for fractional snow cover. We find that, generally, hydrometeorological forecast skill is dependent on the forecast lead time, the memory of the variable within the physical system, and the validation dataset used. Overall, these results benchmark the GEOS-S2S system's ability to forecast HMA hydrometeorology.

54 ENVIRONMENTAL SCIENCES↗

Effects of input gradient regularization on neural networks time-series forecasting of thermal power systems

This study proposes using neural networks, specifically gated recurrent unit (GRU), long-short-term memory (LSTM), and transformer networks, to improve control strategies in a 450 MW coal-fired power plant. However, neural networks face issues of becoming overly dependent on just a few variables to make predictions, which negatively impacts control decisions that rely on the model to determine the value of all manipulated variables. The paper introduces regularization techniques, including noise injection and input gradient regularization, during the training phase. Here, the work presents novel contributions in adapting neural networks to control industrial systems and applying regularization techniques from computer vision to industrial process control. Results demonstrate the effectiveness of input gradient regularization in reducing model dependence on subsets of variables, emphasizing the balance between fidelity and controllability. Further exploration is recommended, including the development of recurrent transformers, closed-loop control testing, and a sensitivity analysis on computer models to provide further insight.

20 FOSSIL-FUELED POWER PLANTS↗

Ensemble cure kinetics network (ECK-Net): A method to derive cure kinetics of thermosetting resin

This paper introduces an Ensemble Cure Kinetics Network (ECK-Net), a neural network (NN)–based framework for modeling the cure kinetics of thermosetting resins within a phenomenological context. ECK-Net replaces traditional analytic models, which require extensive chemical insight and multiple isothermal/non-isothermal experiments, with a data-driven surrogate that maps nonlinear relationships between temperature, degree of cure, and reaction rate from differential scanning calorimetry data. The proposed approach predicts input-dependent kinetic coefficients of a generalized nth-order reaction equation rather than reaction rates directly, enabling a single unified model to represent various epoxy systems without relying on iso-conversional analysis or predefined functional forms. To ensure robustness, multiple independently trained networks under different random initializations are blended through an ensemble strategy, effectively mitigating the stochastic variability inherent to neural networks. The framework is validated using experimental datasets from multiple resin systems, including aerospace-grade materials (Toray 3900-2, Cycom 5320-1, and Hexcel 8552) and a windmill-grade resin (RIMR 035c). The model accurately reproduces the temporal evolution of the degree of cure under manufacturers’ recommended cure cycles across all tested resins systems, yielding Pearson’s correlation coefficients of 0.992, 0.994, 0.993, 0.997, respectively. To demonstrate process-level applicability, the trained network was implemented within the Abaqus environment to simulate out-of-autoclave (OOA) curing process of the CFRP panel composed of Toray T830H-6K/3900-2D prepreg. The simulation results showed excellent agreement with experimental temperature response (maximum peak temperature, simulation: 189.6 °C, experiment: 188.5 °C) and the final degree of cure (simulation: 0.948, experiment: 0.960 ± 0.013), confirming ECK-Net’s capability as a reliable alternative to conventional cure kinetics modeling methods.

Composite curing↗

Thermoelectric air conditioning undergarment for personal thermal management and HVAC energy saving

Personal thermal management systems (PTMS) are highly desirable for improving individual thermal comfort and reducing indoor HVAC energy consumption. Although there have been many attempts in developing PTMS, existing PTMS are generally bulky, heavy or immobile. Besides, previous analysis on additional personal cooling and heating power required for maintaining thermal comfort in changing environments did not consider the individual difference in metabolism and the variability of comfort skin temperature with environmental conditions. In this paper, we report on the development of a novel lightweight (<1 kg) thermoelectric air conditioning undergarment system with a preferred design of branching tubing network (branching angle of 60° and diameter ratio of 0.782) for air distribution. Through thermal manikin and human subject tests, we demonstrated that the novel system is capable of providing a maximum of 15.5 W of personal cooling and 18.1 W personal heating with a coefficient of performance (COP) greater than 0.4, which is sufficient to the expansion of the indoor set-points by at least 2.2 °C (4 °F) on both sides without compromising thermal comfort and the potential saving of about 15% HVAC energy. Furthermore, our study improved the understanding of the required additional personal cooling and heating power for changing environmental conditions, which are essential guidelines for the further development of PTMS.

42 ENGINEERING↗

A Satellite-Derived Upper-Ocean Stratification Data Set for the Tropical North Atlantic With Potential Applications for Hurricane Intensity Prediction

Upper-ocean stratification strongly impacts vertical mixing and the heat flux between the ocean and atmosphere, especially under extreme conditions of tropical cyclones (TCs). Knowledge of prestorm stratification is important for accurate TC intensity prediction. In situ observations of the tropical ocean have significantly increased in the past decade. However, they are still too sparse to resolve ocean stratification variability in near-real time and on small spatial scales. Here, in this study, based on long-term observations and an ocean reanalysis data set from 2004–2017, we investigate the possibility of retrieving upper-ocean stratification from sea surface temperature (SST), sea surface salinity (SSS), and sea surface height (SSH) using a simple regression method. It is found that more than 90% of the mean seasonal cycle and about 30% to 80% of temperature and salinity stratification anomalies can be reconstructed using surface data from either observations or an ocean reanalysis. Simple regression can be used with satellite observations to create a high-resolution, near-real-time-gridded ocean stratification data set that successfully reproduces both the large and mesoscale variability of ocean stratification. When used in a simple expression for TC-induced SST cooling, the satellite-derived stratification shows improvements over an ocean analysis in terms of variance explained of SST cooling, offering promise as a near-real-time indicator of the ocean's impact on TC intensification.

54 ENVIRONMENTAL SCIENCES↗

Hydrobiogeochemical Controls on the Delivery of Dissolved Organic Matter to Boreal Headwater Streams

Understanding controls on the delivery of dissolved organic matter (DOM) from terrestrial to aquatic systems is key for constraining carbon balances and fluxes in boreal headwater catchments. The largest export of DOM generally occurs during intense periods of water delivery to the landscape (e.g., rain events); however, the timing and magnitude of DOM is complicated by geomorphological, hydrometeorological, and biogeochemical variability. To better assess mechanisms controlling the delivery of DOM to boreal streams, DOM was investigated in two morphologically distinct catchments of an experimental forest watershed in western Newfoundland, one dominated by low relief wetlands and the other by steep hillslopes. The DOM was compared during two fall storms of similar magnitude but contrasting moisture conditions during the autumn transition. Variable concentrations and optical character highlighted antecedent conditions are important for the delivery of DOM to boreal streams. Concentration-discharge relationships with stream hydrograph separation analysis suggest preferential flowpaths through shallow mineral horizons as a key pathway for delivery of DOM in the hillslope dominated catchment compared to rapid input from near-stream wetlands in the low relief catchment. Loss of DOM via preferential pathways suggests an additional mechanism for surface soil carbon loss not directly informed by lateral flow and soil hydrology, which has important implications for our conceptualization and modeling of carbon fluxes in boreal systems. These results suggest watershed scale responses of DOM in boreal headwater catchments are complex, and better linkages between catchment scale hydrology and landscape hydrobiogeochemistry are required to constrain terrestrial to aquatic carbon fluxes in boreal landscapes.

54 ENVIRONMENTAL SCIENCES↗

Do automated vehicle (AV) enthusiasts value travel time differently from cautious travelers? an exploration of travelers’ attitudes towards AV

The anticipation of automated vehicles (AV) has resulted in diverse research topics, including understanding the adoption and usage of this emerging mode. Using attitudinal indicators, this paper identifies prevalent types of attitudes towards AV. It then explores the variation of those attitudes across different demographic groups and quantifies the influence of these attitudes on the preference for AV. The paper uses stated preference survey data collected among respondents of Chicago in 2019. The exploratory factor analysis of the attitudinal indicators identifies four factors, namely 'Cautious Travelers', 'Hedonists', 'Utilitarians', and 'Safety-First'. Next, the paper develops an integrated choice and latent variable model to test the influence of these four types of attitudes on the preference for AV. Finally, the results indicate that the cautious travelers associate higher Value of Travel Time (VOTT) to AV and enthusiasts associate lower VOTT to AV than a base model without such attitudes.

33 ADVANCED PROPULSION SYSTEMS↗

Evidence for multi-decadal fuel buildup in a large California wildfire from smoke radiocarbon measurements

In recent decades, there has been a significant increase in annual area burned in California's Sierra Nevada mountains. This rise in fire activity has prompted the need to understand how historical forest management practices affect fuel composition and emissions. Here we examined the total carbon (TC) concentration and radiocarbon abundance (Δ 14 C) of particulate matter (PM) emitted by the KNP Complex Fire, which occurred during California's 2021 wildfire season and affected several groves of giant sequoia trees in the southern Sierra Nevada. During a 26 h sampling period, we measured concentrations of fine airborne PM (PM 2.5 ), as well as dry air mole fractions of carbon monoxide (CO) and methane (CH 4 ), using a ground-based mobile laboratory. We also collected filter samples of PM 2.5 for analysis of TC concentration and Δ 14 C. High correlation among PM 2.5 , CO, and CH 4 time series confirmed that our PM 2.5 measurements captured variability in wildfire emissions. Using a Keeling plot approach, we determined that the mean Δ 14 C of PM 2.5 was 111.6 ± 7.7‰ (n = 12), which was considerably enriched relative to atmospheric carbon dioxide in the northern hemisphere in 2021 (-3.2 ± 1.4‰). Combining these Δ 14 C data with a steady-state one-box ecosystem model, we estimated that the mean age of fuels combusted in the KNP Complex Fire was 40 years, with a range of 29–57 years. These results provide evidence for emissions originating from woody biomass, larger-diameter fine fuels, and coarse woody debris that have accumulated over multiple decades. This is consistent with independent field observations that indicate high fire intensity contributed to widespread giant sequoia mortality. With the expanded use of prescribed fires planned over the next decade in California to mitigate wildfire impacts, our measurement approach has the potential to provide regionally-integrated estimates of the effectiveness of fuel treatment programs.

54 ENVIRONMENTAL SCIENCES↗

Quantifying the impact of climate change on erosion - 20397

The New York State Energy Research and Development Authority (NYSERDA) is the owner of the Western New York Nuclear Service Center (WNYNSC), a 1,351-ha (3,338-ac) site located approximately 48 km (30 mi) south of Buffalo, New York. In 1962, Nuclear Fuel Services, Inc. (NFS) entered into Agreements with the Atomic Energy Commission and New York State to construct the first commercial reprocessing plant of nuclear fuel in the United States. NFS, a private company, built and operated the spent fuel reprocessing plant and waste disposal facilities, processing 640 Mg (metric tons, or 705 short tons) of spent nuclear fuel from 1966 to 1972 under an Atomic Energy Commission license. Nuclear fuel reprocessing operations ended in 1972 and never reopened, leaving behind radioactive and chemical wastes. Erosion can play an important role in the fate and transport of waste at sites where disposal of long-lived waste is anticipated. The statistical characterization of key processes related to erosion is essential to the understanding of site stability through time. One of the key processes governing erosion is extreme precipitation and it is critical that trends in the distribution of extreme precipitation events through time be represented. The evidence of climate change is increasingly well documented and projected impacts on extreme precipitation events should be incorporated in performance assessment studies when relevant. Not evaluating future climate states in a performance assessment is contradictory to good modeling practice. Specifically, excluding climate change limits development of modeling information that could aid in effective decision making. The current climate literature provides both observational evidence and climate model projections of climate trends and/or climate change in the late 20. and early 21. centuries for North America and the northeast United States. In this work, this information was used to assess the impacts of potential changes in climate on erosion processes. The goal was to understand how projections of future climate relate to the performance of the WNYNSC through time. In the first stage of this work multi-temporal historical aerial images were analyzed in conjunction with orthophotography and Lidar data to develop probability distributions for variables representing important erosion processes. The information from these analyses was then used in conjunction with simulated data from the West Valley Erosion Working Group (EWG). Analysis of EWG simulations provides estimates for the change in erosion rates through time that is driven by changes in climate. The time-varying rates of erosion change were applied to the historical aerial imagery data in order to inform time-varying rates of erosion that are driven by changes in climate. Ultimately, this process resulted in the identification of locations for features like gully heads using both the Lidar dataset and projection of the estimated location from the historical aerial photo under consideration back to the Lidar dataset. The distance between the estimate of the location from the historical image and that of the Lidar dataset was the estimated distance the feature has moved. This distance was then divided by the number of years between the Lidar dataset and the year of the historical aerial image of interest to get a rate of movement through time. This was done for several dozen points on each historical aerial image. The analyses from the EWG were used to characterize the relative impact of climate change on the erosion rates. This relative impact was quantified by comparing the LEM based estimates of erosion that were derived using historical climate data with LEM based estimates of erosion that were derived using projections of future climate. The relative increase in the erosion rates was then applied to the historical estimates of erosion derived from the analysis of the historical aerial imagery. This approach was used since the LEM-based estimated of the historical erosion rates from the past century have a bias towards underestimating erosion. This underestimation is hypothesized to be a consequence of an inability of the LEMs to account for the impacts of land-use change (i.e. deforestation) which had been shown to significantly increase erosion in similar Northern hardwood forest ecosystems. In summary, gully head retreat rate and gully widening rate were characterized using statistical probability distributions from the historical aerial image analysis. Simulated data from the EWG was used to estimate the climatically-driven changes in erosion rates through time. The changes through time from the EWG were applied to the gully head retreat rate and gully widening rate characterized using the historical aerial image analysis. This information can be used to inform PPA models to project future risks from a given site. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗