Land Suitability Analysis for Sustainable Urban Development: A Case of Nabatiyeh Region in Lebanon

: Urbanization has led to the degradation of green areas and natural resources, which are critical for preserving biodiversity and natural beneﬁts. Sustainable urban development (SUD) practices aim to balance urbanization with preserving and protecting green land and natural resources to achieve greater sustainability. This study evaluated land suitability for SUD in the Nabatiyeh area of Lebanon using a multi-criteria decision analysis technique, the analytical hierarchy process, combined with a geographic information system. The elevations and slopes were the most important criteria, followed by the distances from urban agglomeration and industrial areas, accessibility, and land cover types, with corresponding weights of 34% and 36%, 12%, 9%, 6%, and 3%, respectively. The study created a suitability map for future SUD. This identiﬁcation assists urban planners and policymakers in making informed decisions regarding the development of urban areas. The study’s ﬁndings highlight the efﬁciency of combining multi-criteria decision analysis with geospatial techniques for urban planning and environmental management in promoting sustainable development and protecting the environment. Future research could explore the possibility of integrating administrative reforms to promote SUD, including improving governance, enhancing the efﬁciency of governmental agencies responsible for urban development, building scientiﬁc capacities, and increasing ofﬁcials’ accountability.


Introduction
Currently, urbanization has reached a high level, and more than half of the world's population lives in urban agglomeration [1]. Generally, this is due to demographic growth and the socio-economic reasons that drive people to move from rural areas to urban agglomeration [2]. Consequently, this accelerated trend towards urbanization leads to serious environmental and ecological problems [2,3], such as air pollution, degradation of green areas, exhaustive use of natural resources, fragmentation of natural spaces, and biodiversity loss [4]. To face these problems, researchers develop several methods to regulate urban development to face unplanned urban growth and maintain sustainability. Thus, the main objective is explained as determining the best geographical locations for urban development and the locations of natural areas that should be preserved [5]. Therefore, effective land cover planning is essential for the identification of the best locations for sustainable urban growth [6,7]. Then, this paper contributes to the emerging literature by taking into account sustainability as an essential issue to be integrated into the land suitability analysis. This issue could represent important guidance for decision-makers, especially when enacting urban masterplans.
This study aims to identify the land suitability for sustainable urban development in the Nabatiyeh region in Lebanon based on an MCDA methodology. The methodology combines the AHP with remote sensing data in a GIS environment. The outputs of this study are important since they could provide insight into the sustainable urban planning procedures in terms of pre-determining new locations for future sustainable urban development in the study area that witnessed the phenomenon of rapid and fragmented urban sprawl. The results will be in the form of a map presenting the most and least-suitable lands for sustainable urban development. This study also aims to show the significance of GIS and remote sensing applications in spatial decision-making for sustainable city development.
It is hypothesized that the availability of infrastructure and amenities and the accessibility to urban agglomerations and activities prevail over other factors affecting the sustainable development of new urban agglomeration. This hypothesis will be tested for validity within the analysis.
With regard to limitations, additional socio-economic factors, such as detailed population and employment densities, were not included in the study since they were not provided by official bodies due to the lack of information/census updates.
The paper is organized as follows: The Section 2 highlights the selection of the study area, the research methodology, and the data collection procedure. The Section 3 presents and discusses the results. In conclusion, the Section 4 summarizes the main findings of this study and presents the proposed topics for future research.

Study Area
In Lebanon, different districts, including the Nabatiyeh region, suffer from unplanned rapid urbanization and have experienced strong land cover dynamism [28]. In a more detailed way, the accelerated unplanned urbanization of the country, and especially in the southern part, including the Nabatiyeh region and its surrounding areas, occurred mainly in the period extending between the years 2005-2013 [29]. This specific period of accelerated and uncontrolled urban growth was spotted due to different political, security, economic, and immigration events of the country, mainly as (1) the Israeli withdrawal from the southern Lebanese region, (2) the conflict with Israel in 2006 and the post-Lebanon-Israel war reconstruction, (3) the influx of money deposits to Lebanese banks in 2008, and (4) the Syrian war and the immigration of Syrian refugees toward Lebanon as indicated by Al-Shaar and Bonin (2021) [30], in addition to limited planning regulations, and a lack of coordination between local and national governance [31]. In this context, it is noteworthy to consider the Israeli withdrawal from most southern Lebanese areas, as these areas returned to the control of the Lebanese state, which enabled the reconstruction and urbanization.
This rapid, extensive, and unplanned urbanization was behind the study area selection. Moreover, few studies examined the urban development in the country of the study area, but none of them discussed the land suitability for sustainable urban development. [32][33][34][35][36]. For instance, El Asmar et al. (2012) and Nader et al. (2008) presented, in their studies, a vision for sustainable urban planning and development in Lebanon that emphasizes community participation, sustainability, and stakeholder collaboration with the help of semi-structured interviews, municipality archive surveys, and open-ended questionnaires [32,33]. Pierre et al. (2013) discussed the difficulties Nabatiyeh, Lebanon, encountered with urban planning and growth both during and after times of violence. The authors highlighted an approach to urban planning and development that is comprehensive, open to participation, and encourages sustainable practices [34]. Tarraf-Najib (2001) investigated the municipality plans and policies for addressing topics, such as urban planning, accessibility to services, and economic development, in the Nabatiyeh area [35]. Makki Haidar (2005) investigated the elements driving rural urbanization, including demographic and economic change and how these factors affect the social and natural environment in Nabatiyeh [36]. The latter studies focus primarily on the relationship between urban development and various factors in Nabatiyeh without discussing the land suitability for sustainable urban development.
The said area comprises three urban villages: Nabatiyeh ElFaouka, Nabatiyeh ElTahata, and Kfar Remmane, with a total area of 29.76 km 2 ( Figure 1). The topographic form of the Nabatiyeh district is characterized by the presence of numerous hills (ranging in altitude from 300 to 950 m).
lighted an approach to urban planning and development that is comprehensive, open to participation, and encourages sustainable practices [34]. Tarraf-Najib (2001) investigated the municipality plans and policies for addressing topics, such as urban planning, accessibility to services, and economic development, in the Nabatiyeh area [35]. Makki Haidar (2005) investigated the elements driving rural urbanization, including demographic and economic change and how these factors affect the social and natural environment in Nabatiyeh [36]. The latter studies focus primarily on the relationship between urban development and various factors in Nabatiyeh without discussing the land suitability for sustainable urban development.
The said area comprises three urban villages: Nabatiyeh ElFaouka, Nabatiyeh ElTahata, and Kfar Remmane, with a total area of 29.76 km 2 ( Figure 1). The topographic form of the Nabatiyeh district is characterized by the presence of numerous hills (ranging in altitude from 300 to 950 m).
The study area is characterized by the presence of various land cover types with predominantly artificial (urban) and agricultural areas.

Developed Methodology
The aim of this study is depicted as conducting a land suitability analysis to define the most suitable sites for future sustainable urban growth. To this end, the adopted technical approach consists of (a) defining the criteria as the driving factors for sustainable urban development, (b) applying the AHP to weight criteria, (c) according to experts' judgments, dividing the maps of each criterion to five classes: highest, high, moderate, low, and lowest suitability, and (d) using both the weights and classes of these criteria to generate a summary weighted overlay map (WOL) showing the overall land suitability.
Sustainability requires considering socio-economic and environmental aspects while choosing well-suited locations for future urban growth. Thus, the land suitability for SUD depends on several factors. For instance, topographic factors, such as slope and elevation, The study area is characterized by the presence of various land cover types with predominantly artificial (urban) and agricultural areas.

Developed Methodology
The aim of this study is depicted as conducting a land suitability analysis to define the most suitable sites for future sustainable urban growth. To this end, the adopted technical approach consists of (a) defining the criteria as the driving factors for sustainable urban development, (b) applying the AHP to weight criteria, (c) according to experts' judgments, dividing the maps of each criterion to five classes: highest, high, moderate, low, and lowest suitability, and (d) using both the weights and classes of these criteria to generate a summary weighted overlay map (WOL) showing the overall land suitability.
Sustainability requires considering socio-economic and environmental aspects while choosing well-suited locations for future urban growth. Thus, the land suitability for SUD depends on several factors. For instance, topographic factors, such as slope and elevation, were found as some of the main driving factors for SUD in cities [19,37]. The accessibility to social activities and economic opportunities, represented by the distance from residential, commercial, and industrial areas, was also found as an essential element in defining the suitable location for SUD [13,38,39]. This requirement is highly provided by compact urban fabrics.
The provision of infrastructural utilities, such as road networks, also represents one of the main factors that affect the location of SUD [19,39].
Additionally, the preservation of natural resources is a crucial issue that should be highly taken into account when planning for sustainable development. Since the land cover distribution affects the environment and ecological protection, determining (a) suitable areas for sustainable urban development and (b) natural spaces that need to be preserved must be carried out [40,41]. Based on this literature review, the criteria used in this study are based on the following factors: (a) slope, (b) elevation, (c) road networks, (d) urban agglomeration, and (e) land cover.
The following graph ( Figure 2) presents the methodology processes to generate the final spatial suitability map for sustainable urban development.
J 2023, 6, FOR PEER REVIEW 5 were found as some of the main driving factors for SUD in cities [19,37]. The accessibility to social activities and economic opportunities, represented by the distance from residential, commercial, and industrial areas, was also found as an essential element in defining the suitable location for SUD [13,38,39]. This requirement is highly provided by compact urban fabrics. The provision of infrastructural utilities, such as road networks, also represents one of the main factors that affect the location of SUD [19,39].
Additionally, the preservation of natural resources is a crucial issue that should be highly taken into account when planning for sustainable development. Since the land cover distribution affects the environment and ecological protection, determining (a) suitable areas for sustainable urban development and (b) natural spaces that need to be preserved must be carried out [40,41]. Based on this literature review, the criteria used in this study are based on the following factors: (a) slope, (b) elevation, (c) road networks, (d) urban agglomeration, and (e) land cover.
The following graph ( Figure 2) presents the methodology processes to generate the final spatial suitability map for sustainable urban development.

Data Collection
The data collection process for this study involved the collection and processing of various types of data from different sources. In particular, the land use/land cover map (LULC) was generated by the National Center for Remote Sensing in Lebanon (NCRS-L), using Spot satellite imagery with a spatial resolution of 1.5 m [42]. The LULC map for the year 2017 was prepared by applying a supervised maximum likelihood classification with reference to the Corine classification system [43]. Visual interpretation was employed to check the accuracy of the classification process. The results show an accuracy of 81% and a kappa value of 0.72 [42]. The LULC is used as a criterion for the MCDA method. The rasterized map for LULC was resampled to 10 m pixel size in ArcGIS.
The digital elevation model (DEM) with a spatial resolution of 30 m was obtained from the USGS Earth Explorer website [44]. This DEM was obtained from the Shuttle Radar Topography Mission (SRTM01) and was utilized to generate the elevation and slope maps using ArcMap 10.8 software. The DEM raster was then resampled to 10 m resolution.
Moreover, road network data were obtained from the OpenStreetMap website [45], and the road network shapefile was extracted to create a Euclidean distance map from the major and local roads. These maps are in the form of raster with 10 m resolution.

Determining Criteria Weights Using the Analytical Hierarchy Process (AHP)
In this phase, the weight ranking of the proposed criteria, based on the AHP, is presented. Defining the pairwise comparison matrix, among the criteria factors, is a prerequisite for the AHP.
The pairwise comparison matrix was defined based on the expert's judgments, where the experts conduct a comparison among all factors (one factor versus another one) according to their relative importance based on the Likert scale from 1 to 9 (please see Appendix A: Table A1). Table 1 presents the criteria for a pairwise comparison matrix, defined according to experts' judgments and statements via interviews [46][47][48][49][50]. Distance from Industrial and commercial areas 0.17 0.14 0.33 1 3 6 Distance from major roads 0.17 0.14 0.33 0.33 1 4 Land cover 0.14 0.14 0.17 0.17 0. 25 1 Based on the pairwise comparison matrix, the calculation of the final weights is conducted following the AHP process steps presented in Appendix A: Figure A1. Table 2 shows the weights of the used criteria. The consistency of the pairwise comparison matrix and the calculated weights argues that the consistency ratio (CR) must be less than 0.1 [51].
Referring to the calculations, consistency is ensured: where RI is the Random Consistency Index, in which the value is given as 1.24 depending on the number of factors [52].

Defining and Classifying Criteria Classes
Sustainability requires taking the socio-economic and environmental aspects into account while identifying suitable locations for future urban growth [53]. In this study, five criteria were selected for evaluating land suitability. These criteria include slopes, elevations, road networks, urban agglomerations, and land cover. As a standardization system for all data, the map of each criterion was reclassified according to a scoring system from 1 to 5, with 1 representing the least suitable and 5 being the most suitable. Defining and classifying the criteria were based on both literature studies and expert judgments, which were obtained via interviews. Hence, several interviews were conducted with (a) real estate developers and contractors [48][49][50]54], (b) urban planning specialists and researchers [46,47], (c) engineers [55][56][57], and (d) local authority members: the mayor of Nabatiyeh region [58] and municipality members [59,60]. Concerning the topographic characteristics, such as slopes and elevations, the urban area should be located in fairly flat lands since landslides may potentially occur more in locations with steep slopes [61]. From an economic point of view, steep slopes are unfavorable for construction because of the elevated costs. In Lebanon, the majority of urban development occurred in the coastal areas in different cities where the elevation is less than 400 m [62]. Real estate developers revealed that they prefer constructing and developing new urban zones in locations in proximity to coastal lines. These areas represent the contiguous expansion of previous urban agglomerations in Lebanon [49,50], and this contiguity to urban agglomerations makes them more attractive for new residents and for companies.
These preferences would, to some extent, (a) ensure the preservation of natural and green areas occurring in higher zones, (b) prevent unorganized urban sprawl, and (c) promote urban sustainability [31]. The urban planning specialists [46,47] and engineers [55][56][57] stated that areas below 400 m are the most suitable areas for urban development in terms of social, environmental, and economic considerations. Additionally, the majority of the Lebanese cities are located near the coastal line, such as the cities of Beirut, Tripoli, Saida, and Tyre. Hence, the closer urbanizations are to coastal areas, the higher accessibility to (a) social amenities, such as health, educational, and entertainment services, and to (b) economic opportunities, such as commercial activities and employment, is provided. Areas with elevations less than 900 m are divided into four classes: (a) Highest suitability for areas with elevations ranging between 0 and 400 m, (b) High suitability for areas with elevations ranging between 400 and 600 m, (c) Moderate suitability for areas with elevations between 600 and 800 m, and (d) Low suitability for areas with elevations between 800 and 900 m [50,58]. Furthermore, the developers [55,56] revealed that areas with elevations more than 900 m are considered the lowest suitable for urban development. The classification of the slope map suggests that areas with slopes over 20% are not suitable for urban growth, as indicated by the experts in urban planning [46,47]; this is in the same line with previous literature findings that reveal that areas with a slope less than 20% are more suitable for urban development [37,63]. Areas with slopes less than 20% are classified into four classes: (a) Highest suitability for slopes ranging between 0 and 5%, (b) High suitability for slopes ranging between 5 and 10%, (c) Moderate suitability for slopes between 10 and 15%, and (d) Low suitability for slopes between 15-20% [48,50]. The provision of major roads and infrastructure contributes to enhancing the accessibility and connectivity between neighboring cities, providing social and economic opportunities for sustainable urban growth [64]. From an economic point of view, and for the sake of sustainable urban development, attempts must encourage urban growth closer to an existing road to reduce the additional cost of new roads. To retain effective connectivity between existing and future built-up areas in the study area, zones located within a distance less than 800 m from major and local roads are considered more suitable for the SUD [46,47,56]. Areas located at distances of more than 800 m are considered the least suitable for sustainable urban development. According to the expert judgments, these areas could be further divided into four classes: (a) Highest suitability for zones located within a distance less than 200 m, (b) High suitability for zones located between 200 and 400 m from roads, (c) Moderate suitability for zones located between 400 and 600 m from roads, and (d) Low suitability for zones located between 600 and 800 m from roads. During the planning phase of new urban development in a specific region, it is essential to consider the connectivity and compactness of different built-up areas. This can help (a) limit the urban sprawl, which is the uncontrolled dispersion of development within rural areas, and (b) benefit from the existing infrastructure and utilities, such as roads, water, drainage, sewerage and electrical networks, and other infrastructural systems [13]. This could reduce the environmental impacts of new developments and make it more sustainable in the long run [46,59]. Additionally, better connectivity can also promote economic growth and improve the quality of life for residents by providing easier access to services and amenities across the entire urban area [47,58]. A compact and well-planned urban area can provide a range of benefits for residents and help overcome the challenges of building sustainable cities [65,66]. In general, urban agglomerations refer to geographic areas consisting of densely populated urban centers and their surrounding suburbs. The national center for remote sensing in Lebanon has defined and divided the urban agglomeration of the study area into two categories: condensed and uncondensed urban agglomeration referring to high-and low-density urban agglomeration [37], as indicated in Figure 3. Jaffal (2022) and Harb (2022) [57,59] specified that new urban agglomeration should be located within a distance of less than 700 m from the closest urban agglomerations to ensure the sustainable aspect of new development. The experts indicated that areas located at distances of more than 700 m are considered the lowest suitable for sustainable urban development. Therefore, areas located within 700 m of urban agglomerations are divided into four classes: (a) Highest suitability for areas located within 100 and 300 m of urban agglomerations, (b) High suitability for areas located within 300 and 500 m of urban agglomerations, (c) Moderate suitability for areas located in distances less than 100 to urban agglomerations, and (d) Low suitability for areas located within 500 and 700 m to urban agglomerations. While certain industrial and commercial activities can generate noise, traffic, and pollution that may be disruptive to residential areas, the main activities of the study area are characterized by light industrial and commercial activities, which do not generate significant noise and pollution and do not also have a negative traffic impact. The industrial While certain industrial and commercial activities can generate noise, traffic, and pollution that may be disruptive to residential areas, the main activities of the study area are characterized by light industrial and commercial activities, which do not generate significant noise and pollution and do not also have a negative traffic impact. The industrial activities are limited to the food industry, textile and furniture manufacturing, and other unlisted industrial services. Commercial activities are only limited to retail services and do not encompass extensive urban logistics.
For instance, having commercial and industrial areas nearby to urban agglomerations is essential for accessibility to employment and opportunities, and to support local businesses, and improve access to goods and services. Additionally, having a mix of land use can create vibrant and vital neighborhoods [59]. Koubeissy (2022) and Ismail [47,58] stated that the most suitable areas for SUD are at a distance less than 900 m from the industrial and commercial areas. These areas are further classified into four classes: (a) Highest suitability for areas located with a buffer distance between 300 and 500 m to urban agglomeration, (b) High suitability for areas located between 500 and 700 m to urban agglomeration, (c) Moderate suitability for areas located in distances less 300 m away from urban agglomeration, and (d) Low suitability for areas located with distances ranging between 700 and 900 m.
The preservation and protection of ecosystems and natural resources are crucial for sustainable urban development. Uncontrolled urbanization and exhaustive development have depleted natural resources, pollution, and environmental degradation. Therefore, it is vital to ensure that future urban development is planned in a sustainable way that preserves the environment and natural resources. Hence, forests, green, and natural areas should be preserved and protected [13,67]. The LULC classification is reclassified depending on the LULC type. Thus, (a) High suitability (value 5) is assigned to unproductive bare soil and grass areas, (b) Moderate suitability (value 3) is assigned to the agricultural areas, and (c) Low suitability (value 1) the wooded and built-up areas [46,47].
The defined criteria classes are presented in Figure 4.    The final phase of the analysis is the generation of the suitability map. This process could be conducted by employing the weighted overlay (WOL) tool in the ArcMap software. This process consists of combining the data of (a) Weights of used criteria and (b) suitability classification maps of these criteria. These weights were derived using the AHP method and are presented in Table 2. Elevation and slope have the most important weights, with corresponding values of 0.34 and 0.36, respectively. The distance from urban agglomeration has moderate importance with a weight value of 0.12.
The least important criteria are the distance from industrial and commercial areas, the distance from main major and local roads, and the land cover types with weight values of 0.09, 0.06, and 0.03, respectively. These weights were then used to combine input classified maps of criteria into a single output layer representing the suitability map using the weighted overlay tool, which is included in ArcGIS software. The calculation consists of (a) multiplying the normalized values of each criterion layer (raster maps) by their respective weights, (b) summing up the resulting values of raster pixels, and (c) rounding the obtained values to have only integer ones.

Results and Discussion
The generated suitability map, presented in Figure 5a, shows the suitability degrees for the SUD. To highlight the suitable/non-suitable areas for future SUD, the existing built-up areas were retained from the analysis, as indicated in Figure 5b.  The final phase of the analysis is the generation of the suitability map. This process could be conducted by employing the weighted overlay (WOL) tool in the ArcMap software. This process consists of combining the data of (a) Weights of used criteria and (b) suitability classification maps of these criteria. These weights were derived using the AHP method and are presented in Table 2. Elevation and slope have the most important weights, with corresponding values of 0.34 and 0.36, respectively. The distance from urban agglomeration has moderate importance with a weight value of 0.12.
The least important criteria are the distance from industrial and commercial areas, the distance from main major and local roads, and the land cover types with weight values of 0.09, 0.06, and 0.03, respectively. These weights were then used to combine input classified maps of criteria into a single output layer representing the suitability map using the weighted overlay tool, which is included in ArcGIS software. The calculation consists of (a) multiplying the normalized values of each criterion layer (raster maps) by their respective weights, (b) summing up the resulting values of raster pixels, and (c) rounding the obtained values to have only integer ones.

Results and Discussion
The generated suitability map, presented in Figure 5a, shows the suitability degrees for the SUD. To highlight the suitable/non-suitable areas for future SUD, the existing builtup areas were retained from the analysis, as indicated in Figure 5b. According to the findings presented in Figure 5b, zones with the highest degree of suitability are covering a small area of 0.11 km 2 , representing only 0.56% of the remaining unbuilt zones of the study area. Conversely, areas with high suitability extend over a larger area of 9.81 km 2 , representing 50.19% of the undeveloped lands. Areas with moderate suitability reach over 6.85 km 2 or 29.38% of the undeveloped lands, whereas low suitability areas cover 0.63 km 2 , accounting for 17.87% of undeveloped lands. The areas with the lowest degree of suitability represent only 1.99% of the undeveloped areas, with a total area of 2.87 km 2 . Consequently, the study implies that future development projects could be implemented, generally, in areas with high and moderate suitability degrees for SUD since these zones have larger spatial extent. Figure 6 shows the percentages of suitability degrees relative to the non-built areas in the case study. According to the findings presented in Figure 5b, zones with the highest degree of suitability are covering a small area of 0.11 km 2 , representing only 0.56% of the remaining unbuilt zones of the study area. Conversely, areas with high suitability extend over a larger area of 9.81 km 2 , representing 50.19% of the undeveloped lands. Areas with moderate suitability reach over 6.85 km 2 or 29.38% of the undeveloped lands, whereas low suitability areas cover 0.63 km 2 , accounting for 17.87% of undeveloped lands. The areas with the lowest degree of suitability represent only 1.99% of the undeveloped areas, with a total area of 2.87 km 2 . Consequently, the study implies that future development projects could be implemented, generally, in areas with high and moderate suitability degrees for SUD since these zones have larger spatial extent. Figure 6 shows the percentages of suitability degrees relative to the non-built areas in the case study. According to the analysis, 64.11% of agriculture areas have the highest and high suitability degrees, with a total value of 5.09 km 2 . While 53.37% of bare soils and grass areas, equivalent to 3.8 km 2 , are considered as high and highly suitable lands. Moreover, 19.77% of the wooded areas, equivalent to 1.03 km 2 , are considered as high and highly suitable lands for SUD. A total of 35.64% of agriculture areas have a moderate suitability degree, with a total value of 2.83 km 2 . While 19.8% of bare soils and grass areas, equivalent to 1.41 km 2 , are considered moderately suitable lands. Moreover, 50.1% of the wooded areas, equivalent to 2.61 km 2 , are considered moderately suitable lands for SUD.
A total of 0.25% of agriculture areas have low and lowest suitability degrees, with a total value of 0.02 km 2 . While 26.83% of bare soils and grass areas, equivalent to 1.41 km 2 , are considered as low and least suitable lands. Moreover, 30.13% of the wooded areas, equivalent to 2.61 km 2 , are considered low and least suitable lands for SUD. It is noteworthy that this result is not in contradiction with the suitability values depicted in Table 3, but then other factors, mainly the slope, the distance from urban agglomeration, and main roads, would further affect this allocation rather than the land cover distribution. Table 4 shows the crosstabulation of low to high lands with land cover types. According to the analysis, 64.11% of agriculture areas have the highest and high suitability degrees, with a total value of 5.09 km 2 . While 53.37% of bare soils and grass areas, equivalent to 3.8 km 2 , are considered as high and highly suitable lands. Moreover, 19.77% of the wooded areas, equivalent to 1.03 km 2 , are considered as high and highly suitable lands for SUD. A total of 35.64% of agriculture areas have a moderate suitability degree, with a total value of 2.83 km 2 . While 19.8% of bare soils and grass areas, equivalent to 1.41 km 2 , are considered moderately suitable lands. Moreover, 50.1% of the wooded areas, equivalent to 2.61 km 2 , are considered moderately suitable lands for SUD.
A total of 0.25% of agriculture areas have low and lowest suitability degrees, with a total value of 0.02 km 2 . While 26.83% of bare soils and grass areas, equivalent to 1.41 km 2 , are considered as low and least suitable lands. Moreover, 30.13% of the wooded areas, equivalent to 2.61 km 2 , are considered low and least suitable lands for SUD. It is noteworthy that this result is not in contradiction with the suitability values depicted in Table 3, but then other factors, mainly the slope, the distance from urban agglomeration, and main roads, would further affect this allocation rather than the land cover distribution.   Table 4 shows the crosstabulation of low to high lands with land cover types. Upon analyzing the land suitability map in Figure 5b and identifying the four indexed zones (A, B, C, and D) as the highly suitable areas for SUD, it is observed that areas with high suitability degrees are characterized by elevations ranging between 400 and 600 m, and by topographic slopes ranging between 0 and 10 degrees.
Zones A and D are located in proximity to the major roadway that links the study area to the principal coastal cities (Saida and Tyre) through the coastal highway. Additionally, the same road connects the study area to prominent mountain villages (Marjeyoun, Hasbaya, and Rachaya) in the opposite direction since the study area serves as an access crossing point to neighboring regions (Figure 7). access crossing point to neighboring regions (Figure 7). Zone B is located in proximity to the main urban agglomeration and to local roads, resulting in higher levels of accessibility, more developed infrastructure, and a concentration of economic activities in addition to social and cultural amenities.
Zone C is located in proximity to industrial and commercial areas and to local roads that provide high accessibility to neighboring villages. Figure 7. The spatial location of the study area relative to the main nearby coastal cities.

Conclusions
The expansion and densification of urban areas can contribute to the degradation of green areas and natural resources. These green spaces perform a critical role in preserving biodiversity and offering various natural benefits. Moreover, urbanization can increase the consumption of natural resources and lands, further impacting the natural environment. However, sustainable urban development practices aim to balance the needs of urbanization with the preservation and protection of green land and natural resources to achieve more sustainability.
It is vital to ensure that the environment is not harmed by future urban development and to keep a planet with sustainable natural resources for future generations.
This study contributes to the emerging literature by integrating sustainability into the land suitability analysis. This could represent an important guidance for decisionmakers, particularly when authorizing strategies and developing zoning and urban masterplans. In detail, this analysis aims to evaluate the land suitability to identify the most Zone B is located in proximity to the main urban agglomeration and to local roads, resulting in higher levels of accessibility, more developed infrastructure, and a concentration of economic activities in addition to social and cultural amenities.
Zone C is located in proximity to industrial and commercial areas and to local roads that provide high accessibility to neighboring villages.

Conclusions
The expansion and densification of urban areas can contribute to the degradation of green areas and natural resources. These green spaces perform a critical role in preserving biodiversity and offering various natural benefits. Moreover, urbanization can increase the consumption of natural resources and lands, further impacting the natural environment. However, sustainable urban development practices aim to balance the needs of urbanization with the preservation and protection of green land and natural resources to achieve more sustainability.
It is vital to ensure that the environment is not harmed by future urban development and to keep a planet with sustainable natural resources for future generations.
This study contributes to the emerging literature by integrating sustainability into the land suitability analysis. This could represent an important guidance for decision-makers, particularly when authorizing strategies and developing zoning and urban masterplans. In detail, this analysis aims to evaluate the land suitability to identify the most suitable areas for sustainable urban development in the Nabatiyeh area, where there was a rise in an unplanned and accelerated phase of urbanization and sprawl during the past 20 years. The adopted methodology for this spatial assessment is based on using a multi-criteria decision analysis technique, the analytical hierarchy process, combined with a geographic information system that includes several suitability assessment criteria.
When determining suitable areas for sustainable urban expansion, it is crucial to consider socio-economic, environmental, and geophysical factors to achieve sustainability.
To this end, the scientific literature and expert judgments have identified several parameters that must be taken into account, including slope, elevation, land cover types, proximity to urban agglomerations, proximity to commercial and industrial zones, and accessibility to main roads. To obtain the data necessary for this assessment, the digital elevation model, land cover map, and road network were collected from different sources: the USGS Earth Explorer website, the National Center for Remote Sensing in Lebanon, and the OpenStreetMap database, respectively. According to expert judgments, and the literature indications in some cases, the analytical hierarchy process indicates that the elevations and slopes are the most important criteria to be considered, with corresponding weights of 34% and 36%, respectively. Additionally, a moderate importance weight value of 0.12 was assigned to the distances from urban agglomeration. Low values of criteria importance were assigned to the distance from industrial and commercial areas, the distance from main major and local roads, and the land cover types with weight values of 0.09, 0.06, and 0.03, respectively. The suitability map for future sustainable urban development has been created and categorized into five classes based on their degree of suitability. The analysis revealed that regions with high suitability accounted for 51% of the total undeveloped area; however, areas with moderate suitability are extended over 29%, and areas with low suitability degrees constitute 20% of the total undeveloped area.
The high and highest suitability areas identified in the land suitability map offer a promising opportunity for sustainable urban development, especially since they take into consideration the region's topography, urban agglomeration, and existing infrastructure. This identification assists urban planners and policymakers in making informed decisions regarding the development of urban areas, ensuring that they are sustainable, efficient, and meet the needs of the population. Regarding the limitations, the neighboring villages were excluded from our study due to insufficient data. However, these villages share similar physical and socio-economic characteristics with the study area; hence, the findings would potentially be the same if these villages were included in the study area.
These findings illustrate the efficiency of combining multi-criteria decision analysis with geospatial techniques for urban planning and environmental management in promoting sustainable development and protecting the environment. This study offers vital information and insights for decision-makers on sustainable urban development using a multi-criteria decision analysis combined with geographic information systems and remote sensing, and it represents a necessary step toward sustainable growth and mitigation measures against natural degradation. With regard to future research, it is recommended in this context to study the possibility of integrating administrative reforms to promote the studies and implementation of sustainable urban development. These studies could be based on some practices of administrative reforms, including improved governance, consolidating, and enhancing the efficiency of governmental agencies and local authorities that are responsible for urban development, building more scientific capacities for the aforementioned entities, and increasing accountability of officials. Moreover, the existing land use and zoning master plan could be updated according to the land suitability map to provide more guidance toward sustainable development. Funding: This research received no external funding.

Data Availability Statement:
The data presented in this study are available on request from the corresponding authors.

Conflicts of Interest:
The authors declare no conflict of interest.

Appendix A
ing-review and editing, W.A.-S., G.F. and A.E.; visualization, A.A.-G. and Y.A.; supervision, W.A.-S. and G.F.; All authors have read and agreed to the published version of the manuscript.
Funding: This research received no external funding.

Data Availability Statement:
The data presented in this study are available on request from the corresponding authors.
Acknowledgments: Sincere gratitude to the CNRS-Lebanon for their technical and administrative support in providing us with the necessary data, which has been instrumental in our research.

Conflicts of Interest:
The authors declare no conflict of interest.
Appendix A Figure A1. The steps followed in the AHP. Figure A1. The steps followed in the AHP.