Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Model selection”

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 181 records · Page 10
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Tracking Soiling Losses: Assessment, Uncertainty, and Challenges in Mapping

Several models have been presented in the recent years to estimate the magnitude of soiling from environmental parameters. However, these models are often based on data from a single site, or at most a few sites, and only limited data are, as of yet, available on their uncertainty. The present work aims to present a first comparative analysis of soiling estimation models, using measured soiling data from various locations in the USA. The study also investigates the impact that the source of the input data can have on the estimation. The results show that the model selection is only one of the factors that can affect the evaluation. Indeed, the use of satellite-derived or ground-mounted particulate matter data can lead to the generation of different soiling maps, with factors greater than 2x between the modeled losses. The current challenges and the unanswered questions that can bias soiling estimation are discussed. Additionally, potential research directions to improve the quality of soiling modeling are identified.

14 SOLAR ENERGY↗
DOE OSTI2025

Water Temperature, Prey Concentration and Salmonid Density Influence Daily Growth of Wild Juvenile Salmonids in Tributaries of the Upper Salmon River, Idaho ( USA )

ABSTRACT Theory, experiments and field studies indicate that the somatic growth rate of freshwater consumers is shaped by the individual, additive and multiplicative effects of multiple factors, including consumer size and condition, temperature, prey resources and biotic interactions. While our understanding of how these factors affect wild populations of freshwater consumers is improving, the topic remains poorly studied, especially with respect to mobile species. Here, we report on an 8‐year, seven‐stream ( n = 49 stream‐year combinations) observational study examining the individual and interactive effects of invertebrate prey concentration (F, mg/m 3 ), mean daily water temperature (T, °C) and juvenile Chinook salmon ( Oncorhynchus tshawytscha ) density (D, fish/100 m 2 ) on summer daily growth rates (%/d) of mobile, anadromous, juvenile Chinook salmon (age‐0+, n = 382) and sub‐yearling (age‐0+, n = 61) and yearling (age‐1+) steelhead trout ( O. mykiss , n = 70) rearing in cold (mean daily summer: 12.1°C, range: 4.2°C–16.7°C) mountain tributaries of the Salmon River basin in central Idaho (USA). AIC c model selection indicated that daily juvenile salmonid growth positively correlated with water temperature, prey biomass concentration, local juvenile Chinook density and the interaction between water temperature and food but with species and age‐specific differences. Water temperature was a covariate in all top‐ranked models, with daily growth (%/day) rate increasing (0.05%–0.23%/d) linearly with mean daily summer water temperature. In addition to a direct positive relationship with daily growth rate, there was evidence that prey concentration positively interacted with water temperature to accelerate daily growth (F × T). The positive relationship between juvenile salmonid daily growth rate and juvenile Chinook density is difficult to explain and could result from confounding factors. The individual success observed in these streams may contribute to population‐level benefits for the focal consumers, as prey‐rich, warm summers may result in larger individuals with higher energy reserves at the end of the summer/autumn growing season, contributing to improved overwinter survival. Our results, taken in combination with evidence from models, experiments and observational studies, have climate change implications. Current and predicted increases in water temperature will necessitate higher rates of prey consumption by aquatic ectothermic consumers to offset accelerated metabolic demands. Thus, to improve the resilience of mobile freshwater consumers in a warming climate, we suggest that natural resource managers not only consider physical and chemical habitat conditions but also biotic conditions, including the spatiotemporal quantity and quality of prey resources.

Kiffney, Peter. M. [Fish Ecology, Northwest Fisher↗
DOE OSTI2023

Rock Valley Direct Comparison Relocation Working Group Location Results and Recommendations

Work accomplished: Collected and compared historic data for the 1993 Rock Valley earthquake sequence; Compared preliminary and prior location work from different location algorithms, phase pick sets, station constellations, and velocity models; Selected a common set of stations that could be used across all location methods for consistency; Reviewed 8 different sets of phase picks and converged on a single, reviewed set of picks for all common stations; Evaluated four pre-existing regional velocity models and incorporated new and preliminary results for five new velocity models that provide information on the very shallow (< 2km) structure near station RTPP; Compared location results from different methods while using the common sets of picks, stations, and velocity models

58 GEOSCIENCES↗
DOE OSTI2022

Sparsifying priors for Bayesian uncertainty quantification in model discovery

We propose a probabilistic model discovery method for identifying ordinary differential equations governing the dynamics of observed multivariate data. Our method is based on the sparse identification of nonlinear dynamics (SINDy) framework, where models are expressed as sparse linear combinations of pre-specified candidate functions. Promoting parsimony through sparsity leads to interpretable models that generalize to unknown data. Instead of targeting point estimates of the SINDy coefficients, we estimate these coefficients via sparse Bayesian inference. The resulting method, uncertainty quantification SINDy (UQ-SINDy), quantifies not only the uncertainty in the values of the SINDy coefficients due to observation errors and limited data, but also the probability of inclusion of each candidate function in the linear combination. UQ-SINDy promotes robustness against observation noise and limited data, interpretability (in terms of model selection and inclusion probabilities) and generalization capacity for out-of-sample forecast. Sparse inference for UQ-SINDy employs Markov chain Monte Carlo, and we explore two sparsifying priors: the spike and slab prior, and the regularized horseshoe prior. UQ-SINDy is shown to discover accurate models in the presence of noise and with orders-of-magnitude less data than current model discovery methods, thus providing a transformative method for real-world applications which have limited data.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Machine Learning-Based Identification of the Interface Regions for Coupling Local and Nonlocal Models

Local-nonlocal coupling approaches provide a means to combine the computational efficiency of local models and the accuracy of nonlocal models. However, the coupling process can be challenging, requiring expertise to identify the interface between local and nonlocal regions. Here, this study introduces a machine learning-based approach to automatically detect the regions in which the local and nonlocal models should be used. The method uses loading functions evaluated at grid points to decide the model selection at those points. Training of the networks is based on datasets provided by classes of loading functions for which reference coupling configurations are computed using accurate coupled solutions, where accuracy is measured in terms of the relative error between the solution to the coupling approach and the solution to the nonlocal model. We study two approaches that vary in data structure. The first, the full-domain input data approach, uses the entire load vector and outputs a complete label vector, performing a global classification. The second, a window-based approach, processes loads into windows and addresses the problem as a node-wise classification where each window's central point is classified individually. The classification problems are solved via deep learning algorithms based on convolutional neural networks. The performance of these approaches is studied on one-dimensional numerical examples using F1-scores and accuracy metrics. Notably, the windowing approach achieves an accuracy of 0.96 and an F1-score of 0.97, highlighting its potential to automate coupling processes effectively and enhance computational efficiency in material science applications.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

A Bayesian approach to evaluation of soil biogeochemical models

Abstract. To make predictions about the carbon cycling consequences of rising global surface temperatures, Earth system scientists rely on mathematical soil biogeochemical models (SBMs). However, it is not clear which models have better predictive accuracy, and a rigorous quantitative approach for comparing and validating the predictions has yet to be established. In this study, we present a Bayesian approach to SBM comparison that can be incorporated into a statistical model selection framework. We compared the fits of linear and nonlinear SBMs to soil respiration data compiled in a recent meta-analysis of soil warming field experiments. Fit quality was quantified using Bayesian goodness-of-fit metrics, including the widely applicable information criterion (WAIC) and leave-one-out cross validation (LOO). We found that the linear model generally outperformed the nonlinear model at fitting the meta-analysis data set. Both WAIC and LOO computed higher overfitting risk and effective numbers of parameters for the nonlinear model compared to the linear model, conditional on the data set. Goodness of fit for both models generally improved when they were initialized with lower and more realistic steady-state soil organic carbon densities. Still, testing whether linear models offer definitively superior predictive performance over nonlinear models on a global scale will require comparisons with additional site-specific data sets of suitable size and dimensionality. Such comparisons can build upon the approach defined in this study to make more rigorous statistical determinations about model accuracy while leveraging emerging data sets, such as those from long-term ecological research experiments.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Carbon- and Oxygen-rich stars in MaStar: identification and classification

ABSTRACT Carbon- and Oxygen-rich stars populating the Thermally Pulsing Asymptotic Giant Branch (TP-AGB) phase of stellar evolution are relevant contributors to the spectra of ∼1 Gyr old populations. Atmosphere models for these types are uncertain, due to complex molecules and mass-loss effects. Empirical spectra are then crucial, but samples are small due to the short (∼3 Myr) TP-AGB lifetime. Here we exploit the vastness of the MaNGA Stellar library MaStar (∼60 000 spectra) to identify C and O-rich type stars. We define an optical colour selection with cuts of (g − r) > 2 and (g − i) < 1.55(g − r) − 0.07, calibrated with known C- and O-rich spectra. This identifies C- and O-rich stars along clean, separated sequences. An analogue selection is found in V, R, I bands. Our equation identifies C-rich and O-rich spectra with predictive performance metric F1-scores of 0.72 and 0.74 (over 1), respectively. We finally identify 41 C- and 87 O-rich type AGB stars in MaStar, 5 and 49 of which do not have a SIMBAD counterpart. We also detect a sample of non-AGB, dwarf C-stars. We further design a fitting procedure to classify the spectra into broad spectral types, by using as fitting templates empirical C- and O-rich spectra. We find remarkably good fits for the majority of candidates and categorize them into C- and O-rich bins following existing classifications, which correlate to effective temperature. Our selection models can be applied to large photometric surveys (e.g. Euclid, Rubin). The classified spectra will facilitate future evolutionary population synthesis models.

Hill, Lewis (ORCID:0000000180528080)↗
DOE OSTI2023

Ex vivo to in vivo model of malignant peripheral nerve sheath tumors for precision oncology

Malignant peripheral nerve sheath tumors (MPNST) are aggressive soft tissue sarcomas that often develop in patients with neurofibromatosis type 1 (NF1). To address the critical need for novel therapeutics in MPNST, we aimed to establish an ex vivo 3D platform that accurately captured the genomic diversity of MPNST and could be utilized in a medium-throughput manner for drug screening studies to be validated in vivo using patient-derived xenografts (PDX). Genomic analysis was performed on all PDX-tumor pairs. Selected PDX were harvested for assembly into 3D microtissues. Based on prior work in our labs, we evaluated drugs (trabectedin, olaparib, and mirdametinib) ex vivo and in vivo. For 3D microtissue studies, cell viability was the endpoint as assessed by Zeiss Axio Observer. For PDX drug studies, tumor volume was measured twice weekly. Bulk RNA sequencing was performed to identify pathways enriched in cells. We developed 13 NF1-associated MPNST-PDX and identified mutations or structural abnormalities in NF1 (100%), SUZ12 (85%), EED (15%), TP53 (15%), CDKN2A (85%), and chromosome 8 gain (77%). We successfully assembled PDX into 3D microtissues, categorized as robust (>90% viability at 48 h), good (>50%), or unusable (<50%). We evaluated drug response to “robust” or “good” microtissues, namely MN-2, JH-2-002, JH-2-079-c, and WU-225. Drug response ex vivo predicted drug response in vivo, and enhanced drug effects were observed in select models. These data support the successful establishment of a novel 3D platform for drug discovery and MPNST biology exploration in a system representative of the human condition.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2024

Learning of networked spreading models from noisy and incomplete data

Recent years have seen a lot of progress in algorithms for learning parameters of spreading dynamics from both full and partial data. Some of the remaining challenges include model selection under the scenarios of unknown network structure, noisy data, missing observations in time, as well as an efficient incorporation of prior information to minimize the number of samples required for an accurate learning. Here, in this work, we introduce a universal learning method based on a scalable dynamic message-passing technique that addresses these challenges often encountered in real data. The algorithm leverages available prior knowledge on the model and on the data, and reconstructs both network structure and parameters of a spreading model. We show that a linear computational complexity of the method with the key model parameters makes the algorithm scalable to large network instances.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

From noise to information: The transfer function formalism for uncertainty quantification in reconstructing the nuclear density

The neutron distribution of neutron-rich nuclei provides critical information on the structure of finite nuclei and neutron stars. Parity violating experiments—such as PREX and CREX—provide a clean and largely model-independent determination of neutron densities. Such experiments, however, are challenging and expensive, which is why sound statistical arguments are required to maximize the information gained. We introduce a new framework, the “transfer function formalism,” aimed at uncertainty quantification, model selection, and experimental design in the context of neutron densities. The transfer functions (TFs) are built analytically by expressing the linear response of the objective function (e.g., χ2) to small perturbations of the data. Using the TF formalism, we are able to analyze the expected overall uncertainty—quantified in terms of bias and variance—of the mean square radius and interior density of 48 Ca and 208 Pb. Using relativistic mean field models as a proxy for the weak-charge density—and assuming that a total of five measurements could be performed on the weak form factor of 48 Ca and 208 Pb—we identify the optimal models and experimental locations that minimize the uncertainty in the extraction of the radius and interior density. We also explore the use of the TF formalism to understand the influence of prior distributions for the model parameters, as well as the optimization of model hyperparameters not constrained by the data. Here,we establish how the choice of experimental locations and the model that is used can have a significant impact on the final uncertainties of the extracted quantities of interest. For challenging experiments such as CREX and PREX, a proper quantification of such uncertainties is critical. We have demonstrated how the TF formalism provides several advantages for this type of analysis.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Local modeling for FRF estimation with noisy input measurements

The frequency response function (FRF) is an essential means by which dynamic systems are qualified. In recent years, local modeling approaches have been extensively researched and shown to significantly outperform traditional FRF estimators. However, the standard local modeling approach assumes a perfectly-known system input, which results in biased FRF estimates in the presence of input noise. This paper derives a simple adjustment that can be used to improve FRF estimation for systems subjected to random excitation with noisy input data. This improvement can be implemented with little modification to standard local modeling algorithms and with little additional computational burden. The adjustment is coupled with a model selection procedure to avoid underfitting and overfitting. In conclusion, the methods presented in this paper are validated on a simulation, and they are shown to reduce bias due to input noise.

47 OTHER INSTRUMENTATION↗
DOE OSTI2020

Progress Toward Simulating Departure from Nucleate Boiling at High-Pressure Applications with Selected Wall Boiling Closures

Recently, a Eulerian-based two-fluid computational fluid dynamics (CFD) framework with a wall heat flux partitioning approach has been intensively investigated for departure from nucleate boiling (DNB) simulation under the U.S. Department of Energy–funded Consortium for Advanced Simulation of Light Water Reactors (CASL) program. Understanding of the DNB characteristics over a range of pressurized water reactor–like operating conditions and accurate prediction of boiling crisis in the nuclear power system have been grand challenges because of the large impact of DNB on reactor safety and operational economics. The ultimate goal of this task in the CASL program is to introduce a robust multiphase CFD–based DNB modeling framework that is capable of characterizing an entire boiling history in which the wall boiling mode experiences the following through multiple stages of heat transfer mode: (1) single-phase convective heat transfer, (2) nucleate boiling heat transfer, and (3) identification of the departure of nucleate boiling. To validate the CASL boiling model, we have benchmarked simulated DNB over three different flow channel configurations (pipe flow, 5 × 5 fuel bundle with mixing vane tests, and 5 × 5 fuel bundle without mixing vane tests) against experimental measurements, and the validation result with open literature is reported. The DNB detection criteria in the simulation are checked by monitoring the peak wall temperature, wall dryout factor, and net energy balance. In addition to the DNB performance test, some preliminary sensitivity results on closure model selection are reported to address the prediction capability of local void profile against measurements. The boiling simulation tested in this study exhibits a maximum deviation of 24% from the measured DNB value in a high-pressure (i.e., 138 bars) subcooled pipe flow test. The ranges of operating conditions are as follows: 1650 to 2650 kg/m 2 ·s for mass flux and 8.5 to 96 K for subcooled inlet temperature. The deviation is even reduced to 7% when the subcooled temperature is less than 40 K. Besides accuracy, base practice guidelines for DNB detection criteria are tested by monitoring three simulation variables: (1) maximum wall temperature, (2) wall dryout factor (i.e., K-value), and (3) energy balance. Numerical robustness of DNB simulation is largely achieved in most of the validation test except for a few high subcooled test cases.

42 ENGINEERING↗
DOE OSTI2021

Improved Treatment of Host-galaxy Correlations in Cosmological Analyses with Type Ia Supernovae

Improving the use of Type Ia supernovae (SNe Ia) as standard candles requires a better approach to incorporate the relationship between SNe Ia and the properties of their host galaxies. Using a spectroscopically confirmed sample of ~1600 SNe Ia, we develop the first empirical model of underlying populations for SNe Ia light-curve properties that includes their dependence on host-galaxy stellar mass; we find a significant correlation between stretch population and stellar mass (99.9% confidence) and a weaker correlation between color and stellar mass (90% confidence). These populations are important inputs to simulations that are used to model selection effects and correct distance biases within the BEAMS with Bias Correction (BBC) framework. Here we improve BBC to also account for SNe Ia-host correlations, and we validate this technique on simulated data samples. Here, we recover the input relationship between SNe Ia luminosity and host-galaxy stellar mass (the mass step, γ) with a bias of 0.004 ±0.001 mag, which is a factor of 5 improvement over previous methods that have a γ bias of ~0.02 ± 0.001 mag. We adapt BBC for a novel dust-based model of intrinsic brightness variations, which results in a greatly reduced mass step for data (γ = 0.017 ± 0.008) and for simulations (γ = 0.006 ± 0.007). Analyzing simulated SNe Ia, the biases on the dark energy equation of state, w, vary from Δw = 0.006(5) to 0.010(5) with our new BBC method; these biases are significantly smaller than the 0.02(5) w bias using previous BBC methods that ignore SNe Ia-host correlations.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Elucidating Biomass-Derived Pyrolytic Lignin Structures from Demethylation Reactions through Density Functional Theory Calculations

Pyrolytic lignin is a fraction of pyrolysis oil that contains a wide range of phenolic compounds that can be used as intermediates to produce fuels and chemicals. However, the characteristics of the raw lignin structure make it difficult to establish a pyrolysis mechanism and determine pyrolytic lignin structures. Herein this study proposes dimer, trimer, and tetramer structures based on their relative thermodynamic stability for a hardwood lignin model in pyrolysis. Different configurations of oligomers were evaluated by varying the positions of the guaiacyl (G) and syringyl (S) units and the bonds βO4 and β5 in the hardwood model lignin through electronic structure calculations. The homolytic cleavage of βO4 bonds is assumed to occur and generate two free radical fragments. These can stabilize by taking hydrogen radicals that may be in solution during the intermediate liquid (pathway 1) formation before the thermal ejection. An alternative pathway (pathway 2) could occur when the radicals use intramolecular hydrogen, turning themselves into stable products. Subsequently, a demethylation reaction can take place, thus generating a methane molecule and new oligomeric lignin-derived molecules. The most probable resulting structures were studied. We used FTIR and NMR spectra of selected model compounds to evaluate our calculation approach. Thermophysical properties were calculated using group contribution methods. The results give insights into the lignin oligomer structures and how these molecules are formed. They also provide helpful information for the design of pyrolysis oil separation and upgrading equipment.

09 BIOMASS FUELS↗
DOE OSTI2023

Chaconne: A Statistical Approach to Nonlocal Compression for Supervised Learning, Semi-Supervised Learning, and Anomaly Detection

This project developed a novel statistical understanding of compression analytics (CA), which has challenged and clarified some core assumptions about CA, and enabled the development of novel techniques that address vital challenges of national security. Specifically, this project has yielded the development of novel capabilities including 1. Principled metrics for model selection in CA, 2. Techniques for deriving/applying optimal classification rules and decision theory to supervised CA, including how to properly handle class imbalance and differing costs of misclassification, 3. Two techniques for handling nonlocal information in CA, 4. A novel technique for unsupervised CA that is agnostic with regard to the underlying compression algorithm, 5. A framework for semisupervised CA when a small number of labels are known in an otherwise large unlabeled dataset. 6. The academic alliance component of this project has focused on the development of a novel exemplar-based Bayesian technique for estimating variable length Markov models (closely related to PPM [prediction by partial matching] compression techniques). We have developed examples illustrating the application of our work to text, video, genetic sequences, and unstructured cybersecurity log files.

99 GENERAL AND MISCELLANEOUS↗
DOE OSTI2026

Scaling open-weight large language models for hydropower regulatory information extraction: A systematic analysis

Information extraction from regulatory and technical documents using large language models (LLMs) involves practical trade-offs between extraction quality and computational cost. We evaluate eight open-weight LLMs spanning 0.6B–70B parameters on hydropower licensing documents and report deployment-oriented evidence under a unified extraction schema and evaluation protocol. Across the model set, we observe clear scale-dependent trends in both baseline extraction quality and the effectiveness of reflective reasoning (self-checking) under our fixed-prompt, no-augmentation setting. Mid-scale models often provide a favorable balance of accuracy and efficiency, whereas the smallest models show limited or inconsistent gains from the reasoning variants tested. Larger models achieve the highest overall F1 scores but incur substantially greater compute and infrastructure requirements. We further find that reliability failure modes can distort conventional metrics in this domain: in particular, high recall can coincide with systematic extraction errors when models fabricate values for fields that are absent from the source text, underscoring the importance of conservative null handling and evidence-grounded evaluation. Overall, our study provides a reproducible resource–performance comparison for open-weight LLM-based extraction in hydropower regulatory documentation and offers practical guidance for model selection under different deployment constraints.

Evaluation protocol↗
DOE OSTI2022

Atacama Cosmology Telescope: Constraints on prerecombination early dark energy

The early dark energy (EDE) scenario aims to increase the value of the Hubble constant (H0) inferred from cosmic microwave background (CMB) data over that found in the standard cosmological model (Λ CDM ), via the introduction of a new form of energy density in the early Universe. The EDE component briefly accelerates cosmic expansion just prior to recombination, which reduces the physical size of the sound horizon imprinted in the CMB. Previous work has found that nonzero EDE is not preferred by Planck CMB power spectrum data alone, which yield a 95% confidence level (C.L.) upper limit f EDE < 0.087 on the maximal fractional contribution of the EDE field to the cosmic energy budget. In this paper, we fit the EDE model to CMB data from the Atacama Cosmology Telescope (ACT) data release 4. We find that a combination of ACT, large-scale Planck TT (similar to WMAP), Planck CMB lensing, and BAO data prefers the existence of EDE at >99.7 % C .L . : f EDE = 0.091$^{+0.020}_{-0.036}$, with H 0 = 70.9$^{+1.0}_{-2.0}$ km/s/Mpc (both 68% C.L.). From a model-selection standpoint, we find that EDE is favored over Λ CDM by these data at roughly 3 σ significance. In contrast, a joint analysis of the full Planck and ACT data yields no evidence for EDE, as previously found for Planck alone. We show that the preference for EDE in ACT alone is driven by its TE and EE power spectrum data. The tight constraint on EDE from Planck alone is driven by its high-ℓ TT power spectrum data. Understanding whether these differing constraints are physical in nature, due to systematics, or simply a rare statistical fluctuation is of high priority. Here, the best-fit EDE models to ACT and Planck exhibit coherent differences across a wide range of multipoles in TE and EE, indicating that a powerful test of this scenario is anticipated with near-future data from ACT and other ground-based experiments.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Explosive Yield Estimation Using Regional Seismic Moment Tensors

Here, we use the Pasyanos and Chiang (2022) data set to calculate the seismic moment M 0 for each explosion and use the measured explosive yield W to validate the W~M 0 relationship in Denny and Johnson (1991; hereafter, DJ91). The M 0 is corrected by transforming to a potency tensor and applying more appropriate near-source geophysical parameter values in the moment estimate. The mean residual between observed and predicted yield is near zero; however, the standard deviation of the residuals results in an F-value (a 95% confidence factor) of about 5. We re-estimate the coefficients in the DJ91 model and find similar values and only a slight improvement in the F-value. Next, we embark on a similar model selection process as DJ91, allowing for non-cube-root yield scaling and other plausible near-source elastic moduli. As was found by DJ91, the yield dependence is not significantly different from unity, and a cube root assumption is valid. Therefore, we yield scale the seismic moment and test the significance of all plausible explanatory variables. Isotropic moment performs better in the response variable than total moment. The preference for isotropic moment could be due to its relationship to volume change, which would be more directly affected by explosive yield. Surprisingly, we find that the overburden pressure, which is a function of depth, is not a significant parameter in the model. We hypothesize that this is due to the competing depth effects on source asymmetry and the incorporation of depth in the Green’s functions used to calculate the seismic moment tensors. Importantly, this emphasizes that only seismic moment tensor-derived moments should be used in these models. After removing insignificant model parameters, we are left with a simple model to predict explosive yield $\widehat{W}$ in kt from isotropic moment M I in N·m, $\widehat{W}$=κ –1.4132 10 0.035626GP M I , in which κ and GP are the near-source bulk modulus and gas porosity in Pa and %, respectively. The F-value for this model is approximately 3.

58 GEOSCIENCES↗