THE MAYAK WORKER DOSIMETRY SYTEM (MWDS-2013): DETERMINATION OF THE INDIVIDUAL SCENARIO OF INHALED PLUTONIUM INTAKE IN THE MAYAK WORKERS

Plutonium in Manhattan Project workers: Using autopsy data to evaluate organ content and dose estimates based on urine bioassay with implications for radiation epidemiology

PURPOSE

Radiation dose estimates in epidemiology typically rely on intake predictions based on urine bioassay measurements. The purpose of this article is to compare the conventional dosimetric estimates for radiation epidemiology with the estimates based on additional post-mortem tissue radiochemical analysis results.

METHODS

The comparison was performed on a unique group of 11 former Manhattan Project nuclear workers, who worked with plutonium in the 1940s, and voluntarily donated their bodies to the United States Transuranium and Uranium Registries.

RESULTS

Post-mortem organ activities were predicted using different sets of urine data and compared to measured activities. Use of urinalysis data collected during the exposure periods overestimated the systemic (liver+skeleton) deposition of 239Pu by 155±134%, while the average bias from using post-exposure urinalyses was -4±50%. Committed effective doses estimated using early urine data differed from the best estimate by, on average, 196±193%; inclusion of follow-up urine measurements in analyses decreased the mean bias to 0.6±36.3%. Cumulative absorbed doses for the liver, red marrow, bone surface, and brain were calculated for the actual commitment period.

CONCLUSION

On average, post-exposure urine bioassay results were in good agreement with post-mortem tissue analyses and were more reliable than results of urine bioassays collected during the exposure.

Lung Cancer in the Mayak Workers Cohort: Risk Estimation and Uncertainty Analysis

The workers at the Mayak nuclear facility near Ozyorsk, Russia are a primary source of information about exposure to radiation at low-dose rates, since they were subject to protracted exposures to external gamma rays and to internal exposures from plutonium inhalation. Here we re-examine lung cancer mortality rates and assess the effects of external gamma and internal plutonium exposures using recently developed Monte Carlo dosimetry systems. Using individual lagged mean annual lung doses computed from the dose realizations, we fit excess relative risk (ERR) models to the lung cancer mortality data for the Mayak Workers Cohort using risk-modeling software. We then used the corrected-information matrix (CIM) approach to widen the confidence intervals of ERR by taking into account the uncertainty in doses represented by multiple realizations from the Monte Carlo dosimetry systems. Findings of this work revealed that there were 930 lung cancer deaths during follow-up. Plutonium lung doses (but not gamma doses) were generally higher in the new dosimetry systems than those used in the previous analysis. This led to a reduction in the risk per unit dose compared to prior estimates. The estimated ERR/Gy for external gamma-ray exposure was 0.19 (95% CI: 0.07 to 0.31) for both sexes combined, while the ERR/Gy for internal exposures based on mean plutonium doses were 3.5 (95% CI: 2.3 to 4.6) and 8.9 (95% CI: 3.4 to 14) for males and females at attained age 60. Accounting for uncertainty in dose had little effect on the confidence intervals for the ERR associated with gamma-ray exposure, but had a marked impact on confidence intervals, particularly the upper bounds, for the effect of plutonium exposure [adjusted 95% CIs: 1.5 to 8.9 for males and 2.7 to 28 for females]. In conclusion, lung cancer rates increased significantly with both external gamma-ray and internal plutonium exposures. Accounting for the effects of dose uncertainty markedly increased the width of the confidence intervals for the plutonium dose response but had little impact on the external gamma dose effect estimate. Adjusting risk estimate confidence intervals using CIM provides a solution to the important problem of dose uncertainty. This work demonstrates, for the first time, that it is possible and practical to use our recently developed CIM method to make such adjustments in a large cohort study.

The Mayak Worker Dosimetry System (MWDS-2013): Internal Dosimetry Results

The distribution of calculated internal doses has been determined for 8043 Mayak Production Associate (Mayak PA) workers. This is a subset of the entire cohort of 25 757 workers, for whom monitoring data are available. Statistical characteristics of point estimates of accumulated doses to 17 different tissues and organs and the uncertainty ranges were calculated. Under the MWDS-2013 dosimetry system, the mean accumulated lung dose was 185 ± 594 mGy (geometric mean = 28 mGy; geometric standard deviation = 9.32; median value = 31 mGy; maximum value = 8980 mGy). The ranges of relative standard uncertainty were from 40 to 2200% for accumulated lung dose, from 25-90% to 2600-3000% for accumulated dose to different regions of respiratory tract, from 13-22% to 2300-2500% for systemic organs and tissues. The Mayak PA workers accumulated internal plutonium lung dose is shown to be close to log normal. The accumulated internal plutonium dose to systemic organs was close to a log triangle. The dependency of uncertainty of accumulated absorbed lung and liver doses on the dose estimates itself is also shown. The accumulated absorbed doses to lung, alveolar-interstitial region, liver, bone surface cells and red bone marrow calculated both with MWDS-2013 and MWDS-2008 have been compared. In general, the accumulated lung doses increased by a factor of 1.8 in median value, while the accumulated doses to systemic organs decreased by factor of 1.3-1.4 in median value. For the cases with identical initial data, accumulated lung doses increased by a factor of 2.1 in median value, while accumulated doses to systemic organs decreased by 8-13% in median value. For the cases with both identical initial data and all of plutonium activity in urine measurements above the decision threshold, accumulated lung doses increased by a factor of 2.7 in median value, while accumulated doses to systemic organs increased by 6-12% in median value.

THE MAYAK WORKER DOSIMETRY SYSTEM (MWDS-2013): TREATMENT OF UNCERTAINTY IN MODEL PARAMETERS

Different dose estimates have been produced for the Mayak PA workforce over recent years (DOSES-2000, DOSES-2005, MWDS-2008). The dosimetry system MWDS-2013 described here differs from previous analyses, in that it deals directly with uncertainty in the assumed model parameters. This paper details the way in which uncertainty is dealt with within MWDS-2013 to produce the final output represented by a multiple hyper-realisation of organ doses. More specifically, the paper describes: Application of the WeLMoS method to calculate Bayesian posterior probability distributions of organ doses.Extension of the WeLMoS method for dealing with multiple intake regimes.How shared and unshared parameters are dealt with using a multiple realisation method.A practical algorithm for the generation of multiple hyper-realisations.How to deal with uncertainty in the intake and the intake regime. The resulting multiple hyper-realisation contains all of the information required to take account of model parameter uncertainty and the effects of shared and unshared parameters in any epidemiological analysis, which uses this information, although it is acknowledged that in practice, certain data simplifications may be required to make such analyses tractable, and comparable to previous analyses. Such simplifications are outside the scope of this paper.

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.