Preliminary Evaluation of Derived Concentration Guideline Level for Surface Soil at Wolsong NPP Site Using RESRAD-ONSITE Code

: A preliminary derived concentration guideline level (DCGL) was calculated for the soil of the Wolsong nuclear power plant (NPP) (the ﬁrst commercial pressurized heavy water reactor in South Korea) site using the RESRAD-ONSITE computational code. In total, fourteen selected radionuclides were analyzed after considering the preliminary evaluation information of radionuclides observed in the pressure tube specimen in the Wolsong Unit 1 heavy water reactor and previous NPP decommissioning cases. Furthermore, a geological structure model of the Wolsong NPP site was established according to the safety analysis report of the Wolsong NPP. In addition, the distribution coefﬁcients (Kd) of various radionuclides were derived using the JAEA-SDB and the pH information of the groundwater around the Wolsong NPP site. The DCGL for surface soil of the Wolsong NPP site was derived via the application of criteria for the site release to facilitate unrestricted reuse. Moreover, preliminary dose evaluation and relevant analysis were performed according to the Wolsong NPP site resident scenario. The novelty of this study lies in the ﬁrst calculation of the preliminary DCGL values for the case of the pressurized heavy water reactor (Wolsong NPP) site. It is expected that further reliable DCGL results might be achievable if more precise radionuclide information and site-speciﬁc parameters with respect to the Wolsong NPP site are secured and applied in the future.


Introduction
South Korea's first commercial pressurized heavy water reactor (PHWR), Wolsong Unit 1, was permanently shut down in 2019 in the absence of any record of decommissioning heavy water reactors worldwide.Currently, South Korea is exerting efforts in preparation for nuclear power plant (NPP) decommissioning through the establishment of the world's first heavy water reactor decommissioning institute.
Wolsong Unit 1 is a pressurized heavy-water reactor type (Canada Deuterium Uranium, CANDU) and, owing to the lack of research and results of decommissioning such NPPs, extensive prior research is necessary before the NPP decommissioning.The U.S. Multi-Agency Radiation Survey and Site Investigation Manual (MARSSIM) presents detailed standards for designing, planning, sampling, interpreting, and quality assuring for irradiation through consistent standard guidelines for site survey and instructs the inference of the derived concentration guideline level (DCGL) for site reuse [1].In the U.S., 10 CFR Part 20, Subpart E prescribes that the exposure dose to the receptor considering all possible exposure pathways should not exceed 0.25 mSv•year −1 based on the effective dose in case the business operator reuses the site in an unrestricted manner after decommissioning.In the case of South Korea, the Nuclear Safety and Security Commission Notification No. 2016-33 presents the standard of 0.1 mSv•year −1 .To reuse the site after NPP decommissioning, the DCGL must be derived from the evaluation of the exposure dose for the site using an evaluation model [2] based on a computer code such as RESRAD.The RESRAD computational code, developed by the Argonne National Laboratory (ANL) in the U.S., is designed to evaluate the risk and exposure dose of residual radioactive materials and is used for the decommissioning of commercial NPPs in the U.S. under the approval of the U.S. Nuclear Regulatory Commission (NRC).The Rancho Seco [3] and Zion [4] NPPs in the U.S. are typical examples that applied the RESRAD computer code to NPP decommissioning.The Rancho Seco NPP was a pressurized water reactor (PWR) plant in California, USA which started commercial operation in 1975 and has now been fully decommissioned.To determine the list of potential radionuclides in this plant decommissioning case, a list of nuclides reflecting site characteristics was derived by referring to literature such as NUREG/CR-0130 [5], NUREG/CR-3474 [6], and NUREG/CR-4289 [7].Furthermore, the DCGL criteria for surface soil, building, and materials were established via the application of the RESRAD computer code [3].Meanwhile, the Zion NPP was another PWR type reactor that started commercial operation in 1973 and was permanently shut down in 1998 after 15 years.The information regarding the NUREG/CR-3474 [6], 4289 [7], WINCO-1191 [8], and site waste specimens needed to derive a list of site-specific nuclides and the DCGL criteria for site decommissioning were derived by using the RESRAD code [4].
The objective of this study was to obtain a preliminary DCGL for the soil of the Wolsong NPP site, which is a heavy water reactor, using the RESRAD-ONSITE computational code.The target radionuclides were selected by referring to previous studies on NPP decommissioning, and a geological structure model of the Wolsong NPP site was established.Further, the distribution coefficient reflecting the site characteristics was derived using the JAEA-SDB [9] and the pH information of groundwater around the Wolsong NPP site.The preliminary DCGL, satisfying the limit of the unrestricted reuse criteria, was derived based on the various site-specific input parameters as input data of the RESRAD-ONSITE code.Furthermore, the uncertainty of the input parameters for the result of exposure dose was quantified using the probabilistic analysis feature of the RESRAD-ONSITE.In addition, the effects of the variables used in the calculation of the RESRAD code on the result of exposure dose were analyzed by performing a sensitivity analysis.

RESRAD-ONSITE Computational Code
The RESRAD (RESidual RADioactivity) computer code [10] was developed in 1989 by the ANL in the U.S. and is used in various countries, institutions, companies, and research groups.It is used to evaluate the radiation exposure of humans, animals, and plants in the ecosystem by residual radioactive materials and the resulting potential exposure dose and risk.In particular, the RESRAD-ONSITE code is used to calculate the radiation exposure of receptors in contaminated sites.In the U.S., the RESRAD-ONSITE code is widely used to satisfy the site opening criteria considering both social/economic factors and technical validity by the Nuclear Regulatory Commission (NRC) and the Department of Energy (DOE) when decommissioning an NPP.The primary exposure pathways of the RESRAD-ONSITE include external radiation, inhalation, and ingestion.The exposure dose by external radiation is calculated by the dose to which the receptor was exposed at 1 m from the ground and via the application of correction factors such as area, thickness, shape, and shielding of the contaminated zone.Further, the exposure pathways by inhalation primarily include the inhalation of radon decay products and contaminated dust.Finally, the exposure pathways by ingestion include exposure of receptors to radioactive substances in food, water, and soil.

Target Radionuclides
The radionuclides to be analyzed were selected using the radionuclide selection method, applied to the decommissioning of the Rancho Seco [3] and Zion [4] NPPs in the U.S., for 37 radionuclides detected in the Wolsong Unit 1 pressure tube specimens from a study on the evaluation of heavy water reactor source terms [11].
As a CANDU-type heavy water reactor, the Wolsong Unit 1 is composed of 380 nuclear fuel channels and each nuclear fuel channel is comprised of a pressure tube and calandria tube.The calandria tube acts as the pressure and thermal shielding boundary between the heavy water moderator and the coolant channel, whereas the pressure tube supports the nuclear fuel at the pressure boundary and nuclear reactor core to maintain the pressure so that the coolant does not boil.The pressure tube is made of Zr alloy containing 2.5% of Nb, with a total length, inner radius, outer radius, and density of 6.4 m, 5.2 cm, 5.6 cm, and 6.6 g•cm −3 , respectively [11].According to the Wolsong Unit 1 source term evaluation report, 37 radionuclides were derived from the preliminary evaluation of the radioactivity of the pressure tube using the ORIGEN2 code for 9 pressure tube specimens of the Wolsong Unit 1 [11].When considering the License Termination Plan (LTP) of the Rancho Seco power plant, nuclides with a half-life of fewer than 2 years were excluded from the target radionuclides because their effect on the result was considered insignificant [3].Thus, based on ICRP 107 [12], among the 37 radionuclides derived above from the pressure tube of the Wolsong Unit 1, 10 radionuclides were found to exhibit a half-life of fewer than 2 years ( 54 Mn, 90 Y, 108 Ag, 109 Cd, 109m Ag, 125m Te, 137m Ba, 182 Ta, 192 Ir, and 192m Ir), while the remaining 27 radionuclides exhibited a half-life of 2 years or longer.Therefore, this study also excluded radionuclides with a half-life of fewer than 2 years from the evaluation by applying the sample case of the Rancho Seco LTP.
Further, a study related to the decommissioning of the Zion NPP in the U.S. reported a list of relative fractions of radionuclides based on the attenuation-corrected radioactivity concentration for the list of nuclides detected in 19 samples collected from the Zion NPP site [4].Consequently, on analyzing the 27 target radionuclides while referring to the literature [4] related to the decommissioning of the Zion NPP and the result of the preliminary evaluation of the radioactivity of pressure tube specimens from the Wolsong Unit 1, 12 radionuclides ( 3 H, 14 C, 55 Fe, 60 Co, 63 Ni, 90 Sr, 93 Zr, 93m Nb, 94 Nb, 108m Ag, 121m Sn, and 125 Sb) were found to exhibit a half-life of 2 years or longer with a relative fraction of 0.01% or higher.The sum of their relative fractions was ca.99.42%.In contrast, another 15 radionuclides yielded a relative fraction of less than 0.01% with the sum of their relative fractions being ca.0.02%.Furthermore, the sum of relative fractions of the other 10 radionuclides with a half-life of fewer than 2 years was calculated to be ca.0.56% [11].In this study, the 25 radionuclides with a relative fraction of less than 0.01% were excluded from evaluation owing to the small sum of their fractions, which exerts a negligible effect on the result.Furthermore, in addition to the 12 radionuclides for analysis, the radionuclides 134 Cs and 137 Cs, which were detected in the previous study on the Kori Unit 1, Rancho Seco, and Zion NPPs, were included; thus, 14 radionuclides were finally selected as the target of analysis, and they are listed in Table 1 [3,4,13].

Geological Structure Model
The Wolsong NPP site is located on the southeastern coast of South Korea and the geological structure of the site comprises sediments of the Cretaceous rocks strongly affected by volcanic activities, volcanic sediments, faults, or intrusions [14].The rocks are mainly composed of argillite, quartz andesite, quartz diorite, granitic rocks, etc. Further, the thicknesses of the sedimentary layers distributed in the Wolsong NPP site have been reported to be in the range of 1-5 m (Wolsong unit 2), 6-15 m (Wolsong unit 3), and 6-12 m (Wolsong unit 4) [14].In the present work, the thickness of the contaminated zone was assumed as 0.15 m by referring to the case of the Rancho Seco NPP, and the thickness of the unsaturated zone was set to be 1-15 m, considering the thickness of the sedimentary layer of Wolsong NPP site (see Figure 1).Furthermore, the geological features of the contaminated and unsaturated zones were assumed to be sand, soil, loam, etc., whereas that of the saturated zone was assumed to be granitic rocks, based on the typical geological features of South Korea [15].

Geological Structure Model
The Wolsong NPP site is located on the southeastern coast of South Korea and the geological structure of the site comprises sediments of the Cretaceous rocks strongly affected by volcanic activities, volcanic sediments, faults, or intrusions [14].The rocks are mainly composed of argillite, quartz andesite, quartz diorite, granitic rocks, etc. Further, the thicknesses of the sedimentary layers distributed in the Wolsong NPP site have been reported to be in the range of 1-5 m (Wolsong unit 2), 6-15 m (Wolsong unit 3), and 6-12 m (Wolsong unit 4) [14].In the present work, the thickness of the contaminated zone was assumed as 0.15 m by referring to the case of the Rancho Seco NPP, and the thickness of the unsaturated zone was set to be 1-15 m, considering the thickness of the sedimentary layer of Wolsong NPP site (see Figure 1).Furthermore, the geological features of the contaminated and unsaturated zones were assumed to be sand, soil, loam, etc., whereas that of the saturated zone was assumed to be granitic rocks, based on the typical geological features of South Korea [15].

Distribution Coefficient
The Japan Atomic Energy Agency (JAEA) provides the JAEA sorption database (JAEA-SDB) for utilization in the disposal safety evaluation modeling of radioactive waste [9].The distribution coefficient (Kd) of radionuclide is a major variable used to predict the migration and retardation behavior of radionuclides in the groundwater system.The Kd values used in this study were selected from the JAEA-SDB considering the geological structure model of the Wolsong NPP site that was assumed above and the mean pH value of the groundwater around the Wolsong NPP site.The pH values of groundwater were obtained from the National Groundwater Information Center (https://www.gims.go.kr;

Distribution Coefficient
The Japan Atomic Energy Agency (JAEA) provides the JAEA sorption database (JAEA-SDB) for utilization in the disposal safety evaluation modeling of radioactive waste [9].The distribution coefficient (Kd) of radionuclide is a major variable used to predict the migration and retardation behavior of radionuclides in the groundwater system.The Kd values used in this study were selected from the JAEA-SDB considering the geological structure model of the Wolsong NPP site that was assumed above and the mean pH value of the groundwater around the Wolsong NPP site.The pH values of groundwater were obtained from the National Groundwater Information Center (https://www.gims.go.kr; accessed on 31 August 2021) which provides comprehensive groundwater information on South Korea.Consequently, the pH of the natural groundwater of Gyeongju City, which is near the Wolsong NPP site, 5.83-8.40,was used to derive the Kd values.
Table 2 shows the Kd values for 14 target radionuclides derived from the JAEA-SDB.Since the JAEA-SDB does not provide the Kd value of Ag for granitic rock that satisfies the above pH condition, the present work replaced the Kd value of Ag with that of Cs for granitic rock by referring to the literature [16].a : Kd value of Cs is applied based on the chemical analogy; b : deterministic value.
Further, there is a lack of data regarding Fe and C for granitic rock in the JAEA-SDB.Thus, it was assumed that the Kd values of Fe for unsaturated zone and granitic rock were identical to each other.Whereas for C, the Kd value suggested in the literature [16] was selectively applied; however, a deterministic single value was input, instead of a distribution function.

Input Parameters
The input parameters database of the RESRAD-ONSITE code used in the present work for the Wolsong NPP was established by referring to overseas NPP decommissioning cases [3,4], previous studies [17][18][19], and the site-specific variables of the Wolsong NPP site [10,14,[20][21][22][23][24][25][26].Consequently, based on the geological structure model of the Wolsong NPP site that was assumed above, the hydrological input parameters for each zone were taken from the NUREG/CR-6697 [17].The area of the contaminated zone, which is a characteristic variable of the Wolsong NPP site, was applied based on a single unit of the Wolsong NPP, and the thickness of the contaminated zone was set to 0.15 m by referring to the case of the Rancho Seco NPP [3,14].Furthermore, with respect to the density of the saturated zone, the averaged data for quartz andesite, quartz diorite, and andesite, reported as the main bedrocks in the Wolsong Final Safety Analysis Report (FSAR), were employed.The hydraulic conductivity was input by referring to the data of a previous study on field permeability tests in the Wolsong NPP site [14,20].Moreover, for the water transport model, the non-dispersion model proposed by default in the RESRAD-ONSITE code for sites with a contaminated zone larger than 1000 m 2 was applied [10].Furthermore, the input parameters related to diet in a contaminated site, assuming resident exposure scenario, were input by referring to the domestic annual consumption reports of South Korea [21][22][23][24][25].In addition, regarding the precipitation and wind speed, the mean values of climate data measured in Gyeongju city from 2011 to 2020 were employed [26].The detailed input parameters used in this study are listed in Tables 3 and A1 in Appendix A. Table A1 presents various input parameters used in the RESRAD-ONSITE code taken from various research works.The data obtained from NUREG/CR-7267 [19] were alternatively used in case of the absence of site-specific parameters available for the Wolsong NPP site.

Probabilistic Analysis
The RESRAD-ONSITE code provides probabilistic analysis features for calculating the results from the input parameters that contain uncertainty.Following various parameters being input as distribution functions, the correlations of the results derived from them can be quantified.Subsequently, the assigned values for specific parameters can be derived and reflected in the result.
In this study, the non-linear monotonic relationship between the input parameters and the exposure dose results was estimated by referring to the LTP case of the Rancho Seco NPP.In addition, the partial rank correlation coefficient (PRCC), which provides a unique contribution to the exposure dose result, was used.If the absolute value of the PRCC exceeded 0.25, the corresponding input parameter was treated as a sensitive parameter.The uncertainty of the input parameter can be quantified by replacing the deterministic value corresponding to the 75% quartile if the PRCC value is positive, and the 25% quartile if it is negative, with the input parameter.The list of sensitive parameters for which the absolute value of PRCC exceeded 0.25 in this study is listed in Table A2 in Appendix A. Consequently, the deterministic value for each radionuclide was derived and the preliminary DCGL was calculated.

Results
NUREG/CR-1757 recommends the use of the peak of the mean dose to derive the DCGL value.According to the previous studies [2,13,27], the DCGL formula can be expressed as follows: Table 4 shows the preliminary DCGL results derived with the assigned value for each radionuclide through probabilistic analysis (see Table A2).As presented in Table A2, the highest PRCC values were mostly due to the external gamma shielding factor, density of contaminated zone, runoff coefficient, distribution coefficients, evapotranspiration coefficient, etc.The PRCC and assigned values will be further discussed in the following chapter.The regulatory dose limit applied in this study was 0.25 mSv•year −1 according to the 10 CFR Part 20, subpart E in the U.S., and 0.1 mSv•year −1 according to the Nuclear Safety and Security Commission Notification No. 2016-33 in South Korea.On deriving the preliminary DCGLs for 14 target radionuclides, the radionuclide 93m Nb yielded the largest DCGL value, followed by 55 Fe, 121m Sn, and 63 Ni.Those large DCGL values of several radionuclides were seemingly caused by their large Kd values (e.g., log Kd saturated zone = 2.93 ± 0.48 for Nb, 3.52 ± 0.70 for Sn, 2.82 ± 0.65 for Ni, respectively) inducing remarkable retardation of migration along with relatively short half-life (e.g., t 1/2 = 2.7 years for 55 Fe and 16 years for 93m Nb).
The result of the preliminary DCGL of the Wolsong NPP site soil derived in this study was compared with other studies [13,[28][29][30][31] as presented in Table 5.The authors note that the regulatory dose limit was unified to the U.S. criterion of 0.25 mSv•year −1 to establish the same criterion for the comparison among various studies.The comparisons indicated relatively similar results for most of the radionuclides.In particular, the results of major radionuclides, such as 14 C, 60 Co, 134 Cs, 137 Cs, 63 Ni, and 90 Sr, were relatively consistent with the values derived from Kori-1 [13].Furthermore, there were no data in the other decommissioning cases for comparing the DCGL values for 93m Nb and 93 Zr.However, the radionuclide 93m Nb showed a significantly higher DCGL value compared to the IAEA standard.As explained above, this is considered to be the effect of the considerably high Kd value of 93m Nb, hindering the spread of the contamination, together with a relatively short half-life.However, although the Kd values of 93m Nb and 93 Zr are similar to each other, the DCGL of 93 Zr obtained in the present work indicated a relatively decreased tendency in comparison to the IAEA standard.It might be caused by a considerably longer half-life of 93 Zr (t 1/2 = 1.53 × 10 6 years) compared with that of 93m Nb (t 1/2 = 16 years).A longer half-life means a larger number of radionuclides required to provide the same radioactivity.Due to the limited retention capacity of contaminated, unsaturated, and saturated zones, the DCGL of 93 Zr cannot be as high as 93m Nb.Nevertheless, it is expected that a more reliable DCGL value of various radionuclides could be derived by using the Kd value directly obtained from the practical adsorption test for Wolsong NPP site soil in the future.

Discussion
The sensitivity parameters and assigned values were calculated with PRCC based on various site-specific input parameters listed in Table A1 in Appendix A. For preliminary calculation of the timing for site release for unrestricted reuse after NPP decommissioning, the total exposure dose was derived by assuming the concentration in the contaminated soil of various radionuclides to 0.037 Bq•g −1 (1 pCi•g −1 ), the same initial condition employed in the previous works [13,32], and is presented in Figure 2. Consequently, the timing when the regulatory limit of 0.1 mSv•year −1 (allowing unrestricted reuse of the site after decommissioning the NPP in South Korea according to Article 4 of the Nuclear Safety and Security Commission Notification No. 2016-33) is satisfied was derived to be 4.07 years after the NPP decommissioning and a final investigation of the site condition.
missioning the NPP in South Korea according to Article 4 of the Nuclear Safety and Security Commission Notification No. 2016-33) is satisfied was derived to be 4.07 years after the NPP decommissioning and a final investigation of the site condition.In addition, this study evaluated the effect of a change of input parameters on the final exposure dose result by using the sensitivity analysis feature of the RESRAD-ON-SITE code.To that end, the sensitivity parameters for the 14 radionuclides were identified by applying the PRCC.Consequently, eight representative sensitivity parameters were derived as listed in Table 6.The sensitivity parameter that yielded the highest PRCC was the external gamma shielding factor, followed by the density of the contaminated zone, runoff coefficient, distribution coefficient (Kd), and evapotranspiration coefficient.Moreover, all the eight sensitivity parameters had positive correlations, and the deterministic values corresponding to the 75% quartile were derived accordingly.Furthermore, an analysis to evaluate the relative effect of changes in individual sensitivity parameters on the unrestricted site reuse criteria was also performed.To that end, the change of the site reuse criteria according to the change of input parameters was evaluated while increasing and decreasing the values of each sensitivity parameter by a factor of two.However, the input values for runoff coefficient and evapotranspiration coefficient among the sensitivity parameters were limited to the range of 0 and 1.Thus, the effect on In addition, this study evaluated the effect of a change of input parameters on the final exposure dose result by using the sensitivity analysis feature of the RESRAD-ONSITE code.To that end, the sensitivity parameters for the 14 radionuclides were identified by applying the PRCC.Consequently, eight representative sensitivity parameters were derived as listed in Table 6.The sensitivity parameter that yielded the highest PRCC was the external gamma shielding factor, followed by the density of the contaminated zone, runoff coefficient, distribution coefficient (Kd), and evapotranspiration coefficient.Moreover, all the eight sensitivity parameters had positive correlations, and the deterministic values corresponding to the 75% quartile were derived accordingly.Furthermore, an analysis to evaluate the relative effect of changes in individual sensitivity parameters on the unrestricted site reuse criteria was also performed.To that end, the change of the site reuse criteria according to the change of input parameters was evaluated while increasing and decreasing the values of each sensitivity parameter by a factor of two.However, the input values for runoff coefficient and evapotranspiration coefficient among the sensitivity parameters were limited to the range of 0 and 1.Thus, the effect on site reuse criteria was evaluated while increasing or decreasing the runoff coefficient by 60% and the evapotranspiration coefficient by 45%.Table 7 summarizes the analysis results for all the eight sensitivity parameters and the relevant changes in the criteria for site release for unrestricted reuse.The result of sensitivity analysis showed that the parameters that have a major impact on the year for site release to facilitate unrestricted reuse were the density of contaminated zone, evapotranspiration coefficient, and external gamma shielding factor.As shown in Figure 3a, if the input for the density of the contaminated zone through sensitivity analysis is doubled, the year for site release for unrestricted reuse was increased to 11.2 years, which is 2.8 times longer than the standard value of 4.07 years, because of the distinct increase in the total exposure dose.This tendency might be caused by the largest impact of external radiation among various exposure pathways on the total exposure dose, together with the increased number of radionuclides contained in the soil due to the increased density of the contaminated zone (see Figure 3b).
Furthermore, the groundwater contamination in the NPP site was mainly caused by the leaching of radionuclides from the water infiltration.This is simulated in the RESRAD-ONSITE code as follows [10]: where L i = leach rate for radionuclide i (year −1 ) R(t) = source strength at time t (pCi•year −1 ) ρ = bulk density of the contaminated zone (kg•m −3 ) A = area of the contaminated zone (m 2 ) T(t) = thickness of the contaminated zone at time t (m) S i (t) = average concentration of the i-th principal radionuclide in the contaminated zone available for leaching at time t (Bq•kg −1 ) zone available for leaching at time t (Bq•kg ) According to Equation ( 2), the leach rate (Li) has an inverse relationship wit density of the contaminated zone (ρ).In other words, with an increase in the dens the contaminated zone, the leaching of radionuclides to groundwater decreases relat Thus, the radionuclides stay longer in the contaminated zone, which increases the d external exposure to receptors, and consequently, the total exposure dose increas well.Further, the relationship between the evapotranspiration coefficient and the in tion rate is expressed in the RESRAD-ONSITE code as follows [10]: Here, if the evapotranspiration coefficient increases, the infiltration rate of the tration water decreases relatively.Thus, similar to the case of the density of contami zone, the transport of radionuclides is relatively delayed.Consequently, there is a crease in the external exposure to receptors which increases the total exposure dose results of this study showed that when the evapotranspiration coefficient was incr by 1.45 times, the year for site release for unrestricted reuse increased to approxima According to Equation ( 2), the leach rate (L i ) has an inverse relationship with the density of the contaminated zone (ρ).In other words, with an increase in the density of the contaminated zone, the leaching of radionuclides to groundwater decreases relatively.Thus, the radionuclides stay longer in the contaminated zone, which increases the direct external exposure to receptors, and consequently, the total exposure dose increases as well.
Further, the relationship between the evapotranspiration coefficient and the infiltration rate is expressed in the RESRAD-ONSITE code as follows [10]: where Here, if the evapotranspiration coefficient increases, the infiltration rate of the infiltration water decreases relatively.Thus, similar to the case of the density of contaminated zone, the transport of radionuclides is relatively delayed.Consequently, there is an increase in the external exposure to receptors which increases the total exposure dose.The results of this study showed that when the evapotranspiration coefficient was increased by 1.45 times, the year for site release for unrestricted reuse increased to approximately 8 years from the reference case of 4.07 years because of a significant increase in the total exposure dose.
The external gamma shielding factor is also considered to be a major factor that significantly affects the evaluation of the receptors' radiation exposure to radionuclides remaining in the contaminated zone.Among the parameters related to external radiation, which is a major exposure path, the external gamma shielding factor is used to derive the environmental transport factor (ETF), which is expressed as follows: where ETF(t) = environment transport factor (dimensionless) FO = occupancy and shielding factor (dimensionless) FS(t) = shape factor (dimensionless) FA(t) = nuclide-specific area factor (dimensionless) FCD(t) = depth-and-cover factor (dimensionless) f otd = fraction of a year spent outdoors, on site (dimensionless) f ind = fraction of a year spent indoors, on site (dimensionless) F sh = indoor shielding factor for external gamma (dimensionless) Thus, according to Equations ( 4) and ( 5), the environmental transport factor (ETF) increases together with the external gamma shielding factor.Consequently, the exposure dose that receptors receive via the external radiation pathway also increases.Therefore, it can be evaluated as a major factor that increases the year for site release for unrestricted reuse.

Conclusions
The preliminary DCGLs for the soil of the Wolsong NPP site, which is a pressurized heavy water reactor, were derived in this study according to the methodology suggested by the MARSSIM.The decommissioning cases of overseas NPPs such as Rancho Seco and Zion NPPs were considered, and various literature was referred to as well for the evaluation.The target radionuclides for calculating the DCGL were selected based on the data of radionuclides identified in the pressure tube specimens of Wolsong Unit 1 that have been reported in the literature.The input parameters were selected considering the site-specific variables of the Wolsong NPP site and employed in the RESRAD-ONSITE code.In addition, input data such as distribution coefficients were derived and applied to simulate the migration and retardation behavior of radionuclides using the JAEA-SDB [9].Consequently, the DCGLs for the Wolsong NPP site soil with respect to the resident scenario were preliminarily derived by using the input parameter database built in the present work.A comparison of the results with that of previous studies showed that they were relatively consistent with each other.
To preliminarily derive the year for site release for unrestricted reuse, the inventory of 14 target radionuclides was set to 0.037 Bq•g −1 (1 pCi•g −1 ) in the contaminated zone, and their total exposure dose was calculated accordingly.The result showed that the exposure dose of 0.1 mSv•year −1 for site reuse was satisfied at 4.07 years after the decommissioning of the Wolsong NPP and the final investigation of the site condition.Furthermore, as a result of sensitivity analysis to identify input parameters that affect the calculation of the year for site release for unrestricted reuse, the density of contaminated zone, evapotranspiration coefficient, and external gamma shielding factor were found to be major sensitivity parameters.These parameters were evaluated to have a significant impact on the external radiation pathway by radionuclide numbers remaining in the contaminated zone.
As the first study that evaluated the DCGL for the pressurized heavy water reactor, Wolsong NPP, the input parameter database and calculation model developed in the present work are expected to be useful as a preliminary study that can be referenced when decommissioning heavy water reactors in South Korea in the future.Furthermore, it is foreseen that the economic and efficient site release standards for pressurized heavy water reactors can be established through the derivation of further reliable exposure dose results if the uncertainties in radionuclide inventory, hydrogeological parameters, and other input parameters are reduced through practical experiments on soil samples in the NPP site, along with additional field surveys in the future.For example, previous research work has confirmed the considerable relationship between the distribution coefficient and exposure dose results in the safety assessment model of nuclear waste repositories [33].In the former investigation, a relatively increased distribution coefficient provided a remarkably delayed release rate of radionuclides owing to the enhanced retention capacity of geologic barriers and presented a relatively delayed peak time and decreased peak value of exposure dose.Since the distribution coefficient is an element-specific and conditional parameter, depending on various hydro-geochemical conditions, securing more precise distribution coefficient data for various radionuclides through a practical field test is expected to provide a significant contribution to further reliable model calculation.

Figure 1 .
Figure 1.Geological structure model of the Wolsong NPP site.

Figure 1 .
Figure 1.Geological structure model of the Wolsong NPP site.

Figure 2 .
Figure 2. Total exposure dose as a function of time after the NPP decommissioning calculated with RESRAD-ONSITE code for Wolsong NPP site.

Figure 2 .
Figure 2. Total exposure dose as a function of time after the NPP decommissioning calculated with RESRAD-ONSITE code for Wolsong NPP site.

Figure 3 .
Figure 3. (a) Sensitivity analysis result for the density of contaminated zone and (b) relevant sure pathways for the reference case.Minor pathways are not shown.

Figure 3 .
Figure 3. (a) Sensitivity analysis result for the density of contaminated zone and (b) relevant exposure pathways for the reference case.Minor pathways are not shown.

Table 1 .
Target radionuclides for the derivation of the preliminary DCGL.

Table 2 .
Distribution coefficients of various radionuclides.

Table 3 .
List of input parameters related to the characteristics of the Wolsong NPP site.

Table 4 .
DCGL results for various radionuclides were calculated with two different regulatory dose limits.

Table 5 .
Comparison of DCGL results (Bq•g −1 ) for various radionuclides derived with a regulatory dose limit of 0.25 mSv•year −1 .
a : Connecticut Yankee.

Table 6 .
Input parameters derived from the sensitivity analysis.

Table 6 .
Input parameters derived from the sensitivity analysis.

Table 7 .
Year for site release for the unrestricted reuse of Wolsong NPP site calculated based on the sensitivity analysis for various input parameters.
a : Total exposure dose is less than 0.1 mSv•year −1 at t = 0 year.

Table A2 .
Partial rank correlation coefficient (PRCC) and assigned value for each radionuclide derived in this study.