Practical Method Proposed to Estimate Silting’s Rate in Small and Hillside Dams

The construction of small dams and hillside reservoirs is a common practice for dealing with erratic and unevenly distributed rainfalls in arid and semi arid areas. Moreover, these small hydraulics structures are recommended as a solution to adapt to climate change in rural areas. They are fundamental for the sustainability of the water resources. However, these dams are potential sinks for upstream sediment. Therefore, it would be useful to monitor and measure the amount of sediments in small dams because it is necessary to understand siltation and better functioning of these vital structures. This control and these measurements require the application of the bathymetry or topography surveys which is too expensive for the responsible of these small dams. In this context, this paper proposes a simple and low cost method called: Triangulation of the estimated silted level. This method is based on direct measurement of the accessible silted level in dam. This new method is compared with another method based on the initial Elevation-Capacity curve of dam, usually used when the bathymetry is not attainable. The results of both of them have been verified with those of the bathymetry and electrical tomography methods and validated the proposal method as the most approximate. Therefore, this method is recommended as a new direct and simple technique to evaluate the amount of silting in small dams and hillside reservoirs.


Introduction
Since their conception in the United States, improving their technique by Italians and, due to their low cost and their adaptation to marginal and mountainous areas, the strategy of building hillside and small dams was quickly

Direct Measurement Techniques Commonly Used for Calculate Rate's Siltation
Direct measurement of silting's rate is based on visiting the study area and survey the current state without using the modeling and empirical formulas.

Bathymetry Survey
The bathymetric approach is based on a simple comparison of reservoirs morphology at two different time periods, first at the time of the construction of the dam and second at the time of the survey. Bathymetric map involves two types of necessary measures: the position and depth measurements [14] [15] [35].

Remote Sensing
It based on the repetition of similar satellite images at regular intervals of 2 to 5 years on areas of dams. The remote sensing technique helps to calculate the present storage capacity, which is compared with the original capacity and reduction in capacity over a period of time is attributed to silt deposition. Remote sensing technique for estimating sedimentation is useful only for lakes and reservoirs where rise and fall in water level is a regular phenomenon [9] [36].

Calculating the Volume of Siltation
We can find different methods to operate the results of the bathymetry survey and remote sensing. All of them are based on the old methods for calculating volumes between profiles raised.

Kolmogorov Method
This method is based on the average between two sections survey profiles Pa, Pb and a fictitious profile Pc, it is explained on Figure 1 and Equations (1), (2) and (3) [14].

The Method of Average Height
As its name suggests, this method is based on the average silted heights (Equation (4)) on surveys profiles Pa and Pb and the surface that separates them Sab which is shows on Figure 2 [14].

The General Method
This method is based on silted sections on the surveys profiles but also takes into account the inclination of sections Pa and Pb on Figure 3 [15].
Now, with the development of topographic software and remote sensing, surveying of the silted areas are automatically integrated into specialized software to calculate directly and easily the volume of sediment in dams.

Technique to Estimate the Silting Volume: The Initial Elevation-Area-Capacity Curve of a Dam
In fact, the area-capacity curve is necessary for defining the storage capacity of a reservoir. An area capacity curve is obtained by planimetering the area enclosed within each contour line in the reservoir area before its filling (see Figure 4). The graphical plotting of area and capacity curves rates the surface area and the storage capacity of the reservoir to the elevation of the water surface (see Table 1 and Figure 5).    To calculate the silting volume in small dams when the survey bathymetry is not possible technically and financially, the manager of these structures applies generally the initial elevation-area-capacity curve of the dam. One or two employees surveyed the level of silt just on the embankment, and then we projected this level on the elevation-capacity curve which informed us about the volume of "silt". But this process assumes that the silting is considered like water. It means that the silting move horizontally on the surface of the dam's basin which is so doubting.

The New Proposal Method: Triangulation of the Estimated Silting Level
The application of bathymetric surveying and Remote Sensing in small dams and hillside lakes is not regularly possible. Therefore, a simple and less expensive method is necessary to estimate silting and monitoring its development in small dams. The proposed method consists to build an estimated digital elevation model based on TIN (triangulated irregular network). A TIN is a vector-based representation of the physical land surface or lake bottom, made up of irregularly distributed nodes and lines with three-dimensional coordinates (x, y, and z) that are arranged in a network of no overlapping triangles.
To prepare this estimated digital elevation model for silted surface in dam, we had surveyed a few accessible points on dam's surface: Figure 6 shows the three surveyed points on the embankment and other accessible points on the surface of silt, although, summer is recommended to survey silting level because of the decrease of water's level.
Actually, with these few points, we can establish an estimated digital elevation model, and then we can superpose on it the horizontal plan corresponding to the level of normal capacity of the dam which is a data. The difference of the horizontal plan and the estimated digital elevation model represents the new capacity of the dam after the siltation (see Figure 7 and Figure 8). Thus; the volume of the sediment is the difference of the initial capacity and the new calculated capacity.

Study Area
The study area consists on seven dams located on the watershed of Souss Massa in center of Morocco (see Figure 9).
The climate is arid to semi-arid. The proximity of the ocean and the influence of the cold current of the Canary Islands mitigate the climate of the area, the mountain barrier of the Anti-Atlas contributes to its protection against the winds from the south. Rainfall is unevenly distributed on the area. The average is about 280 mm and 265 mm in Souss Massa.
The water surface is limited and very irregular. Flows of rivers are erratically within the year. The runoff is only for short periods when it is flooding; sometimes it is so fast and violent. The water potential of the region is estimated at 600 Mm 3 /year. The minimum and maximum flows recorded are 35 million m 3 (1960 to 1961) and 2160 million m 3 (1962-1963) (see Figure 10).
The population is currently around 1,680,000 habitants; nearly 53% are in rural areas which have prompted the country, and some NGO (non-governmental organization) to move towards building small dams as a way to satisfy water requirements for rural communities. Today, there are more than 24 hillside dams in the region of   Souss Massa and it is planned to build more, due to the national program reaching 1000 dams and hill lakes throughout the country. Table 2 shows some characteristics of the studied dams and their watersheds.

The Fieldwork
The equipment required is:   The level of the silt throughout the perimeter of the lake and some accessible points of the dam's area (dry area) are surveyed by GPS. These measured elevations are corrected by the level of the water which has anaccurate elevation surveyed by the wire-weight gauge. We use this correction because of the inaccuracy of the GPS at elevations.

Data Processing
The process of calculating the volume of sediments on dams consists on seven steps below: a) Import of the surveyed points by Gps on the software like Covadis or surfer b) Addition of the three points surveyed on the embankment. c) Building the contour representing the level of the useful dam's capacity. d) Establishment of the digital elevation model with the surveyed points e) Superposition of the horizontal plan corresponding to the useful capacity f) Attribution the elevation of the horizontal plan using the initial curve of elevation capacity of the dam. g) Calculate the existing volume between the horizontal plan and the digital elevation model (see Figure 11).
Finally, the calculated volume is the new useful capacity of dam after siltation, and then we can deduce the volume of silting by the difference of the new and initial useful capacity.
The application of this technique is very affordable compared to other known reliable methods. Nevertheless, it is imperative to check the purpose of validation.

Control and Validation of the Proposed Method
The control and validation of the proposed method involved the comparison of its results with those of the following methods: -The bathymetry survey in the three dams Timicha, Taguenza and Assgherkiys (Figure 12) -The topographic survey in Tildi (dry lake) Figure 11. Steps to calculate the sediment volume. -The electrical tomography survey in Tigmi (dry lake). The above methods are reliable and made on five of seven dams studied, this control gives us a fairly acceptable margin for comparison, control and validation of the new proposed method.

The Bathymetry Survey in the Three Dams Timicha, Taguenza and Assgherkiys
The bathymetry survey is made by the Agency hydraulic basin Souss Massa Daraa. Data of these surveys are retrieved and processed by the Surfer software to deduce the volume of the silting (Figure 12).

The Topography Survey in Tildi Dam
For this dam, it was dry, so we didn't need bathymetry method; the topography survey using a total station was enough.
Surveys are imported from the total station to a computer, and saved in the specific format available for reading points on the software (Covadis, Surfer or other) to build a plan, showed on Figure 13.

The Electrical Tomography Survey in Tigmi Dam
This method is mainly used for the detection of noble metals, because it is too expensive. But due to the efforts of the team in our laboratory, we had the chance to apply this method in the Tigmi dam which was dry.
In fact, this method is based on the method of electrical resistivity which involves the introduction of an artificially generated electric direct current generally of low frequency (I) according to the source in four electrodes [37].
The implementation of four profiles in the dam, allowed having the resistivity and thickness of each geological bed in underground toward the profiles. Figure 14 and Figure 15 show the position of surveyed profiles and measures.    Table 3 and Figure 16 show the measurement results of siltation by the proposed new method Tesl (Triangulation of the estimated silting level), the method of the initial elevation capacity curve and the results of validation methods.

Results and Discussions
The comparison of the new method results and those of the initial curve commonly used, required the calculation of differences versus validation methods results. These later are shown in Table 4. Table 4 shows that the differences between the results obtained by the initial elevation capacity curve method, and those obtained by the validation method (bathymetry) is so significant, the initial elevation capacity curve method underestimates the volume up to 40%. This underestimation decreases if the siltation rate is important. For example, Tildi dam is silted over 90% and showed less difference. This may be explained by silting transport which was by homogeneous decantation and suspension over the entire surface of a lake that induces planar topography of the silt's surface. This sort of phenomenon is observed in the watersheds with friable lithology as marly limestone (case of Tildi and Tigmi dams). Contrariwise siltation of Assgherkiss and Timichalakes is the result thrusting transport which is generally low and induces the formation of deltas at the entrance of the lake and therefore an irregular surface.
The differences between the results obtained by the new method (Tesl) and those given by the validation methods are acceptable with a slight overestimation does not exceed 2.41% maximum. These results appear very encouraging and very reliable especially because there is no difference between lakes where sedimentation occurs a flat surface and those receiving the materials transported by thrusting and induces an irregular surface.
Generally, the deposition begins at the entrance of the dam by large material by thrusting and forms a delta. Afterward, the smaller materials and fine particles are carried by density currents through the embankment. Figure 17 shows the process of sedimentation of a dam [38].
The patterns bellow on Figure 18, show the portion of the silted volume by the initial curve method and by the new proposed method. It is noted that the initial method ignores the deltas created by silting at the entrance of the reservoir and the slope of the reservoir.
The new proposed method gives results closer to those of validation methods that are more accurate. This new method is also applicable on both of planer and irregular surfaces, then it can be considered technically a reliable method and deserves to be applied instead of the commonly used method of the initial elevation capacity curve.
The advantage of the new method is that its cost ten to twenty times less than the bathymetry survey. Therefore, it is recommended for users of small dams and hillside dams who generally have limited financial and technical resources.
The new technique has a scientific contribution since it allows the estimation of sediment transport in the sub watersheds and used in the screening of future sites of small dams and hillside lakes.

Conclusions
Small dams and hillside lakes which are usually located in fragile environments with low economic activity play a vital role in the sustainable development of these areas. Arid and semi-arid areas are most affected by recurrent drought periods, hence followed by the construction of these dams for water conservation. However, these small hydraulic structures are greatly threatened by siltation, which limits their durability and makes their aim uncertain. In fact, the scourge of the dams silting is not new. It tormented the managers of these structures due to antiquity and will continue to do it. This is a natural phenomenon; therefore, the prevention remains the most efficient way to mitigate it. The measurement of siltation in dams is an important step to understand and master its operation. This direct measurement of siltation is not always easy and attainable to the responsible of small dams and hillside lakes. In this context, we would like to propose a new method of measuring siltation, less expensive than those currently used. This new method is validated by the more precise method as bathymetry survey and others recognized methods in five hillside dams. Thus, this method could be recommended for the measurement of siltation in the small dams and hillside lakes especially by the NGOs and communities who have limited technical and financial resources.
The adoption of this method will allow studying a larger number of small dams and hillside lakes sites. The measures will help understand the specificity of each sub-basin, then to identify the significant indicators of siltation. Finally, we can use these detected indicators for studying the proposal best sites for the future small dams.