SEARCH · Engineering Papers
Results for “model 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.
Why Does Amazon Precipitation Decrease When Tropical Forests Respond to Increasing CO 2 ?
Earth system models predict a zonal dipole of precipitation change over tropical South America, with decreases over the Amazon and increases over the Andes. Much of this has been attributed to the physiological response of the rainforest to elevated CO 2 , which describes a basin-wide reduction in stomatal conductance and transpiration. While robust in Earth system model experiments, details of the underlying atmospheric mechanism—specifically how it evolves in the context of land-atmosphere interaction and the diurnal cycle—are unresolved. We investigate this using idealized model simulations and find that within 24 hr of a CO 2 increase, changes occur over the Amazon that engender synoptic timescale feedbacks. Decreased evapotranspiration from the rainforest throttles near-surface moisture, inducing a drier, warmer, and deeper boundary layer. Above this, enhanced turbulent diffusivity increases vapor in the lower free troposphere. Together, these processes reduce convective activity and cause immediate decreases in Amazon rainfall. Over the synoptic timescale, these changes leave behind lower tropospheric moisture, which is advected westward by the background jet and increases Andean precipitation. This produces a dipole of precipitation change consistent across global and regional models as well as parameterized and resolved convection, though details are sensitive to model topography and boundary layer formulation. The mechanism reported here stresses the importance of fast timescale processes affecting stability over a period of hours that can influence longer-term vegetation-climate interactions. These results help clarify the Amazon's physiological response to rising CO 2 and provide insight into possible causes of historical model biases and end-of-century uncertainty in this region.
A Vertically Resolved Canopy Improves Chemical Transport Model Predictions of Ozone Deposition to North Temperate Forests
Abstract Dry deposition is the second largest tropospheric ozone (O 3 ) sink and occurs through stomatal and nonstomatal pathways. Current O 3 uptake predictions are limited by the simplistic big‐leaf schemes commonly used in chemical transport models (CTMs) to parameterize deposition. Such schemes fail to reproduce observed O 3 fluxes over terrestrial ecosystems, highlighting the need for more realistic treatment of surface‐atmosphere exchange in CTMs. We address this need by linking a resolved canopy model (1D Multi‐Layer Canopy CHemistry and Exchange Model, MLC‐CHEM) to the GEOS‐Chem CTM and use this new framework to simulate O 3 fluxes over three north temperate forests. We compare results with in situ measurements from four field studies and with standalone, observationally constrained MLC‐CHEM runs to test current knowledge of O 3 deposition and its drivers. We show that GEOS‐Chem overpredicts observed O 3 fluxes across all four studies by up to 2×, whereas the resolved‐canopy models capture observed diel profiles of O 3 deposition and in‐canopy concentrations to within 10%. Relative humidity and solar irradiance are strong O 3 flux drivers over these forests, and uncertainties in those fields provide the largest remaining source of model deposition biases. Flux partitioning analysis shows that: (a) nonstomatal loss accounts for 60% of O 3 deposition on average; (b) in‐canopy chemistry makes only a small contribution to total O 3 fluxes; and (c) the CTM big‐leaf treatment overestimates O 3 ‐driven stomatal loss and plant phytotoxicity in these temperate forests by up to 7×. Results motivate the application of fully online vertically explicit canopy schemes in CTMs for improved O 3 predictions.
Searching for beyond the Standard Model physics using the improved description of 100 Mo $2\nu \beta \beta$ decay spectral shape with CUPID-Mo
The current experiments searching for neutrinoless double-β ($0\nu \beta \beta$) decay also collect large statistics of Standard Model allowed two-neutrino double-β ($2\nu \beta \beta$ ) decay events. These can be used to search for Beyond Standard Model (BSM) physics via $2\nu \beta \beta$ decay spectral distortions. 100 Mo has a natural advantage due to its relatively short half-life, allowing higher $2\nu \beta \beta$ decay statistics at equal exposures compared to the other isotopes. We demonstrate the potential of the dual read-out bolometric technique exploiting a 100 Mo exposure of 1.47 kg years, acquired in the CUPID-Mo experiment at the Modane underground laboratory (France). We set limits on $0\nu \beta \beta$ decays with the emission of one or more Majorons, on $2\nu \beta \beta$ decay with Lorentz violation, and $2\nu \beta \beta$ decay with a sterile neutrino emission. In this analysis, we investigate the systematic uncertainty induced by modeling the $2\nu \beta \beta$ decay spectral shape parameterized through an improved model, an effect never considered before. This work motivates searches for BSM processes in the upcoming CUPID experiment, which will collect the largest amount of $2\nu \beta \beta$ decay events among the next-generation experiments.
Institutional Transformation (IX) 2.5 Building Module Help Manual
The Institutional Transformation (IX) building module is a software tool developed at Sandia National Laboratories to evaluate energy conservation measures (ECMs) on hundreds of DOE-2 building energy models simultaneously. In IX, ECMs can be designed through parameterizing DOE-2 building models and doing further processing via visual basic for applications subroutines. IX provides the functionality to handle multiple building models for different years, which enables incrementally changing a site of hundreds of buildings over time. It also enables evaluation of the effects of changing climate, comparisons between data and modeling results, and energy use of centralized utility buildings (CUBs). IX consists of a Microsoft Excel® user interface, Microsoft Access® database, and Microsoft Excel® CUB build utility whose functionalities are described in detail in this report. In addition to descriptions of the user interfaces, descriptions of every ECM already designed in IX is included.
Is solid copper oxalate a spin chain or a mixture of entangled spin pairs?
Macroscopic assemblies of interacting spins give rise to a broad spectrum of behaviors determined by the spatial arrangement of the magnetic sites and the electronic interactions between them. Compounds of copper (II), in which each copper carries spin 1/2, exhibit a vast variety of physical properties. For antiferromagnetically coupled spin sites, there are two limiting scenarios: spin chains in which the spins can exhibit a long-range order or a mixture of dimers in which the spins within each pair are entangled but do not communicate with the spins from other dimers. In principle, the two types can be distinguished on the basis of experimental observations and modeling using empirically parameterized effective Hamiltonians, but in practice, ambiguity may persist for decades, as is the case for copper oxalate. Here we use high-level ab initio calculations to establish the validity of the nearest-site Heisenberg model and to predict the interaction strength between the magnetic sites. Finally, the computed magnetic susceptibility provides an unambiguous interpretation of magnetic experiments performed throughout half a century, clearly supporting the infinite spin-chain behavior of solid copper oxalate.
A Community Atmosphere Model with Superparameterized Clouds
In 1999, National Center for Atmospheric Research (NCAR) scientists Wojciech Grabowski and Piotr Smolarkiewicz created a "multiscale" atmospheric model in which the physical processes associated with clouds were represented by running a simple high-resolution model within each grid column of a lowresolution global model. In idealized experiments, they found that the multiscale model produced promising simulations of organized tropical convection, which other models had struggled to produce. Inspired by their results, Colorado State University (CSU) scientists Marat Khairoutdinov and David Randall created a multiscale version of the Community Atmosphere Model (CAM). They removed the cloud parameterizations of the CAM, and replaced them with Khairoutdinov's high-resolution cloud model. They dubbed the embedded cloud model a "super-parameterization," and the modified CAM is now called the "SP-CAM." Over the next several years, many scientists, from many institutions, have explored the ability of the SP-CAM to simulate tropical weather systems, the day-night changes of precipitation, the Asian and African monsoons, and a number of other climate processes. Cristiana Stan of the Center for Ocean-Land-Atmosphere Interactions found that the SP-CAM gives improved results when coupled to an ocean model, and follow-on studies have explored the SP-CAM's utility when used as the atmospheric component of the Community Earth System Model. Much of this research has been performed under the auspices of the Center for Multiscale Modeling of Atmospheric Processes, a National Science Foundation (NSF) Science and Technology Center for which the lead institution is CSU.
Predicting future coexistence in a North American ant community
Global climate change will remodel ecological communities worldwide. However, as a consequence of biotic interactions, communities may respond to climate change in idiosyncratic ways. This makes predictive models that incorporate biotic interactions necessary. We show how such models can be constructed based on empirical studies in combination with predictions or assumptions regarding the abiotic consequences of climate change. Specifically, we consider a well-studied ant community in North America. First, we use historical data to parameterize a basic model for species coexistence. Using this model, we determine the importance of various factors, including thermal niches, food discovery rates, and food removal rates, to historical species coexistence. We then extend the model to predict how the community will restructure in response to several climate-related changes, such as increased temperature, shifts in species phenology, and altered resource availability. Interestingly, our mechanistic model suggests that increased temperature and shifts in species phenology can have contrasting effects. Nevertheless, for almost all scenarios considered, we find that the most subordinate ant species suffers most as a result of climate change. More generally, our analysis shows that community composition can respond to climate warming in nonintuitive ways. For example, in the context of a community, it is not necessarily the most heatsensitive species that are most at risk. Our results demonstrate how models that account for niche partitioning and interspecific trade-offs among species can be used to predict the likely idiosyncratic responses of local communities to climate change.
Atmospheric processing of iron in mineral and combustion aerosols: development of an intermediate-complexity mechanism suitable for Earth system models
Abstract. Atmospheric processing of iron in dust and combustion aerosols is simulated using an intermediate-complexity soluble iron mechanism designed for Earth system models. The solubilization mechanism includes both a dependence on aerosol water pH and in-cloud oxalic acid. The simulations of size-resolved total, soluble and fractional iron solubility indicate that this mechanism captures many but not all of the features seen from cruise observations of labile iron. The primary objective was to determine the extent to which our solubility scheme could adequately match observations of fractional iron solubility. We define a semi-quantitative metric as the model mean at points with observations divided by the observational mean (MMO). The model is in reasonable agreement with observations of fractional iron solubility with an MMO of 0.86. Several sensitivity studies are performed to ascertain the degree of complexity needed to match observations; including the oxalic acid enhancement is necessary, while different parameterizations for calculating model oxalate concentrations are less important. The percent change in soluble iron deposition between the reference case (REF) and the simulation with acidic processing alone is 63.8 %, which is consistent with previous studies. Upon deposition to global oceans, global mean combustion iron solubility to total fractional iron solubility is 8.2 %; however, the contribution of fractional iron solubility from combustion sources to ocean basins below 15∘ S is approximately 50 %. We conclude that, in many remote ocean regions, sources of iron from combustion and dust aerosols are equally important. Our estimates of changes in deposition of soluble iron to the ocean since preindustrial climate conditions suggest roughly a doubling due to a combination of higher dust and combustion iron emissions along with more efficient atmospheric processing.
Sensitivity analysis of a parameterization of coagulation in an aqueous transport model
A three-dimensional aqueous transport model, ALGE, was developed at the Savannah River National Laboratory (SRNL) as a tool for emergency response. Recently, coagulation and break-up of sediment and particulate (contaminant adhered to sediment) matter has been added as a new model capability, representing both processes as functions of salinity and turbulence. A sensitivity analysis was performed by comparing time series of sediment and particulate concentrations at various locations when the model’s parameters (particle size and density) are altered. The results show significant differences between experiments in both freshwater and saltwater systems parameters vary (P-values < 0.05). From this sensitivity study, we have concluded the importance of incorporating both salinity and turbulent effects, as they have a large role in the physical process of coagulation and breakup of sediment. By accounting for these physical processes, we have made the ALGE model a more robust tool for emergency response.
Time–series Analysis of Volume Change at Brady Hot Springs, Nevada, USA using Geodetic Data from 2003 – 2018
Brady Hot Springs geothermal field has exhibited subsidence, as measured by interferometric synthetic aperture radar (InSAR). Previous studies have examined both the temporal evolution of the deformation from 2004 through 2016 and the spatial extent of the deformation, directly relating the observed subsidence to volumetric changes below the surface. We extend the modeling at Brady to analyze a data set of interferometric pairs spanning from the end of 2003 through 2018. We examine spatial and temporal trends in the observed deformation by time–series analysis of each of the 1656 cubic voxels in a parameterized elastic dislocation model to identify areas where the subsurface volume changes as a function of time. Joint time–series analysis of Global Positioning System and InSAR pairs confirm significant changes in rates of volume change during time intervals when well operations were varied. Here, the rate of subsidence increases with increased injection, consistent with the identification of thermal contraction of the rock matrix as the dominant driving mechanism. Conversely, the modeled volume increases when pumping ceases, suggesting thermal expansion of the rock matrix.
A Shallow‐Deep Unified Stochastic Mass Flux Cumulus Parameterization in the Single Column Community Climate Model
Explore the source record for details and available documents.
Soil type and texture impacts on soil organic carbon accumulation in a sub-tropical agro-ecosystem
Soil organic carbon (C) plays a fundamental role in tropical and subtropical soil fertility, agronomic productivity, and soil health. As a tool for understand ecosystems dynamics, mathematical models such as Century have been used to assess soil's capacity to store C in different environments. However, as Century was initially developed for temperate ecosystems, several authors have hypothesized that C storage may be underestimated by Century in Oxisols. We tested the hypothesis that Century model can be parameterized for tropical soils and used to reliably estimate soil organic carbon (SOC) storage. The aim of this study was to investigate SOC storage under two soil types and three textural classes and quantify the sources and magnitude of uncertainty using the Century model. The simulation for SOC storage was efficient and the mean residue was 10 Mg C ha -1 (13%) for n = 91. However, a different simulation bias was observed for soil with <600 g kg -1 of clay was 16.3 Mg C ha -1 (18%) for n = 30, and at >600 g kg -1 of clay, was 4 Mg C ha -1 (5%) for n = 50, respectively. The results suggest a non-linear effect of clay and silt contents on C storage in Oxisols. All types of soil contain nearly 70% of Fe and Al oxides in the clay fraction and a regression analysis showed an increase in model bias with increase in oxides content. Consequently, inclusion of mineralogical control of SOC stabilization by Fe and Al (hydro) oxides may improve results of Century model simulations in soils with high oxides contents
Three-dimensional potential energy surface of Ar–CO
A three-dimensional intermolecular potential energy surface of the Ar–CO complex has been determined by fitting most of the previously reported spectroscopic data, where observed transition frequencies by microwave, millimeter-wave, submillimeter-wave, and infrared spectroscopy were reproduced simultaneously within their experimental accuracies. A free rotor model Hamiltonian considering all the freedom of motions for an atom-diatom system was applied to calculate vibration-rotation energies. A three-dimensional potential energy surface obtained by ab initio calculations at the CCSD(T)-F12b/aug-cc-pV5Z level of theory was parameterized by a model function consisting of 46 parameters. They were used as initial values for the least-squares analysis of the experimental data. A total of 20 parameters were optimized to reproduce all the spectroscopic data.
Assessment and Simulation of Particulate Transport for Delivery of Solid Amendments into the Subsurface: FY25 Status Report
For particulate-based amendments to be viable for field-scale remediation at the Hanford Site (e.g., the 200 DV-1 Operable Unit), amendment particles need to be delivered a reasonable distance away from an injection well to provide cost effective in situ treatment. Field-scale particle transport models can estimate spatial deposition of amendment particles in the subsurface, which is critical for developing an overall remediation strategy. However, field-scale particle simulations are currently limited due to insufficient simulation capability and a lack of experimental data to validate and parameterize particle transport models. During this fiscal year, the following progress has been made toward a field-scale particle transport modeling evaluation: (1) in addition to the two particle transport models implemented last FY, four additional particle transport models have been implemented within PFLOTRAN; (2) a Python-based pre-screening tool was finalized, enabling users to quickly estimate the particle radius of influence (ROI) for any given particle-amendment system; and (3) an initial compatibility assessment was completed using both particle transport simulations deployed through the pre-screening tool, in conjunction with general guidelines to (a) identify the most suitable amendment particle sizes for various Hanford sediments and (b) evaluate amendment-delivery fluid compatibility. The preliminary compatibility assessment revealed that for amendment delivery success to the various Hanford target formations, amendment particle sizes will likely need to be smaller than the amendment sizes tested in the DV-1 treatability study. It is recommended that amendment particles be decreased in size, or alternative smaller size amendments be obtained from the manufacture, prior to any further experimental testing. Also, preliminary testing suggests that xanthan gum may be the most broadly compatible delivery fluid. Planned laboratory experiments will be instrumental in validating and refining the PFLOTRAN particulate transport model formulations, ultimately enabling predictive capabilities to facilitate the design of field-scale amendment delivery systems. This work consists of acquiring new theoretical or experimental knowledge. The information associated with this report should not be used as design input or operating parameters without additional qualification.
Perspective: Ab initio force field methods derived from quantum mechanics
It is often desirable to accurately and efficiently model the behavior of large molecular systems in the condensed phase (thousands to tens of thousands of atoms) over long time scales (from nanoseconds to milliseconds). In these cases, ab initio methods are difficult due to the increasing computational cost with the number of electrons. A more computationally attractive alternative is to perform the simulations at the atomic level using a parameterized function to model the electronic energy. Many empirical force fields have been developed for this purpose. However, the functions that are used to model interatomic and intermolecular interactions contain many fitted parameters obtained from selected model systems, and such classical force fields cannot properly simulate important electronic effects. Furthermore, while such force fields are computationally affordable, they are not reliable when applied to systems that differ significantly from those used in their parameterization. They also cannot provide the information necessary to analyze the interactions that occur in the system, making the systematic improvement of the functional forms that are used difficult. Ab initio force field methods aim to combine the merits of both types of methods. The ideal ab initio force fields are built on first principles and require no fitted parameters. Ab initio force field methods surveyed in this perspective are based on fragmentation approaches and intermolecular perturbation theory. This perspective summarizes their theoretical foundation, key components in their formulation, and discusses key aspects of these methods such as accuracy and formal computational cost. The ab initio force fields considered here were developed for different targets, and this perspective also aims to provide a balanced presentation of their strengths and shortcomings. Finally, this perspective suggests some future directions for this actively developing area.
Impacts of Vaccination and Severe Acute Respiratory Syndrome Coronavirus 2 Variants Alpha and Delta on Coronavirus Disease 2019 Transmission Dynamics in Four Metropolitan Areas of the United States
To characterize Coronavirus Disease 2019 (COVID-19) transmission dynamics in each of the metropolitan statistical areas (MSAs) surrounding Dallas, Houston, New York City, and Phoenix in 2020 and 2021, we extended a previously reported compartmental model accounting for effects of multiple distinct periods of non-pharmaceutical interventions by adding consideration of vaccination and Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2) variants Alpha (lineage B.1.1.7) and Delta (lineage B.1.617.2). For each MSA, we found region-specific parameterizations of the model using daily reports of new COVID-19 cases available from January 21, 2020 to October 31, 2021. In the process, we obtained estimates of the relative infectiousness of Alpha and Delta as well as their takeoff times in each MSA (the times at which sustained transmission began). The estimated infectiousness of Alpha ranged from 1.1x to 1.4x that of viral strains circulating in 2020 and early 2021. The estimated relative infectiousness of Delta was higher in all cases, ranging from 1.6x to 2.1x. The estimated Alpha takeoff times ranged from February 1 to February 28, 2021. The estimated Delta takeoff times ranged from June 2 to June 26, 2021. In conclusion, estimated takeoff times are consistent with genomic surveillance data.
Beyond crystal balls: crosscutting solutions in global health to prepare for an unpredictable future
Efforts in global heath need to deal not only with current challenges, but also to anticipate new scenarios, which sometimes unfold at lightning speed. Predictive modeling is frequently used to assist planning, but outcomes depend heavily on a subset of critical assumptions, which are mostly hampered by our limited knowledge about the many factors, mechanisms and relationships that determine the dynamics of disease systems, by a lack of data to parameterize and validate models, and by uncertainties about future scenarios. We propose a shift from a focus on the prediction of individual disease patterns to the identification and mitigation of broader fragilities in public health systems. Modeling capabilities should be used to perform “stress tests” on how interrelated fragilities respond when faced with a range of possible or plausible threats of different nature and intensity. This system should be able to reveal crosscutting solutions with the potential to address not only one threat, but multiple areas of vulnerability to future health risks. Actionable knowledge not based on a narrow subset of threats and conditions can better guide policy, build societal resilience and ensure effective prevention in an uncertain world.