Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “model interpretability”

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 397 records · Page 22
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Stable isotopes contain substantial additive information about terrestrial carbon and water cycling

Abstract Stable isotope ratios of H ( δ 2 H ), O ( δ 18 O), and C ( δ 13 C) are linked to key biogeochemical processes of the water and carbon cycles; however, the degree to which isotope-associated processes are reflected in macroscale ecosystem flux observations remains unquantified. Here through formal information assessment, new measurements of δ 13 C of net ecosystem exchange ( NEE ) as well as δ 2 H and δ 18 O of latent heat ( LH ) fluxes across the United States National Ecological Observation Network (NEON) are used to determine conditions under which isotope measurements are informative of environmental exchanges. We find all three isotopic datasets individually contain comparable amounts of information about NEE and LH fluxes as wind speed observations. Such information from isotope measurements, however, is largely unique. Generally, δ 13 C provides more information about LH as aridity increases or mean annual precipitation decreases. δ 2 H provides more information about LH as temperatures or mean annual precipitation decreases, and also provides more information about NEE as temperatures decrease. Overall, we show that the stable isotope datasets collected by NEON contribute non-trivial amounts of new information about bulk environmental fluxes useful for interpreting biogeochemical and ecohydrological processes at landscape scales. However, the utility of this new information varies with environmental conditions at continental scales. This study provides an approach for quantifying the value adding non-traditional sensing approaches to environmental monitoring sites and the patterns identified here are expected to aid in modeling and data interpretation efforts focused on constraining carbon and water cycles’ mechanisms.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Neural networks for estimation of divertor conditions in DIII-D using C III imaging

Deep learning approaches have been applied to images of C III emission in the lower divertor of DIII-D to develop models for estimating the level of detachment and magnetic configuration (X-point location and strike point radial location). The poloidal distance from the target to the C III emission front is used to represent the level of detachment. The models perform well on a test dataset not used in training, achieving $F_1$ scores as high as 0.99 for detachment state classification and root mean squared error (RMSE) as low as 2cm for front location regression. Predictions for shots with intermittent reattachment are studied, with class activation mapping used to aid in interpretation of the model predictions. Based on the success of these models, a third model was trained to predict the X-point location and strike point radial position from C III images. Though the dataset covers only a small range of possible magnetic configurations, the model shows promising results, achieving RMSE around 1cm for the test data.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Scientific machine learning for closure models in multiscale problems: A review

Here, closure problems are omnipresent when simulating multiscale systems, where some quantities and processes cannot be fully prescribed despite their effects on the simulation's accuracy. Recently, scientific machine learning approaches have been proposed as a way to tackle the closure problem, combining traditional (physics-based) modeling with data-driven (machine-learned) techniques, typically through enriching differential equations with neural networks. This paper reviews the different reduced model forms, distinguished by the degree to which they include known physics, and the different objectives of a priori and a posteriori learning. The importance of adhering to physical laws (such as symmetries and conservation laws) in choosing the reduced model form and choosing the learning method is discussed. The effect of spatial and temporal discretization and recent trends toward discretization-invariant models are reviewed. In addition, we make the connections between closure problems and several other research disciplines: inverse problems, Mori-Zwanzig theory, and multi-fidelity methods. In conclusion, much progress has been made with scientific machine learning approaches for solving closure problems, but many challenges remain. In particular, the generalizability and interpretability of learned models is a major issue that needs to be addressed further.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Modeling Electron Acceleration and Transport in the Early Impulsive Phase of the 2017 September 10th Solar Flare

The X8.2-class limb flare on 2017 September 10 is among the best studied solar flare events owing to its great similarity to the standard flare model and the broad coverage by multiple spacecraft and ground-based observations. These multiwavelength observations indicate that electron acceleration and transport are efficient in the reconnection and flare looptop regions. However, there lacks a comprehensive model for explaining and interpreting the multi-faceted observations. In this work, we model the electron acceleration and transport in the early impulsive phase of this flare. We solve the Parker transport equation that includes the primary acceleration mechanism during magnetic reconnection in the large-scale flare region modeled by MHD simulations. We find that electrons are accelerated up to several MeV and fill a large volume of the reconnection region, similar to the observations shown in microwaves. The electron spatial distribution and spectral shape in the looptop region agree well with those derived from the microwave and hard X-ray emissions before magnetic islands grow large and dominate the acceleration. Future emission modelings using the electron maps will enable direct comparison with microwave and hard X-ray observations. These results shed new light on the electron acceleration and transport in a broad region of solar flares within a data-constrained realistic flare geometry.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Multi-Omics Reveals Temporal Scales of Carbon Metabolism in Synechococcus Elongatus PCC 7942 Under Light Disturbance

Central carbon metabolism in model cyanobacteria involves multiple pathways to adapt to energy-light limitations across diel cycles. However, the success in mechanistic modeling for phenotypic prediction of the protein regulators in the metabolic state depends on capturing the vast possibilities emerging from multiple regulatory pathways in complex biological processes. Here, we developed a physics-informed machine learning approach based on energy-landscape concepts to predict regulatory proteins responding to cyclic circadian and unforeseen light perturbations in cyanobacterial metabolic networks. Our approach provides interpretable de novo models for inferring gene expression dynamics from Synechococcus elongatus over diel cycles and using redox proteome analysis to distinguish immediate light-responsive elements from circadian-regulated processes in carbon metabolism pathways. We identified distinct temporal signatures with the analysis of the redox proteome: there was an immediate shift in cysteine redox states accompanied by a limited change in protein abundance under constant illumination and after 2 hours of darkness. This discovery indicates that the generation of reductants coordinates photoinduced electron transport with redox metabolic pathways in two discernable molecular mechanisms: fast redox-based protein modifications occur immediately after the light disturbance, followed by slow transcriptional regulations across networks. This temporal regulation reveals how metabolic networks integrate rapid light responses with programmed circadian rhythms to maintain cellular homeostasis under the light-energy limitations over the diel cycle.

Biomolecular & subcellular processes↗
DOE OSTI2024

Interpretable Net Load Forecasting Using Smooth Multiperiodic Features

We consider the problem of forecasting net load over a horizon such as one day, using a trailing window of past net load values as well as date and time. We focus on three variations on this problem: point forecasts, marginal quantile forecasts, and generating conditional samples of the future value. We propose a method that relies on linear regression using some custom engineered time-based features to capture multiple periodicities, such as daily, weekly, and seasonal, and their interactions. Our proposed models are readily interpretable, and rely on efficient and reliable convex optimization [1] to fit. We illustrate our method on four years worth of hourly net load data, comparing predictions made with various subsets of the features.

Ogut, Mehmet G↗
DOE OSTI2020

Prediction of Hydroxymethylfurfural Yield in Glucose Conversion through Investigation of Lewis Acid and Organic Solvent Effects

Hydroxymethylfurfural (HMF) is one of the important renewable platform compounds that can be obtained from biomass feedstocks through glucose conversion catalyzed by Brønsted and Lewis acids. However, it is challenging to enhance the HMF yield due to side reactions. In this study, a systematic approach combining theory and experiment was performed to investigate the influence of Lewis acids and organic solvents on the HMF yield. For the Lewis acid effect, a relationship between chemical hardness and experimental HMF yields was found in the rate-limiting step of glucose-to-fructose isomerization for six metal chlorides; HMF production was promoted when the metal chloride and a substrate had a similar chemical hardness. To study the organic solvent effect, a multivariate model was developed based on the insights gained from the mechanistic study of fructose dehydration, to predict HMF yields in a given water-organic cosolvent system. It showed a reliable accuracy in evaluating HMF yields with a mean absolute error (MAE) of 3.0% with respect to experimental HMF yields for 13 solvents, and also predicted HMF yields with a MAE of 10.7% for four new solvents. Chemical interpretation of the model revealed that it is desirable to use a solvent capable of stabilizing the carbocation intermediates with low proton transfer activity and high hydrogen bond basicity, to maximize the HMF yield. This multivariate model informs experimentalists about rational selection of solvents with very low computational costs needed to calculate only six variables for each solvent. It can be expanded to other catalytic systems such as heterogeneous Brønsted–Lewis bifunctional catalysts and enables optimization of reaction conditions to obtain other useful platform molecules through biomass conversion.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2023

Carbon isotope fractionation by an ancestral rubisco suggests that biological proxies for CO 2 through geologic time should be reevaluated

The history of Earth’s carbon cycle reflects trends in atmospheric composition convolved with the evolution of photosynthesis. Fortunately, key parts of the carbon cycle have been recorded in the carbon isotope ratios of sedimentary rocks. The dominant model used to interpret this record as a proxy for ancient atmospheric CO 2 is based on carbon isotope fractionations of modern photoautotrophs, and longstanding questions remain about how their evolution might have impacted the record. Therefore, we measured both biomass (ε p ) and enzymatic (ε Rubisco ) carbon isotope fractionations of a cyanobacterial strain ( Synechococcus elongatus PCC 7942) solely expressing a putative ancestral Form 1B rubisco dating to ≫1 Ga. This strain, nicknamed ANC, grows in ambient pCO 2 and displays larger ε p values than WT, despite having a much smaller ε Rubisco (17.23 ± 0.61‰ vs. 25.18 ± 0.31‰, respectively). Surprisingly, ANC ε p exceeded ANC ε Rubisco in all conditions tested, contradicting prevailing models of cyanobacterial carbon isotope fractionation. Such models can be rectified by introducing additional isotopic fractionation associated with powered inorganic carbon uptake mechanisms present in Cyanobacteria, but this amendment hinders the ability to accurately estimate historical pCO 2 from geological data. Understanding the evolution of rubisco and the CO 2 concentrating mechanism is therefore critical for interpreting the carbon isotope record, and fluctuations in the record may reflect the evolving efficiency of carbon fixing metabolisms in addition to changes in atmospheric CO 2 .

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2025

Interpretation of Ion Irradiation and Neutron Irradiation Damage in Additively Manufactured 316 Stainless Steel using Multiscale Modeling

The accelerated adoption of nuclear energy necessitates advanced manufacturing technologies, such as additive manufacturing, to meet heightened supply chain requirements and support innovative reactor technologies. Due to the unique microstructural characteristics of additively manufactured materials under distinct solidification conditions, comprehensive evaluation of their performance in reactor environments is essential. The Advanced Materials and Manufacturing Technologies program under the Department of Energy's Office of Nuclear Energy focuses on understanding the irradiation performance and damage evolution of laser powder bed fusion 316 stainless steel, with an emphasis on integrating ion and neutron irradiation data to accelerate the development and qualification of materials for advanced nuclear reactor applications. While ion irradiation is a cost- and time-effective method, modeling and simulation are required to interpret the data for the broader range of irradiation conditions encountered in advanced reactors. In fiscal year 2025, integrated multiscale modeling and simulations were conducted to assess irradiation damage in additively manufactured 316 stainless steel. Key outcomes include predictions of chromium enrichment at grain boundaries, nickel enrichment at dislocation cell walls and void surfaces, and heterogeneous void evolution under ion and neutron irradiation conditions. Cluster dynamics simulations revealed the coarsening of voids at high irradiation temperatures and the suppression of void growth by high network dislocation density, while also demonstrating significant growth and coarsening of voids and self-interstitial atom loops at low dose rates. Machine learning-accelerated atomistic simulations highlighted the impact of the local environment and chromium concentration on vacancy diffusivity, providing key insights on the influence of composition on void swelling and radiation-induced segregation. Additionally, molecular dynamics simulations demonstrated the presence of defect production bias and a significant effect of carbon content on defect cluster behavior. These combined efforts aim to predict the performance of additively manufactured materials under various reactor conditions, supporting their qualification for nuclear reactor applications by interpreting ion irradiation data. This report underscores the potential of integrated multiscale modeling to analyze ion irradiation data in the effort to accelerate the qualification of additively manufactured materials for nuclear reactor components.

316 stainless steel↗
DOE OSTI2025

Charge Radii of Neutron-Rich Scandium Isotopes and the Seniority Symmetry in the 0⁢𝑓 7/2 Shell

Nuclear charge radii of neutron-rich 47–49 Sc isotopes were measured using collinear laser spectroscopy at CERN-ISOLDE. The new data reveal that the charge radii of scandium isotopes exhibit a distinct trend between 𝑁 = 20 and 𝑁 = 28, with 41 Sc and 49 Sc isotopes having similar values, mirroring the closeness of the charge radii of 40 Ca and 48 Ca. Theoretical models that successfully interpret the radii of calcium isotopes could not account for the observed behavior in scandium radii, in particular the reduced odd-even staggering. Remarkably, the inclusion of the new 49 Sc radius data has unveiled a similar trend in the charge radii of 𝑁 = 28 isotones and 𝑍 = 20 isotopes when adding the protons atop the 48 Ca core and the neutrons atop the 40 Ca core, respectively. We demonstrate that this trend is consistent with the prediction of the seniority model.

Nuclear charge distribution↗
DOE OSTI2026

Sim-to-real supervised domain adaptation for radioisotope identification

Machine learning has the potential to improve the speed and reliability of radioisotope identification using gamma spectroscopy. However, meticulously labeling an experimental dataset for training is often prohibitively expensive, while training models purely on synthetic data is risky due to the domain gap between simulated and experimental measurements. In this research, we demonstrate that supervised domain adaptation can substantially improve the performance of radioisotope identification models by transferring knowledge between synthetic and experimental data domains. We consider two domain adaptation scenarios: (1) a simulation-to-simulation adaptation, where we perform multi-label proportion estimation using simulated high-purity germanium detectors, and (2) a simulation-to-experimental adaptation, where we perform multi-class, single-label classification using measured spectra from handheld lanthanum bromide (LaBr) and sodium iodide (NaI) detectors. We begin by pretraining a spectral classifier on synthetic data using a custom transformer-based neural network. After subsequent fine-tuning on just 64 labeled experimental spectra, we achieve a test accuracy of 96% in the sim-to-real scenario with a LaBr detector, far surpassing a synthetic-only baseline model (75%) and a model trained from scratch (80%) on the same 64 spectra. Furthermore, we demonstrate that domain-adapted models learn more human-interpretable features than experiment-only baseline models. Overall, our results highlight the potential for supervised domain adaptation techniques to bridge the sim-to-real gap in radioisotope identification, enabling the development of accurate and explainable classifiers even in real-world scenarios where access to experimental data is limited.

Lalor, Peter W.↗
DOE OSTI2025

Decoding substrate specificity determining factors in glycosyltransferase-B enzymes – insights from machine learning models

Substrate specificity is an essential characteristic of any enzyme's function and an understanding of the factors that determine this specificity is crucial for enzyme engineering. Unlike the structure of an enzyme which is directly impacted by its sequence, substrate specificity as an enzyme attribute involves a rather indirect relationship with sequence as it also depends on structural aspects that dictate substrate accessibility and active site dynamics. In this study, we explore the performance of classifier-based machine learning models trained on curated sequence and structural data for a class of glycosyltransferases (GTs), namely GT-Bs, to understand their substrate specificity determining factors. GTs enable the transfer of sugar moieties to other biomolecules such as oligosaccharides or proteins and are found in all kingdoms of life. In plants, GTs participate in the biosynthesis of plant cell wall biopolymers (e.g.: hemicelluloses and pectins) and are an integral part of the enzymatic machinery that enables the storage of carbon and energy as plant biomass. To elucidate the substrate specificity of uncharacterized GT-Bs, we constructed multi-label machine learning models (Support Vector Classifier, K-Nearest Neighbors, Gaussian Naïve-Bayes, Random Forest) that incorporate both sequence and structural features. These models achieve good predictive accuracies on test datasets. However, despite our use of structural information, we highlight that there is further scope for improvement in training these models to draw interpretable relationships between sequence, structure and substrate specificity determining motifs in GT-Bs.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2026

Attention-based explainability for structure–property relationships

Machine learning methods are emerging as a universal paradigm for constructing correlative structure–property relationships in materials science based on multimodal characterization. However, this necessitates the development of methods for the physical interpretability of the resulting correlative models. Here, we demonstrate the potential of attention-based neural networks for revealing structure–property relationships and the underlying physical mechanisms, using the ferroelectric properties of PbTiO3 thin films as a case study. Through the analysis of attention scores, we disentangle the influence of distinct domain patterns on the polarization switching process. The attention-based Transformer model is explored both as a direct interpretability tool and as a surrogate for explaining representations learned via unsupervised machine learning, enabling the identification of physically grounded correlations. We compare attention-derived interpretability scores with classical SHapley Additive exPlanations analysis and show that, in contrast to applications in natural language processing, attention mechanisms in materials science exhibit high efficiency in highlighting meaningful structural features.

Slautin, Boris [Independent Researcher]↗
DOE OSTI2020

Ordovician-Cambrian Units: Hierarchical Evaluation of Geologic Carbon Storage Resource Estimates

The Indiana Geological and Water Survey (IGWS) led subtask 1.1 to assess the regional distribution and estimate the storage capacity of Ordovician-Cambrian stratigraphic units located within the partnership region. A comprehensive data set of wireline logs and petrophysical information was used to generate these interpretations. These data include core analysis for porosity and permeability, mercury injection capillary pressure (MICP), and existing well data including location and stratigraphic information. This report includes storage resource estimates (SREs) for three potential storage reservoirs(limestone and dolostone from the Upper Ordovician Trenton Limestone/Black River Group and equivalent units, the Middle Ordovician St. Peter Sandstone, and primary target reservoir rocks of the Lower Ordovician and Upper Cambrian Knox Supergroup and equivalent units) calculated using six methodologies: (1) a fixed value of porosity of 10 percent in all units evaluated; (2) a unique average porosity (per well) from wireline-derived porosity (neutron, sonic, and/or density porosity for each unit); (3) porosity values from core analysis; (4) a depth-dependent porosity model (Knox Supergroup only); (5) porosity based on a model based on petrophysical facies; and (6) SREs using National Energy Technology Laboratory’s CO2 Storage prospeCtive Resource Estimation Excel aNalysis (CO2-SCREEN beta V2). All methods used the same values for thickness for each unit. However, the areal extent of each assessment was limited by the data available for each method. Estimated volumes were calculated in 1-by-1 kilometer grid cells and summarized as county and total stratigraphic unit volumes. The resultant SREs mass are displayed using boxplots, which allow for comparing data statistics (mean values and variability) between methods. Differences observed in SRE results from the six methods are mainly attributable to differences in the data and conceptual models used to interpret or estimate porosity in each method. Based on this systematic variability between methods, it is inferred that methods 1, 4, and 6 are best used for regional-scale reconnaissance estimates of storage capacity while methods 2, 3, and 5 are more appropriate for local scales where more data is required. All estimates are data-density dependent and different methods require different amounts of data for reasonable assessments. ArcMap 10.5.1 software was used to portray SREs to help visualize spatial variance of estimates for each methodology, and more importantly, to highlight those areas having the greatest total storage potential estimates.

01 COAL, LIGNITE, AND PEAT↗
DOE OSTI2021

Machine Learning Emulation of Spatial Deposition from a Multi-Physics Ensemble of Weather and Atmospheric Transport Models

In the event of an accidental or intentional hazardous material release in the atmosphere, researchers often run physics-based atmospheric transport and dispersion models to predict the extent and variation of the contaminant spread. These predictions are imperfect due to propagated uncertainty from atmospheric model physics (or parameterizations) and weather data initial conditions. Ensembles of simulations can be used to estimate uncertainty, but running large ensembles is often very time consuming and resource intensive, even using large supercomputers. In this paper, we present a machine-learning-based method which can be used to quickly emulate spatial deposition patterns from a multi-physics ensemble of dispersion simulations. We use a hybrid linear and logistic regression method that can predict deposition in more than 100,000 grid cells with as few as fifty training examples. Logistic regression provides probabilistic predictions of the presence or absence of hazardous materials, while linear regression predicts the quantity of hazardous materials. The coefficients of the linear regressions also open avenues of exploration regarding interpretability—the presented model can be used to find which physics schemes are most important over different spatial areas. A single regression prediction is on the order of 10,000 times faster than running a weather and dispersion simulation. However, considering the number of weather and dispersion simulations needed to train the regressions, the speed-up achieved when considering the whole ensemble is about 24 times. Ultimately, this work will allow atmospheric researchers to produce potential contamination scenarios with uncertainty estimates faster than previously possible, aiding public servants and first responders.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

A Spatiotemporal-Aware Weighting Scheme for Improving Climate Model Ensemble Predictions

Multimodel ensembling has been widely used to improve climate model predictions, and the improvement strongly depends on the ensembling scheme. In this work, we propose a Bayesian neural network (BNN) ensembling method, which combines climate models within a Bayesian model averaging framework, to improve the predictive capability of model ensembles. Our proposed BNN approach calculates spatiotemporally varying model weights and biases by leveraging individual models' simulation skill, calibrates the ensemble prediction against observations by considering observation data uncertainty, and quantifies epistemic uncertainty when extrapolating to new conditions. More importantly, the BNN method provides interpretability about which climate model contributes more to the ensemble prediction at which locations and times. Thus, beyond its predictive capability, the method also brings insights and understanding of the models to guide further model and data development. In this study, we design experiments using an ensemble of CMIP6 climate model simulations to illustrate the BNN ensembling method's capability with respect to prediction accuracy, interpretability, and uncertainty quantification (UQ). We demonstrate that BNN can correctly assign larger weights to the regions and seasons where the individual model fits the observation better. Moreover, its offered interpretability is consistent with our understanding of localized climate model performance. Additionally, BNN shows an increasing uncertainty when the prediction is farther away from the period with constrained data, which appropriately reflects our trustworthiness of the models in the changing climate.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Effective point-charge analysis of crystal fields: Application to rare-earth pyrochlores and tripod kagome magnets 𝑅 ⁢3⁢ Mg⁢ 2 ⁢Sb⁢ 3⁢ O⁢ 14

An indispensable step to understand collective magnetic phenomena in rare-earth compounds is the determination of spatially anisotropic single-ion properties resulting from spin-orbit coupling and crystal field (CF). The CF Hamiltonian has a discrete energy spectrum—accessible to spectroscopic probes such as neutron scattering—controlled by a number of independent parameters reflecting the point symmetry of the magnetic sites. Determining these parameters in low-symmetry systems is often challenging. Here, we describe a general method to analyze CF excitation spectra using adjustable effective point-charges. We benchmark our method to existing neutron-scattering measurements on pyrochlore rare-earth oxides and obtain a universal point-charge model that describes a large family of related materials. We adapt this model to the newly discovered tripod kagome magnets (𝑅 3 ⁢Mg 2 ⁢Sb 3 ⁢O 14 , 𝑅 = Tb, Ho, Er, Yb) for which we report broadband inelastic neutron-scattering spectra. Analysis of these data using adjustable point-charges yields the CF wave functions for each compound. From this, we calculate thermomagnetic properties that accurately reflect our measurements on powder samples and predict the effective gyromagnetic tensor for pseudospin degrees of freedom—a crucial step to understand the exotic collective properties of these kagome magnets at low temperature. We present further applications of our method to other tripod kagome materials and triangular rare-earth compounds 𝑅⁢MgGaO 4 (𝑅 =Yb, Tm). Overall, this study establishes a widely applicable methodology to predict CF and single-ion properties of rare-earth compounds based on interpretable and adjustable models of effective point charges.

Crystal-field theory↗
DOE OSTI2025

A mixture parameterized biologically based dosimetry model to predict body burdens of polycyclic aromatic hydrocarbons in developmental zebrafish toxicity assays

Polycyclic aromatic hydrocarbons (PAHs) are a group of environmental toxicants found ubiquitously as complex mixtures in human-impacted environments. Developmental zebrafish exposures have been used widely to study PAH toxicity, but most studies report nominal exposure concentrations. Nominal exposure concentrations can be unreliable dose metrics due to differences in toxicant bioavailability resulting from disparate exposure methodologies and chemical properties. Toxicokinetic modeling can predict toxicant tissue doses to facilitate comparison between exposures of different chemicals, methodologies, and biological models. We parameterize a biologically based dosimetry model for developmental zebrafish toxicity assays for 9 PAHs. The model was optimized with measurements from media, tissue, and plastic plate walls throughout a static developmental exposure to a mixture of 10 PAHs of high abundance within the Portland Harbor Superfund Site. Plate binding, volatilization, zebrafish permeability, and tissue—media partitioning coefficients vary widely between PAHs. Model predictions accounted for 83% and 54% of 48 hpf body burdens within a factor of 2 resulting from exposures to mixtures and individual PAHs, respectively. Accounting for solubility significantly improves model performance. Competition for active sites in metabolizing enzymes may change biotransformation kinetics between individual PAH and mixture exposures. Area under the curve estimations of concentrations in zebrafish resulted in altered hazard rankings from nominal exposure concentrations. Future work will be oriented to generalizing the model to other PAHs. This PAH dosimetry model improves the interpretability of developmental zebrafish toxicity assays by providing time-resolved body burdens from nominal exposure concentrations.

63 RADIATION, THERMAL, AND OTHER ENVIRON. POLLUTAN↗