12.2 Turbulence Ensemble




Another component of the uncertainty in the concentration calculation is the contribution of random motions by atmospheric turbulence. The model already computes this turbulence when computing particle dispersal. However, normally we would try to release a sufficiently large number of particles to insure that each simulation gives similar results. In the turbulence ensemble approach, we reduce the particle release number and run multiple simulations, each with a different random number seed, and then examine the variations between simulations.

  1. Before running a simulation, we should delete all the ensemble files left over from the previous section. There is a Special Runs / Ensemble / cleanup tab that opens the menu. There are no options, just press Execute Script and all ensemble related files will be deleted from the working directory. Only the model output concentration files with the root name set in the Grid menu with the .000 suffix will be deleted.

  2. Start by retrieving the saved ensemble_control.txt and ensemble_setup.txt settings into the GUI menu. For this simulation, from the Setup Run menus, change the release height back to 10 m and change the output file from ensemble to a unique name such as ensturb.

  3. After the Setup Run and the Advanced Configuration menus have been retrieved, open the Advanced / Configuration Setup / Concentration / Menu #4 and change the particle release per cycle from 5000 to 1852. Although the turbulence ensemble can be composed of any number of members, to parallel the meteorological grid ensemble, 27 variations are also run in this case. The 1852 is determined by dividing 50000 (from the original simulation) by the number of members to maintain about the same total particle number when computing the mean concentration.

  4. After saving the change, press the Special Runs / Ensemble / Turbulence menu tab. There are no data entry options, just a prompt menu to ask if you really want to continue with the ensemble turbulence calculation. The calculation will start and as each member calculation is completed, a message is added to the simulation log and finally a completion message when member 27 finishes.

  5. Following the same procedure as previously with the meteorological ensemble, select the Display / Ensemble / View Map menu tab. This step is required to display probability maps or box plots because it invokes the pre-processor step that creates the probability files from the individual concentration simulations. Display any map and then open the box plot menu and enter the same position 42.25 -78.80, sampler 510 (Little Valley, NY), used in the previous section.

  6. The resulting box plot is similar to the previous example, but with a somewhat lower range between low and high concentrations; note that the ordinate scale is different. For instance, the time period with the initial measured concentration (26 18), shows the 90th percentile concentration uncertainty range (5th to 95th) is about a factor of four rather than a factor of 10 from the meteorological grid ensemble.

The approach shown here, using the model's random number generator, provides some guidance as to how much variability might be introduced by atmospheric turbulence. A similar approach would be to set the namelist parameter KRAND=4, which uses a different method to set a variable seed for the random numbers. The krand approach can be used with any concentration model executable. It is not limited to the variance version of the code.