...
In this section, examples are given with simplified rank distributions in order to demonstrate the rank-mean and rank-std computation for the expected anomaly and uncertainty category calculations. Based on these examples, the users can have a feel on how the rank-mean and rank-std values change with the changing underlying distributions. In addition, the impact of the 0-value in the climatology and in the ensemble forecasts is also demonstrated for different severity of the 0-value problem, with the complexity of these dry cases.
In these simplified examples, a fixed portion of the climatological and/or ensemble forecast distribution is 0, while the non-zero ensemble members have just one rank for the very dry cases and 5 for the non-zero cases, with members having the same rank in each group. This way, the computation methodology can be demonstrated in a simple way that is easier to interpret.
The non-zero value examples highlight that shifting the same rank distribution 'up' or 'down' (i.e. for wetter or drier) does not change the standard deviation (and thus the uncertainty). Also, after 'narrowing' the rank distribution over the same mean value (i.e. making it cluster more), the mean does not change, but the uncertainty drops markedly. Moreover, in a similar manner, by adding extreme members (i.e. near 1 or 100 very extreme members), even if only very few members, the uncertainty will increased quite substantially.
For those forecasts, when some portion or all of the climatological percentiles are 0, it is a general rule of thumb that as the percentage of zero climate percentiles increase, it gets more and more difficult to end up with negative expected forecast anomalies. The lowest possible rank-mean values are going to happen for forecasts with all 0 values, in which case the forecast rank-mean is going to be determined by the size of the 0-value section of the climatology. For example, if the lowest 30% of the climatology is 0, maybe in southern Spain somewhere on a small river catchment, then the rank-mean of the forecast of all a constant 0-value will be about 15, which will put this forecast in the expected anomaly category of 'Low'. But, with 70% of climatology being 0, say further into the Sahara, there the driest possible ensemble forecast of only 0 values is only going to be rank-mean of 35, so with expected anomaly category of 'Bit low'. Drier than this anomaly is simply physically not possible for such a climatological climatologically dry place.
So, for the option of 10% of climatology being 0, the absolute minimum possible forecast rank-mean is 5.5, while for 30% it will be 15.5 and for the totally dry climatology, where all 99 percentiles are zero, the rank-mean will be 50.5. The extent of which the rank-mean of the forecast increases above these minimum rank-mean values, depends on how many of the ensemble members will be non-zero and with which actual rank (determined by the non-zero section of the climatology). For example, one of the most extreme cases is when all 99 climatological percentiles are 0 and none of the ensemble forecast members are 0. For this super unlikely to occur event, the rank-mean of the forecast will always be 100 (and the expected forecast anomaly category 'Extreme high'), regardless of the actual ensemble member values (i.e. how much higher they are than 0). So, even if all forecast ensemble member river discharge values are very low, say from 0.12 to 0.23, the forecast rank-mean will still be 100 and the expected forecast anomaly category 'Extreme high'.
Panel | ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
|
Panel | |||||||
---|---|---|---|---|---|---|---|
| |||||||
In these examples, again for simplicity reasons, the climatological and forecast values will only be in one of 2 categories, either 0-value or non 0-value. This way, the main impact of the 0/non-0 value issue can be demonstrated. In the tables below, the numbers and the related average ranks are given for the two groups of 0 and non-0 ensemble members, with the rank-mean, rank-std and expected anomaly and uncertainty categories determined from those cases. There are 4 tables, with 10% / , 30% / , 70% and 100% of 0 values -value in the climatology (i.e increasingly dry climate). For example, in the 7th row of the 1st table with 10% of 0 in the climatology, 11 ensemble members are 0-value and the remaining 40 are greater than 0. The average rank for the 0-value members are 5.5 (as this is given by the method of handling the 0-value issue with equal representation, explained above), while the average rank for the non-zero members is given as an example of 11. The related rank-mean is then 9.81, making this forecast into the 'Extreme low' expected category while the rank-std is 2.26, with low uncertainty category. These tables demonstrate the complex interaction between the dryness of the climatology and ensemble forecasts, reflected in the forecast rank-mean and rank-std values and the subsequent expected anomaly and uncertainty categories. They also demonstrate, how less likely it becomes to have negative anomalies as the climatology becomes drier and drier. Expand | | ||||||
| |||||||
Number of 0 members | Number of non-0 members | Average rank of 0 members | Average rank of non-0 members | Rank-mean | Forecast anomaly category | Rank-std | Uncertainty category |
0 | 51 | NA (no member to rank) | 11 (the lowest possible rank for a non-zero member if 1-10 percentiles in the climatology are 0) | (0 * 5.5 + 51 * 11)/51 = 11 | Low (10-25) | 0 | Low uncertainty |
0 | 51 | NA | 20 | (0 * 5.5 + 51 * 20)/51 = 20 | Low (10-25) | 0 | Low uncertainty |
0 | 51 | NA | 50 | (0 * 5.5 + 51 * 50)/51 = 50 | Near normal (40-60) | 0 | Low uncertainty |
0 | 51 | NA | 70 | (0 * 5.5 + 51 * 70)/51 = 70 | Bit high (60-75) | 0 | Low uncertainty |
0 | 51 | NA | 100 | (0 * 5.5 + 51 * 100)/51 = 100 | Extreme high (90<) | 0 | Low uncertainty |
11 | 40 | 5.5 | 11 | (11 * 5.5 + 40 * 11)/51 = 9.81 | Extreme low (<10) | 2.26 | Low uncertainty |
11 | 40 | 5.5 | 20 | (11 * 5.5 + 40 * 20)/51 = 16.87 | Low (10-25) | 5.96 | Low uncertainty |
11 | 40 | 5.5 | 50 | (11 * 5.5 + 40 * 50)/51 = 40.40 | Near normal (40-60) | 18.3 | Medium uncertainty |
11 | 40 | 5.5 | 70 | (11 * 5.5 + 40 * 70)/51 = 56.08 | Near normal (40-60) | 26.52 | High uncertainty |
11 | 40 | 5.5 | 100 | (11 * 5.5 + 40 * 100)/51 = 79.61 | High (75-90) | 38.86 | High uncertainty |
21 | 30 | 5.5 | 11 | (21 * 5.5 + 30 * 11)/51 = 8.73 | Extreme low (<10) | 2.70 | Low uncertainty |
21 | 30 | 5.5 | 20 | (21 * 5.5 + 30 * 20)/51 = 14.02 | Low (10-25) | 7.13 | Low uncertainty |
21 | 30 | 5.5 | 50 | (21 * 5.5 + 30 * 50)/51 = 31.67 | Bit low (25-40) | 21.90 | High uncertainty |
21 | 30 | 5.5 | 70 | (21 * 5.5 + 30 * 70)/51 = 43.44 | Near normal (40-60) | 31.74 | High uncertainty |
21 | 30 | 5.5 | 100 | (21 * 5.5 + 30 * 50)/51 = 61.08 | Bit high (60-75) | 46.50 | High uncertainty |
36 | 15 | 5.5 | 11 | (36 * 5.5 + 15 * 11)/51 = 7.11 | Extreme low (<10) | 2.50 | Low uncertainty |
36 | 15 | 5.5 | 20 | (36 * 5.5 + 15 * 20)/51 = 9.76 | Extreme low (<10) | 6.60 | Low uncertainty |
36 | 15 | 5.5 | 50 | (36 * 5.5 + 15 * 50)/51 = 18.58 | Low (10-25) | 20.27 | High uncertainty |
36 | 15 | 5.5 | 70 | (36 * 5.5 + 15 * 70)/51 = 24.47 | Low (10-25) | 29.38 | High uncertainty |
36 | 15 | 5.5 | 100 | (36 * 5.5 + 15 * 50)/51 = 33.29 | Bit low (25-40) | 43.05 | High uncertainty |
51 | 0 | 5.5 | NA (no member to rank) | (51 * 5.5 + 0)/51 = 5.5 | Extreme low (<10) | 0 | Low uncertainty |
Expand | |||||||
| |||||||
Number of 0 members | Number of non-0 members | Average rank of 0 members | Average rank of non-0 members | Rank-mean | Forecast anomaly category | Rank-std | Uncertainty category |
0 | 51 | NA (no member to rank) | 31 (the lowest possible rank for a non-zero member if 1-30 percentiles in the climatology are 0) | (0 * 15.5 + 51 * 31)/51 = 31 | Bit low (25-40) | 0 | Low uncertainty |
0 | 51 | NA | 50 | (0 * 15.5 + 51 * 50)/51 = 50 | Near normal (40-60) | 0 | Low uncertainty |
0 | 51 | NA | 70 | (0 * 15.5 + 51 * 70)/51 = 70 | Bit high (60-75) | 0 | Low uncertainty |
0 | 51 | NA | 100 | (0 * 15.5 + 51 * 100)/51 = 100 | Extreme high (90<) | 0 | Low uncertainty |
11 | 40 | 15.5 | 31 | (11 * 15.5 + 40 * 31)/51 = 27.65 | Bit low (25-40) | 6.37 | Low uncertainty |
11 | 40 | 15.5 | 50 | (11 * 15.5 + 40 * 50)/51 = 42.55 | Near normal (40-60) | 14.18 | Medium uncertainty |
11 | 40 | 15.5 | 70 | (11 * 15.5 + 40 * 70)/51 = 58.24 | Near normal (40-60) | 22.41 | High uncertainty |
11 | 40 | 15.5 | 100 | (11 * 15.5 + 40 * 100)/51 = 81.77 | High (75-90) | 34.75 | High uncertainty |
21 | 30 | 15.5 | 31 | (21 * 15.5 + 30 * 31)/51 = 24.61 | Low (10-25) | 7.62 | Low uncertainty |
21 | 30 | 15.5 | 50 | (21 * 15.5 + 30 * 50)/51 = 35.79 | Bit low (25-40) | 16.97 | Medium uncertainty |
21 | 30 | 15.5 | 70 | (21 * 15.5 + 30 * 70)/51 = 47.55 | Near normal (40-60) | 26.82 | High uncertainty |
21 | 30 | 15.5 | 100 | (21 * 15.5 + 30 * 100)/51 = 65.20 | Bit high (60-75) | 41.58 | High uncertainty |
36 | 15 | 15.5 | 31 | (36 * 15.5 + 15 * 31)/51 = 20.05 | Low (10-25) | 7.06 | Low uncertainty |
36 | 15 | 15.5 | 50 | (36 * 15.5 + 15 * 50)/51 = 25.64 | Bit low (25-40) | 15.71 | Medium uncertainty |
36 | 15 | 15.5 | 70 | (36 * 15.5 + 15 * 70)/51 = 31.52 | Bit low (25-40) | 24.83 | High uncertainty |
36 | 15 | 15.5 | 100 | (36 * 15.5 + 15 * 100)/51 = 40.35 | Near normal (40-60) | 38.50 | High uncertainty |
51 | 0 | 15.5 | NA (no member to rank) | (51 * 15.5 + 0)/51 = 15.5 | Low (10-25) | 0 | Low uncertainty |
Expand | |||||||
| |||||||
Number of 0 members | Number of non-0 members | Average rank of 0 members | Average rank of non-0 members | Rank-mean | Forecast anomaly category | Rank-std | Uncertainty 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 | Low uncertainty |
0 | 51 | NA | 100 | (0 * 35.5 + 51 * 100)/51 = 100 | Extreme high (90<) | 0 | Low uncertainty |
11 | 40 | 35.5 | 71 | (11 * 35.5 + 40 * 71)/51 = 63.34 | Bit high (60-75) | 14.60 | Medium uncertainty |
11 | 40 | 35.5 | 100 | (11 * 35.5 + 40 * 100)/51 = 86.08 | High (75-90) | 26.52 | High uncertainty |
21 | 30 | 35.5 | 71 | (21 * 35.5 + 30 * 71)/51 = 56.38 | Near normal (40-60) | 17.47 | Medium uncertainty |
21 | 30 | 35.5 | 100 | (21 * 35.5 + 30 * 100)/51 = 73.44 | Bit high (60-75) | 31.74 | High uncertainty |
36 | 15 | 35.5 | 71 | (36 * 35.5 + 15 * 71)/51 = 45.94 | Near normal (40-60) | 16.17 | Medium uncertainty |
36 | 15 | 35.5 | 100 | (36 * 35.5 + 15 * 100)/51 = 54.47 | Near normal (40-60) | 29.38 | High uncertainty |
51 | 0 | 35.5 | NA (no member to rank) | (51 * 35.5 + 0)/51 = 35.5 | Bit low (25-40) | 0 | Low uncertainty |
Expand | |||||||
| |||||||
Number of 0 members | Number of non-0 members | Average rank of 0 members | Average rank of non-0 members | Rank-mean | Forecast anomaly category | Rank-std | Uncertainty 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<) | 0 | Low uncertainty |
11 | 40 | 50.5 | 100 | (11 * 50.5 + 40 * 100)/51 = 89.32 | High (75-90) | 20.35 | High uncertainty |
21 | 30 | 35.5 | 100 | (21 * 50.5 + 30 * 100)/51 = 79.61 | High (75-90) | 24.36 | High uncertainty |
36 | 15 | 35.5 | 100 | (36 * 50.5 + 15 * 100)/51 = 65.05 | Bit high (60-75) | 22.55 | High uncertainty |
51 | 0 | 35.5 | NA (no member to rank) | (51 * 50.5 + 0)/51 = 50.5 | Near normal (40-60) | 0 | Low uncertainty |
Example No-1: No zero section in climatology, 5 number/rank groups:
In these simplified examples, a fixed portion of the climatological and/or ensemble forecast distribution is 0, while the non-zero ensemble members have just one other rank for the very dry cases and 5 for the non-zero cases, with members having the same rank in each group. This way, the computation methodology can be demonstrated in a simple way that is easier to interpret.
in a Below, there are examples with simplified rank distributions and specific cases of very dry rivers (0 values) that will demonstrate how the forecast anomaly category and uncertainty category generation work.
For example, when the X% of the climatological distribution is 0, then the average rank of the 0-value ensemble members will always be X/2+0.5 with the even rank representation for the 0-value case explained above (e.g. for 10% being 0 in the climate, the average rank of the 0-value forecast ensemble members will be 5.5).
Example No-2: Lowest 10% of the climatology is 0, 2 number/rank groups:
less likely it becomes to have negative anomalies as the climate becomes drier and drier.
|
...
|
...
|
...
|