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 433 records · Page 24
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Neutrino mass and mass ordering: no conclusive evidence for normal ordering

The extraction of the neutrino mass ordering is one of the major challenges in particle physics and cosmology, not only for its implications for a fundamental theory of mass generation in nature, but also for its decisive role in the scale of future neutrinoless double beta decay experimental searches. It has been recently claimed that current oscillation, beta decay and cosmological limits on the different observables describing the neutrino mass parameter space provide robust decisive Bayesian evidence in favor of the normal ordering of the neutrino mass spectrum. We further investigate these strong claims using a rich and wide phenomenology, with different sampling techniques of the neutrino parameter space. Contrary to the findings of Jimenez et al., no decisive evidence for the normal mass ordering is found. Neutrino mass ordering analyses must rely on priors and parameterizations that are ordering-agnostic: robust results should be regarded as those in which the preference for the normal neutrino mass ordering is driven exclusively by the data, while we find a difference of up to a factor of 33 in the Bayes factors among the different priors and parameterizations exploited here. An ordering-agnostic prior would be represented by the case of parameterizations sampling over the two mass splittings and a mass scale, or those sampling over the individual neutrino masses via normal prior distributions only. In this regard, we show that the current significance in favor of the normal mass ordering should be taken as 2.7σ (i.e. moderate evidence), mostly driven by neutrino oscillation data. Let us stress that, while current data favor NO only mildly, we do not exclude the possibility that this may change in the future. Eventually, upcoming oscillation and cosmological data may (or may not) lead to a more significant exclusion of IO.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Modified gravity constraints from the full shape modeling of clustering measurements from DESI 2024

We present cosmological constraints on deviations from general relativity (GR) from the first-year of clustering observations from the Dark Energy Spectroscopic Instrument (DESI) in combination with other available datasets including the CMB data from Planck with CMB-lensing from Planck and ACT, BBN constraints on the physical baryon density, the galaxy weak lensing and clustering from DESY3 and supernova data from DESY5. We first consider the μ(a,k)–Σ(a,k) modified gravity (MG) parameterization (as well as η(a,k)) in a ΛCDM and a w 0 w a CDM cosmological backgrounds. Using a functional form for time-only evolution gives μ 0 = 0.11 +0.44 -0.54 from DESI(FS+BAO)+BBN and a wide prior on n s . Using DESI(FS+BAO)+CMB+DESY3+DESY5-SN, we obtain μ 0 = 0.05 ± 0.22 and Σ 0 = 0.008 ± 0.045 and similarly μ 0 = 0.02 +0.19 -0.24 and η 0 = 0.09 +0.36 -0.60 , in an ΛCDM background. In w 0 w a CDM we obtain μ 0 = -0.24 +0.32 -0.28 and Σ 0 = 0.006 ± 0.043, consistent with GR, and we still find a preference of the data for a dynamical dark energy with w 0 > -1 and w a < 0. Using functional dependencies in both time and scale gives μ 0 and Σ 0 with a same level of precision as above but other scale MG parameters remain hard to constrain. We then move to binned parameterizations in a ΛCDM background starting with two bins in redshift and obtain, μ 1 = 1.02 ± 0.13, μ 2 = 1.04 ± 0.11, Σ 1 = 1.021 ± 0.029 and Σ 2 = 1.022 +0.027 -0.023 , all consistent with the unity value of GR in the binning formalism. We then extend the analysis to combine two bins in redshift and two in scale giving 8 MG parameters that we find all consistent with GR. We note that we find here that the tension reported in previous studies about Σ 0 being inconsistent with GR when using Planck PR3 data goes away when we use the recent LoLLiPoP+HiLLiPoP likelihoods. As noted in previous studies, this seems to indicate that the tension is indeed related to the CMB lensing anomaly in PR3 which is also resolved when using the recent likelihoods. We then constrain the class of Horndeski theory in the effective field theory of dark energy approach. We consider both EFT-basis and α-basis in the analysis. Assuming a power law parameterization for the EFT function Ω, which controls non-minimal coupling, we obtain Ω 0 = 0.012 +0.001 -0.012 and s 0 = 0.996 +0.54 -0.20 from the combination of DESI(FS+BAO)+DESY5SN+CMB in a ΛCDM background, which are consistent with GR. Similar results are obtained when using the α-basis and assuming no-braiding (α B = 0) giving c M < 1.14 at 95% CL in a ΛCDM background, also in agreement with GR. However, we see a mild yet consistent indication for c B > 0 when α B is allowed to vary which will require further study to determine whether this is due to systematics or new physics.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Spatio-temporal differences in leaf physiology are associated with fire, not drought, in a clonally integrated shrub

Abstract In highly disturbed environments, clonality facilitates plant survival via resprouting after disturbance, resource sharing among interconnected stems and vegetative reproduction. These traits likely contribute to the encroachment of deep-rooted clonal shrubs in tallgrass prairie. Clonal shrubs have access to deep soil water and are typically thought of as relatively insensitive to environmental variability. However, how leaf physiological traits differ among stems within individual clonal shrubs (hereafter ‘intra-clonal’) in response to extreme environmental variation (i.e. drought or fire) is unclear. Accounting for intra-clonal differences among stems in response to disturbance is needed to more accurately parameterize models that predict the effects of shrub encroachment on ecosystem processes. We assessed intra-clonal leaf-level physiology of the most dominant encroaching shrub in Kansas tallgrass prairie, Cornus drummondii, in response to precipitation and fire. We compared leaf gas exchange rates from the periphery to centre within shrub clones during a wet (2015) and extremely dry (2018) year. We also compared leaf physiology between recently burned shrubs (resprouts) with unburned shrubs in 2018. Resprouts had higher gas exchange rates and leaf nitrogen content than unburned shrubs, suggesting increased rates of carbon gain can contribute to recovery after fire. In areas recently burned, resprouts had higher gas exchange rates in the centre of the shrub than the periphery. In unburned areas, leaf physiology remained constant across the growing season within clonal shrubs (2015 and 2018). Results suggest single measurements within a shrub are likely sufficient to parameterize models to understand the effects of shrub encroachment on ecosystem carbon and water cycles, but model parameterization may require additional complexity in the context of fire.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Integrating crop growth models with remote sensing for predicting biomass yield of sorghum

Abstract Plant phenotypes are often descriptive, rather than predictive of crop performance. As a result, extensive testing is required in plant breeding programmes to develop varieties aimed at performance in the target environments. Crop models can improve this testing regime by providing a predictive framework to (i) augment field phenotyping data and derive hard-to-measure phenotypes and (ii) estimate performance across geographical regions using historical weather data. The goal of this study was to parameterize the Agricultural Production Systems sIMulator (APSIM) crop growth models with remote-sensing and ground-reference data to predict variation in phenology and yield-related traits in 18 commercial grain and biomass sorghum hybrids. Genotype parameters for each hybrid were estimated using remote-sensing measurements combined with manual phenotyping in West Lafayette, IN, in 2018. The models were validated in hybrid performance trials in two additional seasons at that site and against yield trials conducted in Bushland, TX, between 2001 and 2018. These trials demonstrated that (i) maximum plant height, final dry biomass and radiation use efficiency (RUE) of photoperiod-sensitive and -insensitive forage sorghum hybrids tended to be higher than observed in grain sorghum, (ii) photoperiod-sensitive sorghum hybrids exhibited greater biomass production in longer growing environments and (iii) the parameterized and validated models perform well in above-ground biomass simulations across years and locations. Crop growth models that integrate remote-sensing data offer an efficient approach to parameterize larger plant breeding populations.

Yang, Kai-Wei↗
DOE OSTI2021

Spatial biases of information influence global estimates of soil respiration: how can we improve global predictions?

Soil respiration (Rs), the efflux of CO2 from soils to the atmosphere, is a major component of the terrestrial carbon cycle, but is poorly constrained from regional to global scales. The global soil respiration database (SRDB) is a compilation of in-situ Rs observations from around the globe that has been consistently updated with new measurements over the past decade. It is unclear whether the addition of data to new versions of the SRDB has produced better-constrained global Rs estimates. We compared two versions of the SRDB (v3.0 n=5173 and v5.0 n=10366) to determine how additional data influenced the global sum and uncertainty of Rs as predicted by a quantile regression forest model. A model parameterized using SRDBv3 yielded a global Rs sum of 88.6 Pg C yr-1, and associated uncertainty of 29.9 (mean absolute error) and 57.9 (standard deviation) Pg C yr-1, whereas parameterization using SRDBv5 yielded 96.5 Pg C yr-1 and associated uncertainty of 30.2 (mean average error) and 73.4 (standard deviation) Pg C yr-1. Empirically estimated global heterotrophic respiration (Rh) from v3 and v5 were 49.9-50.2 (mean 50.1) and 53.3-53.5 (mean 53.4) Pg C yr-1, respectively. Interestingly, SRDBv5’s inclusion of new data from underrepresented regions (e.g., Asia, Africa, South America) thus resulted in overall higher model uncertainty. The SRDBv5 is still biased towards northern latitudes and temperate zones, so we tested that an optimized global distribution of Rs measurements resulted in a global sum of 96.4 ± 21.4 Pg C yr-1 and lower model uncertainty. These results support current global estimates of Rs but highlight spatial inconsistencies that influence model parameterization and interpretation that should be addressed to improve regional- to global-scale predictions.

Stell, Emma↗
DOE OSTI2021

Aerosol–Ice Formation Closure: A Southern Great Plains Field Campaign

Prediction of ice formation in clouds presents one of the grand challenges in the atmospheric sciences. Immersion freezing initiated by ice-nucleating particles (INPs) is the dominant pathway of primary ice crystal formation in mixed-phase clouds, where supercooled water droplets and ice crystals coexist, with important implications for the hydrological cycle and climate. However, derivation of INP number concentrations from an ambient aerosol population in cloud-resolving and climate models remains highly uncertain. We conducted an aerosol–ice formation closure pilot study using a field-observational approach to evaluate the predictive capability of immersion freezing INPs. The closure study relies on collocated measurements of the ambient size-resolved and single-particle composition and INP number concentrations. The acquired particle data serve as input in several immersion freezing parameterizations, which are employed in cloud-resolving and climate models, for prediction of INP number concentrations. We discuss in detail one closure case study in which a front passed through the measurement site, resulting in a change of ambient particle and INP populations. We achieved closure in some circumstances within uncertainties, but we emphasize the need for freezing parameterization of potentially missing INP types and evaluation of the choice of parameterization to be employed. Overall, this closure pilot study aims to assess the level of parameter details and measurement strategies needed to achieve aerosol–ice formation closure. Furthermore, the closure approach is designed to accurately guide immersion freezing schemes in models, and ultimately identify the leading causes for climate model bias in INP predictions.

58 GEOSCIENCES↗
DOE OSTI2020

Relationships between Immersion Freezing and Crystal Habit for Arctic Mixed-Phase Clouds—A Numerical Study

Abstract The number concentration of ice particles in Arctic mixed-phase clouds is a major controlling factor of cloud lifetime. The relationships between ice nucleation mode and ice crystal habit development are not yet constrained by observations. This study uses a habit-predicting microphysical scheme within a 3D large-eddy simulation model to evaluate the relationship between immersion freezing and ice habit in a simulated Arctic mixed-phase cloud case. Three immersion freezing parameterizations are considered: a volume-dependent freezing scheme (VF), a parameterization limited to activated droplets (C-AC), and a parameterization limited to coarse aerosol particles (C-CM). Both C-AC and C-CM are based on classical nucleation theory. The freezing rate with VF is found to be greater in downdraft regions than in updraft regions due to the downdraft having a higher number concentration of large droplets. The C-AC cases show active freezing of small droplets near cloud top, whereas in the C-CM cases, mainly the 8–32- μ m-sized droplets freeze in updraft regions near the cloud base. Because the initial crystal size is assumed to affect the axis ratio of hexagonal plates, the VF cases produce crystals with larger axis ratios, resulting in smaller mode radii than the C-AC cases. In all cases, irregular polycrystals dominate near cloud top and a band-like structure develops within the cloud, which qualitatively agrees with previous observations. In the VF and C-CM cases, unactivated large droplets arising from coarse-mode aerosol particles contributed significantly to the freezing rate, producing an important influence on crystal habit.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Convection–Kelvin Wave Coupling in a Global Convection-Permitting Model

A convectively coupled equatorial Kelvin wave (CCKW) was observed over the equatorial Indian Ocean in early November 2011 during the DYNAMO field campaign. This study examines the structure of the CCKW event using two simulations made using the MPAS model: one with 3-km grid spacing without convective parameterization and another with a 15-km grid and parameterized convection. Both simulations qualitatively capture the observed structure of the CCKW, including its vertical tilt and progression of cloud/precipitation structures. The two simulations, however, differ substantially in the amplitude of the CCKW-associated precipitation. While the 3-km run realistically captures the observed modulation of precipitation by the CCKW, the 15-km simulation severely underestimates its magnitude. To understand the difference between the two MPAS simulations regarding wave–convection coupling within the CCKW, the relationship of precipitation with convective inhibition, saturation fraction, and surface turbulent fluxes is investigated. Results show that the 15-km simulation underestimates the magnitude of the CCKW precipitation peak in association with its unrealistically linear relationship between moisture and precipitation. Precipitation, both in observations and the 3-km run, is predominantly controlled by saturation fraction and this relationship is exponential. In contrast, the parameterized convection in the 15-km run is overly sensitive to convective inhibition and not sensitive enough to environmental moisture. The implications of these results on CCKW theories are discussed.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Cloud-Resolving ICON Simulations of Secondary Ice Production in Arctic Mixed-Phase Stratocumuli Observed during M-PACE

Field measurements and modeling studies suggest that secondary ice production (SIP) may close the gap between observed Arctic ice nucleating particle (INP) concentrations and ice crystal number concentrations ni. Here, we explore sensitivities with respect to the complexity of different INP parameterizations under the premise that ni is governed by SIP. Idealized, cloud-resolving simulations are performed for the marine cold air outbreak cloud deck sampled during the Mixed-Phase Arctic Cloud Experiment (M-PACE) with the Icosahedral Nonhydrostatic (ICON) model. The impact of the droplet shattering (DS) of raindrops and collisional breakup (BR) in addition to the existing Hallet–Mossop rime splintering mechanism were investigated. Overall, 12 different model experiments (12-h runs) were performed and analyzed. Despite the considerable amount of uncertainty remaining with regard to SIP mechanisms and their process representation in numerical models, we conclude from these experiments that (i) only simulations where DS dominates the SIP signal (potentially amplified by BR) capture observed ice-phase and liquid-phase cloud properties, and (ii) SIP events cluster around the convective outflow region and are structurally linked to mesoscale cloud organization. In addition, interactions with primary nucleation parameterizations of varied complexity were investigated. Here, our simulations show that (i) a stable long-lived mixed-phase cloud (MPC) can be maintained in the absence of primary nucleation once SIP is established, (ii) experiments using a computationally more efficient relaxation-based parameterization of primary nucleation are statistically invariant from simulations considering prognostic INP, and (iii) primary nucleation at cloud-top controls the areal extent of the mixed-phase cloud region, and reduces SIP efficacy via DS due to increased depletion of cloud liquid throughout the entire cloud column.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Forecasting Remote Atmospheric Responses to Decadal Kuroshio Stability Transitions

Based on observational estimates and global ocean eddy-resolving coupled retrospective initialized predictions, we show that Kuroshio Extension variability affects rainfall variability along the west coast of North America. We show that the teleconnection between the current undulations and downstream rainfall can lead to improved subseasonal to seasonal predictions of precipitation over California, and we demonstrate that capturing these teleconnections requires coupled systems with sufficient ocean resolution (i.e., eddy-resolving), especially over time scales longer than one season. The improved forecast skill is diagnosed in terms of 35 years of retrospective initialized ensemble forecasts with an ocean eddy-resolving and an ocean eddy-parameterized coupled model. Not only does the ocean eddy-resolving model show sensitivity to Kuroshio Extension variability in terms of western North America precipitation, but the ocean eddy-resolving forecasts also show improved forecast skill compared to the ocean eddy-parameterized model. Additionally, the ocean eddy-parameterized coupled model shows no sensitivity to Kuroshio Extension variability. We also find near-decadal variability associated with a progression of a lower-tropospheric height dipole around the North Pacific and how these height anomalies lead to wind-driven Rossby waves that affect the eddy activity in the Kuroshio Extension with a time lag on the order of four years. This decadal-scale variability (~10 years) opens the possibility of multiyear predictability of western North American rainfall.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2024

Sensitivity of Arctic Clouds to Ice Microphysical Processes in the NorESM2 Climate Model

Abstract Ice formation remains one of the most poorly represented microphysical processes in climate models. While primary ice production (PIP) parameterizations are known to have a large influence on the modeled cloud properties, the representation of secondary ice production (SIP) is incomplete and its corresponding impact is therefore largely unquantified. Furthermore, ice aggregation is another important process for the total cloud ice budget, which also remains largely unconstrained. In this study, we examine the impact of PIP, SIP, and ice aggregation on Arctic clouds, using the Norwegian Earth System Model, version 2 (NorESM2). Simulations with both prognostic and diagnostic PIP show that heterogeneous freezing alone cannot reproduce the observed cloud ice content. The implementation of missing SIP mechanisms (collisional breakup, drop shattering, and sublimation breakup) in NorESM2 improves the modeled ice properties, while improvements in liquid content occur only in simulations with prognostic PIP. However, results are sensitive to the description of collisional breakup. This mechanism, which dominates SIP in the examined conditions, is very sensitive to the treatment of the sublimation correction factor, a poorly constrained parameter that is included in the utilized parameterization. Finally, variations in ice aggregation treatment can also significantly impact cloud properties, mainly through their impact on collisional breakup efficiency. Overall, enhancement in ice production through the addition of SIP mechanisms and the reduction in ice aggregation (in line with radar observations of shallow Arctic clouds) result in enhanced cloud cover and decreased TOA radiation biases, compared to satellite measurements, especially during the cold months. Significance Statement Arctic clouds remain a large source of uncertainty in projections of the future climate due to the poor representation of the microphysical processes that govern their life cycle. Ice formation is among the least understood processes. While it is widely recognized that better constraints on primary ice production (PIP) are needed to improve existing parameterizations, we show that secondary ice production (SIP) and ice aggregation can have also a significant impact on ice number concentrations. Constraining ice formation through the addition of missing SIP mechanisms and reducing ice aggregation can improve the representation of the cloud macrophysical properties and enhance total cloud cover in the Arctic region, which in turn contributes to decreased TOA radiation biases in the cold months.

Meteorology & Atmospheric Sciences↗
DOE OSTI2025

Data Files for Runoff Evaluation in an Earth System Land Model for Permafrost Regions

Modeling of hydrological runoff is essential for accurately capturing spatiotemporal feedbacks within the land–atmosphere system, particularly in sensitive regions such as permafrost landscapes. However, substantial uncertainties persist in the terrestrial runoff parameterization schemes used in Earth system and land surface models. This is particularly true in permafrost regions, where landscape heterogeneity is high and reliable observational data are scarce.This data set includes all files that were produced and applied in the paper Runoff Evaluation in an Earth System Land Model for Permafrost Regions [Xiang et al. in review]. The paper is in review as of July 1 2025 in Geoscientific Model Development (GMD). In this study, we evaluate the performance of runoff parameterization schemes in the Energy Exascale Earth System Model (E3SM) land model (ELM). Our proposed framework leverages simulation results from the Advanced Terrestrial Simulator (ATS), which is a physics-rich integrated surface/subsurface hydrologic model that has been successfully evaluated previously in Arctic tundra regions. We used ATS to simulate runoff from 22 representative hillslopes in the Sagavanirktok River basin, located on the North Slope of Alaska, then compared the output with ELM’s parameterized representation of total runoff. This dataset contains 2 figure image files (*.png, *jpg) that describe the study site and methods, as well as folders (Figure*.zip) that contain the associated data files (*.csv, *.dat) and python code notebooks (*.ipynb) for figures 3-7 in the paper. Jupyter notebook (*.ipynb) files that produce the figure files using the associated data files will run within a python environment configured with Jupyter Lab or Notebook packages.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

A bespoke model of Arctic river basins based on hillslope delineation: Model Archive

This dataset is a model archive of the paper A bespoke model of Arctic river basins based on hillslope delineation (in prep), which introduces a watershed decomposition and parameterization method for large scale permafrost hydrology simulation. With this dataset, this study aims to address the research question: whether a computationally efficient hillslope-based modeling framework can reliably simulate discharge at Arctic river-basin scales. This dataset contains model input and output data for five modeling scenarios at a study site located in the Sagavanirktok River basin. The five modeling scenarios include three modeling cases under temperate conditions using full 3D, decomposed 3D, and decomposed 2D modeling strategies; and two modeling cases under actual Arctic conditions with permafrost using full 3D and decomposed 2D modeling strategies. Simulations were performed using the Advanced Terrestrial Simulator (ATS, v1.6 for three temperate scenarios and v1.5 for two Arctic scenarios), a physics-rich integrated surface–subsurface hydrologic model with cryo-hydrology features. For the three temperate models, simulations were conducted for the period of 10/01/1993 - 09/30/2002; and for the two Arctic models, simulations were conducted for the period of 01/01/1994 - 12/31/2002. To facilitate reproducibility of simulations, all datasets are organized hierarchically. The dataset contains: (1) Mesh files (.exo) for full 3D model, decomposed 3D models, and decomposed 2D models, located in huc/190604020802_gauge15906000/mesh/. Mesh files can be visualized through Paraview or read by Python. (2) Climate forcings (.h5) for full 3D model and decomposed 3D/2D models are located in huc/190604020802_gauge15906000/daymet_onePiece/, and huc/190604020802_gauge15906000/vp_pr_revised_daymet_1980_2006_with_wind/ separately. Accessible by Python. (3) Raw measured gage discharge (.csv) from USGS, located in huc/190604020802_gauge15906000/gaged_basin15906000_discharge_usgs/. Accessible by Python. (4) Delineated subdomain raster (.tif) and shape files (.shp), and the final parameterized results (.npy) for decomposed models, located in huc/190604020802_gauge15906000/data_preprocessed-meshing. Accessible by Python. (5) Temperate models are located in nonpermaf_huc190604020802_gauge15906000/, which includes three cases: decomposed 2D models (inside model_0*-hillslope_*), decomposed 3D models (inside model_1*-subcatchment_*), and full 3D model (inside model_2*-onepiece_*). Two step spin-up results (checkpoint_final.h5) are located in model_*1-*_spinup_steadystate and model_*2-*_spinup_cycle, separately, which are used to initialize real transient models. The input files (.xml) and output results (.dat) of the real transient models are located in model_*3-*_transient/. Especially, for two example hillslope models (ID=-11 and 11), additional h5py files are included in model_03-hillslope_transient/hillslope-11/, model_03-hillslope_transient/hillslope11, model_13-subcatchment_transient/subcatchment-11/, model_13-subcatchment_transient/subcatchment/11, respectively, which are used to plot the saturation figure (Figure 5) in the manuscript. Accessible by Python. (6) Arctic models are located in huc190604020802_gauge15906000/, which includes two cases: decomposed 2D models (inside model_04-hillslope_transient), and full 3D model (inside model_05-onepiece_transient_mannp1_ra). Three step spin-up results (checkpoint_final.h5) are located in model_01-column_freezeup/, model_02-column_spinup/, model_03-hillslope_spinup/, respectively, which are used to initialize real 2D transient hillslope models. The input files (.xml) and output results (.dat) of transient 2D hillslope models are located in model_04-hillslope_transient/. The input files (.xml) and output results (.dat) of the full 3D transient model is located in model_05-onepiece_transient_mannp1_ra/. The full 3D transient model is initialized by model_02-column_spinup/. Accessible by Python. (7) The MOSART routed discharge results (.csv) under Arctic conditions is located in huc190604020802_gauge15906000/MOSART/. Accessible by Python. (8) All Python codes (.py) used to parameterize full 3D model to decomposed 2D models are located in script/. These codes fit with watershed workflow (a watershed delineation tool) v1.4 under the branch gaob/v1.4 from https://github.com/gaobhub/watershed-workflow.git.

EARTH SCIENCE > CRYOSPHERE↗
DOE OSTI2021

Material Models and Credibility for System Level Abnormal Mechanical ModSim Applications

The purpose of this document is to provide evidence for assessing the adequacy of parameterized material models for a collection of materials used in a finite element analyses setting. “Adequacy” is relative to the intended use of the material in particular analyses. The intended application of the material models covered within this document is for system level abnormal mechanical solid mechanics analyses. Generally, material model parameterizations should be valid from temperatures of approximately -50 to 70° C, across a range of strain rates, and (depending on details of the parts involved) large deformations. Each material covered in this document is presented in its own chapter with a common format across materials. Model assumptions, limitations, existing validation results, readiness for use with uncertainty quantification, general usage guidance, and failure considerations are all provided along with specific parameterization inputs suitable for the finite element analysis code Sierra/Solid Mechanics.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Integrated parameter and process learning for hydrologic and biogeochemical modules in Earth System Models

Focus area: Primary focal area #2; secondary focal area #3: Learning about parameters and processes of land surface hydrologic and biogeochemical models in Earth System models by integrating machine learning, physics, and big data. Science challenges: How do we maximally leverage big-data observations to improve hydrobiogeochemical process description and parameterization so that such modules more realistically capture hydrologic and vegetation responses and feedbacks under the future climate? For example, how can we leverage physics, limited observations of vegetation and streamflow to better estimate evapotranspiration, and, relatedly, net primary productivity, especially for drought areas? Vegetation plays a critical role in regional and global water cycles; however, existing vegetation models have failed to predict vegetation response to droughts (McDowell & Xu, 2017) , arctic greening (Keenan & Riley, 2018) , and critical transitions between forest and savanna (Hirota et al., 2011) . These studies suggest that when we build process-based models (PBM) parameterized from regional and global plant traits, we tend to poorly describe plant adaptation and local-scale competition processes. The models and their associated parameters assigned for different regions in the world are not capturing essential heterogeneity in vegetation responses at finer spatial scales. Many parameters of the land surface models control hydrology and vegetation dynamics at the same time. The heterogeneity in vegetation response is a function of (i) plant type, (ii) plant size, (iii) competition and succession, (iv) environmental controls, and (v) local variations due to the unique ecological community that are very difficult to describe (e.g., the size of gaps resulting from fire that facilitated the coexistence of pioneering species). In the demographic models, only factors (i) and (iv) were captured, and plant types were generally described only by leaf phenology and climate zones. With current demographic models, we generally consider more traits to define plant types (i) and calibrate these traits to consider factors (ii), (iii) and (iv); however, it is substantially challenging to scale to regional and global simulations due to trait variations across space (Ali et al., 2016). Moreover, it has been noted that hillslope processes, including ridge-to-valley flow and sunny vs. shady slopes are primary organizers of water, energy, and vegetation (Clark et al., 2015; Fan et al., 2019) . Although gradual improvements in the hydrologic model component in earth system models may reduce this error (at a remarkably slow pace), the long-term, gradual impact of hydrology on plant traits are not well captured. Recent work showed that the hydrologic controls exerted by groundwater and lateral flow are primary regulators of rooting depth (Fan et al., 2017) . Such hydrologic controls have seldom been reflected in vegetation model parameterizations.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

A Lagrangian Study of the Transition from Shallow to Deep Convection Using ASR Observations and Les Simulations (Final Report)

With the support from this grant, we have made the following accomplishments. We have found that a formula for fractional entrainment rate, a key aspect for parameterization moist convection, that works for both shallow and deep convection, which provide a promising direction for unified treatment of entrainment for both shallow and deep convection. We also examined the formulation of vertical velocity and showed that the effect of non-hydrostatic pressure and entrainment on vertical velocity are not constants, as assumed in past studies, but can be computed using a simple model for the cloudy updrafts. These results also apply to both deep and shallow convection, which opens the door for a more physically based parameterization of the updraft vertical velocity. We have also conducted extensive analysis of the RACORO aircraft data and have published the results in Wei et al. (2020). We have made extensive effort scrutinizing the data in detail and examining the various hypotheses with the aid of the LES simulations. This was the first application of the modified Paluch diagram on modern aircraft data and the first comprehensive analysis of mixing in clouds using real data (instead of only LES output) in a very long time. These analyses enhanced the science return of the ASR aircraft measurements. This grant also supported a study of the collisions among cold pools and an effort to build parameterizations of moist convection using tools developed in control theory.

42 ENGINEERING↗
DOE OSTI2021

Testing to Evaluate Processes Expected to Occur during MSR Salt Spill Accidents

Obtaining a license for a new nuclear reactor requires the identification and assessment of the potential consequences of specified accident scenarios, which are achieved using accident progression models. Those models need to be parameterized and validated using experimental data, but existing experimental data addressing processes relevant to molten salt reactor accidents are sparse. Specifically, experimental data that quantify the sensitives of important processes to the initial conditions of the spill, the ambient environment, and the containment features are needed to parameterize individual process models. Integrated experiments that simulate accident scenarios are also needed to provide data for model validation, but these experiments will require the use of proven methods to quantify the processes under evaluation. The overarching objectives of this work are to develop the methods for simulating the targeted processes, to determine the effectiveness of the methods in producing the data required for model development, to generate data that can be used to parameterize individual process models, and to provide key insights into the behavior of spilled molten salt that should be considered in models. Experimental methods were designed to quantify aspects of individual processes expected to occur during or after a molten salt spill accident that will affect the fate of spilled molten salt and the radionuclides within. These processes include 1) molten salt spreading and heat transfer, 2) molten salt flowing and freezing in tubing, 3) stainless steel corrosion kinetics in molten salt, and 4) molten salt splashing and aerosol generation. The initial tests described in this report were conducted using eutectic FLiNaK to demonstrate the test methods, the data that are generated, and the analyses of the data to derive values needed for modeling. The primary variables that were tested include initial salt temperature and the presence of volatile surrogate fission products (e.g., cesium and iodine). The developed methods are shown to be effective in quantifying the desired processes and can be applied to study more complex salt compositions of interest to molten salt reactor developers, a wide range of environmental conditions of interest to modelers, and additional variables relevant to salt spill accidents. The developed methods and insights gained from laboratory tests can also be incorporated in future large-scale integrated tests used to simulate molten salt spill accidents.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2022

Bayesian Optimization Framework for Imperfect Data or Models

Conventional Bayesian optimization methods implicitly assume that the data and model being optimized are “perfect.” This assumption leads to inaccurate posterior probability distribution functions (PDFs) when applied to “imperfect” data or models. The new Bayesian optimization framework presented in this report provides a way to parameterize the effect of imperfections usually encountered in a prior PDF of generalized data or a model on the posterior PDF. The effects of imperfections are parameterized by a set of constraints imposed on the posterior expectation values of deviations between the data and the model and on their covariance matrix elements. A particular set of values for these constraints conveys an evaluator’s best estimate of the effect of imperfections on the corresponding posterior expectation values. When a prior PDF of generalized data is assumed to be normal, an expression for a posterior PDF satisfying an arbitrary set of constraints is derived analytically for linear models. An analogous iterative algorithm is given for nonlinear models. The corresponding posterior PDF should be used to estimate any posterior expectation values in the presence of imperfections parameterized by that set of constraints. A posterior PDF of a conventional Bayesian optimization method is recovered analytically when all evaluator-specified constraints are set to zero (i.e., in the absence of any imperfections). The analytical expressions derived in this report for normal PDFs and linear models were verified numerically by a Metropolis–Hastings Monte Carlo method. The methods presented herein could be applied to any kind of data or models, including differential cross-section data or integral benchmark experiments.

97 MATHEMATICS AND COMPUTING↗