From the questions that have been received to AUSTAL2000, the most important were selected and listed on this page along with a reply. The questions are grouped into general questions (A) and specific issues (S).

Table of Contents


A13) How to classify the critics by R. Schenk in journal "Immissionsschutz" 01/2015?

The regulation TA Luft states how the public is to be protected from hazardous air pollution in the context of plant licensing according to the Federal Immission Control Act. Part of the regulation is a dispersion model based on guideline VDI 3945 Part 3, for which the reference implementation AUSTAL2000 has been developed on behalf of the Federal Environmental Agency. Entitled "AUSTAL 2000 ist nicht validiert" (AUSTAL 2000 not validated), R. Schenk provides in journal "Immissionsschutz" 01/2015 several examples to falsify AUSTAL2000 and guideline VDI 3945 Part 3.

A more detailed analysis reveals that the results by AUSTAL2000 are correct and that the contradictions highlighted by R. Schenk are based on fundamental errors in his argumentation. A detailed explanation is provided in "Immissionsschutz" 03/2015 (Alfred Trukenmüller, Wolfgang Bächlin, Wolfram Bahmann, André Förster, Uwe Hartmann, Heike Hebbinghaus, Ulf Janicke, Wolfgang J. Müller, Jost Nielinger, Ralf Petrich, Nicole Schmonsees, Uwe Strotkötter, Thomas Wohlfahrt, Sabine Wurzler: Reply to the critics of Schenk to AUSTAL2000 in Immissionsschutz 01/2015 - AUSTAL2000 is verified and validated; Erwiderung der Kritik von Schenk an AUSTAL2000 in Immissionsschutz 01/2015 - AUSTAL2000 ist verifiziert und validiert).


A12) What must be considered for the spatial average of odor hour frequencies?

The existance of an odor hour is detected by internal comparison of the hourly mean concentration in a given grid cell  with the threshold value. For this purpose the concentration must be spatially resolved  with a sufficient accuracy. For hourly means, in general a finer spatial resolution is required than for long time means. This aspect can be of particular importance at small distances to near-ground point sources. As a rule of thumb for such a case, at least 5 grid meshes should be between the source and the point of interest. If the grid meshes are too large, the hourly concentration plume is not well resolved and the resulting odor hour frequency may be under- or overestimated, depending on the emission rate.

Based on the results of a sufficiently resolved dispersion calculation, larger spatial averages can be derived. For example the auxiliary program A2KArea.jar can be used for this purpose. This distinctive feature of spatially averaging odor hour frequencies is due to the non-linear relationship between odor hour and concentration (yes/no decision: concentration above or below a threshold). For the spatial average of concentration in contrast, the values provided by the dispersion program on a coarse calculation grid can be directly used.


A11) What must be considered for the specification of coordinates?

  1. All coordinate specifications are in meters.
  2. All coordinate specifications except for the reference point must be of absolute value smaller 200000. The reference point must be chosen accordingly.
  3. In order to avoid internal rounding errors, specification of calculation grid (left and lower border, mesh width), of the reference point, and of anemometer position should be made in integer meters.

A10) How are dust components specified and evaluated?

For the substances

  • pm (dust, unspecified),
  • as (arsenic),
  • pb (lead),
  • cd (cadmium),
  • ni (nickel),
  • hg (mercury),
  • tl (thallium),
  • xx (unknown)

the diameter classes 1 to 4 and class "unknown" according to Annex 3 of TA Luft can be applied by appending "-1", "-2", "-3", "-4", or "-u" to the substance name (e.g. xx-1).  The components "1" and "2" define the parts of PM-10. The components "3", "4", and "u" are handled inside the dispersion calculation separately from components "1" and "2" because they have different settling velocities.

The concentration files (e.g. pb-y00a.dmna) contain the concentration sums of components "1" and "2", because the imission values of TA Luft refer to PM-10 (for substances hg and xx, the gaseous component is also included in the sum). The deposition files (e.g. pm-depa.dmna) contain the sum over all components.

If one is interested in the concentration of one of the components "3", "4", and "u", the dispersion calculation can be carried out with the gaseous component "xx" and explicit setting of the sedimentation and deposition velocity using NOSTANDARD options. Example for the output of the concentration for diameter class 3:

    os "NOSTANDARD;Vs=0.04;Vd=0.05"
    xx  100

Note that the sedimentation and deposition velocity specified in the NOSTANDARD options is used for all substances that are included in the dispersion calculation.


S19) What to consider when modifying the settings file?

In general, the settings file (austal2000.settings or austal2000n.settings) should not be modified. If so, carefully consider the modification and check the resulting effects in detail. Special aspects:

  • The evaluation of highest hourly means is implemented for gaseous substances only (according to the requirements of the current TA Luft). Hence the program aborts if it does not encounter a gaseous component during this evaluation. Even if it does, the evaluation is not useful because no automatic addition of the dust components 1 and 2 takes place as it is the case for annual and daily means. If for a substances highest hourly means is desired, the spectrum of allowed components must be set to "0-0".



S18) How should calms be defined in a AKTerm?

According to the specification of the German Weather Service DWD, calms (hours with no wind) are encoded in a AKTerm file by wind direction 0 (wind from north is defined by 360 deg). In addition, the wind speed should be set to 0: first for obvious plausibility, second to allow a consistent treatment of the wind direction redistribution according to paragraph 8.2 of TA Luft, Annex 3. Thus, a wind direction 0 should only occur together with wind speed 0 and vice versa. This is also the convention followed by the DWD. The program Austal2000 detects calms at hand of wind speed 0.


S17) What must be considered when applying externally generated wind fields?

  • The wind fields must be defined on the Arakawa-C net.

  • Before the actual dispersion calculation, any divergence in the provided wind field is internally removed by a simple adjustment ofthe z-components. To keep this modification at a minimum, the provided wind fields should be already sufficiently free of divergence.

  • Buildings: Buildings are resolved on the calculation grid. A cell is interpreted as buildingcell if the vertical component Vs of the wind vector at the bottom of the cell is less or equal -99 m/s. Below a building cell there may beno free cell (this would result either in an error message or in incorrect concentration values).

  • Terrain profile and nested grids: AUSTAL2000 adjusts the z-values specified in files srfa0li.dmna for the two boundary stripes of inner grids at hand of the values ofthe next coarser grid. This "tuning" guarantees a smooth transition between the grids. The resulting values are written to the wind fieldsgenerated by AUSTAL2000. If external fields are provided, they must containthese tuned z-values. AUSTAL2000 checks the z-values of the wind field files with the internally created ones andaborts with an according error message if they do not agree.



S16) Why is a quality level of 2 in some situations not sufficient to faithfully calculate the frequency of odor hours?

A sufficiently large number of simulation particles is required to avoid systematic effects in the calculation of odor hour frequencies (see the discussion in the annex of the program manual). Usually (one or several sources near ground and close together) a quality level of 2 is sufficient for time series calculations. However, if the sources are far apart or if a source, which is dominant with respect to the emission, has almost no contribution to the total concentration near ground due to its construction height and/or plume rise (but nevertheless takes most of particles, which are assigned to the sources according to their emissions), a considerably larger rate may be required. Once the systematic effect is avoided, the listed (absolute) statistical uncertainty is a faithful measure also for odor hours. The quality level in AUSTAL2000 (parameters qs, default value 0) can be set to a value up to 4 in standard operation. Higher values require the NOSTANDARD setting in the option string (parameter os).