Concentration / Special Simulations

Special simulations may require a different executable file, modifications to the Control file that are not supported by the GUI, or interactions with other items under the Advanced Menu tab. More information is provided below for each special simulation. Special model configurations may not be available for all operating systems. Note that most of the special simulations may be run using a single processor system or multiprocessor system (supporting MPI).

Run Matrix

Although the setup of the concentration matrix calculation is similar to that of the trajectory matrix calculation, there is an additional option that can be set in the Advanced Menu configuration tab that changes the nature of the concentration output file to produce a source-receptor matrix. This will be discussed further below. The matrix calculation is a way to set up the CONTROL file for multiple starting locations that may exceed the GUI limit of 6 under the Concentration Setup menu tab. Hundreds or thousands of starting points may be specified. The Run Matrix menu tab first runs a program that reads the CONTROL file with three starting locations and then rewrites the same CONTROL file with multiple locations. The multiple locations are computed from the number of starting points that fall in the domain between starting point 1 and starting point 2, where each new location is offset is the same as that between starting locations 1 and 3. For instance, if the original control file has three starting locations: #1 at 40N, 90W; #2 at 50N, 80W, and #3 at 41N, 89W; then the matrix processing results in a Control file with 121 starting locations, all one degree apart, all between 40N 90W and 50N 80W.

In the normal model execution mode, the concentration contributions from multiple sources are summed on the concentration grid, hence it is not possible to determine the fraction of the material comes from each source location. This can be seen in the illustration below using the above configuration for the first 12 hours of the sample case.

However, if the "Matrix" conversion module is checked in the Advanced Concentration Configuration menu tab, then the multiple source simulation maintains the identity of each source in the concentration output file. The Display Matrix menu tab permits extraction of information for individual sources or receptors from this special concentration output file. The results of the same simulation are shown in the illustration below. In this case the receptor check-box was set and the receptor location was identified as 45.0, -75.0 with normalization. Therefore the graphic illustrates the fractional contribution of each region to the first 12 hour average concentration at the designated receptor.

The following table illustrates the HYSPLIT matrix configuration. Emissions occur from each of N source locations and the receptors represent the concentration grid of M nodes. A single concentration output file is produced where each source contributes to its own concentration grid of M receptors. When selecting a "source" display, the M columns from the source location (row) represent the downwind concentration pattern for that source. When the receptor location (column) is selected, the contours represent the concentrations (each row of that column) contributing to that receptor from each of the source locations. Source and receptor grids should be of comparable resolution.

Run Ensemble

The ensemble form of the model, an independent executable, is similar to the trajectory version of the ensemble. The meteorological grid is offset in either X, Y, and Z for each member of the ensemble. The model automatically starts each member on a single processor in a multi-processor environment or cycles through the simulations on one processor. The calculation offset for each member of the ensemble is determined by the grid factor as defined in the Advanced Concentration Configuration Tab. The default offset is one meteorological grid point in the horizontal and 0.01 sigma units in the vertical. The result is 27 ensemble members for all offsets. The normal Setup Menu tab is used to configure the CONTROL file. Note that if fewer than 27 processors are available, the ensemble configuration menu permits starting the calculation at any ensemble member number within the valid range. Because the ensemble calculation offsets the starting point, it is suggested that for ground-level sources, the starting point height should be at least 0.01 sigma (about 250 m) above ground. The model simulation will result in 27 concentration output files named according the file name setting in the control file "{cdump}.{001 to 027}" with a suffix equivalent to the ensemble member number. On a single processor system, the calculation may take some. The menu will be locked until the simulation has completed. A message file window will open after termination. Computational progress may be monitored by noting the generation of new concentration output and message files with the ensemble number suffix in the /working directory. The concentration output from each member can be displayed through the concentration display menu tab. However, to display the probabilities associated with the multiple simulations, it is necessary to pre-process the data through the Display Ensemble menu tab. Using the default configurations for the sample simulation, the illustration below represents the 90th percentile concentrations aggregating all four output time periods.

For instance the blue contour in this 90th percentile plot represents the region in which only 10% of the ensemble members have air concentrations greater than 10-15.

Run Dust Storm

A model for the emission of PM10 dust has been constructed (Draxler, R.R, Gillette, D.A., Kirkpatrick, J.S., Heller, J., 2001, Estimating PM10 Air Concentrations from Dust Storms in Iraq, Kuwait, and Saudi Arabia, Atmospheric Environment, Vol. 35, 4315-4330) using the concept of a threshold friction velocity which is dependent on surface roughness. Surface roughness was correlated with geomorphology or soil properties and a dust emission rate is computed where the local wind velocity exceeds the threshold velocity for the soil characteristics of that emission cell. A pre-processing program was developed that accesses the HYSPLIT land-use file over any selected domain and modify the input CONTROL file such that each emission point entry corresponds with a "desert" (active sand sheet) land-use grid cell. The original PM10 flux equation was replaced by a more generic relationship (Westphal, D.L., Toon, O.B., Carlson, T.N., 1987. A two-dimensional numerical investigation of the dynamics and microphysics of Saharan dust storms. J. Geophys. Res., 92, 3027-3029).

The dust storm simulation is configured in the same way as the matrix calculation in that it is necessary to define three source locations, the first two representing the limits of the domain, and the third defining the emission grid resolution. The pre-processor then finds all emission points within that domain that have a desert category and modify the CONTROL file accordingly. The dust box must be checked in the advanced configuration menu to compute the PM10 emission rate. As an example, we can configure the model to run the large Mongolian dust storm of April 2001. An animation of the calculation results can be downloaded. To run the same simulation it will be necessary to obtain the first two weeks of northern-hemisphere meteorological analysis data (FNL.NH.APR01.001). A pre-configured CONTROL file (dust_conc) should be retrieved from the working directory. The CONTROL file defines the emission domain by the three starting locations: 35N-90W to 50N-120W with the grid increment to 36N-91W. There is no point in defining an emission grid of less than one-degree resolution because the resolution of the land-use data file is one-degree and the meteorological data is closer to two-degrees. Once the model is setup for the simulation, including the dust check-box in the configuration menu, execute the model from the Special Simulations / Run Dust Storm menu tab. The window shown below will open to indicate the revision of the CONTROL file.

The message indicates that the initial 3-location CONTROL file was reconfigured by the dustbdy program for 105 source locations. That means in the domain specified, 105 one-degree latititude-longitude grid cells were found to have a desert land-use category. If none are found, then the CONTROL file is deleted to prevent model execution. Click on Yes or No to continue - yes just deletes the window. The model execution will then start. PM10 pollutant dust particles are only emitted from those 105 cells where the wind speed exceeds the emission threshold. Therefore it is possible to have simulations with no emissions. An example of the output after 24 hours simulation time is shown in the illustration below.

The concentrations represent a 3-hour average from 21 to 24 hours after the start of the simulations. It is not possible to say exactly from when or where particles are emitted except to note that the 105 potential source locations are shown. The emission algorithm emits in units of grams, but the in configuring the concentration display, the units were changed to ug by defining a conversion factor of 106. Maximum concentrations are on the-order-of 100, but the layer was defined with a depth of 3 km to facilitate comparison with satellite observations. The simulation could be rerun with a smaller layer to obtain concentrations more representative of near-ground-level exposures.



Table of Contents