Engineering Papers⌕ Search

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.

At least 775 records · Page 43

Investigating Uncertainty in Predicting Carbon Dynamics in North American Biomes: Putting Support-Effect Bias in Perspective

A fundamental strategy in NASA's Earth Observing System's (EOS) monitoring of vegetation and its contribution to the global carbon cycle is to rely on deterministic, process-based ecosystem models to make predictions of carbon flux over large regions. These models are parameterized (that is, the input variables are derived) using remotely sensed images such as those from the Moderate Resolution Imaging Spectroradiometer (MODIS), ground measurements and interpolated maps. Since early applications of these models, investigators have noted that results depend partly on the spatial support of the input variables. In general, the larger the support of the input data, the greater the chance that the effects of important components of the ecosystem will be averaged out. A review of previous work shows that using large supports can cause either positive or negative bias in carbon flux predictions. To put the magnitude and direction of these biases in perspective, we must quantify the range of uncertainty on our best measurements of carbon-related variables made on equivalent areas. In other words, support-effect bias should be placed in the context of prediction uncertainty from other sources. If the range of uncertainty at the smallest support is less than the support-effect bias, more research emphasis should probably be placed on support sizes that are intermediate between those of field measurements and MODIS. If the uncertainty range at the smallest support is larger than the support-effect bias, the accuracy of MODIS-based predictions will be difficult to quantify and more emphasis should be placed on field-scale characterization and sampling. This talk will describe methods to address these issues using a field measurement campaign in North America and "upscaling" using geostatistical estimation and simulation.

Dungan, Jennifer L.↗

Advances in parameterization of cumulus convection in numerical models

Observational studies over the last decade have revealed new sources of difficulty in the already complex problem of parameterization of precipitating cumulus convection in large scale numerical weather prediction models. The structure of mesoscale updrafts and downdrafts and the difficulties they present for cumulus parameterization is reviewed. An approach is proposed which follows Kuo (1974) but includes cumulus and mesoscale downdrafts in the vertical distributions of heating and moistening. The development of a realistic life cycle behavior in the numerical prediction of a mesoscale convective complex will be shown.

Molinari, J.↗

Final Technical Report for "High-resolution global modeling of the effects of subgrid-scale clouds and turbulence on precipitating cloud systems"

The Multiscale Modeling Framework (MMF) embeds a cloud-resolving model in each grid column of a General Circulation Model (GCM). A MMF model does not need to use a deep convective parameterization, and thereby dispenses with the uncertainties in such parameterizations. However, MMF models grossly under-resolve shallow boundary-layer clouds, and hence those clouds may still benefit from parameterization. In this grant, we successfully created a climate model that embeds a cloud parameterization (“CLUBB”) within a MMF model. This involved interfacing CLUBB’s clouds with microphysics and reducing computational cost. We have evaluated the resulting simulated clouds and precipitation with satellite observations. The chief benefit of the project is to provide a MMF model that has an improved representation of clouds and that provides improved simulations of precipitation.

47 OTHER INSTRUMENTATION↗

A study of the sources and sinks of methane and methyl chloroform using a global three-dimensional Lagrangian tropospheric tracer transport model

Sources and sinks of methane and methyl chloroform are investigated using a global three-dimensional Lagrangian tropospheric tracer transport model with parameterized hydroxyl and temperature fields. Using the hydroxyl radical field calibrated to the methyl chloroform observations, the globally averaged release of methane and its spatial and temporal distribution were investigated. Two source function models of the spatial and temporal distribution of the flux of methane to the atmosphere were developed. The first model was based on the assumption that methane is emitted as a proportion of net primary productivity (NPP). The second model identified source regions for methane from rice paddies, wetlands, enteric fermentation, termites, and biomass burning based on high-resolution land use data. The most significant difference between the two models were predictions of methane fluxes over China and South East Asia, the location of most of the world's rice paddies, indicating that either the assumption that a uniform fraction of NPP is converted to methane is not valid for rice paddies, or that NPP is underestimated for rice paddies, or that present methane emission estimates from rice paddies are too high.

Taylor, John A.↗

A Hierarchical Evaluation of Regional Climate Simulations

Global climate models (GCMs) are the primary tools for predicting the evolution of the climate system. Through decades of development, GCMs have demonstrated useful skill in simulating climate at continental to global scales. However, large uncertainties remain in projecting climate change at regional scales, which limit our ability to inform decisions on climate change adaptation and mitigation. To bridge this gap, different modeling approaches including nested regional climate models (RCMs), global stretch-grid models, and global high-resolution atmospheric models have been used to provide regional climate simulations (Leung et al. 2003). In previous efforts to evaluate these approaches, isolating their relative merits was not possible because factors such as dynamical frameworks, physics parameterizations, and model resolutions were not systematically constrained. With advances in high performance computing, it is now feasible to run coupled atmosphere-ocean GCMs at horizontal resolution comparable to what RCMs use today. Global models with local refinement using unstructured grids have become available for modeling regional climate (e.g., Rauscher et al. 2012; Ringler et al. 2013). While they offer opportunities to improve climate simulations, significant efforts are needed to test their veracity for regional-scale climate simulations.

climate↗

Determination of the Fermi level position in dilute magnetic Ga{sub 1-x}Mn{sub x}N films

We report on a combined theoretical and experimental determination of the Fermi level position in wurtzite Ga{sub 1−x}Mn{sub x}N films with x=4% and x=10% as grown by molecular beam epitaxy. By means of ellipsometric measurements, the real part of the frequency-dependent conductivity is determined. An electronic model in the framework of the effective bond-orbital model is parameterized in order to theoretically reproduce the measured transport properties. Predictions for the long-wavelength behaviour as a function of the Fermi level are made. The corresponding density of states obtained in this model is in qualitative agreement with first-principle calculations. The absence of a significant experimental peak in the AC conductivity for small frequencies indicates that the Fermi level lies in a gap between two Mn-related impurity bands in the host band gap.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Stoichiometrically coupled carbon and nitrogen cycling in the MIcrobial-MIneral Carbon Stabilization model version 1.0 (MIMICS-CN v1.0)

Explicit consideration of microbial physiology in soil biogeochemical models that represent coupled carbon–nitrogen dynamics presents opportunities to deepen understanding of ecosystem responses to environmental change. The MIcrobial-MIneral Carbon Stabilization (MIMICS) model explicitly represents microbial physiology and physicochemical stabilization of soil carbon (C) on regional and global scales. Here we present a new version of MIMICS with coupled C and nitrogen (N) cycling through litter, microbial, and soil organic matter (SOM) pools. The model was parameterized and validated against C and N data from the Long-Term Inter-site Decomposition Experiment Team (LIDET; six litter types, 10 years of observations, and 13 sites across North America). The model simulates C and N losses from litterbags in the LIDET study with reasonable accuracy (C: R 2 =0.63; N: R 2 =0.29), which is comparable with simulations from the DAYCENT model that implicitly represents microbial activity (C: R 2 =0.67; N: R 2 =0.30). Subsequently, we evaluated equilibrium values of stocks (total soil C and N, microbial biomass C and N, inorganic N) and microbial process rates (soil heterotrophic respiration, N mineralization) simulated by MIMICS-CN across the 13 simulated LIDET sites against published observations from other continent-wide datasets. We found that MIMICS-CN produces equilibrium values in line with measured values, showing that the model generates plausible estimates of ecosystem soil biogeochemical dynamics across continental-scale gradients. MIMICS-CN provides a platform for coupling C and N projections in a microbially explicit model, but experiments still need to identify the physiological and stoichiometric characteristics of soil microbes, especially under environmental change scenarios.

58 GEOSCIENCES↗

A New Theoretical Framework for Parameterizing Nonequilibrium Fractionation During Evaporation From the Ocean

The evaporation model for water isotopes proposed by Craig and Gordon (1965, https://books.google.co.in/books?id=6wIKAQAAIAAJ) is used in most isotope-enabled atmospheric models for the parameterization of nonequilibrium fractionation during evaporation from the ocean. In this model, one of the most uncertain parameters is the nonequilibrium fractionation factor 𝛼 𝑘 . Many isotope models use the formulation of Merlivat and Jouzel (1979, https://doi.org/10.1029/jc084ic08p05029), which parameterizes 𝛼 𝑘 as a function of wind speed and distinguishes between a smooth and a rough regime to account for the effect of ocean waves. The resulting discontinuity in 𝛼 𝑘 between smooth and rough regimes has been disputed by several empirical studies. Here, we present a new approach to parameterizing 𝛼𝑘 by explicitly accounting for the influence of wave drag on the momentum flux near the surface. Following Cifuentes-Lorenzen et al. (2018, https://doi.org/10.1007/s10546-018-0376-0), we add a third wave-induced component to the total momentum flux, in addition to the viscous and turbulent components, and extend the definition of the eddy viscosity to account for the momentum flux due to waves and turbulent dissipation near the surface. The new scheme predicts a slight decrease of 𝛼 𝑘 with wind speed, similar to the smooth-regime parameterization of Merlivat and Jouzel (1979, https://doi.org/10.1029/jc084ic08p05029). This new parameterization is incorporated into the isotope-enabled Community Atmosphere Model, where it improves the correlation of simulated and measured vapor deuterium excess relative to the default version and a version with constant 𝛼 𝑘 , suggesting that it may be used as a valid representation of fractionation during evaporation from the ocean in future isotope models.

54 ENVIRONMENTAL SCIENCES↗

Simulated Hydrological Dynamics and Coupled Iron Redox Cycling Impact Methane Production in an Arctic Soil

The fate of organic carbon (C) in permafrost soils is important to the climate system due to the large global stocks of permafrost C. Thawing permafrost can be subject to dynamic hydrology, making redox processes an important factor controlling soil organic matter (SOM) decomposition rates and greenhouse gas production. In iron (Fe)-rich permafrost soils, Fe(III) can serve as a terminal electron acceptor, promoting anaerobic respiration of SOM and increasing pH. Current large-scale models of Arctic C cycling do not include Fe cycling or pH interactions. Here, a geochemical reaction model was developed by coupling Fe redox reactions and C cycling to simulate SOM decomposition, Fe(III) reduction, pH dynamics, and greenhouse gas production in permafrost soils subject to dynamic hydrology. In this study we parameterized the model using measured CO 2 and CH 4 fluxes as well as changes in pH, Fe(II), and dissolved organic C concentrations from oxic and anoxic incubations of permafrost soils from polygonal permafrost sites in northern Alaska, United States. In simulations of repeated oxic-anoxic cycles, Fe(III) reduction during anoxic periods enhanced CO 2 production, while the net effect of Fe(III) reduction on cumulative CH 4 fluxes depended on substrate C availability. With lower substrate availability, Fe(III) reduction decreased total CH 4 production by further limiting available substrate. With higher substrate availability, Fe(III) reduction enhanced CH 4 production by increasing pH. Our results suggest that interactions among Fe-redox reactions, pH and methanogenesis are important factors in predicting CH 4 and CO 2 production as well as SOM decomposition rates in Fe-rich, frequently waterlogged Arctic soils.

54 ENVIRONMENTAL SCIENCES↗

Centralized Activities in Support of the Radiative Forcing Model Intercomparison Project

A new paradigm in benchmark absorption-scattering radiative transfer is presented that enables both the globally-averaged and spatially-resolved testing of climate model radiation parameterizations using a novel benchmark computational strategy to uncover persistent sources of biases in the aerosol Instantaneous Radiative Effect (IRE). A proof-of-concept is demonstrated with the GFDL AM4 and CESM 1.2.2 climate models. Instead of prescribing atmospheric conditions and aerosols, as in prior intercomparisons, snapshots of the atmospheric state and aerosol optical properties from the participating models are used as inputs to a highly accurate radiation solver to uncover relevant biases unique to each climate model. These diagnostic results show that the models’ aerosol IRE bias is of the same magnitude as the persistent intermodel spread (~1 W/m2), and varies spatially and with intrinsic aerosol optical properties. The findings underscore the significance of native model error analysis and its ability to diagnose global biases, confirming its fundamental value for the Radiative Forcing Model Intercomparison Project (RFMIP).

58 GEOSCIENCES↗

An algorithm to estimate the heating budget from vertical hydrometeor profiles

A simple algorithm to estimate the latent heating of cloud systems from their vertical hydrometeor profiles is proposed. The derivation as well as the validation of the algorithm is based on output generated by a nonhydrostatic cloud model with parameterized microphysical processes. Mature and decaying stages of a GATE squall-type convective system have been tested. The algorithm-derived heating budget is in reasonable agreement with the budget predicted by the cloud model. The input to the proposed algoritm can be obtained from either a rain retrieval technique based on information from multichannel passive microwave signals or a kinematic cloud model based on information from Doppler radar wind fields and radar reflectivity patterns. Such an application would have significant implications for spaceborne remote sensing and the large-scale weather prediction data assimilation problem.

Tao, Wei-Kuo↗

Parameterization of IR cooling in a middle atmosphere dynamics model. I - Effects on the zonally averaged circulation

A computationally efficient two-stream parameterization of IR cooling is developed for a middle atmosphere dynamics model. The parameterization combines the monochromatic feature of an emissivity formulation with the computational simplicity of a single spatial integral to evaluate the monochromatic transfer equation. Calculations show strong radiative control of the mesopause region and the need for substantial deceleration of the mean zonal winds. Calculated temperatures show good agreement in the summer hemisphere where planetary wave activity is negligible, but in the polar night of the winter hemisphere they are about 30 K too cold, which suggests planetary wave heating must make up the deficit. With observed ozone densities, globally averaged radiative equilibrium temperatures in the 65-75 km region are too cold and an additional heat source about 1 K/day is required.

Apruzese, J. P.↗

Error modeling for surrogates of dynamical systems using machine learning: Machine-learning-based error model for surrogates of dynamical systems

A machine learning–based framework for modeling the error introduced by surrogate models of parameterized dynamical systems is proposed. The framework entails the use of high-dimensional regression techniques (eg, random forests, and LASSO) to map a large set of inexpensively computed “error indicators” (ie, features) produced by the surrogate model at a given time instance to a prediction of the surrogate-model error in a quantity of interest (QoI). This eliminates the need for the user to hand-select a small number of informative features. The methodology requires a training set of parameter instances at which the time-dependent surrogate-model error is computed by simulating both the high-fidelity and surrogate models. Using these training data, the method first determines regression-model locality (via classification or clustering) and subsequently constructs a “local” regression model to predict the time-instantaneous error within each identified region of feature space. We consider 2 uses for the resulting error model: (1) as a correction to the surrogate-model QoI prediction at each time instance and (2) as a way to statistically model arbitrary functions of the time-dependent surrogate-model error (eg, time-integrated errors). We then apply the proposed framework to model errors in reduced-order models of nonlinear oil-water subsurface flow simulations, with time-varying well-control (bottom-hole pressure) parameters. The reduced-order models used in this work entail application of trajectory piecewise linearization in conjunction with proper orthogonal decomposition. Moreover, when the first use of the method is considered, numerical experiments demonstrate consistent improvement in accuracy in the time-instantaneous QoI prediction relative to the original surrogate model, across a large number of test cases. When the second use is considered, results show that the proposed method provides accurate statistical predictions of the time- and well-averaged errors.

97 MATHEMATICS AND COMPUTING↗

Assessing the sensitivity of land-atmosphere coupling strength to boundary and surface layer parameters in the WRF model over Amazon

Modeling tools can be used to diagnose regional land-atmosphere (L-A) coupling strength in the absence of sufficient observations, but subject to uncertainties associated with parameters in model physical parameterizations. Different sensitivity analysis (SA) approaches may lead to different conclusions about the underlying sensitivities. In this study, we quantify simulation uncertainties related to parameter perturbations, and use different approaches to conduct parameter SA on the WRF model pertaining to L-A coupling strength for simulations over the Amazon region. A total of twenty parameters from the Yonsei University (YSU) planetary boundary layer (PBL) and the revised MM5 surface layer (SL) schemes were selected in this analysis. Three different SA methods, the Morris One-at-A-Time (MOAT) method, the Multivariate Adaptive Regression Splines (MARS) method, and the Sobol’ method, were employed to analyze seven WRF-simulated variables and five L-A coupling metrics. Results show that 1) parameter perturbations cause large simulation uncertainties which are comparable to those in the observations; 2) three different SA methods give consistent L-A coupling strength outcomes; 3) six out of the twenty parameters contribute 80%–95% of the total variance in the metrics analyzed, and first-order effects dominate over interaction effects; 4) the twelve variables/metrics of interest show similar sensitivity patterns to the selected parameters, which is consistent across all the methods used. Physical mechanisms for how the sensitive parameters act in determining the L-A coupling strength and associated variables also are illustrated. Our results will help quantifying L-A coupling strength and establishing a basis for parameter calibration over the Amazon region.

Wang, Chen↗

A Flexible Parameterization for Shortwave and Longwave Optical Properties of Ice Crystals and Derived Bulk Optical Properties for Climate Models

We provide a parameterization of the extinction efficiency, single scattering albedo and asymmetry parameter of single ice crystals with any combination of particle volume, projected area, component aspect ratio and crystal distortion at any wavelength between 0.2 and 100 mm. The paramerization is an extension of the one previously published by van Diedenhoven et al. (2014, https://doi.org/10.1175/JAS-D-13-0205.1). In addition, the parameterized optical properties are integrated over size distributions yielding bulk extinction efficiencies, single scattering albedos and asymmetry parameters for large ranges of effective radii, particle component aspect ratios and crystal distortion values. The parameterization of single particle optical properties is evaluated with a reference database. The bulk optical properties are evaluated against the ice model selected for the Moderate Resolution Imaging Spectroradiometer (MODIS) Collection 6 products, for which accurate optical properties are available. Mean absolute errors in parameterized extinction efficiency, asymmetry parameter and single scattering albedo are shown to be 0.0272, 0.00890 and 0.00468, respectively for shortwave wavelengths, while they are 0.0641, 0.0368 and 0.0200 in the longwave. Shortwave and longwave asymmetry parameter and single scattering albedos are shown to vary strongly with particle component aspect ratio and distortion, resulting in substantial variation in shortwave fluxes, but relatively small variations in longwave cloud emissivity. The parameterization and bulk optical properties are made publicly available.

Bastiaan Van Diedenhoven↗

Non-standard neutrino interactions mediated by a light scalar at DUNE

Abstract We investigate the effect on neutrino oscillations generated by beyond-the-standard-model interactions between neutrinos and matter. Specifically, we focus on scalar-mediated non-standard interactions (NSI) whose impact fundamentally differs from that of vector-mediated NSI. Scalar NSI contribute as corrections to the neutrino mass matrix rather than the matter potential and thereby predict distinct phenomenology from the vector-mediated ones. Similar to vector-type NSI, the presence of scalar-mediated neutrino NSI can influence measurements of oscillation parameters in long-baseline neutrino oscillation experiments, with a notable impact on CP measurement in the case of DUNE. Our study focuses on the effect of scalar NSI on neutrino oscillations, using DUNE as an example. We introduce a model-independent parameterization procedure that enables the examination of the impact of all non-zero scalar NSI parameters simultaneously. Subsequently, we convert DUNE’s sensitivity to the NSI parameters into projected sensitivity concerning the parameters of a light scalar model. We compare these results with existing non-oscillation probes. Our findings reveal that the region of the light scalar parameter space sensitive to DUNE is predominantly excluded by non-oscillation probes, especially when considering all nonzero parameters simultaneously for DUNE.

Physics↗

A preliminary regional PBPK model of lung metabolism for improving species dependent descriptions of 1,3-butadiene and its metabolites

1,3-Butadiene (BD), a volatile organic chemical (VOC), is used in synthetic rubber production and other industrial processes. It is detectable at low levels in ambient air as well as in tobacco smoke and gasoline vapors. Inhalation exposures to high concentrations of BD have been associated with lung cancer in both humans and experimental animals, although differences in species sensitivity have been observed. Metabolically active lung cells such as Pulmonary Type I and Type II epithelial cells and club cells (Clara cells) 1 are potential targets of BD metabolite-induced toxicity. Metabolic capacities of these cells, their regional densities, and distributions vary throughout the respiratory tract as well as between species and cell types. Here we present a physiologically based pharmacokinetic (PBPK) model for BD that includes a regional model of lung metabolism, based on a previous model for styrene, to provide species-dependent descriptions of BD metabolism in the mouse, rat, and human. Since there are no in vivo data on BD pharmacokinetics in the human, the rat and mouse models were parameterized to the extent possible on the basis of in vitro metabolic data. Where it was necessary to use in vivo data, extrapolation from rat to mouse was performed to evaluate the level of uncertainty in the human model. A kidney compartment and description of downstream metabolism were also included in the model to allow for eventual use of available urinary and blood biomarker data in animals and humans to calibrate the model for estimation of BD exposures and internal metabolite levels. Results from simulated inhalation exposures to BD indicate that incorporation of differential lung region metabolism is important in describing species differences in pulmonary response and that these differences may have implications for risk assessments of human exposures to BD.

1,3-Butadiene↗