Geochemical modeling and isotopic approach for delineating water resources evolution in El Fayoum depression, Egypt

In this study, an attempt was made to understand the hydrogeochemical processes controlling water resources (surface and groundwater) evolution in El Fayoum depression, Egypt. Hydrochemical, inverse geochemical modeling integrated with environmental isotopes was applied. Forty-three surface and groundwater samples were collected and analyzed for major chemistry (cations, anions) and environmental stable isotopes (δ18O, δD). The ionic sequence based on relative molar proportions for cations was Na+ > Ca2+ > Mg2+ > K+ and Na+ > Mg2+ > Ca2+ > K+ for surface and groundwater, respectively and for anions was HCO3− > Cl− > SO42− for all water samples. δ18O and δD relationship plot on a typical evaporation line, enhance potential evaporation of recharging water prior to infiltration. Saturation index estimation revealed that water resources were affected by carbonate and evaporites minerals dissolution and supersaturated with Goethite, Gibbsite and Barite minerals. Inverse geochemical modeling using PHREEQC was used to identify water resources evolution in the study area through two flow paths A, B. The inverse model was constrained so the primary mineral phase including calcite, Aragonite, Anhydrite and Gypsum with Carbon dioxide (gas) are constrained to dissolve until they reach saturation, Whereas H2(g), hydrogen sulphide (gas) and O2(g) are constrained to out gas in all models for flow path (A). For flow path (B), Dolomite, Anhydrite and Gypsum are tend to dissolve with carbon dioxide (gas), H2(g) and O2(g) gases tend to outgassed. Finally, two processes controlling water resource evolution; carbonate and evaporites dissolution and development of cation exchange process through the aquifer system.


Introduction
El Fayoum area is one of most important agricultural lands of Egypt for a long time due to good soil cover and high crop productivity, In addition to it has a beauty nature represented by natural springs, waterfalls and the great natural lake (Qaroun Lake). El Fayoum's water enters through two main gated intakes, one being the natural channel Bahr Youssef, the other Bahr Hassan Wasef a canal dug in 1905 (Wolters et al. 1987). There are two main drains; The Bats in the east and the Wadi in the west. All drainage water discharged into Qaroun Lake. Population growth and increase of human activities cause several problems in El Fayoum area such as water logging and soil salinization; these problems have already emerged after High Dam construction. Mismanagement of agricultural fertilizers and pesticides, reuse of drainage water for irrigation, wastewater disposal and inadequate land use planning had a major impact on groundwater degradation and soil quality in the studied area. Qaroun Lake, a closed saline lake acts as a natural discharging system for agricultural drainage water in El Fayoum depression. Previous studies focused mainly on this lake (Mansour andSidky 2002, 2003;Gupta and Abd El-Hamid 2003) to investigate the influence of those waters on the lake water quality. As the groundwater resource in El Fayoum depression is a Quaternary aquifer that shows increasing signs of quality deterioration (Abdel Kawy and Belal 2012;Gad and El-Zeiny 2016;Gad and Hattab 2019) attributed to irrigation return and anthropogenic activities. The quality of water is a function of its physiochemical parameters that reflect groundwater characteristics and factors that controlling the quality of water (Hussein et al. 2017;Mgbenu and Egbueri 105 Page 2 of 15 2019). Chemical composition of groundwater is controlled by many factors include composition of precipitation, geological structure, mineralogy of aquifers, and geochemical processes within the aquifer. The interaction of all factors leads to various water facies. Usually, major ion studies are used to define hydrochemical facies of waters and the spatial variability that can provide insight into aquifer heterogeneity and connectivity (Murray, 1996;Rosen and Jones 1998). A number of geochemical models have been developed to investigate geochemical reactions occurring along flow paths within groundwater systems. Geochemical modeling of aqueous systems is based on the assumption of chemical equilibrium among the dissolved species (Plummer 1992). Hydrogeochemical processes including dissolution, precipitation, ion-exchange, sorption and desorption together with the residence time occurring along the flow path, control the variation in chemical composition of groundwater (Apodaca et al. 2002), that can be modeled by inverse geochemical models. Inverse geochemical modeling attempts to account for the chemical changes that occur as water evolves along a flow path. Assuming two water analyses, that represent starting and final water composition along the flow path within a groundwater system, inverse modeling is used to calculate the moles of minerals and gases that must enter or leave solution to account for the differences in composition. The use of geochemical models is increasing in addressing groundwater quality problems involving geochemistry (Slimani et al. 2015;Berihu et al. 2017;Liu et al. 2020). Environmental stable isotopes (δO 18 , δD) modified by meteoric processes provide relevant and useful tracers to identify flow regimes and sources of groundwater recharge (Leontiadias et al. 1988). Frequently, groundwater retains its stable isotopic signatures unless diluted or mixed with waters of different isotopic compositions (Fontes 1980;Gat 1981). Therefore, water from different sources or those exposed to different processes such as evaporation and / or mixing often acquire identifiable isotopic composition which can serve as natural tracers. In this paper, physical, hydrogeochemical information integrated with isotopic data were incorporated into the inverse geochemical modeling to delineate water evolution in El Fayoum depression. Hence, the objectives of this study were: (1) identify various factors affecting water (surface and groundwater) hydrogeochemical processes to understand and distinguish between the water-rock interaction and anthropogenic influences (2) characterize and assess recharge areas and the origin of surface and groundwater in the study area (3) Understand the geochemical evolution of water resources as an essential part in water management and remediation processes with a better knowledge of how the groundwater system functions.

Geology and hydrogeology
El Fayoum area is a green, fertile, near-circular depression and it has a unique physical and environmental setting. El Fayoum area occupies one of the spectacular depressions in Western Desert of Egypt, about 90 km southwest of Cairo between latitudes 29° 00′ and 29° 40′ N and longitudes 30° 20′ and 31° 05′ E (Fig. 1b). It is surrounded by desert terrain except for the southeastern side, where it is connected to the Nile Valley by the Bahr Youssef canal which derives its water from Al Ibrahimya Canal. The depression is limited at the north-western side by a shallow brackish lake, called Birket Qarun (or Lake Qarun), with an area of about 200 km 2 and forming the deepest part at about 45 m below the mean sea level (msl). The total surface area of El Fayoum depression is approximately 6000 km 2 . About 1800 km 2 of the depression is filled with rich alluvial sediments and more than 1580 km 2 is cultivated. This area represents more than 5% of the total cultivated area of Egypt. It is an area of intensive agriculture, reliant solely on irrigation from the Nile River. El Fayoum depression has a special irrigation system due to the nature of its land which slopes downward from south to north. The wastewater flows from south to north towards the slope where it pours into Qarun and Wadi Al Rayan Lakes, the main sources of water drainage for the agricultural lands in the depression (UNDP 2003). Population of El Fayoum is about 2 million and the density of population is about 1000 people per km 2 . The overall climate of this area is characterized by hot, long and dry summers and warm, short winters with scarce precipitation. Also, great temperature differences between summer and winter and between day and night characterize this belt. The temperature ranges between 46.7 and 35.1 °C in summer months and from 20 to 21.5 °C during winter months. The average annual rainfall is around 10.3 mm which is relatively low. The total annual evaporation intensity reaches 2296 mm/ year, and the annual mean of relative humidity is 51.6%.
El Fayoum basin was formed initially during the Jurassic Period, probably by wind erosion in the desert, and its current shape created from tectonic subsidence that terminated in the late Eocene Epoch. The basin then subsided relative to the Nile, allowing the river to break through in flood and deposit fertile alluvial sediments. Rock units that exposed in El-Fayoum region range from Eocene to Quaternary Ages as shown in Fig. 2c. In the study area, the subsurface stratigraphic column is capped by the Pleistocene deposits that are widely distributed over the entire area of El Fayoum area. These deposits are mainly of fluvio-lacustrine origin. The Quaternary lacustrine deposits dominate in the south of Qaroun Lake, which consist of clay stone, gypsum and calcareous materials intercalated with ferruginous sandy silt (Redwan and Elhaddad 2017). The Quaternary sediments overlay the Middle Eocene deposits, which are composed of limestone, marl, sand, and gypsiferous shale that outcrop around the depression (Metwaly et al. 2010). The Middle Eocene deposits overlay the Oligocene deposits, which are composed of fluvio-marine variegated sands and sandstone, with alternating shale, marls, and calcareous beds. The Lower Miocene sediments are composed of continental sandstones and conglomerates, while sediments of Pliocene age are composed of calcareous sandstones and Limestones. Several distortion lines, such as faults, breaks, and some folds striking N-S and NW-SE, that affect El Fayoum Depression are dominant at the peripheries of the depression (El- Sayed et al. 2015). Based on a detailed analysis of the available lithological data, geological map and lithological succession were constructed (Fig. 2b, c).
The surface water hydrology system is represented by the Nile River, irrigation canals, and drains. The area gets irrigation water from the River Nile and the main irrigation branches Bahr Youssef and Bahr Hassan Wassef canals to the depression. Generally, three groundwater aquifers were distinguished in El Fayoum Depression ( Mohamed et al. 2015) namely, the Quaternary aquifer, Eocene aquifer, and Nubian sandstone aquifer. Quaternary aquifer is a perched water table aquifer that is maintained by infiltration from canals, drains, and irrigated fields. Many of the springs in El Fayoum depression have their sources from the Quaternary aquifer. Beneath this aquifer, the Eocene aquifer, that is composed of limestone and marls, has high salinity and is therefore unsuitable for potable or irrigation use. The Nubian aquifer is composed of sandstones and shales of Cretaceous age and it was very deep aquifer (> 1000 m), not exposed in the area of study but reported in the subsurface at Wadi El Rayan and in numbers of oil wells drilled in some localities close to the study area. The Nubian formation is generally characterized by relatively high hydrostatic pressure (Himida and Abdel Baki 1983). The Quaternary aquifer has its maximum thickness (about 60 m) at the center of the depression while, the minimum thickness is less than 1 m along the depression peripheries due to local faulting (Ahmed, 2012). The groundwater depths varied from few centimeters to 9.1 m, with groundwater level ranges from + 20 m (amsl) in the area near Nile River to-50 m (bmsl) near Qarun Lake. In addition, the ground-water flows from SE to N and NW in concordance with the general slope of the depression topography in Fig. 1c that shows the land slopes from Lahun to Lake Qarun which occupies the deepest part in El Fayoum depression (Gad and El Hattab 2019). The Quaternary aquifer is thought to be recharged by the seepage from surface irrigation canals, drains, and Nile River. In addition percolation from the underlying fractured limestone of the Eocene aquifer through hydraulic connection (Doerfliger et al. 1999). This study focused on the shallow water-bearing formations of the Quaternary aquifer that directly related to the objective of this investigation.

sample collection and analysis
A set of representative surface water samples (19 drainage water and 13 canals) and 12 groundwater samples were analyzed for physical measurements including electrical conductivity (EC), pH, temperature, total dissolved solids (TDS), and dissolved oxygen (DO) using portable probes. Chemical parameters comprising major cations (Ca, Mg, Na, K) and anions (Cl, SO 4 , HCO 3 ) concentrations were also analyzed. All the wells located in El Fayoum depression had been in use during the time of sampling (Fig. 1b). The depths of the selected wells ranged from 8 to 13 m below the surface. The samples were collected in 1-L narrow neck pre-washed polyethylene bottles. Analysis of the water samples was carried out following the methods described in APHA (1995). Total hardness (TH) as CaCO 3 and Ca 2+ were analyzed using standard EDTA. Mg 2+ was calculated by taking the differential value between TH and Ca 2+ concentrations. Na + and K + were measured using a flame photometer. Total alkalinity and CaCO 3 , CO 3 2− , and HCO 3 − were estimated by titrating with HCl. Cl − was determined by standard Hg (NO 3 ) 2 titration. SO 4 2− and NO 3 − were analyzed using  Beadnell (1905) and Said (1962). c Lithological succession with the location of proposed shallow aquifer at the center of the depression. Redrawn and modified after Metwaly et al.( 2010) UV/visible spectrophotometer. All parameters are expressed in milligrams per liter and milliequivalents per liter. Data quality was assessed using the charge balance between the difference of cations and anions (expressed in meq/l) divided by their summation according to the following equation: With an acceptable range of ± 5 (Hem 1991) that confirm the water quality assessment. Environmental stable isotope (δ 18 O and δD in ‰) samples were collected according to the procedure described in Clark and Fritz (1997). The analysis was done using Picarro L2120-I Isotopic Water Liquid and Vapour Analyzer. Both results were expressed in δ‰ (IAEA 1995) with analytical error of ± 0.1 and ± 1‰ for δ 18 O and δD respectively. The stable isotope data are reported in Vienna-Standard Mean Ocean Water (V-SMOW) in ‰ using δ notation according to the equation: R sample and R standard are the measured isotopic ratios ( 18 O/ 16 O) and ( 2 H/ 1 H) of sample and standard material.
Ground surface elevation and the location of the sampling sites were recorded using a global positioning system (GPS) manufactured by GARMIN. This was supported by a topographic sheet made available by the Egyptian Survey Department. Measured and estimated groundwater variables were analyzed within ArcGIS 10.2.2 package. All the chemical and Environmental isotope analyses were carried out in the Central Laboratory of Stable Isotope Hydrology at Egyptian Nuclear and Radiological Regulatory Authority (ENRRA).

Geochemical modeling
The saturation index is a vital geochemical parameter in the fields of hydrogeology and geochemistry, often useful for identifying the existence of some common minerals in the groundwater system (Deutsch 1997). In this present study, saturation indices (SIs) was calculated in terms of the following equation (Lloyd and Heathcote 1985): where IAP is the relevant ion activity product, which can be calculated by multiplying the ion activity coefficient γ i and the composition concentration m i , and K s (T) is the equilibrium constant of the reaction using the measured sample temperature. When the groundwater is saturated with some minerals, SI equals zero; positive values of SI represent oversaturation, and negative values show undersaturation (Appelo and Postma 1994;Drever 1997). Inverse geochemical modeling with the thermodynamic program PHREEQC (Parkhurst and Appelo 1999) was used to calculate the (1) Σ (Cations − Anions)∕Σ(Cations + Anions) × 100, (3) SI = log(IAP∕K s (T), geochemical processes responsible for groundwater evolution or changes in quality. To deduce geochemical reactions that accounts for the change in chemical composition of groundwater along the flow path, both aqueous species and saturation indices were calculated. The flow paths were selected on the basis of hydrogeology, hydraulic head variation within the flow system and geographic location in the study area. Minerals used in inverse geochemical modeling are limited to those known to occur in the study area. The mole transfer of minerals and gases reported in mmol kg −1 of water. Positive values of phase mole transfers indicate dissolution of minerals or gases, and negative values represent precipitation or outgassing. The pH values of surface water samples (canals, drains) ranged from 6.16 to 6.7 with mean value of 6.38 and 5.91-7.12 respectively. For groundwater samples the pH values ranged from 5.65 to 6.51 with mean value of 6.09 the data were reflective of slightly acidic to slightly neutral conditions due to that the selected wells were shallow with depth to water ranged from 8 to 13 m, the dissolved carbonates is predominately in HCO 3 − form. In surface water, EC and TDS values exceeded those of groundwater may be due to anthropogenic pollution inputs. The EC of surface water (canals and drains) varied from 606 to 2636 μs/cm, with a mean value of 1108 μs/cm and 565 to 3756 μs/cm with a mean value of 1660 μs/cm respectively. In groundwater samples EC values ranged from 162 to 4426 μs/cm, with a mean value of 1582 μs/cm. TDS values of groundwater samples ranged from 104 to 2832 mg/l exceeding (WHO 2017) standard limits for 50% of collected samples. TDS values ranged from 388 to 1687 mg/l and 361-2431 mg/l in canals and drains samples exceeding (WHO 2017) limits for drinking purposes with 23% and 42% respectively according to (Table 2). According to Konikow and Reilly (1999) and Rhoades et al. (1992) salinity classifications about 42%, 38% of groundwater samples and canal water were slightly saline, while about 42%, 54% of groundwater and canal water were classified as fresh respectively. The relative abundance of major cations in shallow groundwater is Na + > Mg 2+ > Ca 2+ > K + (on a molar basis) and HCO 3 − > Cl − > SO 4 2− for anions. For surface water samples (canals and drains) the predominant cations are Na + > > Ca 2+ > Mg 2+ > K + (on molar basis) and HCO 3 − > Cl − > SO 4 2− for anions. The maximum Ca 2+ and Mg 2+ concentrations of 144 mg/l and 168 mg/l respectively in groundwater samples are however, higher than WHO (2017) standards limits by 50% and 58% respectively. The prevailing common lithology is the old lacustrine deposits which are present within the depression and consist mainly of clay, fine to medium coarse sand with considerable amounts of calcareous and gypsiferous materials. The mean sodium concentrations is 159 mg/l is within the standard limits but potassium concentration is 17.6 mg/l in groundwater samples that exceeds WHO (2004) standards while, for irrigation canals the mean sodium and potassium concentrations are within WHO standards limits While, for drains Na, K exceeds the international standard as illustrated in Table 2. The presence of HCO 3 − ions in shallow groundwater samples is derived from carbon dioxide in the atmosphere, soils and by dissolution of carbonate rocks. Bicarbonate ion represents the first dominant anion in the study area. Most Cl − in the groundwater related to lacustrine deposits found in quaternary aquifer of El Fayoum depression. The most extensive and important occurrences of sulphate ions in the investigated area are sedimentary rocks such as gypsum (CaSO 4 ·2H 2 O) and anhydrite (CaSO 4 ) that give the chance for active dissolution, leaching and ion -exchange processes. This dominates in the central and western parts of the depression. Further addition of sulphate to groundwater comes from breakdown of organic matter in the soil,  and from addition of leachable sulphate in fertilizers of intensively cultivated area in the depression. Hypothetical salt combination could classify the hydrochemical data into three groups as follow:
As indicated by the spatial distribution maps of salinity (Fig. 3), The TDS increases from the peripheries that dissected by drains towards the central part of the study area, where the seepage from drains leads to groundwater salinization. Higher TDS values may be related to water-rock interaction, upward recharge from the underlying fracture limestone of Eocene aquifer, irrigation return flow and seepage from the surrounded drains (El-Sheikh 2004). Low TDS values may be attributed to refreshment from irrigation canals.

Processes controlling major solutes distribution
The rock dominance of the major ion chemistry in the study area provides an insight of chemical weathering in the aquifer, since weathering of different parent rocks Fig. 3 Total dissolved solids distribution maps for a groundwater, b irrigation canals, c drainage water (e.g., carbonates, silicates and evaporites) yield different combinations of dissolved cations and anions to solution. The Gibbs diagram is widely used to study the relationship between water chemistry and aquifer lithology (Gibbs 1970). The functional sources of the dissolved chemical constituents can be divided into three distinct fields in this diagram, namely rainfall dominance, evaporation dominance and rock-weathering dominance areas. The ratios of Cl/ (Cl + HCO 3 ) and (Na + K)/(Na + K + Ca) as a function of TDS were drawn in Fig. 4. The majority of the samples fell in evaporation-dominated field demonstrating the  importance of irrigation return flow on collected surface and groundwater samples chemistry. According to (Table 3) for calculated ionic ratios and ( Fig  5a) the Na + /Cl − showed clear enrichment of Na + in the surface (canals and drains) and groundwater samples. The Na + / Cl − ratios in the study area varied from 1 to 3.3 with an average value of 1.6 in groundwater and for surface water, it varied from 1.2 to 2.2 in canals and 1.2 to 3.4 in drains. In general, the Na + / Cl − molar ratio should be approximately equal to 1, if the halite dissolution is the only source. The majority of the samples (Fig  5a) having a higher Na + /Cl − ratio suggested another possible source of Na + may be ion-exchange or silicate weathering be responsible for the exceeding of Na + concentration in collected water samples. The Ca 2+ /Mg 2+ ratio suggests that dissolution of calcite and dolomite in the aquifer. Results of surface and groundwater samples ratios (Fig 5b) ranged between 0.49 and   (Fig.5d) shows that all surface and groundwater samples are distributed close to or below 1:1 line due to an excess of (SO 4 2− + HCO 3 − ) over (Ca 2+ + Mg 2+ ) an indication of that the dissolution of calcite, dolomite, anhydrite and gypsum are the dominant reactions and ion exchange tends to be predominant process in the system (Mclean and Jankowski 2000).
Chadha's diagram can be used to study the various hydrochemical processes responsible for water system evolution such as base cation exchange, mixing of natural waters and saline water (Chadha 1999). Results of hydrochemical analysis plotted on this diagram (Fig 5e) show that some groundwater samples can be related to field 1of (Ca-Mg-HCO 3 water type) and drains are related to field 2 of (Ca-Mg-SO 4 /Cl mixed water) and the rest of samples related to field 3 (Na-Cl water type). To check the hypothesis of cation exchange is one of geochemical process that may control water chemistry, The relation of (Ca 2+ + Mg 2+ ) − (SO 4 2− + HCO 3 − ) vs. (Na + -Cl − ) meq/l was illustrated in Fig 5f can better indicate the involvement of cation exchange between groundwater and its host environment. Negative values of (Ca 2+ + Mg 2+ ) − (SO 4 2− +HCO 3 − ) are usually used to demonstrate the involvement of cation exchange, inducing the ions absorbed on the surface of fine-grained materials of the aquifer to be replaced by the ions in the solutions. The ratio of this trend in surface and groundwater samples was linear with a slope of − 1.33 and the correlation coefficient was 0.97. As the cation exchange between Na + and Ca 2+ or Mg 2+ will lower the concentration of Ca 2+ dissolved in water, this was illustrated in the relation between Ca 2+ and SO 4 2− in Fig 5c.

Geochemical modeling
In order to evaluate possible hydrochemical reactions along the flow paths, inverse mass-balance models were developed using the PHREEQC program (Parkhurst and Appelo 1999). Three geochemical processes may contribute to generation of solutes in groundwater including evaporation, carbonate dissolution /precipitation and silicate weathering (Garrels and MacKenzie 1971). Statistical summary of thermodynamic speciation calculations for surface and groundwater samples using PHREEQC program were illustrated in Table 4. The results indicated that the surface (canals, drains) and groundwater samples of the Quaternary aquifer Two flow paths were modeled based on Hydrochemical data as illustrated in Fig. 1b, the flow path (A) at which the less mineralized (recharge water) was canal water (sample 20) which has TDS value of 388 mg/l and the moderate mineralized (transition water) one was drainage water (sample 6) which has TDS value of 2431 mg/l. The flow path (B) at the moderate mineralized drainage water (sample 6) also modeled with the more mineralized (Final water) groundwater sample (sample 10) which has TDS value of 2833 mg/l. The inverse modeling in PHREEQC takes into account uncertainty limits that are constrained to satisfy the mole balance for each element and valence state as well as the charge balance for each solution within the simulation. The simulations were constrained within the specified uncertainty limits (0.025 or 0.05) as default for all model runs. Potential phases in the inverse modeling were selected based on (precipitation/dissolution) results from hydrochemical data and saturation index calculations for the modeled samples (Table 5). The results of inverse geochemical modeling were illustrated in Table 5 as shown below: 26 and 14 models were resulted for flow paths A, B respectively with selection done based on the lowest sum of residuals and uncertainty limits. The inverse model was constrained so the primary mineral phase including calcite, aragonite, anhydrite and gypsum with carbon dioxide (gas) are constrained to dissolve until they reach saturation, Whereas H 2(g) , hydrogen sulphide (gas) and O 2(g) are constrained to out gassed in all models for flow path (A). For flow path (B), Dolomite, Anhydrite and Gypsum are tend to dissolve with Carbon dioxide (gas), H 2(g) and O 2(g) gases tend to precipitate. The Quaternary deposits cover wide stretches of the investigated area are distinguished into Holocene deposits and Pleistocene deposits as mentioned in El-Sheikh (2004).The old lacustrine deposits (45 m thick) formed the Pleistocene deposits composed mainly of clay fine to medium coarse sand with considerable amounts of calcareous and gypsiferous materials. As carbonate minerals (calcite, dolomite) dissolution for both flow paths A, B contribute to Ca 2+ , Mg 2+ and HCO 3 − and Gypsum dissolution leads an increase in Ca 2+ and SO 4 2− to the system. Additional increase in Ca 2+ leads to dissolution of more calcite that increase CO 2 , which leads to slightly lower pH as illustrated in Table 1. The out gassing of hydrogen sulphide is occurring in flow path A, which an indication of SO 4 2− reduction. Sulphate reduction occurs extensively in natural systems that involves the consumption of substantial amounts of hydrogen ions and production of HS − at certain pH. As drainage water contains organic carbon or H 2 that is oxidized while sulphate is reduced (Canfield 2001) according to the following equations: The hydrogen ions are consumed and bicarbonate is produced, which leads to an increase in PH as shown in Table 1 for drainage water samples. Field measured ORP (oxidation − reduction potential) for drainage water have ranged between (− 64 to 611 mv) while it was ranged between (227-767 mv) and (125-551 mv) for groundwater and irrigation canals respectively. (El Mahllawy et al. 2013) studied the composition and activation aspects of El Fayoum clays for use as a drilling fluid, the mineralogical analysis of the studied samples using XRD (X-Ray Diffraction) were essentially montmorillonite, kaolinite and illite as clay minerals in a descending order of abundance. Quartz and calcite minerals were also found in minor amounts as non-clay components.
Inverse model calculations were conducted using PHREEQC v.214.03 computer program (Parkhurst and Appelo 1999). Thermodynamic database used: phreeq.dat values are in mol/kg H 2 O. Positive (mass entering water) and negative (mass leaving water) phase mole transfers indicate dissolution and precipitation, respectively, as mmol kg −1 .

Isotopic composition of surface and groundwater samples
Stable isotopes composition (δ 18 O and δD) of 44 surface and groundwater samples were measured to provide basic information on the origin and sources of recharge for the Quaternary aquifer in the studied area. The stable isotope composition of surface and groundwater samples collected together with those of the d-excess values is compiled in  (Fig. 6) (Fig. 6) illustrate that all collected water samples are located within an evaporation line given by a regression equation of slope 6.63 and intercept of 3.98. An enrichment occurs for both δ 18 O and δD values compared with Nile River revealing the occurrence of intensive evaporation process in addition to irrigation return phenomena with an increase of salinity as shown in Fig. 7 for TDS in (mg/l) and δ 18 O in ‰.
TDS values in Fig. 7 show an increment compared with slightly constant values of δ 18 O for collected water samples. This may be due to evaporation and dissolution processes that affect surface and groundwater samples in El Fayoum depression.

Conclusions
The results of this study showed that integration of hydrochemical data, inverse geochemical modeling coupled with isotopic data can help to understand the hydrogeochemical processes controlling water chemistry evolution in El Fayoum depression. The relative abundance of major cations in shallow groundwater is Na + > Mg 2+ > Ca 2+ > K + (on molar basis) and HCO 3 − > Cl − > SO 4 2− for anions. For surface water samples (canals and drains) the predominant cations are Na + > > Ca 2+ > Mg 2+ > K + (on molar basis) and HCO 3 − > Cl − > SO 4 2− for anions. Gibbs plot that illustrate Rock-water interaction which affects Quaternary aquifer system indicates that most water samples fall in evaporation dominant zone. As Nile River being the only recharge source for the depression via Bahr Youssef, the groundwater in the study area is characterized by Na-Ca-Mg/HCO 3 type in 92%, while irrigation canals and drainage water have this type 77%, 58% respectively. The concentration of (SO 4 2− + HCO 3 − ) exceeds (Ca 2+ + Mg 2+ ), which is the result of calcite, dolomite in the system. Gypsum dissolution contributes to increase of both Ca 2+ and SO 4 2− in the aquifer, which leads to a slightly low PH and supersaturation or near equilibrium of calcite. Na + is comparatively higher than Cl − for surface and groundwater samples suggested cation-exchange process that  greatly affects groundwater ion compositions through interaction between water and aquifer materials. Saturation Index using PHREEQC software results indicated that the surface (canals, drains) and groundwater samples were-undersaturated (dissolution is required to achieve equilibrium) with minerals and gases in order of Al (OH) 3 > partial pressure of CO 2(g) > Calcite CaCO 3 > Aragonite CaCO 3 > Gypsum CaSO 4 ·2H 2 O > Anhydrite CaSO 4 > Dolomite CaMg (CO 3 ) 2 . Negative values of PCO 2(g) may reflect the recharge from underlying aquifer and seepage from drains. Surface and groundwater samples have a tendency to precipitate the minerals/gases in order of Hematite Fe 2 O 3 > Goethite FeOOH > Gibbsite Al (OH) 3 > Barite BaSO 4 . The inverse geochemical modeling demonstrated that dissolution of carbonate (calcite, dolomite) and evaporites (gypsum, Anhydrite) minerals with outgassing of hydrogen sulphide as a result of organic matter reduction and microbial processes. Cation-exchange of Ca 2+ over Na + also important geochemical processes in the aquifer system of the area. Environmental stable isotope data revealed that all surface and groundwater samples have enriched values of δ 18 O and δD compared with Nile River revealing the occurrence of intensive evaporation process in addition to irrigation return phenomena that effects on salinity as illustrated for TDS in (mg/l) and δ 18 O in ‰ relationship. Finally, from literature review and field observations, an environmental protection is needed through continuous monitoring of water resources quantity and quality with the possibility of contamination removal and construct sanitation network extends to all regions of El Fayoum depression.
Funding Open access funding provided by The Science, Technology & Innovation Funding Authority (STDF) in cooperation with The Egyptian Knowledge Bank (EKB).
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http:// creat iveco mmons. org/ licen ses/ by/4. 0/.