The uncertainty of the concentration calculation can be addressed in a similar way in which we earlier addressed trajectory uncertainty, by slightly offsetting the meteorological data to test the sensitivity of the advection calculation to the gradients in the meteorological data fields. Start by retrieving captex_control.txt and captex_setup.txt files.
- Recall from the trajectory discussion that the internal grid offset is in the horizontal and vertical. Therefore, to take full advantage of the ensemble, the release height should be changed to 200 m from the Setup Run / Starting Locations menu. Also to reduce the number of output frames, open the grid menu and change the sampling start time to 83 09 26 06 00, rename the output file to ensemble, then save to exit. To further speed up the calculation, we will only compute results through 18 UTC 27 September. Therefore, change the run duration from 68 to 49 hours. For future reference, save the CONTROL file as ensemble_control.txt.
- Open the Advanced / Configuration Setup / Concentration / Menu #4. Because the GUI script will cycle through 27 simulations, we can speed up the calculation by reducing the particle release number from 50000 to 5000. Also the particle position output files are not required. For future reference, save the namelist file as ensemble_setup.txt and then save to close all menus.
- To start the model simulation (this calculation may take some time: ~30 min), press the Special Runs / Ensemble / Meteorology menu tab and accept the prompts to continue and when finished, close the ensemble complete menu. The calculation results are saved to the \working directory and the base output file name, ensemble is appended with a simulation sequence number from .001 to .027 representing different meteorology grid offsets.
- To analyze the simulation results, a special pre-processing program needs to be run. Press the Display / Ensemble / View map menu tab. This opens up the ensemble display menu. The upper part of the menu is used to change how the input data files are processed. In the middle section, change the output selection option Concentration at Percentile from its default 95% to 50%. This means that the output maps will show the median concentration of all the ensemble members. There is one graphic for each time period (starting 26 06) but only the 4th time period is shown here. If you set the time aggregation period to some value greater than 1, then the ensemble results will represent multiple time periods and fewer graphics are produced.
- The View Map processing of the ensemble output files leaves a variety of different files (prob??, cmax??, cmean) in the working directory. These files were previously displayed as maps (using concplot) or they can be displayed as box plots. Start from the Display / Ensemble / box plot tab. This opens a simple menu where you can set the location of the data extraction point. Use the same location, sampler 510 (Little Valley, NY), that we examined in the concentration utility tutorial by entering 42.25 -78.80, then Execute Display.
- The resulting box plot shows the ensemble concentration predictions at Little Valley for each 3-hr averaging period during the simulation. The box bottom and top shows the 25th and 75th percentile concentrations, the single lines the 10th and 90th percentiles, and the circles the 5th and 95th percentile concentrations. The plus symbol and middle line show the mean and median concentrations. Note there are no unit conversion options for any of these programs so the results are all shown as g m-3. An additional graphic with the same axis labels is also created that shows the individual member concentrations indicated by member number.
- To relate the box plot with the original graphic of median concentrations, the median box plot concentration on the 26th between 15 and 18 is perhaps 1x10-9. Box plots are labeled according to the end time of the sample. Examining the graphic at 42.25 -78.80 shows the location to be very close to either the blue (10-10) or yellow (10-9) contours. Your results may differ from those shown here depending upon how many particles were released in your simulation.
- To view how the measured data fall within the ensemble estimates, we can overlay the data points on the boxplot. In the concentration utility tutorial we created a file with just the measurements from station 510 for display purposes. Unfortunately, the box plot program has no overlay capability, but we can manually (e.g. Paint) add the measurements to the box plot graphic after converting the numbers from pg back to g. The initial measured concentration corresponds with the highest predictions, smallest uncertainty range, and is close to both the ensemble mean and median concentrations. Model predictions drop off much more quickly than the measurements which is not accounted for by the meteorological grid ensemble. Note that 10-11 is approximately the detection threshold for the measured concentrations. Values lower than this are effectively zero. The ensemble also provides an estimate of the potential maximum concentration.
- For later reference, examine the box plot again at station 510 and note the concentration range (5%-95%) for the initial time period of the highest predicted concentrations (26 18). It is approximately 2x10-10 to 2x10-9 or about a factor of 10. The range is much larger for many of the other time periods because a significant number of members predicted zero concentrations.
There are many causes of concentration uncertainty, such as how well the meteorological data represent the true flow field, errors in the parameterization of physical processes in the model, and even the effects of atmospheric turbulence not accounted for by the model. In this case we looked at the issue of how well discrete data points represented the flow field. In the subsequent sections we will explore some of the other factors.