Search NASA⌕ Search

SEARCH · Search NASA

Results for “Parameter estimation using data”

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 793 records · Page 44
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

System and method for modeling and forecasting cyclical demand systems with dynamic controls and dynamic incentives

Systems and methods for modeling and forecasting cyclical demand systems in the presence of dynamic control or dynamic incentives. A method for modeling a cyclical demand system comprises obtaining historical data on one or more demand measurements over a plurality of demand cycles, obtaining historical data on incentive signals over the plurality of demand cycles, constructing a model using the obtained historical data on the one or more demand measurements and the incentive signals, wherein constructing the model comprises specifying a state-space model, specifying variance parameters in the model, and estimating unknown variance parameters.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2021

Sensor enabled data-driven predictive analytics for modeling and control with high penetration of DERs in distribution systems

The electric power grid is undergoing a tremendous transformation due to the increasing penetration of renewable energy resources beginning with wind and more recently with the distributed energy resources (DERs) such as solar and battery storage. DERs have dramatically changed the role of the distribution systems in the overall power grid, and they are expected to contribute a significant portion of power generation in the future. If current trends for DERs continue, system operation and control will need to change dramatically for improved grid reliability and resiliency. As renewable resources increase in penetration, new and challenging operational, planning, and design problems are expected to emerge. Some of the key challenges that arise in the planning and operation of the future grid are: 1) Quantifying the impact of high DER penetration in distribution systems on bulk grid behavior over multiple time scales. 2) Identifying whether a particular DER configuration/settings have a large impact on the overall grid behavior. These challenges can be addressed in an offline manner using detailed T&D grid models and they can also be addressed in an online manner using sensor measurements. In particular, the advancement and planned growth in sensor technology in power grid over various voltage levels provide us with a unique opportunity to tackle these challenges from a data analytic perspective without needing detailed T&D grid models. A few questions that naturally arise when addressing the challenges from DERs using sensor data are: 1) How can we use limited sensor measurements to monitor & control voltage stability and small signal stability of the bulk system? 2) How can we ensure that the developed data analytic methods are robust to data availability and quality issues? 3) How can we compute the developed analytics in a scalable manner using streaming measurements? In this project, we addressed the aforementioned challenges arising from DERs and answered the questions raised above on how to effectively use the sensor measurements to enhance the reliability and performance of the electric grid. Thus, the overarching goal of this project is to develop effective reduced/representative system models from data that make the computational complexity sufficiently manageable so as to be useful to simulate, analyze, and even control complex non-linear power systems dynamics with large penetrations of DERs. In order to achieve the objective, the project team established a four-fold technical approach 1) Formulated a combined transmission-distribution co-simulation framework for data generation and validation, 2) Derived reduced/representative models of power systems based on data-driven methods for efficient computation and appropriate representation of system behavior, 3) Developed data driven characterization of power system behavior based on transfer operator theory, machine learning and optimization for model estimation, 4) Incorporated a scalable data management and processing architecture using distributed Kafka streaming applications that coordinate input data streams to the developed data analytics. The key accomplishments of the project are: 1) Development of a scalable multi-timescale T&D co-simulation framework (both for steady state and for dynamic co-simulation) using commercial solvers (PSSE and GridLAB-D). The steady-state T&D co-simulation interface is shared with our industry partner (PJM). 2) A structured reduced order dynamic model of distribution systems that can represent partial motor stalling along with a systematic procedure to derive the model parameters. 3) A PMU based online method to monitor, localize and mitigate fault-induced delayed voltage recovery using DER reactive support and load control in distribution systems. 4) Development of linear operator based robust methodologies for dynamic state estimation, uncertainty quantification, system identification and trajectory prediction for power system dynamics. 5) An adaptive damping control for utilizing wind energy resources to provide oscillation damping and system stability. 6) Implementation of Kafka-based framework for efficient processing of streaming data using Linux-based local virtual environment.

DER integration↗
DOE OSTI2016

Observation and integrated Earth-system science: A roadmap for 2016–2025

This report is the response to a request by the Committee on Space Research of the International Council for Science to prepare a roadmap on observation and integrated Earth-system science for the coming ten years. Its focus is on the combined use of observations and modelling to address the functioning, predictability and projected evolution of interacting components of the Earth system on timescales out to a century or so. It discusses how observations support integrated Earth-system science and its applications, and identifies planned enhancements to the contributing observing systems and other requirements for observations and their processing. All types of observation are considered, but emphasis is placed on those made from space. The origins and development of the integrated view of the Earth system are outlined, noting the interactions between the main components that lead to requirements for integrated science and modelling, and for the observations that guide and support them. What constitutes an Earth-system model is discussed. Summaries are given of key cycles within the Earth system. The nature of Earth observation and the arrangements for international coordination essential for effective operation of global observing systems are introduced. Instances are given of present types of observation, what is already on the roadmap for 2016–2025 and some of the issues to be faced. Observations that are organized on a systematic basis and observations that are made for process understanding and model development, or other research or demonstration purposes, are covered. Specific accounts are given for many of the variables of the Earth system. The current status and prospects for Earth-system modelling are summarized. The evolution towards applying Earth-system models for environmental monitoring and prediction as well as for climate simulation and projection is outlined. General aspects of the improvement of models, whether through refining the representations of processes that are already incorporated or through adding new processes or components, are discussed. Some important elements of Earth-system models are considered more fully. Data assimilation is discussed not only because it uses observations and models to generate datasets for monitoring the Earth system and for initiating and evaluating predictions, in particular through reanalysis, but also because of the feedback it provides on the quality of both the observations and the models employed. Inverse methods for surface-flux or model-parameter estimation are also covered. Reviews are given of the way observations and the processed datasets based on them are used for evaluating models, and of the combined use of observations and models for monitoring and interpreting the behaviour of the Earth system and for predicting and projecting its future. A set of concluding discussions covers general developmental needs, requirements for continuity of space-based observing systems, further long-term requirements for observations and other data, technological advances and data challenges, and the importance of enhanced international co-operation.

Earth system science↗
DOE OSTI2019

Residential retrofits at scale: opportunity identification, saving estimation, and personalized messaging based on communicating thermostat data

Basic insulation and heating system retrofits of existing homes could achieve annual energy savings of up to $4–5 billion in the USA. However, current US utility energy efficiency (EE) programs are costly and challenging to scale. Customer acquisition occurs primarily through energy bill mailers, mass media, and on-line advertising that lack specificity about particular home retrofit opportunities, expected energy savings, and cost-effectiveness. Specific retrofit opportunities are identified via on-site home energy assessments that can be inconvenient to homeowners, expensive, and of variable accuracy. Our paper discusses using communicating thermostat data to significantly increase the customer uptake of energy conservation measures (ECMs) by identifying homes with the most significant retrofit opportunities, estimating post-retrofit energy savings, and formulating home-specific outreach. We extended our previously developed gray-box building model to identify physical building parameters corresponding to the target retrofit opportunities, i.e., whole-home R value, air leakage characteristics, and heating system efficiency. Here, the estimated R values and air leakage characteristics compare favorably with the ground truth. The validated algorithms can calculate home-specific energy savings estimate for each ECM, and the algorithm outputs can be used by utility EE programs to formulate home-specific retrofit offers.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2023

Robust Statistical Approach for Determination of Graphite Nitridation Using Bayesian Model Comparison

A better estimation of surface reaction efficiency of semiconductor-grade graphite with atomic nitrogen, as well as the calibration error are calculated using Bayesian updating based on experimental data. Compared with a conventional deterministic model, the stochastic model approach is a powerful tool in the sense that the model is capable of taking into account underlying error correlations among the data quantities. In this paper, we investigate four different stochastic models (called “stochastic system model classes” herein) corresponding to different descriptions of modeling and measurement error structures, given one deterministic physical model. These stochastic system model classes differ in the covariance matrix structure that is used in the uncertainty model to represent uncertainties associated with the physical model and experimental measurements. For each model class, Bayesian inference is used to estimate the posterior probabilities of the physical model parameters as well as of the stochastic model parameters. Model comparison and selection are then applied based on two measures including Bayesian evidence and Bayesian information criterion, as well as the deviance information criterion. Both measures suggest the stochastic model class, which considers that a correlation between errors in two data quantities among different data points is the most plausible. With the stochastic model class, the range of uncertainty in surface reaction efficiency is estimated to be about two orders of magnitude at [Formula: see text].

Engineering↗
DOE OSTI2016

Improved accuracy of quantitative parameter estimates in dynamic contrast-enhanced CT study with low temporal resolution

Purpose: A previously proposed method to reduce radiation dose to patient in dynamic contrast-enhanced (DCE) CT is enhanced by principal component analysis (PCA) filtering which improves the signal-to-noise ratio (SNR) of time-concentration curves in the DCE-CT study. The efficacy of the combined method to maintain the accuracy of kinetic parameter estimates at low temporal resolution is investigated with pixel-by-pixel kinetic analysis of DCE-CT data. Methods: The method is based on DCE-CT scanning performed with low temporal resolution to reduce the radiation dose to the patient. The arterial input function (AIF) with high temporal resolution can be generated with a coarsely sampled AIF through a previously published method of AIF estimation. To increase the SNR of time-concentration curves (tissue curves), first, a region-of-interest is segmented into squares composed of 3 × 3 pixels in size. Subsequently, the PCA filtering combined with a fraction of residual information criterion is applied to all the segmented squares for further improvement of their SNRs. The proposed method was applied to each DCE-CT data set of a cohort of 14 patients at varying levels of down-sampling. The kinetic analyses using the modified Tofts’ model and singular value decomposition method, then, were carried out for each of the down-sampling schemes between the intervals from 2 to 15 s. The results were compared with analyses done with the measured data in high temporal resolution (i.e., original scanning frequency) as the reference. Results: The patients’ AIFs were estimated to high accuracy based on the 11 orthonormal bases of arterial impulse responses established in the previous paper. In addition, noise in the images was effectively reduced by using five principal components of the tissue curves for filtering. Kinetic analyses using the proposed method showed superior results compared to those with down-sampling alone; they were able to maintain the accuracy in the quantitative histogram parameters of volume transfer constant [standard deviation (SD), 98th percentile, and range], rate constant (SD), blood volume fraction (mean, SD, 98th percentile, and range), and blood flow (mean, SD, median, 98th percentile, and range) for sampling intervals between 10 and 15 s. Conclusions: The proposed method of PCA filtering combined with the AIF estimation technique allows low frequency scanning for DCE-CT study to reduce patient radiation dose. The results indicate that the method is useful in pixel-by-pixel kinetic analysis of DCE-CT data for patients with cervical cancer.

61 RADIATION PROTECTION AND DOSIMETRY↗
DOE OSTI2017

Bayesian K-SVD Using Fast Variational Inference

Recent work in signal processing in general and image processing in particular deals with sparse representation related problems. Two such problems are of paramount importance: an overriding need for designing a well-suited overcomplete dictionary containing a redundant set of atoms— i.e., basis signals—and how to find a sparse representation of a given signal with respect to the chosen dictionary. Dictionary learning techniques, among which we find the popular K-singular value decomposition algorithm, tackle these problems by adapting a dictionary to a set of training data. A common drawback of such techniques is the need for parameter-tuning. In order to overcome this limitation, we propose a fullyautomated Bayesian method that considers the uncertainty of the estimates and produces a sparse representation of the data without prior information on the number of non-zeros in each representation vector. We follow a Bayesian approach that uses a three-tiered hierarchical prior to enforce sparsity on the representations and develop an efficient variational inference framework that reduces computational complexity. Furthermore, we describe a greedy approach that speeds up the whole process. Lastly, we present experimental results that show superior performance on two different applications with real images: denoising and inpainting.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

A Machine Learning–Based Tire Life Prediction Framework for Increasing Life of Commercial Vehicle Tires

In the commercial freight industry, tire retreading decisions are often conservative due to limited knowledge of a tire’s remaining service life. This practice leads to increased costs and material waste. This paper proposes a machine learning–based approach for estimating tire casing life and retreadability, focusing on usage data rather than wear information. This approach could extend the tire’s lifespan and reduce landfill waste. Data integration from diverse tire casing measurement sources presents challenges, including imbalanced removal data. Our methodology addresses these challenges by using historical inspection, telematics, and finite element modeling (FEM) datasets. We introduce “Tire Casing Energy” as a comprehensive usage input and apply a Variance-Reduction Synthetic Minority Oversampling Technique (VR-SMOTE) for data imbalance rectification. A random forest model is used to estimate the state of the tire casing and the casing removal probability, with Bayesian optimization applied for hyperparameter tuning, enhancing model accuracy. Here, the proposed prediction framework is able to differentiate different truck fleets and tire locations based on their usage parameters. With the aid of this machine learning model, the importance and sensitivity of different tire usage parameters can be obtained, which is beneficial to maximize tire life.

Data balancing↗
DOE OSTI2014

The apogee red-clump catalog: Precise distances, velocities, and high-resolution elemental abundances over a large area of the Milky Way's disk

The Sloan Digital Sky Survey III's Apache Point Observatory Galactic Evolution Experiment (APOGEE) is a high-resolution near-infrared spectroscopic survey covering all of the major components of the Galaxy, including the dust-obscured regions of the inner Milky Way disk and bulge. Here we present a sample of 10,341 likely red-clump stars (RC) from the first two years of APOGEE operations, selected based on their position in color-metallicity-surface-gravity-effective-temperature space using a new method calibrated using stellar evolution models and high-quality asteroseismology data. The narrowness of the RC locus in color-metallicity-luminosity space allows us to assign distances to the stars with an accuracy of 5%-10%. The sample extends to typical distances of about 3 kpc from the Sun, with some stars out to 8 kpc, and spans a volume of approximately 100 kpc{sup 3} over 5 kpc ≲ R ≲ 14 kpc, |Z| ≲ 2 kpc, and –15° ≲ Galactocentric azimuth ≲ 30°. The APOGEE red-clump (APOGEE-RC) catalog contains photometry from the Two Micron All Sky Survey, reddening estimates, distances, line-of-sight velocities, stellar parameters and elemental abundances determined from the high-resolution APOGEE spectra, and matches to major proper motion catalogs. We determine the survey selection function for this data set and discuss how the RC selection samples the underlying stellar populations. We use this sample to limit any azimuthal variations in the median metallicity within the ≈45° azimuthal region covered by the current sample to be ≤0.02 dex, which is more than an order of magnitude smaller than the radial metallicity gradient. This result constrains coherent non-axisymmetric flows within a few kiloparsecs from the Sun.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Thermodynamic characterization of synthetic lead-arsenate apatites with different halogen substitutions

Abstract Thermodynamic parameters have been measured for synthetic analogs of the mimetite-group minerals Pb5(AsO4)3X (X = OH, Cl, Br, I) belonging to the apatite supergroup. Phases precipitated from aqueous solutions under ambient conditions with well characterized structures and compositions were studied. For each phase, dissolution enthalpy was experimentally determined by oxide melt drop solution calorimetry in a molten solvent of sodium molybdate (3Na2O·4MoO3) at 976 K. The enthalpy of formation from the elements ΔHf,elo was calculated using thermochemical cycles and was −3030.6 ± 11.5, −3026.6 ± 15.8, −2967.6 ± 25.0, and −2993.1 ± 12.2 kJ/mol for Pb5.00(AsO4)3.00OH0.86(CO3)0.07, Pb5.00(AsO4)3.00Cl0.80(CO3)0.10, Pb5.00(AsO4)3.00Br0.80(CO3)0.10, and Pb5.00(AsO4)3.00I0.45OH0.35(CO3)0.10, respectively. These ΔHf,elo values exhibit typical trends for apatites: they increased (were less negative) with the increasing molar mass and ionic radius of X and decreased with the electronegativity and ionization energy of X. The compilation and comparison of data for Ca-, Pb-, P-, and As-apatites revealed correlations indicating that thermodynamic enthalpic stability is largely influenced by chemical factors (e.g., differences in electronegativities of the elements, ionization energy, or ionic characteristics of the bonds) and to a lesser extent by physical and geometric parameters in the crystal structure related to the mass and size of the X anion. Using the correlations, it was possible to estimate the value of hitherto unknown ΔHf,elo for Pb5(AsO4)3F, −3144.3 ± 66.5 kJ/mol. The observed relationships apply to the entire apatite supergroup and can be used to predict the values of ΔHf,elo for phases that have not been studied experimentally. The new data on environmentally significant phases will contribute to the modeling of mineral-water interactions, particularly for potential use in the remediation of soils and wastes contaminated with Pb and As and in the immobilization of radioactive waste containing I-129.

Geochemistry & Geophysics↗
DOE OSTI2021

Mathematical modeling and parameter estimation of MFI membranes for para/ortho‐xylene separation

Abstract A mathematical model of a membrane with a thin, oriented, and selective MFI layer, which also includes contributions from defects, pore blockages, support layer, and external mass transfer, was formulated based on the Maxwell‐Stefan equations. It was validated using reported (Kim et al., Angew Chemie Int Ed, 2018, 57:480–485; Jeon et al., Nature, 2017, 543:690–694) para/ortho‐xylene separation data from five MFI membranes. The diffusivities of the xylenes were considered to be the same for all membranes, while the thickness and the density of defects and pore blockage were treated as unique properties of each membrane, and their contributions were estimated by fitting the model to the corresponding separation data. The effects of these properties and the role of permeate pressure on the separation performance were subsequently analyzed. The proposed modeling, parameter estimation, and analysis framework allow one to quantitatively interpret the variation of separation performance, to understand separation bottlenecks, and to provide guidance for designing membranes with desired performance.

Liu, Jingjun↗
DOE OSTI2015

Rayleigh-Taylor-instability evolution in colliding-plasma-jet experiments with magnetic and viscous stabilization

The Rayleigh-Taylor instability causes mixing in plasmas throughout the universe, from micron-scale plasmas in inertial confinement fusion implosions to parsec-scale supernova remnants. The evolution of this interchange instability in a plasma is influenced by the presence of viscosity and magnetic fields, both of which have the potential to stabilize short-wavelength modes. Very few experimental observations of Rayleigh-Taylor growth in plasmas with stabilizing mechanisms are reported in the literature, and those that are reported are in sub-millimeter scale plasmas that are difficult to diagnose. Experimental observations in well-characterized plasmas are important for validation of computational models used to make design predictions for inertial confinement fusion efforts. This dissertation presents observations of instability growth during the interaction between a high Mach-number, initially un-magnetized plasma jet and a stagnated, magnetized plasma. A multi-frame fast camera captures Rayleigh-Taylor-instability growth while interferometry, spectroscopy, photodiode, and magnetic probe diagnostics are employed to estimate plasma parameters in the vicinity of the collision. As the instability grows, an evolution to longer mode wavelength is observed. Comparisons of experimental data with idealized magnetohydrodynamic simulations including a physical viscosity model suggest that the observed instability evolution is consistent with both magnetic and viscous stabilization. These data provide the opportunity to benchmark computational models used in astrophysics and fusion research.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2015

Combining of different data pools for calculating a reliable POD for real defects

Real defects are essential for the evaluation of the reliability of non destructive testing (NDT) methods, especially in relation to the integrity of components. But in most of the cases the amount of available real defects is not sufficient to evaluate the system. Model-assisted and transfer functions are one way to handle that challenge. This study is focused on a combination of different data pools to create a sufficient amount of data for the reliability estimation. A widespread approach for calculating the Probability of Detection (POD) was used on a radiographic testing (RT) method. The highest contrast to noise ratio (CNR) of each indication is usually selected as the signal in the 'â vs. a' (signal-response) approach for RT. By combining real and artificial defects (flat bottom holes, side drill holes, flat bottom squares, notches, etc) in RT the highest signals are close to each other, but the process of creating and evaluating real defects is much more complex. The solution is seen in the combination of real and artificial data using a weighted least square approach. The weights for real or artificial data were based on the importance, the value and the different detection behavior of the different data. For comparison, the alternative combination through the Bayesian Updating was also applied. As verification, a data pool with a large amount of real data was available. In an advanced approach for evaluating the digital RT data, the size of the indication (perpendicular to the X-ray beam) was introduced as additional information. The signal now consists of the CNR and the area of the indication. The detectability is changing depending on the area of the indication, a fact that was ignored in the previous POD calculations for RT. This points out that a weighted least square approach to pool the data might no longer be adequate. The Bayesian Updating of the estimated parameters of the relationship between the signal field (the area of the indication) and the geometry of the defects is seen as the appropriate model to combine the different defect types in a useful and meaningful way. This work was carried out together with the Finnish company for spent nuclear fuel and waste management - Posiva Oy. The digital RT is one of the NDT methods that might be used for the inspection of the weld of the copper canister to be used for the spent nuclear fuel in the Scandinavian concept of final disposal.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2021

Validation of the Community Land Model Version 5 over the Contiguous United States (CONUS) using in situ and remote sensing data sets

The Community Land Model (CLM) is an effective tool to simulate the biophysical and biogeochemical processes and their interactions with the atmosphere. Although CLM Version 5 (CLM5) constitutes various updates in these processes, its performance in simulating energy, water and carbon cycles over the Contiguous United States (CONUS) at scales which land surface changes and hydrometeorological and hydroclimatological applications are more locally relevant is yet to be assessed. In this study, we conducted three simulations at 0.125? during 1979-2018 over the CONUS using different configurations of CLM, namely CLM5-biogeochemistry (CLM5BGC), CLM4.5BGC, and CLM5-satellite phenology (CLM5SP). We validated and compared their simulations against multiple remote-sensed and in-situ datasets. Overall, the parametric and structural updates (e.g., carbon cost for nitrogen uptake, variable soil thickness, dry surface layer) in CLM5 improve its ability in capturing terrestrial biogeochemical dynamics. The low evapotranspiration in CLM5BGC is associated with biases in simulating vegetation phenological characteristics rather than soil water limitations. The mismatch between CLM5BGC-simulated peak leaf area index and reference data can be attributed to CLM5BGC's inability in simulating phenology of trees and grasses. The differences between CLM-simulated irrigation and reference estimates can be attributed to differences between processes represented in models and in reality, and uncertainties in input and validation datasets. Evaluation against observations at small catchments suggest that hydrologic parameters needed to be calibrated to improve simulations of runoff, especially subsurface runoff. Additional efforts are needed to incorporate spatially-distributed plant phenology and physiology parameters and regional-specific agricultural management practices (e.g., planting, harvest).

Cheng, Yanyan↗
DOE OSTI2015

A cosmological exclusion plot: towards model-independent constraints on modified gravity from current and future growth rate data

Most cosmological constraints on modified gravity are obtained assuming that the cosmic evolution was standard ΛCDM in the past and that the present matter density and power spectrum normalization are the same as in a ΛCDM model. Here we examine how the constraints change when these assumptions are lifted. We focus in particular on the parameter Y (also called G{sub eff}) that quantifies the deviation from the Poisson equation. This parameter can be estimated by comparing with the model-independent growth rate quantity fσ{sub 8}(z) obtained through redshift distortions. We reduce the model dependency in evaluating Y by marginalizing over σ{sub 8} and over the initial conditions, and by absorbing the degenerate parameter Ω{sub m,0} into Y. We use all currently available values of fσ{sub 8}(z). We find that the combination Y-circumflex =YΩ{sub m,0}, assumed constant in the observed redshift range, can be constrained only very weakly by current data, Y-circumflex =0.28{sub −0.23}{sup +0.35} at 68% c.l. We also forecast the precision of a future estimation of Y-circumflex in a Euclid-like redshift survey. We find that the future constraints will reduce substantially the uncertainty, Y-circumflex =0.30{sub −0.09}{sup +0.08} , at 68% c.l., but the relative error on Y-circumflex around the fiducial remains quite high, of the order of 30%. The main reason for these weak constraints is that Y-circumflex is strongly degenerate with the initial conditions, so that large or small values of Y-circumflex are compensated by choosing non-standard initial values of the derivative of the matter density contrast. Finally, we produce a forecast of a cosmological exclusion plot on the Yukawa strength and range parameters, which complements similar plots on laboratory scales but explores scales and epochs reachable only with large-scale galaxy surveys. We find that future data can constrain the Yukawa strength to within 3% of the Newtonian one if the range is around a few Megaparsecs. In the particular case of f(R) models, we find that the Yukawa range will be constrained to be larger than 80 Mpc/h or smaller than 2 Mpc/h (95% c.l.), regardless of the specific f(R) model.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2019

Impacts of Spatial Heterogeneity and Temporal Non-Stationarity on Intensity-Duration-Frequency Estimates—A Case Study in a Mountainous California-Nevada Watershed

Changes in extreme precipitation events may require revisions of civil engineering standards to prevent water infrastructures from performing below the designated guidelines. Climate change may invalidate the intensity-duration-frequency (IDF) computation that is based on the assumption of data stationarity. Efforts in evaluating non-stationarity in the annual maxima series are inadequate, mostly due to the lack of long data records and convenient methods for detecting trends in the higher moments. In this study, using downscaled high resolution climate simulations of the historical and future periods under different carbon emission scenarios, we tested two solutions to obtain reliable IDFs under non-stationarity: (1) identify quasi-stationary time windows from the time series of interest to compute the IDF curves using data for the corresponding time windows; (2) introduce a parameter representing the trend in the means of the extreme value distributions. Focusing on a mountainous site, the Walker Watershed, the spatial heterogeneity and variability of IDFs or extremes are evaluated, particularly in terms of the terrain and elevation impacts. We compared observations-based IDFs that use the stationarity assumption with the two approaches that consider non-stationarity. The IDFs directly estimated based on the traditional stationarity assumption may underestimate the 100-year 24-h events by 10% to 60% towards the end of the century at most grids, resulting in significant under-designing of the engineering infrastructure at the study site. Strong spatial heterogeneity and variability in the IDF estimates suggest a preference for using high resolution simulation data for the reliable estimation of exceedance probability over data from sparsely distributed weather stations. Discrepancies among the three IDFs analyses due to non-stationarity are comparable to the spatial variability of the IDFs, underscoring a need to use an ensemble of non-stationary approaches to achieve unbiased and comprehensive IDF estimates.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2016

DREISS: Using State-Space Models to Infer the Dynamics of Gene Expression Driven by External and Internal Regulatory Networks

Gene expression is controlled by the combinatorial effects of regulatory factors from different biological subsystems such as general transcription factors (TFs), cellular growth factors and microRNAs. A subsystem’s gene expression may be controlled by its internal regulatory factors, exclusively, or by external subsystems, or by both. It is thus useful to distinguish the degree to which a subsystem is regulated internally or externally–e.g., how non-conserved, species-specific TFs affect the expression of conserved, cross-species genes during evolution. We developed a computational method (DREISS, dreiss.gerteinlab.org) for analyzing the Dynamics of gene expression driven by Regulatory networks, both External and Internal based on State Space models. Given a subsystem, the “state” and “control” in the model refer to its own (internal) and another subsystem’s (external) gene expression levels. The state at a given time is determined by the state and control at a previous time. Because typical time-series data do not have enough samples to fully estimate the model’s parameters, DREISS uses dimensionality reduction, and identifies canonical temporal expression trajectories (e.g., degradation, growth and oscillation) representing the regulatory effects emanating from various subsystems. To demonstrate capabilities of DREISS, we study the regulatory effects of evolutionarily conserved vs. divergent TFs across distant species. In particular, we applied DREISS to the time-series gene expression datasets of C. elegans and D. melanogaster during their embryonic development. We analyzed the expression dynamics of the conserved, orthologous genes (orthologs), seeing the degree to which these can be accounted for by orthologous (internal) versus species-specific (external) TFs. We found that between two species, the orthologs have matched, internally driven expression patterns but very different externally driven ones. This is particularly true for genes with evolutionarily ancient functions (e.g. the ribosomal proteins), in contrast to those with more recently evolved functions (e.g., cell-cell communication). This suggests that despite striking morphological differences, some fundamental embryonic-developmental processes are still controlled by ancient regulatory systems.

Biochemistry & Molecular Biology↗
DOE OSTI2021

On the estimation of boundary layer heights: a machine learning approach

Abstract. The planetary boundary layer height (zi) is a key parameter used in atmospheric models for estimating the exchange of heat, momentum, and moisture between the surface and the free troposphere. Near-surface atmospheric and subsurface properties (such as soil temperature, relative humidity, etc.) are known to have an impact on zi. Nevertheless, precise relationships between these surface properties and zi are less well known and not easily discernible from the multi-year dataset. Machine learning approaches, such as random forest (RF), which use a multi-regression framework, help to decipher some of the physical processes linking surface-based characteristics to zi. In this study, a 4-year dataset from 2016 to 2019 at the Southern Great Plains site is used to develop and test a machine learning framework for estimating zi. Parameters derived from Doppler lidars are used in combination with over 20 different surface meteorological measurements as inputs to a RF model. The model is trained using radiosonde-derived zi values spanning the period from 2016 through 2018 and then evaluated using data from 2019. Results from 2019 showed significantly better agreement with the radiosonde compared to estimates derived from a thresholding technique using Doppler lidars only. Noteworthy improvements in daytime zi estimates were observed using the RF model, with a 50 % improvement in mean absolute error and an R2 of greater than 85 % compared to the Tucker method zi. We also explore the effect of zi uncertainty on convective velocity scaling and present preliminary comparisons between the RF model and zi estimates derived from atmospheric models.

54 ENVIRONMENTAL SCIENCES↗