Documentation of the Evaluation of CALPUFF and Other Long ...

Documentation of the Evaluation of CALPUFF and Other Long ... Documentation of the Evaluation of CALPUFF and Other Long ...

20.04.2013 Views

The puff splitting flag (item 2) is turned on using the IRESPLIT input option. IRESPLIT consists of 24 values corresponding to the hour of the day with values that are either 0 or 1, where 1 turns on the puff splitting flag for puffs. Once the puff splitting flag is turned on, it remains on until the puff splitting occurs at which point the puff splitting flagged is turned off until it is turned back on again by IRESPLIT. The default setting for IRESPLIT is to have all hours zero except for setting hour 17 to 1. The reasoning behind this is to invoke puff splitting in the evening when a nocturnal inversion occurs and there is a decoupling of the winds between the nocturnal inversion layer and above the nocturnal inversion (i.e., the residual “mixed layer”). Setting IRESPLIT to all zeros will result in the puffs never performing vertical puff splitting and setting IRESPLIT to all ones will result in the puff splitting flag always turned on and puffs will always split when the other three criteria for vertical puff splitting are met. The default value for the previous hours minimum mixing height value (item 3) is ZISPLIT = 100 m. This minimum value is used to assure that the current mixing height is not negligible. The ratio of the previous hours mixing height to maximum mixing height encountered by the puff (item 4) is controlled by the ROLDMAX parameter with a default value of 0.25. When vertical puff splitting occurs in CALPUFF, the number of puffs that the puff is split into is controlled by the NSPLIT parameter that has a default value of 3. Horizontal puff splitting occurs when the puff concentrations are above a minimum value (CNSPLITH), the puff has a minimum width that is defined by its sigma‐y in grid cell units (SYSPLITH) and the minim puff elongation rate (SYSPLITH per hour) is above a SHSPLITH factor. The default minimum concentration is CNSPLITH = 10 ‐7 g/m 3 (0.1 µg/m 3 ). Default SYSPLITH value is 1.0 and default SHSPLITH factor is 2.0. When horizontal puff splitting occurs in CALPUFF the number of puffs the puff is split into is controlled by the NSPLITH parameter that has a default of 5. Eight CALPUFF puff splitting sensitivity tests were conducted, which are defined in Table 6‐3. When vertical and horizontal puff splitting occurs in CALPUFF, the default number of puffs to split into was used in the CALPUFF sensitivity tests (i.e., NSPLIT = 3 and NSPLITH = 5). The NOSPLIT sensitivity test set MSPLIT = 0 so no vertical or horizontal puff splitting was allowed to occur. The DEFAULT puff splitting turned on puff splitting (MSPLIT = 1) but only turned on the vertical puff splitting flag at hour 17 every day. Whereas, the ALLHRS sensitivity test made sure that the vertical puff splitting flag was turned on all the time (i.e., IRESPLIT = 24*1) removing criteria 2 from the vertical puff splitting requirement. The ZISPLIT sensitivity test set ZISPLIT to zero thereby removing criteria 3 in the vertical puff splitting, as well as requirement 2 (like ALLHRS). ROLD relaxed the minimum ratio of the previous hours to maximum mixing height for vertical puff splitting from 0.25 to 0.50. The SYS sensitivity test allows horizontal puff splitting to occur more frequently by allowing puff splitting to occur with a puff sigma‐y value is greater than SYSPLITH values of 0.1 (2.6 km) versus the default 1.0 (36 km) value. The last sensitivity test combines the ROLD and SYS sensitivity tests. 131

Table 6‐3. Summary of CALPUFF puff splitting sensitivity tests performed using the ETEX database. Sensitivity Test MSPLIT NSPLIT IRESPLT ZISPLIT ROLDMAX NSPLITH SYSPLITH CNSPLITH NOSPLIT 0 NA NA NA NA NA NA NA DEFAULT 1 3 Hr 17=1 100 0.25 5 1.0 10 ‐7 ALLHRS 1 3 24*1 100 0.25 5 1.0 10 ‐7 CNSMIN 1 3 24*1 100 0.25 5 1.0 10 ‐20 ZISPLIT 1 3 24*1 0 0.25 5 1.0 10 ‐20 ROLD 1 3 24*1 0 0.50 5 1.0 10 ‐20 SYS 1 3 24*1 0 0.25 5 0.1 10 ‐20 SYSROLD 1 3 24*1 0 0.50 5 0.1 10 ‐20 Figure 6‐20 displays the spatial model performance statistics for the CALPUFF puff splitting sensitivity tests. The DEFAULT, ALLHRS and CNSMIN CALPUFF sensitivity tests obtained the exactly same model performance statistics indicating that CALPUFF model performance was not affected by the IRESPLT and CNSMIN puff splitting parameters. There are some small difference in the spatial model performance statistics for the other CALPUFF puff splitting sensitivity tests with the ROLD parameter having the biggest effect when changed from 0.25 to 0.50 that improved model performance a couple of percentage points for the FMS, POD and TS spatial statistics but degraded the FAR spatial statistic by several percentage points. 16% 14% 12% 10% 8% 6% 4% 2% 0% 14% 12% 10% 8% 6% 4% 2% 0% NOSPLIT NOSPLIT DEFAULT Figure of Metric in Space (FMS) (Perfect = 100%) ALLHRS CNSMIN ZISPLIT ROLD Probability of Detection (POD) (Perfect = 100%) DEFAULT ALLHRS CNSMIN ZISPLIT ROLD SYS SYS SYSROLD SYSROLD Figure 6‐20. Spatial model performance statistics for the CALPUFF puff splitting sensitivity tests. 132 64% 63% 62% 61% 60% 59% 58% 57% 56% 12% 10% 8% 6% 4% 2% 0% NOSPLIT NOSPLIT DEFAULT DEFAULT False Alarm Rate (FAR) (Perfect = 0%) ALLHRS CNSMIN ZISPLIT Threat Score (TS) (Perfect = 100%) ALLHRS CNSMIN ZISPLIT ROLD ROLD SYS SYS SYSROLD SYSROLD

Table 6‐3. Summary <strong>of</strong> <strong>CALPUFF</strong> puff splitting sensitivity tests performed using <strong>the</strong> ETEX<br />

database.<br />

Sensitivity<br />

Test MSPLIT NSPLIT IRESPLT ZISPLIT ROLDMAX NSPLITH SYSPLITH CNSPLITH<br />

NOSPLIT 0 NA NA NA NA NA NA NA<br />

DEFAULT 1 3 Hr 17=1 100 0.25 5 1.0 10 ‐7<br />

ALLHRS 1 3 24*1 100 0.25 5 1.0 10 ‐7<br />

CNSMIN 1 3 24*1 100 0.25 5 1.0 10 ‐20<br />

ZISPLIT 1 3 24*1 0 0.25 5 1.0 10 ‐20<br />

ROLD 1 3 24*1 0 0.50 5 1.0 10 ‐20<br />

SYS 1 3 24*1 0 0.25 5 0.1 10 ‐20<br />

SYSROLD 1 3 24*1 0 0.50 5 0.1 10 ‐20<br />

Figure 6‐20 displays <strong>the</strong> spatial model performance statistics for <strong>the</strong> <strong>CALPUFF</strong> puff splitting<br />

sensitivity tests. The DEFAULT, ALLHRS <strong>and</strong> CNSMIN <strong>CALPUFF</strong> sensitivity tests obtained <strong>the</strong><br />

exactly same model performance statistics indicating that <strong>CALPUFF</strong> model performance was not<br />

affected by <strong>the</strong> IRESPLT <strong>and</strong> CNSMIN puff splitting parameters. There are some small<br />

difference in <strong>the</strong> spatial model performance statistics for <strong>the</strong> o<strong>the</strong>r <strong>CALPUFF</strong> puff splitting<br />

sensitivity tests with <strong>the</strong> ROLD parameter having <strong>the</strong> biggest effect when changed from 0.25 to<br />

0.50 that improved model performance a couple <strong>of</strong> percentage points for <strong>the</strong> FMS, POD <strong>and</strong> TS<br />

spatial statistics but degraded <strong>the</strong> FAR spatial statistic by several percentage points.<br />

16%<br />

14%<br />

12%<br />

10%<br />

8%<br />

6%<br />

4%<br />

2%<br />

0%<br />

14%<br />

12%<br />

10%<br />

8%<br />

6%<br />

4%<br />

2%<br />

0%<br />

NOSPLIT<br />

NOSPLIT<br />

DEFAULT<br />

Figure <strong>of</strong> Metric in Space (FMS)<br />

(Perfect = 100%)<br />

ALLHRS<br />

CNSMIN<br />

ZISPLIT<br />

ROLD<br />

Probability <strong>of</strong> Detection (POD)<br />

(Perfect = 100%)<br />

DEFAULT<br />

ALLHRS<br />

CNSMIN<br />

ZISPLIT<br />

ROLD<br />

SYS<br />

SYS<br />

SYSROLD<br />

SYSROLD<br />

Figure 6‐20. Spatial model performance statistics for <strong>the</strong> <strong>CALPUFF</strong> puff splitting sensitivity<br />

tests.<br />

132<br />

64%<br />

63%<br />

62%<br />

61%<br />

60%<br />

59%<br />

58%<br />

57%<br />

56%<br />

12%<br />

10%<br />

8%<br />

6%<br />

4%<br />

2%<br />

0%<br />

NOSPLIT<br />

NOSPLIT<br />

DEFAULT<br />

DEFAULT<br />

False Alarm Rate (FAR)<br />

(Perfect = 0%)<br />

ALLHRS<br />

CNSMIN<br />

ZISPLIT<br />

Threat Score (TS)<br />

(Perfect = 100%)<br />

ALLHRS<br />

CNSMIN<br />

ZISPLIT<br />

ROLD<br />

ROLD<br />

SYS<br />

SYS<br />

SYSROLD<br />

SYSROLD

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!