Search NASA⌕ Search

SEARCH · Search NASA

Results for “Bayesian sensitivity 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 145 records · Page 8
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

High-z galaxies with JWST and local analogues – it is not only star formation

ABSTRACT I present an analysis of the JWST NIRSpec data of SMACS 0723 released as Early Release Observations. As part of this three new redshifts are provided, bringing the total of reliable redshifts to 14. I propose a modification to the direct abundance determination method that reduces sensitivity to flux calibration uncertainties by a factor of ∼3 and show that the resulting abundances are in good agreement with Bayesian photoionization models of the rest-frame optical spectrum. I also show that 6355 is most likely a narrow-line active galactic nucleus (AGN) with $M_*\lt 10^9\, \mathrm{M}_\odot$ at z = 7.66 and argue that 10 612 might also have an AGN contribution to its flux through comparison to photoionization models and low-redshift analogues. Under the assumption that the lines come from star formation, I find that the galaxies have gas depletion times of ∼107 yr, comparable to similar galaxies locally. I also identify a population of possibly shock-dominated galaxies at z < 3, whose near-IR emission lines plausibly come nearly all from shocks and discuss their implications. I close with a discussion of the potential for biases in the determination of the mass–metallicity relation using samples defined by detected [O iii]4363 and show using low-z galaxies that this can lead to biases of up to 0.5 dex with a systematic trend with mass.

Brinchmann, Jarle (ORCID:0000000343598797)↗
DOE OSTI2025

Bayesian inference of nuclear-matter density from proton scattering

Background: Proton elastic scattering at intermediate energy is widely employed as a tool for determining the matter radius of atomic nuclei. Here, the sensitivity of the approach relies on high-resolution measurements at small scattering angles and low-momentum transfer. Under these conditions, the Glauber multiple scattering theory accurately describes the proton-nucleus elastic cross section. Purpose: Investigate the sensitivity of the Glauber multiple scattering theory to uncertainties associated with input parameters such as the nuclear-matter density distribution and nucleon-nucleon data. Method: A joint Bayesian inference was performed using 12 angular distributions of elastic scattering at different energies on 58 Ni, 90 Zr, and 208 Pb targets. A Metropolis-Hastings algorithm was implemented to make an uncertainty quantification analysis for the input parameters used in the Glauber multiple scattering theory. Results: The experimental cross sections were fitted simultaneously using a joint Bayesian inference approach. Posterior probability density distributions of 42 input parameters were obtained from the analysis. A moderate correlation between the nuclear density parameters and the nucleon-nucleon cross sections was found. This correlation impacts the extraction of the nuclear-matter radius. Conclusions: The present analysis provided a consistent method for extracting the nuclear-matter density distribution of 58 Ni, 90 Zr, and 208 Pb from data across different incident energies. Due to the correlation of the nucleon-nucleon cross sections with the other input parameters, a constrained Bayesian inference using free nucleon-nucleon cross section data was performed. The nuclear-matter radii obtained from the analysis are in good agreement with multiple results reported in the literature.

190 ≤ A ≤ 219↗
DOE OSTI2018

Bayesian classification of falls risk

Background. Prior research in falls risk prediction often relies on qualitative and/or clinical methods. There are two challenges with these methods. First, qualitative methods typically use falls history to determine falls risk. Second, clinical methods do not quantify the uncertainty in the classification decision. In this paper, we propose using Bayesian classification to predict falls risk using vectors of gait variables shown to contribute to falls risk. Research Questions: (1) Using a vector of risk ratios for specific gait variables shown to contribute to falls risk, how can older adults be classified as low or high falls risk? and (2) how can the uncertainty in the classifier decision be quantified when using a vector of gait variables? Methods. Using a pressure sensitive walkway, biomechanical measurements of gait were collected from 854 adults over the age of 65. In our method, we first determine low and high falls risk labels for vectors of risk ratios using the k-means algorithm. Next, the posterior probability of low or high falls risk class membership is obtained from a two component Gaussian mixture model (GMM) of gait vectors, which enables risk assessment directly from the underlying biomechanics. We classify the gait vectors using a threshold based on Youden's J statistic. Results. Through a Monte Carlo simulation and an analysis of the receiver operating characteristic (ROC), we demonstrate that our Bayesian classifier, when compared to the k-means falls risk labels, achieves an accuracy greater than 96% at predicting low or high falls risk. Significance. Our analysis indicates that our approach based on a Bayesian framework and an individual's underlying biomechanics can predict falls risk while quantifying uncertainty in the classification decision.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2025

Demonstrate new plasticity models for doped UO 2 that capture dislocation mechanisms

In light water reactors, fuel vendors are investigating the use of dopants to modify the properties of UO 2 pellets, with the goal of improving pellet-cladding mechanical interactions during operation. Dopants are expected to ‘soften’ the pellets; that is, the doped pellets have higher plastic deformation than conventional UO 2 . This leads to a reduction in the severity of mechanical pellet-cladding interactions, helping to reduce the hoop strain on the cladding. By minimizing the strain exerted by the pellet on the cladding, it is anticipated that cladding performance under accident conditions can be enhanced (i.e., lowering the risk of burst during a LOCA). Dopants such as chromium (Cr) promote grain growth during pellet fabrication, leading to larger grains; therefore, understanding the link between chemistry, microstructure and mechanical deformation (enhanced creep rates) behavior of UO 2 is critical to helping operators further substantiate the benefits of doping UO 2 . Historically, the nuclear energy industry has relied on empirical models to make assessments of performance. Compared to empirical models, mechanistic physics-based models provide benefits, such as, fewer data points for validation and better extrapolation where experimental data is scarce or non-existent. In this report, Bayesian inference techniques have been applied to a previously developed lower length-scale-informed diffusional creep model. The objective is to i) infer lower-length-scale parameter distributions from available experiment and then ii) determine the uncertainties in the measurable quantity (in this case creep rates) after propagating the inferred lower length scale parameter uncertainties. The approach requires many evaluations of the model, which becomes computationally insurmountable; therefore, a neural-network model is trained to data obtained by sampling the full model over the most important parameters. This neural-network is then used in the Bayesian inference approach to determine probability distributions in the parameter values that represent the uncertainty in the model given what is known from the experiments (posterior). A significant reduction compared to conservative initial (prior) uncertainties is achieved through inference against the experimental data, demonstrating the efficacy of this approach. Furthermore, by accounting for uncertainties in the experimental conditions and sample non-stoichiometry, it is possible to resolve apparent discrepancies in experimental measurements within a self-consistent grain boundary (Coble) creep model that is sensitive to chemistry. This work has been written up and submitted to Nuclear Technology for a special issue on accelerated fuel qualification (AFQ). This uncertainty quantification (UQ) work not only improves the diffusional model, while accounting for uncertainty, but also establishes a framework which can readily be applied to the mechanistic models of dislocation deformation developed in this study. The most likely values from the Bayesian analysis are incorporated into our UO 2 diffusional creep model and a lower length scale-informed irradiation UO 2 creep mechanistic model to generate a dataset. This dataset has been provided to our INL collaborators for training an artificial neural network surrogate model, which will be implemented in the BISON fuel performance code to assess how the results differ from those currently obtained using a fully empirical model and that of using the nominal (uncalibrated) atomic scale parameters in our mechanistic model. Plastic deformation (creep and glide) in UO 2 is a complex phenomenon, governed by multiple underlying processes such as local defect concentrations, applied stresses, and microstructural characteristics. Consequently, there is a need for a meso-scale model with polycrystalline resolution capable of extrapolating to large grain sizes applicable to doped UO 2 , where data is limited and the model can help bridge the knowledge gap. By integrating atomistic data into the polycrystal LApx code, it becomes possible to predict dislocation climb and glide plasticity that simple analytical models cannot accurately represent. The application of atomic-scale data within LApx demonstrated the importance of climb and glide mechanisms in reproducing high-stress UO 2 behavior. Behaviors such as this are crucial to capture and implement in BISON, as parts of the fuel pellet can reach temperatures where glide can occur before pellet cracking. This model which captures dislocation based mechanisms for UO 2 is then used to stand up the doped model accounting for larger grain sizes. It was found that larger grain sizes can lead to enhanced deformation rates in the glide regime, and therefore can help with the pellet cladding mechanical interaction. Therefore if the fuel pellet reaches conditions (stress/temperature) where glide is active, the enhanced creep rates for larger grains in the glide regime (doped UO 2 ) can help with pellet cladding mechanical interactions. Plastic deformation in UO 2 involves multiple mechanisms, including diffusional creep, dislocation climb, and glide. This milestone contains two parts: (1) UQ of a pre-existing lower length scale informed mechanistic diffusional creep model, and (2) development of a new LApx based model for dislocation-mediated creep mechanisms in UO 2 , with application to large-grain doped UO 2 .

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2020

Revisiting Bayesian constraints on the transport coefficients of QCD

Multistage models based on relativistic viscous hydrodynamics have proven successful in describing hadron measurements from relativistic nuclear collisions. These measurements are sensitive to the shear and the bulk viscosities of QCD and provide a unique opportunity to constrain these transport coefficients. Bayesian analyses can be used to obtain systematic constraints on the viscosities of QCD, through methodical model-to-data comparisons. In this manuscript, we discuss recent developments in Bayesian analyses of heavy ion collision data. Furthermore, we highlight the essential role of closure tests in validating a Bayesian analysis before comparison with measurements. We discuss the role of the emulator that is used as proxy for the multistage theoretical model. We use an ongoing Bayesian analysis of soft hadron measurements by the JETSCAPE Collaboration as context for the discussion.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2026

A Framework for Parametric and Predictive Uncertainty Quantification in the E3SM Land Model: Assessing Site and Observable Generalizability

Quantifying parametric uncertainty using observations from individual sites provides a critical foundation for Earth system modeling, serving as a necessary first step before scaling up to regional or global applications. This study introduces a novel computational framework designed to enhance model predictability by reducing parametric uncertainty and assessing site and observable generalizability using various observational constraints. The framework integrates five components: Model Simulation, Statistical Emulation, Global Sensitivity Analysis (GSA), Model Calibration, and Model Prediction. Using the E3SM land model, we simulated site-level land-atmosphere carbon and energy fluxes from 2003 to 2007 across five evergreen needleleaf FLUXNET sites, perturbing 26 vegetation-related model parameters. Gaussian process emulators were employed to expedite GSA and model calibration. Four critical parameters that strongly influence selected land-atmosphere fluxes were identified by GSA. Bayesian approaches were used to infer parameter probability distributions leveraging synthetic data and FLUXNET observations. The results reveal that posterior parameter distributions vary significantly across different sites and observables within the same plant functional type. Probabilistic predictions indicate that parameters calibrated at one site can enhance predictive accuracy at other sites, although site heterogeneity may sometimes outweigh parametric uncertainty. Additionally, the probabilistic predictions demonstrate that calibration for one variable can also improve predictability for other variables, thereby maximizing predictive capabilities with limited observations. This framework provides a powerful approach for reducing parametric uncertainty in Earth system models and deepening our understanding of carbon dynamics and energy cycles. Its adaptability makes it a valuable tool for broader applications in Earth system modeling.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Improved Limit on Neutrinoless Double-Beta Decay in Te 130 with CUORE

We report new results from the search for neutrinoless double-beta decay in $^{130}$Te with the CUORE detector. This search benefits from a four-fold increase in exposure, lower trigger thresholds and analysis improvements relative to our previous results. We observe a background of $(1.38\pm0.07)\cdot10^{-2}$ counts$/($keV$\cdot$kg$\cdot$yr$)$ in the $0\nu\beta\beta$ decay region of interest and, with a total exposure of 372.5 kg$\cdot$yr, we attain a median exclusion sensitivity of $1.7\cdot10^{25}$ yr. We find no evidence for $0\nu\beta\beta$ decay and set a $90\%$ CI Bayesian lower limit of $3.2\cdot10^{25}$ yr on the $^{130}$Te half-life for this process. In the hypothesis that $0\nu\beta\beta$ decay is mediated by light Majorana neutrinos, this results in an upper limit on the effective Majorana mass of 75-350 meV, depending on the nuclear matrix elements used.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2020

Evaluation of a pair-based, joint-likelihood association approach for regional infrasound event identification

SUMMARY A Bayesian framework for the association of infrasonic detections is presented and evaluated for analysis at regional propagation scales. A pair-based, joint-likelihood association approach is developed that identifies events by computing the probability that individual detection pairs are attributable to a hypothetical common source and applying hierarchical clustering to identify events from the pair-based analysis. The framework is based on a Bayesian formulation introduced for infrasonic source localization and utilizes the propagation models developed for that application with modifications to improve the numerical efficiency of the analysis. Clustering analysis is completed using hierarchical analysis via weighted linkage for a non-Euclidean distance matrix defined by the negative log-joint-likelihood values. The method is evaluated using regional synthetic data with propagation distances of hundreds of kilometres in order to study the sensitivity of the method to uncertainties and errors in backazimuth and time of arrival. The method is found to be robust and stable for typical uncertainties, able to effectively distinguish noise detections within the data set from those in events, and can be made numerically efficient due to its ease of parallelization.

58 GEOSCIENCES↗
DOE OSTI2021

O 16 O 16 collisions at energies available at the BNL Relativistic Heavy Ion Collider and at the CERN Large Hadron Collider comparing α clustering versus substructure

Collisions of . light and heavy nuclei in relativistic heavy-ion collisions have been shown to be sensitive to nuclear structure. With a proposed 16 O 16 O run at the CERN Large Hadron Collider (LHC) and at the BNL Relativistic Heavy Ion Collider (RHIC) we study the potential for finding α clustering in 16 O. Here we use the state-of-the-art iEBE-VISHNU package with 16 O nucleonic configurations from ab initio nuclear lattice simulations. This setup was tuned using a Bayesian analysis on p Pb and PbPb systems. We find that the 16 O 16 O system always begins far from equilibrium and that at LHC and RHIC it approaches the regime of hydrodynamic applicability only at very late times. Finally, by taking ratios of flow harmonics we are able to find measurable differences between α-clustering, nucleonic, and subnucleonic degrees of freedom in the initial state.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2018

Interatomic Potentials Models for Cu-Ni and Cu-Zr Alloys

This study explores a Bayesian calibration framework for the RAMPAGE alloy potential model for Cu-Ni and Cu-Zr systems, respectively. In RAMPAGE potentials, it is proposed that once calibrated potentials for individual elements are available, the inter-species interactions can be described by fitting a Morse potential for pair interactions with three parameters, while densities for the embedding function can be scaled by two parameters from the elemental densities. Global sensitivity analysis tools were employed to understand the impact each parameter has on the MD simulation results. A transitional Markov Chain Monte Carlo algorithm was used to generate samples from the multimodal posterior distribution consistent with the discrepancy between MD simulation results and DFT data. For the Cu-Ni system the posterior predictive tests indicate that the fitted interatomic potential model agrees well with the DFT data, justifying the basic RAMPAGE assumptions. For the Cu-Zr system, where the phase diagram suggests more complicated atomic interactions than in the case of Cu-Ni, the RAMPAGE potential captured only a subset of the DFT data. The resulting posterior distribution for the 5 model parameters exhibited several modes, with each mode corresponding to specific simulation data and a suboptimal agreement with the DFT results.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2018

Current and future constraints on Higgs couplings in the nonlinear Effective Theory

Here, we perform a Bayesian statistical analysis of the constraints on the nonlinear Effective Theory given by the Higgs electroweak chiral Lagrangian. We obtain bounds on the effective coefficients entering in Higgs observables at the leading order, using all available Higgs-boson signal strengths from the LHC runs 1 and 2. Using a prior dependence study of the solutions, we discuss the results within the context of natural-sized Wilson coefficients. We further study the expected sensitivities to the different Wilson coefficients at various possible future colliders. Finally, we interpret our results in terms of some minimal composite Higgs models.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2012

Measurement of the electron antineutrino mass in tritium beta decay in the Troitsk nu-mass experiment

The results obtained in the Troitsk nu-mass experiment by measuring the electron-antineutrino mass in tritium beta decay are presented. The facility used consists of a gaseous windowless tritium source and an electrostatic electron spectrometer involving an adiabatic magnetic collimation. Runs in which measurement conditions were reliably established were thoroughly selected in analyzing data obtained from 1994 to 2004. All known systematic effects were taken into account. For the square of the electron-antineutrino mass, the treatment of measured spectra yielded the following result: m{sub {nu}}{sup 2} = -0.67 {+-} 1.89{sub stat.} {+-} 1.68{sub syst.} eV{sup 2}. The use of the Bayesian method and the Feldman-Cousins unified approach made it possible to obtain the following upper limits on the mass: m{sub {nu}} < 2.12 eV (at a 95% C.L.; Bayesian method) and m{sub {nu}} < 2.05 eV (at a 95% C.L., Feldman-Cousins method). At the same time, an estimation of the sensitivity limit without allowance for negative values of the square of the mass leads to m{sub {nu}} < 2.2 eV (at a 95% C.L.). Measured spectra were analyzed for the possible existence of an additional structure (step) in the electron spectrum near the boundary energy. The conclusion drawn from this analysis was that, within the existing statistical errors, there are no reasons for introducing such a feature.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2014

MOA-2011-BLG-262Lb: A sub-Earth-mass moon orbiting a gas giant primary or a high velocity planetary system in the galactic Bulge

We present the first microlensing candidate for a free-floating exoplanet-exomoon system, MOA-2011-BLG-262, with a primary lens mass of M {sub host} ∼ 4 Jupiter masses hosting a sub-Earth mass moon. The argument for an exomoon hinges on the system being relatively close to the Sun. The data constrain the product M{sub L} π{sub rel} where M{sub L} is the lens system mass and π{sub rel} is the lens-source relative parallax. If the lens system is nearby (large π{sub rel}), then M{sub L} is small (a few Jupiter masses) and the companion is a sub-Earth-mass exomoon. The best-fit solution has a large lens-source relative proper motion, μ{sub rel} = 19.6 ± 1.6 mas yr{sup –1}, which would rule out a distant lens system unless the source star has an unusually high proper motion. However, data from the OGLE collaboration nearly rule out a high source proper motion, so the exoplanet+exomoon model is the favored interpretation for the best fit model. However, there is an alternate solution that has a lower proper motion and fits the data almost as well. This solution is compatible with a distant (so stellar) host. A Bayesian analysis does not favor the exoplanet+exomoon interpretation, so Occam's razor favors a lens system in the bulge with host and companion masses of M{sub host}=0.12{sub −0.06}{sup +0.19} M{sub ⊙} and m{sub comp}=18{sub −10}{sup +28} M{sub ⊕}, at a projected separation of a{sub ⊥}=0.84{sub −0.14}{sup +0.25} AU. The existence of this degeneracy is an unlucky accident, so current microlensing experiments are in principle sensitive to exomoons. In some circumstances, it will be possible to definitively establish the mass of such lens systems through the microlensing parallax effect. Future experiments will be sensitive to less extreme exomoons.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2012

Global Fits of the Minimal Universal Extra Dimensions Scenario

In theories with Universal Extra-Dimensions (UED), the {gamma}{sub 1} particle, first excited state of the hypercharge gauge boson, provides an excellent Dark Matter (DM) candidate. Here we use a modified version of the SuperBayeS code to perform a Bayesian analysis of the minimal UED scenario, in order to assess its detectability at accelerators and with DM experiments. We derive in particular the most probable range of mass and scattering cross sections off nucleons, keeping into account cosmological and electroweak precision constraints. The consequences for the detectability of the {gamma}{sub 1} with direct and indirect experiments are dramatic. The spin-independent cross section probability distribution peaks at {approx} 10{sup -11} pb, i.e. below the sensitivity of ton-scale experiments. The spin-dependent cross-section drives the predicted neutrino flux from the center of the Sun below the reach of present and upcoming experiments. The only strategy that remains open appears to be direct detection with ton-scale experiments sensitive to spin-dependent cross-sections. On the other hand, the LHC with 1 fb{sup -1} of data should be able to probe the current best-fit UED parameters.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

Investigating Kinetic Mechanisms of Soot Formation in Plasma Pyrolysis of Methane via Active Learning (Final Technical Report)

Plasma pyrolysis of methane is an effective route for zero-carbon hydrogen production. Yet, soot generated from pyrolysis of hydrocarbons is detrimental to the climate and human health. There is ample experimental and theoretical evidence that suggests polycyclic aromatic hydrocarbons (PAHs) are the molecular precursors to soot particles. The reaction pathways of PAH formation are intricately dependent on a multitude of process parameters, whose kinetic mechanisms are not well-understood in plasma pyrolysis. This project aims to leverage advances in the kinetic modeling of soot formation in combustion, as well as in surrogate modeling and active learning, to systematically investigate the effects of process parameter on the kinetics of PAH formation in plasma pyrolysis of methane. To this end, we propose to use the PAH formation kinetics model developed by the PPPL/PU group based on the well-established ABF and HACA mechanisms, coupled with low-temperature plasma models. We will develop an active learning (AL) framework based on Bayesian optimization to systematically and data-efficiently explore the complex and multivariable parameter space of plasma pyrolysis in order to quantify the effects of plasma and feed parameters on the ABF and HACA kinetic pathways. AL is the branch of machine learning concerned with systematically querying samples from a system (experimental or computational) to train a data-driven model that maps design parameters to a performance criterion. We will use the data generated via AL to perform global sensitivity analysis, combined with uncertainty quantification, to elucidate the impact of different reaction pathways on minimizing formation of soot precursors. This study will result in an improved understanding of kinetics of PAH formation in plasma pyrolysis and can pave the way for more advanced mechanistic studies (e.g., soot nucleation mechanisms). Additionally, the findings will be useful for establishing practical strategies for increasing the pyrolysis efficiency and producing high-grade carbon for synthesis of nanomaterials.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2010

Uncertainty quantification given discontinuous climate model response and a limited number of model runs.

Uncertainty quantification in complex climate models is challenged by the sparsity of available climate model predictions due to the high computational cost of model runs. Another feature that prevents classical uncertainty analysis from being readily applicable is bifurcative behavior in climate model response with respect to certain input parameters. A typical example is the Atlantic Meridional Overturning Circulation. The predicted maximum overturning stream function exhibits discontinuity across a curve in the space of two uncertain parameters, namely climate sensitivity and CO2 forcing. We outline a methodology for uncertainty quantification given discontinuous model response and a limited number of model runs. Our approach is two-fold. First we detect the discontinuity with Bayesian inference, thus obtaining a probabilistic representation of the discontinuity curve shape and location for arbitrarily distributed input parameter values. Then, we construct spectral representations of uncertainty, using Polynomial Chaos (PC) expansions on either side of the discontinuity curve, leading to an averaged-PC representation of the forward model that allows efficient uncertainty quantification. The approach is enabled by a Rosenblatt transformation that maps each side of the discontinuity to regular domains where desirable orthogonality properties for the spectral bases hold. We obtain PC modes by either orthogonal projection or Bayesian inference, and argue for a hybrid approach that targets a balance between the accuracy provided by the orthogonal projection and the flexibility provided by the Bayesian inference - where the latter allows obtaining reasonable expansions without extra forward model runs. The model output, and its associated uncertainty at specific design points, are then computed by taking an ensemble average over PC expansions corresponding to possible realizations of the discontinuity curve. The methodology is tested on synthetic examples of discontinuous model data with adjustable sharpness and structure.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Which Combinations of Environmental Conditions and Microphysical Parameter Values Produce a Given Orographic Precipitation Distribution?

This study applies an idealized modeling framework, alongside a Bayesian Markov chain Monte Carlo (MCMC) algorithm, to explore which combinations of upstream environmental conditions and cloud microphysical parameter values can produce a particular precipitation distribution over an idealized two-dimensional, bell-shaped mountain. Simulations focus on orographic precipitation produced when an atmospheric river interacts with topography. MCMC-based analysis reveals that different combinations of parameter values produce a similar precipitation distribution, with the most influential parameters being relative humidity (RH), horizontal wind speed (U), surface potential temperature (θ sfc ), and the snow fall speed coefficient (A s ). RH, U, and A s exhibit interdependence: changes in one or more of these factors can be mitigated by compensating changes in the other(s) to produce similar orographic precipitation rates. The results also indicate that the parameter sensitivities and relationships can vary for spatial subregions and given different environmental conditions. In particular, high θ sfc values are more likely to produce the target precipitation rate and spatial distribution, and thus the ensemble of simulations shows a preference for liquid precipitation at the surface. The results presented here highlight the complexity of orographic precipitation controls, and have implications for flood and water management, observational efforts, and climate change.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗