A Methodology to Assess and Evaluate Rainwater Harvesting Techniques in (Semi-) Arid Regions

Arid and semi-arid regions around the world face water scarcity problems due to lack of precipitation and unpredictable rainfall patterns. For thousands of years, rainwater harvesting (RWH) techniques have been applied to cope with water scarcity. Researchers have used many different methodologies for determining suitable sites and techniques for RWH. However, limited attention has been given to the evaluation of RWH structure performance. The aim of this research was to design a scientifically-based, generally applicable methodology to better evaluate the performance of existing RWH techniques in (semi-) arid regions. The methodology integrates engineering, biophysical and socio-economic criteria using the Analytical Hierarchy Process (AHP) supported by the Geographic Information System (GIS). Jessour/Tabias are the most traditional RWH techniques in the Oum Zessar watershed in south-eastern Tunisia, which were used to test this evaluation tool. Fifty-eight RWH locations (14 jessr and 44 tabia) in three main sub-catchments of the watershed were assessed and evaluated. Based on the criteria selected, more than 95% of the assessed sites received low or moderate suitability scores, with only two sites receiving high suitability scores. This integrated methodology, which is highly flexible, saves time and costs, is easy to adapt to different regions and can support designers and decision makers aiming to improve the performance of existing and new RWH sites.


Introduction
Aridity and climate change are the major challenges faced by farmers who rely on rainfed farming [1]. Especially in arid regions, farmers are faced with low average annual rainfall and variability in temporal and spatial distribution. In order to increase the availability of water for crop production and cattle grazing, inhabitants of dry areas have constructed and developed several types of Rain Water Harvesting techniques (RWH). RWH is a method for inducing, collecting, storing and conserving local surface runoff for agriculture in arid and semi-arid regions [2]. RWH is a likely viable option to increase water productivity at the production system level [3]. RWH and management techniques have a significant potential for improving and sustaining the rainfed agriculture in the region [4]. In fact, a wide variety of micro-catchment, macro-catchment and in situ RWH techniques are available in arid and semi-arid regions. The indigenous techniques, or those modified by the indigenous RWH practices, are more common and widely accepted by smallholder farmers than the others [5]. Throughout history, archaeological evidence has revealed RWH sites that were implemented in Jordan, the Al-Negev desert, Syria, Tunisia and Iraq. The earliest signs of RWH are believed Several types of RWH exist in the study area to satisfy water requirements for agriculture and ground water recharge. The most common RWH systems in the region are jessour and tabias; spreading of flood water and groundwater recharge structures in the wadi beds are applied too [23]. To test the RWH evaluation tool, three representative sub-catchments were selected based on four criteria. i.
Representative of the geographic distribution of our watershed; one located in upstream another in the midstream and one in downstream. ii.
Representative of the different types (jessour and tabias), scale (small and large) and age of RWH systems (new and old). iii.
Source and destination of collected rainwater for each sub-catchment. iv.
Accessibility; easy to access physically and acceptance of the local people.
These three sub-catchments are located in the downstream (Sub-catchment 1), middle (Sub-catchment 2), and upstream (Sub-catchment 3) of Oum Zessar watershed as shown in Figure 1. Each jessr (singular of jessour) or tabia consists of three parts: the impluvium or catchment area providing the runoff water; the terrace or cultivation area where the runoff water is collected and crops or trees are grown; and the dyke, which is a barrier to catch water and sediment. Each dyke has a spillway (menfes if the spillway is located on one or both sides and masref if the spillway is located in the middle of the dyke) to regulate water flow between dykes (see Figure 2). Several types of RWH exist in the study area to satisfy water requirements for agriculture and ground water recharge. The most common RWH systems in the region are jessour and tabias; spreading of flood water and groundwater recharge structures in the wadi beds are applied too [23].
To test the RWH evaluation tool, three representative sub-catchments were selected based on four criteria. i Representative of the geographic distribution of our watershed; one located in upstream another in the midstream and one in downstream. ii Representative of the different types (jessour and tabias), scale (small and large) and age of RWH systems (new and old). iii Source and destination of collected rainwater for each sub-catchment. iv Accessibility; easy to access physically and acceptance of the local people.
These three sub-catchments are located in the downstream (Sub-catchment 1), middle (Sub-catchment 2), and upstream (Sub-catchment 3) of Oum Zessar watershed as shown in Figure 1. Each jessr (singular of jessour) or tabia consists of three parts: the impluvium or catchment area providing the runoff water; the terrace or cultivation area where the runoff water is collected and crops or trees are grown; and the dyke, which is a barrier to catch water and sediment. Each dyke has a spillway (menfes if the spillway is located on one or both sides and masref if the spillway is located in the middle of the dyke) to regulate water flow between dykes (see Figure 2). Water 2016, 8,198 4 of 23

General Description of the RWH Evaluation Decision Support Tool
This research aims to develop a more comprehensive and relevant evaluation tool for RWH structures. To achieve this goal, we developed a simple and robust assessment tool for the evaluation of RWH sites (structures) which is inexpensive, simple to apply, reliable and flexible with different criteria and easy to adapt to various RWH techniques and regions. The Analytical Hierarchy Process (AHP) forms the base for this tool.
The AHP is a multi-criteria decision making method, providing a structured technique for organizing and analyzing complex decisions, based on mathematics and expert knowledge [24]. It was developed by Thomas Saaty in the 1970s and, since then, has been applied extensively in different disciplines. The main principle of AHP is representing the elements of any problem hierarchically to show the relationships between each level. The uppermost level is the main goal (objective) for resolving a problem and the lower levels are made up of the most important criteria that are related to the main objective. Pairwise comparison matrixes are constructed and scaled in preference from 1to 9 for each level. Then, the consistency of each matrix is checked through the calculation of a consistency ratio (cr). The cr should be smaller or equal to 10% [25]. The weight for each criterion and the cr are determined, then all matrixes are solved.

Methodology Overview
AHP is particularly useful in multi-index evaluation and consists in our RWH evaluation tool of the following steps: i.
Describe the main objective of the intervention; ii.
Identify the biophysical, engineering (technical) and socio-economical main and sub-criteria; iii.
Develop a decision hierarchy structure; iv.
Collect and process the data for each sub-criteria; v.
Classify the values for each sub-criteria in terms of suitability classes; vi.
Apply the pairwise comparison matrix to identify priorities (weights) for each criterion; vii.
Check the results with the stakeholders; and ix.
Decide based on conclusions and recommendations 2.3.1. Description of the Main Objective of the Intervention In our case study, the main objective is to collect and store runoff water during the rainy season to enable farmers to grow profitable crops and mitigate drought spells in arid and semi-arid regions.

General Description of the RWH Evaluation Decision Support Tool
This research aims to develop a more comprehensive and relevant evaluation tool for RWH structures. To achieve this goal, we developed a simple and robust assessment tool for the evaluation of RWH sites (structures) which is inexpensive, simple to apply, reliable and flexible with different criteria and easy to adapt to various RWH techniques and regions. The Analytical Hierarchy Process (AHP) forms the base for this tool.
The AHP is a multi-criteria decision making method, providing a structured technique for organizing and analyzing complex decisions, based on mathematics and expert knowledge [24]. It was developed by Thomas Saaty in the 1970s and, since then, has been applied extensively in different disciplines. The main principle of AHP is representing the elements of any problem hierarchically to show the relationships between each level. The uppermost level is the main goal (objective) for resolving a problem and the lower levels are made up of the most important criteria that are related to the main objective. Pairwise comparison matrixes are constructed and scaled in preference from 1 to 9 for each level. Then, the consistency of each matrix is checked through the calculation of a consistency ratio (cr). The cr should be smaller or equal to 10% [25]. The weight for each criterion and the cr are determined, then all matrixes are solved.

Methodology Overview
AHP is particularly useful in multi-index evaluation and consists in our RWH evaluation tool of the following steps: In our case study, the main objective is to collect and store runoff water during the rainy season to enable farmers to grow profitable crops and mitigate drought spells in arid and semi-arid regions.

Identification of the Main and Sub-Criteria
This step formulates the set of criteria for the assessment based on the main objective. All major aspects should be represented, but the set should be as small as possible (simple and flexible). In addition to engineering (technical) aspects, social and economic aspects should also be included. Furthermore, the set of criteria has to be operational (e.g., measurable) and not redundant (the set should not count an aspect more than once).
In this study, we looked for criteria that represent the key parameters affecting the performance of RWH interventions and which could be applied to different sites and techniques. The parameters we were concerned with were based on the general definition of RWH, i.e., a method for inducing, collecting, storing and conserving local surface runoff for agriculture in arid and semi-arid regions [2], and information found in literature studies. The main selected criteria and sub-criteria are shown in Figure 3, and reflect the following questions: After assignment of scores, the weight for each criterion was determined by applying AHP with the pairwise comparison matrix. Pairwise comparison concerns the relative importance of two criteria involved in determining the suitability for a given objective. A pairwise matrix is first made for the main decision criteria being used. Other pairwise matrixes are created for additional criteria levels. The comparison and rating between two criteria are conducted using a 9-point continuous  Sub-criteria were chosen based on the relation with the main criteria (above), field investigations, expert discussions and literature studies.

Development of the Decision Hierarchy Structure
In this step, the main criteria and sub-criteria are arranged in a multilevel hierarchical decision structure. In this study case, the objective of the RWH (jessour and tabias) represents the first level. The second level contains the main criteria for the assessment. These criteria define the aspects by which the intervention is assessed e.g., how it fits within the local conditions (climate, drainage length and landscape), functionality and reliability based on the engineering design, and socio-economic aspects. The sub-criteria used to measure the performance of each main criterion are represented in the third level. Figure 3 shows the structure of the applied methodology for our case study.

Collection and Processing of the Data for Each Sub-Criteria
The definition, data collection, field measurements, storage and processing of data, as well as the calculations used for each criterion is explained in detail in the Section 2.4.

Classification of the Values for Each Sub-Criteria in Terms of Suitability Classes
Due to the variety of measurements and scales for the different criteria, a comparable scale between criteria must be identified before applying AHP tools. For instance, rainfall depth is measured in mm while soil texture is measured by the percentage of clay content. Therefore, the selected criteria were re-classified into five suitability classes, namely, 5 (very high suitability), 4 (high suitability), 3 (medium suitability), 2 (low suitability), and 1 (very low suitability). For example, suitability Class 3 is considered to be acceptable performance, while suitability Class 1 means that the RWH does not work well and that one or all criteria that caused this insufficient performance need improvement. Table 1 shows the scores assigned based on discussions and consultations with experienced people and information found in the literature.

Application of Pairwise Comparison Matrix to Identify Priorities (Weights) for Each Criteria
After assignment of scores, the weight for each criterion was determined by applying AHP with the pairwise comparison matrix. Pairwise comparison concerns the relative importance of two criteria involved in determining the suitability for a given objective. A pairwise matrix is first made for the main decision criteria being used. Other pairwise matrixes are created for additional criteria levels. The comparison and rating between two criteria are conducted using a 9-point continuous scale, the odd values 1, 3, 5, 7, and 9 correspond respectively to equally, moderately, strongly, very strongly and extremely important criteria when compared to each other. The even values 2, 4, 6 and 8 are intermediate values [26]. During pairwise comparison, criteria were rated based on the literature review, information from the field survey and discussions with stakeholders and experts. The final weight calculation requires the computation of the principal eigenvector of the pairwise comparison matrix to produce a best-fit set of weights. The consistency of each matrix, which shows the degree of consistency that has been achieved by comparing the criteria, was checked through the calculation of consistency ratio (cr). The cr should be smaller or equal to 10%, otherwise they are judged as not consistent enough to generate weights and, therefore, have to be revised and improved [25]. Table 1. Classification, suitability levels and scores for each criterion for assessment of existing RWH sites in arid and semi-arid regions. Each value, class and score were rated based on the literature review, information from the field survey and discussions with stakeholders and experts. Storage capacity ratio (-), the ratio between the total volume of water inflow and existing storage capacity. The ratio that is close to one is ranked as highly suitable.

Criteria
Over requirement (too large a storage capacity area) <0. Structure dimensions ratio (-), the ratio between the required design height and the existing height of dykes or barriers for each RWH structure. The ratios that are close to one are ranked as highly suitable Reliability ratio (-), the ratio between the total demand and the total supply of water. High suitability scores for the ratio are close to one Sufficient (required water is largely less than supply) <0. 35  To find out the final weight for each criterion and the cr, we solved the pairwise matrixes mathematically. The results of the main criteria from the pairwise comparison and the final weight are presented in the results section.
In this study, two methods were applied. In the first, the hierarchy structure consists of all three levels; the objective, main criteria (5 criteria) and sub-criteria (11 criteria). In the second method, the hierarchy structure consists of just two levels: the objective and the sub-criteria (11 criteria). By applying these two methods, the understanding of the relation between each criterion and its reflection on the main objective becomes much clearer, and they confirm the flexibility of AHP to adopt different criteria on multi-levels. Moreover, this will give an insight into whether there are any mistakes and how they will be distributed or fixed, and gives more reliability and confidence in our methodology for adoption in different regions and/or for different criteria.

Calculation of the RWH Performance (Suitability)
The next step in the assessment methodology is the calculation of the overall suitability for each RWH site. The overall RWH suitability was calculated by applying the following formula: where: S: suitability; W: weight of criteria i; X: score of criteria i; i, n: number of criteria The overall suitability will be classified also from 1 to 5, namely, 5 (very high suitability), 4 (high suitability), 3 (medium suitability), 2 (low suitability) and 1 (very low suitability).

Discussion of the Results with Stakeholders
It is important to check the results with the stakeholders, including the preliminary conclusions and recommendations. If felt that something is missing or has changed, additional measurements or recalculation with different weights might be necessary. Thereafter, results have been presented again to the local stakeholders for discussion and approval.

Decision Making Based on Conclusions and Recommendations
The main results of the assessment will give insight into if and how a RWH structure can be improved to increase its performance. Once there is general agreement on the results between stakeholders and scientists, a well-founded decision can be made on what structure needs to be improved for better performance of the RWH system.

Data Collection
Different data sources were used. Meteorological as well as other biophysical data, was collected from the Institute des Régions Arides (IRA) in Tunisia. Field measurements were carried out in the Wadi Oum Zessar during the period from December 2013 through March 2014. An open structure interview was made with key stakeholders (41 landowners and farmers) and discussions with people working and having experience with RWH (15 experts), particularly the engineers from the Regional Department of Agriculture in Medenine. A pairwise matrix was established and the relative weights for each criterion and suitability rank for classes are assigned as shown in Table 1. GIS was also applied to extract data that are needed in our methodology. All collected and measured data were stored and processed using Excel software.

Climate and Drainage Data
Rainfall Rainfall is one of the major components in any RWH system, with the magnitude of rainfall playing a significant role in assessing the RWH suitability for a given area. In arid and semi-arid regions, rainfall varies greatly in time and space. RWH systems can only function if there is sufficient rainfall in the catchment area to be stored somehow. Average monthly rainfall for the period 1979-2004 was collected from IRA for 7 meteorological stations in the Wadi Oum Zessar watershed, namely Ben Khedache, Toujan Edkhile, Allamat, Koutine, Sidi Makhlouf, Ksar Hallouf and Ksar Jedid. The rainfall amount in the three test sub-catchments was determined by applying the Inverse Distance Weight (IDW) function from ArcGIS 10.0 to interpolate the data from these stations. The rainfall depth data was then reclassified and scored as shown in Table 1. Areas with high annual rainfall are ranked as highly suitable.

Drainage Length
Since RWH interventions (especially jessour and tabias) are located on the hydrographic network and their location is influenced by topography, the distance from the water course has a significant role in the assessment of RWH performance. In this study, the distance from a RWH site to the drainage networks is used to represent the runoff suitability. By determining the location of the furthest point contributing to runoff [31], the drainage system was classified to each of the RWH sites (short distance means fewer water losses). The distances from the water courses to each dyke were measured using Google earth image and ArcGIS software.

Storage Capacity
One of the main principles of RWH is storing water to mitigate drought effects in dry seasons. Technically, the volume of water harvested and the amount retained over a reasonable duration of time is one indicator of the performance of RWH.
Potential runoff (V 1 in m 3 ) from a catchment area was calculated by: C: The mean annual runoff coefficient (-); equal (0.18) based on the simulations done by Schiettecatte et al. [32]. P: The mean annual precipitation (mm) A: The catchment area (m 2 ) The total volume of water inflow (V i ) is, therefore: where V 2 (m 3 ) is the overflow from upstream dyke(s) and V 3 (m 3 ) is the volume of rainfall onto the storage area. During the field measurements, the retention area and maximum potential depth of water (height of spillway) were measured with GPS and measuring tape. Then, the existing storage volumes were calculated (by multiplying the retention area by spillway height). Finally, the ratio between the total volume of water inflow (V i ) and existing storage capacity were calculated and scored. If the ratio, for example, equals 1-2, it means that the total inflow volume will be similar to the storage capacity or there is excess water that will be an overflow to the downstream. Therefore, the ratios that are close to one are ranked as highly suitable (Table 1).

Structure Dimensions
The dimensions of RWH structures are very important for achieving stability, controlling flood hazard and water supply. Furthermore, the primary goal of a structure is to harvest water for irrigation crops; the secondary goal is for flood protection. In this study, we assessed the existing height of dykes or barriers for each RWH structure and then compared this with the theoretical (required) design height.
The existing dyke's height for each site was measured in the field. The total volume of water that could be collected behind each dyke was calculated as noted in the previous section. The effective dyke height was calculated using this information. The free board, the vertical distance between the top of the dam and the full supply level, was calculated using standard dam design principles and added to the effective dyke height to determine the theoretical design height for each site. The ratio between existing and design dyke height was calculated and scored, as shown in Table 1.

Catchment to Cropping Area
To provide sufficient water to the crops, the terrace area should be not too large and the impluvium area should be enough. Therefore, an optimal ratio between impluvium area and terrace area has to be found. Depending on effective rainfall and runoff rates, the ratio between the catchment (impluvium) and cropping (terrace) area (Ca/C) can be determined. According to Schiettecatte et al. [32], the minimum ratio (Ca/C) "impluvium area/terrace area" (design) can be calculated by: where WR is the annual crop water requirement, P is the average annual precipitation (mm) for the period 1979-2004, and C is the average annual runoff coefficient (0.18) of dry soil and wet soil which was measured by Schiettecatte et al. [32]. For olive trees, the WR is estimated to be 500 mm¨y´1 [32]. Catchment area (impluvium) and cropping area were delineated with GPS in the field, and the areas were calculated using ArcGIS. At the end, the CCR ratio between the design and existing "impluvium area/terrace area" were calculated and scored.

Site Characteristics
Soil Texture Soil texture is a very important factor in selecting, designing and assessing the performance of RWH. Soil texture affects both the infiltration rate and surface runoff. The textural class of a soil is determined by the percentage of sand, silt and clay. Soil texture also determines the rate at which water drains through a saturated soil; for instance, water moves more freely through sandy soils than it does through clayey soils. High infiltration rates such as with sandy soil are not suitable for RWH structure. Clay soils have a greater water holding capacity than sandy soils, therefore, soil with high water holding capacity are more suitable for RWH. Indeed, Mbilinyi et al. [33] and others conclude that clay soil is best for water storage due to its low permeability and ability to hold the harvested water.
In this research, the terrace area was sampled at different sites (based on the size of terrace area, 1-3 samples for each site) and at depths of up to 1.3 m. The samples were taken to the IRA laboratory and analyzed. The clay contents (%) were measured, rated and classified into five suitability classes, as shown in Table 1.

Soil Depth
Soil should be deep enough to allow excavation to the prescribed depth for RWH, to ensure both adequate rooting development and storage of the harvested water. Critchley and Siegert [20] and Kahinda et al. [1] used soil depth as one criterion for selecting potential sites for RWH. Both soil depth and soil texture determine the total soil water storage capacity, which controls the availability of water for crops during the dry periods [9]. We measured soil depth in the field using a steel bar hammered into the ground until it could go no further and by checking the soil levels between two successive terraces. Then, soil depth data were categorized and classified into five suitability classes, as shown in Table 1.

Slope
Slope is also a major factor in site selection, implementation and assessment of RWH. It plays a significant role in runoff and sedimentation quantity, the speed of water flow and quantity of material required to construct the dyke structure (dyke's height).
Using DEM (30 m resolution) and ArcGIS 10.0, the slope was extracted for each catchment area and reclassified. Due to the large variety of slope values between jessour and tabias, different suitability classes were used for each type as shown in Table 1.

Structure Reliability
The relation between the demand and supply of water (reliability) is a good indicator of the performance of a RWH structure. Based on the function (purpose) of each technique, the demand for each RWH site was calculated. In our case, the main purpose of RWH is for on-site crop production.
The total demand was calculated by estimating the crop water requirements (evapotranspiration ET c ) plus losses to downward percolation, based on the field measurements by Schiettecatte et al. [32] in the same watershed.
The total demand " ET c`D ownward percolation Schiettecatte et al. [32] applied the Penman-Monteith method to calculate potential evapotranspiration (PET) and used data from the meteorological station at Medenine to calculate the average PET values over the period 1985-1995.
The maximum crop evapotranspiration (ET c ) was calculated by: where k c is the crop coefficient. Table 2 above shows the values for PET, ET c and k c . The infiltration ratios were used to calculate the downward percolation based on the soil texture results, as shown in Table 3 [34]. From the relation between storage capacity and total runoff volume from Equation (2), the total potential volume of supply water was calculated. Reliability was calculated as the ratio between total demand and the total supply of water for each site.

Socio-Economic Criteria
The success of an intervention depends not only on technical aspects but also on how well it fits within the stakeholder's social context and the economic benefit it provides him/her. Bamne [35], Al-Adamat [27] and Nasr [36] argued that one of the main reasons we do not use RWH sufficiently in the Middle East and North Africa is insufficient knowledge of the socio-economic contexts. There are several socio-economic criteria such as ownership, family size, education etc., and to identify good indicators for socio-economic conditions in relation to the functioning of these RWH systems is much more difficult than the biophysical ones. In this case study based on the literature studies and expert discussion, we are using distance to the settlements and cost per cubic meter of water as the socio-economic criteria influencing how suitable the intervention is for the main stakeholders.

Distance to Settlements
Since the local community is targeted in this study, the distance to the settlements is an important parameter in the design, selection and assessment of the RWH suitability [6]. We assumed that the distance to their home would influence the way they manage this system. Each farmer has scattered farming fields at a radius of about 20 m-1 km from his house. Therefore, it is very logical that the closer the field, the easier are the maintenance operations, particularly in the mountain zones where transportation is difficult. The distance for each site was measured using the image from Google earth and the ArcGIS program. Thereafter, as with other criteria, the values were reclassified and scored.

Cost Per Cubic Meter of Water
Cost plays a significant role in the design and assessment of RWH sites. In order to assess the cost effectiveness of each structure, the establishment and annual maintenance costs for each site were calculated. The actual costs for each structure were not available; the main problem with the jessour and tabia is that they do not have fixed designs (different shapes and sizes). Therefore, it is difficult to calculate the exact cost for each structure. Thus, the costs have been estimated by using the best available resources. The cost for each jessr or tabia was calculated based on the World Overview of Conservation Approaches and Technologies (WOCAT) database [30] and interviews with the local farmers. The costs for each jessr/tabia include the establishment and maintenance cost per year. The establishment costs consist of dyke construction, plantations, spillway construction for jessour and diversion channels and terracing for tabia. The maintenance costs consist of crop and tree maintenance, dyke and spillway maintenance, repairs and reconstruction. The overall costs for jessre per year are 3000 US$ for establishment and 900 US$ for maintenance. Whereas, 670 and 200 US$ for establishment and maintenence for tabia per year, respectively [30]. Based on the field measurements, the length for each jessr/tabia was measured and then the cost for each meter length of jessr/tabia was estimated.
These costs are similar to the values that were discussed with local farmers. The volume of collected water in each storage area and maintenance and construction costs of the jessour/tabias were used to calculate the cost per cubic meter of water, which was then classified and scored.

Application of the Assessment Tool for Different Test Sub-Catchments
We first tested our methodology on a catchment that has only one type of RWH structure. Sub-catchment one has just 17 tabias and no jessour and a total area of about 20 ha. It is located in the downstream area of the Oum Zessar watershed, as shown in Figure 1.
To further validate the methodology and criteria, we applied it on the other two sub-catchments, which have different characteristics. The second sub-catchment is located in the middle of Wadi Oum Zessar and has 16 RWH structures, 9 tabias followed (downstream) by 7 jessour, and a total area of about 19 ha. Sub-catchment three is located in the upstream part of Wadi Oum Zessar, with 25 RWH-8 jessour followed by 17 tabias-and a total area of about 45 ha.

Results
All the collected data for each site were stored and analyzed in Excel. The results for each criterion were then classified according to the five classes as defined in Table 1. Figure 4 shows the scores percentages (5 scores) of each sub-criteria (11 criteria) for all 58 sites. The rainfall criterion got a score 3 in all sites since there was no big difference in rainfall pattern nor amount (175-185 mm¨y´1) in the three sub-catchments due to the relatively small area. The criteria related to the design structure, like dimensions, storage capacity, CCR, drainage flow and costs got a high percentage of scores of 1 in many sites. More details about suitability and scores for the three sub-catchments are explained in the following sections. farmers. The volume of collected water in each storage area and maintenance and construction costs of the jessour/tabias were used to calculate the cost per cubic meter of water, which was then classified and scored.

Application of the Assessment Tool for Different Test Sub-Catchments
We first tested our methodology on a catchment that has only one type of RWH structure. Sub-catchment one has just 17 tabias and no jessour and a total area of about 20 ha. It is located in the downstream area of the Oum Zessar watershed, as shown in Figure 1.
To further validate the methodology and criteria, we applied it on the other two sub-catchments, which have different characteristics. The second sub-catchment is located in the middle of Wadi Oum Zessar and has 16 RWH structures, 9 tabias followed (downstream) by 7 jessour, and a total area of about 19 ha. Sub-catchment three is located in the upstream part of Wadi Oum Zessar, with 25 RWH-8 jessour followed by 17 tabias-and a total area of about 45 ha.

Results
All the collected data for each site were stored and analyzed in Excel. The results for each criterion were then classified according to the five classes as defined in Table 1. Figure 4 shows the scores percentages (5 scores) of each sub-criteria (11 criteria) for all 58 sites. The rainfall criterion got a score 3 in all sites since there was no big difference in rainfall pattern nor amount (175-185 mm y −1 ) in the three sub-catchments due to the relatively small area. The criteria related to the design structure, like dimensions, storage capacity, CCR, drainage flow and costs got a high percentage of scores of 1 in many sites. More details about suitability and scores for the three sub-catchments are explained in the following sections.

AHP & Suitability
During pairwise comparison, criteria were rated based on the literature review, interviews with key stakeholders, field survey information and discussions with people working and having experience with RWH, as shown in Table 4. For instance, the reliability and socio-economic criteria have similar relative importance to the main objective of the RWH system, as shown in this Table, and each of them has 1 as a relative importance rate.

AHP & Suitability
During pairwise comparison, criteria were rated based on the literature review, interviews with key stakeholders, field survey information and discussions with people working and having experience with RWH, as shown in Table 4. For instance, the reliability and socio-economic criteria have similar relative importance to the main objective of the RWH system, as shown in this Table, and each of them has 1 as a relative importance rate. A pairwise matrix was established and the relative weights for each criterion and suitability rank for classes are assigned as shown in Figure 5 and Table 1. The climate and rainfall criteria received the highest weights in both methods (three levels and two levels AHP). The values for each criterion were calculated and reclassified based on the 5 suitability classes and Equation (1) was applied to get the final suitability score for each site.  A pairwise matrix was established and the relative weights for each criterion and suitability rank for classes are assigned as shown in Figure 5 and Table 1. The climate and rainfall criteria received the highest weights in both methods (three levels and two levels AHP). The values for each criterion were calculated and reclassified based on the 5 suitability classes and equation 1 was applied to get the final suitability score for each site.  Table 5 shows measurements and scores for each criterion for the tabias receiving the highest (9 and 14) and lowest (10 and 15) suitability scores when AHP Method 1 was applied (before applying Equation 1).   Table 5 shows measurements and scores for each criterion for the tabias receiving the highest (9 and 14) and lowest (10 and 15) suitability scores when AHP Method 1 was applied (before applying Equation (1)). Figure 6A shows the overall suitability scores and the suitability score for each criterion based on Method 1 (three levels) after applying Equation (1). The highest overall score was 3.32 (medium suitability) for tabia 9, whereas the lowest score was 2.04 (low suitability) in tabia 10.

Test Results Sub-Catchment 1
Design criteria (structure dimensions, storage capacity and catchment area to cropping area) are playing a significant (negative) role in the overall RWH suitability for most of the tabias in Sub-catchment 1. These sites scored the lowest on design criteria, resulting in the low overall performance of these RWH sites. This result confirmed the observations of performance in the field.
A possible reason for the poor design is a lack of selection procedure for suitable RWH sites in combination, in this case, with structures built without a proper engineering design. Figure 6B shows the suitability scores for each criterion without multiplying by the weights. Table 5. The measurements and scores for each criterion (indicator) for the tabias receiving the highest (9 and 14) and lowest (10 and 15) suitability scores in Sub-catchment 1, when AHP Method 1 was applied (before applying Equation (1)).  Figure 6A shows the overall suitability scores and the suitability score for each criterion based on Method 1 (three levels) after applying Equation 1. The highest overall score was 3.32 (medium suitability) for tabia 9, whereas the lowest score was 2.04 (low suitability) in tabia 10. Design criteria (structure dimensions, storage capacity and catchment area to cropping area) are playing a significant (negative) role in the overall RWH suitability for most of the tabias in Sub-catchment 1. These sites scored the lowest on design criteria, resulting in the low overall performance of these RWH sites. This result confirmed the observations of performance in the field.
A possible reason for the poor design is a lack of selection procedure for suitable RWH sites in combination, in this case, with structures built without a proper engineering design. Figure 6B shows the suitability scores for each criterion without multiplying by the weights.
In Method two (two levels), the pairwise matrix was applied directly on the sub-criteria. Table 6 shows the overall suitability and the suitability for each criterion for the highest (9 and 14) and lowest (10 and 15) scoring tabias using this method. Once again, the design criteria of dimension and storage capacity had a significant negative impact on the difference between the high-scoring and low-scoring tabias. However with Method 2, CCR did not stand out as a differentiating factor, but reliability and cost did. In Method two (two levels), the pairwise matrix was applied directly on the sub-criteria. Table 6 shows the overall suitability and the suitability for each criterion for the highest (9 and 14) and lowest (10 and 15) scoring tabias using this method. Once again, the design criteria of dimension and storage capacity had a significant negative impact on the difference between the high-scoring and low-scoring tabias. However with Method 2, CCR did not stand out as a differentiating factor, but reliability and cost did. Table 6. The overall suitability and the suitability for each criterion for the highest (9 and 14) and lowest (10 and 15) scoring tabias in Sub-catchment 1, according to Method 2 and after applying Equation (1).

Test Results Sub-Catchments 2 and 3
The suitability scores for each criterion and overall from applying Method 1 (three levels) in Sub-catchments 2 and 3 are shown in Figure 7. The socio-economic criteria played a significant role in the assessment methodology here, especially for jessour in these sub-catchments (8-16 in Sub-catchment 2 and 1-8 in Sub-catchment 3) because of the high cost of implementing and maintaining the RWH compared with the relatively small area and low quantity of water retained behind the dykes. Moreover, these techniques are most common in this region especially in the mountain areas. They seem to be the most suitable techniques to mitigate flood hazard, additionally, the stakeholders consider them to be part of their heritage. Table 6. The overall suitability and the suitability for each criterion for the highest (9 and 14) and lowest (10 and 15) scoring tabias in Sub-catchment 1, according to Method 2 and after applying Equation (1).

Test Results Sub-Catchments 2 and 3
The suitability scores for each criterion and overall from applying Method 1 (three levels) in Sub-catchments 2 and 3 are shown in Figure 7. The socio-economic criteria played a significant role in the assessment methodology here, especially for jessour in these sub-catchments (8-16 in Sub-catchment 2 and 1-8 in Sub-catchment 3) because of the high cost of implementing and maintaining the RWH compared with the relatively small area and low quantity of water retained behind the dykes. Moreover, these techniques are most common in this region especially in the mountain areas. They seem to be the most suitable techniques to mitigate flood hazard, additionally, the stakeholders consider them to be part of their heritage.  Table 7 shows the individual criteria and overall suitability scores for the highest and lowest scoring sites in Sub-catchments 2 and 3 after applying method 2.  Table 7 shows the individual criteria and overall suitability scores for the highest and lowest scoring sites in Sub-catchments 2 and 3 after applying method 2. Catchment to cropping areas ratio (CCR) has a significant effect on overall suitability scores in Sub-catchment 2, whereas in Sub-catchment 3 there was not a difference in CCR between the high and low scoring structures. Moreover, slope played an important role in the overall scores in Sub-catchment 3 but not in Sub-catchment 2 ( Table 7).

Comparison of Methods 1 and 2
A comparison between the two methods of applying AHP (three and two levels structure) in our methodology is shown in Figure 8. Although the results are very similar, Method 2 gives a slightly higher score for the jessour in Sub-catchment 2 (jessour 10-16) and Sub-catchment 3 (jessour 1-8). Table 7. The individual criteria and overall suitability scores for the highest and lowest scoring sites in Sub-catchments 2 and 3 after applying Method 2.
Catchment to cropping areas ratio (CCR) has a significant effect on overall suitability scores in Sub-catchment 2, whereas in Sub-catchment 3 there was not a difference in CCR between the high and low scoring structures. Moreover, slope played an important role in the overall scores in Sub-catchment 3 but not in Sub-catchment 2 ( Table 7).

Comparison of Methods 1 and 2
A comparison between the two methods of applying AHP (three and two levels structure) in our methodology is shown in Figure 8. Although the results are very similar, Method 2 gives a slightly higher score for the jessour in Sub-catchment 2 (jessour 10-16) and Sub-catchment 3 (jessour 1-8).
The consistency of each matrix was calculated using the consistency ratio (cr). For the main criteria matrix in Method 1 cr was 2.9% and for the second method cr was 2.4%. The principles of AHP call for the cr to be smaller or equal to 10%, therefore the cr values were acceptable.
These results suggest that both methods are good and easy to adapt to different criteria, thus researchers can apply either of the two methods.  The consistency of each matrix was calculated using the consistency ratio (cr). For the main criteria matrix in Method 1 cr was 2.9% and for the second method cr was 2.4%.
The principles of AHP call for the cr to be smaller or equal to 10%, therefore the cr values were acceptable.
These results suggest that both methods are good and easy to adapt to different criteria, thus researchers can apply either of the two methods.

Results Validation with the Stakeholders
Based on our discussions with farmers and data collection from literature, we assessed the performance of existing RWH with the evaluation tool. Then, the preliminary results were checked with our field observations and discussed with local farmers and experts. For instance, the RWH sites which scored 2 or lower (low suitability) had been abandoned and or most of their trees were dead. Whereas the sites that scored around 3 (medium suitability) showed well-maintained structures with healthy trees.

Discussion
Fifty-eight RWH sites (44 tabias and 14 jessr) in three sub-catchments were assessed and evaluated on their technical and economic performance as well on social aspects. Using our methodology, 65% of the assessed sites scored around 3 (medium suitability), 31% of the RWH sites got scores of about 2 (low suitability), and only 4%, two sites, scored 4 (high suitability). These results very accurately represent the real performance of each site-both overall and at individual criteria level based on the comparison of our observations and discussion with local farmers and experts. This suggests that the methodology developed is a valid way to assess the performance of RWH structures.
The percentage of each score for each criterion in all sites was shown in (Figure 4). Rainfall had the same score (Score 3) in all sites because of there was no big difference in rainfall pattern nor amount in three sub-catchments. This means the rainfall indicator has no significant impact on overall suitability between sites in our case study, but it can be very important in the comparison between sites in the larger areas [21] with a significant difference in rainfall. Moreover, significantly low score percentages were obtained by the design criteria, drainage length and cost, which was Score 1. For example, drainage length scored 1 for 48% of all sites. That means the distance between watercourses and RWH structures is big and the score would have been higher if these structures were built closer to the watercourse. If the RWH structures were located much closer to the watercourses, the contribution of drainage length to the overall RWH suitability would have been higher for our case study. Therefore, drainage flow has a significant impact on the performance of the RWH, which is not always the case for other types of RWH such as ponds, terraces, etc.
It is interesting to note that although the weight for climate criteria was higher than that for site characteristics criteria, 30% and 26% respectively ( Figure 5),the latter received the highest scores in most of the sites in all three sub-catchments (Figures 6 and 7). This indicates that the sites are generally well selected for their purpose, and the site characteristics criteria had more impacts on the performance of RWH than other criteria such as climate, drainage and structure design. These results are similar to other studies, such as Al-Adamat [6] and Mbilinyi [16], who concluded that site characteristics are the most important criteria to be considered for design and implementation of RWH techniques.
Where RWH performance (suitability) was low, it was in most cases related to a shortcoming in the engineering design, lack of proper maintenance and the high cost of the water storage. The low performance of these RWH sites was confirmed by getting low scores of these criteria, as shown in Figure 4. The evaluation using our methodology clearly shows which criteria should be addressed to improve the performance of, for example, RWH structure design and storage capacity criteria. Due to the small storage area relative to the dyke size, the cost per cubic meter of water, especially in the jessour, was very expensive-such as jessour 10 and 15 in Sub-catchment 1. These results confirm that water harvesting structures with small storage capacity can ultimately be more expensive than large structures, as shown by Lasage, R., & Verburg, P.H. [4]. Therefore, if farmers can improve the dyke design and storage capacity area by following some basic engineering principles such as increasing storage area, constructing a regular spillway and providing periodic maintenance, they will be able to collect more water with less cost and keep the structure working for a longer period of time. Another example is the ratio between catchment size and cultivated area. Where this is not suitable, such as structures 11 and 21 in Sub-catchments 2 and 3, respectively, RWH structure performance can be improved by adapting the cultivated area to the effective area where the water is stored and adapting the crop type or cropping density (which determines the water requirement) to the amount of water stored.
In our methodology, two methods were applied (three levels and two levels of AHP hierarchy structure), and the results for both approaches were very similar. The consistency ratio for both methods was also similar and strong. Therefore, both methods are valid and provide reliable results. Both methods are simple to apply and easy to adapt the criteria in case of different RWH techniques and/or regions in order to cater to stakeholders' objectives. While either method can be used, it is recommended to apply Method 1 (three levels). In Method 1, the impact of any errors in scores (from expert opinion or calculations) will be reduced through the two-step calculation.
In most previous studies, the number of criteria are limited and are aimed primarily at the selection of suitable locations for RWH [1,14] and do not consider other factors or performance over time. In addition, many of those studies were mainly desktop studies using GIS and RS, without including stakeholders' objectives and constrains. Our study showed that socio-economic aspects play an important role in RWH suitability and performance. Thus, the inclusion of such criteria as occurs in our methodology is very important to the goal of meaningful information for improving current RWH effectiveness as well as planning for future structures.
A key precondition for the methodology was that it can be widely applied for different RWH techniques in different regions. In this regard, the structure of the methodology allows it to be easily adapted and applied to different RWH techniques and social-economic settings by simply changing the criteria selected. In addition, the case study showed that it is very possible to select criteria that are easy to assess and still provide accurate results without the need for complex analysis. This keeps the time investment and costs required within reasonable limits.
While Al-Adamat 2008 [6], Jabr 2005 [37] and Mbilinyi 2005 [33] showed that MCA provides a rational, objective and non-biased method for identifying suitable RWH sites, our study demonstrates that combining MCA and expert opinion in a consistent way allows assessment and evaluation of RWH techniques beyond simply site selection. Site conditions and RWH structure performance are likely to change over time, especially in light of predicted climate change. Therefore, a methodology such as ours, which allows evaluation of the performance of current and potential RWH projects, and identification of necessary improvements, is of great value.
An important consideration in the application of our methodology that warrants mention is the establishment of the scores/weighting for each criterion. As this depends on expert opinion [24,27], it is essential to use several experts and take into consideration their area of specialty when analyzing and using their inputs.

Conclusions
An evaluation and decision support methodology/tool was developed and tested for assessment of the overall performance of existing RWH and criteria affecting that performance. A single-objective AHP supported by GIS was put to the test in the Oum Zessar watershed of south-eastern Tunisia to assess the performance of 58 RWH structures (jessour/tabias) in three main sub-catchments. Engineering (Technical), biophysical and socio-economic criteria were determined, weighted and assessed in this study with input from experts and stakeholders. The main conclusions are: (a) The methodology provides an accurate evaluation of RWH performance when compared with the field investigations; (b) The methodology provides a good insight into where in the system improvements are needed for a better performance; (c) In the case study, most sites showed low suitability scores for the criteria structure design, drainage flow and cost, which resulted in a low score on the overall performance of RWH; (d) Site characteristics criteria (both overall and individual criterion) play a more important role in the overall suitability than other criteria; In addition, the methodology can be used to pre-evaluate potential new RWH projects, increasing the chances for good long-term performance. This case study application of our methodology confirmed that it is a highly flexible and applicable tool for the evaluation and improvement of RWH structures, and can employ many different, important and easy to access criteria and indicators in the assessment of different RWH techniques. The time and cost required in using this methodology are also low, making it accessible to the local RWH managers/communities.
To further validate the applicability of the methodology, it needs to be tested in different regions and with different RWH techniques. Moreover, the criteria related to socio-economic suitability/ performance (i.e., ownership, family size, etc.) deserve further investigation. These suggestions will increase the reliability and applicability of our methodology so that it can be used for assessing the performance of existing and new planned RWH structures in any region. This new, scientifically-based evaluation and decision support tool provides a basis on which designers and decision makers can build efficient RWH systems to meet the objectives and needs of the communities in water-scarce regions.