Hysplit: Radionuclide Applications
Hysplit: Radionuclide Applications
Hysplit: Radionuclide Applications
Radionuclide Applications
Roland Draxler
NOAA Air Resources Laboratory
HYSPLIT Overview
HYSPLIT is not just a transport and dispersion model but a
complete system for computing trajectories, dispersion, and
deposition
Can be applied to different air quality problems
Specific applications depend upon configuration of input
files and the use of pre- and post-processing executables
Executable library contains over 100 applications
This presentation is a highly technical focus on configuring
the model for radiological applications
This presentation was developed as a reference document
• Radionuclide mass is
assigned to each particle
• Mass is decayed after it is
released
• A user defined three-
dimensional grid covers
the domain
• Particle masses are
summed in each grid cell
• The mass sum is divided
by the grid cell volume to
obtain air concentration
control
namelist
meteorology HYSPLIT CONCPLOT
binary air
concentration
and deposition
CON2STN
TIMEPLOT
C2DATEM
All the following examples
have been configured for
demonstration purposes
and computational speed!
TIMEPLOT
-iTokai.txt model predictions input text file
-sJAEA_C137.txt measured data file
Model
Measured
1 1
CPAR CPAR
5.887E+13 0.0
51.0 0.0
TG_{MMDDHH}
DG_031500 xxxx
mergelist
fdnpp.bin xxxx xxxx xxxx xxxx xxxx xxxx
4-character pollutant
ID label in output file
control
namelist HYSPLIT binary C2ARRAY
meteorology TCSOLVE
TG_{MMDDHH}
1.30E-12 9.79E-14 0.00E+00 0.00E+00 0.00E+00 0.00E+00 0.00E+00 0.00E+00 0.00E+00 0.00E+00 0.00E+00 0.00E+00 0.00E+00 0.00E+00 0.00E+00 0.00E+00 0.00E+00 1.03E+04
1.67E-13 5.85E-13 5.41E-13 0.00E+00 0.00E+00 0.00E+00 0.00E+00 0.00E+00 0.00E+00 0.00E+00 0.00E+00 0.00E+00 0.00E+00 0.00E+00 0.00E+00 0.00E+00 0.00E+00 1.10E+04
2.06E-13 6.65E-13 1.50E-12 1.47E-12 0.00E+00 0.00E+00 0.00E+00 0.00E+00 0.00E+00 0.00E+00 0.00E+00 0.00E+00 0.00E+00 0.00E+00 0.00E+00 0.00E+00 0.00E+00 1.87E+05
0.00E+00 5.27E-15 2.65E-15 1.06E-13 6.08E-14 0.00E+00 0.00E+00 0.00E+00 0.00E+00 0.00E+00 0.00E+00 0.00E+00 0.00E+00 0.00E+00 0.00E+00 0.00E+00 0.00E+00 3.83E+03
4.72E-15 1.03E-14 0.00E+00 1.33E-14 0.00E+00 0.00E+00 0.00E+00 0.00E+00 0.00E+00 0.00E+00 0.00E+00 0.00E+00 0.00E+00 0.00E+00 0.00E+00 0.00E+00 0.00E+00 4.05E+02
0.00E+00 6.54E-15 0.00E+00 2.62E-15 1.39E-15 0.00E+00 0.00E+00 0.00E+00 0.00E+00 4.35E-17 0.00E+00 0.00E+00 0.00E+00 0.00E+00 0.00E+00 0.00E+00 0.00E+00 3.55E+02
0.00E+00 0.00E+00 0.00E+00 0.00E+00 5.80E-14 1.39E-13 3.29E-13 1.64E-13 4.00E-15 5.43E-15 2.25E-14 0.00E+00 0.00E+00 0.00E+00 0.00E+00 0.00E+00 0.00E+00 5.81E+03
0.00E+00 0.00E+00 0.00E+00 0.00E+00 0.00E+00 1.51E-14 5.88E-14 3.52E-14 2.95E-15 0.00E+00 0.00E+00 0.00E+00 0.00E+00 0.00E+00 0.00E+00 0.00E+00 0.00E+00 3.36E+02
TCSOLVE
-ic2array.csv input CM array (D = Bq/m3)
-u0.001 convert R from mBq to Bq
-otcsolve.txt source term solution in Bq
Alternate Approach:
• Edit the CM file into a square matrix
• Use EXCEL functions MINVERSE() and MMULT()
activity.txt file
RNUC Hr= 0.00 sec Bq rem/h Bq/m3 rem/h Bq/m2 rem/Bq
1.0 95 Nb 3.02000E+06 1.62000E+13 1.26000E-08 2.62000E-10 2.38000E-06
51.0 110 Ag 2.18000E+07 6.48000E+12 4.57000E-08 9.29000E-10 2.38000E-06
00 00 00 00 00 132 Te 2.82000E+05 1.62000E+16 3.36000E-09 7.63000E-11 7.38000E-07
NGAS 131 I 6.94656E+05 8.10000E+15 6.08400E-09 1.31040E-10 7.38000E-07
1.0 133 I 7.49000E+04 8.10000E+15 9.94000E-09 2.22000E-10 7.38000E-07
51.0 133 Xe 4.52995E+05 3.66000E+17 5.00400E-10 0.00000E+00 0.00000E+00
00 00 00 00 00 134 Cs 6.50000E+07 8.10000E+14 2.54000E-08 5.33000E-10 4.63000E-07
1 137 Cs 9.52093E+08 8.10000E+14 3.34080E-11 1.07640E-12 4.63000E-07
CONTROL file
Cloud Ground
mBq/m3
mBq/m3
R/h
R/h
mBq/m3 Total
CONTROL
39.7697 -86.1498 2000.0
39.7697 -86.1498 3000.0
39.7697 -86.1498 4000.0 CONTROL and EMISSIONS
39.7697 -86.1498 5000.0
39.7697 -86.1498 6000.0
06
for a nuclear detonation
0
10000.0
YYYY MM DD HH DURATION(hhhh) #RECORDS
1
EMITIMES
YYYY MM DD HH MM DURATION(hhmm) LAT LON HGT(m) RATE(/h) AREA(m2) HEAT(w)
C:/hysplit4/working/
2010 05 19 09 9999 18
hysplit.t06z.namf
3 2010 05 19 09 00 0006 39.7697 -86.1498 1000.0 0.2 0.0 0.0
P005 2010 05 19 09 00 0006 39.7697 -86.1498 1000.0 0.1 0.0 0.0
0.0 2010 05 19 09 00 0006 39.7697 -86.1498 1000.0 0.0 0.0 0.0
0.1 2010 05 19 09 00 0006 39.7697 -86.1498 2000.0 0.5 0.0 0.0
00 00 00 00 00
2010 05 19 09 00 0006 39.7697 -86.1498 2000.0 0.1 0.0 0.0
P010
0.0 2010 05 19 09 00 0006 39.7697 -86.1498 2000.0 0.0 0.0 0.0
0.1 2010 05 19 09 00 0006 39.7697 -86.1498 3000.0 1.2 0.0 0.0
00 00 00 00 00 2010 05 19 09 00 0006 39.7697 -86.1498 3000.0 0.2 0.0 0.0
P015 3 2010 05 19 09 00 0006 39.7697 -86.1498 3000.0 0.1 0.0 0.0
0.0 1.0 1.0 1.0
0.44 0.0 0.0 0.0 0.0
2010 05 19 09 00 0006 39.7697 -86.1498 4000.0 2.5 0.0 0.0
0.1
00 00 00 00 00 0.0 0.0 0.0 2010 05 19 09 00 0006 39.7697 -86.1498 4000.0 0.2 0.0 0.0
0.0
1 0.0
2010 05 19 09 00 0006 39.7697 -86.1498 4000.0 0.1 0.0 0.0
0.0 0.0 1.0 1.0 1.0 2010 05 19 09 00 0006 39.7697 -86.1498 5000.0 2.5 0.0 0.0
0.01 0.01 0.23 0.0 0.0 0.0 0.0
0.0 0.0 0.0
2010 05 19 09 00 0006 39.7697 -86.1498 5000.0 0.2 0.0 0.0
3.0 3.0 0.0 2010 05 19 09 00 0006 39.7697 -86.1498 5000.0 0.1 0.0 0.0
./ 0.0
cdump 1.0 1.0 1.0
2010 05 19 09 00 0006 39.7697 -86.1498 6000.0 1.7 0.0 0.0
2 0.15 0.0 0.0 0.0 0.0 2010 05 19 09 00 0006 39.7697 -86.1498 6000.0 0.2 0.0 0.0
0.0 0.0 0.0
0 100 0.0
2010 05 19 09 00 0006 39.7697 -86.1498 6000.0 0.1 0.0 0.0
00 00 00 00 00 0.0
00 00 00 00 00
00 01 00
8/30/2019 Air Resources Laboratory 43
Fission and Dose References for CON2REM
• The radionculide yield per 100 fissions was taken from: T.R. England
and B.F. Rider, Los Alamos National Laboratory, LA-UR-94-3106;
ENDF-349 (1993) digitally tabulated and available from
http://ie.lbl.gov/fission.html
• Not part of the HYSPLIT distribution, but an activity.txt file can be
created for any time after the time of fission using England and Rider
• Radionuclide inventories were computed using the rate of 1.45E+23
fissions per kT for nuclear detonations.
• For thermal reactions we assume 3000 MW-hours equals 2.58 kT
• The external dose rate for cloud- and ground-shine is computed from
the factors given by Eckerman K.F. and Leggett R.W. (1996)
DCFPAK: Dose coefficient data file package for Sandia National
Laboratory, Oak Ridge National Laboratory Report ORNL/TM-13347
for %%H in (1409 1412 1415 1418 activity.txt: note that all emissions = 1.0
1421 1500 1503 1506 1509 1512 1515 Mass Nucl T1/2 Emissions Cloudshine Groundshine Inhalation
Hr= 0.00 sec Bq rem/h Bq/m3 rem/h Bq/m2 rem/Bq
1518 1521 1600 1603 1606 1609) do ( 0001 NB95 3.02000E+06 1.00000E+00 1.26000E-08 2.62000E-10 2.38000E-06
0002 A110 2.18000E+07 1.00000E+00 4.57000E-08 9.29000E-10 2.38000E-06
0003 T132 2.82000E+05 1.00000E+00 3.36000E-09 7.63000E-11 7.38000E-07
CON2REM 0004 I131 6.94656E+05 1.00000E+00 6.08400E-09 1.31040E-10 7.38000E-07
0005 I133 7.49000E+04 1.00000E+00 9.94000E-09 2.22000E-10 7.38000E-07
-iCG_03%%H 0006 X133 4.52995E+05 1.00000E+00 5.00400E-10 0.00000E+00 0.00000E+00
0007 C134 6.50000E+07 1.00000E+00 2.54000E-08 5.33000E-10 4.63000E-07
-oDG_03%%H 0008 C137 9.52093E+08 1.00000E+00 3.34080E-11 1.07640E-12 4.63000E-07
-s1 match to input 0009 B140 1.10160E+06 1.00000E+00 2.90520E-09 6.84000E-11 1.03000E-07
0010 L140 1.44979E+05 1.00000E+00 3.99600E-08 7.77600E-10 1.07000E-07
-t0 no decay
CONCSUM
-iDG_03%%H • The 4-character radionuclide label in the
-oSG_03%%H -l activity.txt file must match the input
-pDOSE label output
)
species ID in the CG_* input file created
by the command line options of
CONMERGE -imergelist.txt -ofdnpp.bin condecay.
mBq/m3
R/h