Revealing the geomorphic parameters of neotectonics based on geological and seismic data at Wadi Beda basin, NW Gulf of Suez, Egypt

The current study used remote sensing, GIS, and seismic data to determine the geological and morphometric parameters of the Wadi Beda Basin area, which is located in the northwestern part of the Gulf of Suez Province. The geomorphic parameters related with the development of neotectonics are deduced using the Shuttle Radar Topography Mission (SRTM), DEM-derived catchments, and drainage networks of the research area. The most important geomorphic parameters applied to the current study area are valley floor width to height ratio (Vf), stream length gradient index (Sl), transverse topographic symmetry index (T), hypsometric curve and integral (HI), and drainage basin shape (Bs). Based on the obtained geomorphic parameter results, the Wadi El-Noqra and Wadi Umm Athala basins exhibit moderate tectonic activity and moderate symmetry and are considered to be at their mature stage. The main basins of Wadi Umm Elda and Wadi Beda are moderately symmetrical and reflect the youth stage. However, the Bs index suggests that the Beda main basin has moderate tectonic activity, whereas Wadi Umm Elda has strong tectonic activity. By integrating the geomorphic and geological data, it is possible to determine that the majority of the study area (83.4%) is categorized by a moderate risk. Low-risk and high-risk zones, on the other hand, make up 11.1% and 5.5%, respectively, of the study area. The majority of the wadis in the Wadi Beda Main Basin are controlled by NW-trending faults, and the perturbation sites are concentrated there where these faults cross with E-W-trending faults, making them high-risk areas.


Introduction
Currently, the northwestern Gulf of Suez Province is observing accelerated, sustainable development and many national projects implementation. Many essential infrastructures such as highways, ports, power plants, and urban extension are involved in those projects. Concurrently, the different risks, which may threat those developing projects, dictate conducting hazard assessment to mitigate its impacts on such projects.
There are different definitions of neotectonics (e.g., Mörner 1978;Pavlides 1989;AGI 2009), but the most common accepted one between researchers is the study of young tectonic movements. Furthermore, there is a disagreement on exactly how young they must be considered as neotectonics. Hancock and Williams (1986) proposed that neotectonics began when the present-day stress field of a region was started. Generally, it started between the Oligocene and Miocene periods, which matches Vita-deformation of specific area but also the tectonic activity level (Keller 1986;Keller and Pinter 1996;Demoulin 1998). The most important geomorphic parameters include the ratio of valley floor width to valley height (V f ), the stream length-gradient index (S L ), the transverse topographic system (T), the hypsometric integral (H I ), the drainage basin shape (B s ), and the asymmetric factor (A f ) (e.g., El-Hamdouni et al. 2008;Buczek and Górnik 2020).
Mainly, the present study aims to investigate, analyze, and integrate the geomorphic and geological parameters and, hence, assess the neotectonics of Wadi Beda Main Basin and its sub-basins; Wadi Umm Athala, Wadi El-Noqra, and Wadi Umm Elda (Fig. 1). Such assessment sheds light on the hazardous sites that must be taken into consideration during the future development and urban planning of the Northwestern Gulf of Suez Province.
Several investigations have been carried out either at the whole northwestern Gulf of Suez Province or at a specific part of it. The morphometric analysis which was carried out by Arnous et al. (2011) strongly supported the high probability of flash flooding in different sites of the northwestern Gulf of Suez Province. Geriesh et al. (2004) suggested a suitable management system to control its huge runoff quantities based on an estimation of the flood risks of Wadi Ghoweiba Basin. El-Behiry et al. (2006) deciphered the flash flood risk zones and established an effective hazard mitigation scheme based on the analysis of runoff hazard of the Ain Sukhna industrial area and delineated the drainage networks for predicting surface water runoff. El-Rayes et al. (2009) proposed a system composed of a group of cisterns, retention, and reservoir dams in addition to mitigation canals to control the huge runoff quantities of water at Wadi Hagul Basin. Seleem and Aboulela (2011) evaluated deformation and surface break at the Wadi Hagul area and gave a new insight for understanding and assessment of the seismic hazard. Arnous et al. (2011) concluded that the seismic activity in the northwestern Gulf of Suez region markedly increases from north to south. Hegazi et al. (2013) revealed that the Bitter Lakes area, including Wadi Beda Basin, display four alternating NW-trending seismic zones characterized by unique motion. Two of these seismic zones (III and IV) are displayed at the present study area (Fig. 2). El-Rayes et al. (2015) address the correlation between seismological data and the geographic distribution of geothermal fields in the neighboring South Sinai. They concluded that the distribution of recent earthquake epicenters over South Sinai reveals three broad active trends (N 6° E, N 40°-75° E, and N 50°-60° W), which are consistent with the long fault extensions and include the well-known hot springs in the region. Omran et al. (2021) examine the connection between seismic activity and the potential of rock fall in the nearby South Sinai. They demonstrate that the possible occurrence of rock fall is highly sensitive to the depth of seismic activity, geology, and lineament density.
These previously published investigations have been focused either on seismic, tectonic, runoff, or flash flood parameters. None of them has taken into consideration the integration between the different geomorphic and geological parameters in order to evaluate the neotectonics at the study area.
The goal of the current study is to examine, analyze, and integrate geological and geomorphic parameters in order to assess the neotectonics of the W. Beda Main Basin and its sub-basins, W. Umm Athala, W. El-Noqra, and W. Umm Elda. Such a neotectonic assessment sheds light on the hazardous sites that must be addressed for the future development and urban planning of the Northwestern Gulf of Suez Province.

General geology
Wadi Beda main Basin belongs to the Northwestern Gulf of Suez Province, and it encompasses three sub-basins: the Wadi Umm Athala, Wadi El-Noqra, and Wadi Umm Elda. The Wadi Beda main Basin is bounded by longitudes 32° 22′ and 32° 55′ E and latitudes 29° 41′ and 29° 59′ N, covering an area of about 800 km 2 (Fig. 1).
Many highlands are encountered in the present study area such as Gebel Ataqa block (900 m a.s.l.) to the northeast, and El-Galala El-Bahariya Plateau (1332 m a.s.l.) to the south as well as the Wadi Beda main basin (Fig. 1), and several sub-basins, e.g., Wadi El-Noqra, Wadi Umm Athala, and Wadi Umm Elda.
The exposed rocks at the study area belong to Phanerozoic Era, where they are ranging in age from Eocene to Recent time (Fig. 3). Structurally, the Northwestern Gulf of Suez Province, including the present study area, represents a series of alternating positive and negative fault-blocks which are bounded by ENE-trending faults (Abd El- Rahman and El-Etr 1979). Three of these blocks are identified at the study area; the two positive blocks (topographically high) are El-Galala El-Bahariya Plateau and Gebel Ataqa, whereas the negative block (topographically low) is Wadi Akheider (Fig. 1). Gamal et al. (2021) revealed that competent rock units (Middle Eocene limestone) of Gebel Akheider forming horsts and grabens They are categorized under three main sets of faults; these are NW-, E-W-, and NNW-trending sets. The NW-trending set is the most prominent one in the study area ( Fig. 3) (Abdallah 1993;Youssef and Abdallah 2003).

Materials and methods
Microwave data such as Shuttle Radar Topography Mission (SRTM) images provide a rapid and low-expensive means to analyze and interpret regional geomorphic features (Grohmann et al. 2007). These data are commonly used in accomplishment of digital terrain models, such as shaded relief (Chorowicz et al. 1998;Jordan et al. 2005;Nureddin et al. 2009;and Abdel Kareem et al. 2011). They also provide crucial information for morphotectonic analysis such as geomorphic parameters, drainage asymmetry and networks, and stream orders.
In the current investigation, the SRTM images have been used in producing DEM, contour, and relief maps and extracting drainage networks of the study area. Furthermore, a 3D-perspective with image draping could be useful for the morphotectonic investigations and to sharpen the topography, particularly for cases where the topographic data and optical images are integrated (Le Turdu et al. 1995;Jordan et al. 2005).  The (V f ) index was calculated as follows: where V fw is the width of the floor, E ld and E rd are the elevations of left and right divides, respectively, looking downstream, and E sc is the elevation of the floor.

Stream length-gradient index (S L )
The (S L ) index was computed based on Hack's Eq. (1973) as follows: where ΔΗ/ΔL is the stream gradient at a specific reach (point) in the channel, ΔΗ is the change in elevation of the reach, ΔL is the length of the reach, and L is the channel length from the divide to the midpoint of the channel reach. (1)

Transverse topographic symmetry index (T):
The (T) index was calculated regarding the larger axis of the basin as follows: where D a is the distance from the basin midline to the active meander-belt midline, and D d is the basin divide.

Hypsometry index (H i )
The (H i ) index was calculated relied on Pike and Wilson's Eq. (1971) as follows: where h max , h min , and h mean are the maximum, the minimum, and the mean elevation, respectively.

Drainage basin shape (B s )
The (B s ) index, based on Ramirez-Herrera (1998), can be formulated as  Table 1 The calculations of (V f ) index of sub-basins Values in bold indicate profile numbers as illustrated in figure (9) Sub-basin Profile no where B l is the length of the basin measured from the headwater to the outlet, and B w is the basin width in the widest point of the basin.

Structural lineaments
Lineaments have been extracted from the enhanced STRM images, analyzed, and integrated with lineament parameters (length, number, and directional families). Linear-features extraction is executed by manual lineaments extraction, where the confirmed structural lineaments are extracted from different enhanced satellite images by using visual interpretation. The resulted structural lineaments map of the study area is shown in Fig. 7.

Seismicity data
The epicenters, which have been detected at the study area during the last century (1900-2003) (Earthquake catalog of Egypt, 1900Egypt, -2003, are displayed on the map, which is shown in Fig. 8. The uniformity and incompleteness of seismological data has been carried out based on an equation proposed by Al-Ibiary (2000): where N is the number of observation, I i indications of the seismic intensity at K th location for i th earthquake, and I i,max is the observed maximum intensity for the i th earthquake.

Results and discussion
The attained results of geomorphic parameters are discussed in this section. The sub-basins of Wadi Beda Basin along with the transverse topographic profiles on which the measurements carried out are shown in Fig. 9. An example of the extracted topographic profiles and the structural cross-sections of the different sub-basins (based on field mapping), the topographic profiles of Wadi El-Noqra sub-basin are shown in Fig. 10.

The geomorphic and geological parameters
Generally, the distribution map of the V f index shows a tendency to increase from the upstream to the downstream (Fig. 11) Athala sub-basin (0.43-11.95) are representing ideal cases for the tendency to be gradually-increased from upstream to downstream. The V f values of those sub-basins start low at their upstream and then gradually increase and become higher near the mountain fronts indicating broad valley floors and lateral cutting of the streams. The V f index was originally used to distinguish V-shaped valleys from U-shaped valleys (Bull and McFadden 1977). Wadi Umm Elda is U-shaped at the upstream, whereas as it becomes V-shaped at the downstream. On the other hand, Wadi El-Noqra and Wadi Umm Athala have V-shaped valleys and are considered areas of active uplift and deep, linear stream incision. Wadi Beda main has a wider and flat floor across its course. The relatively high V f values indicate tectonically stable areas where stream valley bottom tend to be flat and wider and eroding laterally into the adjacent hill slopes. The V-shaped valleys are due to the sculpting of streams along the NW-trending faults related to the tectonics of Gulf of Suez Rift (Fig. 3), as an example of these valleys is Wadi El-Noqra which is pursuing its course along one of the rift-related graben (Fig. 12). The careful examination of the geometry of that graben shows the presence of slightly tilted fault-blocks toward the SW direction, and hence, slightly shifting of valley course toward the same direction.
The S L index investigates the relationships between tectonic activity, rock resistance, and topography expressed by sudden changes in channel slope. The differences in the S L values indicate changes in the longitudinal profile of a stream and hence reflect the possible tectonic activity. Baruah et al. (2020) concluded that both stream length gradient index (SL) and steepness index (SI) have their common genetic relationship between the surface uplift and response of any channel gradient to it.
The distribution map of the S L index (Fig. 13) shows a tendency to increase from the upstream to downstream. Based on the attained results, S L parameters change between 23.84 and 287.02 (Table 2). The streams display a smoothly changing, concave longitudinal profile, except some locations that are tectonically perturbed (Fig. 14). The S L values are low at upstream and suddenly became high at sites where streams crossing fault and then gradually increase to reach the maximum value at the downstream. Such perturbations in stream profiles could be attributed to tectonic activity, namely, E-W-trending normal faulting (Fig. 3).
Transverse topographic symmetry factor (T) is used for determining the possible tilt direction (e.g., El-Hamdouni et al. 2008), where the perfectly symmetric basin has a value of transverse topographic symmetry (T) equals zero, as the asymmetry increases, (T) increases and approaches the value of one. The distribution map of the factor (T) shows that the Wadi Beda is moderately asymmetric (Fig. 15), except for some highly asymmetrical parts of the streams, in which the factor approaches one (0.93) ( Table 3). An exception to these results is encountered in the middle part of Wadi El-Noqra and downstream stream of Wadi Umm Athala where (T) values approach zero, indicating symmetry (Table 3) (Fig. 15).
The asymmetry of the sub-basins could be interpreted in terms of the regional tectonic setting of the North-western Gulf of Suez Province. Where tilted fault-blocks style is prevailing, in which blocks between roughly parallel faults are subjected to tilting. The tilting produces mountains and intervening valleys. Three tilted fault blocks are exposed in the study area; these are Gebel Ataqa, Gebel Abu Treifiya-Khahaliya-Umm Zeita, and Gebel Qattamiya, which are separated by three main wadis; Wadi Hagul, Wadi Beda, and Wadi Ghoweiba (Fig. 3). These three tilted fault blocks are composed of Eocene to Miocene rocks are affected by E-W elongated belts of left-stepped en echelon normal faults. These belts act as transfer zones between NW oriented normal faults synchronous with faults of the same trend in the Suez rift (Moustafa and Abd-Allah 1992;Moustafa et al. 1998).  Keller and Pinter (2002), which means that the watershed is susceptible to erosion (equilibrium or mature stage). However, the values of 0.52 and 0.54 obtained for Wadi Beda basin and Wadi Umm Elda subbasin indicate that they are near to the lower limit of the youth stage, in other words, they are highly susceptible to erosion.
The (B s ) values range from 1.84 to 2.89 (Table 5), which seems to have moderate tectonic activity, except for Wadi Umm Elda which is displaying relatively higher value. High values of B s are associated with elongated basins, generally associated with relatively higher tectonic activity. Low values of B s indicate a more circular-shaped basin, generally associated with low tectonic activity (e.g., El-Hamdouni et al. 2008;Buczek and Górnik 2020). The elongated drainage basins are considered to be relatively young and are commonly encountered in active tectonic areas and could be primarily due to the downcutting of the stream (Bull and McFadden 1977).
The directional set parameters were treated independently and each direction is given a definite weight value, based on its relative abundance in the study area. Four directional families are recorded, and they are NW, NNW, WNW, and NE. Furthermore, the weighting values assigned to each direction are as follow: NW (15%), NNW (15%), WNW (5%), NE (5%), NS (2.5%), NNE (2.5%), ENE (2.5%), and EW (2.5%). The sum of weighted values of the three lineament parameters, length, number and directional groups, is considered the photo-lineaments factor values that are displayed as an integrated colored map (Fig. 16). The photolineament factor map exhibits areas of high anomalies that are varying in extensions and size and are most likely underlain by highly fractured rocks. Areas of high anomalies are located at three sites, two of them locate at the upstream Wadi Umm Athala, and the third one at the northeast border   of the Wadi Beda main basin. The low anomalies over Gebel El-Ramliya, Gebel Akheider can be interpreted, in terms of the three lineament parameters, as areas characterized by short, low abundant and low-weighted lineaments. The resulted seismic intensity map of the study area is shown in Fig. 17. An examination of the map reveals the presence of large areas with high to very high seismic intensity. Gebel El-Noqra, Gebel El-Ramliya, Gebel Akheider, and Gebel Umm Zeita are the locations of these areas. The research area's reminder has low and moderate seismic intensities.

The integrated thematic map
The map assigned weights and the geomorphic potentiality values of Wadi Beda Basin are listed in Table 6. The thematic map of Wadi Beda Main Basin (Fig. 18) was produced by integrating the different geomorphic and geological parameters, based on these assigned weights and potentialities. Careful inspection of the integrated thematic map discloses that the majority of the study area exhibits a moderate-active zone (83.4%), whereas some areas display lowactive zone (11.10%), and other areas display high-active zone (5.5%) (Table 7). Specifically, three definite high-active spots could be recognized; these are.

Gebel Umm Zeita and downstream of Wadi Umm Athala 2. Gebel Akheider and downstream of Wadi Akheider 3. Gebel El-Ramliya and downstream of Wadi El-Noqra
Based on the weighted basin and channel geomorphometric parameters (circularity, basin length, basin area, basin perimeter, bifurcation ratio, drainage density, and stream frequency), the Wadi Beda Main Basin was given a moderate rank of flood potentiality (Arnous et al. 2022). This moderate rank indicates that Wadi Beda Main Basin is not capable to collect high quantity of runoff water, and hence, a low quantity of water moves out the basin, which is lowrisk Basin.
The high-active areas are characterized by the presence of E-W belt of left steeped en-echelon normal faults and WNW to NW-oriented normal faults and very high seismic intensities (Figs. 3 and 7). The E-W faults have a left-stepped en-echelon arrangement belts which are proposed to overlie E-W deep-seated faults. Moreover, these deep-seated faults have been rejuvenated by right-lateral wrenching during the early phases of opening of the Suez rift (Early Miocene)   (Moustafa and Abdallah 1991;Abdallah 1993). Aeromagnetic data confirmed that the Gulf of Suez region appears to be a small down-lifted block dragged down between the larger faults (Abuelnaga et al. 2019). It is clear that the main tectonic trends, which have affected the study area, are the NW-trending Faults and E-Wtrending faults. Hegazi et al. (2013) correlated these trends with the well-known tectonic events affecting Egypt. The NW-trending faults are related to the Gulf of Suez rifting (Oligocene -Early Miocene). Whereas the E-W-oriented faults were mostly originally formed in Late Triassic-Early Jurassic during the rifting of the African-Arabian plate away from the Eurasian plate. These faults were reactivated to deform the upper Cretaceous rocks and continued to affect the overlying Tertiary rocks. The intersections of these two trends represent sites of high clustering of epicenters. The majority of drainage lines in the Wadi Beda Main Basin are controlled by NW-trending faults, and the sites of perturbation are localized at the intersections with the E-W-trending faults, which make them high-risky zones.
According to the integrated geo-environmental hazards map of the northwestern Gulf of Suez area (Al-Moazamy 2018), the Wadi Beda Main Basin has flash flood and seismo-tectonic threats (Fig. 19). The unpredictable nature of flash floods makes detecting, analyzing and predicting them a unique problem. The best mitigation approach for flash flood threats is the construction of retardation dams, the digging of artificial basins to collect floodwater, or the

Summary and conclusion
The present study aims at investigate, analyze and integrate the geomorphic and geological parameters, and hence assess the tectonic activity of Wadi Beda Main Basin and its sub-basins; Wadi Umm Athala, Wadi El-Noqra, and Wadi Umm Elda.
The Shuttle Radar Topography Mission (SRTM) images have been used in producing DEM, contour, and relief maps and extracting drainage networks a 3D-perspective model of the study area. The geomorphic parameters include valley floor width to height (V f ) index, stream length-gradient index (S L ), transverse topographic symmetry index (T), hypsometry index (H i ), and drainage basin shape (B s ). Whereas the geological parameters are the structural lineaments, and in addition to seismicity data.
The distribution maps of V f index, S L index, factor (T), photo-lineament factor, and seismic intensity have been integrated to produce a thematic map of Wadi Beda Main Basin. This thematic map shows that the majority of the study area exhibits a moderate-active zone (83.4%), whereas some areas display low-active zone (11.10%), and other areas display high-active zone (5.5%). Specifically, three definite high-active spots could be recognized; these are Gebel Umm Zeita and downstream of Wadi Umm Athala, Gebel Akheider and downstream of Wadi Akheider and Gebel El-Ramliya and downstream of Wadi El-Noqra. These highactive areas are characterized by the presence of E-W belt of left steeped en-echelon normal faults and WNW to NWoriented normal faults and very high seismic intensities. It is clear that the main tectonic trends, which have affected the study area, are the NW-trending faults related to Gulf of Suez Rift (Oligocene -Early Miocene), and E-W-trending faults mostly originally formed in Late Triassic-Early Jurassic during the rifting of the African-Arabian plate away from the Eurasian plate. The intersections of these two trends represent sites of high clustering of epicenters and represent highly-active-zones.
Author contribution A.E.: data analysis and interpretation of results and drafting, reviewing, and submitting the original manuscript. A.M.: field data collection, analysis and interpretation of results, and sharing the first draft. M.O.: data integration and GIS modeling and sharing the first draft of manuscript. A.A.: collecting, analyzing, and interpreting field data, as well as drafting figures.
Funding Open access funding provided by The Science, Technology & Innovation Funding Authority (STDF) in cooperation with The Egyptian Knowledge Bank (EKB).
Data availability All data generated or analyzed during this study are included in this published article.

Declarations
Ethics approval and consent to participate Not applicable.

Conflict of interest The authors declare no competing interests.
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/.