Search NASA⌕ Search

SEARCH · Search NASA

Results for “biokinetic modeling”

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.

14 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Mathematical solutions in internal dose assessment: A comparison of Python-based differential equation solvers in biokinetic modeling

Abstract In biokinetic modeling systems employed for radiation protection, biological retention and excretion have been modeled as a series of discretized compartments representing the organs and tissues of the human body. Fractional retention and excretion in these organ and tissue systems have been mathematically governed by a series of coupled first-order ordinary differential equations (ODEs). The coupled ODE systems comprising the biokinetic models are usually stiff due to the severe difference between rapid and slow transfers between compartments. In this study, the capabilities of solving a complex coupled system of ODEs for biokinetic modeling were evaluated by comparing different Python programming language solvers and solving methods with the motivation of establishing a framework that enables multi-level analysis. The stability of the solvers was analyzed to select the best performers for solving the biokinetic problems. A Python-based linear algebraic method was also explored to examine how the numerical methods deviated from an analytical or semi-analytical method. Results demonstrated that customized implicit methods resulted in an enhanced stable solution for the inhaled 60 Co (Type M) and 131 I (Type F) exposure scenarios for the inhalation pathway of the International Commission on Radiological Protection (ICRP) Publication 130 Human Respiratory Tract Model (HRTM). The customized implementation of the Python-based implicit solvers resulted in approximately consistent solutions with the Python-based matrix exponential method ( expm ). The differences generally observed between the implicit solvers and expm are attributable to numerical precision and the order of numerical approximation of the numerical solvers. This study provides the first analysis of a list of Python ODE solvers and methods by comparing their usage for solving biokinetic models using the ICRP Publication 130 HRTM and provides a framework for the selection of the most appropriate ODE solvers and methods in Python language to implement for modeling the distribution of internal radioactivity.

61 RADIATION PROTECTION AND DOSIMETRY↗
DOE OSTI2023

An age- and sex-specific biokinetic model for radon *

Publication 137 of the International Commission on Radiological Protection (ICRP) describes a biokinetic model for radon used to derive dose coefficients for occupational intake of radon isotopes. The model depicts transfer of inhaled or ingested radon to blood, exchange of radon between blood and tissues, and gradual loss of radon from the body based on physical laws governing transfer of a non-reactive and soluble gas between materials. Here, this paper describes an age- and sex-specific variation of that model developed for use in an upcoming ICRP series of reports on environmental intake of radionuclides by members of the public titled ‘Dose Coefficients for Intakes of Radionuclides by Members of the Public’. The proposed model modifies the model structure and transfer coefficients presented in Publication 137 to allow more realistic dosimetric treatment of bone marrow and breast and expands the model to address pre-adult ages, based on the physical principles used in the development of the model of Publication 137 together with anatomical and physiological changes occurring during human development.

61 RADIATION PROTECTION AND DOSIMETRY↗
DOE OSTI2026

Statistical Uncertainty of Inhalation Dose Coefficients: Impact of Particle Deposition in ICRP 66 Human Respiratory Tract Model

Inhaled radioactive materials can pose a long-term health concern, as the material can be incorporated into the body’s metabolic pathways and remain in organs and tissues for extended durations. During the retention period, the radioactive material may localize in a source organ and irradiate adjacent target organs and tissues. Distribution of these materials changes over time, requiring biokinetic modeling to evaluate their movement through various tissues and organs. The evolving distribution depends on multiple inputs characterizing the inhaled material, such as particle size and size distribution, particle density, aspect ratio, specific radionuclide, the chemical form, and solubility. In addition, biological parameters such as breathing rate, breathing type (nasal or nasal/oral), respiratory system morphometry, tidal volume, functional residual capacity, and anatomical dead space all influence material transport. These aerosol properties and physiological characteristics of the respiratory tract jointly define a range of initial conditions that influence the time-dependent distribution of radioactive material. To evaluate both uncertainty in the initial conditions of inhalation exposure and the final output (committed effective dose) from biokinetic models, a Python-based software tool, Radiological Exposure Dose Calculator (REDCAL), was developed to propagate uncertainty within the human respiratory tract model. Focusing on deposition fraction uncertainty, the primary objective was to characterize the initial activity distribution across respiratory regions as a function of anticipated particle sizes and distributions. The impact of the deposition fraction uncertainty was propagated to committed effective dose coefficients for selected radionuclides in a companion publication. For each particle size, a lognormal distribution, characterized by its geometric mean as defined within ICRP Publication 66, serves as the basis for introducing uncertainty into the physical processes governing deposition in various lung regions. Finally, this study addresses the deposition process and examines how uncertainty in deposition mechanisms affects activity distribution in the airways, ultimately presenting the expected range and standard deviation of deposited activity as a function of particle size.

International Commission on Radiological Protectio↗
DOE OSTI2023

Dose reconstruction for plutonium-239 intakes at the Rocky Flats Plant

Purpose: The Rocky Flats (RF) Plant was a weapons manufacturing facility that operated from the early 1950s to 1989. Its primary missions were the production of plutonium (Pu) pits for thermonuclear weapons and the processing of retired weapons for Pu recovery. The purpose of this study was to estimate radiation doses to a cohort of 4499 RF workers from an intake of 239 Pu, the primary plutonium isotope handled at the site. Materials and methods: The latest biokinetic models of the International Commission on Radiological Protection, or site-specific variations of those models, were used to estimate 239 Pu intakes for each worker based on model fits to bioassay data often coupled with lung measurements. Results: Urinary excretion and lung retention data for most 239 Pu intakes could be fit reasonably well by a mixture of Pu dioxide and moderately soluble material. For some workers, better fits were obtained by application of other absorption types including Type S, 239 Pu nitrate, or pure 239 Pu dioxide, or by assuming intake via a wound. The lungs typically received the highest tissue doses, with fifty-year committed equivalent doses in the range of 0.5–1 Sv for 275 workers, 1–5 Sv for 115 workers, and greater than 5 Sv for 12 workers. Conclusions: RF was a unique site regarding a large number of lung measurements available for determining the appropriate absorption types for inhaled material. Further, this provided higher confidence in reconstructed 239 Pu doses than is generally gained from urinary data alone.

61 RADIATION PROTECTION AND DOSIMETRY↗
DOE OSTI2024

Separation of Lipoproteins for Quantitative Analysis of 14 C-Labeled Lipid-Soluble Compounds by Accelerator Mass Spectrometry

To date, 14 C tracer studies using accelerator mass spectrometry (AMS) have not yet resolved lipid-soluble analytes into individual lipoprotein density subclasses. The objective of this work was to develop a reliable method for lipoprotein separation and quantitative recovery for biokinetic modeling purposes. The novel method developed provides the means for use of small volumes (10–200 µL) of frozen plasma as a starting material for continuous isopycnic lipoprotein separation within a carbon- and pH-stable analyte matrix, which, following post-separation fraction clean up, created samples suitable for highly accurate 14 C/ 12 C isotope ratio determinations by AMS. Manual aspiration achieved 99.2 ± 0.41% recovery of [5- 14 CH 3 ]-(2R, 4'R, 8'R)-α-tocopherol contained within 25 µL plasma recovered in triacylglycerol rich lipoproteins (TRL = Chylomicrons + VLDL), LDL, HDL, and infranatant (INF) from each of 10 different sampling times for one male and one female subject, n = 20 total samples. Small sample volumes of previously frozen plasma and high analyte recoveries make this an attractive method for AMS studies using newer, smaller footprint AMS equipment to develop genuine tracer analyses of lipophilic nutrients or compounds in all human age ranges.

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

Modeling Plutonium Decorporation in a Female Nuclear Worker Treated with Ca-DTPA after Inhalation Intake

The present work models plutonium (Pu) biokinetics in a female former nuclear worker. Her bioassay measurements are available at the US Transuranium and Uranium Registries. The worker was internally exposed to a plutonium-americium mixture via acute inhalation at a nuclear weapons facility. She was medically treated with injections of 1 g Ca-DTPA on days 0, 5, and 14 after the intake. Between days 0 and 20, fecal and urine samples were collected and analyzed for 239 Pu and 241 Am. Subsequently, she was followed up for bioassay monitoring over 14 y, with additional post-treatment urine samples collected and analyzed for 239 Pu. The uniqueness of this dataset is due to the availability of: (1) both early and long-term bioassay data from a female with plutonium intake; (2) data on chelation therapy for a female; and (3) fecal measurement results. Chelation therapy with Ca- and/or Zn-salts of DTPA is known to aid in reducing the internal radiation dose by enhancing the excretion of plutonium and americium from the body. Such enhancement affects plutonium biokinetics in the human body, posing a challenge to the internal dose assessment. The current radiation dose assessment practice is to exclude the data affected by Ca-DTPA from the analysis. The present analysis is the first to explicitly model the chelation-affected bioassay data in a female by using a newly developed chelation model. Thus, the bioassay data collected during and after the Ca-DTPA administrations were used for biokinetic modeling and dose assessment. The Markov Chain Monte Carlo method was used to investigate model parameter uncertainty, based on the bioassay data and assumed prior probability distributions. A χ 2 /nData (number of data points) ≈ 1 was observed in this study, which indicates self-consistency of the data with the model. Results of this study show that the worker’s 239 Pu intake was 12 Bq, with a committed effective dose to the whole-body of 1.2 mSv and a committed equivalent dose to the bone surfaces, liver, and lungs of 37.8, 9.1, and 0.8 mSv, respectively. This study also discusses the worker’s dose reduction due to chelation treatment.

61 RADIATION PROTECTION AND DOSIMETRY↗
DOE OSTI2026

Statistical Uncertainty of Inhalation Dose Coefficients in Consequence Management: Propagated Dose Uncertainty in ICRP 66 Human Respiratory Tract Model

Reference inhalation dose models rely on deterministic biokinetics and reference computational phantoms, limiting their applicability to the variability present in population-specific exposures encountered in emergency response scenarios. Here, this study introduces REDCAL, a Python-based computational framework developed to propagate uncertainty in inhalation dose coefficients using the International Commission on Radiological Protection (ICRP) Publication 66 Human Respiratory Tract Model. REDCAL integrates ICRP deposition and clearance models, systemic biokinetics, and governing physics principles, and leverages Sandia National Laboratories’ Dakota toolkit for uncertainty quantification via Latin Hypercube Sampling. REDCAL was validated against DCAL, with biokinetic retention results differing by less than 1% and effective dose coefficients by less than 2% across all tested radionuclides. Stochastic sampling introduced variability in dose coefficients, with geometric standard deviations (GSD) in committed effective dose coefficients (CEDC) ranging from 1.0 to 1.5, based on lognormal distribution fits. Analysis demonstrated that variations in the activity median aerodynamic diameter (AMAD) notably influenced the computed CEDC values. Smaller particles (<1 µm) increased doses by 20–30% due to deeper lung deposition and prolonged retention for alpha emitting radionuclides, such as 241 Am and 239 Pu. Radionuclides with fast clearance, such as 133 I, demonstrated a dose reduction exceeding 50%, as AMAD increased beyond 5 µm due to upper airway deposition and rapid mucociliary clearance. The greatest GSD among the radionuclides reported in this study was for 241 Am. In most cases, the largest GSDs in the CEDC were associated with larger particle sizes, an expected outcome, as ICRP Publication 66 defines GSD in particle size as a function of AMAD, resulting in an extended tail of the lognormal distribution. The findings support improved inhalation dose assessments and enhance consequence management strategies for the U.S. Federal Radiological Monitoring and Assessment Center by quantifying uncertainty in dose coefficients and strengthening decision-making for emergency response scenarios.

Biokinetic Modeling↗
DOE OSTI2024

Analysis of biokinetic parameters reveals patterns in mercury accumulation across aquatic species

Mercury (Hg) is a potent neurotoxicant and poses a risk to human health through the ingestion of Hg-contaminated fish. Mercury, especially in its organic form methylmercury (MeHg), biomagnifies up food chains such that even small aqueous concentrations of Hg can result in significant concentrations of total Hg in fish. Understanding the ecological and human health risks associated with Hg and MeHg exposure requires an understanding of the factors that affect its bioaccumulation in aquatic species. We compiled estimates of three biokinetic parameters: uptake rate (k u ), assimilation efficiency (AE), and efflux rate (k e ). These parameters describe contaminant uptake from aqueous (k u ) and dietary (AE) exposure and the rate of excretion (k e ). We found parameter values for 38 and 34 different species of fish and aquatic invertebrates, respectively, and collected 502 parameter values in total. Here, we used a machine learning technique to establish the relationships between experimental and physiological variables and these parameter values. We found differences in which variables were associated with biokinetic parameter values for fish and aquatic invertebrates. The form of Hg was the most impactful variable, influencing values of all parameters except k u for invertebrates, for which aqueous exposure time was the only significant predicator variable. The parameter k e were the only values significantly influenced by more than one variable, with water type (freshwater, brackish, or marine), organism weight, and form of Hg significantly impacting parameter values for fish and/or invertebrates. To our knowledge, this study represents the most extensive review of biokinetic parameters of Hg and MeHg accumulation in aquatic organisms. Environmental parameters found to significantly impact Hg and MeHg bioaccumulation in past studies were not identified as important in our analyses across aquatic ecosystems and species. Our dataset and analysis reveal novel patterns that may help us better understand and manage Hg bioaccumulation.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS1979

Man-machine analysis of translation and work tasks of Skylab films

An objective approach to determine the concurrent validity of computer-graphic models is real time film analysis. This technique was illustrated through the procedures and results obtained in an evaluation of translation of Skylab mission astronauts. The quantitative analysis was facilitated by the use of an electronic film analyzer, minicomputer, and specifically supportive software. The uses of this technique for human factors research are: (1) validation of theoretical operator models; (2) biokinetic analysis; (3) objective data evaluation; (4) dynamic anthropometry; (5) empirical time-line analysis; and (6) consideration of human variability. Computer assisted techniques for interface design and evaluation have the potential for improving the capability for human factors engineering.

Hosler, W. W.↗
DOE OSTI2023

Chelation Modeling of a Plutonium-238 Inhalation Incident Treated with Delayed DTPA

This work describes an analysis, using a previously established chelation model, of the bioassay data collected from a worker who received delayed chelation therapy following a plutonium-238 inhalation. The details of the case have already been described in two publications. The individual was treated with Ca-DTPA via multiple intravenous injections and then nebulizations beginning several months after the intake and continuing for four years. The exact date and circumstances of the intake are unknown. However, interviews with the worker suggested that the intake occurred via inhalation of a soluble plutonium compound. The worker provided daily urine and fecal bioassay samples throughout the chelation treatment protocol, including samples collected before, during, and after the administration of Ca-DTPA. Unlike the previous two publications presenting this case, the current analysis explicitly models the combined biokinetics of the plutonium-DTPA chelate. Further, using the previously established chelation model, it was possible to fit the data through optimizing only the intake (day and magnitude), solubility, and absorbed fraction of nebulized Ca-DTPA. This work supports the hypothesis that the efficacy of the delayed chelation treatment observed in this case results mainly from chelation of cell-internalized plutonium by Ca-DTPA (intracellular chelation). It also demonstrates the validity of the previously established chelation model. As the bioassay data were modified to ensure data anonymization, the calculation of the “true” committed effective dose was not possible. However, the treatment-induced dose inhibition (in percentage) was calculated.

63 RADIATION, THERMAL, AND OTHER ENVIRON. POLLUTAN↗
DOE OSTI2025

Uncertainty Analysis of Inhalation Dose Coefficients for Nuclear Incident Response

This study addressed the need to characterize variability in inhalation dose coefficients due to uncertainties in respiratory tract deposition, systemic biokinetics, and physiological parameters. A Python-based implementation of the International Commission on Radiological Protection Publication 66 Human Respiratory Tract Model was developed called the Radiological Exposure Dose Calculator (REDCAL) to propagate parameter uncertainty.

61 RADIATION PROTECTION AND DOSIMETRY↗
DOE OSTI2023

Thermal adaptation of soil microbial growth traits in response to chronic warming

ABSTRACT Adaptation of soil microbes due to warming from climate change has been observed, but it remains unknown what microbial growth traits are adaptive to warming. We studied bacterial isolates from the Harvard Forest Long-Term Ecological Research site, where field soils have been experimentally heated to 5°C above ambient temperature with unheated controls for 30 years. We hypothesized that Alphaproteobacteria from warmed plots have (i) less temperature-sensitive growth rates; (ii) higher optimum growth temperatures; and (iii) higher maximum growth temperatures compared to isolates from control plots. We made high-throughput measurements of bacterial growth in liquid cultures over time and across temperatures from 22°C to 37°C in 2–3°C increments. We estimated growth rates by fitting Gompertz models to the growth data. Temperature sensitivity of growth rate, optimum growth temperature, and maximum growth temperature were estimated by the Ratkowsky 1983 model and a modified Macromolecular Rate Theory (MMRT) model. To determine evidence of adaptation, we ran phylogenetic generalized least squares tests on isolates from warmed and control soils. Our results showed evidence of adaptation of higher optimum growth temperature of bacterial isolates from heated soils. However, we observed no evidence of adaptation of temperature sensitivity of growth and maximum growth temperature. Our project begins to capture the shape of the temperature response curves, but illustrates that the relationship between growth and temperature is complex and cannot be limited to a single point in the biokinetic range. IMPORTANCE Soils are the largest terrestrial carbon sink and the foundation of our food, fiber, and fuel systems. Healthy soils are carbon sinks, storing more carbon than they release. This reduces the amount of carbon dioxide released into the atmosphere and buffers against climate change. Soil microbes drive biogeochemical cycling and contribute to soil health through organic matter breakdown, plant growth promotion, and nutrient distribution. In this study, we determined how soil microbial growth traits respond to long-term soil warming. We found that bacterial isolates from warmed plots showed evidence of adaptation of optimum growth temperature. This suggests that increased microbial biomass and growth in a warming world could result in greater carbon storage. As temperatures increase, greater microbial activity may help reduce the soil carbon feedback loop. Our results provide insight on how atmospheric carbon cycling and soil health may respond in a warming world.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2024

Distribution of plutonium and radium in the human heart

Since 1968, the United States Transuranium and Uranium Registries (USTUR) has studied the biokinetics and tissue dosimetry of uranium and transuranium elements in nuclear workers. As part of the USTUR collaboration with the Million Person Study of Low-Dose Health Effects, radiation dose to different parts of the human heart is being estimated for workers with documented intakes of 239 Pu or 226 Ra. The study may be expanded for workers with intakes of 238 U and other radionuclides. The distribution of radionuclides, expressed in terms of concentration (Bq per kg of tissue) serves as an important parameter for estimating radiation dose. Based on available organs from workers who donated their bodies or tissues for research, nine undissected hearts were selected: seven from USTUR registrants with plutonium exposure (males) and two individuals with radium intakes (female and male). For the plutonium workers, estimated 239 Pu systemic deposition ranged from <74 Bq to 1765 Bq. Estimated 226 Ra ‘initial systemic intakes’ were 10.1 MBq and 14.8 kBq for the female patient and male worker, respectively. Organ dissection was based on a heart model published by Borrego et al (2019 J. Radiol. Prot. 39 950–65). This model includes nine cardiac substructures: aorta, left main coronary artery, left atrium, left anterior descending artery, left circumflex artery, left ventricle, right atrium, right coronary artery, and right ventricle. In addition, heart valves, fat attached to epicardium, fluids, and a coronary bypass graft were collected resulting in 111 samples that are currently undergoing radiochemical analyses and mass-spectrometric measurements. The 239 Pu and 226 Ra evaluations are not completed. The results of this study are intended to support radiation worker health studies by improving associated dosimetric and epidemiological models.

USTUR↗