Design Rainfall Data and Analysis
Design Rainfall Data and Analysis
Obs. withObs.
missing
without
datamissing
Minimum
data Maximum Mean Std. deviation
0 28 26.700 115.000 48.607 17.711
G(Critical value)p-value Step
2.876 0.000 1
Z-scores
4
1
Z-score
0
1
-1
-2
-3
-4
Observations
Outlier test for Deg
Summary statistics:
Variable
Std. deviation Var1
Grubbs test for outliers / Two
G (Observed
G (Critical
p-value (Tw
alpha
Test interpretation:
H0: There is no outlier in the
Ha: The minimum or maximu
As the computed p-value is lo
null hypothesis H0, and accep
Outliers:
Var1
110.000
Z-scores:
Value
34.200
41.800
110.000
45.000
83.500
52.000
25.000
58.500
40.200
42.500
45.100
54.200
58.000
48.100
43.300
65.000
30.000
24.000
36.000
50.800
57.200
39.100
55.500
50.000
40.200
40.200
Values displayed in bold are o
Z-score
0
1 2 3 4 5 6 7 8 9
-1
-2
-3
-4
Outlier test for Degahabur
XLSTAT 2019.2.2.59614 - Grubbs test for outliers - Start time: 5/10/2023 at 2:33:19 PM / End time: 5/10/2023 at 2:33:20 PM / Microsoft
Data: Workbook = Design rainfall data and analysis.xlsx / Sheet = Max. Daily RF of stations / Range = 'Max. Daily RF of stations'!$G$3:$G
Alternative hypothesis: Two-sided
Significance level (%): 5
Iterations: Maximum: 1
Summary statistics:
Observations
Obs. withObs.
missing
without
datamissing
Minimum
data Maximum Mean Std. deviation
26 0 26 24.000 110.000 48.823 17.806
Grubbs test for outliers / Two-tailed test:
3.436
2.841
0.002
0.05
Test interpretation:
H0: There is no outlier in the data
Ha: The minimum or maximum value is an outlier
As the computed p-value is lower than the significance level alpha=0.05, one should reject the
null hypothesis H0, and accept the alternative hypothesis Ha.
Z-score
-0.821
-0.394
3.436
-0.215
1.947
0.178
-1.338
0.543
-0.484
-0.355
-0.209
0.302
0.515
-0.041
-0.310
0.909
-1.057
-1.394
-0.720
0.111
0.470
-0.546
0.375
0.066
-0.484
-0.484
Values displayed in bold are outliers
Z-scores
4
1
Z-score
0
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26
-1
-2
-3
-4
Observations
023 at 2:33:20 PM / Microsoft Excel 16.014332 XLSTAT 2019.2.2.59614 - Grubbs test for outliers - Start time
Daily RF of stations'!$G$3:$G$28 / 26 rows and 1 column Data: Workbook = Design rainfall data and analysis.xlsx / She
Alternative hypothesis: Two-sided
Significance level (%): 5
Iterations: Maximum: 1
Summary statistics:
Variable Observations
Obs. withObs.
missing
without
datamissing data
Var1 21 0 21
Grubbs test for outliers / Two-tailed test:
G (Observed 2.416
G (Critical 2.734
p-value (Tw 0.194
alpha 0.05
Test interpretation:
H0: There is no outlier in the data
Ha: The minimum or maximum value is an outlier
As the computed p-value is greater than the significance leve
null hypothesis H0.
Z-scores:
Value Z-score
34.500 -0.932
42.200 -0.574
33.700 -0.969
67.100 0.582
106.600 2.416
39.000 -0.723
70.800 0.754
31.300 -1.081
47.300 -0.338
64.500 0.461
35.800 -0.872
36.000 -0.862
41.500 -0.607
62.000 0.345
45.500 -0.421
49.800 -0.222
42.300 -0.570
87.600 1.534
49.600 -0.231
60.900 0.294
98.000 2.017
Values displayed in bold are outliers
Z-scores
3
1
Z-scores
3
Z-score
0
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20
-1
-2
-3
Observations
bs test for outliers - Start time: 5/10/2023 at 3:03:30 PM / End time: 5/10/2023 at 3:03:31 PM / Microsoft Excel 16.014332
ll data and analysis.xlsx / Sheet = Max. Daily RF of stations / Range = 'Max. Daily RF of stations'!$K$3:$K$23 / 21 rows and 1 column
cores
cores
10 11 12 13 14 15 16 17 18 19 20 21
servations
cel 16.014332
3 / 21 rows and 1 column
Station Name :Jigjiga
1 Annual Maximum Daily Rainfall Series and Quality of Data
SUM 1246.00
MEAN 46.15
STANDARD DEVIATION 12.24
SKEWNESS COEFFICIENT 0.642
Km 3.25 79.10
Gumbel-Powell Method
return period T yrs 20.00
Standard variate, Sn 1.10
standard mean, Yn 0.53
yt 2.97
frequency factor , KT 1.87
Y=Xmean + Z*XStandared deviation 68.99
Y=Ymean + Z* y σ 1.838
Design Rainfall, X50 = 10Y 68.80
2.5 Normal
Design Period, T 20.00
Probability, P 0.05
K=(Cs/6) 0.000
W=(Ln(1/P )) 2 0.5
2.45
KT=(w((2.515517+0.802853*w+0.010328*w )/(1+1.432788*w+0.189269*w +0.001308*w )))
2 2 3
1.65
X50=Xmean + KT*σx 66.29
2.6 D-indexes Take at least above Xmean for the value of XI.
Normal
Rank XI PR w
1 79.10 0.032 2.6207
2 62.00 0.065 2.3413
3 61.90 0.097 2.1612
4 61.20 0.129 2.0237
5 59.00 0.161 1.9103
6 54.90 0.194 1.8123
Sum
Sum/Xm
Log Pearson
Normal Type III
Rank XI XI -'XI' XI -'XI'
1 79.10 10.312 6.497
2 62.00 2.738 4.548
3 61.90 0.170 0.939
4 61.20 1.203 1.100
5 59.00 0.741 1.103
6 54.90 1.836 1.133
Sum 17.000 15.319
Sum/Mean 0.368 0.332
Point Rainfall 66.29 68.94
Design Point Rainfall = 73.355
Upper confidence limit(Pu)@95%CL 87.454
Lower confidence limit(Pl)@95%CL 59.255
Logarithmic plotting
Rank Value/Yo/ (Yo-Ym)2 (Yo-Ym)3 position
1 1.8982 0.0617342 0.0153387 0.04
2 1.7924 0.0203573 0.0029046 0.07
3 1.7917 0.0201577 0.0028620 0.11
4 1.7868 0.0187796 0.0025735 0.14
5 1.7709 0.0146747 0.0017777 0.18
6 1.7396 0.0080747 0.0007256 0.21
7 1.7340 0.0071042 0.0005988 0.25
7 1.7340 0.0071042 0.0005988 0.25
9 1.7300 0.0064419 0.0005170 0.32
10 1.7059 0.0031529 0.0001770 0.36
11 1.6812 0.0009940 0.0000313 0.39
12 1.6739 0.0005871 0.0000142 0.43
13 1.6665 0.0002824 0.0000047 0.46
14 1.6532 0.0000122 0.0000000 0.50
15 1.6385 0.0001260 -0.0000014 0.54
16 1.6064 0.0018776 -0.0000814 0.57
17 1.6031 0.0021686 -0.0001010 0.61
18 1.6021 0.0022708 -0.0001082 0.64
19 1.6010 0.0023756 -0.0001158 0.68
20 1.5798 0.0048901 -0.0003420 0.71
21 1.5740 0.0057277 -0.0004335 0.75
22 1.5599 0.0080651 -0.0007243 0.79
23 1.5079 0.0201234 -0.0028546 0.82
24 1.5065 0.0205085 -0.0029370 0.86
25 1.4955 0.0237679 -0.0036643 0.89
26 1.4829 0.0278353 -0.0046440 0.93
27 1.4265 0.0498189 -0.0111197 0.96
44.54 0.34 0.0010
1.65
0.11 0.1142
0.028 0.0278
44.88 10.53 0.97 0.86
80
70
60
50
40
80
70
60
ok
50
40
30
20
10
0
0.01 0.10 1.00
25 100
83.84 99.67
79.23 92.86
79.72 94.84
76.76 86.97
79.86 93.65
75.58 82.76
Remark:The value of b =3/8 for all distributions except Gumbel whose b=11/25
But in this case m/(N+1) probability plotting equation is used
Normal Log Nor
KT XI' XI-'XI' PR w
1.849 68.79 10.312 0.032 2.62068205
1.518 64.74 2.738 0.065 2.341298795
1.300 62.07 0.170 0.097 2.161191762
1.131 60.00 1.203 0.129 2.023705929
0.989 58.26 0.741 0.161 1.910261392
0.865 56.74 1.836 0.194 1.812306671
17.000
0.36838
One-tail
a = 0.005
a = 0.010
a = 0.025
a = 0.050
a = 0.100
1.00
or all distributions except Gumbel whose b=11/25
ability plotting equation is used
Log Normal
KT Y XI' XI-'XI' PR
1.84900948637617 1.86084787314786 72.5851657935149 6.515 0.032
1.5182361346521 1.82307744983464 66.5391808159651 4.539 0.065
1.30034247479262 1.79819655865261 62.834267798039 0.934 0.097
1.1310494079582 1.77886528312814 60.0987284076957 1.101 0.129
0.989128202840989 1.76265955202417 57.8974654213409 1.103 0.161
0.864749550378005 1.74845697394414 56.0346900399631 1.135 0.194
15.327
0.33212
CALCULATION SUMMARY
DFs 20-year Design RainfD-index statistic
Normal 66.29 0.37
Log Normal 68.80 0.33
Pearson Type III 68.25 0.36
Log Pearson Type III 68.94 0.33
Gumbel EVI 73.35 0.37
Gumbel-Powell 68.99 0.51
Ven Te Chow 99.27 1.48
Stochastic method 195.29 7.85
obs lag xi(xi+1) xi2 xi3 xi4
53.7 3324.03 3831.61 237176.659 14681235.1921
50.8 2727.96 2883.69 154854.153 8315668.0161
36.3 1844.04 2580.64 131096.512 6659702.8096
46.4 1684.32 1317.69 47832.147 1736306.9361
54.9 2547.36 2152.96 99897.344 4635236.7616
39.9 2190.51 3014.01 165469.149 9084256.2801
37.5 1496.25 1592.01 63521.199 2534495.8401
62 2325 1406.25 52734.375 1977539.0625
54.2 3360.4 3844 238328 14776336
31.3 1696.46 2937.64 159220.088 8629728.7696
38 1189.4 979.69 30664.297 959792.4961
40 1520 1444 54872 2085136
61.9 2476 1600 64000 2560000
40.4 2383.6 3481 205379 12117361
61.2 2472.48 1632.16 65939.264 2663946.2656
43.5 2662.2 3745.44 229220.928 14028320.7936
45 1957.5 1892.25 82312.875 3580610.0625
26.7 1201.5 2025 91125 4100625
48 1281.6 712.89 19034.163 508212.1521
#REF! #REF! 2304 110592 5308416
30.4 #REF! #REF! #REF! #REF!
40.1 1219.04 924.16 28094.464 854071.7056
47.2 1892.72 1608.01 64481.201 2585696.1601
54.2 2558.24 2227.84 105154.048 4963271.0656
0 0 2937.64 159220.088 8629728.7696
#REF! #REF! #REF! #REF!
R S2 S3 S4
R(bar) #REF!
var('R) #REF!
std #REF!
u #REF!
Two-tail
2.576 a = 0.01
2.326
1.960 a = 0.05 a 5%
1.645 a = 0.10 a/2 0.025
1.282 1-a/2 0.975
Since |u| =0.693 is less than u0.025 = 1.96, the data show independent and indicates no trend at 5% significance level.
Validation @
acceptable (5%)
acceptable (5%)
acceptable (5%)
no higher outlier (10%)
no lower outlier(10%)
acceptable (5%)
Log person Type III
w K Z KT y XI'
2.6207 0.0004 1.8490 1.8499 1.8610 72.6026
2.3413 0.0004 1.5182 1.5187 1.8231 66.5478
2.1612 0.0004 1.3003 1.3006 1.7982 62.8386
2.0237 0.0004 1.1310 1.1312 1.7789 60.1004
1.9103 0.0004 0.9891 0.9891 1.7627 57.8973
1.8123 0.0004 0.8647 0.8647 1.7484 56.0333
1.4 test for data trend absence using M-K
M-K test
Time obs
1 1989 61.9
2 1990 53.7
3 1991 50.8
4 1992 36.3
5 1993 46.4
6 1994 54.9
7 1995 39.9
8 1996 37.5
9 1997 62
10 1998 54.2
11 2000 31.3
12 2001 38
13 2002 40
14 2003 59
15 2004 40.4
16 2005 61.2
17 2006 43.5
18 2007 45
19 2008 26.7
20 2009 48
21 #REF! #REF!
22 2011 30.4
23 2012 40.1
24 2013 47.2
n S
13 Variance
z
One-tail
Two-tail a = 0.005
ucrit a = 0.010
1.96 a = 0.025
a = 0.050
a = 0.100
Since z =11.958 is less than z0.025 = 1.96, the
Person Type III
XI-'XI' PR w K Z
6.497 0.032 2.620682050339 0.019 1.849009486
4.548 0.065 2.341298795082 0.019 1.518236135
0.939 0.097 2.161191761884 0.019 1.300342475
1.100 0.129 2.023705928916 0.019 1.131049408
1.103 0.161 1.910261391565 0.019 0.989128203
1.133 0.194 1.812306671211 0.019 0.86474955
15.319
0.33196
absence using M-K
S
#REF!
#REF!
#REF!
#REF!
#REF!
#REF!
#REF!
#REF!
#REF!
#REF!
#REF!
#REF!
#REF!
#REF!
#REF!
#REF!
#REF!
#REF!
#REF!
#REF!
#REF!
-4
-5
-6
#REF!
268.7
#REF!
Two-tail
2.576 a = 0.01
2.326
1.960 a = 0.05
1.645 a = 0.10
1.282
Since z =11.958 is less than z0.025 = 1.96, the data has trend at 5% significance level.
1
YT -3.314076
QYT -262.1434
YT2 10.983098
Q 79.1
erson Type III Gumbel-Powell Gumbel
KT XI' XI-'XI' PR KT XI' XI-'XI' PR
1.89445062428146 69.34452602 9.755 0.042 2.011 70.78 8.324 0.032
1.54238677894908 65.03372167 3.034 0.083 1.454 63.95 1.949 0.065
1.31280998427967 62.22269561 0.323 0.125 1.120 59.86 2.040 0.097
1.13571794749214 60.05431299 1.146 0.167 0.877 56.89 4.314 0.129
0.9881172466782 58.24703357 0.753 0.208 0.684 54.52 4.480 0.161
0.859403193276086 56.67100948 1.771 0.250 0.521 52.53 2.368 0.194
16.782 23.475
0.36365 0.50869
1.5 test for data homogeneity and stationarity using M-W
M-W test
26
27
28
29
R 91
p 5
q 8
N 13
V 76
W -36
U -36
U(bar) 20 Two-tail
sumT 0 a 5% ucrit
a/2 0.025 1.96
var(U) 46.666666667 1-a/2 0.975
|u| 8.20
ucrit @5% 1.96
Since |u| =8.07 is not less than u0.025 = 1.96, the data is not considered homogeneous and stationary at 5% significance level.
2 3 4 5 6 7 8 9 10
-2.602232 -2.177463 -1.869825 -1.62602322 -1.422286 -1.24589932 -1.245899 -0.9473544235 -0.816824
-161.3384 -134.785 -114.4333 -95.93536999 -78.08351 -67.5277433 -67.52774 -50.872932544 -41.49465
6.7716122 4.741345 3.4962445 2.6439515126 2.0228979 1.552265125 1.5522651 0.8974804038 0.6672012
62 61.9 61.2 59 54.9 54.2 54.2 53.7 50.8
a -2.6147
b -26.11145
Gumbel EVI ven Te Chow methoStochastic method PMF
YT KT XI' XI-'XI' Xsim Xsim Xobs - Xsim
3.418 2.6296362 78.35 0.754 53.981225 25.12 172.43707 93.34
2.708 1.9827843 70.43 8.426 53.096007 8.90 141.43208 79.43
2.285 1.5975978 65.71 3.810 52.179721 9.72 123.29533 61.40
1.979 1.3192503 62.30 1.102 51.230109 9.97 110.4271 49.23
1.738 1.0991981 59.61 0.607 50.244655 8.76 100.44572 41.45
1.537 0.9157969 57.36 2.462 49.220545 5.68 92.290343 37.39
17.160 68.148 362.228
0.37184 99.270578 1.47672 195.28668 7.84923 254.30265
1.6 test for data trend absence using Spearman's Rank-Correlation Method
tleft -2.05
t=tright 2.05
Since tl=1.257 is less than t0.975 = 2.052 & greater than t0.025=-2.052 , the data has no trend at 5% significance level.
11 12 13 14 15 16 17 18 19 20
-0.695167 -0.580505 -0.471358 -0.366513 -0.264936 -0.165703 -0.067948 0.0291892 0.126615 0.2253515
-33.36801 -27.39983 -21.871 -16.49308 -11.52471 -6.6944 -2.724707 1.1675694 5.0519373 8.5633565
0.483257 0.3369859 0.222178 0.1343317 0.070191 0.0274575 0.0046169 0.000852 0.0160314 0.0507833
48 47.2 46.4 45 43.5 40.4 40.1 40 39.9 38
d at 5% significance level.
21 22 23 24 25 26 27 28 29 30
0.3266343 0.4320714 0.5439315 0.6657298 0.8036112 0.9704218 1.2036341
12.248785 15.68419 17.514594 21.369927 25.153029 29.500822 32.13703
0.1066899 0.1866857 0.2958615 0.4431962 0.6457909 0.9417184 1.448735 0 0 0
37.5 36.3 32.2 32.1 31.3 30.4 26.7
sum
#N/A #N/A #N/A #N/A #N/A #N/A #N/A #N/A #N/A -14.55282
#N/A #N/A #N/A #N/A #N/A #N/A #N/A #N/A #N/A -1025.826
#N/A #N/A #N/A #N/A #N/A #N/A #N/A #N/A #N/A 40.743724
#N/A #N/A #N/A #N/A #N/A #N/A #N/A #N/A #N/A 1246
XLSTAT 2019.2.2.59614 - Grubbs test for outliers - Start time: 18/09/2023 at 16:43:25 / End time: 18/09/2023 at 16
Data: Workbook = Design rainfall data and analysis .xlsx / Sheet = Design Rainfall_Jigjiga / Range = 'Design Rainfall_J
Alternative hypothesis: Two-sided
Significance level (%): 5
Iterations: Maximum: 1
Summary statistics:
Variable Observations
Obs. withObs.
missing
without
datamissing
Minimum
data Maximum Mean Std. deviation
Var1 27 0 27 26.700 79.100 46.148 12.244
G (Observed 2.691
G (Critical 2.859
p-value (Tw 0.103
alpha 0.05
Test interpretation:
H0: There is no outlier in the data
Ha: The minimum or maximum value is an outlier
As the computed p-value is greater than the significance level alpha=0.05, one cannot reject the
null hypothesis H0.
Z-scores:
Value Z-score
61.900 1.286
53.700 0.617
50.800 0.380
36.300 -0.804
46.400 0.021
54.900 0.715
39.900 -0.510
37.500 -0.706
62.000 1.295
54.200 0.658
31.300 -1.213
38.000 -0.665
40.000 -0.502
59.000 1.050
40.400 -0.469
61.200 1.229
43.500 -0.216
45.000 -0.094
26.700 -1.588
48.000 0.151
30.400 -1.286
40.100 -0.494
47.200 0.086
54.200 0.658
32.100 -1.147
32.200 -1.139
79.100 2.691
Values displayed in bold are outliers
Z-scores
3
1
Z-score
0
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27
-1
-2
-3
Observations
/ End time: 18/09/2023 at 16:43:28 / Microsoft Excel 16.05173
ga / Range = 'Design Rainfall_Jigjiga'!$D$6:$D$32 / 27 rows and 1 column
Thieson polygon
Area of catchmnt in km2 under the M.stations
Jigjiga Degahabur Gode
Km2 4280 0 0
% 1 0 0
Design RF from frequency analysis (mm) 68.94 0.00 0.00
Design Rf for the whole catchment (mm) 68.94068
Area reduction factor 0.5613686050758
Design rainfall after area reduction factor (mm) 38.701133057505
768
0.00512
0.512