Next Article in Journal
Bioremediation of Waste Water to Remove Heavy Metals Using the Spent Mushroom Substrate of Agaricus bisporus
Previous Article in Journal
The Threshold Effect of Environmental Regulation, FDI Agglomeration, and Water Utilization Efficiency under “Double Control Actions”—An Empirical Test Based on Yangtze River Economic Belt
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Statistical Study of Rainfall Control: The Dagum Distribution and Applicability to the Southwest of Spain

1
Engineering Projects Department, Universidad de Extremadura, Avda. Elvas s/n, 06071 Badajoz, Spain
2
Graphical Expression Department, Universidad de Extremadura, Avda. Elvas s/n, 06071 Badajoz, Spain
3
Infrastructure Department, Diputación de Badajoz, 06002 Badajoz, Spain
*
Author to whom correspondence should be addressed.
Water 2019, 11(3), 453; https://doi.org/10.3390/w11030453
Submission received: 1 January 2019 / Revised: 19 February 2019 / Accepted: 22 February 2019 / Published: 4 March 2019
(This article belongs to the Section Hydrology)

Abstract

:
It is of vital importance in statistical distributions to fit rainfall data to determine the maximum amount of rainfall expected for a specific hydraulic work. Otherwise, the hydraulic capacity study could be erroneous, with the tragic consequences that this would entail. This study aims to present the Dagum distribution as a new statistical tool to calculate rainfall in front of frequent statistical distributions such as Gumbel, Log-Pearson Type III, Gen Extreme Value (GEV) and SQRT-ET max. The study was performed by collecting annual rainfall data from 52 meteorological stations in the province of Badajoz (Spain), using the statistical goodness-of-fit tests of Anderson–Darling and Kolmogorov–Smirnov to establish the degree of fitness of the Dagum distribution, applied to the maximum annual rainfall series. The results show that this distribution obtained a flow 21.92% greater than that with the traditional distributions. Therefore, in the Southwest of Spain, the Dagum distribution fits better to the observed rainfall data than other common statistical distributions, with respect to precision and calculus of hydraulics works and river flood plains.

1. Introduction

River flooding and hydrologic studies are carried out ensuring that the waters will reach a certain maximum level during statistic rainfalls with a certain return period. The possible level of flooding and the waterworks are designed by means of different classical statistical distributions applied to rainfall, using a series of maximum annual recorded rainfall data.
Some models of future scenarios suggest that climate change will involve a significant modification in the distribution of extreme rainfall intensity [1].
There have been several significant flooding events in the southwest of Spain in the last few decades. Particularly in Badajoz (Spain) in the year 1997, there were even casualties as a result of the river Rivillas breaking its banks. Therefore, undoubtedly another solution must be sought to provide a better fit to the historical rainfall data than those options currently available [2].
According to the scientific literature the statistical distributions most commonly used in Europe and Spain are: The Gumbel distribution, developed by the German mathematician Gumbel [3] and later applied to hydrology [4]; the Log-Pearson Type III distribution [5,6] put forward by several authors for use in hydrology [7,8]; extreme values distribution (GEV) [9]; and the SQRT-ET max. distribution [10,11], which best fits the characteristics of Spanish rainfall.
The method most commonly used in Spain for the calculation of the flood flow (m3/s) for a specific return period with the corresponding calculation changes is regulated by Spanish Highway Instruction Norm 5.2-IC [12]. For the calculation of the flow this regulation proposes the use of the Rational Method [13], which is based on the fact that the duration of the rainfall is greater than or equal to the concentration time in the basin and, therefore, it determines maximum flow rate. The rational method is used in Spain for small basins. For large ones, with concentration times of over 6 hours or more than 1 km2 (and up to 3000 km2), the modified method [8] was used, referred to by hydrologists such as Ferrer (1993) [14]. Currently, with the development of the computational capacity, more accurate hydro-meteorological models are widely used (e.g., software Hydrologic Modeling System, HEC-HMS). In other words, both models are in use depending on the basin size or predictable concentration times of more than 6 hours.
Panagoulia and Caroni (2011) [15] investigated the fit of generalized extreme value (GEV) distributions to maximum precipitation over the Mesochora mountainous catchment in central-western Greece. They observed that the stationary model was adequate for historical data on annual daily maxima for 1972–1992 and also for 1 × CO2 for the period 1961–2000. However, the 2 × CO2 series for 2061–2100 required a cubic time trend in location to obtain a satisfactory fit.
Lang et al. (1999) [16] reviewed tests and methods useful for modeling the process of over-threshold values, the choice of the threshold level, the verification of the independence of the values and the stationarity of the process. The independence criteria for the extraction of the peak-over-threshold events were recently modified by Onyutha (2017, 2019) [17,18] in terms of the independency ratio, threshold, and the inter-event time. Despite the desirable properties of the peak-over-threshold in the analyses of extreme events, this study focused on the annual maximum series.
It must be pointed out that these distributions are limited, as they must rely on rainfall for the desired return period. Additionally, there are scarce statistical distributions currently used in hydrology. They only were conceived at the beginning of the last century and scientific methods have, without doubt, moved on since then. For example, it seems that: (i) the distribution of Gumbel in high return periods underestimates precipitation values, and creates uncertainty about the safety of works designed with this criterion; (ii) the Log-Pearson Type III distribution does not underestimate values, but has three estimable parameters; (iii) the SQRT-ET max distribution is perhaps more precise and easier to apply, although it also has estimable parameters; (iv) the GEV distribution takes into account the extreme values maximum annual rainfall of one day but not the intermediate values. Other distributions, though not considered in this study, are Pareto [19] and Generalized Logistic [20].
However, the Dagum distribution [21] presents all the requisites to be used in hydrological studies. For example, it is a continuous distribution of extreme and non-negative values of the continuous variable. The Dagum distribution is used in different scientific fields, such as econometrics, economic theory (particularly for model ranges as diverse as wealth distribution), the banking sector (by using the distribution methods of aggregate losses) and the insurance sector (applied to analyze solvency) [22].
This distribution began to be used mostly in latitudes different from Europe for the calculation of rainfall, Mielke and Johnson [23] verified that a distribution similar to that of Dagum is possible in meteorology and hydrology. Later we found studies by Alam S. et al. (2014) [24], Alam and Rahat (2015) [25] in rivers of Dhaka (Bangladesh), Crabbe (2014) [26] in climate change systems in Charlottetown (Canada) and Mayooran and Laheetharan (2014) [27] in Colombo (Sri Lankan), possibly as a consequence of the existing problems of rainfall and floods in those countries.
On the other hand, in the East of Spain, Pérez-Sánchez and Senent-Aparicio (2018) [28], studied 29 meteorological stations between 1993 and 2013, for a given daily rainfall, simulating scenarios —for all those stations—with several distributions: Burr, Dagum, error, generalized extreme value, generalized logistic, generalized Pareto, Gumbel Max, inverse Gaussian, Johnson SB, Log-Logistic, Log-Pearson 3, Triangular, Weibull, and Wakeby. They observed that only the series of annual maximum dry spells offer a good adjustment for all the weather stations, showing that the Wakeby distribution presented the best results. Finally, Domma and Condino (2017) [29] carried out a simulation study that shows the good performance of the maximum likelihood estimators for finite samples.
However, Mayooran and Laheetharan (2014) [30] used the same form of Dagum as in this work and compared it to other 44 different distributions. The parameters of the selected probability distributions were used to generate random numbers for both actual and estimated maximum daily precipitation. Log-Pearson 3 and Burr (4P) were found to be the best-fit probability model for the annual period and the first inter monsoon study period, respectively [31].
The transmuted Dagum model provides the broader range of hazard behavior than the Dagum model [32]. The parameters of the new model are estimated by maximum likelihood using Newton Raphson approach and the information matrix and confidence intervals are also obtained. Other simulation results showed that both the corrected Akaike information criterion and Bayesian information criterion (BIC) always detected nonstationary, but the BIC selected the correct model more often except in very small samples [33]. Simulation studies indicated that the bias corrected and accelerated (BCa) method is best overall for the extreme percentiles that are often the focus of interest [34].
Despite all the above, it is not understandable that such a small number of distributions are used in professional practice in Europe and especially in Spain, considering that there are other efficient distributions in the field of hydrology.
To demonstrate the effectiveness of the Dagum distribution, the adjustment of the statistical distribution to the observed maximum annual rainfall values will be confirmed using the Anderson–Darling [35] and Kolmogorov–Smirnov goodness-of-fit tests and comparing them with the other distributions.
A statistical distribution must provide as good fit to the rainfall data as possible, since the better the fit, the more precise the value for the calculated rainfall. It can be used in the sizing of waterworks and flood plains.
This study intends to introduce the Dagum distribution as a new statistical tool to calculate rainfall because it fits better to the observed rainfall data by testing it with a dataset from Spain.

2. Materials and Methods

Firstly, a review of the statistical distributions used in hydrology studies will be performed, mainly in Europe and Spain, including the distribution of Dagum, with its fundamental characteristics.
The demonstration of the validity of this distribution in the field of civil engineering is analyzed to check whether the Dagum distribution provides a better fit, according to the goodness-of-fit tests, to the maximum annual rainfall distributions than the fit given by the commonly used distributions in the province of Badajoz, using real maximum annual rainfall data from the meteorological stations in that province.

2.1. Statistical Distribution Functions

The current method used in Spain to find out the flow rate that allows to dimension waterworks or calculate a flood plain of a river is the Rational Method [13], (except for large basins size). This method transforms the statistical rainfall associated to a certain return period (mm) to a flow rate (L/s):
Q = c × I × A
where c is a constant called runoff coefficient; I is the maximum intensity of precipitation and A is the area of the basin.
According to Instruction 5.2-IC, variable I can be obtained from the IDF curves (intensity, duration, frequency), with the following expression:
I i = P d 24 × ( I i I d ) 28 0 , 1 D 0 , 1 0 , 4
As it can be verified, the correct estimation of Pd is very important, since it is the maximum daily precipitation obtained through the series of daily rainfall recorded in rainfall stations.
The rainfall associated with the return period is currently calculated using the statistical distributions commonly used in hydrology: Gumbel, Log-Pearson Type III, SQRT-ET max and GEV. Thus, starting from a historical record obtained from the rainfall stations the maximum value of rainfall associated with a certain return period is determined (frequency).
The use of one of these statistical distributions is nearly always found in hydrological studies [12].

2.1.1. The Gumbel Distribution

According to Gumbel [36], the density function is:
F ( x ) = e e α ( x μ )
where x is the value of the random variable, F(x) represents the probability that rainfall is less than or equal to x, and α and μ are parameters of the fit which depend on the mean and standard deviation of the variable yi, which in turn depends on the sample size.
y i = L n   L n   ( N + 1 i ) α = S N S x μ = x ¯ y N ¯   S X S N
Being y ¯ and SN the mean and standard deviation of the variable yi, respectively, x ¯   and Sx the mean and standard deviation respectively of the sample from the data of daily maximum values of annual rainfall:
x ¯ = i = 1 N x i N S x = i = 1 N ( x x ¯ ) 2 N 1 2    y N ¯ = i = 1 N y i N    S N = i = 1 N ( y i y N ¯ ) 2 N 2
The return period, T(x), is related to the distribution function, F(x), by Equation (6):
T ( x ) = 1 1 F ( x )
After entering sample values, the analytical expression is reached Equation (7), after clearing x:
x = P 24 T = P d = μ 1 α [ L n ( L n T ( x ) 1 T ( x ) ) ]
where the expected daily maximum precipitation Pd is obtained for a given return period T(x).

2.1.2. The Log-Pearson Type III Distribution

The Log-Pearson density function is:
f ( x ) = λ β ( y ε ) β 1 e λ ( y ε ) x Γ ( β )
being y = log(x) and Γ(β) the function Gamma and e the Euler’s number.
β, λ and ε are the parameters of form, scale and position respectively and are drawn from the equation (9–11):
β = ( 2 C s ) 2
C s = n ( n 1 ) × ( n 2 ) × ( S log ( x ) ) 3 × ( log ( x ) log ( x ) ¯ ) 3
λ = S y β
The density function of this distribution is not integral, so it is resolved by parametric methods.
P 24 T = 10 log ( x ) + K × S log ( x )

2.1.3. The Distribution of SQRT-ET Max

The density function is as follows:
F ( x ) = exp [ k × ( 1 + α × x ) × exp ( α × x ) ]
where F(x) is the probability that the value will be less than x, and k and α are parameters to be estimated that depend on the mean and typical deviation of the data series.

2.1.4. The Gen Extreme Value (GEV)

The density function [10] is:
f ( x ) = { 1 σ e ( 1 + k z ) 1 / k ( 1 + k × z } ) 1 k × z , k 0 ,   f ( x ) = 1 α e e z , k = 0
where z = xμ/α, and k, μ and α are the parameters of form, position and scale.

2.1.5. The DAGUM Distribution

The Dagum distribution has long been used in different fields such as economics, econometrics [37] and social sciences. However, there are just a few applications found in hydrology. The importance of using this distribution in hydrology is both its adoptability to extreme data and similar capability to traditional distributions.
In probability theory, statistics and econometrics, the Dagum distribution is a continuous distribution with a probability distribution defined on real positive numbers. The Dagum distribution arose from several variations of a new model in the size distribution of personal incomes and is associated above all with the study of incomes. This distribution can be used for three parameters (Type I) and for four parameter (Type II). The density function is defined by:
F ( x ) = α k ( x γ β ) α k 1 β ( 1 + ( x γ β ) α ) k + 1
and the distribution function
f ( x ) = ( ( 1 + x y β ) α ) k
where k is a continuous shape parameter (k > 0), a is continuous shape parameter (a > 0), β is a continuous scale parameter (β > 0) and γ is a continuous location parameter (γ = 0 yields the three-parameter Dagum distribution) (γx ≤ ¥). Figure 1 shows the density function of Dagum.

2.2. Tests of Goodness

For the adjustment of the distributions the EasyFit software [38] was used, which adjusts the probability laws to the rainy series and allows performing the goodness-of-fit tests by the Kolmogorov–Smirnov, Anderson–Darling and chi-square methods.

2.2.1. The Kolmogorov–Smirnov test

The Kolmogorov–Smirnov test [39] is a nonparametric, single sample, bi-sample and continuous test that proves particularly useful for large samples and is therefore optimal for the study [40].
The Kolmogorov–Smirnov test considers two hypotheses:
H0: F(X) = Fs(X), H1: F(X) ≠ Fs(X)
where F(X) is the distribution function to be studied, and Fs(X) is the probability or theoretical proportion of values that must be less than or equal to x assuming the proposed hypothesis to be true.
Sample: n independent observations.
Contrast statistic: │Dn│ = maximum │Fs (X) − S (X) │

2.2.2. The Anderson–Darling Test

The Anderson–Darling test has been widely used in hydrology due to its reliability in comparison with other tests and its common use in samples with pronounced tails. This test is very interesting compared to commonly used tests when faced with a variety of hydraulic engineering alternatives [41].
The Anderson–Darling test [42] uses the following formulation as a test statistic:
A2 = −NS
S = i = 1 N ( 2 i 1 ) N L n ( F ( x ) i + L n ( 1 F ( x ) N + 1 i ) )
where A2 is the test statistic, N is the sample size, and F(x) is the frequency.

2.3. Case Study

The present study was performed using annual rainfall data from 52 meteorological stations in Badajoz, provided by the Spanish Meteorological Institute.
An exhaustive study was made using the statistical goodness-of-fit tests of Anderson–Darling and Kolmogorov–Smirnov to establish the degree of fitness of the Dagum distribution applied to the maximum annual rainfall series, and thus, be able to compare the fits of this distribution to those of the classical statistical distributions such as Gumbel, SQRT-ET max, Log-Pearson type III and the GEV.
Adjustment tests are widely used in hydrology due to the high degree of precision they provide when reflecting the fit of the statistical distribution to the rainfall data series available.
For the fit of the distributions, the software Easyfit was used. On introducing the maximum annual rainfall data, the program gives the degree of fitness of each statistical distribution according to the two goodness-of-fit tests used.
The 52 meteorological stations used in the analysis are shown in Table 1 and Figure 2. These contained maximum and minimum temperatures and daily precipitation for the period between 1990 and 2015. The quality control procedures of the Algorithm Theoretical Basis Document (ATBD) project, developed by the Royal Netherlands Meteorological Institute (KNMI) for the European Climate Assessment & Dataset (ECA&D), have been applied [43]. The blended series passed the standard homogeneity test, the Buishand range test, the Pettitt test and the Von Neumann ratio, as described by Wijngaard et al. [44] and ECA&D. Some series presenting missing values were completed following the recommendations of WMO [45] and Allen et al. [46]. The daily data from each station were processed and analyzed. Coefficients of variation and maximum precipitations at each meteorological station are also shown in Table 1.
The rainfall data from the 52 meteorological stations were introduced in EasyFit statistical program. The function of density was applied to the rainfall histogram for its adjustment. Finally, the goodness-of-fit to the rainfall histograms was studied as a function of the density of each statistical distribution.

3. Results

Figure 3 shows how the Dagum probability distribution fits the rainfall histogram in the town of San Vicente de Alcántara.
Figure 4 shows how the density function fits the cumulative histogram of a maximum annual rainfall series. The density function of the statistical distribution is never going to reproduce the exact values of the histogram as desired that is why the goodness-to-fit tests are used to check which one provides the best fit to the rainfall data series when comparing various statistical distributions. That is to say, the chosen statistical distribution should be the one that fits the rainfall histogram more accurately.
After fitting the statistical functions to the rainfall data, Figure 5 shows graphically how the Dagum distribution is aligned with both the fit of Gumbel distribution and Log-Pearson type III distribution, being difficult to decide the best, since the curves are very close together.
Similarly, Figure 6 shows that it is difficult to determine which of the three distributions from all the cumulative distribution function fits more accurately the histogram of rainfall data in San Vicente de Alcántara.
As previously noted, it is difficult to draw any conclusions. Therefore, it is necessary to apply the goodness-to-fit tests. In this study, we applied the test to four distributions, as shown in Table 2. The distribution of Dagum clearly presents the lowest goodness-to-fit statistics in the San Vicente de Alcántara, Jerez de los Caballeros and Herrera del Duque data, which means that it fits better than the other three to the rainfall data.
These tests were applied to rainfall data from the remaining 51 stations [47], in which the analysis reflects a similar trend, and the statistical distribution Dagum presents lower goodness statistics than the rest of the distributions. Other tests were taken into account to select the best distribution model, such as the corrected Akaike information criterion and the Bayesian information criterion (BIC), but finally the tests of Anderson–Darling and the test of Kolmogorov–Smirnov were chosen.
Subsequently, and to confirm the above results, tests of goodness-of-fit were carried out with a series of statistical distributes (applied both in hydrology and in other disciplines), using the ten stations with the largest sample size (among the 52 stations).
In Figure 7, it can see that Dagum appears as one of the most frequent distributions (within the five best settings), just below the GEV distribution, but above the Gumbel and the Log-Pearson type III distribution.
Table 3 shows that the application of the Dagum Distribution to the cases of Cabeza la Vaca, Monterrubio and Campanario obtain different flows to those of the traditional distributions. For example, Cabeza la Vaca shows 21.92% flow greater than the most commonly used distribution (Gumbel).
It can be deduced that the statistical distributions that provide the greatest rainfall are the Log- Logistic 3P distribution followed by the Dagum distribution. The quantitative differences in relation to the value provided by the Gumbel distribution, the most widespread in studies and projects, are variable and in some cases are even considerable, up to 58% higher in the case of Jerez de los Caballeros. Therefore, in order to be on the safety side, the Log-Logistic 3P distributions and the Dagum distribution, must be used, since you can be completely sure that the flows and precipitations derived from its application will be greater than the results obtained with the distributions of Gumbel, SQRT-ET max and Log Pearson 3. It is important to emphasize that these distributions are the best fit in the Kolmogorov–Smirnov, Anderson–Darling kindness tests and their weighting.

4. Discussion

It is important to denote that the comparison of distributions is complicated. Goodness of fit tests are not very powerful and with the typical sample sizes available in practice it is rarely possible to reject statistically some distribution candidates. Therefore, the comparison must be done on a larger scale [48]. Thus, in this case study, sites located throughout the province considered are sufficiently numerous and evenly distributed to obtain significant results. However, these results cannot be extended elsewhere, that is, the choice of a particular distribution at a given place should be carefully studied and selected. Although some distributions, such as the Gumbel or Log-Pearson type III, have been extensively used in many hydrologic studies and without any additional consideration related to particular conditions of the basins, the inclusion of the spatial factor would reduce the uncertainty concerning the choice of the model [49]. Physical factors such as large-scale meteorological phenomena could create regional probabilities dependencies which have to be accounted for. In consequence, as it was previously indicated, each region or zone should be initially characterized for the choice of the statistical distribution, which better explains the expected rainfall events [50].
As a result of the analyses carried out in the 52 locations throughout the province of Badajoz, the Dagum model scored better than the other models which have been traditionally used in hydrologic and hydraulic works. There are few previous studies where the Dagum distribution had been used for these topics.
The Dagum model was found to overestimate a great number of times when compared to the Gumbel distribution. Therefore, the Dagum distribution seems to be the most recommendable distribution for a conservative design and for to plan accordingly [51].
Because of the ample availability of computers nowadays, many statistical distributions have to be considered when a single-site flood frequency analysis is done [52]. Moreover, as more data are being accumulated since the recent and coming years are providing more new information, new analyses could be performed with regionalized parameters of proven model for each location [53]. In this sense, the consideration of the Dagum model can provide more accurate results in many places of southwestern Spain that those obtained using traditional distributions. The evaluation and simulation of rainfall scenarios indicate that changes in rainfall characteristics have a considerable impact on the built drainage system and that Low Impact Development (LID) practices can adequately control flooding [54].
Future work should aim at verifying the applicability of the Dagum distribution in other regions of southern Europe.

5. Conclusions

On analyzing the maximum annual rainfall data from 52 stations (strategically located throughout the zone) and treating them by using the goodness-of-fit tests of Anderson–Darling and Kolmogorov–Smirnov, it is confirmed that, in addition to the distributions traditionally used in hydrology (such as the Gumbel, Log-Pearson type III and the EVD distribution), there is another statistical distribution, the Dagum, which can be used in hydrology and meets the formulation of extreme values (outliers) and fits better to the rainfall histograms.
Based on the statistical data from the study, it is concluded that the Dagum distribution presents lower statistics in the two goodness-of-fit tests mentioned above and, therefore, adjusts significantly better to the histograms of the maximum annual rainfall data than the commonly used distributions. Particularly, this new statistical distribution is more appropriate to reflect the rainfall regime in Badajoz.
In conclusion, the Dagum statistical distribution is proposed to improve hydrological studies in Badajoz, since the rainfalls given by its density function are more precise (as shown by the goodness-of-fit tests) than the rainfall data calculated through classical statistical distributions. Its use in the professional field would allow for greater flows rates to be considered when designing drainage systems and studies of flooding, thus preventing future possible rainfall damage.

Author Contributions

Conceptualization, F.L.-R.; Data curation, A.G.-C.; Formal analysis, F.M.-G.; Funding acquisition, F.M.-G.; Investigation, F.L.-R., J.G.-S.-C. and A.G.-C.; Methodology, F.L.-R.; Project administration, J.G.-S.-C.; Resources, F.M.-G.; Software, A.G.-C.; Supervision, J.G.-S.-C.; Validation, J.G.-S.-C.

Funding

This research received no external funding.

Acknowledgments

The authors wish to acknowledge to the Junta de Extremadura and the European Social Found (FEDER) for the support of this research work. This study has been carried out through the Research Project GR-18029, GR-18081 and GR-18086 linked to the VI Regional Plan for Research, Technological Development and Innovation from the General Government of Extremadura 2017–2020.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Sanches Fernandes, L.; Pereira, M.-G.; Morgado, S.G.; Macário, E.B. Influence of Climate Change on the Design of Retention Basins in Northeastern Portugal. Water 2018, 10, 743. [Google Scholar] [CrossRef]
  2. Bai, T.; Mayer, A.L.; Shuster, W.D.; Tian, G. The Hydrologic Role of Urban Green Space in Mitigating Flooding (Luohe, China). Sustainability 2018, 10, 3584. [Google Scholar] [CrossRef]
  3. Gumbel, E.J. Multivariate Extremal Distributions; Bull. Inst. International the Statistique: Paris, France, 1960. [Google Scholar]
  4. Chow, V.T.D.R.; Maidment, L.W.M. Hidrología Aplicada; McGraw-Hill: Colombia, Santa Fé de Bogotá, 1994. [Google Scholar]
  5. Ponce, V.M. Engineering Hydrology, Principles and Practices; Prentice Hall: Upper Saddle River, NJ, USA, 1989. [Google Scholar]
  6. Pearson, K. Mathematical Contributions to the Theory of Evolution, XIX: Second Supplement to a Memoir on Skew Variation; Philosophical Transactions of the Royal Society: London, UK, 1916. [Google Scholar]
  7. Ferrer, F.J. Obtención de la lluvia neta según la metodología del Soil Conservation Service; Cedex: Madrid, Spain, 1991. (In Spanish) [Google Scholar]
  8. Témez, J.R. Extended and Improved Rational Method. Version of the Highways Administration of Spain; XXIV Congreso Internacional de la IAHR: Madrid, Spain, 1991. [Google Scholar]
  9. Kotz, S.; Nadarajah, S. Extreme Value Distributions: Theory and Applications; Imperial College Press: London, UK, 2000. [Google Scholar]
  10. Zorraquino, C. La Función SQRT-ET Max; Revista de Obras Públicas: Madrid, Spain, 2004. (In Spanish) [Google Scholar]
  11. de Fomento, M. Máximas Lluvias diarias en la España Peninsular; Ministerio de Fomento: Madrid, Spain, 1999. (In Spanish)
  12. Ministerio de Obras Públicas y Urbanismo. Instrucción de Carreteras, Norma 5.2-IC de Drenaje Superficial. It has been doneMadrid; It has been doneMadrid; Ministerio de Obras Públicas y Urbanismo: Madrid, Spain, 1990. (In Spanish)
  13. Catalá, F. Cálculo de caudales en las redes de saneamiento; Colegio de Ingenieros de Caminos, Canales y Puertos: Madrid, Spain, 1997. [Google Scholar]
  14. Ferrer, F.J. Recomendaciones para el cálculo hidrometeorológico de avenidas; Cedex: Madrid, Spain, 1993. (In Spanish) [Google Scholar]
  15. Panagoulia, D.; Caroni, C. Modelling maximum precipitation in a mountainous area of Greece under global warming, Hydro-climatology: Variability and Change, Proceedings of Symposium J-H02 held during IUGG2011 in Melbourne, Australia. IAHS Publ. 2011, 344, 176–181. [Google Scholar]
  16. Lang, M.; Ouarda, T.B.M.J.; Bobée, B. Towards operational guidelines for over-threshold modeling. J. Hydrol. 1999, 225, 103–117. [Google Scholar] [CrossRef]
  17. Onyutha, C. On rigorous drought assessment using daily time scale: Non-stationary frequency analyses, revisited concepts, and a new method to yield non-parametric indices. Hydrology 2017, 4, 48. [Google Scholar] [CrossRef]
  18. Onyutha, C. Hydrological model supported by a step-Wise calibration against sub-flows and validation of extreme flow events. Water 2019, 11, 244. [Google Scholar] [CrossRef]
  19. Pickands, J. Statistical inference using extreme order statistics. Ann. Stat. 1975, 3, 119–131. [Google Scholar]
  20. Modarres, R.; Sarhadi, A. Statistically-based regionalization of rainfall climates of Iran. Glob. Planet. Change 2011, 75, 67–75. [Google Scholar] [CrossRef]
  21. Dagum, C. Un nuevo modelo de distribución del ingreso personal: Especificación y estimación. Economía Aplicada 1977, 30, 413–437. (In Spanish) [Google Scholar]
  22. Kleiber, C. A guide to the Dagum distributions. In Modeling Income Distributions and Lorenz Curves; Springer: New York, NY, USA, 2008; pp. 97–117. [Google Scholar]
  23. Mielke, P.R., Jr.; Johnson, E.S. Some generalized beta distributions of the second kind having desirable application features in hydrology and meteorology. Water Resour. Res. 1974, 10, 223–226. [Google Scholar] [CrossRef]
  24. Alam, S.; Sabbir, M.; Khan, M. Statistical Characterization of Extreme Hydrologic Parameters for the Peripheral River System of Dhaka City. J. Water Resour. Ocean Sci. 2014, 3, 30–37. [Google Scholar] [CrossRef]
  25. Alam, S.; Khan, S.M.; Rahat, S.H. A Study on Selection of Probability Distributions of Extreme Hydrologic Parameters for the Peripheral River System of Dhaka City. In Proceedings of the Academics World 15th International Conference, Bangkok, Thailand, 28–29 December 2015; pp. 29–34. [Google Scholar]
  26. Crabbe, P.J. Conceptual Framework for a Systems Dynamics Adaptation Model to Climate Change for Charlottetown, P.E.I., Canada. J. Glob. Resour. 2016, 2, 1–11. [Google Scholar]
  27. Mayooran, T.; Laheetharan, A. The Statistical Distribution of Annual Maximum Rainfall in Colombo District. Sri Lankan J. Appl. Stat. 2014, 15, 107–130. [Google Scholar] [CrossRef]
  28. Pérez-Sánchez, J.; Senent-Aparicio, J. Analysis of meteorological droughts and dry spells in semiarid regions: A comparative analysis of probability distribution functions in the Segura Basin (SE Spain). Theor. Appl. Climatol. 2018, 133, 1061. [Google Scholar] [CrossRef]
  29. Domma, F.; Condino, F. The Beta-Dagum distribution: Definition and properties. Commun. Stat. Theory Methods 2013, 42, 4070–4090. [Google Scholar] [CrossRef]
  30. Cheng, Y.; Wang, R. A novel stormwater management system for urban roads in China based on local conditions. Sustain. Cities Soc. 2018, 39, 163–171. [Google Scholar] [CrossRef]
  31. Douka, M.; Karacostas, T. Statistical analyses of extreme rainfall events in Thessaloniki, Greece. Atmos. Res. 2018, 208, 60–77. [Google Scholar] [CrossRef]
  32. Shahzad, M.N.; Asghar, Z. Transmuted dagum distribution: A more flexible and broad shaped hazard function model. Hacet. J. Math. Stat. 2016, 45, 227–244. [Google Scholar] [CrossRef]
  33. Panagoulia, D.; Economou, P.; Caroni, C. Stationary and non-stationary generalized extreme value modeling of extreme precipitation over a mountainous area under climate change. Environmetrics 2014, 25, 29–43. [Google Scholar] [CrossRef]
  34. Caroni, C.; Panagoulia, D.; Economou, P. Non–stationary modelling of extremes of precipitation and temperature over mountainous areas under climate change. In Proceedings of the International Conference in Current Topics on Risk Analysis: ICRA6 and RISK, Barcelona, Spain, 26–29 May 2015; pp. 203–209. [Google Scholar]
  35. Anderson, T.W.; Darling, D.A. Asymptotic theory of certain “goodness-of-fit” criteria based on stochastic processes. Ann. Math. Stat. 1952, 23, 193–212. [Google Scholar] [CrossRef]
  36. Gumbel, E.J. Les valeurs extrêmes des distributions statistiques. Ann. Inst. Henri Poincaré 1935, 5, 115–158. [Google Scholar]
  37. Champernowne, D.G. The graduation of income distributions. Econ. J. Econom. Soc. 1952, 20, 591–615. [Google Scholar] [CrossRef]
  38. MATHWAVE. Available online: http://www.mathwave.com/easyfit-distribution-fitting.html (accessed on 15 December 2018).
  39. Kolmogorov, A.N. Foundations of the Theory of Probability; Chelsea Publishing Company: New York, NY, USA, 1950. [Google Scholar]
  40. Miller, I.; Freund, J.E. Probabilidad y estadística para ingenieros; Reverté: Madrid, Spain, 1969. [Google Scholar]
  41. Shorack, G.R.; Wellner, J.A. Empirical Processes with Applications to Statistics; Wiley: New York, NY, USA, 1986. [Google Scholar]
  42. Anderson, T.W.; Darling, D.A. A Test of Goodness-of-Fit. J. Am. Stat. Assoc. 1954, 49, 765–769. [Google Scholar] [CrossRef]
  43. ECA&D. Algorithm Basis Document (ATBD); Royal Netherlands Meteorological Institute KNMI: De Bilt, The Netherlands, 2013. [Google Scholar]
  44. Wijngaard, J.B.; Tank, A.K.; Können, G.P. Homogeneity of 20th century European daily temperature and precipitation series. Int. J. Climatol. 2003, 23, 679–692. [Google Scholar] [CrossRef]
  45. WMO. Guide to Climatological Practices (WMO No. 100), 3rd ed.; WMO: Geneva, Switzerland, 2011. [Google Scholar]
  46. Allen, R.G.; Pereira, L.S.; Raes, D.; Smith, M. Crop Evapotranspiration. Guidelines for Computing Crop Water Requirements; FAO Irrigation and Drainage Paper, Nº 56; FAO: Rome, Italy, 1998. [Google Scholar]
  47. García-Conde, A.J. Análisis de distribuciones estadísticas alternativas a las tradicionales para la optimización de los caudales de cálculo empleados en los estudios de hidrología. Ph.D. Thesis, University of Extremadura, Badajoz, Spain, 2013. (In Spanish). [Google Scholar]
  48. Bobée, B.; Cavadis, G.; Ashkar, F.; Bernier, J.; Rasmussen, P. Towards a systematic approach to comparing distributions used in flood frequency analysis. J. Hidrol. 1993, 142, 121–136. [Google Scholar] [CrossRef]
  49. Haktanir, T.; Horlacher, H.B. Evaluation of various distributions for flood frequency analysis. Hidrol. Sci. J. Sci. Hydrol. 1993, 38, 15–32. [Google Scholar] [CrossRef] [Green Version]
  50. Courty, L.G.; Rico-Ramirez, M.A.; Pedrozo-Acuña, A. The Significance of the Spatial Variability of Rainfall on the Numerical Simulation of Urban Floods. Water 2018, 10, 207. [Google Scholar] [CrossRef]
  51. Candelario-Garrido, A.; García-Sanz-Calcedo, J.; Reyes Rodríguez, A.M. A quantitative analysis on the feasibility of 4D planning graphic systems versus conventional systems in building projects. Sustain. Cities Soc. 2017, 35, 378–384. [Google Scholar] [CrossRef]
  52. Panagoulia, D.; Economou, P.; Caroni, C. Quantifying uncertainty in modelled estimates of annual maximum precipitation: Confidence intervals, European Geosciences Union. Precipitation uncertainty and variability: Observations, ensemble simulation and downscaling. Geophys. Res. Abstr. 2016, 18, 17–22. [Google Scholar]
  53. Carretero-Ayuso, M.J.; García-Sanz-Calcedo, J. Analytical study on design deficiencies in the envelope projects of healthcare buildings in Spain. Sustain. Cities Soc. 2018, 42, 139–147. [Google Scholar] [CrossRef]
  54. Zhihua, Z.; Xiaohong, C. Evaluating the Effects of Low Impact Development Practices on Urban Flooding under Different Rainfall Intensities. Water 2017, 9, 548. [Google Scholar] [CrossRef]
Figure 1. Density function of the Dagum distribution as a function of its parameters.
Figure 1. Density function of the Dagum distribution as a function of its parameters.
Water 11 00453 g001
Figure 2. Map of the spatial distribution of the meteorological stations. “Coordinate system is ETRS89/UTM zone 30” N.
Figure 2. Map of the spatial distribution of the meteorological stations. “Coordinate system is ETRS89/UTM zone 30” N.
Water 11 00453 g002
Figure 3. Dagum probability density distribution function in San Vicente de Alcántara for γ = 0.
Figure 3. Dagum probability density distribution function in San Vicente de Alcántara for γ = 0.
Water 11 00453 g003
Figure 4. Cumulative distribution function for the Dagum distribution applicable to data from San Vicente de Alcántara.
Figure 4. Cumulative distribution function for the Dagum distribution applicable to data from San Vicente de Alcántara.
Water 11 00453 g004
Figure 5. Probability density function of Gumbel, Log-Pearson type III and Dagum Distributions. Comparison of the fit among the three different distributions applied to the histogram of the data from San Vicente de Alcántara.
Figure 5. Probability density function of Gumbel, Log-Pearson type III and Dagum Distributions. Comparison of the fit among the three different distributions applied to the histogram of the data from San Vicente de Alcántara.
Water 11 00453 g005
Figure 6. Cumulative distribution functions of Gumbel, Log-Pearson type III and Dagum distributions, applied to the histogram of data from San Vicente de Alcántara.
Figure 6. Cumulative distribution functions of Gumbel, Log-Pearson type III and Dagum distributions, applied to the histogram of data from San Vicente de Alcántara.
Water 11 00453 g006
Figure 7. Number of events in the different distributions according to the goodness-of-fit test at the locations for the ten largest simple sizes.
Figure 7. Number of events in the different distributions according to the goodness-of-fit test at the locations for the ten largest simple sizes.
Water 11 00453 g007
Table 1. Meteorological stations used in the analysis.
Table 1. Meteorological stations used in the analysis.
Meteorological StationsCV (%)Pmax (mm)Meteorological StationsCV (%)Pmax (mm)
1San Vicente Alcántara19.6112.627GranjaTorrehermosa25.9116.5
2Alburquerque22.3104.228Maguilla27.990.3
3La Roca de la Sierra27.8140.629Peraleda del Zaucejo30.275.3
4Badajoz (Sagrajas)29.5102.230Puebla del Prior27.6181.8
5Montijo (Institute)32.5135.831Puebla de la Reina22.3173.1
6Olivenza30.5126.332Monterrubio19.698.7
7Cheles19.779.133Puerto Hurraco22.6106.2
8Alconchel25.9121.334Castuera27.495.4
9Villanueva del Fresno22.2116.235Quintana de Serena25.4102.3
10Valencia del Mombuey19.9188.636Valle de la Serena29.689.7
11Zahinos30.1163.337Alange22.8119.5
12Higuera de Vargas25.6217.438Manchita23.6109.3
13Jerez de los Caballeros26.8171.439Guareña28.4116.5
14Barcarrota 27.5157.140Aceuchal26.892.5
15La Albuera31.2146.641Mérida30.9107.8
16Talavera la Real29.8116.542Santa Amalia26.1143.1
17Fregenal de la Sierra19.8108.243La Coronada20.9135.4
18Segura de León18.5182.144Campanario19.999.1
19Cabeza la Vaca16.9164.945Acedera22.676.8
20Fuente de Cantos19.6124.546Orellana la Sierra28.5148.3
21Puebla del Maestre28.6134.247Casas de Don Pedro30.292.2
22Casas de Reina25.8149.848Capilla/Baterno29.4122.3
23Villagarcía de la Torre27.6162.849Siruela25.4145.5
24Berlanga28.9134.950Herrera del Duque27.2130.3
25Valverde de Lerena30.2141.651Villarta de Montes22.1140.1
26Azuaga27.4150.652Helechosa25.698.5
Table 2. Goodness-of-fit test results.
Table 2. Goodness-of-fit test results.
DistributionSan Vicente AlcántaraJerez de los CaballerosHerrera del Duque
Tests of Anderson–DarlingTest of Kolmogorv–SmirnovTests of Anderson–DarlingTest of Kolmogorov–SmirnovTests of Anderson–DarlingTest of Kolmogorov–Smirnov
Dagum0.04900.17460.04210.09140.03910.1257
Log-Pearson III0.05840.19470.06240.17580.59580.1985
Gumbel0.06030.25530.05870.14040.06870.0984
GEV0.05640.18750.05900.06540.05930.1751
Table 3. Maximum daily precipitation for different distributions.
Table 3. Maximum daily precipitation for different distributions.
DistributionCabeza la VacaMonterrubioCampanario
P24T(mm)Variation Respect GumbelP24T(mm)Variation Respect GumbelP24T(mm)Variation Respect Gumbel
GEV158.470.43%98.23−12.96%89.03−0.15%
Dagum192.3721.92%123.049.28%112.7726.47%
Log-Logistic 3P188.9719.76%124.049.90%116.1730.29%
Fechet 3P151.64−3.90%110.94−1.69%85.35−4.28%
Pearson 5 3P146.19−7.35%100.3811.05%86.47−3.02%
Gumbel157.79-112.85-89.16-
SQRT-ET max164.904.51%118.805.27%93.454.81%
Log Pearson Type III145.01−8.10%97.79−13.34%84.97−4.71%
Regional analysis172.009.00%118.004.56%98.009.91%

Share and Cite

MDPI and ACS Style

López-Rodríguez, F.; García-Sanz-Calcedo, J.; Moral-García, F.J.; García-Conde, A.J. Statistical Study of Rainfall Control: The Dagum Distribution and Applicability to the Southwest of Spain. Water 2019, 11, 453. https://doi.org/10.3390/w11030453

AMA Style

López-Rodríguez F, García-Sanz-Calcedo J, Moral-García FJ, García-Conde AJ. Statistical Study of Rainfall Control: The Dagum Distribution and Applicability to the Southwest of Spain. Water. 2019; 11(3):453. https://doi.org/10.3390/w11030453

Chicago/Turabian Style

López-Rodríguez, Fernando, Justo García-Sanz-Calcedo, Francisco J. Moral-García, and Antonio J. García-Conde. 2019. "Statistical Study of Rainfall Control: The Dagum Distribution and Applicability to the Southwest of Spain" Water 11, no. 3: 453. https://doi.org/10.3390/w11030453

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop