Enhanced groundwater availability through rainwater harvesting and managed aquifer recharge in arid regions

Climate change in desert areas and semi-arid watersheds may offer a promising solution for the water scarcity problem that Bedouins and local inhabitants face. This study investigated the integrated water resources management in arid and semi-arid regions using rainwater harvesting in combination with the managed aquifer recharge (RWH-MAR) technique. The study also used recharge wells and storage dams to achieve the sustainability of groundwater supplies in the context of climate change and management of the flow to the Gulf of Suez. Therefore, different return periods of 10, 25, 50, and 100 years were considered for the annual flood volume resulting from those watersheds. Moreover, hydrologic modeling was carried out for the El Qaa plain area, South Sinai, Egypt, using the Watershed Modeling System (WMS) and the groundwater modeling of SEAWAT code. Our findings show that for different scenarios of climate change based on return periods of 10, 25, 50, and 100 years, the aquifer potentiality reached 24.3 MCM (million cubic meters) per year, 28.8 MCM, 36.7 MCM, and 49.4 MCM compared to 21.7 MCM at 2014 with storage of groundwater ranges 11.8%, 32.1%, 69%, and 127.4%, respectively. These findings have significant implications for the system of RWH-MAR and groundwater sustainability in El Qaa Plain, South Sinai. The RWH-MAR proved to be an effective approach that can be applied in different water-stressed and arid regions to support freshwater resources for sustainable future development and food security, as well as protect communities from extreme flash flood events.


Introduction
Limited water resources are increasingly strained by growing demands for agriculture, drinking water, and industry, particularly in areas already coping with drought.This situation is further exacerbated by the impacts of climate change, which disrupts traditional watershed patterns and reduces water availability (Abd-Elaty et al., 2022a;Ibrahim et al. 2023).Flash floods (FF) are critical natural disasters worldwide, influencing people's lives, the economy, and infrastructure, especially in arid or semi-arid zones; watershed management and planning are essential for reducing flood damage (Akib Jabed et al. 2020;Jonkman 2005;Prama et al. 2020;Abd-Elaty et al. 2023).Also, FF is one of the most strenuous challenges we face in the twenty-first century (Jha et al. 2012).The Intergovernmental Panel on Climate Change (IPCC, 2014) published that climate change exacerbates many regions due to increasing the recurrence and magnitude of FF (Sharma and Ravindranath 2019;Abd-Elaty et al. 2024).
Rainwater harvesting and artificial recharge are among the most promising alternatives for supplying freshwater in the face of increasing water scarcity and escalating demand (Misra 2019).The best management of FF promotes fresh groundwater storage, preventing FF hazards and keeping the soil from drifting.Also, infiltration through the basin streams from the rainwater harvesting systems is the major source of artificial groundwater recharge and management for the present and future demands during the dry season (Mostafa et al. 2019;Zhang et al. 2020).Artificial recharge is a source and technique for sustainable groundwater resources in arid and semi-arid areas.The potential artificial recharge by flooding harvesting is better for mitigating the water shortage from over-pumping (Ibn Ali et al. 2017).Also, the aquifer storage and recovery (ASR) uses the same pumping wells to recharge the aquifers in wet and low-demand seasons from the precipitation and abstraction in dry seasons.This methodology has been widely used as an alternative solution for surface water storage in aquifers to reduce the construction cost of these structures and storage in surface reservoirs ( Jeong et al. 2018).ASR technique stores the additional water through deep injection in coastal aquifers.The stored water is then recovered using the same recharge wells to provide the community's water supply during the next dry or high-demand years (Alqahtani et al. 2021).Also, ASR describes intentional banking and treatment of water in aquifers (Dillon 2005).
The interactions between fresh groundwater and surface water are effective for water demand.This sharing between the two systems can increase efficiency and improve water shortages and environmental conditions.Increasing the abstraction well rates due to overpopulation and climate change are expanding the water demands and the abstraction from the aquifers.Moreover, fresh groundwater resources provide about 31% of the world's drinking water (Bear et al. 1999;El Shinawi et al. 2022).Eldeeb (2019) conducted a field study to investigate the interaction between surface water and groundwater due to the construction of New Dairut Group Regulators in Assuit, Egypt.The field measurements showed that a confined aquifer is underlying the study area, and the two systems have a hydraulic connection.Abd-Elaty and Zelenakova (2022b) used different techniques and strategies for coastal deep and shallow aquifers' freshwater sustainability in high-aridity regions using the hydraulic method.The results showed that the hydraulic method is positive and affects the mitigation of saline water intrusion in shallow and deep coastal aquifers.
The flood's runoff losses are compound infiltration and evaporation.The infiltration is estimated by the Horton Equation, Green-Ampt model, Phi index method, initial and constant loss, and the Soil Conservation Service (SCS) curve number (CN) method (SCS 1972).Smemoe et al. (2004) tested the CN and the Green-Ampt methods.The study indicated that the Green-Ampt method performed better than the CN method.Chahinian et al. (2005) compared the performance and consistency of four infiltration models.The results indicated that Morel-Seytoux's model was better than the Philip, Horton, and SCS methods.Also, the consistency of Horton's method was better than Philip's.The flood runoff estimation uses different methods, including kinematic wave, nonlinear kinematic wave routing, synthetic unit hydrographs, EPA SWMM/nonlinear reservoir, 2D diffuse wave, Wallingford, large catchment, SPRINT, and SWMM.Sherman (1932) introduced the concept of unit hydrograph (UH), which is the runoff for unit depth of excess rainfall based on a uniform-intensity storm.Morad (2000) used the HEC-1 code to estimate the potential of water resources in the Wadi Sudr basin, Sinai, Egypt.The results showed that the storm events and the generated runoff hydrographs were deduced to construct the interrelationship of the basin hydrologic characteristics.Al-Salamah et al. (2011) estimated its excessive use's impact on Saq aquifer depletion in Buraydah in the Kingdom of Saudi Arabia.The results concluded that the pumping from the aquifer would result in significant cones of depression if the existing excessive pumping rates prevailed.It has been recommended that the management scheme be adopted to protect groundwater resources.Elewa and Qaddah (2011) integrated enhanced Thematic Mapper Plus (ETM +) images, GIS, a watershed modeling system (WMS), and weighted spatial probability modeling (WSPM) to investigate the groundwater potential in Sinai, Egypt.The study showed that Sinai has moderate groundwater potentiality for about 33,120 km 2 , representing 52% of its total area.Ahmed (2012) applied gravity and magnetic data over the El Qaa plain by using 381 ground gravity stations and 800 lines of aeromagnetic data to assess the hydro-geologic settings, aquifer structures extension, geometry, and water storage in El Qaa plain.Aboelkhair et al. (2020) studied the integration of geospatial techniques for mapping groundwater potentialities in El Qaa Plain, Sinai, Egypt; the results showed that remote sensing data along with geospatial techniques could provide a powerful tool for groundwater probabilities in arid lands.They thus can be applied in regions with similar conditions, such as the Middle East countries.Yousif and Hussien (2020) increased the opportunities for groundwater potential in Sharm El Sheikh, South Sinai, Egypt.They mitigated the flooding risk, enhancing the environment under arid conditions through groundwater recharge.Elbarbary et al. (2021) applied the HEC-HMS package to forecast the runoff in Wadi El-Meleiha, Sinai, Egypt.The study showed that accuracy would be improved by accounting for the excess rainfall for the spatial and temporal variation in predicting direct runoff hydrographs.The current study uses the Watershed Modeling System (WMS) and variable density groundwater model using SEAWAT.It investigates the impact of rainwater harvesting recharge by flash flooding for sustainable groundwater supplies in the desert of El Qaa plain, Sinai, Egypt, due to climate changes for different recurrence intervals of 10, 25, 50, and 100 years.The study proposed rainwater harvesting using recharge wells and storage dams for the basin water resources management and decreasing the flow to the Gulf of Suez.
This study assesses the potential of using rainwater harvesting in combination with managed aquifer recharge (RWH-MAR) in the desert of El Qaa Plain, Sinai, Egypt.This technique holds promise for enhancing freshwater supplies and securing water resources for the sustainable development of arid coastal regions, especially in drought seasons.Also, the RWH-MAR prevents rainwater from evaporating or flowing into the sea unused, maximizing water capture and protecting the communities from the risk of flood hazards.

Hydraulic aquifer recharge
This study uses two techniques for rainwater harvesting: artificial recharge using recharge wells and natural recharge using a storage dam.The method is required to increase the aquifer freshwater storage to overcome over-pumping problems and climate change in arid and semi-arid regions.
Artificial recharge consists of a positive or pressure barrier; the coastal aquifers are artificially recharged to increase fresh groundwater's hydrostatic force and increase the piezometric heads to overcome the saltwater intrusion.Generally, this method minimizes flash flood risk, increases fresh groundwater heads and storage, relieves drawdown by overpumping, improves water quality, and suppresses the saline water body (Yomoto and Islam 1992;Tanim and Goharian 2021;).Ghazaw (2009) investigated how damages from runoff can be reduced and how extra runoff can be used to recharge.The study results will be helpful in runoff control, mitigation of damages from runoff, and appropriate storage and recharge of runoff in various areas of Saudi Arabia.
Natural recharge this technique feeds the aquifers with excess surface water by constructing dams and weirs to manage the runoff flowing directly to the sea and flood protection (Fig. 1a).This methodology aims to help infiltrate the retained water in the soil to increase the volume of fresh groundwater storage and help in SWI management.However, this solution is better for unconfined aquifers requiring high top layer permeability; depending on top layer properties, the aquifer recharging process could take a long time.The construction and maintenance costs of required structures for this method are very high.This method is unsuitable for confined and deep aquifers (Abbas et al. 2020).
The recharge well is used for low flood volume with a 1 m diameter and 100 m depth, while the storage dams are used for large flooding volume.The well recharge is suitable and operational in certain hydrogeological settings, for the aquifers do not get the natural recharge because of the confining layers of low permeability (Hernandez 2014); the recharge well is presented in Fig. 1b.The recharge dams for the groundwater occur through the base of the pond and downwards to the underlying aquifer.The water in a reservoir is released over several days/weeks so it can gradually infiltrate through the riverbed downstream.The most effective recharge dams have silt removal mechanisms to allow the release of silt to maintain the infiltration rate (Standen et al. 2020) (Fig. 1c).
The recharge type was selected based on flood volume and aquifer type.Using a storage recharge dam with high permeability in an unconfined aquifer is good.In contrast, the recharge wells are good in semi-confined and confined aquifers for low permeability.In this study area, the Wadi El-Aawag watershed has high permeability, and the aquifer is unconfined, so the surface runoff will easily infiltrate groundwater.

Wadi El-Aawag Watershed
Wadi El-Aawag watershed, one of the largest basins in southwestern Sinai, flows into the Gulf of Suez drainage system.It locates between latitudes 28° 15ʹ, 28° 45ʹN and longitudes 33° 20ʹ, 34° 00ʹ E. It occupies an area of about 1960 km 2 , locally named Sahl El Qaa (Fig. 2).This area was selected as a promising zone in the Sinai for future development and tourism activities, increasing the number of inhabitants and expanding land use (El-Rawy et al. 2020).This area of the Wadi El-Aawag watershed is bounded in the North and East South by mountainous hard rocks while in the South by the Red Seas and Gulf of Suez at the west.The main water supply for the city of El Tor depends on the surface and groundwater recharge by the Wadi El-Aawag Watershed.
The watershed's trunk channel generally extends about 58 km from the north and northeast to southwest.Topographic elevations varied between 1 m and 2566 m (Fig. 3a).This plan is a promising area for future land use planning and agricultural activities for shoreline extension to about 180 km with suitable soils of low, sandy, and gravelly surface, which differ significantly from the adjacent eastern high mountainous hard rocks of South Sinai.
The watershed varied in slope degree from a gentle slope (< 6) to a very steep slope in the eastern part of the watershed (67.7 m/m) (Fig. 3b).The watershed is covered by four land use /land cover units.The built area class is represented mainly by the area close to El Tur City at the Wadi outlet.Cultivated lands are occupied by the downstream area (El Qaa Plain area), which is represented by many agricultural farms irrigated by the water extracted from the El Qaa Plain aquifer.The bare lands occupy most of the watershed, while rangelands are observed on the eastern and western sides.The land use and land cover (LULC) is significant for mapping the variation in soil types and anthropogenic activities in the plain area, which is critical for runoff infiltration (Fig. 3c).The figure shows the five classes of LULC: basement rocks, old alluvial deposits, recent and wadi deposits, sedimentary rocks, urban & roads, and vegetation (Yousif et al. 2020).

Hydrological model
The hydrological frequency analysis was carried out using the HYFRAN-Plus software V2.1 in the Wadi El-Aawag watershed to estimate the return periods (RP).Also, the precipitation in El Tur and Saint Catherine stations from 1981 to 2022, where the minimum, average, maximum, and total annual rainfall was 0.38 mm, 7.18 mm, and 30.01 mm for El Tur station, while the Saint Catherine station reached 0.41 mm, 33.73 mm, and 8.92 mm, respectively, as presented in Fig. 4.
The recurrence intervals were estimated for different distributions.Log-Pearson type III developed a good fit for the

Wadi El-Aawag Watershed Geology and Hydrogeological setting
Wadi El-Aawag watershed geology predominantly comprises Precambrian and Cambrian rocks (Fig. 5).The Quaternary deposits are mainly from the El Qaa Plain, a coastal plain that is a promising aquifer for groundwater storage.
Wadi El-Aawag consists of twelve sub-watersheds from the north to South.The El Qaa Plain has a high mountainous range in the eastern part; these mountains play a vital role in the coastal plain development.Many drain of the watersheds plain (wadies) originated from these granitic mountainous massive and formed the Gulf of Suez drainage system in the study area.These wadies are buried within El Qaa Plain during the Quaternary times and become active during the rainy periods (Eltahan et al. 2021;Gaafar et al. 2018).These separate drainage systems are suitable for delivering mountain sediments to the El Qaa Plain and recharging groundwater during flooding.This plan attains the sixth order following the main geologic structures and rock contacts.
Figure 5a presents the El Qaa Quaternary aquifer system on the eastern coast of the Suez Gulf in the South of the study area.Also, it is bounded from the east by the high Precambrian crystalline basement rocks and from the west by the sedimentary carbonate rocks (Gabel Qabeliat).It is formed mainly of alluvial deposits and detrital materials detached from the surrounding basement rocks.It is characterized by high to moderate productivity.The Quaternary aquifer system in the El Qaa Plain is characterized by high storage capacity.It receives considerable amounts of recharge from the El-Aawag Watershed.This Quaternary aquifer in El Qaa Plain is the main source of fresh groundwater in South Sinai for El Tur and Sharm El Sheikh Cities.Lithologically, the aquifer consists of three layers, which are hydraulically connected.The aquifer lithological consists of sand, gravel, silt, and rock fragments detached from the basement rocks, varying between 700 m and 1000 m (Abbas et al. 2020), as shown in Fig. 5b.

El-Quad watershed modeling description
The Watershed Modeling System (WMS) is software for the El Qaa Plain to carry out hydrologic and hydraulic modeling of a watershed.The software modules include: The terrain module uses digital elevation models (DEM) data for floodplain delineation and mapping.This stage was applied for morphometric analysis of watersheds using DEM with 90 m; the study area resolution was obtained from the SRTM (Shuttle Radar Topography Mission data).Also, the geological and topographic maps were used to understand the different unit distributions for land use and soil type in the study area.
Drainage module is used for drainage basins and streams.The drainage network extraction for El Qaa Plain was performed through its sub-modules using the TOpographic PArameteriZation program (TOPAZ) after the watershed was delineated and subdivided into sub-watersheds.The drainage basin's characteristics, including areas, slopes, and maximum flow distance, were estimated.
The hydrologic modeling module uses Hydrologic Engineering Center-Hydrologic Modeling System (HEC-HMS) software to determine the rainfall-runoff relationships based on watershed characteristics by simulating the hydrological processes for basins and sub-basins to estimate the hydrographs for different watersheds.The HEC-1 can be used to develop the hydrological model by simulating a single watershed or a system of multiple hydrologically connected watersheds; it requires different The SCS curve number method is developed in the current simulation to estimate the losses.Also, the water quantities and inflow rate calculations were calculated using the SCS unit hydrograph approach.Furthermore, the hydraulic modeling module used the Hydrologic Engineering Center-River Analysis System (HEC-RAS) to design the streams and detention structures (El-Fakharany and Mansour 2021; Gülbaz et al. 2019).The mean annual rainfall of El Qaa Plain was used in performing the WMS 8.0 © from the meteorological stations' El Tor using different recurrence intervals.Also, the losses of the SCS-CN method used 88, the hydrograph approach was applied using SCS, UH, and the hydrograph was estimated.The flood volumes in the current study area of El Qaa Plain reached 89,892 m 3 at an average annual rainfall of 17.1 mm in 2014, which was used to calibrate the groundwater model.
Figure 6a presents the stream order for the El Aawag watershed, which has the sixth order.Table 1 and Fig. 6b show the characteristics of the El Aawag watershed and its sub-basins for the basin area, slope, average overflow, length, perimeter, shape factor, sinuosity factor, average basin elevation, maximum flow distance and slope, and maximum stream length and slope.

Surface runoff estimation
This method has Clark Unit Hydrograph, Snyder UH, Nash UH, and SCS UH.The watershed area was delineated using the SCS method as the following (USDA-TR55 1986; Gitika and Saikia 2014): The concentration time, defined as the runoff time required to travel from the most distant point to the basin's outlet point, is used in flood assessment.Kirpich equation (Chow et al. 1988) is used to compute the time of concentration described as follows: where S = potential maximum retention (mm), CN = curve number, P = rainfall for different return periods (cm), Q = the depth of direct runoff (mm), T c is the time of concentration in hours; L is the length of the overland flow (m); S is the average overland slope in (m/m).

El Qaa groundwater modeling description
El Qaa aquifer was simulated using SEAWAT code by variable density coupling.It is the modified MODFLOW (Langevin et al. 2003) and MT3DMS (Zheng and Wang 1999) into SEAWAT_V4 (Langevin 2009).The variable density groundwater flow equation is solved using the The groundwater flow model is solved using the finitedifferential method by dividing the system into a grid of cells.Each cell has a single node point where the groundwater heads are calculated (McDonard 1988).

Model design and geometry
The groundwater conceptual model was developed based on topography using DEM files, geology data, and hydrogeological studies.The Quaternary aquifer of El Qaa consists of gravel, sands, clays, and sandy clays.It is considered the main sediments composing the plain area (Sultan et al. 2009).The aquifer system is divided vertically into three hydrological units north of El Qaa Plain.The first aquifer is the upper unit by the lithological composition of gravels with a sand layer.It represents a limited area in the W. El-Aawag Watershed near El Tur City.Layer #2 is the middle unit by a soil of sand and gravel with clay lenses.It is considered the main groundwater unit with a 60 km length and 10 km width, occupying about 600 km 2 based on the hydrogeological maps of Sinai (JICA 1999).Its bottom is 20 m to over 100 m below sea level (masl).The deepest place is about 15 km northeast of El Tur City, representing the lower unit by sands and gravels with silt and clay.As previously discussed, these hydrological units are hydraulically connected (Fig. 3a).
The model has a rectangular dimension of 60,324 m in width and 65,380 m in length, representing an active aquifer dimension of 66 km (max.length) and 24 km (max.width).The El Qaa aquifer model includes many rows and columns, 200 and 200, respectively, using cell dimensions of 160,000 m 2 .

Boundary conditions and hydraulic parameters
Figure 7a presents the boundary conditions for the El Qaa aquifer, which are essential for a flow model construction, reflecting the relationship between the studied aquifer (4)  2).The aquifer hydraulic parameters were based on the data collected from JICA (1999), the Water Resource Research Institute (WRRI), and other studies using pumping tests for the production wells.The hydraulic conductivity of the eastern margin reached about 50 m day −1 with sandy gravel and gravelly sand.At the same time, the coastal area consists mainly of limestone and has a lower hydraulic conductivity of about 5 m day −1 .The transmissivity varies between 106 m 2 day −1 and 2150 m 2 day −1 , with an average of 768 m 2 day −1 .The specific storage ranges between 0.33-6.11,with an average value of 2.64.The effective porosity is 43.1% to Sultan et al. (2009).The three aquifers have hydraulic gradients of about 0.00395, 0.006, and 0.000775 (Elbeih 2015).

Abstraction and recharge
The abstraction rates for production wells in the El Qaa aquifer increase with increasing agricultural activities.All these wells depend basically on the El Qaa Quaternary aquifer system.Wells' yield ranges between 20 m 3 hr −1 and 15 m 3 hr −1 , with an average value of 56.7 m 3 hr −1 .The groundwater abstraction from the El Qaa aquifer reached 8550 m 3 day −1 , 13040 m 3 day −1 , 10270 m 3 day −1 , and 8820 m 3 day −1 in the years 1971, 1984, 1987, 1990, and 1992, respectively (El-Rawy et al. 2020), and reached 17,000 m 3 day −1 , 18,000 m 3 day −1 in 2004 and 2008, respectively (Elbeih 2015;Masoud 2015).The amount of recharge from the watersheds on the eastern side is about 5.9 × 10 6 m 3 year −1 (Masoud 2015).The recharge was estimated to be 16,840 m 3 day −1 in the year 2000.
Also, the aquifer received a modern annual recharge of 39 × 10 6 m 3 year −1 in 2020.The renewable water in the El Qaa aquifer ranges between 20 and 70% of the total aquifer storage (Yousif and Hussien 2020).The current model used an abstraction rate of 23,136 m 3 day −1 in 2014, while the recharge was 59,507 m 3 day −1 , as presented in Fig. 7b.

Model scenarios
The groundwater model was simulated for the seven scenarios for RP by 10, 25, 50, and 100 years applying Rainwater Harvesting combination with the Managed Aquifer Recharge (RWH-MAR) technique for the El Qaa plan using WMS software and the groundwater flow model of SEAWAT for rainfall events of 18.80 mm, 30.10 mm, 40.70 mm, and 53.40 mm respectively.

Results
The analysis and rainfall results, the watershed hydrology, and the groundwater flow model were presented in this study for El Qaa Plain.Also, the impact of climate change on the watershed basin in high-stress regions for different RP of 10, 25, 50, and 100 years were considered.

Estimation of watershed hydrograph for different recurrence periods
Table 3 presents 3.These results suggest that this freshwater water would have been lost without using rainwater harvesting in such desert regions as the El Aawag watershed.

Management of watershed for different sub-basins
Figure 8a presents intensity-duration-frequency (IDF) and the infiltration rates for the El Qaa plan.The initial and constant infiltration rates reached 1.23 mm hr −1 and 0.09 mm hr −1 , respectively.The rainwater harvesting plan for the El-Aawag watershed is a chance to increase the water infiltration in the El Qaa groundwater system    (Elewa et al. 2023).The hydrological model was carried out to estimate groundwater harvesting for nine subbasins.Table 4 shows sub-basin results for the RP of 10, 25, 50, and 100 years.Sub-basins # 1, #2, #3, #4, #5, #6, #7, and #8 were selected for storage recharge dam.In contrast, sub-basin # 9 was chosen to recharge the well due to the low volume of water harvesting at this subbasin #9 (Table 5).The selection of the dams was based on a recharge pond or reservoir because the top layer is Quaternary and consists of gravel, sands, clays, and sandy clays, which have high permeability with an average value of 23.9 m day −1 .Figure 8b presents the location of the protection structures using nine recharge dams and one recharge well.

El Qaa groundwater modeling results
Model calibration is essential for predicting groundwater levels under different scenarios.This step includes comparing the simulated groundwater system with the observed levels for a well.The current model calibration process has been performed through several trials by adjusting the hydraulic conductivity based on steady-state nonlinear conditions.
Figure 9a shows the difference in groundwater heads between the calculated heads (MODFLOW) and observed values for 22 observation wells by field investigation.The groundwater head residuals between the simulated wells ranged from -0.05 m to 1.19 m.The mean and absolute  residual reached -0.11 m and 0.75 m, while the standard error was 0.17 m.The root-mean-square error (RMSE) was 0.786 m with a normalization RMSE of 2.61%.The calibrated groundwater heads and the directions of groundwater flow for the El Qaa aquifer, South Sinai, Egypt, are presented in Fig. 9b and c, respectively.The general flow direction is observed from the east to the west in the southern half of the study area.In contrast, it is from northwest to east-west.The piezometric level varies between 5 m (a.s.l) in the southwestern part, near El Tur City, to 30 m (a.s.l) in the northern part of the El Qaa aquifer.The average velocity of groundwater reached 0.13 m day −1 .

El Qaa groundwater recharge using RWH-MAR
The calibrated groundwater model was tested using the system of RWH-MAR at RP of 10, 25, 50, and 100 years in the El Qaa aquifer to protect the city at different recurrence intervals.Figure 10 presents groundwater flow at each interval.It indicates that increasing the recurrence interval led to increased groundwater heads differences (GWHD), reaching 0.65 m, 9.5 m, 23 m, and 39 m at RP of 10, 25, 50, and 100 years, respectively.Thus, it is imperative not to ignore the impact of return periods when developing flash flood-prone areas.This thoughtful consideration is paramount, especially in safeguarding the freshwater resources vital to coastal regions and fortifying the critical infrastructures that underpin the resilience of these cities.By carefully accounting for return periods, urban planners and developers can implement sustainable measures that mitigate the immediate risks of flash floods and contribute to these vulnerable regions' long-term ecological and infrastructural sustainability.This holistic approach ensures that urban growth harmoniously aligns with the natural environment, fostering resilience and adaptability in unpredictable climatic events.

Discussion
The water harvesting model used changes in climatic variables to simulate future water availability.El Qaa Plain, in South Sinai, Egypt, is the largest groundwater aquifer in Sinai.It is considered the most suitable area for agricultural development and the main drinking water source for the major South Sinai cities (Mills and Shata 1989).Figure 11 shows the results of the El Aawag study area for  To mitigate the adverse impacts of escalating flooding, exacerbated by the effects of climate change, one viable strategy involves the integration of rooftop rainwater harvesting with managed aquifer recharge (Baptista et al. 2023).This approach aims to proactively address the escalating challenges of climate-induced flooding, presenting a sustainable and effective solution to safeguard water resources and mitigate the consequences of changing environmental conditions (Khanal et al. 2020).The of Rainwater Harvesting in combination with the Managed Aquifer Recharge (RWH-MAR) approach was adopted using storage dams and recharge wells in the study area to manage the flow to the flow to the Gulf of Suez and increase the aquifer potentiality for the desert region of El Aawag plain, Sinai, Egypt, applying eight recharge dams and one recharge well (Nachson et al. 2022).The aquifer potentiality reached 24.28 MCM, 28.82 MCM, 36.71 MCM, and 49.38 MCM, compared with 21.72 MCM in 2014 for the recurrence interval of 10, 25, 50, 100 years, respectively.
The study findings for using the system of RWH-MAR method agree with Missimer et al. (2015), who studied the aquifer recharge by the design of RWH-MAR in wadi dams using the groundwater model of MODFLOW; the results showed that about 80% of the water was lost by evaporation, while the using artificial recharge by wells improved the aquifer storage.The existing wells could store up to 1,000 m 3 day −1 under gravity-feed conditions and up to 3,900 m 3 day −1 with the shut-in of the well to produce a pressurized system.Gado and El-Agha (2020) studied the feasibility of rainwater harvesting for sustainable water management in urban areas of Egypt.The results indicated that the implementation of such RWH system has a significant impact on the regional water cycle, where the effective infiltration coefficient increased from 10% without RWH to 75% with RWH.Also, the volume of runoff decreased using RWH by around 82%. Arezoo et al. (2021) optimized the urban stormwater control strategies and assessed aquifer recharge through dry wells in an urban watershed; the results showed that the combining implementation of dry wells (DWs), bio-retention cells (BCs), and permeable pavement (PP) would increase infiltration by 19%, 15.6%, and 14%.Also, the runoff reduction rate in the presence of DWs would rise by 11.7%, 7.0%, and 6.1%.Recently, Szabó et al. (2023) studied the long-term potential and environmental impacts of rooftop rainwater harvesting coupled with shallow well infiltration, which is a local scale inexpensive solution that could contribute to easing water shortage in Danube-Tisza Interfluve (Hungary).Baptista et al. (2023) studied rooftop water harvesting for managed aquifer recharge and flood mitigation in João Pessoa, northeast Brazil.The results showed that the solutions provided mean annual estimates of aquifer recharge between 57 m 3 year −1 and 255 m 3 year −1 from 2004 to 2019, respectively.
The present study necessitates a thorough exploration through a comprehensive cost-benefit analysis and feasibility study concerning the construction of dams and recharge wells.This examination is crucial to accurately assess the cost-tobenefit ratio.Furthermore, addressing additional limitations, such as the assumption of uniform rainfall distribution over the basin, is imperative.It is recommended that the existing archival data be augmented by expanding the network of rainfall stations.This enhancement will contribute to refining the estimations of aquifer storage across each basin, thereby fortifying the robustness and precision of our findings.

Conclusions
Climate change in desert regions (arid and semi-arid regions) adversely impacts water supplies in regions with limited water resources.This study investigated the integration of For management plans and the integration between the surface and groundwater technique using rainwater harvesting, a comprehensive cost-benefit analysis and feasibility study are required to construct dams and recharge wells.Ultimately, the efficacy of the RWH-MAR system underscores its potential as a valuable tool.This effective approach holds considerable promise for policymakers in safeguarding urban areas against flash flooding, augmenting freshwater storage in arid regions, alleviating water scarcity, and promoting sustainable aquifer management.

Fig. 2
Fig. 2 Study site location and drainage network of W. Elawag watershed

Fig. 4
Fig. 4 Annual rainfall depths in the El Qaa aquifer

Fig. 5
Fig.5The study area: a simplified geological map of W. El-Aawag Watershed after Conoco (1987) and b conceptual (schematic cross section of El Qaa Plain after JICA(1999)

Fig. 6
Fig. 6 The El Aawag watershed: a basin's stream order and b its characteristics and sub-basins the WMS results for the El Aawag watershed considering 10, 25, 50, and 100 years of RP.The results of peak discharge reached 7.77 m 3 s −1 , 100.78 m 3 s −1 , 281.23 m 3 s −1 , and 587.10 m 3 s −1 for the different RP compared with 2.49 m 3 s −1 in 2014.In contrast, the flood annual volume reached 328,428 m 3 , 4,872,530 m 3 , 12,760,832 m 3 , and 25,432,076 m 3 at RP of 10, 25, 50, and 100 years, respectively, compared with 89,892 m 3 in 2014, as shown in Table

Fig. 7
Fig. 7 El Qaa Quaternary aquifer: a model geometry and boundary conditions and b average groundwater extraction through the years

Fig. 8
Fig. 8 El Aawag basin: a the location of storage dams and recharge well used for the sub-basins watershed, and b rainfall intensity and infiltration rates

Fig. 9
Fig. 9 Groundwater model for a observed and calculated heads, b flow heads, and c flow directions

Table 1
Characteristics of El Aawag watershed

Table 2
Calibrated hydraulic parameters and boundary conditions used in the El Qaa groundwater model

Table 3
Runoff discharge for W.
infiltration, water harvesting recharge, and aquifer potentiality for each return period.The infiltration of rainfall through the aquifer reached 23.95MCM, 23.95 MCM, 23.95 MCM, and 23.95 MCM  at RP of 10, 25, 50, and 100 years, respectively, compared  with 21.72 MCM in 2014.The hydrology losses and the aquifer recharge by the natural infiltration by 23.95 MCM before the runoff and the flash flood water go to the Gulf of Suez.Moreover, the system of RWH-MAR reached 0.33 MCM, 4.87 MCM, 12.76 MCM, and 25.43 MCM compared with 0.00 MCM in 2014, while the aquifer potentiality reached 24.28MCM, 28.82 MCM, 36.71MCM, and  49.38 MCM compared with 21.72 MCM at 2014.This represents a percentage change of the aquifer potentiality by + 11.78%, + 32.70%, + 69.02%, and + 127.36% at RP of 10, 25, 50, and 100 years, respectively, compared with 2014.The aquifer potentiality variation (P-P 2014 ) / P 2014 was calculated, where P 2014 is the aquifer potentiality in 2014, and P is the aquifer potentiality for the given scenario.
Infiltration Harvesting recharge Aquifer potential Fig.11El Qaa aquifer potential for different recurrence intervals Rainwater Harvesting combination with the Managed Aquifer Recharge (RWH-MAR) technique for the sustainability of groundwater supplies in the aquifer of El Qaa Plain using the Watershed Modeling System (WMS) and SEAWAT software under the impact of climate change.The groundwater is recharged from precipitation and runoff from the adjacent watershed of El Awag.Results indicated that the annual rainfall reached 18.8 mm year −1 , 30.1 mm year −1 , 40.7 mm year −1 , and 53.40 mm year −1 at return periods of 10, 25, 50, and 100 years, respectively.Also, the annual flood volume reached 328,428 m 3 , 4,872,530 m 3 , 12,760,832 m 3 , and 25,432,076 m 3 , respectively.The integrated water resources management was achieved in El Qaa plain by recharging the aquifer.This was done by proposing nine storage dams and one recharge well based on the RWH-MAR quantity and soil infiltration capacity of the sub-basins system.The rainwater for the Gulf of Suez was managed.The annual potentiality of groundwater resources in the aquifer reached 24.28MCM, 28.82 MCM, 36.71MCM, and  49.38 MCM for return periods of 10, 25, 50, and 100 years,  respectively, compared with 21.72 MCM in 2014.