Uncertainties in internal doses calculated for Mayak workers--a study of 63 cases

The Importance and Quantification of Plutonium Binding in Human Lungs

Epidemiological studies have shown that the main risk arising from exposure to plutonium aerosols is lung cancer, with other detrimental effects in the bone and liver. A realistic assessment of these risks, in turn, depends on the accuracy of the dosimetric models used to calculate doses in such studies. A state-of-the-art biokinetic model for plutonium, based on the current International Commission on Radiological Protection biokinetic model, has been developed for this purpose in an epidemiological study involving the plutonium exposure of Mayak workers in Ozersk, Russia. One important consequence of this model is that the lung dose is extremely sensitive to the fraction (fb) of plutonium, which becomes bound to lung tissue after it dissolves. It has been shown that if just 1% of the material becomes bound in the bronchial region, this will double the lung dose. Furthermore, fb is very difficult to quantify from experimental measurements. This paper summarizes the work carried out thus far to quantify fb. Bayesian techniques have been used to analyze data from different sources, including both humans and dogs, and the results suggest a small, but nonzero, fraction of < 1%. A Bayesian analysis of 20 Mayak workers exposed to plutonium nitrate suggests an fb between 0 and 0.3%. Based on this work, the International Commission on Radiological Protection is currently considering the adoption of a value of 0.2% for the default bound fraction for all actinides in its forthcoming recommendations on internal dosimetry. In an attempt to corroborate these findings, further experimental work has been carried out by the US Transuranium and Uranium Registries. This work has involved direct measurements of plutonium in the respiratory tract tissues of workers who have been exposed to soluble plutonium nitrate. Without binding, one would not expect to see any activity remaining in the lungs at long times after exposure since it would have been cleared by the natural process of mucociliary clearance. Further supportive study of workers exposed to plutonium oxide is planned. This paper ascertains the extent to which these results corroborate previous inferences concerning the bound fraction.

The Mayak Worker Dosimetry System (MWDS-2013): A Bayesian Analysis to Quantify Pulmonary Binding of Plutonium in Lungs Using Historic Beagle Dog Data

The revised human respiratory tract model, published in Part 1 of the International Commission on Radiological Protection's (ICRP) report on Occupational Intakes of Radionuclides (OIR), includes a bound fraction, fb, to represent radionuclides that have become chemically bound in the lungs following dissolution of particulates in lung fluid. Bound radionuclides are not subject to particle transport clearance but can be absorbed to blood at a rate, sb. The occurrence of long-term binding of plutonium can greatly increase lung doses, particularly if it occurs in the bronchial and bronchiolar regions. However, there has been little evidence that currently supports the existence of a long-term bound state for plutonium. The present work describes the analysis of measurements of lung data obtained from a life span study of Beagle dogs that were exposed by inhalation to different concentrations of plutonium-239 (239Pu) nitrate aerosol at Pacific Northwest Laboratories, USA. The data have been analysed to assess whether a bound state was required to explain the data. A Bayesian approach was adopted for the analysis that accounts for uncertainties in model parameter values, including uncertainties in the rates of particle transport clearance. Furthermore, it performs the analysis using two different modelling hypotheses: a model based on the current ICRP human respiratory tract model and its treatment of alveolar particle transport clearance; and a model of particle transport clearance that is based on the updated model developed by ICRP to calculate dose coefficients for the OIR. The current model better represents clearance in dogs at early times (up to 1 year following intake) and the latter better represents retention at greater times (>5 years following intake). The results indicate that a long-term bound fraction of between 0.16 and 1.1%, with a mean value of between 0.24 and 0.8% (depending on the model) is required to explain the data.

The Mayak Worker Dosimetry System (MWDS-2013): Plutonium Binding in the Lungs-An Analysis of Mayak Workers

Estimates of plutonium lung doses from urine bioassay are highly dependent on the rate of absorption from the lungs to blood assumed for the inhaled aerosol. Absorption occurs by dissolution of particles in lung fluid followed by uptake to blood. The latter may occur either rapidly or dissolved ions may first become temporarily bound within airway tissue. The presence of long-term binding can greatly increase lung doses, particularly if it occurs in the bronchial and bronchiolar regions. Analyses of autopsy data from Beagle dogs and USTUR Case 0269, obtained following exposure to plutonium nitrate, suggest that a small fraction of 0.2-1.1 and 0.4-0.7%, respectively, of plutonium becomes permanently bound within the lungs. The present work performs a further analysis using autopsy data of former plutonium workers of the Mayak Production Association to determine values of the bound fraction that are supported by these data. The results suggest a bound fraction value of 0-0.3%. The results also indicate that the Mayak worker population median values of the particle transport clearance parameters from the alveolar-interstitial region are largely consistent with expected values, but suggest the rate from the alveolar region to the interstitium may be lower than initially thought.

The Mayak Worker Dosimetry System (Mwds-2013): Plutonium Dissolution in The Lungs-An Analysis of Mayak Workers

Lung doses resulting from inhalation of plutonium aerosols are highly dependent on the assumed rate of particle clearance, which occurs by two competing processes: (1) particle transport clearance to the alimentary tract and to the thoracic lymph nodes and (2) clearance to systemic tissues, which occurs by dissolution of particles in lung fluid followed by uptake to blood, which is a process collectively known as absorption. Unbiased and accurate estimates of the values of lung absorption parameters are required to obtain reliable estimates of lung dose, particularly those inferred from urine bioassay. Parameter values governing the rate of absorption are best estimated from data, such as autopsy measurements of plutonium in the lungs and systemic tissues, which directly relate to the exposed workers of interest. However, because the mathematical models that determine clearance from the lungs and systemic tissues are complex and consist of many parameters, estimates of model parameter values are subject to significant uncertainties. With this in mind, this paper uses a Bayesian approach to estimate one of the most important dissolution parameters: the slow rate of dissolution. This is estimated for both plutonium nitrate and plutonium oxide bearing aerosols in the lungs of former workers of the Mayak Production Association. A value of 2.6 × 10-4 d-1 is estimated for plutonium nitrates, and 4.7 × 10-5 d-1 for plutonium oxides.

THE MAYAK WORKER DOSIMETRY SYSTEM (MWDS 2013): A RE-ANALYSIS OF USTUR CASE 0269 TO DETERMINE WHETHER PLUTONIUM BINDS TO THE LUNGS

Radionuclides in ionic form can become chemically bound in the airways of the lungs following dissolution of inhaled particulates in lung fluid. The presence of long-term binding can greatly increase lung doses from inhaled plutonium, particularly if it occurs in the bronchial and bronchiolar regions. However, the only published evidence that plutonium binding occurs in humans comes from an analysis of the autopsy and bioassay data of United States Transuranium and Uranium Registries Case 0269, a plutonium worker who experienced a very high (58 kBq) acute inhalation of plutonium nitrate. This analysis suggested a bound fraction of around 8 %, inferred from an unexpectedly low ratio of estimated total thoracic lymph node activity:total lung activity, at the time of death. However, there are some limitations with this study, the most significant being that measurements of the regional distribution of plutonium activity in the lungs, which provide more direct evidence of binding, were not available when the analysis was performed. The present work describes the analysis of new data, which includes measurements of plutonium activity in the alveolar-interstitial (AI) region, bronchial (BB) and bronchiolar (bb) regions, and extra-thoracic (ET) regions, at the time of death. A Bayesian approach is used that accounts for uncertainties in model parameter values, including particle transport clearance, which were not considered in the original analysis. The results indicate that a long-term bound fraction between 0.4 and 0.7 % is required to explain this data, largely because plutonium activity is present in the extra-thoracic (ET2), bronchial and bronchiolar airways at the time of death.

An intake prior for the Bayesian analysis of plutonium and uranium exposures in an epidemiology study

In Bayesian inference, the initial knowledge regarding the value of a parameter, before additional data are considered, is represented as a prior probability distribution. This paper describes the derivation of a prior distribution of intake that was used for the Bayesian analysis of plutonium and uranium worker doses in a recent epidemiology study. The chosen distribution is log-normal with a geometric standard deviation of 6 and a median value that is derived for each worker based on the duration of the work history and the number of reported acute intakes. The median value is a function of the work history and a constant related to activity in air concentration, M, which is derived separately for uranium and plutonium. The value of M is based primarily on measurements of plutonium and uranium in air derived from historical personal air sampler (PAS) data. However, there is significant uncertainty on the value of M that results from paucity of PAS data and from extrapolating these measurements to actual intakes. This paper compares posterior and prior distributions of intake and investigates the sensitivity of the Bayesian analyses to the assumed value of M. It is found that varying M by a factor of 10 results in a much smaller factor of 2 variation in mean intake and lung dose for both plutonium and uranium. It is concluded that if a log-normal distribution is considered to adequately represent worker intakes, then the Bayesian posterior distribution of dose is relatively insensitive to the value assumed of M.

Application of Bayesian inference to the bioassay data from long-term follow-up of two refractory PuO2 inhalation cases

The dominant contribution to the uncertainty in internal dose assessment can often be explained by the uncertainty in the biokinetic model structure and parameters. The International Commission on Radiological Protection (ICRP) is currently updating its biokinetic models, including the Human Respiratory Tract Model (HRTM). Gregoratto et al. (2010) proposed a physiologically-based particle transport model that simplifies significantly the representation of particle clearance from the alveolar interstitial region. Bayesian inference using the Weighted Likelihood Monte-Carlo Sampling (WeLMoS) method is applied to the bioassay and autopsy data from the U.S. Transuranium and Uranium Registries' (USTUR) tissue donors 0202 and 0407 exposed to "high fired," refractory PuO2 aerosols in order to examine the applicability of the revised model and to estimate the uncertainties in model parameters and the lung doses as expressed by the posterior probability distributions. It is demonstrated that, with appropriate adjustments, the Gregoratto et al. particle transport model can describe situations involving exposure to highly insoluble particles. Significant differences are observed in particle clearance pattern characteristics to these two individuals' respiratory systems. The respiratory tract of registrant 0202 was most likely compromised by his prior occupational exposure to coal dust, smoking habit, and chronic obstructive pulmonary disease, while donor 0407 was a non-smoker and had no prior history of lung disorder. However, the central values of the particle transport parameter posterior distributions for both cases are found to be still within the 68% probability range for the inter-subject variability derived by Gregoratto et al. PuO2 particles produced by the plutonium fire were extremely insoluble, with about 99% absorbed into blood at a rate of approximately 4.8 × 10 d (Case 0202) and 5.1 × 10 d (Case 0202). When considering this type of plutonium material, doses to other body organs are small in comparison to those to tissues of the respiratory tract. More than 95% of the total committed weighted equivalent dose is contributed by the lungs.

A method for calculating Bayesian uncertainties on internal doses resulting from complex occupational exposures

Estimating uncertainties on doses from bioassay data is of interest in epidemiology studies that estimate cancer risk from occupational exposures to radionuclides. Bayesian methods provide a logical framework to calculate these uncertainties. However, occupational exposures often consist of many intakes, and this can make the Bayesian calculation computationally intractable. This paper describes a novel strategy for increasing the computational speed of the calculation by simplifying the intake pattern to a single composite intake, termed as complex intake regime (CIR). In order to assess whether this approximation is accurate and fast enough for practical purposes, the method is implemented by the Weighted Likelihood Monte Carlo Sampling (WeLMoS) method and evaluated by comparing its performance with a Markov Chain Monte Carlo (MCMC) method. The MCMC method gives the full solution (all intakes are independent), but is very computationally intensive to apply routinely. Posterior distributions of model parameter values, intakes and doses are calculated for a representative sample of plutonium workers from the United Kingdom Atomic Energy cohort using the WeLMoS method with the CIR and the MCMC method. The distributions are in good agreement: posterior means and Q(0.025) and Q(0.975) quantiles are typically within 20 %. Furthermore, the WeLMoS method using the CIR converges quickly: a typical case history takes around 10-20 min on a fast workstation, whereas the MCMC method took around 12-72 hr. The advantages and disadvantages of the method are discussed.

Uncertainty analysis of doses from ingestion of plutonium and americium

Uncertainty analyses have been performed on the biokinetic model for americium currently used by the International Commission on Radiological Protection (ICRP), and the model for plutonium recently derived by Leggett, considering acute intakes by ingestion by adult members of the public. The analyses calculated distributions of doses per unit intake. Those parameters having the greatest impact on prospective doses were identified by sensitivity analysis; the most important were the fraction absorbed from the alimentary tract, f(1), and rates of uptake from blood to bone surfaces. Probability distributions were selected based on the observed distribution of plutonium and americium in human subjects where possible; the distributions for f(1) reflected uncertainty on the average value of this parameter for non-specified plutonium and americium compounds ingested by adult members of the public. The calculated distributions of effective doses for ingested (239)Pu and (241)Am were well described by log-normal distributions, with doses varying by around a factor of 3 above and below the central values; the distributions contain the current ICRP Publication 67 dose coefficients for ingestion of (239)Pu and (241)Am by adult members of the public. Uncertainty on f(1) values had the greatest impact on doses, particularly effective dose. It is concluded that: (1) more precise data on f(1) values would have a greater effect in reducing uncertainties on doses from ingested (239)Pu and (241)Am, than reducing uncertainty on other model parameter values and (2) the results support the dose coefficients (Sv Bq(-1) intake) derived by ICRP for ingestion of (239)Pu and (241)Am by adult members of the public.

Plutonium worker dosimetry

Epidemiological studies of the relationship between risk and internal exposure to plutonium are clearly reliant on the dose estimates used. The International Commission on Radiological Protection (ICRP) is currently reviewing the latest scientific information available on biokinetic models and dosimetry, and it is likely that a number of changes to the existing models will be recommended. The effect of certain changes, particularly to the ICRP model of the respiratory tract, has been investigated for inhaled forms of (239)Pu and uncertainties have also been assessed. Notable effects of possible changes to respiratory tract model assumptions are (1) a reduction in the absorbed dose to target cells in the airways, if changes under consideration are made to the slow clearing fraction and (2) a doubling of absorbed dose to the alveolar region for insoluble forms, if evidence of longer retention times is taken into account. An important factor influencing doses for moderately soluble forms of (239)Pu is the extent of binding of dissolved plutonium to lung tissues and assumptions regarding the extent of binding in the airways. Uncertainty analyses have been performed with prior distributions chosen for application in epidemiological studies. The resulting distributions for dose per unit intake were lognormal with geometric standard deviations of 2.3 and 2.6 for nitrates and oxides, respectively. The wide ranges were due largely to consideration of results for a range of experimental data for the solubility of different forms of nitrate and oxides. The medians of these distributions were a factor of three times higher than calculated using current default ICRP parameter values. For nitrates, this was due to the assumption of a bound fraction, and for oxides due mainly to the assumption of slower alveolar clearance. This study highlights areas where more research is needed to reduce biokinetic uncertainties, including more accurate determination of particle transport rates and long-term dissolution for plutonium compounds, a re-evaluation of long-term binding of dissolved plutonium, and further consideration of modeling for plutonium absorbed to blood from the lungs.

Uncertainties of Mayak urine data

A method of parameterising the likelihood functions representing the uncertainty of Mayak plutonium urine bioassay measurements is described. The Poisson-lognormal model is assumed and data from 63 cases (1,087 urine measurements in all) are used to empirically determine the lognormal normalisation uncertainty, given the measurement uncertainties obtained from count quantities. An outlier-insensitive procedure is used to fit the cumulative probability distribution of scaled deviations in order to determine the normalisation uncertainty. The natural logarithm of the geometric standard deviation of the total normalisation uncertainty is found to be 0.34 including a measurement component estimated to be 0.2.

IMPDOS (improved dosimetry and risk assessment for plutonium-induced diseases): internal dosimetry software tools developed for the Mayak worker study

A collection of software tools developed for the Mayak worker study is described. IMPDOS is software for modelling, data analysis, and activity and dose calculations using the bioassay and postmortem data from Mayak workers provided by Southern Urals Biophysics Institute. The capabilities include: parameter fitting of data for individual cases, Bayesian dose calculations using the fit results for collections of cases with extensive data as a biokinetic prior, and database storage of results for retrieval, analysis and interpretation.