Floods and Hydrograms of Floods of Rivers in Arid Zones of the Mediterranean, Case of the Kingdom of Morocco ()
1. Introduction
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 [1].
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 [2].
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 [3]. The damage caused by the torrential rains of 2008 and 2009 on the road network (roads and civil engineering works) represents more than 2.9 billion dirhams [4]. Development of flood protection projects, sizing of road crossing structures and spillways etc. is always linked to the qualitative assessment of extreme flows (floods) and their flood hydrographs.
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 [5]. The flood index is generally a stable value for each watershed (it represents the inter-annual average flood) and is easy to relate to the pluviometric and physical parameters of the watershed [6] [7].
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 [6] [8]. In the absence of data, the maximum annual flood can be evaluated using simple formulas including meteorological, physiographic, morphological data, etc. of the watershed [6] [7] [9].
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 [5] [9].
2. Material and Methods
2.1. Material
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 (Table 1), they are distributed almost uni-formly over all the hydrographic basins of Morocco (Tensift, Moulouya, Sebou, Sous massa, Oum Er Rebia, Guir, Ziz Ghéris and Mediterranean coastal basins). The observation period varies from one station to another, it varies on average between 15 and 20 years (1970-2010).
2.2. Methods
2.2.1. Peak Flood Flows
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 [9]. From a statistical point of view, the region is said to be homogeneous, if the observations of the different sites in this region come from the same parent population [5].
Table 1. Number of watersheds according to their areas.
There are several methods for delimiting homogeneous regions, the most used are based on geographic criteria, physiographic and/or climatic criteria and hydrological criteria [8].
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 (Figure 1), where: P is average annual rainfall.
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 [7]:
(1)
where:
is a scale factor representing the average of the maximum annual floods, it is also called the flood index;
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.
Figure 1. Relation between the specific flow q0 (m3/s∙km2) and the area of the watershed S (km2) in different climatic zones.
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 [10].
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:
(2)
where:
is standardized frequency distribution;
is return period flood flow T;
is flood index.
The mathematical expression
for the three homogeneous regions is the same, it has the following form:
(3)
For the three homogeneous regions, the value of the parameter ξ is close to unity
. 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 [6] propose a simple method to determine the flood index (maximum annual average flood Q0) according to the area of the watershed S and its geographic position for a slice of rain given by the following relation:
(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 (Table 2).
The formula of R. Hazen and D. Lazarvic for calculating the flood index Q0 therefore requires corrections.
2.2.2. Flood Hydrographs
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:
Table 2. Average flood flows of river (wadi) R’dat and Oued Zat.
· 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.
Figure 2 present an example of the unit hydrogram determination of the flood that occurred on October 23, 2010, on the Ghdat wadi (Sidi Rahal).
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 Table 3.
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 Figure 3, where: Qt is flood flow at time t; Qp is peak flood flow;
is the rise time.
Figure 2. Dissociation of the flood hydrograph in surface and underground flow.
Table 3. The ordinates (in % of the total runoff volume) of the distribution diagram of the unit flood hydrographs at the Tahanaout station (Oued Rheraya).
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 [3]:
(5)
The analysis of this mathematical formulation by setting
and
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
(Figure 4).
3. Result and Discussion
3.1. Peak Flows
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
Figure 4. Distribution of the adimensional hydrograph as a function of σ.
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:
(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 (Table 4).
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 (Figure 5).
The standardized regional frequency distribution model for the three homogeneous regions of Moroccan territory is as follows:
(7)
The parameter λ is tabulated according to the climatic zone and the average annual rainfall (Table 4).
Table 4. Values of the regional parameters β, n and λ.
Figure 5. Comparison between the calculated and observed flood index values.
The final regional frequency model for calculating peak flows in Moroccan territory is given by the following form:
(8)
or:
(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 (Table 4).
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.
Table 5 presents some examples of calculation of the flood flows of some rivers in the Moroccan territory.
The values of the flows calculated for wadi Ourika (Aghbalou) deviate appreciably from the values provided by the Agency of the Tensift Hydraulic Basin
Table 5. Calculation of peak flows of some wadis in Moroccan territory.
(ABHT) [11]. The difference is certainly due to the update of the frequency study of flows by incorporating the deadly and catastrophic flood (1030 m3/s) which was recorded on August 17, 1995, this flood does not come from the same population (natural flood flows), it is the result of an aggravated natural phenomenon (sudden rupture of the enormous temporary natural dam formed by a large stone plug upstream of the Ourika valley) [12] [13].
By the way, the hydrological characteristics of this catastrophic flood (Table 6) are very different from the characteristics of other floods recorded before and after 1995 [14].
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.
3.2 Flood Hydrographs
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
(Figure 4).
· The shape parameter σ = 0.4 is fixed (Figure 4). For other values of σ, the hydrographs will not have the meaning of a normalized dimensionless hydrograph.
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:
(10)
Table 6. Characteristics of certain floods of the Ourika river at Aghbalou [14].
where:
= flood flow in m3/s at time t in hours;
= peak flow in m3/s, it is determined by the formula (8);
= flood rise time in hours;
= shape parameter of the flood hydrograph;
The flood rise time
is generally comparable to the concentration time
for short-term showers (less than 24 hours) [15]. Which are characteristic of arid and semi-arid climates.
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 [16]. With a slight modification which takes into account the specificity of the Moroccan relief, it is of the following form:
(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 (Table 7).
The shape parameter k of the flood hydrograph is given by the following relation [16]:
(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 (Figure 6).
Table 7. Comparison of calculated and observed concentration time values.
4. Application
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 [17].
1) The average annual flood flow (flood index) by the formula (6):
2) The peak flow of the 100-year flood by the formula (8):
3) The flood hydrogram by (10):
4) The rise time tm by the formula (11):
5) The form coefficient k by the formula (12):
6) The equation of the flood hydrograph is therefore:
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 Figure 7.
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) (Table 5) than that of Fès and being in the same homogeneous region, the value of the flood calculated for wadi Fès can be considered reasonable.
Figure 7. Hydrograph of the centennial flood (Q100) of river Fès (Morocco).
5. Conclusions
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.