Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “uncertainty 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 505 records · Page 28
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Statistical analysis of non-Maxwellian electron distribution functions measured with angularly resolved Thomson scattering

Angularly resolved Thomson scattering is a novel extension of Thomson scattering, enabling the measurement of the electron velocity distribution function over many orders of magnitude. Here, details of the theoretical basis of the technique and the instrument designed for this measurement are described. Angularly resolved Thomson-scattering data from several experiments are shown with descriptions of the corresponding distribution functions. A reduced model describing the distribution function is given and used to perform a Monte Carlo analysis of the uncertainty in the measurements. The electron density and temperature were determined to a precision of 12% and 21%, respectively, on average, while all other parameters defining the distribution function were generally determined to better than 20%. It was found that these uncertainties were primarily due to limited signal to noise and instrumental effects. Furthermore, measurements with this level of precision were sufficient to distinguish between Maxwellian and non-Maxwellian distribution functions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Towards a self-consistent analysis of the anisotropic galaxy two- and three-point correlation functions on large scales: application to mock galaxy catalogues

We establish a practical method for the joint analysis of anisotropic galaxy two- and three-point correlation functions (2PCF and 3PCF, respectively) on the basis of the decomposition formalism of the 3PCF using tripolar spherical harmonics. We perform such an analysis with MultiDark-Patchy mock catalogues to demonstrate and understand the benefit of the anisotropic 3PCF. We focus on scales above $80\, h^{-1}\, {\rm Mpc}$, and use information from the shape and the baryon acoustic oscillation (BAO) signals of the 2PCF and 3PCF. We also apply density field reconstruction to increase the signal-to-noise ratio of BAO in the 2PCF measurement, but not in the 3PCF measurement. In particular, we study in detail the constraints on the angular diameter distance and the Hubble parameter. We build a model of the bispectrum or 3PCF that includes the non-linear damping of the BAO signal in redshift space. We carefully account for various uncertainties in our analysis including theoretical models of the 3PCF, window function corrections, biases in estimated parameters from the fiducial values, the number of mock realizations to estimate the covariance matrix, and bin size. The joint analysis of the 2PCF and 3PCF monopole and quadrupole components shows a $30{{\ \rm per\ cent}}$ and $20{{\ \rm per\ cent}}$ improvement in Hubble parameter constraints before and after reconstruction of the 2PCF measurements, respectively, compared to the 2PCF analysis alone. Furthermore, this study clearly shows that the anisotropic 3PCF increases cosmological information from galaxy surveys and encourages further development of the modelling of the 3PCF on smaller scales than we consider.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Uncertainty Quantification of Calculated Temperatures for the AGR 5/6/7 Experiment

This report documents the quantification of uncertainty of the calculated temperature data for the Advanced Gas Reactor (AGR) 5/6/7 fuel irradiation experiment conducted in the Advanced Test Reactor at Idaho National Laboratory in support of the Advanced Reactor Technologies? research and development program. Recognizing uncertainties inherent in physics and thermal simulations of the AGR 5/6/7 capsules, the results of the numerical simulations are used in combination with statistical analysis methods to improve qualification of measured data. The calculated fuel temperatures for AGR tests are also used for validation of the fission product transport and fuel performance simulation models. These crucial roles of the calculated fuel temperatures in ensuring achievement of the AGR experimental program objectives require accurate determination of the model temperature uncertainties. This report covers temperature uncertainty results for each of the five AGR 5/6/7 capsules. To quantify the uncertainty of calculated temperatures determined using the ABAQUS finite element heat transfer code, this study identifies and analyzes model parameters of potential importance to the calculated temperatures of fuel compacts and thermocouples. The selection of input parameters for uncertainty quantification is based on the ranking of their influences upon temperature predictions. Thus, selected input parameters include those with high sensitivity and those with the largest uncertainty. Propagation of model parameter uncertainty and sensitivity is then used to quantify the overall uncertainty of calculated temperatures. Measurement uncertainty, analysis of modeling assumptions, and expert judgment are used as the basis to quantify the uncertainty range for selected input parameters. The input uncertainties are dynamic, accounting for the effect of unplanned events and changes in thermal properties of capsule components over extended exposure to high temperatures and fast neutron irradiation. The sensitivity analysis performed in this work went beyond the traditional local sensitivity. Using experimental design, analysis of pairwise interactions of model parameters was performed to establish sufficiency of the time dependent first order (linear) expansion terms in constructing the temperature response surface. To achieve completeness, uncertainty propagation made use of pairwise noise correlations of model parameters. Furthermore, using an interpolation scheme over the input parameter domain, the analysis obtains time dependent sensitivity over the test campaign duration. This allows computation of uncertainty for the calculated fuel temperatures and the calculated graphite temperatures at thermocouple locations during the entire irradiation period.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2020

Key Nuclear Data Impacting Reactivity in Advanced Reactors

Advanced reactor concepts currently being developed throughout the industry are significantly different from light water reactor (LWR) designs with respect to geometry, materials, and operating conditions, and consequently, with respect to their reactor physics behavior. Given the limited operating experience with non-LWRs, the accurate simulation of reactor physics and the quantification of associated uncertainties are critical for ensuring that advanced reactor concepts operate within the appropriate safety margins. Nuclear data are a major source of input uncertainties in reactor physics analysis. As part of an ongoing project at Oak Ridge National Laboratory (ORNL), the effects of nuclear data uncertainties on key figures of merit associated with advanced reactor safety are being assessed for selected advanced reactor technologies. Key nuclear data relevant for reactor safety analysis for each selected advanced reactor technology were identified, and their impact on important key figures of merit was assessed. Available advanced reactor specifications were reviewed, results from studies performed at ORNL and other research institutions were consulted, and available evaluated nuclear data libraries were analyzed. This report summarizes the key nuclear data for nuclides in the fuel, as well as other significant data, including scattering and neutron capture in various materials for the moderator, coolant, and structure of the considered advanced reactors. For the considered advanced reactors that use low-enriched uranium (LEU) fuel, results from LWR studies provided insight into relevant nuclear data given the lack of available studies specifically addressing these new systems. The major nominal missing data that were identified consist of thermal scattering data and 135m Xe cross section data for molten salt reactor (MSR) analysis. The identified major gaps with respect to nuclear data uncertainties are missing uncertainties of thermal scattering data for high temperature gas-cooled reactors and moderated MSR systems, and incomplete uncertainties on angular distributions in particular for fast spectrum systems, such as sodiumcooled fast reactors, fast molten salt reactors, and heat pipe reactors. Furthermore, it was found that special attention should be paid to cross section and uncertainty differences between different evaluated nuclear data library releases, because significant differences in nuclear data that can lead to major differences in reactivity calculations were found, even for well-known nuclides.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2022

Physics-informed machine learning with differentiable programming for heterogeneous underground reservoir pressure management

Abstract Avoiding over-pressurization in subsurface reservoirs is critical for applications like CO $$_2$$ 2 sequestration and wastewater injection. Managing the pressures by controlling injection/extraction are challenging because of complex heterogeneity in the subsurface. The heterogeneity typically requires high-fidelity physics-based models to make predictions on CO $$_2$$ 2 fate. Furthermore, characterizing the heterogeneity accurately is fraught with parametric uncertainty. Accounting for both, heterogeneity and uncertainty, makes this a computationally-intensive problem challenging for current reservoir simulators. To tackle this, we use differentiable programming with a full-physics model and machine learning to determine the fluid extraction rates that prevent over-pressurization at critical reservoir locations. We use DPFEHM framework, which has trustworthy physics based on the standard two-point flux finite volume discretization and is also automatically differentiable like machine learning models. Our physics-informed machine learning framework uses convolutional neural networks to learn an appropriate extraction rate based on the permeability field. We also perform a hyperparameter search to improve the model’s accuracy. Training and testing scenarios are executed to evaluate the feasibility of using physics-informed machine learning to manage reservoir pressures. We constructed and tested a sufficiently accurate simulator that is 400 000 times faster than the underlying physics-based simulator, allowing for near real-time analysis and robust uncertainty quantification.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Quantifying replication through repeated analysis of UVM-A, a liquid reference material for cosmogenic 10 Be and 26 Al studies

In this study, in-situ produced cosmogenic nuclide sample preparation and analysis is a complicated, multi-step process with numerous possible sources of error, many of which have not been robustly quantified. Here, we use a liquid reference material (UVM-A) to test whether the analytic precision associated with individual isotopic ratios (n = 96) measured at two Accelerator Mass Spectrometry facilities (Lawrence Livermore National Laboratory and Purdue Rare Isotope Measurement Laboratory) accurately represents the range of values determined for aliquots prepared in a single laboratory over a period of five years. We find that 10 Be/ 9 Be ratios measured at the different accelerator facilities (n = 26 and n = 70) have statistically indistinguishable central tendencies, matching to within 0.1 %. Based on this large dataset, we suggest preliminary consensus values for UVM-A of 1.45 ± 0.06 x 10 -13 for 10 Be/ 9 Be (mean, one standard deviation, n = 96) and 4.47 ± 0.22 x 10 -13 for 26 Al/ 27 Al (mean, one standard deviation, n = 27). For both 10 Be/ 9 Be and 26 Al/ 27 Al, the relative standard deviations of the measured ratios are similar to the reported analytic uncertainties of sample measurement, suggesting that sample preparation introduced little if any additional scatter beyond the uncertainty of isotopic analysis. This dataset demonstrates that robust sample preparation and measurement can generate isotope ratio data reproducible at the level of counting statistics.

58 GEOSCIENCES↗
DOE OSTI2022

Interpreting machine learning prediction of fire emissions and comparison with FireMIP process-based models

Annual burned areas in the United States have increased 2-fold during the past decades. With more large fires resulting in more emissions of fine particulate matter, an accurate prediction of fire emissions is critical for quantifying the impacts of fires on air quality, human health, and climate. This study aims to construct a machine learning (ML) model with game-theory interpretation to predict monthly fire emissions over the contiguous US (CONUS) and to understand the controlling factors of fire emissions. The optimized ML model is used to diagnose the process-based models in the Fire Modeling Intercomparison Project (FireMIP) to inform future development. Results show promising performance for the ML model, Community Land Model (CLM), and Joint UK Land Environment Simulator-Interactive Fire And Emission Algorithm For Natural Environments (JULES-INFERNO) in reproducing the spatial distributions, seasonality, and interannual variability of fire emissions over the CONUS. Regional analysis shows that only the ML model and CLM simulate the realistic interannual variability of fire emissions for most of the subregions (r >0.95 for ML and r =0.14~0.70 for CLM), except for Mediterranean California, where all the models perform poorly (r =0.74 for ML and r <0.30 for the FireMIP models). Regarding seasonality, most models capture the peak emission in July over the western US. However, all models except for the ML model fail to reproduce the bimodal peaks in July and October over Mediterranean California, which may be explained by the smaller wind speeds of the atmospheric forcing data during Santa Ana wind events and limitations in model parameterizations for capturing the effects of Santa Ana winds on fire activity. Furthermore, most models struggle to capture the spring peak in emissions in the southeastern US, probably due to underrepresentation of human effects and the influences of winter dryness on fires in the models. As for extreme events, both the ML model and CLM successfully reproduce the frequency map of extreme emission occurrence but overestimate the number of months with extremely large fire emissions. Comparing the fire PM 2.5 emissions from the ML model with process-based fire models highlights their strengths and uncertainties for regional analysis and prediction and provides useful insights into future directions for model improvements.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

AATTENUATION—The Atmospheric Attenuation Model for CSP Tower Plants: A Look-Up Table for Operational Implementation

Attenuation of solar radiation between the receiver and the heliostat field in concentrated solar power (CSP) tower plants can reduce the overall system performance significantly. The attenuation varies strongly with time and the average attenuation at different sites might also vary strongly from each other. If no site specific attenuation data is available, the optimal plant design cannot be determined and rough estimations of the attenuation effect are required leading to high uncertainties of yield analysis calculations. The attenuation is caused mainly by water vapor content and aerosol particles in the lower atmospheric layer above ground. Although several on-site measurement systems have been developed during recent years, attenuation data sets are usually not available to be included during the plant project development. An Atmospheric Attenuation (AATTENUATION) model to derive the atmospheric transmittance between a heliostat and receiver on the basis of common direct normal irradiance (DNI), temperature, relative humidity, and barometric pressure measurements was developed and validated by the authors earlier. The model allows the accurate estimation of attenuation for sites with low attenuation and gives an estimation of the attenuation for less clear sites. However, the site-dependent coefficients of the AATTENUATION model had to be developed individually for each site of interest, which required time-consuming radiative transfer simulations, considering the exact location and altitude, as well as the pre-dominant aerosol type at the location. This strongly limited the application of the model despite its typically available input data. In this manuscript, a look-up table (LUT) is presented which enables the application of the AATTENUATION model at the site of interest without the necessity to perform the according complex radiative transfer calculations for each site individually. This enables the application of the AATTENUATION model for virtually all resource assessments for tower plants and in an operational mode in real time within plant monitoring systems around the world. The LUT also facilitates the generation of solar attenuation maps on the basis of long-term meteorological data sets which can be considered during resource assessment for CSP tower plant projects. The LUTs are provided together with this manuscript as supplementary files. The LUT for the AATTENUATION model was developed for a solar zenith angle (SZA) grid of 1°, an altitude grid of 100 m, 7 different standard aerosol types and the standard AFGL atmospheres for mid-latitudes and the tropics. The LUT was tested against the original version of the AATTENUATION model at 4 sites in Morocco and Spain, and it was found that the additional uncertainty introduced by the application of the LUT is negligible. With the information of latitude, longitude, altitude above mean sea level, DNI, relative humidity (RH), ambient temperature (Tair), and barometric pressure (bp), the attenuation can be now derived easily for each site of interest.

14 SOLAR ENERGY↗
DOE OSTI2025

Summary of the 5th IAEA technical meeting on fusion data processing, validation and analysis (FDPVA)

The purpose of the 5th International Atomic Energy Agency technical meeting on fusion data processing, validation and analysis (FDPVA) (Ghent University, Ghent, Belgium, 12–15 June 2023) was to provide a platform during which a set of topics relevant to FDPVA were discussed with the view of meeting the needs of next step fusion devices such as ITER. The validation and analysis of experimental data obtained from diagnostics used to characterize fusion plasmas are crucial for a knowledge-based understanding of the physical processes governing the dynamics of these plasmas. This paper presents the recent progress and achievements in the domain of plasma diagnostics data analysis and synthetic diagnostics reported at the meeting, including concept description of new devices; fusion databases; integrated data analysis; inverse problems; uncertainty propagation, verification and validation; probabilistic methods and machine learning. The relevant results underline trends observed in the current major fusion confinement devices.

fusion databases↗
DOE OSTI2026

Lattice Relaxations in Bixbyite‐Structured Rare‐Earth Sesquioxides

This report examines atomic relaxations in rare‐earth sesquioxides with the cubic bixbyite crystal structure. In particular, the experimentally determined unit cell atom positions in bixbyite compounds are compared to “ideal” bixbyite atom positions, where the concept of ideality is based on similarities between the bixbyite and fluorite crystal structures. Analyses indicate that all rare‐earth sesquioxide compounds deviate from ideality in similar ways and that there is a trend wherein the magnitudes of these atomic deviations decrease with increasing cationic size. The deviations from an ideal, fluorite‐like structure produce highly distorted coordination polyhedra and variations in metal‐oxygen bond lengths. A modified bond valence sum analysis revealed that the first nearest neighbor cation‐anion bond length variations lead to large (20%) variations in cation‐anion bond strength. This report also provides a qualitative analysis regarding the accuracy of atomic relaxed deviation parameters, based on x‐ray diffraction measurements. The analysis indicates that uncertainties in the accuracy of measured oxygen relaxed deviation parameters in bixbyite sesquioxides are likely much greater than the published uncertainties.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Computationally inexpensive part-scale thermal history of additive friction-stir deposition

This study presents an analytical model for steady-state power generation and tool heat loss in additive friction-stir deposition (AFSD), developed to enable part-scale thermal simulation while remaining computationally inexpensive. The model predicts total generated power, yielding 3.7–4.7 kW across deposition temperature setpoints of 400–460 °C for the deposition of AA6061 with a Be-Cu tool. This corresponds to 90–95% of the reported spindle power. Tool heat loss is experimentally determined by calibrating a steady-state energy balance between the generated power, the substrate-deposition thermal gradient, and a temperature dependent tool heat loss term: q tool (T) = a + b (T - 400°C) with a = 2.7 x 10 6 Wm -2 and b = 9.5 x 10 3 Wm -2 K -1 . The calibration indicates that about 69% of the generated heat is conducted into the tool for this configuration, which is much higher than previously reported. The calibrated heat-source is implemented in finite element software (Adamantine) to simulate the transient thermal history of a 100 cm 3 representative build in 8 min on a standard desktop (at 0.635 mm build-height resolution). For the first three layers, the substrate temperatures between simulation and experiment are within 10% mean absolute percentage error. Sensitivity analysis indicates that uncertainties in average deposition temperature and deformation localization (stir-zone geometry, depth, and spatial dependance of strain-rate and flow stress) dominate model variance, motivating additional experimental verification.

Additive Friction-Stir Deposition↗
DOE OSTI2025

Revisiting a minimally destructive analytic approach for determining electrochemical kinetic parameters: Measuring aluminum corrosion across a wide pH range based on the Butler-Volmer equation

Here, this study revisits the three-point sampling of the simplified Butler-Volmer equation to address the limitations of strong potentiodynamic polarization, which can introduce irreversible damage and uncertainty in corrosion analysis. The method extracts electrochemical kinetic parameters while minimizing polarization effects, evaluates noise sensitivity relative to overpotential, and accounts for errors from signal noise, OCP drift, ohmic resistance, and mass-transfer constraints. Verified against the Tafel extrapolation method for aluminum corrosion across a wide pH range, this low-polarization approach enables accurate evaluations with specific error estimates, offering a robust alternative to linear polarization resistance methods that assume constant Tafel slopes.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Detailed temperature diagnostic studies for Radishock and COAX experiments

Spectral diagnostics provide a powerful probe of high energy-density physics experiments. By shining an x-ray source on a target, absorption features can be used to determine accurate temperature profiles of that target material. Many studies produce a single temperature/density measurement by fitting these observed spectra. This paper demonstrates how, by leveraging detailed simulations, we can not only measure the average temperature and density, but the full density and temperature profiles. To do so, we must conduct a careful analysis of the uncertainties in the diagnostic measurement. Here we discuss the characteristics and associated uncertainties of the spectral diagnostic used in the COAX, Radishock and OuTi experiments, ultimately demonstrating how these detailed studies increase the potential of this powerful probe.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

A Primer on Dose-Response Data Modeling in Radiation Therapy

An overview of common approaches used to assess a dose response for radiation therapy–associated endpoints is presented, using lung toxicity data sets analyzed as a part of the High Dose per Fraction, Hypofractionated Treatment Effects in the Clinic effort as an example. Each component presented (eg, data-driven analysis, dose-response analysis, and calculating uncertainties on model prediction) is addressed using established approaches. Specifically, the maximum likelihood method was used to calculate best parameter values of the commonly used logistic model, the profile-likelihood to calculate confidence intervals on model parameters, and the likelihood ratio to determine whether the observed data fit is statistically significant. The bootstrap method was used to calculate confidence intervals for model predictions. Correlated behavior of model parameters and implication for interpreting dose response are discussed.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2026

Creating high-precision reference gas standards of 85Kr for groundwater age-dating

Absolute gas counting (AGC) was applied to two gas blends of 85Kr in argon-methane (P10) counting gas to establish a high-precision specific activity (Bq/cm3) reference value for characterizing 85Kr detection efficiency for groundwater age dating measurements. The AGC or length-compensated technique has been utilized by the metrology community for decades and is an accepted method for developing radioactive gas standards. The AGC capability at Pacific Northwest National Laboratory (PNNL) uses a set of nine unequal-length proportional counters with precisely-measured internal volumes, and a gas loading system with high-precision pressure and temperature sensors. A series of AGC measurements were collected at multiple pressures to determine the inverse pressure relationship (1/P) for 85Kr and define a wall-effect correction that accounts for events decaying into the detector wall and not depositing sufficient energy in the gas to be detected. In addition to the wall-effect, two additional corrections were evaluated and are discussed in detail. Specifically, the threshold effect which accounts for events deposited below the analysis threshold and a detection efficiency as a function of detector volume effect that was observed during analysis. A robust uncertainty model was developed using the Guide to the expression of Uncertainty in Measurements (GUM) approach. The combination of carefully scrutinized correction factors, precise measurements of pressure, temperature and detector volume, and robust counting statistics resulted in the determination of high-precision specific activity values with 0.50% or less total combined uncertainty for two Kr-in-P10 reference gas standards (KP10) that will enable new groundwater age-dating measurements at PNNL.

Kr-85↗
DOE OSTI2023

Contextually aware roadside radiation measurement testbed

Here we demonstrate a contextually aware multimodal roadside radiation measurement detection testbed for traffic monitoring applications in nuclear nonproliferation. Many variables in traffic such as vehicle or cargo size, mass, speed, shape, and distance of closest approach can have significant impacts on the radiation measured from a vehicle-transported radiation source. These factors can lead to uncertainties in the analysis of the radiation source, especially for lower-strength radiation sources of interest. Our testbed, known as the Multimodal Measurement System (MMS) uses non-radiation sensors including magnetometers, geophones, radiofrequency receivers, cameras, and LiDAR to extract contextual information about vehicles passing by the system. These contextual data can then be fused with data from radiation measurements to increase the system’s sensitivity and accuracy in nuclear threat detection applications. This work describes the instrumentation of the MMS and its data acquisition pipeline. Furthermore, we describe the pre-analysis performed on the raw multimodal data streams for data fusion, and the high-level machine learning analyses for detection and characterization. The variety of sensors within the MMS provides a valuable testbed that can be used to identify the combinations of contextual sensors that provide the greatest improvements to radiation source detection and characterization within the restrictions for various proliferation detection applications. The MMS is also modular so that additional combinations of sensors can be explored in the future.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2026

A unifying equation for fermentation sustainability across the titer-rate-yield landscape

Industrial fermentation is central to the sustainable production of fuels and chemicals, yet commercial viability of emerging technologies hinges on improving fermentation titer, rate, and yield (TRY). How these metrics shape system cost remains difficult to generalize due to complex interactions among feedstocks, fermentation, separations, catalytic upgrading, waste management, and facility design. Here, we systematically map theoretical fermentation performance spaces (formed by all potential TRY combinations) for 32 representative biomanufacturing facilities—spanning distinct choices for feedstocks, fermentation regimes and products, separations, and catalytic upgrading—by simulating and evaluating them (via techno-economic analysis, TEA) under uncertainty (600,000 Monte Carlo simulations) and across TRY combinations (7500 TRY combinations for each of 32 configurations). Across this wide design and thermodynamic simulation space, we find the relationship between fermentation TRY and system cost is captured by a simple, generalizable mathematical equation (R 2 of 0.992 − 1.000 across our simulations; 0.954 − 1.000 when validated against prior studies that used different tools). We use this equation to elucidate key drivers that shape cost sensitivity to fermentation performance, generating widely applicable insights. By demonstrating a unifying relationship governs the impact of fermentation on biomanufacturing economics, this work establishes a foundation for agile, holistically predictive, resource-efficient strategies to prioritize fermentation research and development needs and accelerate commercialization of emerging biomanufacturing technologies.

applied mathematics↗
DOE OSTI2022

Projected U.S. drought extremes through the twenty-first century with vapor pressure deficit

Global warming is expected to enhance drought extremes in the United States throughout the twenty-first century. Projecting these changes can be complex in regions with large variability in atmospheric and soil moisture on small spatial scales. Vapor Pressure Deficit (VPD) is a valuable measure of evaporative demand as moisture moves from the surface into the atmosphere and a dynamic measure of drought. Here, VPD is used to identify short-term drought with the Standardized VPD Drought Index (SVDI); and used to characterize future extreme droughts using grid dependent stationary and non-stationary generalized extreme value (GEV) models, and a random sampling technique is developed to quantify multimodel uncertainties. The GEV analysis was performed with projections using the Weather Research and Forecasting model, downscaled from three Global Climate Models based on the Representative Concentration Pathway 8.5 for present, mid-century and late-century. Results show the VPD based index (SVDI) accurately identifies the timing and magnitude short-term droughts, and extreme VPD is increasing across the United States and by the end of the twenty-first century. The number of days VPD is above 9 kPa increases by 10 days along California’s coastline, 30–40 days in the northwest and Midwest, and 100 days in California’s Central Valley.

54 ENVIRONMENTAL SCIENCES↗