An integrated approach to identify suitable areas for built-up development using GIS-based multi-criteria analysis and AHP in Siliguri planning area, India

Unplanned and haphazard built-up growth due to immense pressure of population is a common problem for maximum cities of developing countries which paves to create urban sprawl by capturing spacious areas of natural environment. For this, selection of suitable places for built-up development by maintaining ecological balance has become inevitable part for proper urban planning. Siliguri planning area, which has witnessed a tremendous growth in built-up features after the 1990, has been chosen for this work to find out best places for future built-up development with minimum effects on different natural features and also calculates the amount of different types of suitability which will be transferred from different land use and land cover categories in near future which is a novelty of this work. The work has been executed through GIS-based multi-criteria analysis, and a total of nine criteria have been taken to locate best places for built-up development. Analytic hierarchy process or AHP has been used to calculate weights of each criterion by using pair wise comparison matrix. Final site suitability map has been prepared with four different classes: high suitable, moderate suitable, less suitable and restricted area and the percentage of land is high in moderate suitable category that is 52.33% which covers the intermediary portion, mainly wastelands and agricultural land of the study area with higher facility of accessibility. Results from this work would be helpful for the planners, stakeholders, policy makers by identifying best locations for the intensive developmental projects in near future.


Introduction
Urban growth is a serious and alarming issue which can be defined by change of spatial features as well as demographic structures and also refers to the increase importance of towns and cities for the concentration of population [1,2]. Simply, it is an expansion of urbanized land category that is built up by urban extension [3]. Most of the cities around the world have been facing this problem of urban growth [4] though it is more severe in developing countries like India [5,6] due to the excessive growth of urban population and economic growth which stimulates people to migrate from neighbouring rural areas to settle there and simultaneously increases households and other impervious surfaces which in a single term built up mostly in an unplanned way and leads to decrease the spacious areas of natural environment like, forests, agricultural land and also rural lands [7,8]. These unplanned, unorganized and haphazard patches of built-up growth pave the way to create urban sprawl [9].
Due to this unplanned growth of built up or urban sprawling and unprecedented population growth and also to sustain ecological balance, it is necessary to find out suitable areas for the new built-up development [10] and proper use of different land categories. Suitability analysis by the help of GIS platform is a process which aims to find out the best locations or areas for the proposed work while considering the ecological sustainability [4] and works according to the need of stakeholders [11]. Land suitability analysis for future urban growth is considered to be the best and effective method to find out suitable areas which works with different parameters or criteria that influence on one to settle in a place and which have some certain weightage [12]. It has become significant because this technique helps to reduce the time as well as effort to manually find out a suitable place and also filtering out the least suitable or not suitable places for the construction of new built up [13]. As this analysis is based on different criteria or alternatives and selection of criteria is based on human judgement, so it is called multi-criteria decision-making (MCDM) analysis and this MCDM analysis is capable of providing different tools for complex planning in terms of strutting priorities, weighing and selecting alternatives and also allocate resources across competing activities [14]. This MCDM can be classified into two categories, based on the number of alternatives under consideration and these are, Multi-Attribute Decision Making or MADM and Multi-Objective Decision Making or MODM [15][16][17]. MADM is mainly used for suitability analysis of a particular facility, and alternatives are predefined and limited in number [18] like, to find out a place for solid waste disposal [19,20], for locating new hospitals [21], etc. But in the MODM technique, alternatives are not predefined and are numerous in number [18]. In multi-criteria decision-making technique, weight or scoring of a particular criterion plays a significant role as they provide the importance of a particular criterion [22]. To determine the weights of criteria, different models can be grouped into two and these are subjective model and objective model. In subjective approaches, decision makers or experts on that respective field give their opinion on the significance of criteria and the weights of criteria are calculated based on the information gathered from the experts or decision makers, whereas in objective methods, weights of criteria are determined based on the information contained in a decision-making matrix by applying certain mathematical models and it neglects different opinions of experts or decision makers [23]. Some well-known objective approaches are, Entropy method, Criteria Importance Through Interaction Correlation or CRITIC method, FANMA method, etc., and some best subjective methods are Decision-making Trial and Evaluation Laboratory or DEMATEL method, Step-wise Weight Assessment Ratio Analysis or SWARA method, Best Worst Method or BWM, Full Consistency Method or FUCOM, Analytic Hierarchy Method or AHP. DEMATEL has been used in different kinds of studies [24,25], but its main disadvantage is unable to validate the result, i.e. the measurement of inconsistency and for this DEMATEL method is often used to find out the weights of criteria with combination of Analytic Network Process or ANP [26] which can slightly reduce its weakness. SWARA method has also simplicity of calculation with few numbers of steps, but like the DEMATEL method, it also does not have the ability to determine the degree of consistency of the pair wise comparison [23]. In recent, AHP, BWM, LBWA and FUCOM have been increasingly applied in different works because of the ability of determining the degree of consistency. Among these methods, required number of pair wise comparison in FUCOM is minimum, because it requires only (n-1) of criteria, whereas AHP requires (n(n−1)/2) criteria. But the further process of calculation in FUCOM method is complicated. BWM method requires only (2n−3) number of pairwise comparison which is more than the half number of pairwise comparison that requires in AHP and also more than the FUCOM. However, a large number of comparisons in the pairs of criteria define the limitations for solving nonlinear model make the application of the BWM significantly more complex [27]. Level-Based Weight Assessment or LBWA model also allows for the calculation of weight coefficients with a minimum number of criteria in pairwise comparison, that is (n−1) as FUCOM method [27]. Another advantages of LBWA method are the algorithm does not become more complex when criteria increases and weightage values are obtained by simple mathematical calculations which eliminates inconsistencies of expert preferences, which are tolerated in BWM method [27]. Beside these advantages, this method needs to emphasize the flexibility in terms of additional corrections of weight coefficient values by the elasticity coefficient which allows decision makers to further adjustment in weight coefficient values in accordance with their own preferences. Saaty [28] proposed the Analytic Hierarchy Process or AHP to find out the weights of a particular criterion in multi-criteria analysis by pair wise comparison matrix table, and it is based on hierarchical structure, and for this each and every criterion can be focussed in a better and transparent way [16]. It is a widely used method along with GPS and remote sensing in GISbased MCA technique for proper urban land use planning and management [18,29]. AHP can set the criteria of land use suitability planning in a structured way based on available data, and it can standardize the criteria with simple mathematical calculations [18], and it is much easier for the social science researchers as it follows few numbers of steps to calculate weights of a criterion and does not involve any complex mathematical calculation [16]. Finally, GIS-based AHP can overlay all the standardized criteria that help to make rankings on a locational basis [18]. Based on these considerations, this AHP method has been chosen for this work to calculate weights of different criterion.
In recent years, Multi-Attribute Decision-Making (MADM) analysis and Analytic Hierarchy Process or AHP have been used in various literatures for the demarcation of suitable places for new built-up development. For example, Aburas et. al. [4] explained a review on different literature which deals with land suitability analysis for further urban growth by AHP and GIS techniques. They categorized all criteria into five classes and presented the whole methodology step by step. Ullah and Mansourian [18] used multi-criteria decision-making model and AHP for the evaluation of suitable land for urban land use planning. They took different factors, such as physical, accessibility and also community factors for the analysis. They also described that the justification of taking AHP method for the weightage calculation of these factors. Weldu and Deribew [30] developed a suitability map for the identification of potential sites for the housing development by using multi-criteria evaluation technique along with AHP and GIS. They took a total of ten criteria and presented the whole calculation of AHP in a lucid way. Chougale et. al. [10], Kumar and Shaikh [13] also used GIS-based multicriteria evaluation technique along with AHP method for the identification of suitable sites for the urban development and presented the whole calculations of AHP in an efficient way. For the pairwise comparison of AHP, it is necessary to calculate the weightage on the basis of information that gathered from different experts as AHP falls in the subjective model, but most of the authors have ignored this step and they put the value of Saaty's nine point scale according to their opinions. This present work has been carried out using AHP method and has been followed the decisions that have been put by the experts. This work has been carried out in Siliguri which has witnessed a tremendous growth in built-up features due to huge migration in the partition time and after the partition it got its importance due to its location as it serves as a transit point for the entire north eastern part of the country and some neighbouring countries and also as an important economic centre. For this, this city attracts population from different states for last four decades. According to census of India, 2001 as a total of 139,418 persons from 1991 to 2001 have migrated from different states and districts of West Bengal which was 106,098 from 1981 to 1991 (Census of India, 1991, 2001, Migration table, D series data). This huge growth of population give rises to unplanned built-up growth in and around the city and transforms different productive land into built-up and impervious land which need immediate actions to control this unscientific and haphazard growth. No such work to overcome this problem and also to identify the best locations for scientifically built-up development has not been carried out for this area. So, for the sustainable urban development or built-up growth, it is needed to improve the current pattern and for this reason, and major thrust or aim of this work has been given to find out best places for the built-up development with the minimum effects on different natural features using GIS-based multi-criteria analysis and analytic hierarchy process in and around the city Siliguri (Fig. 1). This work has been done with the consultation of different experts who are responsible for this type of developmental works and based on their advices and thinking, and this work has been preceded further which is a novelty of this work. This work also calculates the amount of different types of suitability which will be transformed from different land use and land cover categories in near future. At last, the calculated result of highly and moderate suitable places has been verified by physically. The whole paper has been arranged in a following way. At first, criteria have been chosen and for the criteria weight, AHP method has been followed in where pair wise comparison table has been filled up with the consultation of experts. Data for each criterion map have been collected from different repositories and prepared in GIS platform, and then calculated weights have been assigned in each criterion map. Final suitability map has been prepared by integrating all criteria maps by using weighted overlay tool in GIS platform (Fig. 2).

Study area
Siliguri, one of the most rapidly developing and fastest growing metropolis, is located in the northern part of West Bengal and foothills of Darjeeling Himalaya. It was a small village with less than 800 populations in 1901 and turned into a municipality town in 1949 with a population of 32,480. It is because of the partition of India after independence and huge in-migration from the Bangladesh (erstwhile East Pakistan), Nepal, southern parts of Bhutan and adjacent states of it. Due to different developmental projects and business activities, it increases its population by natural way and by pulling from neighbouring areas; it got a status of Municipal Corporation in the year 1994. In the year 1980, Siliguri Jalpaiguri Development Authority (SJDA) was set up by the government of West Bengal with an aim of planning and development in Siliguri and Jalpaiguri planning area. Due to rapid growth of population and built-up area that have been experienced by Siliguri after 1980, a need for the control of growth and proper land use was felt by the authority and for this Outer Development Planning Area (ODPA) was formed in the year 1992 for three separate area, Siliguri, Jalpaiguri and Naxalbari; but from the year 2001, this ODPA term has changed and a new term Land Use and Development Control Plan or LUDCP has introduced. Out of these three LUDCP or planning area, Siliguri Planning Area has been selected for this work as this area has witnessed huge changes in built-up pattern, which are not in a proper planned way. At present, this area covers Siliguri Municipal Corporation, twelve census towns and eighty eight villages with a total population of 963,498.

Selection of criteria
Appropriate and significant criteria selection is the first step to calculate the suitable area properly for any purpose though the final selection of criteria also depends on the availability of data [4]. Selection of criteria for the suitable area of built-up growth varies from researcher to researcher, yet all the criteria can be grouped into four major factors, i.e. accessibility, physical, socio-economic and environmental factor. Table 1 shows the factors and criteria that are used by different researchers in previous studies and based on this table and also synthesizing different literatures [11,12,18,30] most significant and frequent used criteria have been taken for this work.
Accessibility and physical factors are most determining factors followed by environmental and socio-economic factors as these have been used by maximum researchers. Three accessibility factors (distance from road, distance from airport and distance from railway station), three environmental factors (land use land cover types, distance from river and distance from canal), two physical factors (slope and geological formation) and one socioeconomic factor (population density) have been taken in this work for the identification of suitable site for built-up growth ( Table 2).   Step-wise flow chart of GIS-AHP-based methodology to identify suitable places for built-up development

Calculating weightage by AHP
AHP or Analytic Hierarchy Process is a well-known and universally used subjective method to determine weights of criteria in multi-criteria decision-making technique [22]. It can be used to find out weights for qualitative as well as qualitative data [40]. It is based on questionnaire survey by expert. Two of them are professor and assistant professor of university and specialized in urban geography and regional planning and they have the experience of doing different future prospective projects for this area under the Siliguri Jalpaiguri Development Authority, one is real estate developer, who already have done several projects related to apartments and housing complexes, and rest two are local people who want to buy residential houses in this region. They are asked to give a score on the basis of Saaty's 1 to 9 ratio scale [38] through pair wise comparison.
The comparison between two criteria is made by judgement of an expert that represents how much more one criterion influences over another [11]. If two criteria influence same or have equal importance, a score of 1 have to give and a score of 9 is given when the influence of one criterion is extreme over another criterion in the comparison table [22]. After normalizing each value by dividing the actual value to its sum of column value, weights are derived by arithmetic mean method (Table 3). After calculating weights, it is necessary to calculate consistency ratio (CR) to check the degree of consistency of judgements that has made by any expert [30]. Saaty [41] has provided an index to justify that the comparison matrix is consistent or not. The formula of getting CR value is,  Very strong to demonstrated importance An activity is favoured very strongly over another, its dominance demonstrated in practice 8 Very, very strong 9 Extreme importance The evidence favouring one activity over another is of the highest possible order of affirmation Reciprocals of above If activity i has one of the above nonzero numbers assigned to it when compared with activity j, then j has the reciprocal value when compared with i A reasonable assumption

Rationals
Ratio arising from the scale If consistency were to be forced by obtaining n numerical values to span the matrix where, CR is the consistency ratio, CI is the ratio of consistency index and RI is the random inconsistency index of randomly generated pair wise comparison matrix of order 1 to 10 calculated by approximating random indices using a sample of 500 [39]. In the formula of consistency ratio (CR), CI is calculated by the following equation.
where, CI is the consistency index, λ is simply the average of consistency vector, n refers to the number of total criteria. The consistency ratio has been designed in such a way that, if the value is < 10%, then the matrix is said to be a sufficiently consistent and the matrix is accepted, but if it is more than 10%, then it is an indication of inconsistent judgement (Table 4).

Preparing criteria maps
After selecting criteria for suitable site of future urban area, different thematic maps have been prepared. Total nine thematic maps have been prepared in a GIS platform as the selected criteria are nine which are mostly in vector data format.

Rasterization of criteria maps
As the mostly maps are in a vector data format, so it is the next step to convert the format of vector maps into raster for further analysis.

Integration of criteria maps and final suitability map
After converting into raster data format and assigning weights, these criteria maps are then integrated by using weighted overlay tool in GIS platform to get the final suitability map of future urban development.

Results and discussion
As almost all the cities in developing countries are expanding haphazardly, so the identification of suitable sites for proper sustainable urban growth has become a fundamental and important work for these cities. The further analysis for suitable site has been carried out by using nine parameters, and the amalgamation of these nine criteria has helped to find out different suitable places for future built up.

Accessibility factor
Smooth and efficient movement and better connectivity with the entire city and other places are the prime need for an individual to set up new construction either for residence or commercial purposes. In this regard, proximity to road and railway play an important role though proximity to airport also influences people who frequent use this service. Construction of new road is expensive, and government provides new roads to a place which has an ample number of settlements, so people wants to settle in a place which has already better road connectivity.
Similarly as the number of railway stations, metro stations, airports in a city are fixed, so people tend to set up their dwellings or commercial unit close to these facilities for easily accessible.
In this work, to find out the accessibility in the entire area, major roads like, Asian highway, national highways, state highways and important village and municipal roads which are more than 15 feet width and has connection with different places, six railway stations and one airport have been digitized from the Google earth. Now to find out better accessibility and also to find out the best locations to the existing roads, railway stations and airport, buffer zone maps by taking 250 m as interval has been prepared (Fig. 3) and their importance for builtup growth has been identified by ranking.

Environmental factors
Environmental factors play a significant role for identifying suitable places for any developmental work. Though it is not possible all the time, but mostly all the constructional projects have been implemented by preserving the natural environment. In this work, three environmental factors, such as, land use land cover types, distance from river and distance from canal, have taken to identify the suitable places.
Land use land cover map has been prepared from the Landsat 8 data of the year 2019. It has been categorized into built-up land, water bodies, vegetation, agricultural land and wasteland based on NRSA manual, 2006. In the further analysis, existing built-up land has been considered as restricted category because no such new construction can be done in future, though vertical growth can take place but this work mainly tries to find out the places which are vacant. Similarly, water bodies are also categorized as restricted land, whereas wasteland and agricultural land have been given the priority for locating suitable places of built up. Siliguri and the surrounding study area lie in the piedmont of Darjeeling Himalaya, so the numbers of rivers are naturally more and are perennial in character. This area also receives huge rainfall during the monsoon, which makes these rivers more furious, but the possibility of occurrence of flood is less in here because this area located in the piedmont region with an average height of 150 m. Two major rivers Mahananda and Balason are flowing through this area. For identifying the best location, buffer zones have created with an interval of 250 m from existing water bodies with a ranking from 1 to 5. The area which are just beside the river ranked as 5 which means very less suitable, and the area located beyond 1000 m ranked as 1 means very high suitable because people wants to reside in a relatively higher area with a minimum distance from river to avoid the flood problem. Similar ranking pattern has assigned to the layers of canal map, which is located to the southern part of the study area (Fig. 4).

Physical factors
Physical factors, like elevation and slope, are the basic parameters for setting up settlements mainly in the hilly terrain. Steep slope lands are not suitable for any kind of construction because it needs more construction costs and also risky. Therefore, gentle slopes are suitable for construction. Though in plain region these parameters have negligible influences on constructing new built up. For the major built-up projects like, shopping malls, hotels, townships, high rise buildings and flats, geological set up of that area should be examined properly and it is necessary for this study area because this area falls in the high risk zone of earthquake according to the seismic zones of India and this area has shown a tremendous growth in high rise buildings, townships and shopping malls.
The whole study area has been categorized into five slope zones starting from less than 10 to more than 25. Classified map of slope shows that, maximum area falls in the less than 10 slope zone, whereas only banks of the river show more than 25 slope zone. Less than 10 slope zone are considered as gentle slope [42], and it is more suitable for the new construction. So this category has been assigned as rank 1 means very high suitable and more than 25 as 5 means very less suitable. Geological data have been collected from Geological Survey of India which defines that five types of formation have been seen in this area and these are Present day deposits, Jalpaiguri formation, Baikunthapur formation, Duars formation and Shaugaon formation. Baikunthapur formation has covered the most portion of this area, and it consists of very fine white sand inter layered with silty clay and overlain by dark grey to thick silty loam which is ideal for new construction and for this it has been ranked as 1, whereas present day deposits are the youngest formation and mainly found in the river bed consists of sediments which is not suitable for the construction of new built up and hence, it is ranked as 5 means unsuitable (Fig. 5).

Socio-economic factors
Human being always finds a place to reside which is socially favourable for them. Most of the time, same social group of people resides closely to each other for the similarity of culture. But it cannot be possible all the time in rapidly growing cities. In these type of cities, maximum people wants to buy a plot at a cheaper rate which mostly lies away from the city core, whereas few people wants to buy in the city for any price to avail the civic facilities. But for the overall people, it can be said that, land price is a determinant factor to set up new construction of housing. Beside these, population density, housing density are also important factors which play a significant role in new built-up development, because most of the people wants to settle in a low or medium dense area than in a clumsy environment.
This study area shows a huge growth of population in last two decades. It was 510,137 in the year 1991 which has turned to almost double in the year 2011, that is 963498 according to Census of India, 2011. Administrative unit wise population density varies from 44 persons per sq. km. in Tarabarir chat to 16,394 persons in Matigarahat, whereas household density ranges from 13 houses per sq. km. to 3112 houses which defines a huge variation and unbalanced distribution of built up in the entire area. For this work, the whole area is divided into five groups by population density (Fig. 6) which starts from less than 1000 persons to more than 4000 persons. As most people wants to reside in a low dense area, so less than 1000 category has been assigned as rank 1 means very high suitable and more than 4000 as rank 5 means very less suitable for creating new built up (Fig. 6).

Pair wise comparison matrix and normalized pair wise comparison matrix
As we have taken nine criteria for the analysis, so, a total of 36 questions have been formed in the questionnaire. Then, the questionnaires have been filled up by five aforesaid experts following Saaty's 9 point scale. Most appropriate view has been taken for the further analysis. In the pair wise comparison matrix table, 81 boxes are there and we have filled up 36 boxes according to the value that we get from the questionnaire. Diagonal value always will be 1 because the influence between the same criteria will  (Table 5). Normalized pair wise comparison matrix table is the second step after pair wise matrix table to get the final weightage of each criterion. Normalized value has been calculated by dividing each value of pair wise matrix table with their respective total column value. Then to find out weight of each criterion, arithmetic mean method has been used row wise (Table 6). Road is the highest influencing criterion as the weightage value is 29.7% followed by airport (18.6%), railway station (13.5%) and land use land cover (13.2%), whereas geology has the least influence over people to set up their dwellings as it has scored only 1.7% weightage. It is also clear from the normalized pair wise matrix table that, accessibility factor is the significant in the decision making of new built-up construction ( Table 6).

Calculation of the consistency ratio
It is necessary to check the degree of consistency of the judgements. Before going to calculate the value of CR, the value of consistency index (CI) and λ should be calculated. Table 7 shows the calculation of consistency vector.

Calculation of lambda
The result of consistency ratio (CR) is less than 10%, so it can be said that the judgement is consistent and it is allowed for the further analysis.

Analysing final suitability map
After calculating final weights of each criterion and preparing thematic and buffer maps in GIS platform, all map layers except land use land cover and slope map are then converted to raster data format by using conversion tool for further suitability analysis. These maps are then integrated, and reclassification of the sub-layers of each map has been done as reclassification tool is employed to convert 5on a scale [43] for land suitability (Table 8), and using weighted overlay analysis tool, suitable sites for future built-up development have been delineated. Final suitability map has been classified into four groups, restricted or not suitable, highly suitable, moderate suitable and less suitable. Table 9 shows the area that has fallen in each zone of suitability. About 31% area that falls on restricted zone means no such new construction can be done in this zone which comprises existing built-up land and water bodies. Highly suitable zone is relatively less among different zones, which is only 1.76% or 4.58 sq. km. because a little amount of land has fulfilled entire criteria which also demonstrate that this area is now at best position by its available built-up features. Further construction in this area will harm the productive features like agricultural land, planted vegetated area mainly tea plantation area because this entire area is mostly dominated by moderate suitable zone for built-up development, and this category mostly covered by agricultural land. It can be predicted that if this present trends of growth of population and built-up continue, then in future maximum portion of agricultural land, mainly tea planted area will be vanished. Major roads are relatively less than the requirement and railway stations also few in number and not distributed equally to the entire region, which is also a reason for having less amount of highly suitable area and high amount of moderate suitable area. Maximum villages of Bagdogra area have shown highly and moderate suitable CI = ( − n) ∕ (n − 1) = (10.13 − 9) ∕ (9 − 1) = 0.14 CR = CI ∕ RI * 100 = 0.14∕1.45 * 100 = 9.82 % Table 7 Computation of consistency vector LULC Land use land cover, PD Population density, RV River, RD Road, RS Railway station, AP Airport, SP Slope, GF Geological formation, CN Canal Criteria Weighted sum vector  Bagdogra are nearer to the Bagdogra railway station and the national highways, but these are too much congested by its population and different impervious surfaces which makes these places moderately suitable for new built-up construction. But for the villages that are in the moderately suitable category are far apart from the existing facilities, like, railway stations, roads and are nearer to the forest and river, but the density of built-up category and population are moderate, which makes these villages moderately suitable for the new built-up construction. Maximum villages and census towns of Matigara block are in the category of moderately suitable. Reasons behind this are villages and census towns which are nearer to the national highway are already developed by impervious categories and the density of these areas are too high, that is more than 4000 persons per sq. km. and the villages which are at the intermediary part or in the edge of this block, lacks to the road facility and railway stations and also covered by tea gardens or other agricultural lands which makes these villages moderate or less suitable for the new built-up construction (Fig. 7). Land use land cover and proximity to road are two important factors to delineate suitable places for builtup development. Table 10 shows the distribution of area of each suitability zone that falls in different LULC types. About 99% area of water bodies has fallen in restricted zone as this category of LULC has been considered as builtup free zone and it should not be captured. About 80% of agricultural land has fallen in the moderate suitable zone and only 3.5% area of it is in the highly suitable zone. Same thing has happened in case of wasteland categories. In this, about 68% area is in the moderate suitable zone and only 1.3% area is in the highly suitable zone. From Table 11, it can be said that in the first zone, that is the area which is in the less than 250 m from road, about 55 and 3.4% area of it falls in the moderately suitable and highly suitable category, respectively. These two areas can be given priority to set up new construction either for residential or for commercial purposes. About 41% area of this zone is in the restricted category as this area is already developed by built-up land. A substantial area of next two zones of roads that are 250 to 500 m and 500 to 750 m falls in the moderately suitable category. This area can be used mainly for residential uses, as the land price would be cheaper than the first zone. From the fourth zone, the amount of suitable area decreases because of increasing distance.

Conclusion
In this present growing world, haphazard and unscientific rapid growth of built up is common for maximum cities in the developing nations. So, suitability analysis or suitable site selection for proper urban development has become necessary to overcome the problem of haphazard growth within limited area. To solve this, GIS-based multi-criteria analysis with analytic hierarchy process is simple and cost efficient technique which can identify various suitability categories within a selected region. Siliguri and its surrounding region also face the same situation due to unprecedented growth of population, and for this, this area has been chosen to find out best places for future built-up development with the consultation of experts and also calculates the amount of different types of suitability lands which will be transferred from different land use and land cover categories in near future which is a newness of this work. This study reveals that only 1.76% area has been identified as highly suitable area which has located sparsely with a great distant from the core city, whereas maximum area has been categorized as moderately suitable land which is 52.33% for new construction. The maximum portion of highly suitable area will be come from agricultural land that is 4.08 km 2 and the area which is beside the existing road that is 4.37 km 2 which reveals that most of the peoples tend to reside just beside the road for easily accessible. Moderately suitable category will also come from agricultural land (79.8%) followed by wasteland that is 68.2% and for the road criteria, maximum moderate suitable part will come from 250 to 500 m category (64%) followed by < 250 m category that is 54.8%. This area will be mostly useful for residential purposes as the price of this category will be less than the highly suitable category. Though AHP method has been used in different literatures for the calculation of weightage of different criteria, but it has some limitations which have been experienced while doing this work. It is almost impossible to do with AHP when the number of criteria is above nine because for the pairwise comparison it requires n(n-1)/2 number of comparison. So, we have to limit our criteria at nine, though it can be solved by dividing the criteria into sub-criteria, but after dividing it complicates the model. Another problem is that, we cannot involve more decision makers or experts for the fulfilling of pairwise comparison matrix, as more experts can make the matrix more complicated in assigning weightage values. So we have bound to take only five experts of different fields in this work. Beside these, this type of regional planning work would be helpful for the planners, policy makers, stakeholders, etc., which allows them to identify areas where intensive developmental projects can be implemented in a lucid and efficient way.
Acknowledgement The authors are thankful to the Editor-in-chief and three anonymous reviewers for their constructive suggestions on the manuscript. First author is also thankful to his friend Mr. Kunal Chakraborty for his constant support and encouragement for working on this topic Author contributions AS contributed to conceptualization, methodology, draft writing, software. RR contributed to conceptualization, supervision, formal analysis, editing. Data availability All the data and materials that have been used in this article are openly available in different repository.

Compliance with ethical standards
Conflict of interest We do not have any conflict of interest regarding the research, authorship and publication of this article.
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creat iveco mmons .org/licen ses/by/4.0/.