SEARCH · Search NASA

Results for “Model fitting”

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 1,045 records · Page 58

All sources NASA NTRS DOE OSTI

DOE OSTI2012

EVOLUTION IN THE VOLUMETRIC TYPE Ia SUPERNOVA RATE FROM THE SUPERNOVA LEGACY SURVEY

We present a measurement of the volumetric Type Ia supernova (SN Ia) rate (SNR{sub Ia}) as a function of redshift for the first four years of data from the Canada-France-Hawaii Telescope Supernova Legacy Survey (SNLS). This analysis includes 286 spectroscopically confirmed and more than 400 additional photometrically identified SNe Ia within the redshift range 0.1 {<=} z {<=} 1.1. The volumetric SNR{sub Ia} evolution is consistent with a rise to z {approx} 1.0 that follows a power law of the form (1+z){sup {alpha}}, with {alpha} = 2.11 {+-} 0.28. This evolutionary trend in the SNLS rates is slightly shallower than that of the cosmic star formation history (SFH) over the same redshift range. We combine the SNLS rate measurements with those from other surveys that complement the SNLS redshift range, and fit various simple SN Ia delay-time distribution (DTD) models to the combined data. A simple power-law model for the DTD (i.e., {proportional_to}t{sup -{beta}}) yields values from {beta} = 0.98 {+-} 0.05 to {beta} = 1.15 {+-} 0.08 depending on the parameterization of the cosmic SFH. A two-component model, where SNR{sub Ia} is dependent on stellar mass (M{sub stellar}) and star formation rate (SFR) as SNR{sub Ia}(z) = A Multiplication-Sign M{sub stellar}(z) + B Multiplication-Sign SFR(z), yields the coefficients A = (1.9 {+-} 0.1) Multiplication-Sign 10{sup -14} SNe yr{sup -1} M{sup -1}{sub Sun} and B = (3.3 {+-} 0.2) Multiplication-Sign 10{sup -4} SNe yr{sup -1} (M{sub Sun} yr{sup -1}){sup -1}. More general two-component models also fit the data well, but single Gaussian or exponential DTDs provide significantly poorer matches. Finally, we split the SNLS sample into two populations by the light-curve width (stretch), and show that the general behavior in the rates of faster-declining SNe Ia (0.8 {<=} s < 1.0) is similar, within our measurement errors, to that of the slower objects (1.0 {<=} s < 1.3) out to z {approx} 0.8.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2012

THERMAL PROPERTIES OF A SOLAR CORONAL CAVITY OBSERVED WITH THE X-RAY TELESCOPE ON HINODE

Coronal cavities are voids in coronal emission often observed above high latitude filament channels. Sometimes, these cavities have areas of bright X-ray emission in their centers. In this study, we use data from the X-ray Telescope (XRT) on the Hinode satellite to examine the thermal emission properties of a cavity observed during 2008 July that contains bright X-ray emission in its center. Using ratios of XRT filters, we find evidence for elevated temperatures in the cavity center. The area of elevated temperature evolves from a ring-shaped structure at the beginning of the observation, to an elongated structure two days later, finally appearing as a compact round source four days after the initial observation. We use a morphological model to fit the cavity emission, and find that a uniform structure running through the cavity does not fit the observations well. Instead, the observations are reproduced by modeling several short cylindrical cavity 'cores' with different parameters on different days. These changing core parameters may be due to some observed activity heating different parts of the cavity core at different times. We find that core temperatures of 1.75 MK, 1.7 MK, and 2.0 MK (for July 19, July 21, and July 23, respectively) in the model lead to structures that are consistent with the data, and that line-of-sight effects serve to lower the effective temperature derived from the filter ratio.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2021

Surrogate modelling the Baryonic Universe II: On forward modelling the colours of individual and populations of galaxies

ABSTRACT Among the properties shaping the light of a galaxy, the star formation history (SFH) is one of the most challenging to model due to the variety of correlated physical processes regulating star formation. In this work, we leverage the stellar population synthesis model fsps, together with SFHs predicted by the hydrodynamical simulation IllustrisTNG and the empirical model universemachine, to study the impact of star formation variability on galaxy colours. We start by introducing a model-independent metric to quantify the burstiness of a galaxy formation model, and we use this metric to demonstrate that universemachine predicts SFHs with more burstiness relative to IllustrisTNG. Using this metric and principal component analysis, we construct families of SFH models with adjustable variability, and we show that the precision of broad-band optical and near-infrared colours degrades as the level of unresolved short-term variability increases. We use the same technique to demonstrate that variability in metallicity and dust attenuation presents a practically negligible impact on colours relative to star formation variability. We additionally provide a model-independent fitting function capturing how the level of unresolved star formation variability translates into imprecision in predictions for galaxy colours; our fitting function can be used to determine the minimal SFH model that reproduces colours with some target precision. Finally, we show that modelling the colours of individual galaxies with per cent-level precision demands resorting to complex SFH models, while producing precise colours for galaxy populations can be achieved using models with just a few degrees of freedom.

Chaves-Montero, Jonás (ORCID:0000000295534261)↗

DOE OSTI2017

A nondestructive method to estimate the chlorophyll content of Arabidopsis seedlings

Chlorophyll content decreases in plants under stress conditions, therefore it is used commonly as an indicator of plant health. Arabidopsis thaliana offers a convenient and fast way to test physiological phenotypes of mutations and treatments. But, chlorophyll measurements with conventional solvent extraction are not applicable to Arabidopsis leaves due to their small size, especially when grown on culture dishes. We provide a nondestructive method for chlorophyll measurement whereby the red, green and blue (RGB) values of a color leaf image is used to estimate the chlorophyll content from Arabidopsis leaves. The method accommodates different profiles of digital cameras by incorporating the ColorChecker chart to make the digital negative profiles, to adjust the white balance, and to calibrate the exposure rate differences caused by the environment so that this method is applicable in any environment. We chose an exponential function model to estimate chlorophyll content from the RGB values, and fitted the model parameters with physical measurements of chlorophyll contents. As further proof of utility, this method was used to estimate chlorophyll content of G protein mutants grown on different sugar to nitrogen ratios. Our method is a simple, fast, inexpensive, and nondestructive estimation of chlorophyll content of Arabidopsis seedlings. This method lead to the discovery that G proteins are important in sensing the C/N balance to control chlorophyll content in Arabidopsis.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2016

West Flank Coso FORGE: Magnetotelluric Inversion

The Coso Magnetotelluric (MT) dataset of which the West Flank FORGE MT data is a subset, was collected by Schlumberger / WesternGeco and initially processed by the WesternGeco GeoSolutions Integrated EM Center of Excellence in Milan, Italy. The 2011 data was based on 99 soundings that were centered on the West Flank geothermal prospect. The new soundings along with previous data from 2003 and 2006 were incorporated into a 3D inversion. Full impedance tensor data were inverted in the 1-3000 Hz range. The modelling report notes several noise sources, specifically the DC powerline that is 20,000 feet west of the survey area, and may have affected data in the 0.02 to 10 Hz range. Model cell dimensions of 450 x 450 x 65 feet were used to avoid computational instability in the 3D model. The fit between calculated and observed MT values for the final model run had an RMS value of 1.807. The included figure from the WesternGeco report shows the sounding locations from the 2011, 2006 and 2003 surveys.

15 GEOTHERMAL ENERGY↗

DOE OSTI2019

Probing Galaxy assembly bias in BOSS galaxies using void probabilities

ABSTRACT We measure the void probability function (VPF) of galaxies in the Baryon Oscillation Spectroscopic Survey (BOSS). The VPF provides complementary information to standard two-point statistics in that it is sensitive to galaxy bias in the most extreme underdensities in the cosmic web. Thus, the VPF is ideal for testing whether halo occupation of galaxies depends on large-scale density, an effect known as galaxy assembly bias. We find that standard HOD model – one parameterized by halo mass only – fit only to the two-point function, accurately predicts the VPF. Additionally, for HOD models where density dependence is explicitly incorporated, the best-fitting models fit to the combination of the correlation function and the VPF have zero density dependence. Thus, galaxy assembly bias is not a strong source of systematic uncertainty when modelling the clustering of massive galaxies.

Walsh, Kilian↗

DOE OSTI2012

FERMI-LAT OBSERVATIONS OF THE DIFFUSE {gamma}-RAY EMISSION: IMPLICATIONS FOR COSMIC RAYS AND THE INTERSTELLAR MEDIUM

The {gamma}-ray sky >100 MeV is dominated by the diffuse emissions from interactions of cosmic rays with the interstellar gas and radiation fields of the Milky Way. Observations of these diffuse emissions provide a tool to study cosmic-ray origin and propagation, and the interstellar medium. We present measurements from the first 21 months of the Fermi Large Area Telescope (Fermi-LAT) mission and compare with models of the diffuse {gamma}-ray emission generated using the GALPROP code. The models are fitted to cosmic-ray data and incorporate astrophysical input for the distribution of cosmic-ray sources, interstellar gas, and radiation fields. To assess uncertainties associated with the astrophysical input, a grid of models is created by varying within observational limits the distribution of cosmic-ray sources, the size of the cosmic-ray confinement volume (halo), and the distribution of interstellar gas. An all-sky maximum-likelihood fit is used to determine the X{sub CO} factor, the ratio between integrated CO-line intensity and H{sub 2} column density, the fluxes and spectra of the {gamma}-ray point sources from the first Fermi-LAT catalog, and the intensity and spectrum of the isotropic background including residual cosmic rays that were misclassified as {gamma}-rays, all of which have some dependency on the assumed diffuse emission model. The models are compared on the basis of their maximum-likelihood ratios as well as spectra, longitude, and latitude profiles. We also provide residual maps for the data following subtraction of the diffuse emission models. The models are consistent with the data at high and intermediate latitudes but underpredict the data in the inner Galaxy for energies above a few GeV. Possible explanations for this discrepancy are discussed, including the contribution by undetected point-source populations and spectral variations of cosmic rays throughout the Galaxy. In the outer Galaxy, we find that the data prefer models with a flatter distribution of cosmic-ray sources, a larger cosmic-ray halo, or greater gas density than is usually assumed. Our results in the outer Galaxy are consistent with other Fermi-LAT studies of this region that used different analysis methods than employed in this paper.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2018

Ion mass dependence of irradiation-induced damage accumulation in KTaO 3

Damage production and amorphization resulting from the interaction of medium-energy (from 40 to 480 keV) noble-gas ions (from He to Kr) with potassium tantalate (KTaO 3 ) are determined using ion channeling measurements. A disorder accumulation model has been fit to the maximum damage concentration versus ion fluence to extract the cross sections for direct-impact and defect-stimulated amorphization, and the results indicate that defect-stimulated amorphization is the dominant mechanism. These cross sections exhibit a strong dependence on the calculated cross sections for displacing lattice atoms, indicating a dominant contribution of nuclear interactions to the defect production and amorphization processes under the irradiation conditions used in this study. Furthermore, these experimental findings, along with the model fits, suggest that the difference in recoil spectra between He and the other heavier ions may be the main driving force for the decreased damage efficiency observed for He ions, which results in a reduced rate of damage accumulation.

36 MATERIALS SCIENCE↗

DOE OSTI2023

Formal Experimentation and Analysis of Handheld RFID Readers as a Tool for Nuclear Material Accounting

Commercial, off-the-shelf Radio Frequency Identification (RFID) systems have been successfully deployed for inventory tracking in numerous industries, but their viability in the tracking of complex environments containing nuclear material is less understood. Of primary interest in this setting is a RFID tracking system for nuclear material accounting which can reliably identify as many tags as possible with high accuracy, while also reducing an operator’s exposure to radiation. In this work, we develop a formal statistical approach to identify relevant handheld RFID reader settings which optimize tagging performance. To achieve this goal, we design a full factorial split-plot experiment for a static shelf configuration scene with 50 randomly placed nuclear material containers affixed with RFID tags. We use Bayesian inference to fit a second-order response surface model which expresses the probability of a successful match for each container as a function of the experimental factors. Such effects are allowed to vary across individual containers and the containers’ population in its entirety to estimate overall effects. Uncertainties of estimates and predictions are quantified via their corresponding posterior distributions. Following extensive model checking and validation, the fitted model is used to identify experimental factors which maximize matching probabilities at both the container-level and for the full shelf configuration scene. We also analyze sensitivity of performance to relevant factors.

25 ENERGY STORAGE↗

DOE OSTI2016

Exploring the spatial variability of soil properties in an Alfisol Catena

Detailed digital soil maps showing the spatial heterogeneity of soil properties consistent with the landscape are required for site-specific management of plant nutrients, land use planning and process-based environmental modeling. We characterized the short-scale spatial heterogeneity of soil properties in an Alfisol catena in a tropical landscape of Sri Lanka. The impact of different land-uses (paddy, vegetable and un-cultivated) was examined to assess the impact of anthropogenic activities on the variability of soil properties at the catenary level. Conditioned Latin hypercube sampling was used to collect 58 geo-referenced topsoil samples (0–30 cm) from the study area. Soil samples were analyzed for pH, electrical conductivity (EC), organic carbon (OC), cation exchange capacity (CEC) and texture. The spatial correlation between soil properties was analyzed by computing crossvariograms and subsequent fitting of theoretical model. Spatial distribution maps were developed using ordinary kriging. The range of soil properties, pH: 4.3–7.9; EC: 0.01–0.18 dS m –1 ; OC: 0.1–1.37%; CEC: 0.44– 11.51 cmol (+) kg –1 ; clay: 1.5–25% and sand: 59.1–84.4% and their coefficient of variations indicated a large variability in the study area. Electrical conductivity and pH showed a strong spatial correlation which was reflected by the cross-variogram close to the hull of the perfect correlation. Moreover, cross-variograms calculated for EC and Clay, CEC and OC, CEC and clay and CEC and pH indicated weak positive spatial correlation between these properties. Relative nugget effect (RNE) calculated from variograms showed strongly structured spatial variability for pH, EC and sand content (RNE < 25%) while CEC, organic carbon and clay content showed moderately structured spatial variability (25% < RNE < 75%). Spatial dependencies for examined soil properties ranged from 48 to 984 m. The mixed effects model fitting followed by Tukey's post-hoc test showed significant effect of land use on the spatial variability of EC. Our study revealed a structured variability of topsoil properties in the selected tropical Alfisol catena. Except for EC, observed variability was not modified by the land uses. Investigated soil properties showed distinct spatial structures at different scales and magnitudes of strength. Our results will be useful for digital soil mapping, site specific management of soil properties, developing appropriate land use plans and quantifying anthropogenic impacts on the soil system.

Geo-statistical analysis↗

DOE OSTI2017

Evaluating penalized logistic regression models to predict Heat-Related Electric grid stress days

Understanding the conditions associated with stress on the electricity grid is important in the development of contingency plans for maintaining reliability during periods when the grid is stressed. In this paper, heat-related grid stress and the relationship with weather conditions is examined using data from the eastern United States. Penalized logistic regression models were developed and applied to predict stress on the electric grid using weather data. The inclusion of other weather variables, such as precipitation, in addition to temperature improved model performance. Several candidate models and datasets were examined. A penalized logistic regression model fit at the operation-zone level was found to provide predictive value and interpretability. Additionally, the importance of different weather variables observed at different time scales were examined. Maximum temperature and precipitation were identified as important across all zones while the importance of other weather variables was zone specific. The methods presented in this work are extensible to other regions and can be used to aid in planning and development of the electrical grid.

Bramer, L. M.↗

DOE OSTI2012

THE INITIAL-FINAL MASS RELATION AMONG WHITE DWARFS IN WIDE BINARIES

We present the initial-final mass relation derived from 10 white dwarfs in wide binaries that consist of a main-sequence star and a white dwarf. The temperature and gravity of each white dwarf were measured by fitting theoretical model atmospheres to the observed spectrum using a {chi}{sup 2} fitting algorithm. The cooling time and mass were obtained using theoretical cooling tracks. The total age of each binary was estimated from the chromospheric activity of its main-sequence component to an uncertainty of about 0.17 dex in log t. The difference between the total age and white dwarf cooling time is taken as the main-sequence lifetime of each white dwarf. The initial mass of each white dwarf was then determined using stellar evolution tracks with a corresponding metallicity derived from spectra of their main-sequence companions, thus yielding the initial-final mass relation. Most of the initial masses of the white dwarf components are between 1 and 2 M{sub Sun }. Our results suggest a correlation between the metallicity of a white dwarf's progenitor and the amount of post-main-sequence mass loss it experiences-at least among progenitors with masses in the range of 1-2 M{sub Sun }. A comparison of our observations to theoretical models suggests that low-mass stars preferentially lose mass on the red giant branch.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2016

WE-H-BRA-07: Mechanistic Modelling of the Relative Biological Effectiveness of Heavy Charged Particles

Purpose Uncertainty in the Relative Biological Effectiveness (RBE) of heavy charged particles compared to photons remains one of the major uncertainties in particle therapy. As RBEs depend strongly on clinical variables such as tissue type, dose, and radiation quality, more accurate individualised models are needed to fully optimise treatments. MethodsWe have developed a model of DNA damage and repair following X-ray irradiation in a number of settings, incorporating mechanistic descriptions of DNA repair pathways, geometric effects on DNA repair, cell cycle effects and cell death. Our model has previously been shown to accurately predict a range of biological endpoints including chromosome aberrations, mutations, and cell death. This model was combined with nanodosimetric models of individual ion tracks to calculate the additional probability of lethal damage forming within a single track. These lethal damage probabilities can be used to predict survival and RBE for cells irradiated with ions of different Linear Energy Transfer (LET). ResultsBy combining the X-ray response model with nanodosimetry information, predictions of RBE can be made without cell-line specific fitting. The model’s RBE predictions were found to agree well with empirical proton RBE models (Mean absolute difference between models of 1.9% and 1.8% for cells with α/β ratios of 9 and 1.4, respectively, for LETs between 0 and 15 keV/µm). The model also accurately recovers the impact of high-LET carbon ion exposures, showing both the reduced efficacy of ions at extremely high LET, as well as the impact of defects in non-homologous end joining on RBE values in Chinese Hamster Ovary cells.ConclusionOur model is predicts RBE without the inclusion of empirical LET fitting parameters for a range of experimental conditions. This approach has the potential to deliver improved personalisation of particle therapy, with future developments allowing for the calculation of individualised RBEs. SJM is supported by a Marie Curie International Outgoing Fellowship from the European Commission’s FP7 program (EC FP7 MC-IOF-623630)

60 APPLIED LIFE SCIENCES↗

DOE OSTI2018

Luminescence methodology to determine grain-boundary, grain-interior, and surface recombination in thin-film solar cells

We determine the grain-boundary (GB) recombination velocity, S GB , and grain-interior (GI) lifetime, t GI , parameters in superstrate CdS/CdTe thin-film solar cell technology by combining cathodoluminescence (CL) spectrum imaging and time-resolved photoluminescence (TRPL) measurements. We consider critical device formation stages, including after CdTe deposition, CdCl 2 treatment, and Cu diffusion. CL image analysis methods extract GB and GI intensities and grain size for hundreds of grains per sample. Concurrently, a three-dimensional CL model is developed to simulate the GI intensity as a function of t GI , S GB , grain size, and the surface recombination velocity, S surf . TRPL measurements provide an estimate of S surf for the CL model. A fit of GI intensity vs. grain size data with the CL model gives a self-consistent and representative set of S GB and t GI values for the samples: S GB (t GI ) = 2.6 x 10 6 cm/s (68-250 ps), S GB (t GI ) = 4.1 x 10 5 cm/s (1.5-3.3 ns), and S GB (t GI ) = 5.5 x 10 5 cm/s (1.0-3.8 ns) for as-deposited, CdCl 2 -treated, and CdCl 2- A nd Cu-treated samples, respectively. Thus, we find that the CdCl2 treatment both helps to passivate GBs and significantly increase the GI lifetime. Subsequent Cu diffusion increases GB recombination slightly and has nuanced effects on the GI lifetime. Finally, as a partial check on the S GB and t GI values, they are input to a Sentaurus device model, and the simulated performance is compared to the measured performance. The methodology developed here can be applied broadly to CdTe and CdSeTe thin-film technology and to other thin-film solar cell materials including Cu(In 1-x Ga x )Se 2 , Cu 2 ZnSnS 4 , and perovskites.

14 SOLAR ENERGY↗

DOE OSTI2014

CHARA/MIRC observations of two M supergiants in Perseus OB1: Temperature, bayesian modeling, and compressed sensing imaging

Two red supergiants (RSGs) of the Per OB1 association, RS Per and T Per, have been observed in the H band using the Michigan Infra-Red Combiner (MIRC) instrument at the CHARA array. The data show clear evidence of a departure from circular symmetry. We present here new techniques specially developed to analyze such cases, based on state-of-the-art statistical frameworks. The stellar surfaces are first modeled as limb-darkened disks based on SATLAS models that fit both MIRC interferometric data and publicly available spectrophotometric data. Bayesian model selection is then used to determine the most probable number of spots. The effective surface temperatures are also determined and give further support to the recently derived hotter temperature scales of RSGs. The stellar surfaces are reconstructed by our model-independent imaging code SQUEEZE, making use of its novel regularizer based on Compressed Sensing theory. We find excellent agreement between the model-selection results and the reconstructions. Our results provide evidence for the presence of near-infrared spots representing about 3%-5% of the stellar flux.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2025

Revenue-Maximizing Shared Parking and Electric Vehicle Charging Management in Multi-Unit Dwellings

In urban areas, searching for parking and electric vehicle (EV) charging can result in cruising, congestion, and environmental externalities. Recognizing the business opportunity of offering private parking and charging infrastructure access within multi-unit dwellings (MUDs) during daytime, we model a shared parking and EV charging management system. We maximize the revenue of MUD charging hubs in mixed land use, catering to public demand. Our approach accounts for the objectives of the two stakeholders involved: a demand model is fitted on the choices of EV charging users, and the supply model optimizes the allocation of parking and charging requests in an MUD parking lot. A binary integer linear programming model for the allocation of parking and charging spaces with a rolling horizon is integrated with matching rules that handle both parking and charging requests. In our numerical experiments in a neighborhood of Chicago, Illinois, we estimate the performance of the MUD parking and charging system with metrics that include revenue, number of matchings, and utilization rates. At any given time, MUDs with lower prices attract more charging requests, particularly those of longer duration, resulting in higher revenue and greater charging utilization. Dynamic pricing facilitates a more equitable distribution of requests; as MUD parking lots reach capacity and their fees increase, other MUDs become more competitive, attracting additional requests. Comparing our method against first-come-first-served and optimal-solution benchmarks, we demonstrate our model’s effectiveness in dynamically managing mixed parking and charging demand in MUD charging hubs.

electric vehicle, multi-unit dwelling, charging in↗

DOE OSTI2026

Semi-Analytical Hierarchical Bayesian Inference of Nonlinear Model Structure in Stochastic Dynamics: Applied to Compartmental Models of Infectious Diseases

A Bayesian computational framework for parsimonious inference in stochastic nonlinear dynamical systems is presented. This framework enables the concurrent estimation of system states, time-varying parameters, time-invariant parameters, and the optimal sparsity structure of the model parameters. Because differential equation-based models are often simplified mechanistic or phenomenological representations, robust inference from noisy measurement data requires explicit treatment of model error and uncertainty. Model error and time-varying parameters can be represented as random processes, enabling inference while making minimal assumptions about the underlying sources of discrepancy and variability. Adopting stochastic differential equation representations affords the model significant flexibility, but can also render it susceptible to overfitting during statistical inversion, where the inferred model may track noise rather than the underlying signal. To alleviate the effects of overfitting and to enable the discovery of the optimal sparse representation of the time-invariant parameters, a Bayesian sparse learning algorithm is embedded within the framework. This sparse learning framework adopts an approximate hierarchical Bayesian setting defined by a series of semi-analytical expressions. The model structure inference framework is validated using a stochastic compartmental model for tracking and forecasting active cases of an infectious disease. Compartmental models describe population-level infectious disease dynamics through interactions among population fractions grouped by disease state. Mathematically, such models consist of a system of coupled ordinary differential equations. This example adopts an expressive compartmental model that includes multiple possible interactions between disease states, motivated by early uncertainty surrounding COVID-19 reinfection dynamics and their implications for long-term epidemic forecasting. The sparse learning exercise permits the inference of a priori unknown epidemiological dynamics from simulated public health data, discovering the nested compartmental model that optimizes the trade-off between average data-fit and model complexity. It is shown that inducing sparsity among the model parameters eliminates redundant interactions between compartments, equivalently revealing the optimal coupling structure between differential equations.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2018

Explosion Source Models and the Scattering Origin of Regional Phases from SPE Phase 1 Coda Spectral Ratios

The Source Physics Experiment(SPE)explosions provide comprehensive ground truth information,and nearby recordings that are rarely available in typical monitoring scenarios. Here, we look for relationships with local and regional distance data to understand source and propagation effects in an effort to improve monitoring at a distance. We examine source spectral ratios, which can be measured to high precision by taking advantage of many SPE and independently operated stations, as well as the redundancy available from coda waves. Classical source models (Mueller and Murphy, 1971, MM71; Denny and Johnson, 1991, DJ91)predict the source ratios poorly; however, a hybrid model(MM71 with DJ91 cavity radii)performs better, and by varying cavity and elastic radii we obtain models that fit to within measurement precision. Based on this, the best prediction model for the SPE suite is the hybrid model with cavity radius decreased by 95%, and elastic radius decreased by 75%. This works best for mid-range scaled depths of burial (SPE-1, 2, 3, and 5). We do not guarantee transportability to other geologies,depth and scaled depth ranges. Finally, we observe a distinct spectral modulation at 6-9Hzthat is not predicted by classical models, and is most likely the result of interference of short period surface waves(Rg)with those produced by spallation of near surface layers.The same modulationis observed for compressional (P) and shear(S) waves at distance, indicating the origin of these waves as scattered Rg.The scattering process must be accounted for quantitatively when using distant signals to monitor treaty compliance over broad regions.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗