Land-Use Suitability Assessment Using Delphi and Analytical Hierarchy Process (D-AHP) Hybrid Model for Coastal City Management: Kuala Terengganu, Peninsular Malaysia

: Since at least half of the world’s population resides and works within coastal land, the coastal zone processes and resource management is of great economic and social importance. One of the fundamental issues for coastal city planners, researchers, managers, and engineers is the coastal city land-use suitability. Land-use suitability is the ability of a given type of land to support a deﬁned use. Rapid urbanization and consequent haphazard growth of cities result in deterioration of infrastructure facilities, loss of agricultural land, water bodies, open spaces, and many micro-climatic changes. Hence, accurate data on coastal city hazards are essential and valuable tools for coastal planning and management, sustainable coastal development, coastal environment conservation, selection of a site for coastal city structures, and coastal resources. In this investigation, the Delphi and Analytical Hierarchy Process (D-AHP) Hybrid model and Geographic Information System (GIS) technique for Coastal Land-Use Assessment (CLUA) are mapped to detect the most suitable and unsuitable areas in the Kuala Terengganu coastal zone. Furthermore, this research offered information not only on the present urban land-use trend and established amenity status in Kuala Terengganu, but also on the suitability of land for the potential establishment of urban facilities for improved urban planning and appropriate decision-making. Using the D-AHP Hybrid model and GIS tool for coastal city management is broadly practical for government, policymakers, and planners to appropriately strategize and plan for the future of coastal cities in Malaysia and other analog coastal cities around the world.


Introduction
Land suitability assessment is a land evaluation approach that identifies the key qualitative and quantitative evaluation limiting criteria [1]. Climate, hydrology, topography, vegetation, and soil characteristics are examined in qualitative land suitability assessments, whereas the results are more thorough, and yield is calculated in quantitative assessments [2]. The systematic examination of land and water potential, land-use alternatives, and economic and social circumstances to choose and implement the optimal land-use options is known as land suitability assessment [3].
Coastal city lands have relatively stable or consistent cyclic behavior as part of the Earth's surface [4]. Coastal city land-use suitability assessment for agriculture, settlement, housing, and industry depends on its elements. However, these elements have been depleted, especially in the last few centuries [5]. Coastal city lands refer to the areas foretasted respectively −1 to 32% and +0.6 to 3.4 • C in 60 years, whiles, the rise of sea level due to global warming is forecasted roughly 1. 30-9.40 mm/year in 100 years in Peninsular Malaysia [26].
Ref. [27] observed that in East Coast of Malaysia, apart from erosion, accretion also happens within the neighborhood of a 25-km stretch of coastline, especially throughout the northeast monsoon season. However, in the coastlines of Kelantan Terengganu, little effort had been made to explain this erosional discrimination feature of the northeast monsoon waves. Kuala Terengganu's climate is dominated by the northeast monsoon that is prevalent from November to March, and the yearly temperature ranges from 25.6 • C to 27.8 • C [28]. The observed average sea-level rate near the Malaysian coast was estimated between 2.7-7.0 mm/year based on satellite altimetry data from 1993 to 2010. In Peninsular Malaysia, the anticipated sea-level increase for the year 2100 was about 2. 50-5.20 mm/year, with the low-lying northeast and west coast portions of the peninsular (Kelantan and Kedah) having the highest value [29]. According to the [30] study, around 3.3% of the 1963 km of coastline is classified as severely vulnerable, with the northern reaches of Kedah's shoreline and the southern lengths of Terengganu's beachfront being among these locations. The study revealed a significant increase in the trend of sea-level rise in the past 5 years, compared to 20 years ago. The National Coastal Erosion Study (NCES) found that of a total of 4809 km of the coastal area of Peninsular Malaysia, roughly 1360 km, covering 28% of the coastline currently threatened with the erosion of varying degrees of intensity, threatening the economic values of human activities [31].
The Malaysia coastal land erosion problem arises mainly due to the dynamic process, human activities. Evidence from past studies indicated spatial variations of coastal lands and erosion, and accretion in Terengganu coastal area disclosed the geomorphological changes. Kuala Terengganu is one of the areas that was endangered by serious erosion and land loss, and it became necessary to monitor its coastal land [32]. Conversely, another factor that also causes coastal city erosion and land-use change are human activities such as settlements, industrial development, agriculture, deforestation, and coastal protection in the coastal area within watersheds, river catchments, and offshore [33]. The combination of human activities and natural forces, particularly development in coastal lands usually worsen coastal city erosion [34]. The effects of coastal city erosion are generally because of inadequate foresight in the planning and location of infrastructure or facilities [35]. The continuous location of facilities in coastal areas that are vulnerable to erosion could aggravate the problem of coastal land erosion except for developers appropriately consider the causes and effects of erosion [36]. In the areas that have developed the central threat is on the houses or residential structures as well as the local access roads linking those houses several housing and residential land might be lost if the erosion processes are left uncontrolled. The likelihood of an occurrence of intensity equal to or greater than a category value occurring during a period of exposure is described as a hazard. Subsidence, sea-level rise, erosion are examples of environmental and climatic hazards (H) in the research region [37].
Finally, risk (R) is defined as the likelihood of a loss impacting an element occurring because of a high-intensity event. Risk is calculated as a function of hazard and vulnerability and was assessed using the equation R = (H, V) [38]. The likelihood of a disaster and the financial worth of potential losses or consequences are frequently the only factors considered in risk assessments. These findings are used to determine if preventing certain calamities is economically possible [39]. Risk assessment is an important step in the risk management of land use, and many researchers have provided their valuable work in this field during recent decades [40].
The assessment of coastal land suitability is analogous to the selection of a land-use position and attempts to map a suitability index for the whole research area [41]. It is necessary to choose the most suitable coastal land suitability evaluation algorithm for current and future land-use planning [42]. One of the most significant prerequisites for coastal land-use decisions is the creation of suitability maps [43]. It is important to consider many conditions and requirements related to coastal land and its uses.
Social and economic influences are far more manipulated than bio-physical and environmental factors through human interventions [44]. With many of the cities clustered near the sea, there are long lands in Peninsular Malaysia, and over 90% of the coast is among the easily erodible [45].
Some essential economic sectors in Malaysia are confronted with adverse effects of climate change and such sectors include industries, energy, transport, agriculture, forestry, public health and waste sector as well as water and coastal resources. Therefore, rapid responses via adaptation and mitigation approaches are required to curtail the risks, possible losses as well as maximize the beneficial opportunities. The Malaysia government has taken diverse programs and actions to reduce the effect of climate change in all important sectors [46].
The East Coast Economic Region (ECER) in Malaysia has a widespread and diverse coastal city zone, and as with several coastal areas of the globe, would also be influenced by future human-oriented climate change especially via increased coastal land erosion [47]. In essence, low-lying areas that have large socio-economic activities and populations are in danger of being inundated. For instance, the low-lying Kuala Terengganu area has more susceptibility to coastal erosion throughout the monsoon season especially [48].
This study selects Kuala Terengganu as the study area due to its location which represents an active tourism gateway to the ECER. There are several distinctive tourism attractions especially in the areas of mainland coastal and sustainable island tourism, ecotourism, urban tourism, and its typical cultural and heritage tourism [49]. The ECER Malaysia was established as an economic corridor with the objective of creating the socio-economic transformation of the east coast of Peninsular Malaysia. Thus, ECER covers the states of Terengganu, Kelantan, Pahang as well as the and the districts of Johor and Mersing [50].
Therefore, assessing the weight of each factor influencing the suitability of coastal land is a crucial step in evaluating the suitability of coastal land use for hazard evaluation. Because for the interest groups concerned, each element has its degree of value or weight, it can be difficult to weigh them [51].
Land suitability assessment using a Geographic Information System (GIS) has been widely used in the real world, for example, mapping [52], urban growth [53], managing transportation issues [54], Environmental Impact Analysis [55], risk and disaster management [56], flood risk areas [57], Natural Resources Management [58], Geology [59], Planning and Community Development [60], Irrigation Water Management [61], regional planning and land-use suitability for agricultural activities [62]. GIS tools, according to [63], were proven to be a quick tool for CLUA, especially in partial data settings. The GIS-based Multi-Criteria Decision Analysis (MCDA) is the most widely used approach in various literature to assess land-use suitability for future coastal city sustainable development MCDA [64]. This method uses several approaches for creating decision problems as well as weighing and ranking various choices [65]. To establish the weighting and score of various criteria, researchers employed a variety of methods [66].
Multi-Criteria Decision Analysis (MCDA) is used to help in complex decision-making, particularly in uncertain situations where there exist a limited practicable set of options, and to select the best alternative according to the scores of many characteristics and its multi-attribute decision-making problem [67]. MCDA offers greater-supported methods to compare project alternatives based on decision matrices. It also gives structured approaches to incorporate the opinions of project stakeholders in the ranking of alternatives [68]. One vital merit of MCDA is its ability to bring out the similarities and possible conflict areas among stakeholders in group decision-making that enhance a thorough understanding of the values of others [69]. The Analytical Hierarchy Process (AHP) is a multi-criteria analysis approach extensively used to aid individual and group decision-making across the world [70]. The AHP remains one of the broadest techniques in the MCDA setting [71]. AHP is frequently used in conjunction with other techniques, with writers estimating the weight of the criterion using AHP. Fuzzy AHP is one of the most widely used methods of MCDA, with many articles devoted to a comparison of classical AHP with fuzzy AHP [72]. The current method of MCDA used specifically for environmental sciences includes the AHP, best-worst method (BWM, level-based weight assessment (LBWA), and full consistency method (FUCOM) [73]. The best-worst method (BWM) is a MCDA approach for weighing a collection of options against a set of decision criteria [74]. The BWM is based on a pairwise comparison of the choice index in a methodical manner [75]. The new levelbased weight assessment (LBWA) approach allows experts from many areas to collaborate to define the relationships between indicator and make reasonable decisions [76]. The approach may be used in real-world situations in specialized decision-making support systems, as well as in virtual settings for alternative conflict settlement [77]. The complete consistency method (FUCOM) is a novel methodology for weighing literature criteria [78]. This approach is a semi-objective assessment method that lowers the number of indexes comparisons inside each other and optimizes the criterion weights with a few comparisons using an optimization algorithm [79].
There are many advantages of AHP techniques, namely it is relatively easy to use, can integrate qualitative and subjective factors, tolerates speedy replanning, obtain the greatest judgmental decision estimate from numerous credit analysts, and offers a method to assess the consistency of these judgments [80]. The most fascinating aspect remains its ability to include judgments and personal values reasonably therefore making a complex decision problem to be structured hierarchically [81]. This viewpoint was corroborated by [82] who submitted that AHP is the frequently and extensively applied as well as fastest emergent decision analytic approaches. AHP has been used in numerous complex decisions especially in the field of environment, namely the challenges of management of natural resources [83], Integrated Coastal Zone Management (ICZM) [84], environmental issues [85], social science, and technology [86], economic studies [87], risk management process [88], land-use planning [89], urban planning [90].
The Delphi technique might be used to improve the quality of expert views used in the AHP approach [91]. The Delphi technique has been applied in numerous studies in the fields of environmental studies, project management, and industrial engineering [92]. Therefore, this technique is suitable as a research instrument especially when the knowledge of a phenomenon or problem is incomplete [93].
Since the 1990s, coastal city planners have significantly increased their attention to the incorporation of the MCDA approach with GIS to solve the problems of spatial planning [94]. The ability of GIS to handle spatial facets of risk assessment has improved its' application in the index-based assessment for prioritization and selection of possible appropriate and inappropriate areas. This is since most of the conditions for CLUA are spatial data [95]. GIS tools were a speedy tool for coastal city land-use suitability, especially incomplete data situations [96]. In support of currently available literature [97], we believe AHP is the most suitable technique that widely used MCDA methods in decision support system processes which includes multiple actors, scenarios, and GIS layers classes [98].
The overall goal of this research is to create risk and land used suitability maps for hazard management in Kuala Terengganu, a coastal city on Peninsular Malaysia's east coast. We have developed a CLUA model by integrating GIS tools and the D-AHP Hybrid model for coastal city land use. The combination of the GIS technique and the D-AHP model is a powerful approach that used risk assessment and land-use suitability in coastal city areas [99].

Study Area
Terengganu is a constitutive state of federal Malaysia (Peninsular Malaysia) and is bordered in the northwest by Kelantan, the southwest by Pahang, and the east by the South China Sea. Kuala Terengganu which is located at the mouth of the broad Terengganu River is the largest city in Terengganu state. The eastern part of Kuala Terengganu where faces the South China Sea is characterized by sandy coastal beaches that cover the entire stretch of both parts ( Figure 1). The study area is selected along the coastline of Kuala Terengganu which extends from Merang to Rusila for about 68.84 km. In this area, the coastal area was characterized by the northeast monsoon, with semidiurnal tides and an average annual temperature of 25 • C to 27 • C. The shoreline in this area is categorized as a critical region by the local authorities. It is mainly due to the dynamic physical process in the area, sea-level change, and human activities.

Study Area
Terengganu is a constitutive state of federal Malaysia (Peninsular Malaysia) and is bordered in the northwest by Kelantan, the southwest by Pahang, and the east by the South China Sea. Kuala Terengganu which is located at the mouth of the broad Terengganu River is the largest city in Terengganu state. The eastern part of Kuala Terengganu where faces the South China Sea is characterized by sandy coastal beaches that cover the entire stretch of both parts ( Figure 1). The study area is selected along the coastline of Kuala Terengganu which extends from Merang to Rusila for about 68.84 km. In this area, the coastal area was characterized by the northeast monsoon, with semidiurnal tides and an average annual temperature of 25 °C to 27 °C. The shoreline in this area is categorized as a critical region by the local authorities. It is mainly due to the dynamic physical process in the area, sea-level change, and human activities.

Data Used
The study collects two kinds of data, namely primary and secondary data, for the CLUA framework (Table 1). Thus, the D-AHP model's primary data were obtained via the Expert Choice Matrix (ECM) and were given to 12 specialists from various agencies or departments in Malaysia (DID, SMART, NAHRIM). Thus, the ECM was run for five-month (for three runs). The study conducts the first run between 14

CLUA Model Framework
The D-AHP Hybrid model's semi-quantitative architecture has improved system consistency as well as weight computation for layers and classes [100]. The D-AHP method is organized hierarchically at diverse levels with each level containing a finite amount of decision components. The general goal of CLUA is represented by the upper level of the hierarchy whereas one or more intermediate levels symbolize the decision GIS layers and classes. Hence, the lower level comprises all potential options (Zones) as alternatives. Applies the D-AHP model to assess the vulnerability, hazard layers, and classes on coastal city region's erosion that necessitated semi-quantitative analysis of the coastal region's information. The semi-quantitative design uses two GIS layers for the measurement of erosion risk, namely vulnerability and hazard layers. Therefore, D-AHP's first step is to create a hierarchy by shaping the problem into its basic elements. This model ascribes significance weights to the variables with the use of pairwise comparison at every hierarchy level. The class of the risk appraisal process is that the pairwise compression of the layers at a hierarchy level with the objective at the next level. CLUA components are based on their common features and the development of a hierarchical model that has diverse levels [101].
Twelve phases are involved in the CLUA framework. Phase 1 of the D-AHP Hybrid structure is associated with the configuration of the decision problem into hierarchical coastal land use. Thus, vulnerability and hazard maps in semi-qualitative hierarchy design could involve vulnerability and hazard Index ra-tios. The AHP semi-quantitative design, vulnerability, and Hazard maps may encompass vulnerability and hazard GIS layers and ratios ( Figure 2). ISPRS Int. J. Geo-Inf. 2021, 10, x FOR PEER REVIEW 8 of 37 Phase 1 of the D-AHP Hybrid structure is associated with the configuration of the decision problem into hierarchical coastal land use. Thus, vulnerability and hazard maps in semi-qualitative hierarchy design could involve vulnerability and hazard Index ratios. The AHP semi-quantitative design, vulnerability, and Hazard maps may encompass vulnerability and hazard GIS layers and ratios ( Figure 2).  Phase 2 comprises an empirical study under the Delphi method and expert choice matrix to enhance gathered essential layers and classes from a review of the literature and to examine other layers emanating from the opinion of the expert panel.
Phase 3 is associated with the collection of input data by pairwise comparisons of decision components as well as obtaining the judgment scales. Having arranged the problem in a hierarchical technique, the subsequent step is to appraise the relative significance of the emanating and classes to the general goal. The study applies two designs expert choice matrix with the first design being a 5-point expert choice matrix of which the user enables the selection of the best and essential layers from the expert. The second design signifies a 9-point expert choice matrix which is used for the semi-quantitative data layers. By using pairwise relative comparison, the preference between layers and GIS classes with a 9-point rating scale would be asked for deriving weights by the AHP model. The relative significance of every variable is evaluated through the application of Saaty's nine scales. For the semi-quantification process in AHP, the pairwise comparisons of layers and classes use a scale of unconditional judgments generally known as Saaty's nine scales, which specifies the degree at which one item prevails over the other regarding the certain feature. This scale changes semi-quantification assessments into statistical values from 1 to 9 which is used for semi-quantitative GIS classes to complete the pairwise comparison matrix [102].
Phase 4 is concerned with the evaluation of the weight consistency of comparisons. This step conducted the calculations to determine the maximum eigenvalue, Consistency Index (CI), Consistency Ratio (CR), and normalized values for every layer and class [103]. The AHP determines the general consistency of judgments through a CR Equation (1). As noted by some studies [104,105], the judgmental matrix consistency can be accomplished through the examination of total CR, the ratio of CI and RI denoted as follows: where RI signifies the random consistency index of a reciprocal matrix that was randomly created from the 9-point scale, with reciprocals forced [106] (the average CI of 500 randomly filled matrices). A CR below 10% suggests that the matrix could be considered to be possessing an acceptable consistency [107]. Thus, [108] computed the random indices as indicated in Table 2. The scores that attained a certain level of consistency were accepted as demonstrated by a CI calculated by [109] Equation (2). He proposed a CI that is denoted as follows: The CR allows determining the number of errors that were generated when offering expert judgments. The calculated CI can be compared with the created paired comparison matrix RI to examine CR ( Table 2). The CR indicated whether the judgment weights of the decision-maker were accepted. A value of CR below 0.1 reveals that the errors are relatively small and considered acceptable. This is because human judgments do not necessarily need to be always consistent since some inconsistencies may be introduced due to the nature of the scale used [110].
Phase 5 concerns local weights and consistency of comparisons using expert choice software for the calculation of the D-AHP weight. This step calculates the local weights of the components by using expert choice software. Phase 6 deals with the aggregation of weights through several levels to obtain the ultimate weights of GIS layer classes. After obtaining the local weights of components of diverse levels as highlighted in the preceding step, they are combined or aggregated to obtain the ultimate weights of the decision options. Thus, the layer class weights known as local priorities are seen as decision components in the decision process second step. The decision-maker must offer these preferences or priorities by pairwise comparisons concerning the weights and scores. The weights (w j ) and scores S ij values are stimulated from these comparisons and denoted in a decision table.
Phase 7 of the AHP combines or aggregates all preferences from the decision table with a weighted summation of the kind [111] and creates the local priorities through all layers to determine the global preference or priority [112] Equation (3). The ultimate GIS layer classes weights in semi-quantitative data would be calculated using the formula given as follows: where Z I represents the global priority of the GIS layer classes, S ij signifies local priority while w j denotes the weight of the layer classes j. Thus, global priorities (Z I ) attained in this research are ultimately used to normalize through the division of the score of every layer class solely by the layer classes score and choose the optimum classes in each layer. Phase 8 To produce the vulnerability and hazard map, all shapefile layers are transformed to raster layers and each vulnerability and hazard layer is computed and multiplied with D-AHP weights in Arc Map.
For Phase 9, 10 and 11, a seven-class method was used to categorize the assessment findings of hazard, vulnerability, and risk maps to derive the land-use suitability region from the risk map. When executing the categorization technique, assigning rank to the alternative of the values (1) Extreme low (null); (2) Very low (not subjected areas); (3) low (potential areas); (4) medium (fragile areas); (5) high (critical areas); (6) Very high (degraded areas); (7) Extreme high (No data) are proportionally assigned to give more sensitive vulnerability, hazard, and risk levels in the coastal city region [113,114].
In this research vulnerability, hazard and risk map the classified into seven classes because, according to several Malaysian reports (DID, 1985, and NAHRIM 2010), the Kuala Terengganu coastline area has extreme erosion, the low and high extreme erosion classes are significant in determining the best suited and unsuitable areas. on the other hand, Kuala Terengganu coastline, offers a broad range of tourism attractions, including mainland coastal and sustainable island tourism, ecotourism, urban tourism, traditional culture and historical tourism, and traditional culture and historical tourism. As a result, we need an accurate and significant classification of vulnerability, hazard, and risk in cities and coastal areas, which may help planners and politicians make better decisions.
Phase 9 vulnerability assessment that refers to a susceptible element's vulnerability to a hazard, such as human populations, physical structures, economic assets, or sensitive environmental components. Most of them have quantitative measurements in place, such as the number of lives lost or the monetary value of serious hazard [115].
The D-AHP weights have been calculated and multiplied with each vulnerability layer, using GIS tools Equation (4).
where (V) represents vulnerability, (L v ) signifies vulnerability layer and (W vlc ) denotes the weight of the vulnerable layer class. Phase 10 hazard is a potentially dangerous event that might threaten the lives of people or damage to property or the environment. Severity, length, amount, intensity, and frequency are all factors to consider. A bigger threat is represented by a higher numeric number [116].
The capacity of GIS tools to incorporate the D-AHP weight, which is connected to each hazard layer, has been demonstrated in this section of the research Equation (5).
where (H) signifies hazard, (L h ) represents hazard layer and (W hlc ) denotes the weight of the hazard layer class. Phase 11 The possibility of a calamity is referred to as risk. On a small range, it can be measured as a weighted rank [117].
The development of the CLUA followed a combination of the above model's Equation (6) in the GIS [118]. Thus, the hazard index was multiplied by the vulnerability index in GIS to generate the risk map as indicated s below: where (R) signifies the risk layer of CLUA, (H) denotes the hazard layer, and (V) symbolizes the vulnerability layer. The variability of catastrophes is referred to as the hazard potential factor. In general, as catastrophe intensity and frequency rise, so does the severity of the damage and loss. All persons and things potentially threatened by the hazard factors, such as local inhabitants, animals, residences, crops, and so on, are referred to as impacted objects. An area with a higher property density that is exposed to risk factors will have a higher loss potential.
Phase 12 This project calculates and generates a land-use suitability map from a risk map with seven classes in a GIS that depicts the sensitivity and suitability in the coastal city area. The three varieties of low-risk classes on the risk map are suitable area, whereas the three types of high-risk classes are unsuitable areas. There are many null lands and deforestation areas in the coastal city region, particularly in Kuala Terengganu. Therefore, the seven rankings developed for this study are beneficial for determining the optimum land-use suitable location on a small regional scale with considerable coastal erosion. The final risk maps from GIS for the Kuala Terengganu coastal zone used to create two-class maps in suitable and unsuitable for six zones.

Analysis of Vulnerability Maps on the Local Scale
The vulnerability assessment for the semi-quantitative model design refers to a region with the specific relative weakness of exposure to natural disasters and their ability to handle them when they occur. Based on the results of the model from experts' opinion, the study constructed, classified seven coastal vulnerability layers into two domains, namely environmental layers and human activity layers. The environmental layers domain refers to the soil, slope, geology, digital elevation map, river, and land use. In the human activity layers, the population density, built-up, and road network were taken into consideration. Conversely, another factor that also causes shoreline change and coastal city erosion is human activity in the coastal area within watersheds, river catchments, and offshore [116]. The combination of human activities and natural forces, particularly development in the coastal areas usually worsen coastal erosion [117].
The environmental layers classes are weighted due to serious coastal city erosion problems. The computed weights for each environmental layer's classes will be used in the GIS analysis for an estimate of hazard, vulnerability, and risk area on the local scale Table 3.  Figure 3 shows the final results of the vulnerability analysis for the 3 vulnerability erosion zones of the Kuala Terengganu coastal areas that were selected and choose this area base on literature and some government agency reports. Expert knowledge was used to rank vulnerability parameters for each zone from extremely high (7)

Analysis of Hazard Maps on a Local Scale
The hazard assessment was done with the D-AHP analysis of the Kuala Terengganu and the GIS map of different classes is presented in Figure 4. The GIS layers and D-AHP weights are illustrated in Table 4. In Batu Rakit (Zone 2) with 11,500 ha, 6.2% (710 ha) of the area is considered to be an extremely highly vulnerable, 14.16% (1678 ha) is seen as very highly vulnerable, 38.7% (4452 ha) is highly vulnerable, 13.7% (1580 ha) represents moderate vulnerable, 10.9% (1254 ha) signifies low vulnerable, 6.4% (739 ha) is considered to be very low vulnerable and 9.5% (1088 ha) represents extremely low vulnerable ( Figure 3).

Analysis of Hazard Maps on a Local Scale
The hazard assessment was done with the D-AHP analysis of the Kuala Terengganu and the GIS map of different classes is presented in Figure 4. The GIS layers and D-AHP weights are illustrated in Table 4.        Figure 4 presents the hazard potential of 3 hazard erosion zones of Kuala Terengganu. The three zones that have the highest rate of hazard, namely Zone 2 (Batu Rakit), Zone 3 (Kuala Nerus), Zone 5 (Pengadang Buluh) and are also presented. In Batu Rakit (Zone 2) with a total area of 11,619 ha, 6.06% (704.56 ha) is an extremely high hazard area, 6.11% (710 ha) is a very high hazard area, 15.99% (1858 ha) is a high hazard area, 28.77% (3342 ha) is moderate hazard area, 31.35% (3642 ha) is low hazard area, 5.61% (652 ha) is very low hazard area and 6.11% (710 ha) is an extremely low hazard area. Kuala Nerus (Zone 3) with a total area of 10,590 ha has 1.6% (173 ha) as an extremely high hazard area, 1.7% (176 ha) as a very high hazard area, 9.9% (1047 ha) as a high hazard area, 34.5% (3655 ha) as moderate hazard area, 44.7% (4732 ha) as low hazard area, 7.1% (748 ha) as very low hazard area and 0.5% (58 ha) as an extremely low hazard. Finally, Pengadang Buluh (Zone 5) which covers a total area of 9023 ha has an extremely high hazard area of 20 ha (0.2%), a very high hazard area of 122 ha (1.3%), a high hazard area of 1591 ha (17.6%), moderate hazard area of 2784 ha (30.9%), low hazard area of 3515 ha (39%), very low hazard area of 773 ha (8.6%) and extremely low hazard area of 219 ha (2.4%) (Figure 4).
The result shows that there is an extreme flood in Zones 3, 2, and 5 as well as a river. Zones 2 and 3 are in low-lying and low elevation in the DEM layer in GIS. There are shoreline changes and coastal erosion in transects 10, 11, 12, and 13 with the rates of 6.44, 5.90, 6.92, and 3.95 m/year, respectively, which are located in Zones 2 and 3. Pengadang Buluh (Zone 5) has the lowest hazard with a rate of 19.2% compared to Zones 3 and 2. This could be due to the construction of some coastal protections, namely brake water and a jetty in this area. 3) with a risk rate of 57.9%, and Batu Rakit (Zone 2) with a risk rate of 35.5%. A change in the vulnerability of the area would also change the rate of risk. This shows that Zone 5 has the highest potential to be significantly affected by the negative impacts of climate change hazards.

Assessment of Land-Use Suitability Area
The result of land-use suitability and unsuitability assessment for Kuala Terengganu shows that among three zones (2, 3, and 5), Zone 2 has the most suitable area, followed by

Assessment of Land-Use Suitability Area
The result of land-use suitability and unsuitability assessment for Kuala Terengganu shows that among three zones (2, 3, and 5), Zone 2 has the most suitable area, followed by Zones 3 and 5. Out of 8524 ha land area in Batu Rakit (Zone 2), about 52.1% (4444 ha) is a suitable area and around 47.9% (4080 ha) is an unsuitable area ( Figure 6). Out of 9847 ha total land area in Kuala Nerus (Zone 3), about 38.8% (3818 ha) is considered to be a suitable area 61.2% (6029 ha) is an unsuitable area. Out of 6708 ha of the total land of Pengadang Buluh (Zone 5), only 29.8% (2001 ha) is a suitable area while about 70.2% (4706 ha) is an unsuitable area ( Figure 6). It shows that Zone 5 is less suitable for construction, agriculture, and tourism activities compared to other zones. The suitability maps provide the most satisfactory means of conveying the result of the assessment to the user in a summary form. Overall, the northeast and southeast parts of the study area had been categorized as suitable areas, because these regions are less vulnerable to hazard factors such as a river, flood, coastal hazards, and human activity compared to the middle part of the study area. The matching of current land-use type and land suitability classification can support and contribute information for site selection and planning for the government. For example, some of the government agencies, universities, and international airports, as well as build-up areas near the coast, are in highly vulnerable and risk zones. Therefore, before undertaking any plan and development in the coastal area, the suitability or unsuitability of the land for a specific purpose should be taken into consideration. Suitability analysis allows regional planners to make timely and more informed decisions. Furthermore, it determines how the land would be used in the future. Hazards and risks affect the future development of land. Hence, it is important to include hazard and risk information in the land-use planning process, particularly in the areas facing natural hazards.

Discussion
The coastal city erosion vulnerability control and management plan is a continuing plan requiring federal government actions that the nature of this action will change with time but will change gradually through able to be foreseen future. The federal government actions are reactionary, and they must produce management, decisions, structural control works to protect facilities developed over past years. Additionally, they produce guidance to ensure that disastrous erosion risk consequences do not threaten facilities developed in future years.
Classification is an informed judgment on the environmental, socio-economic, and physical impacts of coastal city risk assessment. That judgment depends on the integration of the many technical analyses and the varied technical information comprising. The rate of erosion is the base for determining when and how seriously, a value will be impacted. Values considered are those generally affected by erosion in the Kuala Terengganu coastal area. To achieve these ends, the eroding areas are classified and grouped into three classes (Table 5): Housing ( Critical erosion risk areas, the anticipated socio-economic impacts indicate that the social and economic costs of erosion may soon be unacceptable. Possibility studies to more precisely define and quantify those costs and the costs of control so that protection investment decisions can be made for affected facilities are needed for the Kuala Terengganu area. Those decisions have already been made for six areas-Merang, Batu Rakit, Kuala Nerus, Manir, Pengadang Buluh, Rusila in this research. Kuala Nerus and Pengadang Buluh in Kuala Terengganu are recreational resorts near Gong Merbau and Left Bank. K. Ibai. The recreational facility is a sandy beach and is quite well known with excellent name recognition. Resort housing and food stalls support the beach. Erosion has already had a high physical impact; the remaining beach is only 30 or 40 m wide. There are alternate beach areas in the vicinity and recreation could relocate to those without great inconvenience. However, Kuala Nerus and Pengadang Buluh in Kuala Terengganu are the beaches people are accustomed to using, it has established support facilities, and it is well known. These considerations combine to indicate a moderate demographic impact. Similar reasoning gives it a moderate economic impact assessment. Damages have already occurred; the beach has been diminished and a swimming pool has been destroyed. This gives the feasibility study need considerable urgency and places it in priority. The feasibility study will compare groins, breakwaters, Jetty, and beach restoration singly and in combination, to determine the most appropriate protection plan. Kuala Nerus and Pengadang Buluh (Zone 3 and Zone 5) is a big town that has a big river (Kuala Terengganu River) crossing it. This area is a low-lying and low elevation area in Digital Elevation Map (DEM) layer in GIS. The area with the highest vulnerability was identified and could serve as a guide for adaptive management planning in the future. The Nerus, Sekati, Manir, Kepung, Telemung, Cepuh, Pertang, Tersat, and Berang, Pur, Pauh, Pueh, Dura, and Sala are the primary tributaries of the Kuala Terengganu River basin, with a total catchment area of roughly 5000 km 2 . These rivers flow through a variety of socioeconomically structured regions, including agriculture, tourism, farming, aquaculture, commercial industries, urban and rural settlements, reserves, and forests. The purpose of studying rivers is to develop a field of knowledge that includes topics such as floods, erosion, sedimentation, water quality, and other studies of streams that are often researched today. In this research river class has the highest weighting of 0.281 and is very sensitive to erosion in the environment layer, whereas the geology class has the lowest loading of 0.028. The land-use class is the most sensitive and significant to erosion in the environment layer, with the greatest weighting of 0.160 after the river class. Land-use activities trigger land erosion rates and create sediment yields due to the reduction of agricultural land during historical events. Land use is a type of synthetic land surface modification that has a significant impact on environmental function and aquatic life. The kind and variables of land use will be used to evaluate the influence of river corrosion and erosion on the catchment grassland and vegetation cover. As a result, Idle Grassland, which has the greatest value of 0.233 in the land-use criterion, is essential to erosion in the Kuala Terengganu coastline area.
The soil class is sensitive and crucial to erosion in the environment layer, with the greatest weighting of 0.094, and the water body (River and coastal area) as a sub-criterion has the highest weighting of 0.299. Soil erosion is one of the factors that contribute to increased total suspended solid concentrations. This is because land erosion is one of the most difficult problems to predict, especially in tropical areas with a lot of rain. If the sediment lifting process cannot be solved, it will cause a slew of issues in the river and coastal area care and management. Furthermore, human activities such as municipal, industrial, aquaculture, and aquaculture led to the occurrence of total suspended solid concentrations. This is evident when the overall suspended solid concentration in the Kuala Terengganu River and coastal area is high.
This circumstance demonstrated the expansion of land-use development such as residential, industrial, and tourism structures surrounding the Kuaka Terengganu River and coastal area, which raises the risk of soil erosion.
Right Bank Kg. Merang, Kg. Seberang Takir, Left Bank of Kuala Marang in Kuala Terengganu are a fishing village. The primary access road to the village has been destroyed by erosion. There is an alternative path, but it passes through the village and is much less convenient than the primary road. The impact on transportation is considered moderate. Threatened houses are temporary type structures of relatively low economic value and road loss forces all traffic to and from the village to detour and experience a delay. The sand beach that once fronted the village has been eliminated by erosion. Because it had the only minor capability for recreation, the physical impact is had the only minor capability for recreation, the physical impact is moderate. The threatened dislocation of so many families and the loss of the primary access road combine to considerably affect and disrupt the established communal life of the village near Kuala Terengganu River and coastal area. Because high demographic impact and moderate economic impact are probable within 3 to 5 years, the need for the feasibility study is urgent.
Acceptable erosion areas distinguished from other eroding areas because no significant social or economic losses from erosion are expected within the next 10 to 15 m years (Table 5). Consequently, feasibility studies are not needed in the foreseeable future. Furthermore, feasibility studies to determine erosion protection needs should not ever be needed because such needs should not be allowed to develop. Unacceptable losses are not expected either because activities in the area are already adequately protected or because the area is undeveloped, not because it is free of erosion. Any development in 3 class areas must properly consider the causes and considers the causes and consequences of erosion and incorporate appropriate control measures (Table 5).

Conclusions
D-AHP Hybrid model was used to develop the CLUA analysis. Understanding and determining land-use change causal indicators is critical for developing the coastal city hazard model in the research region. Other significant variables include choosing the best local agency in the country or an international agency in another country, as well as the number and quality of experts needed to build the expert choice matrix and the best MCDA model. The criteria, sub-criteria, and alternatives were chosen based on information gathered from the literature research, different Malaysian publications such as NAHRIM (2010) and DID (2010), and foreign papers such as the Ministry of Environment (2008) and ETC CCA (2011). Others include fieldwork studies and ECM development by specialists from Malaysian government agencies such as DID, SMART, and NAHRIM.
As a result, several theme data layers relating to erosion cause vulnerability and hazard indicators, such as topography, river, slope, soil, geology, land use, road, build-up, population, flood, sea-level residual, shoreline change, wind, and rainfall, must be organized. These indexes are categorized as preparation indexes categories that affected coastal land-use change and erosion and they are considered to be responsible for vulnerability and hazards in the research location. One of the most significant aspects of this study was the selection of the best and most important criteria and sub-criteria using the AHP Hybrid design and Delphi technique to provide a useful outcome.
The understanding and determination of erosion causative layers are crucial for the study area to develop the CLUA model. Although many layers could be important concerning the land-use change occurrence for the specific coastal city region, the same factors may not be important for other regions. These layers are categorized as the preparatory layers categories that affected the land-use change directly or indirectly and they are considered to be responsible for land-use change and erosion occurrence in the study area. One of the most significant aspects of this research was the selection of the best and most important layers and classes using AHP design for coastal erosion and the Delphi technique to achieve a useful result. Combines approach through two components and presents the findings and contributions. These were the design of environmental indicator consisting of semi-quantitative data including hazard and vulnerability index in a GIS environment with a D-AHP Hybrid model and generated new coastal land-use change layers in the GIS environment. To estimate erosion risks, a good erosion model and accurate mapping of erosion-prone areas are needed. The application of D-AHP Hybrid and GIS techniques produce the final maps for important and sensitive zones of the study area. Thus, a risk assessment was performed for the vulnerability of coastal hazards due to environmental and climate change factors. The D-AHP Hybrid model was used to investigate the potential coastal city land-use change and erosion vulnerability and risk. The output of D-AHP was used as an input layer in GIS to map the hazard, vulnerability, and risk in the area. From the final results of GIS maps, we found high and low-risk areas and suitable and unsuitable areas for each of the six zones. We conclude that D-AHP and GIS techniques can improve the decision-making process.
The definition of suitable and unsuitable areas would be very important for the government to recognize which area is facing high risk or low risk for coastal land construction. For example, results show that UMT is one of the vulnerable areas to erosion and the international airport in Kuala Terengganu is at moderate risk because the airport building is in a low-lying and flooding area.
Hence the government could adopt the new policies to protect this area against coastal erosion by constructing jetties and seawalls to protect the inland area against wave action and prevent coastal erosion. The results may help managers to make better decisions and new management, as well as enable the government to plan how to solve the problem in this area. Policymakers can evaluate the sustainability of cities or near the sea, regarding how to identify the best suitable and unsuitable areas for land-use planning.
The results indicate that for land-use change and erosion vulnerability assessment, the river criterion has a highly significant effect, and the geology criterion has less effect for erosion. Conversely, built-up in human activity criterion has a high effect and road criterion has less effect on erosion. Zone 5 has highly vulnerable erosion for the coastal and city area. The shoreline change (coastal hazard) and rainfall (climate change) have a highly significant effect on erosion. Zone 5 is highly vulnerable to shoreline change, and Zones 5 and 6 are highly vulnerable to rainfall. Vulnerability is significantly affected by hazards in the coastal area of Kuala Terengganu. The Northern and Southern regions of Kuala Terengganu coastal area are comparatively less vulnerable to erosion.
Since this research has been conducted using limited coastal city data and maps, therefore, it is clear that it can only give a general view of environmental hazards for the Kuala Terengganu coastal area. However, as one of the first overall coastal management and planning studies in the Kuala Terengganu coastal area, it could indicate ways of conducting future related research. The definition of suitable and unsuitable areas would be very important for the government to recognize which area is facing high risk or low risk for coastal construction.
There are several constraints for integrated environmental geospatial modeling in coastal cities areas such as Kuala Terengganu. They are as follows: One of the main limitations in using MCDA models is the inadequacy of sources of the environmental criteria and sub-criteria for erosion on coastal areas. The limitation of a fundamental database for the environmental and coastal city in the study area. For environmental indexes, we estimated the semi-quantitative measurements from the experts in some Malaysian agencies. For environmental and climate data, requesting data and collected data were one of the main problems. They were collected from different departments and local agencies in Malaysia. Another limitation was the lack of frequently measured shoreline data sources with appropriate accuracy and time intervals as well as coastal engineering data. The short history of reliable shoreline imagery data caused long-term shoreline study to be impossible for this research. It was difficult and, in some cases, impossible to find identifiable and consistently defined shoreline indicators in aerial photographs in specific coastal sites. Since this study has been conducted using limited coastal city data and maps, therefore, it is clear that it can only give a general view of environmental hazards for the Kuala Terengganu coastal area. However, as one of the first overall coastal management and planning studies in the Kuala Terengganu coastal area, it could indicate ways of conducting future related research.
According to the findings of this study, the research indicates the need for future research for the following issues: Results of this study show that integration of MCDA models and GIS can be applied for erosion risk assessment in the Kuala Terengganu coastal area. Our study is the first attempt to use the integration of these methods in a coastal city risk assessment modeling that used different causative factors. Although the models can demonstrate the potential of using these integration methods and provide useful results in the study area, further research and more application of this methodology are required to test the transferability of the models in several criterion and alternatives with similar and different coastal hazard conditions.
For future studies, coastal environmental conservation can aggregate the environmental sustainability and vulnerability, because the hazard and vulnerability such as erosion and flood can affect the conservation and sustainability in coastal city areas. This research can also assist coastal managers, planners, and developers to identify threatened areas using the geospatial model and improve coastal land-use management. The future sustainability of the coastal city systems can then be projected using a time-series data map and satellite images for the identified indicators, the projected sustainability will be more informative. Therefore, decision-making information can be successfully derived and assist government policymaking to establish stronger sustainable development of the coastal zones. Using the spatial data and digital representation of rainfall, land use, topography, soils and other factors by GIS can help to better define the hazard and vulnerability and finally risk of erosion area in the coastal city zone. The use of these maps (suitable and unsuitable maps) enables managers and planners to improve the current policies and management of land use. It should be noted that the erosion criteria of this study should be changed in another study area because these indicators in Kuala Terengganu coastal city may not be applicable for other regions.