Moroccan southeast areas have limited water ressources, vulnerable to climate change and characterized by a significant spatio-temporal variability. In response, to ensure the availability of water for local comunity, it is advised to develop some alternatives that improve the local water resources management throughout these areas. Rainwater harvesting (RWH) has been used widely as an alternative technique towards water scaricity. However, taking into account socio-economic constraints for identifying ecological sites for RWH remains a complex task for water managers. The present study was conducted using WaTEM/SEDEM model, GIS techniques and Fuzzy Analytical Hierarchy Process (FAHP) method to identify important ecological RWH sites. For this purpose, several data sources were employed to generate needful thematic layers. The soil conservation service curve number (SCS-CN) method was utilized for preparing a yearly runoff potential map. Then, the thematic layers were weighted for generating RWH suitability map. The results show that yearly surface runoff ranges from 136 to 500 mm. Moreover, the spatial distribution map of soil erosion of WaTEM/SEDEM shows that Ziz upper watershed can be classed into four classes: (i) slight (very suitable), (ii) moderate (suitable), (iii) high (less suitable) and (iv) severe (not suitable). Aproximatly 76.1% of the study area falls within slight soil erosion class. The GIS tools were used for generating the important ecological RWH suitability map. The GIS’s intersect tool was used to eliminate insuitable soil erosion classes and the buffer tool was used for integrating the socio-economic factors including the distance from residential areas and agriculltural fields. Prediction accuracy of the resultant map of RWH suitability showed the value of the area under the curve (AUC) equal to 59.6% for FAHP method in this study. Ecological RWH map, mainly, introduced western areas and some central parts of Ziz upper watershed as suitable RWH areas. The present study demonstrated that coupling WaTEM/SEDEM model with FAHP method and GIS tools provide a valuable approach for identifying the ecological RWH sites in large semi-arid areas.