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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 541 records · Page 30

All sources NASA NTRS DOE OSTI

DOE OSTI2018

RHIC polarization for Runs 9-17

This note describes the RHIC polarization measurements for use by the collider experiments for the Run91 and later polarized proton running periods. The measurement procedure is outlined and the resulting polarization parameters are defined. The systematic uncertainties for each step of the procedure are discussed and estimated; when possible the uncertainties are evaluated using the present data. Finally the use of the provided results to determine mean polarization and uncertainty for a data set is described. The results used for this are compiled on the web pages linked at https://wiki.bnl.gov/rhicspin/Results; there, for each year the results are at the link 'Results'.

43 PARTICLE ACCELERATORS↗

DOE OSTI2018

Efficient Generalizable Deep Learning

Stochastic optimization is a fundamental field of research for machine learning. Stochastic gradient descent (SGD) and related methods provide a feasible means to train complicated prediction models over large datasets. SGD, however, does not explicitly address the problem of overfitting, which can lead to predictions that perform poorly on new data. This difference between loss performance on unseen testing data verses that of training data defines the generalization gap of a model. We introduce a new computational kernel called Stochastic Hessian Projection (SHP) that uses a maximum likelihood framework to simultaneously estimate gradient noise covariance and local curvature of the loss function. Our analysis illustrates that these quantities affect the evolution of parameter uncertainty and therefore generalizability. We show how these computations allow us to predict the generalization gap without requiring holdout data. Explicitly assessing this metric for generalizability during training may improve machine learning predictions when data is scarce and understanding prediction variability is critical.

42 ENGINEERING↗

DOE OSTI2023

Assessing the nexus between groundwater and solar energy plants in a desert basin with a dual-model approach under uncertainty

Globally, many solar power plants and other types of renewable energy are being located in water-scarce regions. Many projects rely on groundwater resources whose sustainability is uncertain. In the Chuckwalla Basin in California, quantification of recharge and trans-valley underflow is needed to estimate the impacts of solar project withdrawals on the water table. However, such estimates are highly challenging due to data scarcity, heterogeneous soils and long residence times. Conventional assessment employs isolated groundwater models configured with crude and uniform estimates of recharge. Here, we employ a data-constrained surface subsurface processes model, PAWS+CLM, to provide an ensemble of recharges and underflows with perturbed parameters. Then, the Parameter Estimation (PEST) package is used to calibrate MODFLOW aquifer conductivity and filter out implausible recharges. The novel dual-model approach, potentially applicable in other arid regions, can effectively assimilate groundwater head observations, reject unrealistic parameters, and narrow the range of estimated drawdowns. Simulated recharge concentrates along alluvial fans at the mountain fronts and ephemeral washes where run-off water infiltrates. If an evenly distributed recharge was assumed, it resulted in under-estimated drawdown and larger uncertainty bounds. The withdrawals are approaching total inflow, suggesting the system will be nearing, if not exceeding, its sustainable groundwater production capacity, and a boom of such projects will not be sustainable. Especially, the cost/benefit of pumped-storage projects is called into question as the initial-fill phase depletes entire area’s recharge. Our study highlights the stress on groundwater resources of solar development, and that the speed of groundwater recovery does not indicate sustainability. Main point 1: A novel dual model approach, involving an integrated surface/subsurface model and a groundwater parameter-estimation model, was able to better constrain the model. Main point 2: The groundwater system may be nearing, if not exceeding, its sustainable groundwater production capacity and the speed of recovery is not indicative of sustainability. Main point 3: Results from using conventionally-assumed uniform recharge distort calibrated K fields and impacts assessment

14 SOLAR ENERGY↗

DOE OSTI2019

Equation Discovery Using Fast Function Extraction: a Deterministic Symbolic Regression Approach

Advances in machine learning (ML) coupled with increased computational power have enabled identification of patterns in data extracted from complex systems. ML algorithms are actively being sought in recovering physical models or mathematical equations from data. This is a highly valuable technique where models cannot be built using physical reasoning alone. In this paper, we investigate the application of fast function extraction (FFX), a fast, scalable, deterministic symbolic regression algorithm to recover partial differential equations (PDEs). FFX identifies active bases among a huge set of candidate basis functions and their corresponding coefficients from recorded snapshot data. This approach uses a sparsity-promoting technique from compressive sensing and sparse optimization called pathwise regularized learning to perform feature selection and parameter estimation. Furthermore, it recovers several models of varying complexity (number of basis terms). FFX finally filters out many identified models using non-dominated sorting and forms a Pareto front consisting of optimal models with respect to minimizing complexity and test accuracy. Numerical experiments are carried out to recover several ubiquitous PDEs such as wave and heat equations among linear PDEs and Burgers, Korteweg–de Vries (KdV), and Kawahara equations among higher-order nonlinear PDEs. Additional simulations are conducted on the same PDEs under noisy conditions to test the robustness of the proposed approach.

42 ENGINEERING↗

DOE OSTI2021

A Rate Theory Model of Radiation-Induced Swelling in an Austenitic Stainless Steel

Many rate theory models of cavity (void) swelling have been published over the past 50 years, all having the same, or similar, structures. A rigorous validation of the models has not been possible because of the dearth of information concerning the microstructures that correspond with the swelling data. Whereas the lack of microstructure information is still an issue for historical swelling data, in the past 10–20 years data have been published on the evolution of the microstructure (point defect yields from collision cascades, cavity number densities, and dislocation densities/yield strengths) allowing certain gaps in information to be filled when considering historic swelling data. With reasonable estimates of key microstructure parameters, a standard rate theory model can be applied, and the model parameter space explored, in connection with historical swelling data. By using published data on: (i) yield strength as a function of dose and temperature (to establish an empirical expression for dislocation density evolution); (ii) cavity number densities as a function of temperature; and (iii) freely migrating defect (FMD) production as a function of primary knock-on atom (PKA) spectrum, the necessary parameter and microstructure inputs that were previously unknown can be used in model development. This paper describes a rate-theory model for void swelling of 316 stainless steel irradiated in the EBR-2 reactor as a function of irradiation temperature and neutron dose.

36 MATERIALS SCIENCE↗

DOE OSTI2017

Assessing the nexus between groundwater and solar energy plants in a desert basin with a dual-model approach under uncertainty

14 SOLAR ENERGY↗

DOE OSTI2017

Fitting of Failure Rate Data to a Closed Form Solution using Gamma-Poisson Method of Moments

Due to the responsibility and gravity of estimation of nuclear power plant failure rates, it is required that there be consistent and accurate methodology for describing the likelihood of an event occurring. However, the methodology currently employed can vary dependent on the data, as well as subjective conjecture from the expert performing the analysis. The current implementation of the empirical Bayes method to a Gamma-Poisson (GaP) distribution utilizes algorithms to solve for the parameters that do not provide consistent answers. Additionally, to achieve a distribution of the likelihood, the variance of each parameter of the GaP must be determined. There does not exist an exact solution to one of the parameters variance and is typically estimated using techniques like the Kass-Steffey adjustment. Thus a new approach to the problem is proposed, built upon the method of moments for a negative binomial, based upon the work of Bradlow. Using the method of moments approach, we are able to achieve a closed form estimation of the mean and variance for each parameter in the negative binomial distribution. Due to the relationship between the negative binominal and GaP distribution, results can be made between the distributions. The hyper-priors defined assume beta prime distribution which is appropriately informed; the results of this application translate back to the Gamma distribution for easy utilization in SAPHIRE. Additionally several cases are explored using publically available NRC data that considers zero inflated, over dispersed, and Poisson data.

Failure Rate↗

DOE OSTI2018

Resource Analysis for Deep Direct-Use Feasibility Study in East Texas

The National Renewable Energy Laboratory, Southern Methodist University Geothermal Laboratory, Eastman Chemical, Turbine Air Systems, and the Electric Power Research Institute are evaluating the feasibility of using geothermal heat to improve the efficiency of natural gas power plants. The area of interest is the Eastman Chemical plant in Longview, Texas, which is on the northwestern margin of the Sabine Uplift. The study is focused on determining the potential for a geothermal reservoir within a 10 km radius of the site as defined by data from existing geological studies and cross-sections within the depths of 2,100 to 3,400 meters. Wells within a 20 km radius are included for broader geological comparison to determine the heat flow, temperature-at-depth, and oil and gas field porosity and permeability. The geothermal reservoir model is based on the multiple formation top data sources, published literature data, and well log interpretations within the 10 km radius. Area thickness estimates, reservoir extent bounding parameters, potential flow rates, and temperatures are combined to calculate a reservoir productivity index and develop a reservoir production model. Historical fluid volumes production data are used as an independent check for the reservoir productivity index and production model results. The reservoir parameters calculated here are being used for the surface engineering model to determine the economic viability of using geothermal fluids for a deep direct use application at this site. The data files are submitted as separate workbooks in 'content model' format, including: Well Fluid Production, Heat Flow, and Geologic Reservoir.

15 GEOTHERMAL ENERGY↗

DOE OSTI2025

Atacama Cosmology Telescope: Multiprobe cosmology with unWISE galaxies and ACT DR6 CMB lensing

We present a joint analysis of the cosmic microwave background (CMB) lensing power spectra measured from the Data Release 6 of the Atacama Cosmology Telescope (ACT) and Planck PR4, cross-correlations between the ACT and Planck lensing reconstruction and galaxy clustering from unWISE, and the unWISE clustering auto-spectrum. We obtain 1.5% constraints on the matter density fluctuations at late times parametrized by the best constrained parameter combination 𝑆$^{3⁢x⁢2⁢pt}_{8}$ ≡ 𝜎 8 ⁢(Ω 𝑚 /0.3) 0.4 = 0.815 ± 0.012. The commonly used 𝑆 8 ≡ 𝜎 8⁢ (Ω 𝑚 /0.3) 0.5 parameter is constrained to 𝑆 8 = 0.816 ± 0.015. In combination with baryon acoustic oscillation (BAO) measurements we find 𝜎 8 = 0.815 ± 0.012. We also present sound-horizon-independent estimates of the present day Hubble rate of 𝐻 0 = 66.4$^{+3.2}_{−3.7}$ km s −1 Mpc −1 from our large scale structure data alone and 𝐻 0 = 64.3$^{+2.1}_{−2.4}$ km s −1 Mpc −1 in combination with uncalibrated supernovae from Pantheon+. Using parametric estimates of the evolution of matter density fluctuations, we place constraints on cosmic structure in a range of high redshifts typically inaccessible with cross-correlation analyses. Combining lensing cross- and autocorrelations, we derive a 3.3% constraint on the integrated matter density fluctuations above 𝑧 = 2.4, one of the tightest constraints in this redshift range and fully consistent with a Λ cold dark matter (Λ⁢CDM) model fit to the primary CMB from Planck. Finally, combining with primary CMB observations and using the extended low redshift coverage of these combined datasets we derive constraints on a variety of extensions to the Λ⁢CDM model including massive neutrinos, spatial curvature, and dark energy. We find in flat Λ⁢CDM⁢ ∑𝑚 𝜈 < 0.12 eV at 95% confidence using the large scale structure data, BAO measurements from Sloan Digital Sky Survey, and primary CMB observations.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2018

Evaluating meteorological data from weather stations, and from satellites and global models for a multi-site epidemiological study

Highlights: • Earth observation data has potential for use in epidemiology but is seldom validated. • The performance of two EO products varied depending on numerous factors. • Epidemiologists should use their own weather stations to record meteorological data. • Daily estimates of meteorological variables are associated with rotavirus infection. • Choice of climate data source should be determined by the research question. Longitudinal and time series analyses are needed to characterize the associations between hydrometeorological parameters and health outcomes. Earth Observation (EO) climate data products derived from satellites and global model-based reanalysis have the potential to be used as surrogates in situations and locations where weather-station based observations are inadequate or incomplete. However, these products often lack direct evaluation at specific sites of epidemiological interest.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2020

Planck 2018 results - VIII. Gravitational lensing

We present measurements of the cosmic microwave background (CMB) lensing potential using the final Planck 2018 temperature and polarization data. Using polarization maps filtered to account for the noise anisotropy, we increase the significance of the detection of lensing in the polarization maps from 5σ to 9σ. Combined with temperature, lensing is detected at 40σ. We present an extensive set of tests of the robustness of the lensing-potential power spectrum, and construct a minimum-variance estimator likelihood over lensing multipoles 8 ≤ L ≤ 400 (extending the range to lower L compared to 2015), which we use to constrain cosmological parameters. We find good consistency between lensing constraints and the results from the Planck CMB power spectra within the ΛCDM model. Combined with baryon density and other weak priors, the lensing analysis alone constrains σ 8 Ωm 0.25 = 0.589 ± 0.020 (1σ errors). Also combining with baryon acoustic oscillation data, we find tight individual parameter constraints, σ8 = 0.811 ± 0.019, H0 = 67.9-1.3+1.2 km s -1 Mpc -1 , and Ω m = 0.303 -0.018 +0.016 . Combining with Planck CMB power spectrum data, we measure σ 8 to better than 1% precision, finding σ 8 = 0.811 ± 0.006. CMB lensing reconstruction data are complementary to galaxy lensing data at lower redshift, having a different degeneracy direction in σ 8 - Ω m space; we find consistency with the lensing results from the Dark Energy Survey, and give combined lensing-only parameter constraints that are tighter than joint results using galaxy clustering. Using the Planck cosmic infrared background (CIB) maps as an additional tracer of high-redshift matter, we make a combined Planck-only estimate of the lensing potential over 60% of the sky with considerably more small-scale signal. We additionally demonstrate delensing of the Planck power spectra using the joint and individual lensing potential estimates, detecting a maximum removal of 40% of the lensing-induced power in all spectra. The improvement in the sharpening of the acoustic peaks by including both CIB and the quadratic lensing reconstruction is detected at high significance.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2017

Constraining the Magmatic System at Mount St. Helens (2004-2008) Using Bayesian Inversion With Physics-Based Models Including Gas Escape and Crystallization

Physics-based models of volcanic eruptions track conduit processes as functions of depth and time. When used in inversions, these models permit integration of diverse geological and geophysical data sets to constrain important parameters of magmatic systems. We develop a 1-D steady state conduit model for effusive eruptions including equilibrium crystallization and gas transport through the conduit and compare with the quasi-steady dome growth phase of Mount St. Helens in 2005. Viscosity increase resulting from pressure-dependent crystallization leads to a natural transition from viscous flow to frictional sliding on the conduit margin. Erupted mass flux depends strongly on wall rock and magma permeabilities due to their impact on magma density. Including both lateral and vertical gas transport reveals competing effects that produce nonmonotonic behavior in the mass flux when increasing magma permeability. Using this physics-based model in a Bayesian inversion, we link data sets from Mount St. Helens such as extrusion flux and earthquake depths with petrological data to estimate unknown model parameters, including magma chamber pressure and water content, magma permeability constants, conduit radius, and friction along the conduit walls. Even with this relatively simple model and limited data, we obtain improved constraints on important model parameters. We find that the magma chamber had low (<5 wt %) total volatiles and that the magma permeability scale is well constrained at ~10 –11.4 m 2 to reproduce observed dome rock porosities. Here, compared with previous results, higher magma overpressure and lower wall friction are required to compensate for increased viscous resistance while keeping extrusion rate at the observed value.

58 GEOSCIENCES↗

DOE OSTI2016

PBPK modeling of the cis- and trans-permethrin isomers and their major urinary metabolites in rats

Permethrin, a pyrethroid insecticide, is suspected to induce neuronal and hormonal disturbances in humans. The widespread exposure of the populations has been confirmed by the detection of the urinary metabolites of permethrin in biomonitoring studies. Permethrin is a chiral molecule presenting two forms, the cis and the trans isomers. Because in vitro studies indicated a metabolic interaction between the trans and cis isomers of permethrin, we adapted and calibrated a PBPK model for trans- and cis-permethrin separately in rats. The model also describes the toxicokinetics of three urinary metabolites, cis- and trans-3-(2,2 dichlorovinyl)-2,2-dimethyl-(1-cyclopropane) carboxylic acid (cis- and trans-DCCA), 3-phenoxybenzoic acid (3-PBA) and 4′OH-phenoxybenzoic acid (4′-OH-PBA). In vivo experiments performed in Sprague–Dawley rats were used to calibrate the PBPK model in a Bayesian framework. The model captured well the toxicokinetics of permethrin isomers and their metabolites including the rapid absorption, the accumulation in fat, the extensive metabolism of the parent compounds, and the rapid elimination of metabolites in urine. Average hepatic clearances in rats were estimated to be 2.4 and 5.7 L/h/kg for cis- and trans-permethrin, respectively. High concentrations of the metabolite 4′-OH-PBA were measured in urine compared to cis- and trans-DCCA and 3-PBA. The confidence in the extended PBPK model was then confirmed by good predictions of published experimental data obtained using the isomers mixture. The extended PBPK model could be extrapolated to humans to predict the internal dose of exposure to permethrin from biomonitoring data in urine. - Highlights: • A PBPK model of isomers of permethrin and its urinary metabolites was developed. • A quantitative link was established for permethrin and its biomarkers of exposure. • The bayesian framework allows getting confidence interval on the estimated parameters. • The PBPK model can be extrapolated to human and used in a reverse dosimetry context.

60 APPLIED LIFE SCIENCES↗

DOE OSTI2018

Quantifying Model Form Uncertainty of Epidemic Forecasting Models from Incidence Data

Objective We present a mathematical framework for non-parametric estimation of the force of infection, together with statistical upper and lower confidence bands. The resulting estimates allow to assess how well simpler models, such as SEIR, fit the observed time series of incidence data.IntroductionUncertainty Quantification (UQ), the ability to quantify the impact of sample-to-sample variations and model misspecification on predictions and forecasts, is a critical aspect of disease surveillance. While quantifying the impact of stochastic uncertainty in the data is well understood, quantifying the impact of model misspecification is significantly harder. For the latter, one needs a "universal model" to which more restrictive parametric models are compared too.MethodsThis talk presents a useful modeling framework for time series of incidence data from contagious diseases that enables one to identify and quantify the impact of model form uncertainty. Specifically, we propose to focus on estimating the timedependent force of infection. The latter is a universal parameters for all contagious disease model. Using a machine learning technique for estimating monotone functions, i.e., isotonic regression and its variants, one can estimate the force ofinfection without addtional assumptions. We note that most contagious disease model do satisfy this monotonicity assumption, due to a combination of factors: depletion of susceptibles, implementation of mitigation strategies, behavior change, etc. Comparing the resulting "non-parametric" estimate with parametric estimates, obtained by fitting an SEIR for example, can reveal model deficiencies and help quantify model form uncertainties.Finally, we discuss how ideas from "strict bound theory" can be used to develop upper and lower uncertainty bands for force of infection that acknowledge the intrinsic stochasticity in the data.ResultsWe demonstrate the application of the methodology to weekly Influenza Like Illness (ILI) incidence data from France andcompare the results to fitted SIR and SEIR models. This comparison can be seen as a nonparametric goodness of fit test, providing one with tools to do simple model selection.ConclusionsWe present a novel and flexible model to statistically describe the force of infection as a function of time. Comparing the fit to incidence data of that model with the fit of simpler parametric models enables the quantification of model form uncertainty and associated model selection.

Hengartner, Nicholas↗

DOE OSTI2020

Fatigue behavior and modeling for additive manufactured 304L stainless steel: The effect of surface roughness

The fatigue strength of additively manufactured metallic parts in their as-built surface condition is mainly dominated by the surface roughness. Post-processing is often inevitable to reduce surface roughness effects even though post-processing diminishes the main advantage of additive manufacturing, which is net-shaped direct-to-service production. This work investigates the underlying mechanisms responsible for fatigue failure of additively manufactured 304L stainless steel parts in as-built and machined/polished surface conditions. Both strain- and force-controlled, fully reversed fatigue tests were conducted to gain a comprehensive understanding of surface roughness effects on fatigue behavior. The sensitivity to surface roughness is shown to be dependent on the control mode, with stress-based fatigue tests showing greater sensitivity than strain-based fatigue tests. Moreover, the fatigue life estimation for as-built specimens was performed based on surface roughness parameters as well as the fatigue properties of the specimens in machined/polished surface condition of the material without using any fatigue data of specimens in as-built surface condition. Accordingly, the effect of surface roughness on the fatigue behavior could be estimated reasonably well in both strain-life and stress-life approaches.

36 MATERIALS SCIENCE↗

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 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↗