Spatial and seasonal assessment of physico-chemical characteristics of soil in Wadi El-Rayan lakes using GIS technique

Wadi El-Rayan wetland is an inland Egyptian Ramsar site of international importance, due to its rich biodiversity. The wetland constitutes two man-made lakes. The lakes size, particularly the Lower Lake, is in continuous reduction, and the vegetation productivity pattern is changing as result of deteriorating water and soil quality, domestic and industrial activities. This research was carried out to map the spatial and seasonal distribution of certain physico-chemical parameters of geo-referenced soil samples covering the two lakes parts using ordinary kriging GIS technique. The study utilized ArcGIS ver.10.7 software, and 76 soil samples collected from representative sites; 38 samples during winter and 38 samples during summer seasons. Soil data were analyzed statistically and geo-statistically based on best fitted semivariogram model. The results indicated high concentrations of some soil physical and chemical parameters, especially in summer and in the Lower Lake of Wadi El-Rayan. The summary statistics for soil parameters had shown that the median of some soil parameters is higher than the mean which indicated the presence of abnormal data. This finding has been approved through the spatial distribution mapping which clearly showed the differences in spatial and temporal distribution between the upper and lower lakes of Wadi El-Rayan. The study concluded that the generated spatial distribution maps using ordinary kriging method can be used as an effective tool in Wadi El-Rayan wetlands’ soil management.


Introduction
Sustainable wetland management generally requires increased management ability and input [25]. Water, soil and vegetation are among the main components of any wetland ecosystem. Soil quality is a reflection of water quality, and it determines the soil capacity to produce biomass [10]. Soil physical and chemical parameters play an important role in wetlands functions as they influence their vegetation productivity and indirectly their biodiversity. Numerous previous studies concluded that there is a significant correlation between wetlands vegetation productivity and soil quality [23,29,44,55]. However, those parameters are characterized by high spatial and temporal variability depending on climate, season, topography, land use, water quality, anthropogenic activities (e.g., grazing), as well as management. Hence, understanding spatial and temporal variability of soil physical and chemical parameters are essential for devising site specific parameter management strategies with the aim of better wetlands management and increased sustainability in wetlands vegetation productivity.
Geo-statistical methods are used to interpolate the spatial variability of soil parameters and also to predict an unknown soil parameters value at a certain locations. Geostatistical methods are based on [51] low, which states that "everything is related to everything else, but near things are more related than distant things." They have been applied since 1970s, and many methods have been developed through the time such as; Inverse Distance Weighting (IDW) [49], Ordinary Kriging (OK) [12,57], Simple Kriging (SK) [54], Co-Kriging (CK) [37], pedo-transfer functions [43] and artificial neural network [53]. Recently, OK has been most widely used in many studies that map spatial distribution of soil physical and chemical parameters in different soils of Egypt [26,16,17,15,36] due to its efficiency in predicting and interpolating specific soil parameter on spatial distribution maps [50].
Wadi El-Rayan wetlands are part of Wadi El-Rayan Protected Area (Ramsar site) which covers an area of 1759 km 2 [31]. Wadi El-Rayan wetland constitutes two man-made connected lakes formed during 1970s. One of the reasons behind there formation was to slow-down the increase in the water table in El-Fayoum main depression in Egypt. The construction of the two Lakes led to restoration of the severely affected biodiversity of Wadi El-Rayan Protected Area and attraction of large populations of birds (e.g., water fowl) which came after vegetation cover; mainly common reeds; distributed along the shores of the two lakes (Saleh 1987). Wadi El-Rayan lakes get their water through around 9 km length tunnel from El-Wadi drain. The tunnel was receiving about 8 m 3 /s in 1999, and this amount decreased to about 4.5 m 3 /s in 2008 [19]. Based on the reduction in total received drainage per year and the rainfall pattern, besides the fact that the Lower Lake is a closed basin with no outflow, these all led to the accumulation of salts in Wadi El-Rayan Lakes, particularly the Lower Lake [1]. Soil salinity is the most soil limitation factor for vegetation productivity [36]; however, other soil physical and chemical characteristics can also influence it. Wetland's cover classification is an important tool in wetlands management [14]. According to land cover classification performed on selected two Landsat imageries for the years 1990 and 2020, the total area of the wetlands' surface water and vegetation cover are in continuous reduction and open soil is becoming more exposed. The water body class in Wadi El-Rayan wetland has declined by 30.98 km 2 at an annual rate of 1.03 km 2 /year. While, the vegetation cover have lost over 80% of its area over the past three decades (Table 1).
Therefore, diagnosis of soil physical and chemical parameters-related limitations assumes a greater significance to assess the degree of impact of Wadi El-Rayan wetlands water quality on their soil. Assessment of spatial variability of different soil parameters is a viable option to identify the most affected sites. This will enable Wadi El-Rayan wetlands mangers to delineate the most productive sites and identify the critical impact zones.
There have been few studies and insufficient data on the physico-chemical characteristics of the study areas' soils. Those studies include the work of [2], who studied the 16 physical and chemical soil variables in order to evaluate the soil-vegetation relationship in Wadi El-Rayan wetlands and found that soil salinity, pH, calcium carbonate, organic matter and the relative concentration of the major cations were the most important variables affecting the distribution of the vegetation pattern. El-Hennawy [20] studied the soil physical and chemical characteristics of the two lakes and found that the salinity level is particularly higher in the Lower Lake comparing to the Upper Lake. He also revealed that the wetland vegetation productivity in a certain site is correlated with its correspondent environmental variables. Sayed [47] studied the soil chemical properties of Wadi El-Rayan lakes and revealed that chlorides, sulfates, sodium and magnesium were the major ions concentrations recorded during the study and that the Lower Lake recorded higher concentrations than those of the Upper Lake. He also said that the rate of evaporation and the decreased water flow rate from the Upper Lake are the major factors controlling these concentrations. Abdou [4] studied six chemical concentrations of soil samples collected from Wadi El-Rayan Lakes and found that the two lakes soil nitrites content was the highest recorded concentration among the selected samples, and he also attributed his results to non-point sources soil pollutants such as agricultural and industrial wastes and runoff.
However, up to this date, no studies were reported on the spatial and/or temporal distribution of soil physical and chemical parameters of Wadi El-Rayan wetlands. In this study, we attempt to use ordinary kriging method (OK) to address this knowledge gap by mapping variety of soil physical and chemical parameters spatialy according to their standard soil classification. The objectives of this study are, to assess the spatial variations in soil physical and chemical parameters in different seasons as well as study their status.

Description of the study area
The area under investigation belongs to Wadi El-Rayan protected area located in a depression about 140 km southwest of Cairo, in El-Fayoum Governorate, Egypt [31]. It is located between longitude 29°00′00″ and 29°24′11″ E, and latitude 30°00′00″ and 30°34′00″ N. Wadi El-Rayan is located in hyper-arid Saharan climate zone with hot summers and mild winters [56]. The mean elevation of the area is 43 m below sea level. According to metrological data recorded by the nearest station to the protected area in El-Fayoum over the past fifty years, the mean temperature during winter and summer is 28.5 °C and 13.7 °C, respectively, and the average precipitation rate is 10.1 mm per year, whereas the highest rainfall occurs in winter and there is zero rain fall precipitate in summer. The mean annual aridity index is 0.004, which is due to the high evapotranspiration rate throughout the year [46]. The main habitats in protected area are desert, lakes, oasis and wetlands [24]. Hawela and El-Khatib [27], classified the soils in Wadi El-Rayan as Typic Torriorthents, Typic Torripsamments and Typic Salorthids and a number of sub-land types. Wadi El-Rayan Protected Area was declared a Protected Area by Prime Ministerial Decree (No. 943) in1989 according to the law No. 102/1983 of the protected areas in Egypt. In 2012, the protected area was declared one of the four Egyptian RAMSAR sites of international importance [13]. The area has drawn the attention of national and international visitors due to its natural landscapes, biodiversity, sand beaches, bird watching, camping, waterfalls and the Valley of Whales in Wadi El-Hitan which was designated by the UNESCO as the first World Heritage Site of natural category located in Egypt. Wadi El-Rayan Protected area is visited by over 200 thousands visitors annually [24].

Soil sampling
Seventy-six soil samples were collected from 18 representative sites covering Wadi El-Rayan Lakes' parts (upper and lower lakes) during the year of 2019 for two times. Thirty-eight samples were collected during winter, and 38 samples were collected during summer from the same locations. SW-GIS data collector android application used to record a track of each sampling point, which helped us to take the second season samples from the same locations with the same coordinates. All the soil samples were collected at depth of 0-30 cm in duplicate in labeled plastic bags and carried out in the Ecology Laboratory, Department of Botany and Microbiology, Faculty of Science, Al-Azhar University, Cairo, Egypt.

Laboratory analysis
Sand/Silt/Clay Content was determined using hydrometer method according to Estefan et al. [22]. Samples were dried in an oven at 40 °C for 72 hours, sieved with a 2 mm sieve and stored for further analyses. Water Holding Capacity was determined after Estefan et al. [22]. Soil pH was determined in soil extract for each sample (1 g soil: 5 ml distilled water) using digital portable pH meter Adwa® AD11. Organic Carbon (OC) and Organic Matter (OM) were determined using titration method according to Piper [39]. Total Dissolved Salts (TDS) and Electrical Conductivity (EC) were both measured using digital portable TDS meter Adwa® AD201, in soil extract (1 g soil: 5 ml distilled water). Calcium Carbonate (CaCO 3 ) was determined using titration method described by Estefan et al. [22]. Chlorides (Cl) and Bicarbonates (HCO 3 ) were determined by titration method according to Jackson (1967). Sulfates (SO 4 ) were determined using the turbid metric method after Estefan et al. [22]. Calcium (Ca) and Magnesium (Mg) were determined by titration method according to Page, Miller and Keeney (1982). Potassium (K) and Sodium (Na) were determined using Flame photometer at 767nm wavelength and Flame photometer at 589 nm wavelength, respectively, according to Estefan et al. [22]. Nitrates (NO 2 ) and Nitrites (NO 3 ) were determined using UNICO 2100 UV Visible Spectrophotometer according to Greenberg, Clesceri, and Eaton (1995). Phosphorous (P) was determined using UNICO 2100 UV Visible Spectrophotometer according to Allison et al. [5].

Statistical analysis
Minitab ver.18 was used to produce descriptive statistics (mean, median, standard deviation, range, minimum and maximum) for different soil parameters during both winter and summer seasons.

Geo-statistical analysis
Geo-statistical methods, including Ordinary Kriging (OK), semivariogram and cross-validation, were used to assess the spatial and temporal distribution of 17 selected soil physical and chemical parameters. Geo-statistical analysis was done using geo-statistical analyst tool in ArcGIS ver.10.7 software through application of ordinary kriging method.
OK is the major widely used type of kriging in spatial distribution researches [33]. It is expressed as; where z * (u) is the measured value of z at location (u) ; ∖ a is the weight associated with the estimated value of z at location a and N is the number of estimated values used in measurement in the neighborhood of a.
According to ESRI [21], the semivariogram is defined as Where s i and s j are pair of points close to each other in terms of the distance and Z s i − Z s j are the difference in their values. The semivariogram assumes that when the distance between two pair of point is small, the difference in their values is expected to be similar and vice versa.
There are 11 models can be fitted to the experimental semivariogram, from which we chosen three models; circular, exponential and Gaussian models in order to select the best fit model with Mean Standardized (MS) close to zero and Root Mean Square Standardized (RMSS) close to 1 [41]. The cross-validation method was then applied to validate the model parameters using three major parameters; nugget, sill and range, which assessed the selection of most appropriate model to interpolate and predict the different soil physical and chemical parameters in the study area [36].

Results and discussions
Results of soil texture analysis of Wadi El-Rayan wetlands during winter and summer showed that the main soil textures that can be found are mainly sandy loam, loamy sand (distributed around the Upper Lake) and sand (distributed around the Lower Lake), and these results are in accordance with Abbas and Afefe [1].

Descriptive statistics
The median of some soil parameters, sand, water holding capacity, pH, total dissolved salts, sulfates and calcium during winter and sand, pH, calcium carbonate and sodium during summer, were higher than the mean, which indicates that the effects of abnormal data were significant ( Table 2). Soil water holding capacity (WHC) ranged from 10.4-60.56% to 20.7-56.7% during winter and summer, respectively. The mean value of soil WHC during summer was higher during winter, and it showed a trend with soil texture. These results agreed with Bandyopadhyay [7], who said that the soil WHC has direct relation with the soil texture and the higher the soil silt and clay percentages; the higher the soil WHC. During winter, pH values of the soil water extract ranged from 7.2 to 8.1, while during summer they ranged from 7.4 to 7.8. According to Proffitt [40], the mean value of pH of soil extract of Wadi El-Rayan wetlands was strongly alkaline during winter and slightly alkaline during summer, and these agree with the study of El-Hameed et al. [18]. The soil organic carbon (OC) and organic matter (OM) varied from 0.12 to 2.88%, 0.2-4.95% during winter and from 0.12 to 2.28%, 0.2-3.92% during summer, respectively. Organic carbon and organic matter play a very important role in soil fertility in Egypt [3]. According to Proffitt [40] classification, most studied soil samples had low OC and OM percentages. These results agreed with Dregne [11], who reported that soil of arid areas usually have low level of organic matter. Soil total dissolved salts (TDS) and Electrical conductivity (EC) ranged from 186 to 4560 ppm, 0.29-9.34 mS/cm during winter and from 533 to 7370 ppm, 0.83-11.51 mS/cm during summer. According to [42] classification, the mean value of TDS of soil extract of Wadi El-Rayan wetlands was moderately saline during winter and summer and the mean value of EC of soil samples was slightly saline during winter and moderately saline during summer. However, salinity concentration was higher during summer season compared to the winter season; this result is in accordance with [48], who said that the concentration of TDS in hot period (summer) is much higher than cold period (winter). Generally, the high soil salinity concentration could be attributed to reduction rate of drainage water inflow for reclamation projects and the high evaporation rate [19,47]. Soil calcium carbonate content ranged from 0.3-21.9% to 0-22% during winter and summer, respectively. The mean value of soil calcium carbonate during winter (11.64%) was lower than during summer (16.55%). The high percentages of soil calcium carbonate could be attributed to the accumulation of large amounts of shell remains. The results of soil calcium carbonate and pH analysis indicate that the soil is calcareous, and this is according to 58], who said that the calcareous soils normally have pH values between 7.5 and 8.5 and calcium carbonate between 5 and 8% or greater. Generally, the concentration of major anions, sulfates and chloride, reflect that of TDS because they also tend to increase in the direction of flow; however, bicarbonate does not show the same trend relative to other anions. The mean value of soil sulfates content was higher during winter (114 mg/L) than during summer (96.68 mg/L). According to [30] classification, soil sulfates during winter and summer were very high (> 20). Sulfate is an indication of pollution, and this pollution could be due to industrial discharge, treated sewage release and K-sulfates fertilizers.
Results of soil chlorides analysis revealed that the soil chlorides content ranged from 20-700 mg/L to 100-1300 mg/L during winter and summer, respectively. According to [38] guidelines, most of the studied sites ranged from low during winter to moderate during summer and these results are in accordance with [1]. The four most abundant cations in soils are Calcium ( Ca +2 ) , magnesium ( Mg +2 ), Potassium ( K + ) and sodium ( Na + ). The cations zinc ( Zn +2 ), manganese ( Mn +2 ), iron ( Fe +2 ), aluminum ( Al +3 ) and copper ( Cu +2 ) are usually present in amounts that do not contribute significantly to the cation complement. Therefore, it is common practice to measure only the concentration of the four abundant cations [28]. The soil calcium content range was between 200-1600 mg/L and 200-2000 mg/L during winter and summer, respectively. According to [35] classification, the mean values of soil calcium content were low during the two seasons. Soil magnesium content ranged from 120-1080 mg/L to 250-3582 mg/L during winter and summer, respectively. According to [30], the mean value of soil magnesium content ranged from high during winter to very high during summer. Soil potassium content ranged from 484.57-589.53 mg/L to 529-1455.33 mg/L during winter and summer, respectively. The mean values of soil potassium ranged from high during winter (541.3 mg/L) to excessive during summer (946.8 mg/L) [30]. Soil sodium content ranged from 305-1760 mg/L to 602.5-2862.5 mg/L during winter and summer, respectively. The mean values of soil sodium were very high (> 480 mg/L) during the two seasons [35].
According [34], the mean values of soil nitrites and nitrates content were low (< 20 mg/L) and (< 10 mg/L) during winter and summer, respectively. Soil phosphorous content ranged between low and moderate during winter, and low and excessive during summer [34]. The mean values of soil phosphorous content during summer was much higher than during winter, these could be attributed to the accumulation of deposits containing dead planktons which increase the phosphorous percentage in the soil and according to [24], WadiEl-Rayan Protected Area is visited by over 200 thousands visitors annually who come to the area for recreational activities (boats, cafeterias and camping) around the Lower Lake especially during summer season which contribute to the increase in the pollution rate in the lake.

Semivariogram analysis
Semivariogram parameters of the soil physical and chemical data of the winter and summer seasons and their validation information that used to evaluate the goodness of fit are showed in (Table 3). These semivariogram parameters are the nugget value, sill (summation of nugget and partial-sill), nugget/sill ratio, range and lag size. Generally, positive nuggets were observed for all the soil data and this could be explained by short range variability, random and inherent variability and sampling error.
For bicarbonates, the variance of the input data was too small to be calculated; however, the rest of the soil physical and chemical parameters recorded lag size less than half of the range (Table 3). According to [32], the sampling interval (lag size) should not exceed the half of the range. Hence, using this sampling design was appropriate for the current study, and it is foreseeable that good spatial structures will be shown on the interpolated maps. According to [45], the nugget/sill ratio is used to classify the spatial dependence of soil parameters and to quantify short distance autocorrelation. Low (< 25%), Moderate (25-75%) and High (> 75%) nugget/sill ratios indicate strong, Moderate and weak spatial dependence or spatial autocorrelation, respectively [9]. In the present study, the nugget/sill ratio during winter showed that TDS, Cl, SO 4 , Ca, Mg, Na, NO 2 and NO 3 were strong spatially dependent (0.26-20.45 %), and pH, OC, OM, EC and CaCO 3 were moderate spatially dependent (34.62-69.41 %). The nugget/sill ratio during winter also showed that K and P were weak spatially dependent (82.3-95.5 %). The nugget/sill ratio during summer showed that pH, OC, OM, TDS, EC, CaCO3, Cl, SO4, Ca, K , Na and NO 3 were strong spatially dependent (1.72-23.01 %), and WHC, Mg and NO 2 were moderate spatially dependent (25.31-53.68 %). The nugget/ sill ratio during summer also showed that P was weak spatially dependent (93.83%).
The strong and moderate spatially dependence of total dissolved salts, electrical conductivity and their correspondent anions and cations during winter and summer seasons could be attributed to the high evaporation rate and reduction in the inflow rate of the Lower Lake of Wadi El-Rayan wetlands; however, the strong and moderate spatial dependence of nitrates and nitrites during the two seasons could be attributed to the fact that the source of the two Lakes is agricultural drainage water which is rich with fertilizers. The strong and moderate spatially dependence of organic carbon and organic matter during the study period could be attributed to the grazing activities particularly, around the Upper Lake' shores. The weak spatially dependence of phosphorous during winter and summer and weak spatially dependence of potassium during winter season indicated that the spatial patterns of this soil parameters were not influenced by those parameters, however, influenced by factors such as evaporation rate, inflow rate, quality of drainage water and anthropogenic activities such as grazing, aquaculture Fig. 1. Among selected three models, there was no best fitted model for all the soil parameters; however, the model varied according to the soil parameter and season. Exponential model was best fitted to WHC, OC, OM, Ca, Mg, Na and P during winter, and WHC, EC, CL, SO4, Mg, NO 2 during summer. Similarly, circular model was best fitted to pH, TDS, EC, Cl, SO 4 , K and NO 3 during winter, and OC, OM, K, Na and NO 3 during summer. Furthermore, Gaussian model was best fitted to CaCO 3 and NO 2 during winter, and CaCO 3 , Ca and P during summer Fig. 2.

Spatial, seasonal distribution and spectrum
Ordinary kriging method was applied to estimate the value of the different soil physical and chemical parameters at un-sampled sites using three semivariogram models. The soil parameters data showed normal distribution and did not require any transformation before computing the spatial distribution maps. Raster maps (30 * 30 m) of different soil parameters were generated for the study area during winter and summer seasons (Figs. 3, 4, 5, 6 7 and 8). Spectrum of different soil physical and chemical parameters according to some soil interpretation guidelines is shown in Table 4.
Spatial maps of WHC during winter and summer seasons showed that the WHC value was mainly low, which reflects the soil sandy texture nature of Wadi El-Rayan wetlands. The spatial maps of the soil pH showed that the soil pH was alkaline and was significantly higher in winter than summer. This alkalinity may be due to the exchangeable cations, which increase the soil pH. The spatial maps of soil organic carbon showed that the OC value in the Lower Lake ranged from very low during winter to low during summer, while in the Upper Lake it ranged between low and moderate during the two seasons. Similarly, the spatial maps of soil organic matter showed that the OM value Fig. 1 Landsat Scene of Rayan Lakes in Egypt, and sampling points' locations map in the Lower Lake was low during the two seasons, while in the Upper Lake it ranged between low to moderate and did not differ significantly between the two seasons. The spatial maps of TDS and EC showed that the Lower Lake salinity is much higher than the Upper Lake, which showed low to moderate salinity levels during the two seasons. Spatial maps of soil Cl showed that the Lower Lake soil chloride content was higher than the Upper Lake and did not differ significantly between the two seasons. Spatial maps of soil SO 4 showed that soil sulfates content in the two lakes during the two seasons was very high. The spatial maps of the soil exchangeable cations (Mg, K and Na) showed that the soil magnesium, potassium and calcium content of the two lakes ranged between high and very high and did not show significant difference between the seasons; however, spatial maps of soil Ca showed that the soil calcium content in the Lower Lake was higher than the Upper the Lake during the two seasons. The spatial maps of soil NO 2 and NO 3 content showed that there were no significant spatial and temporal deference of soil nitrites and nitrates content between the two lakes and that their values were mainly low. The spatial maps of soil P content showed that the soil phosphorous content in the two lakes was moderate during winter, while it was high during summer.

Cross-validation results
Evaluation of performance of the ordinary kriging method was done through cross-validation procedures using mean prediction error (MPE), mean standardized prediction error (MSPE) and root mean square standardized prediction error (RMSSPE) of the spatial maps of soil parameters (Table 5). Some parameters showed high MPE which explain the high spatial difference of those parameters between the two lakes of Wadi El-Rayan wetlands. This agreed with the overall research findings which showed that soil of the Lower Lake is much deteriorated comparing to the Upper Lake. However, MSPE and RMSSPE values were low as they ranged from (0.075-0.103) to (1.37-0.822), respectively, indicating that OK provided superior prediction and mapping accuracy for Wadi El-Rayan soil parameters.

Conclusion
This study provides the first comprehensive analysis of physical and chemical parameters of Wadi El-Rayan wetlands' soils. The results indicated high concentrations of some soil physical and chemical parameters, especially in summer and in the Lower Lake of Wadi El-Rayan. The summary statistics for soil parameters had shown that the median of some soil parameters is higher than the mean which indicated the presence of abnormal data. This finding has been approved through the spatial distribution mapping which clearly showed the differences in spatial and temporal distribution between the upper and lower lakes of Wadi El-Rayan. Semivariogram models were fit for all soil properties, and the best fitted semivariogram model for each parameter was identified through crossvalidation procedures. Exponential, Gaussian and Circular models performed well in analyzing the spatial variability of the 17 soil physical and chemical parameters contents. A strong and moderate spatial dependence of soil parameters was observed, indicating that the soil parameters were controlled by evaporation rate, inflow rate, quality of drainage water and anthropogenic activities such as grazing, aquaculture. Cross-validation of semivariogram models through ordinary kriging showed that the spatial prediction and interpolation of soil parameters is preferred than assuming the mean of the observed values at any unstamped site. Finally, the generated spatial distribution maps using best fitted semivariogram models and   ordinary kriging method can be used as an effective tool in Wadi El-Rayan wetlands' soil management. Recommendations • The amount of water inflow into the two lakes should be measured and controlled in order to solve the high salinity problem, particularly the Lower Lake. • Regular monitoring of the soil quality of the Wadi El-Rayan Lakes is important for sustainable wetlands management. • Great efforts and cooperation between wetland managers, decision makers and other concerned authorities are needed in order to keep the wadi El-Rayan wetlands' biodiversity.
• Application of GIS and geo-statistical mapping, spatial distribution of soil and water physical and chemical parameters is helpful in wetlands management.

Future work
The authors are currently working on evaluation of vegetation productivity of Wadi El-Rayan and correlation of distinctive productive sites to their correspondence environmental gradients. Furthermore, the authors are also working on monitoring wetland change detection of Wadi El-Rayan Protected Area in Egypt (RAMSAR site) over the past three decades using remote sensing and GIS.
Funding The authors would like to extend their sincere thanks to the Pan African University for Life and Earth Sciences including Health and Agriculture for providing the research grant to the first author.

Data Availability
The datasets used and/or analyzed during the current study are available from the corresponding author on reasonable request.

Compliance with Ethical Standards
Conflict of interest The authors declare that they have no conflict interest.
Ethical approval and consent to participate A research site permission letter was obtained from the concerned institution. An effort was made to conduct the research in an ethical manner.

Consent for publication
All authors agreed and approved the manuscript for publication in Environmental Systems Research.
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/.