The Moroccan territory is often exposed to brutal and destructive floods. The latter is triggered by intense rains and sudden duration (rain showers, thunderstorms). The human and material damage caused by these floods has always been very high. The violence of the flood (1000 m 3/s) of Oued Ourika on Aug 17, 1995, killed 730 people. The material damage caused by torrential rain only for the years 2008 and 2009 on the road network (roads and civil engineering works) represents more than 2.9 billion dirhams. The development of flood protection studies requires the values of the instantaneous frequency flow of a certain probability and the corresponding flood hydrogram, generally, these values are not available in the watersheds studied. The objective of this study is to propose for the Moroccan territory a method of calculating flood flows and hydrographs of floods for stream and wadis ungauged for the benefit of engineers and planners for possible studies of flood protection. The proposed method is based on the delimitation of homogeneous regions of Moroccan territory (three homogeneous regions). The frequency analysis carried out for these three zones made it possible to identify the climatic parameters necessary to calculate the flood index and consequently the frequency floods. The proposed Galton type flood hydrograph includes three parameters namely the rise time, the shape parameter and the flood index. These latter parameters are easily calculated by the formulas proposed.
The Moroccan territory is often exposed to brutal and destructive floods. Often well-defined meteorological situations lead to violent and torrential floods which are generally linked to intense, stormy and localized precipitation. Torrential mountain floods are characterized by their destructive force, they start after short bursts of showers and thunderstorms.
In Ourika, a strong thunderstorm which erupts upstream on the mountain (more than 100 mm in less than an hour) on Aug 17, 1995, led to a very rapid increase in flows which would have reached approximately 1000 m3/s. The violence of this flood killed 730 people and the material damage was very high [
Following the exceptional precipitation from 20 to 30 November 2014 compared to the usual standards in the Tensift basin (136 mm in Yaacoub El Mansour, 360 mm in Armed, 291 mm in Tazitoune and 519 mm in Tourcht), all the rivers in the basin carried exceptional flows. The most spectacular flood flows of November 2014, recorded on Oued Tensift reached 1597 m3/s and 3500 m3/s at the stations of Abadla and Talmest respectively. Such floods have never been observed before in this region [
Between the years (1999-2009), Morocco experienced fifteen major floods, which left 1068 dead and affected more than 146,400 people. The damage caused by these floods was very considerable [
If a long series of observation rates are available, statistical analysis allows this assessment to be made. The construction of a flood hydrogram is based on the transformation of the hyetogram of net rain into a hydrogram of direct flow. In the majority of cases, the hydro-meteorological data necessary for the statistical analysis of flood flows and the construction of their hydrographs are not available for the watershed studied. In this case, the evaluation of these characteristics is carried out using hydrological models, most often, empirical models, established on the basis of other available information (physiographic, hydrological and rainfall data) in the region or in the vicinity of the site studied. However, extrapolating these models to other regions with hydrometeorological contexts different from the one for which they were developed can lead to outliers.
The objective of this study is to propose for the Moroccan territory a method of calculating the flows of floods and hydrographs of the floods for the stream and the wadis ungauged for the benefit of engineers and planners, exercising in the hydraulics. The proposed method is based on frequency analysis and the flood index.
The regional frequency analysis method is applied for homogeneous regions [
The regional frequency analysis of floods allows the estimation of frequency floods from the sole knowledge of the annual maximum flood (flood index). This index can be reasonably estimated from only a few years of measurement [
Concerning the standardized frequency distribution curve, it is generally modeled for each of the regions by means of quantiles calculated by the laws of statistical distribution [
The data on the flood flows and the hydrographs of the watersheds used in this study come from the database of the Agencies of the Hydraulic Basins (ABH) of Morocco. Part of the data is collected from scientific publications and state doctoral theses. The latter is used just to fill the geographic gaps by verifying each time by the method of hydrological analogy the homogeneity of these data.
The number of hydrometric stations is 50 (
In the absence of a more elaborate approach, unknown to date, allowing the estimation of flood flows, the method of regional frequency analysis is used in several countries of the world for the regionalization of flows.
Delimitation of homogeneous regions on Moroccan territory
If a group of rivers or wadis from the same region exhibits similar hydrological behavior, the region is considered to be homogeneous. In this case, it can be studied separately from the other zones [
| Area of the watershed S (km2) | ||||
|---|---|---|---|---|
| <100 | 100 - 1000 | 1000 - 3000 | 3000 - 6000 | >6000 |
| 3 | 18 | 11 | 10 | 8 |
There are several methods for delimiting homogeneous regions, the most used are based on geographic criteria, physiographic and/or climatic criteria and hydrological criteria [
The Analysis of the relationship between the specific flow (q0, m3/s∙km2) and the area of the watershed (S, km2) in the different geographic regions (Central Rif, Western Rif, Eastern Rif, Middle Atlas and Saharan High Atlas) and climatic (humid zone, subhumid zone and semi-arid zone) shows a certain distribution of points around three straight lines (
The three lines expressing the relation q0 = fonct(S) obtained make it possible to delimit three homogeneous regions for the Moroccan territory according to the annual rainfall.
Regional frequency analysis
The basic assumption of regional frequency analysis is that the region must be homogeneous. In this case, all the frequency distributions of the maximum annual floods at each of the N stations i are identical, to a scale factor which is the flood index. It is therefore possible to write [
Q i = μ i K ( T ) , i = 1 , ⋯ , N (1)
where:
μ i is a scale factor representing the average of the maximum annual floods, it is also called the flood index;
K ( T ) is standardized regional frequency distribution.
The frequency analysis of the observation data of flood flows for each homogeneous region was developed on the basis of asymmetric hypergaussian laws (lognormal law, Gumbel law, Pearson III law, Gamma law, etc.). With the exception of a few series, all of these laws can present flood flows in Moroccan territory.
In very irregular regimes such as those observed in the Mediterranean region and in arid and semi-arid climates, several authors recommend calculating flood flows with low probability of occurrence by the log-normal law [
This law was used to calculate the quantiles of the flow rates of rivers in Moroccan territory. The normalization of flood flows was carried out using the standardized frequency distribution:
K ( T ) = Q ( T ) Q 0 (2)
where:
K ( T ) is standardized frequency distribution;
Q ( T ) is return period flood flow T;
Q 0 is flood index.
The mathematical expression K ( T ) = f o n c t ( T ) for the three homogeneous regions is the same, it has the following form:
K ( T ) = ξ exp [ λ ⋅ log ( T ) ] (3)
For the three homogeneous regions, the value of the parameter ξ is close to unity ( ξ ≅ 1 ) . The value of the parameter λ is a function of the annual rainfall.
Determination of the flood flow index
For the Moroccan territory, R. Hazen and D. Lazarvic [
Q 0 = A ( S ) m (4)
The formula (4) provides for two watersheds of the same size, located in the same geographic region but with different annual rainfall, the same flood index Q0. In reality, Q0 is different for these two basins (
The formula of R. Hazen and D. Lazarvic for calculating the flood index Q0 therefore requires corrections.
The determination of the unit hydrogram can be done directly if rain-flood pairs are known, either otherwise by empirical methods or by comparison with similar basins.
The data from the showers and their time distributions making it possible to define the unit hydrogram are not available for the case of the watersheds studied. To allow the selection of the unit hydrographs observed, we based on the following hypotheses:
| River (Station) | Area of the watershed S, km2 | Average annual rainfall P, mm | Flood index Q0, m3/s Measured |
|---|---|---|---|
| R’Dat (Sidi Rahal) | 532 | 460 | 185 |
| Zat (Tafériet) | 528 | 521 | 125 |
· Small watersheds of small sizes are subject to showers of constant intensity and uniformly distributed in space (unitary downpour) and sufficient to generate a unit hydrograph. For arid and semi-arid climates, the uniformity of the downpour can be perfectly obtained on these watersheds.
· Hydrographs distorted compared to the unit hydrogram generally come from localized showers. Flood hydrographs of this type will not be selected.
· Complex flood hydrographs are generally provided by generalized floods, if their flows are large it is possible to dissociate them into several unit hydrographs and therefore increase the information.
· Large watersheds are rarely fully covered by downpours. In this case, we only select flood hydrographs with low probability of occurrence (very high peak flows) which are more likely to be formed by generalized runoff over the entire watershed.
It is obvious that the superimposition of the reduced unit hydrographs of a catchment area, thus obtained for all floods considered unitary provides, by consideration of the median or average ordinates, a typical hydrograph which can be considered as a hydraulic characteristic of this basin. The relative values of the ordinates (in % of the total volume of runoff) of the hydrogram distribution diagram resulting from some unitary rains at the Tahanaout station (Oued Rheraya) are given in
From the selected unit flood flows, a typical hydrograph for each watershed was calculated based on the averages of the established standard hydrographs. This mode of representation is more direct because it avoids going through the volume. The typical hydrographs of certain rivers in Moroccan territory are illustrated in
| Date of flood discharge | Hours | ||||||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|
| 0 | 1 | 2 | 3 | 4 | 5 | 6 | 7 | 8 | 9 | 10 | 11 | 12 | 13 | 14 | |
| Apr. 1, 2002 | 0.0 | 0.0 | 0.0 | 5.3 | 18.6 | 20.7 | 19.1 | 13.6 | 7.4 | 5.2 | 4.1 | 2.7 | 2.1 | 1.0 | 0.0 |
| Aug. 14, 2003 | 0.0 | 0.0 | 0.0 | 0.0 | 17.8 | 37.9 | 23.2 | 16.1 | 2.7 | 0.9 | 0.7 | 0.7 | 0.0 | 0.0 | 0.0 |
| Aug. 17, 1995 | 0.0 | 0.0 | 0.0 | 0.0 | 0.0 | 31.8 | 27.2 | 22.9 | 10.9 | 5.4 | 0.8 | 0.5 | 0.4 | 0.0 | 0.0 |
| Feb. 10, 1987 | 0.0 | 2.2 | 8.0 | 10.9 | 19.6 | 21.8 | 11.8 | 8.7 | 6.0 | 6.0 | 3.1 | 2.0 | 0.0 | 0.0 | 0.0 |
| Average | 0.0 | 0.5 | 2.0 | 4.1 | 14.0 | 28.1 | 20.3 | 15.3 | 6.7 | 4.4 | 2.2 | 1.5 | 0.6 | 0.3 | 0.0 |
| Median | 0.0 | 0.3 | 1.0 | 2.0 | 15.9 | 29.9 | 21.8 | 15.7 | 6.3 | 4.9 | 1.5 | 1.1 | 0.2 | 0.0 | 0.0 |
Analysis of the standard hydrograph established for each watercourse shows that the rising part of the curve is tight on the left that of the descending part and spread out to the right.
We suppose that the analytical form adopted for these standard hydrographs is that of the density function of the log-normal law (Galton law) which is chosen on the basis of the following considerations:
1) law in the form of an asymmetrical bell, spread to the right.
2) law taking only positive values.
3) limited number of parameters constituting the law.
The density of the Galton law is presented in the following form [
f ( x ) = { 0 p o u r x ≤ 0 1 x 2 π σ exp [ − 1 2 ( l n x σ ) 2 ] p o u r x > 0 (5)
The analysis of this mathematical formulation by setting f ( x ) = Q T and x = t / t m shows that the closest distribution of a standard adimensional unit hydrograph or of a standard hydrogram whose volume is equal to the volume arbitrary fixed and whose shape represents on average that of all the hydrographs of runoff, is obtained almost for σ = 0.4 (
The analysis of the relationship between the average specific flow of the flood q0 and the area of the watershed S in the different geographical and climatic regions of the entire Moroccan territory has made it possible to delimit three homogeneous regions of genesis of the rain floods. This relationship also made it possible to define two laws for the Moroccan territory which can be expressed in the following manner:
· For the same area of the watershed, the specific flow of the annual flood q0 increases if the climate changes from a semi-arid region to a humid region.
· For the same climatic zone, the specific flow of the average annual flood q0 decreases with the growth of the area of the watershed.
These two laws confirm once again that the specific flood flow q0 and therefore the flood index (maximum annual average flood flow Q0 is a function not
only of the surface of the basin but also of the pluviometry. The relationship proposed for the Moroccan territory as an alternative to that of R. Hazen and D. Lazarvic links the flood index Q0, the area of the watershed S and the average annual rainfall P received by the latter, it is as follows:
Q 0 = β P ( S ) n (6)
where:
Q0 = annual average maximum flood flow (m3/s);
S = area of the watershed (km2);
P = average annual rainfall received by the watershed (mm);
β and n are parameters, a function of the climatic zone and the average annual rainfall (
The values of the flood index calculated by the formula (6) do not differ much from the values of the observed flood index. The correlation coefficient between the calculated and observed values is greater than 0.8 (
The standardized regional frequency distribution model for the three homogeneous regions of Moroccan territory is as follows:
K ( T ) = exp [ λ ⋅ log ( T ) ] (7)
The parameter λ is tabulated according to the climatic zone and the average annual rainfall (
| N˚ | Climatic zone | Average annual rainfall P, mm | Parameters | ||
|---|---|---|---|---|---|
| β | n | λ | |||
| I | humid | ≥900 mm | 0.0025 | 0.80 | 0.6 |
| II | Subhumid | 900 - 400 mm | 0.0103 | 0.50 | 0.9 |
| III | semi-arid and arid | ≤400 mm | 0.1440 | 0.25 | 1 |
The final regional frequency model for calculating peak flows in Moroccan territory is given by the following form:
Q T = Q 0 exp [ λ ⋅ log ( T ) ] (8)
or:
Q T = β P ( S ) n exp [ λ ⋅ log ( T ) ] (9)
where:
QT = peak flows (m3/s), of return period T (years);
P = average annual rainfall received by the watershed (mm);
S = area of the watershed (km2);
The regional parameters β, n and λ are tabulated according to the climatic zone and the average annual rain P received by the basin (
The characterization of the performance of the established regional model and the precision of the uncertainties linked to the use of the model are carried out by comparing the values of the flood flows observed and calculated by the regionalized model.
The comparison shows a good correlation between these values, it is translated by a correlation coefficient of the order of R = 0.9.
The values of the flows calculated for wadi Ourika (Aghbalou) deviate appreciably from the values provided by the Agency of the Tensift Hydraulic Basin
| River (Station) | Area S, km2 | Annual rain P, mm | Regional parameters ( | Flood flow index Q0, m3/s (6) | Frequency peak flow QT, m3/s C a l c u l a t e d ( 8 ) O b s e r v e d | ||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| β | n | λ | 10 | 20 | 50 | 100 | 1000 | ||||
| Mekhazine Bg. Mekh | 1820 | 1130 | 0.0025 | 0.8 | 0.6 | 1146 | 2088 1288 | 2501 2179 | 3175 3156 | 3804 4416 | 6931 11 , 686 |
| Rheraya (Tahanaout) | 225 | 670 | 0.0103 | 0.5 | 0.9 | 104 | 255 150 | 334 223 | 478 345 | 626 460 | 1540 |
| N’Fis PK89/PK10 | 848 | 700 | 0.0103 | 0.5 | 0.9 | 210 | 517 525 | 670 725 | 969 1000 | 1270 1250 | 3125 |
| Ourika (Aghbalou) | 503 | 600 | 0.0103 | 0.5 | 0.9 | 139 | 341 490 | 447 730 | 639 1160 | 838 1680 | 2062 |
| K’Sob (Adamna) | 1480 | 330 | 0.144 | 0.25 | 1 | 295 | 802 880 | 1084 1190 | 1613 1620 | 2180 1940 | 5925 |
| Moulouya Bg. Homadi | 51,500 | 200 | 0.144 | 0.25 | 1 | 434 | 1179 1500 | 1594 | 2372 | 3206 3000 | 8714 7500 |
(ABHT) [
By the way, the hydrological characteristics of this catastrophic flood (
The relative mean calculation error by the regional model (8) does not exceed ±8%, it can reach or sometimes exceed ±30% in the semi-arid and arid zone.
Usually, the difference between the calculated and observed values increases in the following cases:
· The gradual transition from the wetland to the semi-arid and arid zone.
· The incorrect determination of the average annual rain over the entire catchment area.
The distribution of the Galton law (5) adopted as an analytical representation to derive a unit hydrograph actually has two following drawbacks:
· The peak of the flood is reached well before time t = t m (
· The shape parameter σ = 0.4 is fixed (
To avoid these drawbacks, that is to, say fix the point in the desirable place and allow the hydrograph to marry several forms without deformation, we introduced a slight modification of the formulation of the Galton law and propose the expression (10) which will be called synthetic hydrograph type Galton:
Q ( t ) = Q T ( t t m ) − 0.1 exp [ − 1 2 ( ln ( t t m ) k ) 2 ] (10)
| Date of flood | Peak flow m3/s | Volume 103 m3 | Basic time (h) | Rise time (h) | Peak coefficient |
|---|---|---|---|---|---|
| December 02, 1987 | 650 | 34.36 | 39.5 | 14 | 2.7 |
| July 14, 1989 | 823 | 34.42 | 23 | 15 | 2 |
| August 17, 1995 | 1030 | 3.6 | 4 | 0.16 (≈10 min) | 4.1 |
| October 28, 1999 | 762 | 25.99 | 22 | 9.5 | 2.5 |
where:
Q ( t ) = flood flow in m3/s at time t in hours;
Q T = peak flow in m3/s, it is determined by the formula (8);
t m = flood rise time in hours;
k = shape parameter of the flood hydrograph;
The flood rise time t m is generally comparable to the concentration time t c for short-term showers (less than 24 hours) [
The concentration time for the case of the watersheds of the kingdom of Morocco can be calculated by similarity by the relationship proposed for the Algerian basins [
t c = 1.06 ( S L I ) 0.19 (11)
where:
tc = concentration time in hours;
S = area of the watershed in km2;
L—length of the main river in km;
I—average slope of the river in %.
The relation (11) provides values of the concentration time close to those presented by the ABH (
The shape parameter k of the flood hydrograph is given by the following relation [
k = 0.0102 ( S + 1 ) 0.4 + 0.15 (12)
where: S is area of the watershed in km2.
Verification of the simulation quality by the proposed model was based on the comparison of the hydrographs measured with those calculated by formula (10). The relationship between the calculated and measured values is quite good and is expressed by a correlation coefficient, varying between 0.7 and 0.8 (
| River | Catchment area S, km2 | Length of the main river L, km | Average slope of the river I | Concentration time tc | Reference | ||
|---|---|---|---|---|---|---|---|
| m/m | % | Calculated by (11) Hour or (Hour: min) | Observed (Hour: min) | ||||
| RDAT | 569 | 50 | 0.032 | 3.2 | 6.7 (6 h 42 min) | 5 h 30 min | [ |
| RHERAYA | 225 | 32 | 0.072 | 7.2 | 4.8 (4 h 48 min) | 4 h 30 min | |
| N’FIS | 1270 | 82 | 0.025 | 2.5 | 8.7 (8 h 42 min) | 8 h 30 min | |
| OURIKA | 503 | 45.5 | 0.069 | 6.9 | 5.9 (5 h 54 min) | 5 h 00 min | |
| ZAT | 528 | 55 | 0.048 | 4.8 | 6.4 (6 h 24 min) | 6 h 00 min | |
For the study of protection of the city of Fès against floods, we propose to calculate the peak flow for T = 100 years and its corresponding hydrograph. The basin
data are: the area of the watershed, S = 879 km2; the length of the wadi L = 70 km; the average slope of the wadi I = 1.44%. The average annual rain over the entire basin P = 550 mm. The characteristics of the watershed are taken from reference [
1) The average annual flood flow (flood index) by the formula (6):
Q 0 = β P ( S ) n = 0.0103 × 550 × ( 879 ) 0.5 = 168 m 3 /s
2) The peak flow of the 100-year flood by the formula (8):
Q 100 = Q 0 exp [ λ ⋅ log ( T ) ] = 168 exp [ 0.9 log ( 100 ) ] = 1016 m 3 /s
3) The flood hydrogram by (10):
Q ( t ) = Q 100 ( t t m ) − 0.1 exp [ − 1 2 ( ln ( t t m ) k ) 2 ]
4) The rise time tm by the formula (11):
t m = 1.06 ( S L I ) 0.19 = 1.06 ( 879 × 70 1.44 ) = 8.3 h o u r s
5) The form coefficient k by the formula (12):
k = 0.0102 ( S + 1 ) 0.4 + 0.15 = 0.0102 ( 879 + 1 ) 0.4 + 0.15 = 0.3
6) The equation of the flood hydrograph is therefore:
Q ( t ) = 1016 ( t 8.3 ) − 0.1 exp [ − 1 2 ( ln ( t 8.3 ) 0.3 ) 2 ]
We propose some values of the time (t > 0) before and after the rise time (for example 2, 4, 8.3, 14, 16, 20, 25 et 30 hours) then we calculate the corresponding flow rates (0.015, 57, 1016, 211, 87, 13, 1, et 0 m3/s). The plot of the flood hydrogram of river Fès is presented in
Unfortunately, we do not have flood values measured on wadi Fès to validate the value of the 100-year flood (1016 m3/s) calculated. Referring to the value of the 100-year flood of the Oued N’Fis watershed at PK89/PK10 (1250 m3/s) which is almost the same size (848 km2) with significantly higher rainfall (700 mm) (
The design of flood protection works requires the knowledge of the value of the flood flow, which unfortunately often remains unknown in several regions of Moroccan territory, due to lack of measurement data and/or absence of universal methods calculation. The regionalized model for calculating the peak flow of the flood established for the Moroccan territory makes it possible to provide values characteristic of the flood regime of ungauged rivers, using only readily available characteristics, namely: the area of the watershed, the average annual rainfall received by the latter and some regional parameters, tabulated according to the average annual rainfall.
The synthetic hydrograph type Galton developed will also make it possible to derive flood hydrographs for ungauged basins throughout the Moroccan territory from the knowledge of the following hydrological characteristics: the peak flow, the surface of the basin, the length and the average slope of the river (wadi). The formulas developed above are not very precise but provide representative values for the physiographic and climatic conditions of the watersheds of the kingdom of Morocco. As data is collected, the proposed methodology can be reviewed and corrected later.
The authors declare no conflicts of interest regarding the publication of this paper.
Djamel, B., Houari, Z. and Fares, B. (2020) Floods and Hydrograms of Floods of Rivers in Arid Zones of the Mediterranean, Case of the Kingdom of Morocco. International Journal of Geosciences, 11, 651-666. https://doi.org/10.4236/ijg.2020.1110033