Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Inverse model calibration”

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

Nuclear Data Adjustment for Nonlinear Applications in the OECD/NEA WPNCS SG14 Benchmark—A Bayesian Inverse UQ-Based Approach for Data Assimilation

The Organisation for Economic Co-operation and Development Working Party on Nuclear Criticality Safety has proposed a benchmark exercise to assess the performance of current nuclear data adjustment techniques applied to nonlinear applications and experiments with low correlation to applications. This work introduces Bayesian inverse uncertainty quantification (IUQ) employing scientific machine learning surrogate models as a method for nuclear data adjustments in this benchmark, and compares IUQ to the more traditional methods of generalized linear least squares (GLLS) and Monte Carlo Bayes (MOCABA). Posterior predictions from IUQ showed agreement with GLLS and MOCABA for linear applications. Here, when comparing GLLS, MOCABA, and IUQ posterior predictions to computed model responses using adjusted parameters, we observe that the GLLS predictions failed to replicate the computed response distributions for nonlinear applications, while MOCABA showed near agreement, and IUQ used the computed model responses directly. We also discuss observations on why experiments with low correlation to applications can be informative to nuclear data adjustments and identify some properties useful in selecting experiments for inclusion in nuclear data adjustment. Performance in this benchmark indicates potential for Bayesian IUQ in nuclear data adjustments.

Bayesian calibration↗

Applicaton of radiative transfer theory to microwave transmission medium calibrations

Precise determinations of the transmission medium loss and noise temperature contribution which are important to the performance characterization of low noise microwave receiving systems and thermal noise standards are discussed. Tropospheric loss is frequently inferred from microwave radiometer noise temperature measurements. Interpretation of these measurements requires an inversion of the radiative transfer integral equation. This is inconvenient even with computer techniques. Solutions of a rapidly convergent power series of the radiative transfer equations are presented. This solution is applicable to a low loss medium with either uniform or nonuniform loss distributions. A four layer atmosphere model is investigated to demonstrate the accuracy of the solution relative to the model. Applications include thermal noise standards and single- and dual-frequency water radiometers.

Stelzried, C. T.↗

Bayesian State-Space Modeling Framework for Understanding and Predicting Golden Eagle Movements Using Telemetry Data

Predicting raptor movements through a wind power plant under given atmospheric and topographical conditions is a crucial first step in the overall goal of quantifying the risk of turbine-related collisions and mortalities. Extracting behavioral traits of golden eagles (Aquila chrysaetos) from telemetry data requires the fusion of noisy and sparse movement data (location, heading, velocity) with a stochastic mathematical representation of the eagles' decision-making processes. In this study, we framed this problem in a Bayesian state-space framework where both observations and decision-making are assumed to be stochastic processes connected through hidden states (mode of flight, intent), and the unknown model parameters are assumed to be random variables that are calibrated using the available telemetry data. This framework allowed for rigorous consideration of underlying uncertainties while allowing for both data and prior biological knowledge to contribute to a probabilistic and predictive agent-based movement model. We implemented and applied the Bayesian framework to understand movement behavior of 23 GPS-tagged golden eagles travelling in the western US for years 2019 and 2020. Our preliminary findings show that the Bayesian state-space framework provides a robust inverse modeling apparatus to decode eagle behavioral characteristics from telemetry data. This study was primarily aimed at verifying and validating the framework with selected golden eagle tracks (both long- and short-ranged), with future research aimed at extending the framework to include multi-mode flight, consideration of atmospheric data and uplift mechanisms, eagle-to-eagle interaction, and eagle-to-turbine interaction.

Bayesian modeling↗

The Dark Energy Survey Supernova Program: an updated measurement of the Hubble constant using the inverse distance ladder

We measure the current expansion rate of the Universe, Hubble’s constant $H_0$, by calibrating the absolute magnitudes of supernovae to distances measured by baryon acoustic oscillations (BAO). This ‘inverse distance ladder’ technique provides an alternative to calibrating supernovae using nearby absolute distance measurements, replacing the calibration with a high-redshift anchor. We use the recent release of 1829 supernovae from the Dark Energy Survey spanning $0.01\lt z\lt 1.13$ anchored to the recent baryon acoustic oscillation measurements from Dark Energy Spectroscopic Instrument (DESI) spanning $0.30 \lt z_{\mathrm{eff}}\lt 2.33$. To trace cosmology to $z=0$, we use the third-, fourth-, and fifth-order cosmographic models, which, by design, are agnostic about the energy content and expansion history of the universe. With the inclusion of the higher redshift DESI-BAO data, the third-order model is a poor fit to both data sets, with the fourth-order model being preferred by the Akaike Information Criterion. Using the fourth-order cosmographic model, we find $H_0=67.19^{+0.66}_{-0.64}\mathrm{~km} \mathrm{~s}^{-1} \mathrm{~Mpc}^{-1}$, in agreement with the value found by Planck without the need to assume Flat-$\Lambda$CDM. However, the best-fitting expansion history differs from that of Planck, providing continued motivation to investigate these tensions.

79 ASTRONOMY AND ASTROPHYSICS↗

Radiative Flux Measurements from ARISE: A Comparison with CERES Top-of-Atmosphere Radiative Fluxes

Uncertainty in top-of-atmosphere (TOA) radiation fluxes observations are larger in the Arctic than in other regions. These uncertainties are due to the low sun angles and the highly reflective, anisotropic, and heterogeneous surface conditions. Quantifying, attributing, and reducing Arctic TOA radiative flux uncertainty enables a better understanding of the rapidly changing Arctic. To advance this goal, we compare the Cloud and Earth’s Radiant Energy System (CERES) TOA radiative fluxes with Arctic Radiation-IceBridge Sea and Ice Experiment (ARISE) campaign measurements collected in September 2014. We compare CERES TOA and aircraft radiative flux measurements using two complementary approaches: grid box average fluxes and instantaneously matched footprints. The grid box mean flux comparison indicates an agreement between CERES and aircraft measurements within 2 uncertainty (calibration and inversion) in the longwave for all five grid boxes and for four-of-five grid boxes in the shortwave; shortwave and longwave mean differences are -7.9 and +2.3 Wm‑2, respectively. The comparison of 36 instantaneously matched footprints with aircraft measurements reveals mean differences of -10.5 and 0.4 Wm‑2 in the shortwave and longwave, respectively. To further explore the persistent negative difference in the shortwave, we further quantify the effects of temporal and spatial sampling differences, angular distribution models, and scene identification to CERES-aircraft differences. Our analysis indicates that sampling differences (including scene evolution) account for an additional 1.8 and 1.7% uncertainty in the shortwave and longwave, respectively and indicates no bias. After accounting for this sampling uncertainty, all CERES-aircraft grid box mean fluxes agree within 2 uncertainty. Scene identification errors due to sea ice concentration data set differences exhibit no bias in the shortwave flux difference and indicate the possibility of substantial differences in the CERES fluxes in specific cases with large spatial heterogeneity. Considering the instantaneously matched footprints, we find that the angular distribution models account may account for up to ‑7.3 Wm-2 of the persistent CERES-aircraft shortwave flux difference due to systematic differences in the anisotropy for sea ice partly cloudy scenes. Additional analysis using a special scan model with the CERES FM2 instrument suggests a significant view zenith angle dependence of the CERES fluxes for sea ice partly cloudy scenes where shortwave fluxes systematically decrease with increasing view zenith angle; no dependence is found for other scene types. We conclude that (1) spatial heterogeneity and scene temporal evolution substantially limit our ability to use aircraft measurements to place strong constraints on CERES TOA fluxes and (2) that the representation of anisotropy in sea ice partly cloudy scenes is a significant factor contributing to the persistent negative CERES-aircraft shortwave flux difference in this comparison and require additional data to analysis fully quantify the potential bias.

Patrick C Taylor↗

A comparison between CERES TOA Radiative Fluxes and Airborne Radiative Flux Measurements from ARISE

Uncertainty in top-of-atmosphere (TOA) radiation fluxes observations are larger in the Arctic than in other regions. These uncertainties are due to the low sun angles and the highly reflective, anisotropic, and heterogeneous surface conditions. Quantifying, attributing, and reducing Arctic TOA radiative flux uncertainty enables a better understanding of the rapidly changing Arctic. To advance this goal, we compare the Cloud and Earth’s Radiant Energy System (CERES) TOA radiative fluxes with Arctic Radiation-IceBridge Sea and Ice Experiment (ARISE) campaign measurements collected in September 2014. We compare CERES TOA and aircraft radiative flux measurements using two complementary approaches: grid box average fluxes and instantaneously matched footprints. The grid box mean flux comparison indicates an agreement between CERES and aircraft measurements within 2σ uncertainty (calibration and inversion) in the longwave for all five grid boxes and for four-of-five grid boxes in the shortwave; shortwave and longwave mean differences are -7.9 and +2.3 Wm 2, respectively. The comparison of 36 40 instantaneously matched footprints with aircraft measurements reveals mean differences of -120.25 and -1.00.4 Wm 2 in the shortwave and longwave, respectively. To further explore the persistent negative difference in the shortwave, Wwe further quantify the effects of temporal and spatial sampling differences, angular distribution models, and scene identification to CERES-aircraft differences. Our analysis indicates that sampling differences (including scene evolution) account for an additional 1.8 and 1.7% uncertainty in the shortwave and longwave, respectively and , but indicates no bias. After accounting for this sampling uncertainty, all CERES-aircraft grid box mean fluxes agree within 2σ uncertainty. Scene identification errors due to sea ice concentration data set differences exhibit no bias in the shortwave flux difference and indicate the possibility of substantial differences in the CERES fluxes in specific cases with large spatial heterogeneity. Considering the instantaneously matched footprints, we find that the angular distribution models account may account for up to 7.3 Wm-2 of the persistent CERES-aircraft shortwave flux difference due to systematic differences in the anisotropy for sea ice partly cloudy scenes. Additional analysis using a special programmable scan model with the CERES FM2 instrument suggests a significant view zenith angle dependence of the CERES fluxes for sea ice partly cloudy scenes where shortwave fluxes systematically decrease with increasing view zenith angle; no dependence is found for other scene types. We conclude that (1) spatial heterogeneity and scene temporal evolution substantially limit our ability to use aircraft measurements to place strong constraints on CERES TOA fluxes and (2) that the representation of anisotropy in sea ice partly cloudy scenes is likely a significant factor contributing to the persistent negative CERES-aircraft shortwave flux difference in this comparison and require additional data to analysis fully quantify the potential bias.

Patrick C. Taylor↗

The influence of tillage and fertilizer on the flux and source of nitrous oxide with reference to atmospheric variation using laser spectroscopy

Nitrous oxide (N 2 O) is the third most important long-lived greenhouse gas and agriculture is the largest source of N 2 O emissions. Curbing N 2 O emissions requires understanding influences on the flux and sources of N 2 O. We measured flux and evaluated microbial sources of N 2 O using site preference ( S P ; the intramolecular distribution of 15 N in N 2 O) in flux chambers from a grassland tilling and agricultural fertilization experiments and atmosphere. We identified values greater than that of the average atmosphere to reflect nitrification and/or fungal denitrification and those lower than atmosphere as increased denitrification. Our spectroscopic approach was based on an extensive calibration with 18 standards that yielded S P accuracy and reproducibility of 0.7 ‰ and 1.0 ‰, respectively, without preconcentration. Chamber samples from the tilling experiment taken ~ monthly over a year showed a wide range in N 2 O flux (0–1.9 g N 2 O-N ha -1 d -1 ) and S P (- 1.8 to 25.1 ‰). Flux and S P were not influenced by tilling but responded to sampling date. Large fluxes occurred in October and May in no-till when soils were warm and moist and during a spring thaw, an event likely representing release of N 2 O accumulated under snow cover. These high fluxes could not be ascribed to a single microbial process as S P differed among chambers. However, the year-long S P and flux data for no-till showed a slight direct relationship suggesting that nitrification increased with flux. The comparative data in till showed an inverse relationship indicating that high flux events are driven by denitrification. Corn ( Zea mays ) showed high fluxes and S P values indicative of nitrification ~ 4 wk after fertilization with subsequent declines in S P indicating denitrification. Although there was no effect of fertilizer treatment on flux or S P in switchgrass ( Panicum virgatum) , high fluxes occurred ~1 month after fertilization. In both treatments, S P was indicative of denitrification in many instances, but evidence of nitrification/fungal denitrification also prevailed. At 2 m atmospheric N 2 O S P had a range of 31.1 ‰ and 14.6 ‰ in the grassland tilling and agricultural fertilization experiments, respectively. These data suggest the influence of soil microbial processes on atmospheric N 2 O and argue against the use of the global average atmospheric S P in isotopic modeling approaches.

Environmental Sciences & Ecology↗

Quantification of LEU Holdup using gamma ray imaging and inverse transport solver

Holdup is the residual amount of special nuclear material (SNM) remaining in a processing facility after the bulk materials have been cleaned out. In commercial uranium processing facilities, quantification of holdup is a major challenge because of the highly variable shapes and sizes of the deposits. Any method that attempts to generalize and calibrate deposit shapes in order to quantify holdup will be prone to high uncertainties. Uncertainties on the order of ±50% are typical in holdup results. In international safeguards applications, a ±50% uncertainty can result in a large amount of material unaccounted for (MUF) thereby increasing the difficulty of detecting material diversion and facility misuse. An imaging-based methodology has been developed with the objective of significantly reducing this uncertainty by using the true deposit shape, instead of relying on oversimplified geometric assumptions. The project is a collaboration between ORNL, Y-12, and the University of Tennessee, Knoxville, TN. Uranium sources of known masses were measured using the Germanium Gamma-ray Imager (GeGI), a high-resolution imaging spectrometer, creating a pixelated map for each spectral bin. Two different gamma imaging methods are employed in this work: coded aperture imaging and Compton imaging. A validated MonteCarlo model of the detector has been developed using the GEANT4 code for determining the intrinsic response of the detector, its enclosure, and the coded aperture mask. An inverse transport solver based on the Markov Chain Monte-Carlo approach known as Differential Evolution Adaptive Metropolis (DREAM) is employed to use the measurement data from the image pixels (coded aperture or Compton) to solve for the mass of 235 U in the deposit. A reliable method based on the DREAM solver has been developed to flag the infinite thickness condition of a uranium deposit. The project team is working towards improving the image reconstruction for Compton imaging so that a better localization of the source can be achieved. Besides treating the coded aperture and Compton imaging methods independently, the project is also evaluating a combined method that uses the Compton scatter data from a coded aperture measurement. GEANT4 simulations are being performed to evaluate the combined approach. The impact on the DREAM optimization as the source thickness progressively approaches infinite thickness is being evaluated. A number of uranium sources available at ORNL have been measured, and the DREAM results have been tested and validated for the coded aperture imaging. A similar effort will be carried out to validate the Compton based method once the development of algorithms for better localization are complete. The imaging based quantification is very amenable to unattended monitoring of holdup accumulation at key measurement points. A proof of concept measurement has been completed to demonstrate this capability The current work used the high energy resolution imager GeGI. However, the approach and methodologies are applicable to other imagers such as the cadmium zin telluride (CZT) based imager manufactured by H3D, Inc.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Viscosity estimates for the crust and upper mantle from patterns of lacustrine shoreline deformation in the Eastern Great Basin

The deformed shorelines of Lake Bonneville constitute a classic source of information on lithospheric elastic thickness and upper mantle viscosity. We describe and apply a new model to a recently augmented data set. New data better constrain both the complex spatio-temporal pattern of the lake load and the crustal deformation response to that load. The history of lake level fluctuations has been significantly refined and somewhat modified. This is due to both more radiocarbon dates from within the Bonneville basin and to an improved calibration of the radiocarbon timescale itself. The data which constrain the crustal deformation pattern consist of ages and shoreline elevations from several hundred points which sample three major levels of Lake Bonneville and corresponding elevations from the high stands of three smaller lakes situated to the west of Lake Bonneville. The geometry of the Earth model incorporates an arbitrary number of layers overlying a half-space, and the rheology of each level can accommodate an arbitrary number of Maxwell viscoelastic elements in parallel. The inverse modeling comprises three complementary approaches: for the simplest configurations, we performed a direct search of the parameter space and delineated the irregular boundary of the subspace of acceptable models. For more complex configurations, we constrained the elastic parameters to their seismically determined values and then solved for viscosity versus depth profiles by either expressing the log(viscosity) versus log(depth) profile as a series of specially constructed orhtogonal polynomials, or by allowing each of 8-10 layers (plus the half-space) to have an independently determined viscosity. We found that the data do not strongly support (nor can they conclusively exclude) a more complex rheology than simple Maxwell viscoelasticity. The orthogonal polynomial solution exhibits an essentially monotonic decrease in viscosity with depth.

Bills, Bruce G.↗

A Bayesian Approach to Regional Decadal Predictability: Sparse Parameter Estimation in High-Dimensional Linear Inverse Models of High-Latitude Sea Surface Temperature Variability

Stochastic reduced models are an important tool in climate systems whose many spatial and temporal scales cannot be fully discretized or underlying physics may not be fully accounted for. One form of reduced model, the linear inverse model (LIM), has been widely used for regional climate predictability studies—typically focusing more on tropical or midlatitude studies. However, most LIM fitting techniques rely on point estimation techniques deriving from fluctuation–dissipation theory. In this methodological study we explore the use of Bayesian inference techniques for LIM parameter estimation of sea surface temperature (SST), to quantify the skillful decadal predictability of Bayesian LIM models at high latitudes. We show that Bayesian methods, when compared to traditional point estimation methods for LIM-type models, provide better calibrated probabilistic skill, while simultaneously providing better point estimates due to the regularization effect of the prior distribution in high-dimensional problems. We compare the effect of several priors, as well as maximum likelihood estimates, on 1) estimating parameter values on a perfect model experiment and 2) producing calibrated 1-yr SST anomaly forecast distributions using a preindustrial control run of the Community Earth System Model (CESM). Finally, we employ a host of probabilistic skill metrics to determine the extent to which an LIM can forecast SST anomalies at high latitudes. We find that the choice of prior distribution has an appreciable impact on estimation outcomes, and priors that emphasize physically relevant properties enhance the model’s ability to capture variability of SST anomalies.

54 ENVIRONMENTAL SCIENCES↗

Retrieval of Surface Directional Reflectance Properties Using Ground Level Multiangle Measurements

Knowledge of the directional reflectance properties of natural surfaces such as soils and vegetation canopies is essential for classification studies and canopy model inversion. Atmospheric correction schemes, using various levels of approximation, are described to retrieve surface bidirectional reflectance factors (BRFs) and directional-hemispherical reftectances (albedos)from multiangle radiance measurements taken at ground level. The retrieval schemes are tested on simulated data incorporating realistic surface BRFs and atmospheric models containing aerosols. Sensitivity of the atmospherically corrected BRFs and associated directional-hemispherical reflectances to various aerosol properties and the sun-view geometry is illustrated. A measurement strategy for obtaining highly accurate surface reflectance properties also is examined in the context of instrument radiametric calibration, knowledge of the atmospheric properties, and sun-view angular coverage.

Martonchik, John V.↗

Improving the Quasi‐Biennial Oscillation via a Surrogate‐Accelerated Multi‐Objective Optimization

Accurate simulation of the quasi-biennial oscillation (QBO) is challenging due to uncertainties in representing convectively generated gravity waves. We develop an end-to-end uncertainty quantification workflow that calibrates these gravity wave processes in E3SM for a realistic QBO. Central to our approach is a domain knowledge-informed, compressed representation of high-dimensional spatio-temporal wind fields. By employing a parsimonious statistical model that learns the fundamental frequency from complex observations, we extract interpretable and physically meaningful quantities capturing key attributes. Building on this, we train a probabilistic surrogate model that approximates the fundamental characteristics of the QBO as functions of critical physics parameters governing gravity wave generation. Leveraging the Karhunen–Loève decomposition, our surrogate efficiently represents these characteristics as a set of orthogonal features, capturing cross-correlations among multiple physics quantities evaluated at different pressure levels and enabling rapid surrogate-based inference at a fraction of the computational cost of full-scale simulations. Finally, we analyze the inverse problem using a multi-objective approach. Our study reveals a tension between amplitude and period that constrains the QBO representation, precluding a single optimal solution. To navigate this, we quantify the bi-criteria trade-off and generate a set of Pareto optimal parameter values that balance the conflicting objectives. This integrated workflow improves the fidelity of QBO simulations and offers a versatile template for uncertainty quantification in complex geophysical models.

54 ENVIRONMENTAL SCIENCES↗

Tests of two convection theories for red giant and red supergiant envelopes

Two theories of stellar envelope convection are considered here in the context of red giants and red supergiants of intermediate to high mass: Boehm-Vitense's standard mixing-length theory (MLT) and Canuto & Mazzitelli's new theory incorporating the full spectrum of turbulence (FST). Both theories assume incompressible convection. Two formulations of the convective mixing length are also evaluated: l proportional to the local pressure scale height (H(sub P)) and l proportional to the distance from the upper boundary of the convection zone (z). Applications to test both theories are made by calculating stellar evolutionary sequences into the red zone (z). Applications to test both theories are made by calculating stellar evolutionary sequences into the red phase of core helium burning. Since the theoretically predicted effective temperatures for cool stars are known to be sensitive to the assigned value of the mixing length, this quantity has been individually calibrated for each evolutionary sequence. The calibration is done in a composite Hertzsprung-Russell diagram for the red giant and red supergiant members of well-observed Galactic open clusters. The MLT model requires the constant of proportionality for the convective mixing length to vary by a small but statistically significant amount with stellar mass, whereas the FST model succeeds in all cases with the mixing lenghth simply set equal to z. The structure of the deep stellar interior, however, remains very nearly unaffected by the choices of convection theory and mixing lenghth. Inside the convective envelope itself, a density inversion always occurs, but is somewhat smaller for the convectively more efficient MLT model. On physical grounds the FST model is preferable, and seems to alleviate the problem of finding the proper mixing length.

Stothers, Richard B.↗

A Method to Solve Interior and Exterior Camera Calibration Parameters for Image Resection

An iterative method is presented to solve the internal and external camera calibration parameters, given model target points and their images from one or more camera locations. The direct linear transform formulation was used to obtain a guess for the iterative method, and herein lies one of the strengths of the present method. In all test cases, the method converged to the correct solution. In general, an overdetermined system of nonlinear equations is solved in the least-squares sense. The iterative method presented is based on Newton-Raphson for solving systems of nonlinear algebraic equations. The Jacobian is analytically derived and the pseudo-inverse of the Jacobian is obtained by singular value decomposition.

Samtaney, Ravi↗

Evaluating Soil Moisture Retrievals from ESA's SMOS and NASA's SMAP Brightness Temperature Datasets

Two satellites are currently monitoring surface soil moisture (SM) using L-band observations: SMOS (Soil Moisture and Ocean Salinity), a joint ESA (European Space Agency), CNES (Centre national d'tudes spatiales), and CDTI (the Spanish government agency with responsibility for space) satellite launched on November 2, 2009 and SMAP (Soil Moisture Active Passive), a National Aeronautics and Space Administration (NASA) satellite successfully launched in January 2015. In this study, we used a multilinear regression approach to retrieve SM from SMAP data to create a global dataset of SM, which is consistent with SM data retrieved from SMOS. This was achieved by calibrating coefficients of the regression model using the CATDS (Centre Aval de Traitement des Donnes) SMOS Level 3 SM and the horizontally and vertically polarized brightness temperatures (TB) at 40 deg incidence angle, over the 2013 - 2014 period. Next, this model was applied to SMAP L3 TB data from Apr 2015 to Jul 2016. The retrieved SM from SMAP (referred to here as SMAP_Reg) was compared to: (i) the operational SMAP L3 SM (SMAP_SCA), retrieved using the baseline Single Channel retrieval Algorithm (SCA); and (ii) the operational SMOSL3 SM, derived from the multiangular inversion of the L-MEB model (L-MEB algorithm) (SMOSL3). This inter-comparison was made against in situ soil moisture measurements from more than 400 sites spread over the globe, which are used here as a reference soil moisture dataset. The in situ observations were obtained from the International Soil Moisture Network (ISMN; https:ismn.geo.tuwien.ac.at) in North of America (PBO_H2O, SCAN, SNOTEL, iRON, and USCRN), in Australia (Oznet), Africa (DAHRA), and in Europe (REMEDHUS, SMOSMANIA, FMI, and RSMN). The agreement was analyzed in terms of four classical statistical criteria: Root Mean Squared Error (RMSE),Bias, Unbiased RMSE (UnbRMSE), and correlation coefficient (R). Results of the comparison of these various products with in situ observations show that the performance of both SMAP products i.e. SMAP_SCA and SMAP_Reg is 48 similar and marginally better to that of the SMOSL3 product particularly over the PBO_H2O, SCAN, and USCRN sites. However, SMOSL3 SM was closer to the in situ observations over the DAHRA and Oznet sites. We found that the correlation between all three datasets and in situ measurements is best (R 0.80) over the Oznet sites and worst (R 0.58) over the SNOTEL sites for SMAP_SCA and over the DAHRA and SMOSMANIA sites (R 0.51 and R 0.45 for SMAP_Reg and SMOSL3, respectively). The Bias values showed that all products are generally dry, except over RSMN, DAHRA, and Oznet (and FMI for SMAP_SCA). Finally, our analysis provided interesting insights that can be useful to improve the consistency between SMAP and SMOS datasets.

SMOS↗

Gamma ray Imaging Based Methods to Improve the Accuracy of Uranium Holdup Quantification

Deposits in commercial uranium processing facilities of interest to nuclear safeguards are predominantly low enriched uranium (LEU). The shapes and sizes of holdup deposits can vary significantly, thus making it difficult to model the deposits accurately for calibration purposes. This presents a challenge to currently employed methods such as the Generalized Geometry Holdup (GGH) that rely on quantifying 235U mass by simplifying deposit shapes as a point, a line, or an area. In this work, gamma ray imaging using high energy resolution Germanium Gamma Imagers (GeGI™) is employed to determine the distribution of uranium inside the source containment. Two types of imaging methods are employed; coded-aperture imaging and Compton imaging. Gamma-ray emissions seen in different regions of the image are quantified using an inverse gamma-ray transport solver being developed at ORNL. Using measured data, the inverse method will be used to solve for unknown source parameters such as source matrix thickness, density, and attenuation due to container wall and shielding. The intrinsic efficiency of the detector is determined based on a library of built-in response functions and key dimensions of the detector (e.g., thickness and radius of the detector crystal). Measurements are performed in both imaging modalities using uranium sources of well-known masses and enrichments configured inside mocked-up holdup fixtures. Results from the two imaging modalities are intercompared. The paper presents preliminary results of the imaging measurements and the progress made in the computational methods to quantify uranium mass.

98 NUCLEAR DISARMAMENT, SAFEGUARDS, AND PHYSICAL P↗

A hierarchical Bayesian model for background variation in radiation source localization

In this paper, we apply a new model to account for varying background radiation in radiological source localization. We present a hierarchical Bayesian model that simultaneously infers background and source location parameters without requiring separate estimation of the background radiation at each detector location. We employ a simplified photon transport model to reduce the computational expense of Bayesian model calibration. Furthermore, we demonstrate the model accuracy by localizing a cesium-137 source in a simulated city block, and we analyze experimental field measurements with varying background. In both cases, the model provides sufficient fidelity that we can locate the source while simultaneously estimating background radiation.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Reduced scale stripline platform to extend accessible pressures on the Z machine: Shockless compression of platinum to 650 GPa

Reaching astrophysically relevant high energy density (HED) material states in the laboratory is an ongoing effort at multiple experimental facilities. We have developed a new dynamic compression platform for the Z Pulsed Power Facility that allows for sample sizes 100s of [Formula: see text]m in thickness that accommodate multiple grains in order to fully capture bulk properties, such as material strength. A pair of experiments compressed platinum (Pt) to HED conditions and conventional inverse Lagrangian analysis as well as a recent Bayesian calibration technique were used to determine the principal isentrope to 650 GPa with density uncertainties of <2%. These low uncertainties are calculated for single sample experiments, presenting the possibility of even smaller experimental uncertainties with multiple samples the platform allows. Our new platform extends the accessible Pt ramp pressures achievable on the Z machine to over 80% of the pressure recently achieved using the National Ignition Facility planar Hohlraum platform. This new capability, the next generation evolution of the stripline platform, was made possible by advancements in both our understanding of the Z pulsed power driver and our overall magnetohydrodynamic modeling capabilities.

Porwitzky, Andrew↗