Importance of Beach Evaluation for Integrated Coastal Zone Management A Case Study of Tyre Southern Lebanon

Coastal erosion and accretion are short-term and long-term processes that may cause changes in shoreline and beach profiles. Monitoring beach profile is important to figure out the factors that are responsible for these changes so it helps the stakeholders in planning sustainable beach management that plays a vital role in achieving Integrated Coastal Zone Management (ICZM). For this reason, a beach survey was conducted at the study area to monitor and determine the factors that led to erosion and accretion on beach profile. DEM (Digital Elevation Model) helped on comparing between the elevations of ground control points through four seasons of the survey, which compensated: November 2020, December 2020, February 2021, April 2021 and July 2021. The used methodology of using GPS (Global Positioning System) depended on stop and go methodology within the two locations in the study area: site A (Jal El Bahr north of Tyre) and site B (Bahr Slaiyeb south of Tyre). GIS (Geographic Information System) is another tool that was used to have LC (land cover) map. Meteorological data was extracted from Saksakye weather station north of Tyre. Results of beach survey indicated that the two study areas were exposed to change on profile through the study period especially at Bahr Slaiyeb. These changes in the two study areas are a result of anthropogenic pressure (grooming and cleaning, grazing, human trampling, construction, etc..), vegetation cover and meteorological factors (wind speed, wind direction and precipitation). This database could be used to support in outlining the baseline of strategic beach monitoring and ICZM.


Introduction
Coastal zones are among the world's most productive locations, providing a diverse range of rich habitats and ecosystem services that have drawn humans and human activities for millennia (Tissier, 2020). As a result of increased industry, trade and commerce, tourism, and the population growth and migration, as well as diminished water quality, the world's coastal zone became increasingly stressed (Nayak, 2000). Long-term management tools must be used to improve coastal resource conservation while increasing the sustainable growth of their use (Tissier, 2020).
ICZM is a widely acknowledged management framework for addressing coastal and marine environmental challenges and disputes, as well as achieving sustainable coastal resource use, especially in poor countries (Pollnac and Pomeroy, 2005). The mission of the ICZM is to accomplish sustainable development of coastal and marine areas, which includes preserving important natural processes, promoting economic development, and resolving disputes arising from numerous uses of coastal space (Meliadou et al., 2012). The Mediterranean Sea is dealing with a number of environmental issues, including habitat morphological changes, water deterioration, eutrophication, alien species introduction, shoreline erosion, and increased surface water vulnerability (Maugendre et al., 2015). Because of these realities, an integrated management strategy for this coastal zone is required (Zaucha et al., 2016). ICZM was used in the Mediterranean Sea to attain Good Environmental Status (GEnS) for its ecological systems by 2020 (Knežević and Petović, 2016).
Beaches are already threatened by a variety of human activities. Human changes to sandy shores began at least two centuries ago (Nordstrom, 2000), and are predicted to intensify exponentially over the next few decades (Brown et al., 2008). Evidence for ecological change in beach ecosystems, which are exposed to human pressures at scales and intensities unmatched in history, is accumulating worldwide (Brown and McLachlan 2002;Schlacher et al., 2007). Aside from direct anthropogenic effects on beaches, global climate change is expected to have profound, broad, and long-term ramifications for the world's marine ecosystems (Feagin et al., 2005;Harley et al., 2006). Planners and decision makers need a greater awareness of the instability of sandy beaches and the possible impact of human activities. It is vital to know the shoreline evaluation according to the linked time scale for coastal management purposes (Toure et al., 2019). Furthermore, comprehending coastal erosion necessitates a comprehension of all the elements that interact along the coast, as well as an awareness of various timelines (Marchand et al., 2011). Erosion is a phenomenon that occurs over a wide range of time scales in this context. As a result, in order to combat erosion, a comprehensive view of processes at many scales is required (Fekete et al., 2009).
On sandy coasts, many present and future conservation and mitigation approaches necessitate active threat management. The examination of morphological changes in beach-dune systems are required for the research and management of coastal regions, at least in certain sectors (Pardo-Pascual et al., 2005). Monitoring the behavior of shorelines is of considerable social and economic importance in support of setback planning, hazard zoning, erosion-accretion management, regional sediment budgets, and the establishment and validation of models for coastline changes (Zuzek et al., 2003). Shoreline change assessment is a crucial component of coastal engineering project planning and evaluating, large-scale shoreline change studies are required to determine the magnitude of the issues and important concern areas (Stauble et al., 2003). Monitoring shoreline changes aids in determining the nature and processes that caused the changes in any given location, assessing the human impact, and planning management (Jana et al., 2016). To build sustainable shoreline management plans, it is vital to monitor and understand long-term, seasonal, and short-term shoreline changes (Kankara et al., 2015).
No matter how several initiatives were done by the Lebanese government to achieve ICZM, these initiatives did not reach their aims. Lebanon still lacks an integrated plan or strategy for the administration and protection of its coastal zone. Due to some of the irreversible losses to Lebanon's coastal zone, which is considered a national asset of significant economic importance, the adoption of such a strategy/plan has grown increasingly vital and urgent by the day. Different research confirms that sandy beaches on Tyre as Lebanese beaches have been exposed to the erosion and accretion from Qasmieh till Naqoura (Saad et al., 2021). Despite this fact, there is no studies light upon the factors that lead to this phenomenon, also there is no integrated system of beach management designed to evaluate beach environments from all aspects. This study investigates the beach profile change which aims to establish a base information about human and physical factors that lead to this change. This research is an early stage of the decision making process that helps stakeholders and authorities in the planning for ICZM.

Study Area
Tyre shoreline is located to the south of the Lebanese coast. Tyre is one of the most wellknown Phoenician cities, it is one of the biggest areas of Lebanon's coastal plain, ranging in width from a few hundred meters to three kilometers (Nir, 1996).
Tyre's maritime location is noteworthy; it is located along the Mediterranean coast, approximately 80 kilometres south of Beirut and 26 kilometres north of the international border. Tyre's coastal zone is a densely populated area marked by significant levels of immigration, which has been a feature of this region since decades. A diverse spectrum of economic, commercial, and recreational activities is attracted to this location. Tyre as many Mediterranean cities plays a vital role in touristic sector. In busy summer days, Tyre's public beach can attract up to 20,000 visitors.
Two sandy beaches were selected in this study (Figure 1), the first is Jal El Bahr (Site A) located at north of Tyre city, and the second is Bahr Slaiyeb (Site B) located at the South of Tyre city. Site A is a part of Al Abbasiye and El Bourghliye Beach, where the width of the beach exceeds 50 meters, with fine sand smoothly rising up to the dunes (Kasparek, 2015). The beach profile at Jal El Bahr consists of dunes, berms, wrack line, beach face and swash zone. While a beach is connected to the sea and influenced by waves and tides, dunes are connected to the ground and influenced by wind, they are home to a wide range of highly specialized plants and animals, including endangered and rare species. Sand dunes in this site are characterized by their grey color with well-developed sand dune vegetation. The height of dunes at this site is ranged between 50 cm and 1.5 m.
The swash zone is the flow region that is continuously covered and uncovered by the run-up and backwash of water as waves impose on a beach (Kelly and Dod, 2010). Berms are typically generated between a mean sea level and a maximum wave run-up level due to these sediment transport rates (Suzuki, 2020). The material found in the swash zone, in lines along the foreshore, and occasionally on the back of the beach, especially after storms, is referred to as wrack line zone. Season, coastal landform, offshore substrates (determining algae/seagrass development), currents, tidal forces, wind, and wave action, all influenced the quantity. The beach face is the part that interacts regularly with wave run up and its shape can respond on short timescales, depending on the balance of forces moving sediment (Bujan et al., 2019).
Bahr Slaiyeb is recognized by a spectacular sandy beach, the most famous beach in Lebanon and from the vital beaches in Levantine basin. This site is located in the touristic zone of Tyre Coast Nature Reserve (TCNR). TCNR is a Ramsar site (The Ramsar Convention on Wetlands) that was ratified in 1/ 3/1999 (Law No.23) as well as a nature reserve (The Reserve is part of the city of Tyre, which was designated in 1984 by the United Nations Educational, Scientific and Cultural Organization (UNESCO) as a World Heritage Site). It's a home to a diverse range of flora, animals, and insects, in addition to providing amusement for visitors. It is a nesting habitat for the endangered Loggerhead and Green Sea Turtles, as well as a haven for the Arabian Spiny Mouse and a variety of other significant animals. A profile of Bahr Slaiyeb is formed from the dunes, berms, wrack line, beach face and swash zone as site A (Figure 1). The sand is extremely fine in texture and heavily compressed. Sand dunes are very flat, and they appear to be the remnants of a much larger and higher sand dune complex. There are dune slacks with permanent stagnant water (Kasparek, 2015). Beach are made up mainly of Biodetritiques sands and have a CaCO3 content (Sanlaville, 1977).

Real Time Kinematic (RTK) GPS surveys
Changes in beach shorelines are a complicated process, including oceanographic and coastal processes such as waves, tides, currents, and sediment movements. A range of techniques can be used to measure and develop an understanding of these processes, including the use of mathematical models, aerial photography, and topographic surveys (Binti et al., 2016). Various methods were applied to monitor the beach profile changes, including direct measurement using GPS-RTK (Global Positioning System-Real Time Kinematic) (Ruggiero et al., 2009), Total survey station (Lee et al., 2013) as well as remote sensing methods such as coastal imaging (Holman and Stanley, 2007) airborne LIDAR (light detection and ranging) (Saye et al., 2005) and Unmanned Aerial Vehicles (UAVs) for coastal surveying (Turner et al., 2016). Besides, photogrammetric techniques and terrestrial laser scanning (TLS) were used to derive digital elevation models (DEMs) which were used to improve our understanding of coastal changes at different spatial and temporal scales (James et al., 2013). The use of GPS allows the detection of elevation and geographic location. One-dimensional GPS surveys can be used to monitor the berm, high-water line, or vegetation line quickly and affordably. In geodetic and engineering surveying, relative positioning with Global Positioning System (GPS) carrier phase measurements were widely employed for precise positioning (Chen et al., 2009). Most GPS inaccuracies can be cancelled out or greatly reduced in a relative positioning mode (Olynik et al., 2002). As a result, some sources of GPS inaccuracies, such as the earth tide, phase wind-up effects, and satellite antenna offset, are not critical for relative positioning over a short baseline (Chen et al., 2009). Beach morphology monitoring using RTK GPS surveying techniques was widely accepted as an accurate and efficient means to collect coastal morphology data (Binti et al., 2016).
Due to the advantages of RTK GPS highlighted by previous researchers, the RTK GPS was selected for this study. We used GPS/GNSS (The Global Navigation Satellite System), RTK RECIVER, Full-spectrum of satellite signals, Full connectivity. Kinematic GPS Surveys involve collecting data while the antenna of the mobile receiver is moving. The movement of the antenna can be intermittent: stop and go, or continuous: fully kinematic. We relied on Stop and go methods represented in Figure 2. Typically, a few minutes of data collection at each point were sufficient. A survey that had started from November 2020 till July 2021 as mentioned in Table 1, covers autumn, winter, spring and summer seasons. The user of stop-and-go kinematic surveys are usually on land and looking for a quick survey of a succession of stationary spots (Morton et al., 1993). During the survey, we linked the GPS receiver rover to GPS receiver station set up on (SRIFA STATION) that has stereographic projection. This receiver base measures its coordinates according to satellites and compare it with entered coordinate to calculate the error in the observation. Basing to the measured error, GPS base delivers to GPS rover corrections, and the vertical accuracy of the GPS rover was about ± 3cm. Noticing that many obstacles have been faced during this survey. The first was the lockdown due to Corona Virus, the second was the Lebanese fuel crisis. These issues hindered our track in complete all months of the seasons.

Sections
We divided each site into different sections. At Jal El Bahr beach, we created three sections due to the variety of topographic features on this site. Also Site B was divided into two sections. Survey on different sections was launched from sand dunes till the swash zone. Every coordinate (x,y) represented by code number from 1 to 20 for ex: x= -366266.81 and Y= -94569.11 is a specific ground control point.

Data acquisition
One Landsat image was downloaded from the USGS Earth Explorer (https://earthexplorer.usgs.gov/) for the year of 2021 to the study area, the detailed data are shown in Table 2. This data was used for the performing of LC maps via the GIS (Geographic Information System).

Image classification
Classification errors are widely reported using univariate statistical metrics such as producer, user, and overall accuracies, as well as omission and commission errors. Maximum likelihood classification was chosen to detect a two class (vegetation cover and sand area).
Furthermore, the 'Kappa Coefficient Index,' a discrete multivariate statistical model, has become a commonly used standard tool for CAA (classification accuracy assessment) processes (Jensen, 2005;Congalton and Green, 2009). Producer's Accuracy (PA) is the map accuracy from the point of view of the map producer. It is the number of reference sites classified accurately divided by the total number of reference sites for that class. The minimum acceptable accuracy for a class is 90%. The User's Accuracy (UA) is the accuracy of the attitude of a map user. The User's Accuracy is calculated by taking the total number of truthful classifications for a particular class and dividing it by the row total (Story and Congalton, 1986). The proportion of correctly identified samples in relation to the total number of samples on the map is referred to as overall accuracy. The overall accuracy percentage was calculated as 100 x (Total number of correct samples) / (Total number of samples).
The Kappa coefficient (Cohen, 1960(Cohen, , 1968 is commonly used to evaluate picture classification accuracy. It is mostly based on the difference between the amount of agreement shown and the amount of the agreement that would be expected to occur by chance alone (Thakkar et al., 2017). The formula for computing Cohen's k is expressed as (P0 -Pe) / (1-Pe) where P0 is percent agreement, defined as the proportion of subjects on which the raters agree, and Pe is chance agreement, defined as the proportion of agreement that would be expected by chance.

Meteorological parameters
Meteorological data has been extracted for Tyre coastline during the study period that extended from 8 November 2020 to 25 July 2021, from Saksakiye weather station, it is located north Tyre. The parameters are: high temperature (High Temp), Low temperature (Low Temp), precipitation (Precip), high wind speed (High WS), and wind direction.

Sections
The length of sections at site A are different, section one is 143.5 m, section two is 139.4 m while section three is 131.7 m, the distance between sections one and two is 29.1 m while between two and three is 18.3 m (Figure 3). The length of section one at site B is 199.9 m and section two is 183.5 m, the distance between sections one and two is 90.1 m. This difference in length is due to different topography deviation during survey (Figure 4).

Site A (Jal El Bahr)
DEMs are used to assess coastal processes and coastal morphology changes over various temporal and spatial scales (Li and King, 2007;White and Wang, 2003) and can be used to support coastal management (Davidson et al., 2007). Section one ( Figure 5) indicates a different elevation of ground control points through same month, for example, on November ground control point 2 registered 23.471 m, then it declined 21.115 m at ground control point 16. Depending on this figure, we noticed small modification on the profile of this section between the periods of the study. Section two and section three (Figures 6 and 7) also showed that the elevations of each ground control point registers different values through the same month and during the months. All profiles showed that there were few erosion and accretion processes taking place on site A. The survey recorded high elevation at dunes then it declined at beach face and swash zone.

Site B (Bahr Slaiyeb)
Section one ( Figure 8) indicated high variations between control points through the same month and between the months and periods. In November, the elevation was 22.518 m at control point 1 then it increased to reach 23.069 m at control point 11 and 23.544 m at control point 12, then it decreased to reach 20.477 m at control point 19. We noticed also vital variations in the elevation of ground control points through months. For example, at control point 11 the elevation was 23.06 m in November then it declined to 22.962 m in December then it increased in February to reach 23.027 m. After that the elevation decreased to reach 22.914 m and 21.977 m in April and July respectively. High change was recorded during July as represented in the graph during the period of (April-July), the curve of July declined in obvious rate which led to dramatically change of section one, profile B. Figure 9 showed that the profile of section two site B was exposed to a slight variation, especially in December as represented on graph one (November-December). Beach survey recorded different elevations within ground control points. The elevation was high at flat dunes, then it declined to reach its minimum values at swash zone.

Vegetation cover
Sediment distribution is the outcome of the dune environment's retention capacity, which is influenced by a variety of elements including wind regimes, runoff, coastal currents, wave energy, vegetation coverage, and coastal exploitation (Valentini et al., 2020). Beach erosion is largely prevented by vegetation, plant roots in the sand serve to keep the sand in place and make dunes more stable (Robertson, 2010). Plants contribute to build dunes and improve geomorphological stability, while dunes provide habitat for them. Disturbances related to sand movement, whether in accretion (dune formation) or erosion, have a significant impact on these interactions (Miller, 2015).   The vegetation cover is weak in Bahr Slaiyeb. Vegetation cover at this area includes: European Marram Grass, Phragmites australis, Ammophila arenaria (Kasparek, 2015). This vegetation cover is found on the flat dunes noticing that the bloom of plants differs from season to another depending on the weather and meteorological condition. The two maps represented a green color within the wrackline and this is due to the deposition of organic materials (seagrass, shell etc.) by waves. The overall accuracy of LC is 99.3% of site B (Table 4). The Cohen's k is equal 0.92 for site A and 0.98 for site B. Landis and Koch (1977) proposed benchmarks for interpreting Cohen's k, with 0.8 to 1.0 indicating nearly perfect agreement, 0.6 to 0.8 indicating substantial agreement, 0.4 to 0.6 indicating moderate agreement, 0.2 to 0.4 indicating fair agreement, zero to 0.2 indicating slight agreement, and zero or lower indicating poor agreement.

Wind and precipitation
Depending on Table 5, the temperature recorded moderate level during the seasons of the study period. The maximum average temperature was recorded during summer season 26.82°C while the minimum average recorded in winter season 16.0°C.
Maximum precipitation was recorded in the period that extended from 14 December to 7 December 2021 which registered 120.90 mm. The minimum value of precipitation was recorded in the summer within a period that extended from 5 July 2021 to 25 July 2021 (0 mm).
Maximum high wind speed recorded in a period of (April 5, 2021-April 11, 2021) which registered 11.3 m/s while the minimum wind speed recorded in a period of (July 19, 2021-July 25, 2021) which registered 5 m/s. Different wind direction was recorded during the period of the survey. The direction of the wind was SW in autumn, SSE and S through winter, WSW and NW in April while the wind direction was WSW and W in summer season.

Discussion
The beach profile changes for five sections have been successfully obtained. The research and management of coastal environments, at least in some sectors, requires the investigation of morphological changes in beach-dune systems caused by specific agents such as waves, wind, and human activities, as well as antecedent topography, vegetation, material texture, and so on (Pardo-Pascual et al., 2005).

Site A: Jal El Bahr
Depending on Figure 11, both erosion and accretion were witnessed at section one through the periods. While the accretion is more than erosion at section two and three during the period of November-December and December-February. The value of accretion reached 0.264 m at the ground control point 10 as a maximum value at section one and 0.13 m at section three, ground control point 12. The moderate level of rainfall that ranged between 102.6 mm during November, 150 mm during December and 19 mm during February helped the dunes to keep their stability with the support of vegetation cover that acted as a natural obstruction against winds and protect the dune from a high erosion. Collectively Sections one, two and three observed a rate of erosion during April, we can notice that during the period of February-April the rate of erosion at section two reached -0.184 m (18.4 cm) as a maximum level at ground control point 6 and -0.048 m (-4.8 cm) as minimum at ground control point 3. While the rate of erosion at the same period in section three reached -0.371 m (37.1 cm) as a maximum level at ground control point 13 and -0.026 m (-2.6 cm) as minimum value at ground control point 10. This change in the elevation of sand was happening as a result of rainfall, wind speed and wind direction. During July, the three sections were exposed to accretion at the dunes, beach face and swash zone as a result of stable weather and vegetation cover. Jal El Bahr was exposed to human pressure as a result of many activities like grazing and burning. Grazing included the removal of biomass, reduction fast-growing competitive species (Kooijman, 2004). Livestock may affect both sedimentation and marsh accretion negatively (e.g. Leonard, 1997). Furthermore, livestock caused compaction of the soil by trampling and hence reduced surface elevation change (Elschot et al., 2013;Pagès et al., 2019).
Off-road vehicles (ORVs) were witnessed at Jal El Bahr beach and this has affected also the elevation of sand. ORVs are commonly used on beaches and dunes worldwide (Schlacher et al., 2007) and cause damage that includes: disturbing the physical attributes and stability of dunes and beaches (Priskin, 2003) destroying dune vegetation, leading to lower diversity and less floral ground cover (Groom et al., 2007); and disturbing, injuring or killing beach fauna (Schlacher et al., 2007). These activities will lead gradually to change the profile of Jal El Bahr beach. About 50 touristic tents are installed on site B (Bahr Slaiyeb beach) every summer to receive visitors from Tyre and around. This huge flat beach is most likely the remains of sand dunes which sand was extracted for construction purposes. Human impact (constructions, human trampling, camps, cleaning and grooming, mining) and physical factors at this site led to modification of the beach profile. Referring to Figure 12 (section one), erosion and accretion happened obviously. There is an accretion that occurred between November and December.

Site B: Bahr Slaiyeb
The direction of the wind in relation to the shoreline has a significant impact on sediment delivery to the dune, changing apparent fetch lengths and foredune face slopes, and hence the location and amplitude of deposition (Cohn et al., 2019). During February (Figures 8 and 9) we witnessed an erosion in both sections and this was due to the S and WNS wind (alongshore and offshore transport) that reached 9.4 m/s noticing that the precipitation was 18 mm and this led to erode the sand from the flat dunes to carry it into swash zone. During the period of (December-February) the rate of accretion in the swash zone on section one (ground control point 18 and 19) ranged between 0.094 m (9.4 cm) and 0.095m (9.5 cm) while the rate of accretion on section two (ground control point 15 and 16) at swash zone was 0.116m (11.6 cm) and 0.153m (15.3 cm). Throughout April the value of precipitation was about 82.1 mm and the speed of offshore wind ranged between 9.4 m/s and 11.1 m/s, then this led to decrease the elevation of sand at the two sections, so we noticed erosion at the period of February-April.
Cleaning or grooming is a common practice used on tourist beaches causing heavy pressure (Dugan et al., 2003;Davenport and Davenport 2006), this act involves raking and sieving the sand, often with heavy equipment, to remove macrophyte wrack, litter, and other debris (Kinzelman et al., 2003;Fanini et al., 2005). During February 2021, tar disposal from Israel washed up Tyre beach and this pushed different campaigns to remove tar from the beach by using different kinds of tools like sieves and shovels, and this also led to changes in the elevation of sands. Tyre's public beach hosts approximately 375,000 visitors a monthly average during summer season. Trampling has a seasonal influence on supralittoral fauna at medium-high latitudes, allowing for some recovery during low-use seasons (see e.g., Scapini et al., 2005). Human trampling, tents and camps not only lead to harm fauna and flora but also to affect the elevations of sand.
During summer, the power of waves is low and may reach beach face while in winter the power of waves is large and has high energy moving offshore. This led to having accretion during July (Figures 8 and 9) within the two sections except control point 12 and 13 on section one. The rate of erosion at ground control 12 recorded -2.026 m (the ground control point was 23.544 m during December then it declines to 21.518 m during July) while the rate of erosion at ground control point 13 registered -1.12 m (the ground control point was 22.553 m during December then it declines to 21.441 m during July). This due to the process done during May, where Tyre's municipality extracted sand from section one to fill the water stream that excavated it between section one and section two at the beginning of winter (during November) in order to direct the water that has been collected in the parking towards the sea by creating artificial water path as mentioned on Figure 13). Noticing that the summit on ground control points 10 and 15 is a result of the gliding sand caused by the Bulldozer. Figure 13: Location of parking, water path, section one and section two at Bahr Slaiyeb, derived from Google Earth, 12,2020 Jal El Bar ( Figure 11) as Bahr Slaiyeb (Figure 12), recognized lessen in elevations at swash zone, wrack line and beach face, this was mainly due to longshore, cross-shore, swash and backwash sediment transport that were created by waves and currents. The swash zone is critical to the nearshore sediment budget because it is where shoreline erosion and accretion occur (Puleo et al., 2000). Because of local hydrodynamic pressure and cross-shore and alongshore sediment transport gradients, the sediment bed is active (Lanckriet et al., 2014).

Conclusion
Natural and man-made sources have been identified as the primary causes of ongoing erosion and accretion processes within the two sites of the study area: Jal El Bahr and Slaiyeb Beach. This paper presented an example of how the relative kinematic GPS method can be used to monitor the advance or retreat of coastal areas. Analysis of DEM created the lines of beach profile effectivity, the erosion and accretion processes determined that they have occurred throughout the research period. Despite the fact that the changes may not be noticeable, it is still possible to determine that the coastline exhibits were difference in beach profiles each time data was collected. Essentially, it conveyed the impression that the beach is constantly changing over time. Wind speed, wind direction and precipitation played vital role in the change of elevations through the study period. Human intervention and anthropogenic pressure were found within the two study areas, especially at Bahr Slaiyeb, so it can be predicted that this open touristic zone will always change from month to month, from year to year and this is due to the effects of unrestricted human activities. Vegetation on Jal El Bahr beach makes the dunes more stable than the flat dunes at the Bahr Slaiyeb beach.
Geodatabase represented by physical and human factors helps in perfect understanding of the rate of response of beach profile to these factors. This short-term evaluation of Tyre beach aids in constructing sustainable shoreline management plans.

Conflicts of interest
The authors declare no conflict of interest.