THE HYSPLIT4 SO2/SO4 CHEMISTRY MODULE
Last Revised: 29 August 2003


ABSTRACT

     In 1990, a simple Lagrangian-based sulfur chemistry module 
was developed (Rolph, 1992 and 1993) to be incorporated into the 
HYSPLIT 3.0 model. The sulfur module is nonlinear (in hydrogen 
peroxide) and includes gas- and aqueous-phase oxidation of sulfur 
dioxide (SO2) and dry and wet removal of SO2 and sulfate (SO-) 
particles. The current version is identical to the original, except 
that it has been incorporated into HYSPLIT 4.6. 

Rolph G.D., Draxler R.R. and de Pena R.G. (1993) The use of model-derived
and observed precipitation in long-term sulfur concentration and deposition
modeling.  Atmospheric Environment 27A, No 13, 2017-2037.

Rolph G.D., Draxler R.R. and de Pena R.G. (1992) Modeling sulfur
concentrations and depositions in the United States during
ANATEX. Atmospheric Environment 26, 73-93.


PC VERSION

	The model is run by executing the batch file "run_model",
which runs the executable (hysp_grs.so2) in the \Exec directory.
When the simulation completes, the output file "cavrg" will have been
created containing the hourly average concentration fields for all
species.  Executing "run_plot" will create an extract file (xdump)
with just the hourly sulfate concentrations and then a Postscript file
"concplot.ps" with maps of hourly concentrations. The output files 
are compatible with all standard Hysplit plotting programs. A 
Postscript file viewer or printer is required. The model is 
pre-configured to create output for a three-day test period, June 
24-27, 2003, over the Northeast states, using EDAS archive meteorology.
All required files are provided in the \Fixed directory. The 
configuration (CFG) and CONTROL files have been pre-set to run this 
example case. Other simulations can be run if the appropriate 
meteorological data have been downloaded and the control and 
configuration files modified accordingly.  A familiarity with the 
regular PC version of Hysplit is strongly suggested prior to trying 
run this version of the model.  More detail is provided below about 
the format of some of the input configuration files as well as the 
emissions inventory.  The use of different inventories may require 
reformatting. 


CONTROL file

	The CONTROL file cannot be edited through the standard HYSPLIT
GUI. Manual editing is required. Point source emissions are not 
supported. A gridded emission inventory is required. The emission 
inputs in the CONTROL file have slightly different meanings when using
a gridded inventory.  The emission duration is the same, but the 
emission rate is a multiplier to the gridded emission value.  Normally
it would be set to 1.0. Gridded emissions can produce many particles 
each time step. To minimize the number of particles, one can set the 
emission duration in the control file to a very small value (0.1 hrs)
and set the emission rate multiplyer in the control to 10, thus the 
total emissions are conserved.

	The top of the snapshot concentration grid should match the top
of the model domain to avoid the transport of particles that are not 
chemically transformed. All concentration grids should cover the same
spatial domain. Grid #1 should always be defined as the snapshot grid.
The reactions are determined by the species concentrations on the 
snapshot grid.

	The emission domain is determined by the source lat/lon limits 
set in the CONTROL file.  Note that the source domain must lie entirely within 
the meteorological grid limits, otherwise the results may be un-
predictable. The release height represents the height from the ground 
through which pollutants will be initially distributed. 

Three species must be defined for sulfur chemistry. Only one is defined
in the gridded emissions file.    

3
SO2
10.0 
0.1
00 00 00 00 00
SO4
10.0         
0.1 
00 00 00 00 00
SO4W
10.0         
0.1   
00 00 00 00 00 


The model requires a gridded emissions inventory for SO2, in units of
kg/hour.  Model output units will be ug/m3. The following abbreviations
are used to identify species:

#  ID	Mol WT	

1: SO2	64	- SO2 from emissions inventory
2: SO4	96	- SO4 dry air concentrations computed from model 
3: SO4W	96	- SO4 wet air concentrations computed from model

Wet and dry deposition is computed within the sulfate chemistry
module rather than in the standard Hysplit subroutines.  However,
the dry deposition velocities are computed within Hysplit (from 
parameters set in the control file) and passed to the sulfate module.


SETUP.CFG NAMELIST FILE

	The standard setup.cfg results in the release of one 
particle at each grid point per time step.  If more particles are
to be released (at substantial computational cost) then the NUMPAR
parameter should to increased to the appropriate multiple of the
initial number of non-zero emission grid points as noted in the 
MESSAGE file. As noted above in the discussion of the control 
file, the model is configured to release all particles in one
time step, each hour (hence qcycle=1.0).  About 350000 particles
will be created over the 3 day test simulation.  Longer duration
simulations will result in particles moving off the domain.

INITD=0 		run in 3D particle mode 
KHMAX=9999 		don't drop any old particles
NUMPAR=500 		one particle per point if < number grid pts
MAXPAR=350000           max number of particles in simulation
MAXDIM=3  		number of species on one particle
QCYCLE=1.0 		new emission cycle each hour
DELT=0.0 		automatic time step computation
ICHEM=4			forces conc grid = meteo grid


PRCHEM.DAT file

Contains constants that can be changed to modify the sulfate 
chemistry conversion rates.

0.90      CRITICAL RELATIVE HUMIDITY                    (RHCRIT)
0.015     DAY SO2-SO4/HR                                (KSO2DY)
0.002     NIGHT SO2-SO4/HR                              (KSO2NT)
0.65      DRY TO WET SO4 CONVERSION FRACTION            (D2WFR)
0.90      HYDROGEN PEROXIDE (PPB)                       (H2O2)
0.90      LIQUID WATER CONTENT (G/M3)                   (LWC)

Suggested values for other seasons:

	Winter	Spring	Summer	Autumn
RHCRIT	.90	.90	.90	.90
KSO2DY	.008	.020	.035	.015
KSO2NT	.001	.002	.003	.002
D2WFR	.65	.65	.65	.65
H2O2	0.7	1.5	2.0	1.2
LWC	0.2	1.0	2.2	0.8



EMISSIONS INVENTORY

	The file "emission.txt" contains all the information
that is required to interpred the data in the gridded emission
inventory file. The file's first record contains information 
about the internal grid cell size that is used by the dispersion
model to accumulate the file's emissions. The emission file 
defines the emissions at latititude-longitude points, the values 
at these points are accumulated in an internal grid, the size of 
which is defined on the first record.  This value can be 
arbitrarily changed according to the desired resolution of the 
simulation.  The pollutant puffs are released with an initial 
size comparable to the accumulation cell size.  Because the 
emission file data are remapped to an internal grid, the file 
can consist of emissions data on a regular grid or just a 
collection of individual cells.  The emission rate in the Control
file is used as an additional multiplication factor for the data 
in the emission file.  Also note that previously discussed
particle number restrictions still apply.  The initial number of 
particles are spread out over the duration of the emission and the 
number of grid cells that are defined in the emission domain. 
The format of the emission.txt file is given below:

Record #1
        I4      - Number (n) of pollutant species in file
        I4      - Number of emissions defined for each 24 hour period
        F10.4   - Conversion factor: file units to model units/hour
        2F10.4  - Accumulation cell size (latitude & longitude)

Record #2
        nA4     - Pollutant character identification each pollutant

Record #3
        A       - the /directory/filename of the emission data file


The 1985 NAPAP emissions inventory file was provided by Stu McKeen
(stu@al.noaa.gov) from NOAA's Aeronomy Laboratory. The file 
NOx_VOC_SO2 contains NOx, Volatile Organic Carbon, and SO2 
emissions. The file is an ASCII formatted file od hourly emissions 
for 24472 1/6 deg latitude by 1/4 degree longitude grids between
24.5 and 59.5 degrees north and 52-135 degrees west longitude.

The records were written with the following FORTRAN write statements:
	
	WRITE(10,'(2I4,2F10.4)')II,JJ,XLON,XLAT
	WRITE(10,'(1P,12E10.3)')(NOX(IH),IH=1,24)
	WRITE(10,'(1P,12E10.3)')(VOC(IH),IH=1,24)
        WRITE(10,'(1P,12E10.3)')(SO2(IH),IH=1,24)

II and JJ are indexes used in the programming. XLON (negative is 
west longitude), and XLAT are the longitude and latitude of the 
south-west corner of each grid square. The emissions in British 
tons per hour per grid cell. The indexing of NOX, VOC, and SO2
are:

	 1 = 0000-0100 (GMT) average
 	 2 = 0100-0200 (GMT) average
	....
	24 = 2300-2400 (GMT) average.

Only grids for which nonzero emissions of either NOX, VOC, or SO2
are within the NAPAP inventory are in the file.  All point, mobile
and nonmobile emissions are included, with no distinction in height 
of emissions from point sources.  The NAPAP85 summertime-weekday 
scenario is used. VOC emissions represent all nonmethane hydrocarbons.
That is, CO and CH4 are not included in the emissions.  VOC emissions
are for 80 degrees F.



COMPILATION (for reference only ... source code not available)

	The Hysplit SO2 program is a conditional compilation 
variation of the standard model.  The main concentration program, 
source code, hymodelc.F, is compiled with the parameter -DSO2 to 
create the executable hysp_so2.

The following subroutines are required in addition to the standard
Hysplit library routines:

prchem.f - computes SO2 to SO4 conversion, wet and dry deposition is
           automatically computed regardless of control file settings 

emsgrd.f - in the standard Hysplit library, uncomment the mass to
           volume conversion for ppm output. Conversion  is initiated
           by the specification of molecular weight in the control file

