Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Parameterization”

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 649 records · Page 36

PBX9502 PAGOSA/SURF Calibration for 100 µm Grid Size

There are two methods to calibrate the SURF reactive burn explosives model.The first is to calibrate using Pop plot data. Because there are only two main parameters in SURF, called A and B, there is just enough flexibility in SURF to determine the Pop plot line, i.e., equivalent information to a slope and a single point on the line. However, due to the linearity of the Pop plot being based on logarithmic scales, a slight deviation off the Pop plot (though the fit to the Pop plot might be very good overall) can make considerable difference when comparing with the velocity profile data. The second calibration method thusly focuses strongly on matching entire velocity profile data. This is to say that individual velocity profiles extracted from embedded gages (in the case of PBX 9502, we refer to R. Gustavsen’s experimental data) may fit poorly, while the overall fit to the Pop plot may seem reasonable. In addition, there are other factors that could produce slightly different calibrations for SURF, while leading to an overall acceptable fit to data. One of those may be the physics code used and here we use Pagosa. Another is mesh size dependence. A third one stems from fitting a reactive burn model being done in conjunction with an equation of state (EOS) of the reactants and one for the products – different EOS can result in slightly different parameterizations. Finally, there are several secondary material parameters in SURF, and the overall fit will depend on the entire parameterization. The goal of this report is to clarify some of our previously made statements about SURF calibration, as well as to redo a systematic calibration using 100 µm grid resolution.

42 ENGINEERING↗

Macro-physical Properties of Shallow Cumulus from Integrated ARM Observations (Final Report)

Fair-weather shallow cumuli (ShCu) play an important role in many climate-related processes. Irregular geometry of ShCu and their strong temporal and spatial variability make it challenging to observe ShCu holistically and to represent them correctly in climate models. To improve ShCu parameterizations, information on both vertically and horizontally resolved cloud properties is required. Commonly, the vertically resolved cloud properties are provided by zenith pointing lidar-radar observations with a very narrow field of view (FOV). Thus, these “pencil-beam” properties may not be representative of a larger surrounding area. Limited number of areal-averaged cloud properties, such as fractional sky cover (FSC), are offered typically by wide-FOV observations. The main goal of our project was to integrate advantages of the narrow-FOV (vertical structure of clouds) and wide-FOV (spatial arrangement of clouds) observations for an improved characterization of single-layer ShCu observed at the Atmospheric Radiation Measurement (ARM) Southern Great Plains (SGP) site for an 18-yr period (2000-2017). There are four major accomplishments of our project, First, an updated operational cloud classification for days with ShCu has been suggested and evaluated through a detailed comparison with the manually curated records. Our classification extends successfully the latest ARM cloud type Value Added Product (VAP) based on the Active Remote Sensing of Clouds (ARSCL) cloud product by incorporating both cloud fraction (CF) provided by narrow-FOV ceilometer data and FSC from wide-FOV images offered by a Total Sky Imager (TSI). Moreover, our classification allows one to identify impact of instrumentation changes at the SGP site, namely the transition to KAZRARSCL with the updated cloud radar, on the identification of periods with single-layer ShCu. Second, a new approach that resolves cloud area distributions for a given region (up to 4x4 km 2 ) has been suggested and cloud equivalent diameters (CEDs) have been estimated for the first time. These estimations have been performed over a wide range of cloud sizes (about 0.01–3.5 km) with high temporal resolution (30s) using wide-FOV TSI images and cloud base height (CBH) provided by complementary narrow-FOV lidar measurements. Our simple and computationally inexpensive approach offers a previously unavailable dataset for process studies in the convective boundary layer and evaluation of ShCu parameterizations in cloud-resolving models. Third, a long-term integrated record of ShCu macrophysical properties has been developed. The developed record represents the longest available compilation of events with ShCu and includes (i) a novel visualization of the spatial variability in cloud cover both along- and across-wind directions, (ii) updated estimates of narrow-FOV CF and wide-FOV FSC, (iii) updated narrow-FOV CBH, and (iv) complementary data, such as wind speed and direction from the 915-MHz Radar Wind Profiler (RWP) data. The developed record has been used successfully to assess conventional observational estimates of cloud cover and their sensitivity to the following two factors: (i) instrument-dependent cloud detection and data merging criteria and (ii) FOV configuration. Fourth, co-variability of the ShCu macrophysical properties and environmental parameters has been analyzed for a 3-yr period (2016-2018). Our initial analysis includes diurnal changes of FSCs obtained for clouds with small, moderate and large CEDs and several environmental parameters, such as lifted condensation level (LCL) and mixed layer height (zi). Preliminary results of our analysis suggest that the horizontal extent of ShCu is controlled substantially by the sign and magnitude of difference between these two parameters (zi-LCL): the CED tends to grow with increase of this difference (zi exceeds LCL). We have initiated relationships between the ShCu and key atmospheric parameters that control both the development and evolution of ShCu using our new data product, which combines effectively the advantages of narrow-FOV data offered by zenith pointing cloud radars and lidars and wide-FOV TSI images. While the latest instrumentation at the ARM sites may address these challenging relationships in the future, we believe that the historical ARM data at the SGP site has not yet been fully utilized. Overall, our data product can be used by researchers working on a wide range of climate-related projects. These projects may include (i) a comprehensive evaluation of outputs from the Large-Eddy Simulation (LES) and single-column models for their future improvement, (ii) the representativeness of “short-period” results obtained from the previous model and observational studies and (iii) the planning of future field campaigns with focus on improved understanding of the diurnal cycle of cumulus convection.

54 ENVIRONMENTAL SCIENCES↗

Implementation of a Quasi-3D Multiscale Modeling Framework to ACME

Representation of unresolved cloud processes has been one of the most challenging problems in numerical modeling of the atmosphere. During the last decades, considerable efforts have been spent to better represent the collective effects of those processes in terms of resolved scales (see Arakawa 2004 for a review). Still, the progress on this problem has been unacceptably slow as pointed out by Randall et al. (2003). The large part of the difficulty is inherent in the artificial separation of atmospheric processes into the explicitly simulated resolved scales and the highly parameterized subgrid scales. Arakawa (2004) and Arakawa et al. (2011) suggested two plausible approaches to eliminate the problem of artificial separation ofatmospheric processes: “parameterize everything” and “resolve everything”.

58 GEOSCIENCES↗

Use of Remote Sensing and In-situ Observations to Develop and Evaluate Improved Representations of Convection and Clouds for the Accelerated Climate Model for Energy. Final report

Representation of unresolved cloud processes has been one of the most challenging problems in numerical modeling of the atmosphere. During the last decades, considerable efforts have been spent to better represent the collective effects of those processes in terms of resolved scales (see Arakawa 2004 for a review). Still, the progress on this problem has been unacceptably slow as pointed out by Randall et al. (2003). The large part of the difficulty is inherent in the artificial separation of atmospheric processes into the explicitly simulated resolved scales and the highly parameterized subgrid scales. Arakawa (2004) and Arakawa et al. (2011) suggested two plausible approaches to eliminate the problem of artificial separation of atmospheric processes: “parameterize everything” and “resolve everything”.

58 GEOSCIENCES↗

Ice Nucleating Particles, Aerosols and Clouds over the Higher Latitude Southern Ocean

Improved understanding of aerosol-cloud-precipitation interactions is a key goal of the DOE-ASR program. Ice Nucleating Particles (INPs) are rare aerosol particles that trigger ice formation in clouds. By doing so they may initiate precipitation in clouds, and this impacts cloud lifetime. Nowhere may the role of INPs be manifested so dramatically as in the Southern Ocean (SO). This vast region was the focus of the DOE-ARM MARCUS (Measurement of Aerosols, Radiation and CloUds over the Southern Oceans) and MICRE (Macquarie Island Cloud and Radiation Experiment) campaigns. In the SO, a deficit of INPs has been a leading hypothesis for why liquid clouds persist, which underlays the bias of global climate models in predicting excess shortwave radiation reaching the ocean surface between 55°S and Antarctica. In this study, we used the 6-month acquisition of INP data on four MARCUS cruises from Hobart, Tasmania, to Antarctica, and the single-location full annual INP cycle observed in MICRE, to gain a holistic, quantitative picture of the INP number concentrations as a function of temperature and their variability over SO latitudes from 43°S to Antarctica. Augmenting initial processing of filter collections of aerosol particles during MARCUS and MICRE, additional freezing studies following thermal treatment (to removes INPs associated with microbes) and peroxide treatment (to remove all organics), as well as ionic chemistry and total organic carbon analysis, were applied to samples to improve time and space resolution, and differentiate marine from terrestrial contributions. We also performed Next Generation Sequencing to determine that bacterial composition can tag periods of enhanced and degraded marine organic INPs and the relative lack of occurrence of land-sourced bio-particles. Through these and other analyses, INP data were categorized into representative source types (e.g., marine versus terrestrially-influenced). The products of our analyses are being used to parameterize ice formation for use in numerical modeling studies to determine if, and when, ocean-derived INPs control the microphysical composition and radiative balance of SO clouds in a manner not presently captured by global models. The specific objectives and approaches proposed were largely accomplished, including: 1) Additional and value-added chemical and biological analyses of archived samples of aerosols were completed to categorize and quantify INP types and concentrations over the SO. 2) Analyses were completed to place INP data in meteorological and aerosol context. 3) Sea spray aerosol INPs were demonstrated to dominate the SO region in all but episodic events. 4) New parameterizations for INP sources relevant to the SO were constructed. 5) Collaborative modeling studies demonstrating the crucial role of INPs in determining cloud radiative and precipitation properties were conducted with MARCUS and MICRE partners. Collaborative publication submissions in this regard ensued within the two year study. This work will advance the science of interactions of aerosols, clouds and precipitation, to improve their representation in regional and global climate models. Results have informed representation of primary ice crystal nucleation, give inference to the role of secondary processes, and improve investigations of aerosol influences on clouds via ice nucleation over SO high latitudes and similar vast ocean regions. Application will improve representation of such interactions for clouds in both regional and global climate models. The data base and analyses methods applied will continue to serve research studies in the future.

58 GEOSCIENCES↗

Unitary Bethe-Salpeter Methods in Two- and Three-Body Systems

The study of resonances, unstable particles formed in particle scattering, is motivated by the questions about the strong interaction and the composition of hadrons. In order to advance this study in the intermediate energy region, a variety of techniques are used. In this thesis, we present work done on two different analyses that use the principle of unitarity and the Bethe-Salpeter equation to study the resonances Lambda(1405), Sigma(1385), and a1(1260). Firstly, we perform a simultaneous analysis of Sand P-waves of the Strangeness = -1 two-body meson-baryon scattering amplitude using all low-energy data. For the first time, differential cross section data are included for chiral unitary coupled-channel models. From this model, S- and P-wave amplitudes are extracted and we observe both well-known I(JP)=0(1/2^-) S-wave states as well as a new I(JP)=1(1/2^+) state absent in quark models and lattice QCD results. Multiple statistical and phenomenological tests suggest that, while the data clearly require an I =1 P-wave resonance, the new state just accounts for the absence of the decuplet Sigma(1385)3/2+ in the model which is subsequently included in the parameterization. Secondly, we formulate the final state interaction of the a1(1260) resonance decay in a manifestly three-body-unitary parameterization and fit it to the a1(1260) lineshape measured by the ALEPH experiment. Dalitz plots calculated from this fit are presented. The work demonstrates the feasibility to numerically solve a previously derived amplitude and its generalization to isobars with spin and coupled channels. The model is a good test case because it can also be applied to other meson decays including exotic states and modified for the finite-volume problem as it arises in lattice QCD due to its manifest unitarity. The understanding of the resonances being analyzed in both the two-body work and the three-body work gives a better understanding of intermediate energy QCD.

Sadasivan, Daniel↗

Iodine Waste Form Corrosion: PNNL Summary Testing Report

A test method to correlate the laboratory corrosion behavior of iodine-containing waste forms (IWFs) to a conceptual long-term corrosion models is required. Several candidate methods have been evaluated previously, and each method presents unique challenges and limitations. Dynamic leach tests have emerged as the preferred candidate test method to evaluate IWFs and parameterize an eventual long-term corrosion model. This report was prepared in support of this effort to define the test suite to evaluate the corrosion of IWFs, with the goal of defining and parameterizing a conceptual long-term corrosion model. A review was presented highlighting the high degree of variability between the IWF corrosion data that is available in the literature and the need for a standardized test suite. The targeted test suite is composed of a single interval screening test, a dynamic leach test, and post-test solid characterization of the monolithic IWF samples. The data presented serves two purposes: to expand the range of IWFs evaluated with the test suite to assess method variability and to further understand the corrosion of IWFs to best inform conceptual long-term corrosion model development. A series of IWF types were included in the testing, including hot isostatically pressed Ag-mordenite, silver functionalized silica-based aerogels, iodosodalites (hot isostatic pressing, sintered, glass bonded ceramic), iodoapatite, tellurite glass and glass composite materials.

42 ENGINEERING↗

Facilitating better and faster simulations of aerosol-cloud interactions in Earth system models

Focal Area(s): 1. Predictive modeling through the use of AI techniques and AI-derived model components; the use of AI and other tools to design a prediction system comprising a hierarchy of models. 2. Insight gleaned from complex data (both observed and simulated) using AI, big data analytics, and other advanced methods, including explainable AI and physics- or knowledge-guided AI. Science Challenge: One major challenge that Earth system models (ESMs) face in providing credible prediction of the Earth system and its water cycle characteristics (e.g., mean state, variability, and extreme events) is to accurately simulate aerosol-cloud interactions (ACI). The physical, chemical, and dynamical processes affecting ACI are extremely complex and they range from nanoscale to planetary scale. In each model development cycle, scientists spend significant efforts investigating model deficiencies and uncertainties associated with aerosols (e.g., emissions, chemical processes, aerosol microphysics, and transport) and clouds (e.g., macrophysics, microphysics, turbulence, and large-scale circulation) in order to develop improved treatments. However, despite decades of active research, ACI is still a major source of uncertainty in climate projections, even though great progress has been made. Specific scientific challenges include: (i) Parameterizations are developed based on limited data; (ii) The complexity of a parameterization required for accurate predictions is not understood; (iii) Incomplete and unknown physics leads to errors in the fully coupled Earth system; and (iv) Complex physics is computationally too expensive to employ in ESMs.

54 ENVIRONMENTAL SCIENCES↗

Structurally flexible cloud microphysics, observationally constrained at all scales via ML-accelerated Bayesian inference

We discuss the challenge of developing observationally informed parameterizations of microphysics for use at a hierarchy of modeling scales. Our proposed approach is applicable to any domain that suffers from a two-fold parameterization problem, where physical processes are not resolved at the model scale (the first problem), and those processes are uncertain at any scale (the second problem). For such problems, a physical approach facilitates modeling across scales, as well as systematic observational inference accelerated by machine learning (ML) surrogate models, and quantification of physical uncertainties.

54 ENVIRONMENTAL SCIENCES↗

Wind Energy Forecasting with the Weather Research and Forecasting Model

This was a collaborative effort between Lawrence Livermore National Security, LLC as manager and operator of Lawrence Livermore National Laboratory (LLNL) and Siemens Energy, Inc. (Siemens) to develop a wind resource forecasting tool. LLNL was to develop an independent high-resolution mesoscale modeling capability forecasting tool that could be implemented in conjunction with existing wind farm control and monitoring software to provide forecasting of wind resources using local observations of winds and temperature. Research with LLNL’s state-of-the-art large-eddy simulation meteorological prediction model, based on the community WRF model and innovative turbulence parameterizations, would improve that model’s applicability to large wind farms offshore and in complex terrain. The modeling capability would include uncertainty quantification. Finally, the application of the modeling tool and existing global climate change predictions would enable the delineation of the likely effects of climate change on wind resources. Siemens was to provide high time resolution hub-height wind speed and other meteorological data streams, including temperature profiles from wind farms, for LLNL to incorporate into the modeling system, to validate and tune this forecasting model for their locations of interest. These data streams would also be used for longer-term studies of correlations of wind resources to climate oscillations to indicate how long-term climate change trends may affect the available wind resource. Siemens would also provide information and observations of turbine wakes for incorporation into the modeling tool. By implementing state-of-the-art turbulence parameterizations into a simulation model and/or ensembles of simulation models, and by integrating real-time hub height wind speed and other meteorological datastreams from wind farms into that model or ensemble of models, LLNL would develop a forecasting tool that could be implemented by Siemens as an add-on to existing wind farm control and monitoring software to provide owners with useful resource forecasting. The desired outcome was that the accuracy level of the output would be sufficient to substantiate power output commitments. The final deliverable for this work would consist of a document outlining the algorithms and software tools that could be integrated into Siemens Wind Park Supervisor.

17 WIND ENERGY↗

Aerosol Chemical Speciation Monitor (ACSM) Composition-Dependent Collection Efficiency (CDCE) Value-Added Product Report

Aerosol particles influence the Earth’s radiation balance directly by absorbing and scattering light and indirectly by influencing cloud formation, properties, and lifetimes. Measurements of aerosol particle optical properties, mass loading, size distributions, microphysical properties, cloud formation properties, and chemical composition are important for understanding the aerosol life cycle and for validating earth system models that predict these quantities. The U.S. Department of Energy Atmospheric Radiation Measurement (ARM) user facility’s Aerosol Observing System (AOS) is a highly instrumented platform that measures many of these aerosol properties in situ. The Aerodyne aerosol chemical speciation monitor (ACSM) is a baseline instrument deployed in the AOS. The ACSM provides a quantitative measurement of aerosol particle chemical composition for non-refractory (operationally defined as components that evaporate on the 600ºC vaporizer) aerosol components in real time. Standard output is the mass concentration of particulate organics, nitrate, sulfate, chloride, and ammonium. One well-known limitation to the accuracy of the ACSM data is in evaluating the fraction of the ambient aerosol particles that are detected by the instrument. This quantity is referred to as the collection efficiency (CE) and is often less than unity. This is attributed to particles rebounding after impaction onto the heated vaporizer rather than being trapped, volatilized rapidly, and detected. Other factors, such as divergence of the aerosol particle beam, may also impact CE, but the particle rebound effect is the dominant factor. Scientists will sometimes assume CE = 0.5 (i.e., one half of sampled particles are detected) for ambient particles collected during field missions. However, parameterizations have been developed that express CE as a function of the measured chemical composition, referred to as the composition-dependent collection efficiency (CDCE). The physical explanation for a CDCE, supported by laboratory studies, is that rebound from the vaporizer is a function of the particle phase, with liquid-like particles “sticking” to the vaporizer and solid-like or crystalline particles “bouncing” from the vaporizer. Thus, particles with compositions are liquids-like under the measurement conditions (e.g., certain organics, acidic particles, particles enriched in nitrate) have CE close to unity while crystalline or solid particles (e.g., deliquesced ammonium sulfate) have lower CE. This value-added product (VAP) implements the procedure described by Middlebrook et al. (2012) to correct the ACSM data for the composition-dependent collection efficiency. Applying this parameterization improves the accuracy of the ACSM data and brings them into better agreement with other co-located aerosol measurements.

54 ENVIRONMENTAL SCIENCES↗

Assessing the drivers of Isoprene SOA: laboratory studies, field observations and modeling (Final Technical Report)

This report describes the research results and products stemming from the funding award DE-SC0018221. The main goals of this project were to develop a quantitative understanding of the factors driving secondary organic aerosol (SOA) formation from biogenic volatile organic compounds, namely isoprene and monoterpenes, using a combination of laboratory measurements, field data analysis, and modeling. Specifically, we aimed to: 1) develop an observationally constrained volatility distribution of trace gases produced from the oxidation of biogenic hydrocarbons that can explain the formation of SOA from in situ observations made during DOE ASR/ARM field campaigns such as HI-SCALE and BAECC 2) use new and existing (e.g. DOE ASR funded SOAFFEE experiments) laboratory chamber experiments to develop detailed parameterizations of isoprene-derived epxoy diol (IEPOX) reactive uptake, including product branching, product volatility, and solubility constants, as well as SOA formation generally from the formation of highly oxygenated organic molecules (HOM) formed from isoprene and monoterpene photo-oxidation that can be directly incorporated into models. 3) collaborate with PNNL modeling teams to incorporate new parameterizations into detailed box models, such as MOSAIC, and regional or Earth System models such as WRF.

54 ENVIRONMENTAL SCIENCES↗

Ellipsoidal Fitting Methodology for Defect Clusters in Gallium Arsenide

In assessing the initial spatial distribution of defects from neutron or heavy ion irradiation, it is useful to have a reliable, automated, and fast-running tool to evaluate characteristic metrics such as the number of sub-clusters or the overall cluster volume. The latter metric, for instance, can be utilized to estimate a reference neutron fluence level at which inter-cluster interaction effects begin to become significant. This paper details a methodology to fit an arbitrarily complex defect map with a set of ellipsoids (one per identified sub-cluster) in which the constituent defects of a sub-cluster are determined using fuzzy degree-of-membership analysis. Specifically, a parameterized model is developed for point defects in gallium arsenide. Cluster volume calculations based on the model are compared against convex hull and single- ellipsoid representations. Results show that the parameterized sub-cluster model begins to deviate from the two reference models at a recoil energy of about 100 keV in GaAs, with the convex hull and single-ellipsoid representations increasingly overestimating the volume thereafter.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Small field campaign: aerosol – ice formation closure pilot study. Final Report

Prediction of atmospheric ice formation from aerosol particles by heterogeneous nucleation represents one of the grand challenges in atmospheric science. Our insufficient predictive understanding of primary ice formation is the reason that climate models typically do not include heterogeneous ice nucleation with subsequent effects on climate uncertainty. Mixed-phase clouds, where supercooled water droplets and ice crystals coexist play globally an important role regulating climate. This is especially the case for the Arctic region that experiences the greatest warming due to climate change compared to other regions in the world. Immersion freezing initiated by ice-nucleating particles (INPs) in supercooled water droplets is recognized as the dominant primary ice formation pathway in mixed-phase cloud regimes. For this reason, it is crucial to evaluate our capability to predict immersion freezing for a given ambient aerosol population. The goal of this project is conducting a field-based pilot study at the U.S. DOE Atmospheric Radiation Measurement (ARM) user facility at Southern Great Plains (SGP) to evaluate our capability to predict the number concentration of aerosol particles that serve as INPs in the immersion freezing mode. Successful prediction of INP number concentrations is also termed “closure”. This field-observational approach represents a first-of-its kind attempt of an aerosol–ice formation closure study (AEROICESTUDY). Very few closure studies related to INPs have been conducted, and to our knowledge, none using robust size-resolved ambient aerosol composition measurements as a starting point. Achievement of aerosol–ice formation closure relies on our ability to characterize the ambient aerosol population with respect to particles size and composition and to determine INP number concentrations for specified freezing temperatures. This requires numerous online and offline instrumentation resulting in this pilot field campaign being a multi-institutional and community-collaborative effort. We chose the ARM SGP megasite for this first aerosol-ice formation closure pilot study due to its significant measurement capabilities available to obtain detailed physical characterization of the local aerosol population including size distribution, mass loading, and chemical composition of non-refractory aerosol particles. The overall objective of this project is to identify ice nucleation parameterizations that produce the most robust predictions of INP numbers and thus are best suited to be included in cloud and climate models. This objective includes the following goals for this field and laboratory-based project: i) What are the crucial aerosol physicochemical property measurements needed to accurately guide ice nucleation representations in models and long-term INP measurements? ii) What level of parameter details needs to be known to achieve aerosol–ice formation closure? iii) What are the leading causes for climate model bias in INP predictions? We found that the advances in our understanding of immersion freezing garnered over the last 20 years allowed us to yield partial and full closures of atmospheric immersion freezing from ambient aerosol particles. When the aerosol population is physicochemically complex and parameterizations for representative INP types are not yet available, we still struggle to accurately predict INP number concentrations. This project clearly demonstrates that with more laboratory and field measurements that are accompanied by particle composition analysis, the necessary datasets to achieve aerosol–ice formation closure for various locations will emerge, thus providing a robust foundation for guiding the representation of INPs in cloud and climate models.

54 ENVIRONMENTAL SCIENCES↗

FY2022 Status Update: A Probabilistic Model for Stress Corrosion Cracking of SNF Dry Storage Canisters

Understanding the potential risk of stress corrosion cracking of spent nuclear fuel dry storage canisters has been identified as a knowledge gap for determining the safety of long-term interim storage of spent nuclear fuel. To address this, the DOE is funding a multi-lab DOE effort to understand the timing, occurrence, and consequences of potential canister SCC. Sandia National Laboratories has developed a probabilistic model for canister penetration by SCC. This model has been continuously updated at SNL since 2014. Model uncertainties are treated using a nested loop structure, where the outer loop accounts for uncertainties due to lack of data and the inner aleatoric loop accounts for uncertainties due to variation in nature. By separating uncertainties into these categories, it is possible to focus future work on reducing the most influential epistemic uncertainties. Several experimental studies have already been performed to improve the modeling approach through expanded process understanding and improved model parameterization. The resulting code is physics-based and intended to inform future work by identifying (1) important modeling assumptions, (2) experimental data needs, and (3) necessary model developments. In this document, several of the sub-models in the probabilistic SCC model have been exercised, and the intermediate results, as the model progresses from one sub-model to the next, are presented. Evaluating the sub-models in this manner provides a better understanding of sub-model outputs and has identified several unintended consequences of model assumptions or parameterizations, requiring updates to the modeling approach. The following updates have been made, and future updates have been identified.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Toward Better Understanding of Microphysical Processes and Resulting Precipitation Physics: A Merger Of Observations and Cloud Models (Final Report)

A large database of global disdrometer observations and a diverse set of simulations from the Regional Atmospheric Modeling System (RAMS) were compared using principal component analysis (PCA) in order to better understand warm and ice-based precipitation processes. The analysis demonstrated that six distinct precipitation groups (PGs) with common characteristics were revealed in both the observations and model simulations. These PGs were defined by similar co-variability of rain parameters. However, the model showed large concentrations of drops with a given size compared to the observations. Detailed investigations determined the parameterization of drop breakup was forcing the drops to an equilibrium size with a stronger constraint than was indicated by the observations. A series of sensitivity studies permutating the rain shape parameter in the RAMS 2-moment bulk scheme demonstrated how the assumed shape of the rain distribution influences the microphysics as well as precipitation characteristics. Results were compared to the same simulations performed with a bin microphysics scheme where the rain shape parameter is freely evolving. In the bin scheme, a wide range of shape parameters were found, with horizontal and vertical variability, in addition to being a function of storm lifetime. The results of this study highlight the limitations of a fixed assumed rain distribution in 2-moment microphysics schemes. Model simulations were used to probe the microphysical origins of the six PGs. Rain budgets from a supercell case showed that several of the groups hypothesized to be associated with strong ice microphysical processes had significant contributions from melting hail, but the complexity of ice and warm-rain processes in the supercell precluded definitive determination of microphysical origins for many of the PGs. A case study from the Mid-latitude Continental Convection and Clouds Experiment (MC3E) compared the model PGs to disdrometer observations, and demonstrated the spatial and temporal variability of the PGs which is not possible from disdrometer point measurements alone. The analyses performed during this project demonstrated how statistical analysis can provide a robust method for comparing model simulations and observations, which ultimately resulted identification of limitations in some current microphysics parameterizations as well as insights into precipitation variability which ultimately benefits both observations and modeling efforts.

54 ENVIRONMENTAL SCIENCES↗

HPC for Fission Modeling in Support of Nuclear Data [Slides]

The presentation discusses the background of the High-Performance Computing (HPC) division at the Los Alamos National Laboratory (LANL), with emphasis on its computational needs and computational techniques. It also discusses the highlights of the HPC and the results of calculations, including energy sharing parameterized from DFT calculations, parameterization used in CGMF calculations, angular momentum from microscopic calculations, and correlations between FF angular momenta. The presentation concludes with an outlook including what remains to be done in the HPC division and how it will benefit nuclear data, specifically how it will provide a better understanding of the fission process.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Illuminating the pathways to carbon liberation: a systems approach to characterizing the consequential unknowns of carbon transformation and loss from thawing permafrost peatlands (Final Report)

The IsoGenie3 Project delivered new systems-level insights into carbon cycling in thawing permafrost landscapes, with an emphasis on methane and carbon dioxide emissions. From >200 samples from the site collected over a decade, co-analyzed for geochemistry and microbiology, the team recovered ~1,500 assembled microbial genomes and ~1,900 viral population genomes, revealing appreciable genetic novelty - from a new highly abundant bacterial phylum, to novel methane consumers and their activities, to rampant viral novelty. IsoGenie3 linked these organisms to carbon compound transformations (which define the cycling of organic matter in soils, and the loss of the greenhouse gases carbon dioxide and methane), and saw that the microbes at each stage of permafrost thaw had different genetic potential to degrade categories of compounds, expressed that genetic potential differently, and actually transformed carbon compounds into greenhouse gases in different ways. IsoGenie 3 identified that some of the thaw-stage differences were due to plant-microbiome relationships; the plant species across the thaw gradient contributed different carbon compounds into the soil, and hosted distinct microbiota (differing among parts of plants as well as species). Lastly, microbes in the saturated post-thaw conditions appeared likely to contribute to the mobilization and toxification of mercury released during thaw. In parallel with ongoing field sampling and analysis, hypotheses arising from field observations were tested via lab incubation experiments. When communities are taken out of their native habitats, they behave differently, and the team first rigorously quantified the magnitude of this effect on microbiome composition and functional capacity, organic matter composition, and gas production; overall the main system processes were maintained in the lab incubations under the conditions tested. Further, the microbial data could inform geochemical reaction network models of those processes. Then, the team ran experiments with additions of compounds, varying temperature, and “live” vs. “dead” peat (the latter having been gamma irradiated, with a few additional variants to control for methodological artifacts). From these, we (a) determined the importance of plant-derived soluble phenolic compounds in bogs’ extraordinary recalcitrance of organic matter, and carbon gas emissions skewed to carbon dioxide; (b) proposed an abiotic ‘tanning’ mechanism, which could contribute to Sphagnum’s inhibitory effect on anaerobic decomposition through alteration of N availability. IsoGenie3 illuminated longer-term and landscape-scale interactions of permafrost thaw and carbon cycling, advancing knowledge of the drivers of methane dynamics not only across in the permafrost-associated peatland (where hydrology and plant communities dictate microbiomes) but also their interconnected lakes (where sediment carbon quality and resident microbiota are determined by position within lake, and lake features). By leveraging observations of site methane dynamics extending well before this project, the team was able to construct a 44-year portrait of the interplay of permafrost thaw, hydrology, vegetation dynamics, and carbon gas emissions, and the doubling of the fully-thawed fens over this time. From the detailed study of this focal site, IsoGenie3 also aimed to improve model representation of these kinds of sites and processes. To improve predictions of methane transformations, we incorporated acetate and isotope dynamics into the ‘DNDC’ biogeochemistry model. In addition, recovered genomes were grouped into ‘functional groups’, i.e. the genomes that perform a specific function of interest, then used to parameterize maximum growth rate and optimum growth temperature (via signatures in their sequence composition) for the BioCrunch model. The BioCrunch model was then in turn used to test the impact of increasing functional resolution of the microbes, on the carbon gas emissions. Lastly for modeling, the ecosys model was parameterized from the microbial and other data, and used to evaluate drivers of e.g. change in methane emissions. Finally, this project also led to the development of a range of new methods and tools, a new metric of organic matter decomposability, as well as a graph-database solution to multidisciplinary data storage and querying. This project’s ongoing analyses at our focal site also contributed to broader advancements in understanding elements of genetic plasticity and methane metabolism, climate change microbiology and community assembly, global peatland geochemistry and Arctic lakes’ roles in climate feedbacks.

54 ENVIRONMENTAL SCIENCES↗