# Estimation of Intervention Distances for Urgent Protective Actions Using Comparative Approach of MACCS and InterRAS

Authors: Mazzammal Hussain; Salah Ud-Din Khan; Waqar A. Adil Syed; Shahab Ud-Din Khan
Journal: Science and Technology of Nuclear Installations
 (2014)

---

## Abstract

Distances for taking evacuation as a protective measure during early phase of a nuclear accident have been approximated using MELCOR Accident Consequence Code System (MACCS). As a reference data, the source term of Pakistan Research Reactor 1 (PARR-1) and meteorological data of Islamabad, Pakistan, have been considered. Based on comparison with published data and international radiological assessment (InterRAS) code results, it is concluded that MACCS is a rational tool for estimation of urgent protective actions during early phase of nuclear accident by taking into account the variations in meteorological and release concentrations parameters.


---

## Body

## 1. Introduction

A range of probable accidents are associated with nuclear power plants starting with minor incident to immense disaster. Different researchers had used MACCS code to analyze the nuclear power plants accident consequences. Jongtae Jeong and Wondea Jung studied estimation of early health effects for different combinations of release parameters and meteorological data using MACCS code for Younggwang 3 and 4 nuclear power plants in Korea concluding that, with the same amount of radioactive material released to the atmosphere, a large difference in early health effects from case to case was observed [1]. Haste and coworkers attempted to demonstrate a MELCOR-MACCS capability to simulate the accident scenarios for the whole plant, including the containment response and off-site consequences arising from fission product release from the containment. The results provided a good basis for the NPP analysis foreseen [2]. Ke-Shih performed reevaluation of emergency planning zones (EPZs) for nuclear power plants (NPPs) using MACCS2 code and concluded that the radius identified previously is a reasonable conservative value of EPZs for each of the three operating NPPs in Taiwan [3]. Jeong and Ha studied influence of source term release parameters on health effects for Younggwang 3 and 4 nuclear power plants in Korea using MACCS code and concluded that the research work will be very useful for developing strategies for reducing off-site consequences of accident management if they are combined with influence of weather conditions [4]. Thaning and Baklanov considered a simulated accident at a nuclear power plant that could cause a large release of radioactivity into the atmosphere. The consequence analysis was performed using a 3-dimensional mesoscale model and MACCS code [5].


## 2. Research Methodology

### 2.1. Computational Technique: MACCS Code

Sandia National Laboratories developed MELCOR Accident Consequence Code System (MACCS) for the consequence assessment of severe accidents at nuclear power plants. MACCS is organized into three modules, the ATMOS module which performs the atmospheric transport and deposition portion of the calculation, the EARLY module which estimates the consequences of the accident immediately after the accident usually within the first week, and the CHRONC module which estimates the long term consequences of the accident. The flow diagram of MACCS code [6–8] is shown in Figure 1.Figure 1
Flow diagram of MACCS code.In MACCS code, Gaussian plume model has been used for the study of atmospheric dispersion of radioactive material and vertical and crosswind distributions. The plume dimensions are defined in vertical and crosswind direction by the standard deviations (σ
y
,
σ
z) of the normal distribution of material concentration in crosswind and vertical directions. The general form of the Gaussian plume equation is
(1)
χ
(
x
,
y
,
z
)
=
Q
2
π
u
-
σ
y
σ
z
×
exp
⁡
⁡
[
-
1
2
(
y
σ
y
)
2
]
exp
⁡
[
-
1
2
(
z
-
h
σ
z
)
2
]
,
where χ
(
x
,
y
,
z
) is the time integrated air concentration (Bq-s/m3) at the downwind location (
x
,
y
,
z
), Q is the source strength (Bq), u
- is the mean downwind speed, σ
y
and
σ
z are the standard deviations (m) of the normal concentration distribution along crosswind and vertical direction, and h is the release height (m).Equation (1) is not applicable when plume expands vertically and is bounded by mixing layer or by the ground. To solve this problem, ground and mixing layers are considered as totally reflecting boundaries. This is achieved by adding a mirror image source below the ground and above the inversion layer. By considering this effect in (1), the centerline air concentration χ
(
x
=
0
,
y
=
0
,
z
=
H
) and ground concentrations χ
(
x
=
0
,
y
=
0
,
z
=
0
), after time of release to the time at which the concentrations become uniform along vertical direction, are given by
(2)
χ
(
x
,
y
=
0
,
z
)
=
Q
2
π
σ
y
σ
z
u
-
[
exp
⁡
[
-
1
2
(
z
-
H
σ
z
)
2
]
+
exp
⁡
[
-
1
2
(
z
+
H
σ
z
)
2
]
+
∑
n
=
1
5
{
exp
⁡
[
-
1
2
(
z
-
H
-
2
n
L
σ
z
)
2
]
+
exp
⁡
[
-
1
2
(
z
+
H
-
2
n
L
σ
z
)
2
]
+
exp
⁡
[
-
1
2
(
z
-
H
+
2
n
L
σ
z
)
2
]
+
exp
⁡
[
-
1
2
(
z
+
H
+
2
n
L
σ
z
)
2
]
}
]
,
where H
=
(
h
+
Δ
h
) is the height of the plume centerline (m),h = the initial release height of the plume before plume rise (m), Δ
h = the amount of plume rise (m), and L = the height (m) of the inversion layer (mixing height).In MACCS, only first five terms are considered and the rest of the terms are neglected. When a uniform vertical distribution is attained, the following equation is used to calculate centerline air concentration:(3)
χ
(
x
,
y
=
0
,
z
)
=
Q
2
π
u
-
σ
y
L
.
For the estimation of radiation doses during early phase of nuclear emergency, the following pathways are taken into consideration:(i)
external and internal exposure from cloud shine,
(ii)
exposures from ground shine,
(iii)
internal exposure from resuspension inhalation and skin doses from deposited material onto the skin,
(iv)
acute and lifetime doses from early exposures.


### 2.2. Meteorological Conditions

For the execution of MACCS code, the hourly-based data for the year 2010 of Islamabad collected from the ground station (with an approximate height of 10 meters) has been used. As the near surface wind speed increases with altitude, the same data if used for the stake height (61 m) will overestimate the plume rise. This could produce significant underestimation of the radiation doses. To incorporate the effect of wind speed (m/s) with altitude, the following theoretical formula was used to estimate the wind speed at higher altitudes:(4)
u
=
u
o
(
z
z
o
)
p
,
where u = wind speed at height z (m/s), u
o = wind speed at surface (m/s), and p = parameter (dimensionless) that varies with stability class and surface roughness.The wind rose for the reported year is presented in Figure2. The percentage value of wind direction remained 3.7% in north, 8.7% in northeast, 8.4% in east, 13.8% in southeast, 10.8% in south, 38.3% in southwest, 6.4% in west, and 10% in northwest, respectively. Data presented in Figure 2 represents that dominant wind direction was southwest (38.3%) during year 2010. The dominant atmospheric stability class as presented in Figure 3 was stability class-F.Figure 2
Wind rose for sixteen directions.Figure 3
Dominant stability class.


## 2.1. Computational Technique: MACCS Code

Sandia National Laboratories developed MELCOR Accident Consequence Code System (MACCS) for the consequence assessment of severe accidents at nuclear power plants. MACCS is organized into three modules, the ATMOS module which performs the atmospheric transport and deposition portion of the calculation, the EARLY module which estimates the consequences of the accident immediately after the accident usually within the first week, and the CHRONC module which estimates the long term consequences of the accident. The flow diagram of MACCS code [6–8] is shown in Figure 1.Figure 1
Flow diagram of MACCS code.In MACCS code, Gaussian plume model has been used for the study of atmospheric dispersion of radioactive material and vertical and crosswind distributions. The plume dimensions are defined in vertical and crosswind direction by the standard deviations (σ
y
,
σ
z) of the normal distribution of material concentration in crosswind and vertical directions. The general form of the Gaussian plume equation is
(1)
χ
(
x
,
y
,
z
)
=
Q
2
π
u
-
σ
y
σ
z
×
exp
⁡
⁡
[
-
1
2
(
y
σ
y
)
2
]
exp
⁡
[
-
1
2
(
z
-
h
σ
z
)
2
]
,
where χ
(
x
,
y
,
z
) is the time integrated air concentration (Bq-s/m3) at the downwind location (
x
,
y
,
z
), Q is the source strength (Bq), u
- is the mean downwind speed, σ
y
and
σ
z are the standard deviations (m) of the normal concentration distribution along crosswind and vertical direction, and h is the release height (m).Equation (1) is not applicable when plume expands vertically and is bounded by mixing layer or by the ground. To solve this problem, ground and mixing layers are considered as totally reflecting boundaries. This is achieved by adding a mirror image source below the ground and above the inversion layer. By considering this effect in (1), the centerline air concentration χ
(
x
=
0
,
y
=
0
,
z
=
H
) and ground concentrations χ
(
x
=
0
,
y
=
0
,
z
=
0
), after time of release to the time at which the concentrations become uniform along vertical direction, are given by
(2)
χ
(
x
,
y
=
0
,
z
)
=
Q
2
π
σ
y
σ
z
u
-
[
exp
⁡
[
-
1
2
(
z
-
H
σ
z
)
2
]
+
exp
⁡
[
-
1
2
(
z
+
H
σ
z
)
2
]
+
∑
n
=
1
5
{
exp
⁡
[
-
1
2
(
z
-
H
-
2
n
L
σ
z
)
2
]
+
exp
⁡
[
-
1
2
(
z
+
H
-
2
n
L
σ
z
)
2
]
+
exp
⁡
[
-
1
2
(
z
-
H
+
2
n
L
σ
z
)
2
]
+
exp
⁡
[
-
1
2
(
z
+
H
+
2
n
L
σ
z
)
2
]
}
]
,
where H
=
(
h
+
Δ
h
) is the height of the plume centerline (m),h = the initial release height of the plume before plume rise (m), Δ
h = the amount of plume rise (m), and L = the height (m) of the inversion layer (mixing height).In MACCS, only first five terms are considered and the rest of the terms are neglected. When a uniform vertical distribution is attained, the following equation is used to calculate centerline air concentration:(3)
χ
(
x
,
y
=
0
,
z
)
=
Q
2
π
u
-
σ
y
L
.
For the estimation of radiation doses during early phase of nuclear emergency, the following pathways are taken into consideration:(i)
external and internal exposure from cloud shine,
(ii)
exposures from ground shine,
(iii)
internal exposure from resuspension inhalation and skin doses from deposited material onto the skin,
(iv)
acute and lifetime doses from early exposures.


## 2.2. Meteorological Conditions

For the execution of MACCS code, the hourly-based data for the year 2010 of Islamabad collected from the ground station (with an approximate height of 10 meters) has been used. As the near surface wind speed increases with altitude, the same data if used for the stake height (61 m) will overestimate the plume rise. This could produce significant underestimation of the radiation doses. To incorporate the effect of wind speed (m/s) with altitude, the following theoretical formula was used to estimate the wind speed at higher altitudes:(4)
u
=
u
o
(
z
z
o
)
p
,
where u = wind speed at height z (m/s), u
o = wind speed at surface (m/s), and p = parameter (dimensionless) that varies with stability class and surface roughness.The wind rose for the reported year is presented in Figure2. The percentage value of wind direction remained 3.7% in north, 8.7% in northeast, 8.4% in east, 13.8% in southeast, 10.8% in south, 38.3% in southwest, 6.4% in west, and 10% in northwest, respectively. Data presented in Figure 2 represents that dominant wind direction was southwest (38.3%) during year 2010. The dominant atmospheric stability class as presented in Figure 3 was stability class-F.Figure 2
Wind rose for sixteen directions.Figure 3
Dominant stability class.


## 3. Initial Conditions and Assumptions

### 3.1. Source Term and Release Fraction

The source term for the fission product in the reactor core was taken from the international published data [6] for a postulated accidental airborne release from Pakistan Research Reactor 1 (PARR-1), Islamabad, on upgraded power of 10 MW. The fractions of releases were based on USNRC document NUREG-1150. The release fractions of 1, 0.4, 0.3, 0.05, and 0.02 for noble gases, halogens, alkali metals, the tellurium group, and the Ba-Sr group, respectively, have been considered.


### 3.2. Meteorological Options

MACCS code provides different options for the selection of the meteorological data. The following four meteorological options (MO) have been used to estimate the radiation doses:(i)
MO-1: user specified day and hour,
(ii)
MO-2: weather bin sampling,
(iii)
MO-3: constant weather boundary data,
(iv)
MO-4: special case (maximum wind speed).


### 3.3. Release Scenarios

Scenario 1. In this scenario, the releases and meteorological data have been considered at height of 61 meters (stake releases). The total effective dose equivalent (TEDE) for whole body modeled using four different meteorological options is presented in Figure 4.Figure 4
Scenario 1 mean TEDE (7 days).Scenario 2.In this scenario, the releases and meteorological data have been considered at ground level (10 m). The TEDE for whole body modeled by using four different meteorological options is presented in Figure 5.Figure 5
Scenario 2 mean TEDE (7 days).


### 3.4. Comparison of Doses for Constant Weather Condition

The comparison of radiation doses of the scenarios has been made with the literature data (using the meteorological conditions of the site [6]) for 61-meter height using constant weather conditions. The comparison is presented in Figure 6.Figure 6
Comparison of mean TEDE (7 days) for constant weather.


## 3.1. Source Term and Release Fraction

The source term for the fission product in the reactor core was taken from the international published data [6] for a postulated accidental airborne release from Pakistan Research Reactor 1 (PARR-1), Islamabad, on upgraded power of 10 MW. The fractions of releases were based on USNRC document NUREG-1150. The release fractions of 1, 0.4, 0.3, 0.05, and 0.02 for noble gases, halogens, alkali metals, the tellurium group, and the Ba-Sr group, respectively, have been considered.


## 3.2. Meteorological Options

MACCS code provides different options for the selection of the meteorological data. The following four meteorological options (MO) have been used to estimate the radiation doses:(i)
MO-1: user specified day and hour,
(ii)
MO-2: weather bin sampling,
(iii)
MO-3: constant weather boundary data,
(iv)
MO-4: special case (maximum wind speed).


## 3.3. Release Scenarios

Scenario 1. In this scenario, the releases and meteorological data have been considered at height of 61 meters (stake releases). The total effective dose equivalent (TEDE) for whole body modeled using four different meteorological options is presented in Figure 4.Figure 4
Scenario 1 mean TEDE (7 days).Scenario 2.In this scenario, the releases and meteorological data have been considered at ground level (10 m). The TEDE for whole body modeled by using four different meteorological options is presented in Figure 5.Figure 5
Scenario 2 mean TEDE (7 days).


## 3.4. Comparison of Doses for Constant Weather Condition

The comparison of radiation doses of the scenarios has been made with the literature data (using the meteorological conditions of the site [6]) for 61-meter height using constant weather conditions. The comparison is presented in Figure 6.Figure 6
Comparison of mean TEDE (7 days) for constant weather.


## 4. Comparison of MACCS and InterRAS Codes

Comparison of MACCS code with InterRAS code [9, 10] for releases at the height of 61 meters using constant weather conditions has been made. The InterRAS code estimate radiation doses maximum up to 48 hours. The doses estimated using InterRAS code are comparable to doses estimated using MACCS code. The comparison of MACCS and InterRAS code results is presented in Figure 7.Figure 7
Comparison of MACCS and InterRAS codes.


### 4.1. Mean TEDE Trends (One Week)

Mean TEDE for different release heights and meteorological conditions has been analyzed and presented in Figures4–6. Through analysis, it was found that radiation doses reduce exponentially over the distance.It was found that the protective action “evacuation” was required at different distances in different accident situations. The radiation doses fall below the intervention level for evacuation (50 mSv) at 10 kilometers for Scenario 1, 10.75 kilometers for Scenario 2, and 2.75 kilometers for constant weather condition using site specific data.


### 4.2. Mean Distance for Intervention Levels

The intervention level for evacuation (50 mSv) was achieved at different distances for different emergency scenarios and for different meteorological options. Mean distances (Km) for taking intervention (evacuation) for different meteorological options (MO) are presented in Table1.Table 1
Evacuation intervention distance for different MET options.
Scenarios
Intervention (evacuation) distances (Km)
MO-1
MO-2
MO-3
MO-4
Scenario 1
2.00
7.75
10.00
0.25
Scenario 2
2.75
6.75
10.75
4.75
Literature [6]
—
—
2.75
—A comparison of MACCS and InterRAS code output results (mean effective dose equivalent to whole body for two days) for constant weather conditions has been presented in Figure7. From the trend analysis, it is found that the maximum of 5.10
×
10
3 mSv has been observed at mean distance of 0.25 km. At the nearby distances, the doses estimated using InterRAS code are comparable to doses estimated using MACCS code.


## 4.1. Mean TEDE Trends (One Week)

Mean TEDE for different release heights and meteorological conditions has been analyzed and presented in Figures4–6. Through analysis, it was found that radiation doses reduce exponentially over the distance.It was found that the protective action “evacuation” was required at different distances in different accident situations. The radiation doses fall below the intervention level for evacuation (50 mSv) at 10 kilometers for Scenario 1, 10.75 kilometers for Scenario 2, and 2.75 kilometers for constant weather condition using site specific data.


## 4.2. Mean Distance for Intervention Levels

The intervention level for evacuation (50 mSv) was achieved at different distances for different emergency scenarios and for different meteorological options. Mean distances (Km) for taking intervention (evacuation) for different meteorological options (MO) are presented in Table1.Table 1
Evacuation intervention distance for different MET options.
Scenarios
Intervention (evacuation) distances (Km)
MO-1
MO-2
MO-3
MO-4
Scenario 1
2.00
7.75
10.00
0.25
Scenario 2
2.75
6.75
10.75
4.75
Literature [6]
—
—
2.75
—A comparison of MACCS and InterRAS code output results (mean effective dose equivalent to whole body for two days) for constant weather conditions has been presented in Figure7. From the trend analysis, it is found that the maximum of 5.10
×
10
3 mSv has been observed at mean distance of 0.25 km. At the nearby distances, the doses estimated using InterRAS code are comparable to doses estimated using MACCS code.


## 5. Summary and Conclusion

The most frequent wind direction southwest, that is, 38.3%, and the prevailing stability class F, that is, 58.16%, were recorded during year 2010. Intervention distance for evacuation remained in the range 2.0 to 2.75 kilometers for user specified MET data; 6.75 to 7.75 kilometers for weather bin sampling method; 2.75 to 10.75 kilometers for constant weather data; and 0.25 to 4.75 kilometers for maximum wind speed. With site specific meteorological conditions, intervention distance range is very low, that is, 2.0 to 2.75 kilometers, which is very close to the already published value, that is, 1 to 2 kilometers, using “InterRAS” code.The seven-day TEDE has been modeled and compared with the reference levels [11] for evacuation, that is, 50 mSv for not more than seven days. It is found that the protective action “evacuation” was required at different distances in different accident scenarios.MACCS code during early phase of nuclear emergency reasonably estimates the radiation doses by taking into account the variation in meteorological data, release duration, and concentrations. MACCS code could be an efficient tool to be used during early phase of nuclear emergencies if integrated with pre- and postprocessor for handling and presentation of data.



---

*Source: 874134-2014-05-28.xml*