...
The ensemble forecasts have 51 members, which will be assigned an extremity rank each. Using these 51 ranks the forecasts will be put in one of the 7 anomaly categories (as described in Table 1). This is done based on the arithmetic mean of the 51 ensemble member rank values (rank-mean). This rank-mean will also be a number between 1 and 100, but this time a real (not integer) number. If the anomaly is 50.5, that is exactly the normal (median) condition, i.e. no anomaly whatsoever. If the anomaly is below 50.5, then drier than normal conditions are forecast, if above 50.5, then wetter than normal. The lower/higher the anomaly value is below/above 50.5, the drier/wetter the conditions are predicted to be. The lowest/highest possible value is 1/100, if all ensemble members are 1 or /100 (the most extremely dry/wet). Then, based on this rank-mean, we define the dominant anomaly category (one of the 7 categories in Table 1) for the ensemble forecast, by placing the rank-mean into the right categories, as defined in Table 1 above. For example, all rank-mean values from 40.0 to 60.0, interpreted as 40.0<= <60.0, will be assigned to 'Near normal', or category-4.
The ensemble forecast anomaly was not based on the most probable of the 7 anomaly categories, as that would make it prone to jumpiness. For example, in the super uncertain case of 6, 8, 7, 7, 7, 9, 7 members being in each of the 7 anomaly categories, the forecast category (the dominant one) would be the 'High' category (cat-6), as that has the most members (9). However, it is likely that nearby river pixels could easily be only slightly different with 7, 9, 7, 7, 7, 7, 7 members in each category, in which case the dominant anomaly category would be the 'Low' category (cat-2), as now that has the most (again 9) members. It is worth mentioning that very uncertain cases are especially likely to happen at longer ranges. These two forecasts are only slightly different in terms of distribution, but the ensemble forecast anomaly categories would be almost the complete opposite of each other, making the signal look possibly very jumpy geographically. With the mean-rank definition we avoid this and simply assign the 'Near normal' category (cat-4) for both these forecasts, as the mean of the ranks are certainly very close to each other and both will be quite near the median.
Figure 2. Schematic of the forecast extremity ranking of the 51 ensemble members and the 7 anomaly categories in the context of the climatological distribution.
There is a consequence of the 0-value problem over dry or very dry areas (described above), as some or all of the low anomaly signals will not be impossible to occur. Please consider these idealised examples as a help to demonstrate the situation with the 0-value problem. In these, for simplicity reason, a fixed portion of the ensemble members are being 0 and the rest having a simple average rank. This way the rank-mean and dominant anomaly category computation is simpler and easier to interpret.
...
Example No-3: Lowest 70% of the climatology is 0:
Number of 0 members | Number of non-0 members | Average rank of 0 members | Average rank of non-0 members | Rank-mean | Dominant anomaly category |
---|---|---|---|---|---|
0 | 51 | NA (no member to rank) | 71 (the lowest possible rank for a non-zero |
member if 1-70 percentiles in the climatology are 0) | (0 * 35.5 + 51 * 71)/51 = 71 | Bit high (60-75) | |||
0 | 51 | NA | 100 | (0 * 35.5 + 51 * 100)/51 = 100 | Extreme high (90<) |
11 | 40 | 35.5 | 71 | (11 * 35.5 + 40 * 71)/51 = 63.34 | Bit high (60-75) |
11 | 40 | 35.5 | 100 | (11 * 35.5 + 40 * 100)/51 = 86.08 | High (75-90) |
21 | 30 | 35.5 |
71 | (21 * 35.5 + 30 * 71)/51 = 56.38 | Near normal (40-60) | |||
21 | 30 | 35.5 | 100 | (21 * 35.5 + 30 * 100)/51 = 73.44 | Bit high (60-75) |
36 | 15 | 35.5 |
71 | (36 * 35.5 + 15 * 71)/51 = 45.94 | Near normal (40-60) | |||
36 | 15 | 35.5 | 100 | (36 * 35.5 + 15 * 100)/51 = 54.47 | Near normal (40-60) |
51 | 0 | 35.5 | NA (no member to rank) | (51 * 35.5 + 0)/51 = 35.5 | Bit low (25-40) |
Example No-4: All percentiles of the climatology (1-99) are 0:
Number of 0 members | Number of non-0 members | Average rank of 0 members | Average rank of non-0 members | Rank-mean | Dominant anomaly category |
---|---|---|---|---|---|
0 | 51 | NA (no member to rank) | 100 (the lowest possible rank for a non-zero member if 1-99 percentiles in the climatology are 0) | (0 * 50.5 + 51 * 100)/51 = 100 | Extreme high (90<) |
11 | 40 | 50.5 | 100 | (11 * 50.5 + 40 * 100)/51 = 89.32 | High (75-90) |
21 | 30 | 35.5 | 100 | (21 * 50.5 + 30 * 100)/51 = 79.61 | High (75-90) |
36 | 15 | 35.5 | 100 | (36 * 50.5 + 15 * 100)/51 = 65.05 | Bit high (60-75) |
51 | 0 | 35.5 | NA (no member to rank) | (51 * 50.5 + 0)/51 = 50.5 | Near normal (40-60) |
- If all ensemble members are 0, then the ranks will spread evenly from 1 to 10 and the rank-mean will be around 5.5, so in the Extreme low category.
...