UNSCEAR makes use of generic radon fluxes to estimate radon release rates and a generic air dispersion model to estimate the environmental radon concentrations as a function of distance from the site. UNSCEAR then converts the concentrations to population doses using assumed population densities for areas out to a distance of 2000 km and a radon dose conversion factor. Doses are accumulated over an assumed long term exposure period (10 000 years). The results are normalised to a unit amount of electrical energy produced. A summary of UNSCEAR assumptions, some of which were originally derived in Ref 4, is given in Table 2. These UNSCEAR assumptions are discussed below.
Uranium Fuel Requirements
UNSCEAR assumes that 210 tU (250 t U3O8) are required to produce 1 GW-year of electrical energy. Requirements are dependent on the reactor type, ranging from 180 tU for heavy water reactors to 330 tU for Magnox reactors (Ref 1, p105).
Normalised Tailings Surface Area
The basis for the 1 hectare per GW-year value is not given (originally used in Ref 5, p140). However, for perspective, the thickness of tailings, with a density of 1.6 t m-3 and a surface area of 1 hectare, resulting from the production of 210 tU from 1% uranium grade ore, would be about 1.4 m (assuming 92% recovery). This value is inversely proportional to the grade, with 0.3% ore requiring a thickness of about 4.8 m. In practice, however, tailings usually exceed these thicknesses, at least for the sites examined in this study, and therefore the area per unit of electrical energy produced is usually less than that assumed by UNSCEAR. (The radon release rate per unit surface area does increase proportionately with ore grade, but only minimally with increasing thickness beyond the first couple of metres of tailings.)
Radon Release Rate
The UNSCEAR unit radon release rate is based on reported emission rates ranging around 10 Bq m-2s-1, and the assumption that some reasonably impermeable cover would reduce the rate to 3 Bq m-2s-1 (Ref 1, p106; Ref 5, p140). The rate is assumed to be unchanged over at least 10 000 years because of the long radioactive half-life (80 000 years) of Th-230, the precursor of Ra-226.
As described in the next chapter, this rate substantially exceeds the rate expected for most of the tailings sites examined in this study, due to the planned covers and/or saturation of the tailings. The normalised emission rate of 1 (actually 0.946) TBq per year per GW-year is derived from the previous assumptions.
Population Density
The assumed population densities for the reference tailings site of 3 and 25 persons km-2 in the <100 km and 100 to 2000 km regions, respectively (Ref 1, p106), were originally derived from a 1975 study of tailings sites in the United States (Ref 4, p168).
While not unreasonable values for the rural, southwest United States, these densities are significantly greater than the densities derived for the sites examined in this study. The final estimate of the population dose is directly proportional to the assumed population density.
Air Dispersion Factors
The average air dispersion factor at 1 km and its reduction with distance for the model site were derived in part from a Gaussian plume model with a nominal source release height of 10 m and various assumptions about atmospheric conditions (Ref 4 and Ref 6). The factor at 1 km is 3 x 10-6 Bq m-3 per Bq s-1. Beyond 1 km, and it was assumed that the concentrations decreased as (distance) -1.5.
If the assumed population densities and air dispersion factors are combined, and the radon progeny equilibrium factors (see below) are assumed to be constant with distance, then the total population dose is proportional to:
where P1 and P2 are the population densities (persons km-2) in the <100 km and 1002000 km regions respectively. Integrating this equation indicates that about 97% of the estimated population dose is for people living more than 100 km from the site for the population densities assumed by UNSCEAR. For a uniform population density across both regions, about 78% of the population dose would be for people living in the 1002000 km region around the site.
It is not clear if UNSCEAR includes the decay of radon (3.82 day half-life) in its dispersion calculations. At a distance of 2000 km, and with an assumed 2.5 m s-1 average wind speed, the decay would reduce the radon concentrations by about a factor of five.
Dose Conversion Factor
The dose conversion factor of 9 nSv h-1 per Bq m-3 (equilibrium equivalent concentration, EEC) (Ref 1, p54) is based on now superseded dose factors that were derived using a dosimetric approach.
Subsequent to the publication of Ref 1, the International Commission on Radiological Protection (ICRP) published (Ref 7) a revised radon dose conversion factor for members of the public (derived by the ICRP using an epidemiological approach). The factor is 4 mSv per WLM (working level month) and is equivalent to 6.4 nSv h-1 per Bq m-3 (EEC), the value used in this study. (Based on the results from a recent meta-analysis of major epidemiological studies (Ref 8), an even lower dose conversion factor (Ref 9) of less than 2 mSv per WLM may be appropriate.)
Radon Progeny Equilibrium Factors
The radon progeny equilibrium factors (Ref 1, p54) were based on several reported studies relative to typical background conditions, although there is some suggestion from some European and US studies that 0.6, rather than 0.8, might be a more representative outdoor factor. Using the 0.6 value would slightly reduce the time-weighted average equilibrium factor (considering indoor and outdoor occupancy) from 0.48 to 0.44.
These factors refer to typical background conditions. However, when relatively close to a source of radon, such as uranium tailings, the outdoor factor is much lower because there is insufficient time for radon progeny in-growth. For example, for an assumed 2.5 m s-1 average wind speed, the outdoor equilibrium factor at 1 km from a radon source would be about 0.1. (The indoor factor would be dependent on the air exchange rate of the building.)
Collective Dose Factors
The final estimate of the collective dose is directly proportional to the assumed cumulative exposure period of 10 000 years. UNSCEAR (Ref 1, p107) acknowledges that the estimated result of 150 person-Sv per GW-year (the same estimate obtained by Ref 5, p140) is highly dependent on a number of assumptions, including future tailings management practices, and they suggest a range from 1 to 1000 person-Sv per GW-year about their central estimate. The numerical basis for UNSCEAR's quoted range is not given.
Radon Source Terms and Population Densities
For this assessment, radon releases from the major, operating uranium production facilities were considered in the evaluation of the post-decommissioning source terms. The analysis did not consider radon from tailings areas no longer in use or which have previously been decommissioned. Similarly, as done by UNSCEAR, radon from mining and milling operations and from waste rock with residual trace radioactivity were not included in this analysis.
In order to develop a better understanding of the potential source terms, the operators of the major currently producing facilities were sent information requests. The following source term descriptions and population densities for the various sites were based on the responses, and on information supplied by the Uranium Institute and derived from publicly available literature.
The post-decommissioning radon source term estimates are summarised in Table 3. The population density estimates are summarised in Table 4. The bases for these estimates are in given in Ref 2.
Air Concentration Modelling
Modelling of long range transport requires sophisticated models, comprehensive meteorological data and extensive set-up effort. Existing data were available at SENES from previous studies to model long range transport for North American sites and two sites were selected. The first site was applicable to northern Saskatchewan meteorological conditions and the Canadian uranium mill tailings examined in this study, while the second (Mexico) was taken as representative of mid latitude meteorological conditions.
Analysis of the dispersion patterns for these two quite different regions provided information on how dispersion with distance varied. This analysis thus allowed some insight into the variability in the patterns of dispersion with distance that could be anticipated for mine sites at different locations. It was the intention of this analysis to illustrate the possible range of uncertainty associated with using the same air dispersion factors for all sites.
Northern Latitude Site
The dispersion modelling for northern Saskatchewan used information and methodology previously developed for a study of the long-range pollutant transport in North America (Ref 17). The US Environmental Protection Agency (EPA) CALPUFF/CALMET modelling package was used to address long range flow patterns and the earth's curvature. Meteorological conditions were estimated for a 100 km by 100 km grid. Set-up effort is extensive for this type of modelling and, therefore, such detailed modelling could not be applied to other sites within the constraints of this project.
Concentrations from 100 to 2000 km were estimated using the CALPUFF/CALMET model and the existing data available to SENES. The US EPA ISC3 model was used to model concentrations from 1 km out to 100 km. The ISC3 model requires less set-up effort and provides reliable estimates over this range of distances. The modelled source, located at 58N and 103W and at a height of 1 m, had an emission rate of 1.0 Bq m-2s-1 over a 250 000 m² area source for a total release rate of 0.25 MBq s-1.
Theory predicts that concentration drops off quickly, by much more than a factor of r-1, close to the source due to the combination of rapidly increasing area with distance and due to vertical mixing. At larger distances, concentrations tend to drop off by no less than r-1 since the mixing heights limit the vertical dispersion.
It is not clear if the UNSCEAR air dispersion modelling was done accounting for radioactive decay of radon (half-life of 3.82 days). As a result, the radon concentrations may be overestimated for this reason, especially at large distances from the source. For this study, a correction for removal of radon due to radioactive disintegration was developed based on the time required for the radon to reach the location. The average duration of transport from the source to the receptor was approximated by dividing the distance by a 2.5 m s-1 average wind speed. For example, the duration of transport to a receptor 2000 km from the source would be 9.3 days.
The corrected concentration (c) of radon at a receptor was approximated by:
where lambdar is the radioactive decay rate, 2.1 x 10-6 s-1, of radon, and t is the elapsed time between release at the source and arrival at the receptor. The time was estimated by dividing the distance by a speed of 2.5 m s-1 (as estimated from the North American meteorology). The method estimates that radon concentrations at 2000 km distance from the source would be about 19% of modelled values assuming no radioactive decay (about a factor of five lower).
Concentrations were summarised for 16 directions from the source and the mean, maximum, and minimum concentrations are plotted on Figure 1 in comparison to those predicted by UNSCEAR dispersion factors. Concentrations drop off rapidly with distance with mean levels decreasing from about 300 mBq m-3 at 1 km from the source to 0.3 mBq m-3 at a distance of 100 km (a factor of 1000 lower). The mean concentration is lower than 0.001 mBq m-3 at a distance of 2000 km and continues to drop at increasing distances due to both ongoing dilution and the radioactive decay of radon. The ratios of the predictions using the UNSCEAR dispersion model to the predicted mean values in this study range from about two at a distance of 10 km to 17 at a distance of 2000 km.
The incremental radon levels at all distances are much lower than typical outdoor radon concentrations, which are in the order of 10 000 mBq m-3 (Ref 1, p54).
Mid Latitude Site
Radon concentrations were estimated for the 100 km by 100 km North American grid with the source located at 25.5N and 103.0W in Mexico. Although no uranium facility is present at this location, the site was selected to illustrate dispersion characteristics in the mid latitudes of the Northern Hemisphere.
Figure 2 shows the pattern of concentrations with distance for both the northern Saskatchewan and mid latitude locations for distances of 100 to 2000 km from the source in comparison to concentrations predicted using UNSCEAR's dispersion model. Mean concentrations for the mid latitude location are higher (by about a factor of two) than those calculated for the northern Saskatchewan location. (Although not shown, the maximum and minimum concentrations (i.e. upwind or downwind) vary by about a factor of 10 for the mid latitude location.)
Reference Concentrations
Concentrations for the sites examined in this study were estimated based on the 1 to 100 km mean concentrations estimated using for northern latitude (Saskatchewan) meteorology, and the 100 to 2000 km mean concentrations using the mid latitude (Mexico) meteorology. Concentrations were prorated by the ratio of site-specific emission rates to the reference case emission rate, 0.25 MBq s-1, used in the dispersion modelling described previously.
The dispersion factors, and hence concentrations, used by UNSCEAR and those derived in this study for the northern Saskatchewan region differ by factors of two to three in the <100 km region. However, the differences increase with distance, with the difference being up to a factor of about 17 (with UNSCEAR concentrations larger) at 2000 km. This may be due in part to the effect of radon decay that was accounted for in this study and which reduces concentrations by about a factor of five at 2000 km.
Population Dose Estimates
Population exposure estimates were based on multiplying the population size in an area by the average radon concentration in the area. This section describes the population exposure estimates for each mining area that were estimated by multiplying the average (area-weighted) site-specific concentration in the <100 km and 100 to 2000 km regions by the site-specific uniform population density in the regions.
Since concentrations decrease rapidly with distance, the true distribution of population within the area can significantly impact the population exposure. In order to examine the potential impacts of using a uniform rather than a true population distribution, population exposure estimates by distance from the source were investigated in detail for the northern Saskatchewan site. Site specific population distributions were not available for the other sites. In this situation, population exposure estimates were found to be about a factor of three to four lower if the true population distribution with distance were used as compared to the assumption of uniform population density within the two regions, i.e. the <100 km and about three times lower for the 100 to 2000 km regions.
The long term (10 000 year) population dose estimates for the uranium tailings sites examined in this study are summarised in Table 5. The normalised estimates range from 0 to 5.9 person-Sv per GW-year, with an overall 1997 production weighted average of 0.96 person-Sv per GW-year. These estimates may be compared to the UNSCEAR central estimate of 150 person-Sv per GW-year. The way in which the estimates for each site were derived is described in Ref 2.
Background Doses
People living around the eight sites examined in this study will be exposed to background radon, irrespective of the operation of uranium production facilities. For the total of 2.11 x 10 (Ref 8) people around the sites (counting only one in each country in order not to double count people), the total background dose is given by:
The assumed 14.4 Bq m-3 (EEC) average background radon concentration was derived from Reference 1 as follows. The previously noted outdoor concentration of 10 000 mBq m-3 (10 Bq m-3) converts to 8 Bq m-3 (EEC) based on an outdoor radon progeny equilibrium factor of 0.8. The indoor concentration of 40 Bq m-3 converts to 16 Bq m-3 (EEC) based on an indoor equilibrium factor of 0.4 (Ref 1, p54). For 80% indoor occupancy and 20% outdoor occupancy, these concentrations give an overall average of 14.4 Bq m-3 (EEC).
This is probably a conservative (low) dose estimate because background concentrations in areas with uranium deposits are generally higher than typical background levels. The dose of 1.7 x 109 person-Sv per 10 000 years is a factor of 366 000 larger than the total population dose of 4650 person-Sv for all the sites examined in this study. These results are discussed in more detail in the next section.
Comparing UNSCEAR and SENES Estimates
The overall, normalised population dose estimate in the SENES study of 0.96 person-Sv per GW-year of electrical energy produced is about a factor of 150 lower than UNSCEAR's estimate1 of 150 person-Sv per GW-year. The more recent and site-specific information obtained for the SENES study was not available to UNSCEAR, which had to rely on generic radon emission and population data for its estimates.
UNSCEAR (Ref 1) provided a range of 1 to 1000 person-Sv per GW-year about their central estimate of 150 person-Sv per GW-year. The central estimates for every site examined in this study were below the UNSCEAR central estimate, while the central estimates for six of the sites were below UNSCEAR's suggested bottom range estimate of 1 person-Sv per GW-year. Some of the more significant parameters that contribute to these differences in the estimates are discussed below.
Radon Emission Rates
Based on the information summarised in Table 3, the long term emission rate of radon from tailings used by UNSCEAR substantially exceeds the rates expected for the most of the tailings sites examined in this study. While the value of the long term emission rate at any site is certainly speculative (especially to 10 000 years), it is clear that UNSCEAR's central estimate did not account for the current and short term future tailings management practices that would essentially eliminate radon emissions at some sites i.e. saturated, water-covered tailings. The lower estimate of 1 person-Sv per GW-year given by UNSCEAR does indicate that it acknowledges the uncertainty in its estimate (and specifically that its estimate may be too large), although the numerical basis for its estimated range is not given by UNSCEAR.
Normalised Emanating Area of Tailings
The UNSCEAR estimate of 1 hectare per GW-year differs from (overestimates) the surface area of tailings based on more recent data, or, conversely, underestimates the thickness of tailings, and correspondingly the amount of potential electrical energy, associated with the radon emissions. This increases the final estimate of person-Sv per GW-year, especially for the deeper (thicker) tailings areas. Deeply buried tailings essentially contribute zero radon to the environment but the ore associated with those tailings does result in the production of electrical power.
Population Density
The overall population densities assumed by UNSCEAR for its modelled tailings site typically overestimate the current population densities at the sites examined in this study. Data collected for this study shows that for the <100 km region, the ratios of the UNSCEAR estimate of 3 persons km-2 to the site-specific densities range from 0.9 to 88.2. The ratios for the 1002000 km region range from 3.5 to 16.7, whereas the UNSCEAR estimate is 25 persons km-2.
This study also examined the effects of using uniform, rather than actual, population distributions as a function of distance around a tailings site. For the Saskatchewan site, for which data were available to this study, the assumption of uniform (i.e. two region) distributions overestimates the cumulative population dose by about a factor of 2.5. While this difference is not necessarily applicable to other sites, the analysis indicates the magnitude of one of the sources of uncertainty associated with the final dose estimates.
Air Dispersion Factors
As in the SENES study, UNSCEAR makes use of generic air dispersion factors in assessing the dispersion of radon released from the decommissioned tailings. The air dispersion factors derived in this study compare within a factor of about three with the UNSCEAR factors at distances <100 km, but diverge from the UNSCEAR factors, by factors of four to 10 and more, in the 1002000 km region.
The analyses carried out in this study suggest that only part of this difference in the distant region is due to site-specific parameters, since the direction-averaged dispersion factors for the northern latitude (Saskatchewan, Canada) and mid latitude sites (Mexico) differed by less than a factor of two. Not accounting for radon decay with distance also appears to have contributed to UNSCEAR's higher estimate (relative to this study) of the air dispersion factors. If radon decay were to be included in the UNSCEAR values, the air dispersion factors derived in this study would be within about a factor of three of the UNSCEAR values.
Uncertainties
Similar to the UNSCEAR long term population dose estimates, the estimates derived in this study are inherently uncertain. Some of the sources of this uncertainty are qualitatively discussed below.
How Representative are the Sites?
The eight sites examined in this study are currently (1997) responsible for 67% of the world's production of uranium. (As previously noted, the results of the study could be considered to represent the impacts of long term radon emissions based on 80% of current worldwide uranium production if allowance is made for in-situ leach facilities.) However, the conditions at other production facilities may be significantly different than at these major producers, in terms of population densities and especially in terms of overall radon emission rates per unit of potential electrical energy produced. The overall population dose for all existing facilities may therefore be different.
Nevertheless, given that tailings management practices are continually evolving and that the sites examined here may be considered representative of likely future practices, the estimate derived in this study is probably a fair representation of conditions in the foreseeable future.
Other Fuel Cycles
The population dose estimates from both this study and UNSCEAR were based on assumed uranium fuel requirements of 210 tU per GW-year. Other fuel cycles, including reprocessing, could significantly lower fuel requirements and would correspondingly result in lower population doses per GW-year due to radon emissions from uranium tailings.
Population Growth
The population dose estimates were based on current (generally within 10 years) population densities. There is the potential for long term population growth but the rate and numbers over the very long term (10 000 years) is unknown. On a worldwide basis, and considering that the earth's resources are finite, the long term population is unlikely to be more than a factor of five to ten larger than today's population. On the other hand, in the <100 km region, because many of the uranium producing facilities are in remote locations, the local populations may well decrease once those facilities are no longer operating.
All this is clearly speculative, but the uncertainties due to population growth affect both the estimates derived in this study and the UNSCEAR estimate. The estimated population doses from background radiation would also be affected by population growth.
Population Distribution
As discussed above, the assumption of uniform (two region) population densities resulted in nearly a threefold overestimation of population dose for the northern latitude site, based on the actual population distribution. A similar or greater magnitude of overestimation or underestimation could exist at any site. The location of the population relative to wind direction causes larger uncertainties, by factors of five to 15 when comparing maximum (downwind) and minimum (upwind) population doses (Table 5). These uncertainties are equally inherent in the UNSCEAR estimate of population dose.
Air Dispersion Factors
Air dispersion characteristics, especially near-field, can be very site-specific. However, except for sites that might have population centres located generally downwind of the site in the <100 km region, the largest population dose will occur to people living in the 1002000 km region. The air dispersion factors estimated in this study for both a northern and mid latitude site were comparable (within a factor of two), but smaller than the UNSCEAR dispersion factors. If radon decay (up to a factor of about five at 2000 km) were included in the UNSCEAR dispersion factors, the latter would also be comparable (within a factor of about three) to the dispersion factors used in this study.
Overall Uncertainty
The major objective of this study was to estimate the normalised (to unit of power produced) long term population dose due to radon emitted from uranium tailings, using the same methodology as UNSCEAR but based on more site-specific information. A quantitative uncertainty analysis of the results was beyond the scope of work for the study. However, the following qualitative comments are offered.
The contribution to the overall uncertainty due to imprecision in the site-specific information (overall radon emission rate, population densities, and the potential amount of uranium associated with the tailings) is likely to be much smaller than from the uncertainty in other factors that affect the estimation of population dose. Considering the uncertainties in population distributions within each region (<100 km and 1002000 km), meteorological dispersion, the air dispersion model and overall population growth, the overall uncertainty range for sites with specific information is subjectively estimated at about a factor of ten about the central estimate.
A reliable estimate of the population dose for sites not considered in this study (corresponding to about one-fifth of current world production) cannot be made at present because their site-specific conditions (radon emissions, population densities, etc.) were not available.
Other Issues
Perspective on the Estimated Doses
Notwithstanding the uncertainties associated with the analyses undertaken here, it is perhaps informative in terms of perspective to examine the magnitude of the estimated radon exposures and doses.
The total long term (10 000 year) population dose to radon emissions from the uranium tailings sites examined in this study is about 4650 person-Sv. Using the ICRP's nominal probability coefficient of 0.06 total cancers per person-Sv (Ref 18), this population dose converts to about 280 cancers over 10 000 years, or less than three cancers over a typical lifetime. This may be compared to the more than 60 million background cancers expected in the lifetime of the approximately 210 million people living within 2000 km of the sites (assuming a 30% background cancer incidence rate). The 4650 person-Sv estimate is a factor of 366 000 below the background population dose of 1.7 x 109 person-Sv for the same sites.
The area-weighted concentrations in the <100 km region are factors of 200 or more lower than background concentrations. A large fraction (3720 person-Sv or 80%) of the population dose is incurred by people living beyond 100 km of the sites. The area-weighted average radon concentrations in the 1002000 km region around the sites are estimated to range from near zero to about 0.2 mBq m-3. These concentrations are factors of more than 50 000 lower than typical outdoor background concentrations.
Linear No-Threshold Model
The population doses estimated in this study and by UNSCEAR implicitly assume the validity of the linear, no-threshold (LNT) dose response model; that is, the risks of exposure to radiation are assumed to be directly proportional to the dose received, down to zero dose. There is much current discussion of the appropriateness of the LNT model for estimating impacts from doses that are extremely small fractions of natural background radiation. The presence of a dose threshold, even a practical threshold in which competing causes of death defer the risk from radiation beyond the expected lifespan for detrimental effects, or of a hormetic effect, would render any assumed impacts associated with the population doses estimated in this study or by UNSCEAR invalid.
Integration Period
To be comparable to the UNSCEAR estimate, the population doses estimated in this study were integrated over a 10 000 year exposure period. While UNSCEAR chose this period for illustrative purposes, the assumption that conditions for which the dose estimates were derived would be constant over this period is clearly speculative.
However, it is quite likely that total or at least partial cures for some of the potential cancers associated with radiation exposure would have become available in the 10 000 year time period. On this basis, notwithstanding the issues associated with the LNT model noted previously, or from future improvements in medical care of cancer, it is likely the any impacts associated with the population doses derived in this study or by UNSCEAR will be overestimates of the actual (if any) impacts.
Conclusions
The major conclusions from the present study are:
- Based on the site-specific data available to this study, UNSCEAR's generic estimate of radon emissions from uranium tailings sites is too large for the sites examined in this study.
- Relative to the assumption of uniform (two region) population densities, the use of site-specific population distributions can result in significantly different (larger or smaller) population dose estimates.
- The radon concentrations associated with the tailings emissions are extremely small on both a relative (compared to typical background levels) and absolute (in terms of dose and risk) level. In the authors' view, the individual risk of cancer associated with the predicted concentrations is below a level that can be considered completely insignificant and trivial, i.e. de minimis.
- The uncertainties in the estimates provided in this study were reduced relative to the uncertainty in the UNSCEAR estimate because of the availability of more site-specific data on radon emissions and population densities. Further refinement of the analyses would require the use of site-specific meteorological data and population distributions.
- UNSCEAR's central estimate of the long term (10 000 year) population dose due to radon emissions from uranium mill tailings is too large, by about a factor of about 150, based on site-specific data and current and proposed tailings management practices at the sites examined in this study.