15.1 Radioactive Decay and Dose




In this section we're going to assume the release consisted of radioactive particles which decay and deposit. The air concentration and deposition amounts will be converted to dose. Although we will try to make some reasonable assumptions, this is still a simplistic example of what can be a complex modeling problem. The details presented here should just be used as an example in configuring the model and not as a template for simulating a nuclear accident. To avoid unintended options, start by pressing the Reset button on the main menu to clear all previous settings.

  1. The decay of radioactive pollutants is treated the same way as transformation and deposition. The mass on a particle is reduced using the time constant approach:
    • m2 = m1 exp(-βrad Δt)

    In addition to the particle mass, deposited radioactive pollutants also decay, and hence deposition amounts are also reduced each time step. The decay constant for radioactive processes is defined by the half-life (T½):
    • βrad = ln 2 / T½

    In HYSPLIT the decay process starts at the time of the particle release, however in reality the decay starts immediately after the fission product is created. From the perspective of a nuclear accident, the time of release to the atmosphere may not coincide with the time the fission has stopped (reactor shutdown), or the valid time of the emission inventory. In any event, the longer the half-life or the shorter the duration of the release, the less this becomes a problem. This issue will be discussed in more detail in another section.

  2. For the purposes of this example, we will assume a release from a hypothetical nuclear reactor located at the CAPTEX release point, essentially the same configuration we have been using for all the previous examples. Start by retrieving the previously saved captex_control.txt and captex_setup.txt settings into the GUI menu. We will assume that a 3000 MW reactor that has been running for one year created about 1.0E+17 Bq of Cesium-137 (and hundreds of other fission products). We also assume that 10% of the Cesium is released in the accident over a one-hour period.

  3. To configure this release, set the run duration to 25 h, open the Setup Run / Pollutant menu and enter the C137 label, release amount 1.0E+16, and the 1.0 hour duration. Save the changes and open the Setup Run / Grid menu and change the name of the output file to c137.bin, and add the second level 0 for deposition. Then open the Setup Run / Deposition menu and check the Cs-137 radio-button to set the deposition defaults for Cesium. Save the changes and then run the model.

  4. To convert the air concentration and ground deposition values into a number that can more easily be related to health effects, we need to find the dose conversion factors for Cesium-137. For the airborne plume assume 3.4E-11 (REM/hr) / (Bq/m3) and 1.1E-12 (REM/hr) / (Bq/m2) for the deposited material.

  5. Because the air concentrations are three-hour averages and we want to compute the dose in REM, the air concentration dose conversion factor will be 1.0E-10 and if we compute the total deposition over 24 hours, then the deposition dose conversion factor becomes 2.6E-11. Open the Display / Contours menu and enter the dose conversion factors. Also check the total radiobutton to see the accumulated deposition. Before executing the display, the plot labels should be modified to reflect the changes. Just changing the units field in the menu to REM is insufficient because it will incorrectly show as REM/m3.

  6. To make more complex changes to the graphics label, open the Advanced / File Edit / Border Labels menu and enter Cesium 137 for a new plot title, Dose Equivalent for the map type, and REM for the mass units into their respective text entry fields. Leave the other fields blank. Then save the changes and go back and execute the Contour Display menu and advance the frames to the end to see the one-day accumulated ground-shine dose map.

  7. Because Cesium has such a long-half life, the deposited material will continue to be a problem for quite some time. The decay after one year is so small that we can compute the one-year ground-shine by multiplying the dose conversion factor by 8760 rather than 24 (9.6E-09). The accumulated dose becomes 365 times larger. The conversion factor could also be adjusted for the actual Cesium decay.

  8. Although the prospective maximum dose from this simple example appears to be small (1 mR) for the year, other situations for a similar release amount may have different consequences. For instance, the 25 km resolution concentration grid would under-estimate air concentration and deposition near the source. A finer grid is required to assess the near source impacts. In addition, this event had no wet deposition during the early stages of transport, further reducing the impact of the release. To facilitate some additional testing save the configuration to conc_decay_control.txt and conc_decay_setup.txt

The example presented here was intended show how the model can be configured for a radiological simulation. The emission and dose conversion values are only intended as examples and a guideline for an actual accident scenario, which is composed of a complex mix of short- and long-lived radionuclides. Although there are many public sites with radiological monitoring data, frequently these sites only show background values, except when there has been some event.