Landfill Site Selection Using MCDM Methods and GIS in the Sulaimaniyah Governorate, Iraq

A shortage of land for waste disposal is one of the serious problems that faces urban areas in developing countries. The Sulaimaniyah Governorate, located in the north of Iraq, is one of the major cities in the Kurdistan Region of the country, covering an area of 2400 km2 with a population of 856,990 in 2016. Currently, there is no landfill site in the study area that meets scientific and environmental criteria, and inappropriate solid waste dumping is causing negative environmental impacts. The process of landfill site selection is considered a complex process and is restricted by numerous factors and regulations. This paper proposes multi-criteria decision-making (MCDM) methods in a model for landfill site decision. The model assumes the input of two groups of factors that need to satisfy the optimal values of weight coefficients. These groups of constants are natural factors and artificial factors, and they included thirteen selected criteria: slope, geology, land use, urban area, villages, rivers, groundwater, slope, elevation, soil, geology, road, oil and gas, land use, archaeology and power lines. The criteria were used in the geographic information system (GIS), which has a high capacity to process and analyze various data. In addition, multi-criteria decision-making (MCDM) methods followed by the weighted linear combination (WLC) method were used to derive criteria weightings using a matrix of pair-wise comparison. Finally, all the multi criteria decision methods were combined to obtain an intersection of the suitability index map for candidate landfill sites. Seven appropriate sites for landfill were suggested, all of which satisfied the scientific and environmental criteria which were adopted in this study.


Introduction
Municipal solid waste (MSW) is the main by-product of urban life. Approximately 3.5 million tons of MSW are generated daily on a global scale [1]. Population growth, improving standards of living after economic recovery, and industrial activities are all major reasons for a significant increase in solid waste quantity output in urban areas [2].
The problem of solid waste is very serious in developing countries. About 80% of the world population are living in these countries [3]. Due to the situation of these countries, such waste is affecting humans and their environment. For this reason, legislations and laws have been developed since the mid-twentieth century to ensure sustainable development and management of waste in a prudent manner to minimize the effect on human and environment [4]. This article offers a method

Materials and Methods
To present a suitable site selection landfill map, spatial analysis tools in GIS were used to prepare thirteen different criteria as layers covering the study area, with each criterion consisting of a database of a digital map prepared within GIS. All processed datasets used to input the AHP model of landfill dumpsite followed these principal steps: 1 Development of a digital GIS database (spatial information) and creation of mapping limitations using regions of exclusion. 2 Creating appropriate buffer zones or special limitations around significant regions to fit each criterion map. 3 Using MCDA models to specify the most preferred alternative and determine each criterion weighting according to the priority of the study objective. 4 Weightings of sub-criteria based on literature, expert opinion, environmental guidelines and regulations. 5 Overlap of all weighted criteria resulted from (AHP, SRS, RSW and SAW) methods in GIS to prepare a suitable index map and indicate the priority for optimizing sites for landfill siting.

The Hierarchical Scheme for Landfill Siting
The hierarchical framework for decision-making was based on expert views in multi-criteria analysis, prior research and available information in the study area [33]. The hierarchical scheme of the site selection involves three levels as presented in Figure 2.
The first level in the proposed model represents identification of the factors affecting the potential landfill site model, which consists of two categories: natural environmental factors, and artificial factors. The second level contains six primary criteria: land, topographical, hydrological, geology, infrastructure, accessibility and social-cultural factors. The third level involves all thirteen criteria that were used in the present research to determine suitable sites as the optimum solution for landfill solid waste disposal.
The hierarchical scheme was designed according to the relations between each criterion and the decision process for landfill site selection.

Materials and Methods
To present a suitable site selection landfill map, spatial analysis tools in GIS were used to prepare thirteen different criteria as layers covering the study area, with each criterion consisting of a database of a digital map prepared within GIS. All processed datasets used to input the AHP model of landfill dumpsite followed these principal steps: 1.
Development of a digital GIS database (spatial information) and creation of mapping limitations using regions of exclusion.

2.
Creating appropriate buffer zones or special limitations around significant regions to fit each criterion map. 3.
Using MCDA models to specify the most preferred alternative and determine each criterion weighting according to the priority of the study objective. 4.
Weightings of sub-criteria based on literature, expert opinion, environmental guidelines and regulations.

5.
Overlap of all weighted criteria resulted from (AHP, SRS, RSW and SAW) methods in GIS to prepare a suitable index map and indicate the priority for optimizing sites for landfill siting.

The Hierarchical Scheme for Landfill Siting
The hierarchical framework for decision-making was based on expert views in multi-criteria analysis, prior research and available information in the study area [33]. The hierarchical scheme of the site selection involves three levels as presented in Figure 2.
The first level in the proposed model represents identification of the factors affecting the potential landfill site model, which consists of two categories: natural environmental factors, and artificial factors. The second level contains six primary criteria: land, topographical, hydrological, geology, infrastructure, accessibility and social-cultural factors. The third level involves all thirteen criteria that were used in the present research to determine suitable sites as the optimum solution for landfill solid waste disposal. The hierarchical scheme was designed according to the relations between each criterion and the decision process for landfill site selection.
The "river" criterion was adopted in order to protect surface water from organic, inorganic and microbial pollutants and contamination from leachate leaking from the solid waste landfill, which is a very common phenomenon in rapidly developing countries [34,35]. Distances from roads to the selected landfill location must ensure proximity to the road network with a view to facilitating transport and thus reducing the relative costs. Moreover, access to a landfill site should be adequate to ensure that the aesthetics of the region are not adversely affected [16,36]. The criterion of topography (elevation) was selected in respect of the risk of flooding in the area [37]. Steep slopes would make it difficult to transport the waste to the site and cause leachate drainage into flat regions and water bodies that cause contamination, while slopes that are too flat would affect the runoff drainage [16,23,38].
For the "urban area" criterion, an approximate distance of solid waste sites depends on many factors, such as governmental regulations and environment suggestions from scientists. Generally, the approximate distance should be appropriate for urban planning and potential future expansion, as well as ensuring impacts on public health, the aesthetic aspect of the city. Furthermore, a landfill site should be at a proper distance from the boundaries of urban regions due to the economic factors of the transport costs [18,[39][40][41].
The "soil" criterion was used to preserve the study region from contamination by leachate transportation through the unsaturated zone within soil layers. The soil of the study area consists of four types of soil which are known as alluvial and recent deposits [42,43].
The "power lines" criterion was selected to prevent any danger to and from high voltage power lines and to ensure minimal disruption to the infrastructure near the solid waste site [44]. The "archaeological" criterion is the national cultural heritage, which includes numerous archeological and historical sites. In the current study, the criteria "oil and gas field" were adopted to protect important infrastructure, which is considered the most valuable source of national wealth [45,46]. "Geological formations" describes the outcropped and subsurface lithology. In the current study, this criterion was adopted to protect the study area from geological hazards by understanding the engineering properties of the formations for settlement and other large structure building and by finding economic and valuable formations like groundwater and oil accumulation, together with raw material for cement, gypsum and aggregates [42].
The "land use" criterion was adopted according to the land use classification and land specification due to priority of land use and the distance between the disposal site and factories, pasture and forests, also considering the cost of unused land and ownership of agriculture land [5]. The "river" criterion was adopted in order to protect surface water from organic, inorganic and microbial pollutants and contamination from leachate leaking from the solid waste landfill, which is a very common phenomenon in rapidly developing countries [34,35]. Distances from roads to the selected landfill location must ensure proximity to the road network with a view to facilitating transport and thus reducing the relative costs. Moreover, access to a landfill site should be adequate to ensure that the aesthetics of the region are not adversely affected [16,36]. The criterion of topography (elevation) was selected in respect of the risk of flooding in the area [37]. Steep slopes would make it difficult to transport the waste to the site and cause leachate drainage into flat regions and water bodies that cause contamination, while slopes that are too flat would affect the runoff drainage [16,23,38].
For the "urban area" criterion, an approximate distance of solid waste sites depends on many factors, such as governmental regulations and environment suggestions from scientists. Generally, the approximate distance should be appropriate for urban planning and potential future expansion, as well as ensuring impacts on public health, the aesthetic aspect of the city. Furthermore, a landfill site should be at a proper distance from the boundaries of urban regions due to the economic factors of the transport costs [18,[39][40][41].
The "soil" criterion was used to preserve the study region from contamination by leachate transportation through the unsaturated zone within soil layers. The soil of the study area consists of four types of soil which are known as alluvial and recent deposits [42,43].
The "power lines" criterion was selected to prevent any danger to and from high voltage power lines and to ensure minimal disruption to the infrastructure near the solid waste site [44]. The "archaeological" criterion is the national cultural heritage, which includes numerous archeological and historical sites. In the current study, the criteria "oil and gas field" were adopted to protect important infrastructure, which is considered the most valuable source of national wealth [45,46]. "Geological formations" describes the outcropped and subsurface lithology. In the current study, this criterion was adopted to protect the study area from geological hazards by understanding the engineering properties of the formations for settlement and other large structure building and by finding economic and valuable formations like groundwater and oil accumulation, together with raw material for cement, gypsum and aggregates [42].
The "land use" criterion was adopted according to the land use classification and land specification due to priority of land use and the distance between the disposal site and factories, pasture and forests, also considering the cost of unused land and ownership of agriculture land [5]. Groundwater ranges in depth from 19 to 280 m below the ground surface in most regions of the Sulaimaniyah Governorate. The landfill location should be located in an area with a deep groundwater depth in order to decrease the risk of leachate from solid waste and protect water aquifers from pollution [18,34]. Groundwater ranges in depth from 19 to 280 m below the ground surface in most regions of the Sulaimaniyah Governorate. The landfill location should be located in an area with a deep groundwater depth in order to decrease the risk of leachate from solid waste and protect water aquifers from pollution [18,34].

Dataset Criteria Map
In this study, thirteen data sets obtained from international organizations and various government sources were used as input parameters for AHP modelling. The first source is available government data, including urban area, archaeological sites, power lines, oil and gas fields, villages, soil and geological formations. These digital data were converted into format files using GIS.
The river, road, and elevation data from the United States Geological Survey USGS Earth Explorer and the elevation map were downloaded as a digital elevation model that was converted to a slope map. The source of the data is presented in Table 1.
Water level depth data was obtained from the Sulaymaniyah groundwater authority, and 235 wells were added to GIS using the spatial analysis tools "Kriging" to create a map of the groundwater table in the Sulaymaniyah Governorate. The land use classification was prepared using satellite data and processed by remote sensing software (ENVI 5.4).

Criteria Restriction
Specific geographic features created using buffer zones around each criterion using GIS spatial analysis tools. A buffer zone is a neutral area that can be separated by grade to decrease or eliminate the impact of land use practices on delicate regions or natural characteristics. To determine the minimum and maximum distance from landfill locations, each criterion should be characterized by the recommended distance according to the potential environmental risks, human health and excess cost, and by taking into account the requirements of government regulations [47,48]. As shown in Table 2, the buffer areas were created at a distance from each feature to the required criteria based on a literature review.

Sub-Criteria Rating Values
Each criterion was classified into sub-criteria in this research and assigned a suitability rating value from zero to ten [58]. Due to the associated legislation, limitations and regulations, the criteria and variables are determined and assessed for their suitability for the research region, together with experience, information from literature reviews and scientific experts in this field. This was done in order to collect more information about category priorities for the study area and to rate the categories' importance in respect of solid waste site selection.
A number of steps were followed to evaluate the rating value for each criterion and sub-criteria using the GIS spatial analysis tools in a sequence (Buffer, Clip, Extract, Overlay, Proximity, Convert, Reclassify and Map Algebra). Sub-criteria buffer zone and rating values for the input layer are shown in Table 3.  In this revise, a distance of less than 1000 m from the boundaries of a river was given a rating value of 0 to reduce the potential for river contamination from landfill, and a rating value of 10 was provided at distances higher than 1000 m ( Figure 3A). Main roads and highways were integrated into the "roads" criterion. The greatest score of 10 was allocated to buffer areas from roads to landfill locations of 1000-2000 m. Buffer zones of less than 500 m and those of more than 3000 m were given grades of 0 and 3, respectively. Buffer areas of 500-1000 m were allocated a grade of 7 depending on the aesthetic effects on the region, while buffer zones of 2000-3000 m were given a grade of 5 to reduce transportation costs ( Figure 3B).
For the elevation layer, the digital elevation model in the form of the raster elevation map was divided into five categories due to the study area altitude, whereby 1350-2100, 1100-1350, 920-1100, 750-920 and 500-750 above mean sea level (a.m.s.l.) were given ratings of 2, 4, 6, 8 and 10, respectively. In this research, appropriate elevations for the selection of landfill sites were given grading values of 8 and 10 ( Figure 3C).
The layer "urban area" criterion was classified into five categories. Buffer zones of more than 20 km were given the highest score of 10. Buffer zones of 15-20 km, 10-15 km, 5-10 km and less than 5 km were given the scores of 7, 5, 3 and 0, respectively ( Figure 3D). Cover soil in the Sulaimaniyah government is classified into four types of soil. These soil groups are (S33), consisting of Brown Soil, Medium and Shallow Phase over "Bakhtyary Gravel," (S37), consisting of Chestnut Soil, Deep phase, (S38), consisting of Rough Broken and Stony Land, and (S39), defined as Rough Mountainous Land. These soil types were ratings of 4, 6, 8 and 10 respectively ( Figure 3E).
Four categories have been identified for Land Slope, less than 5 • , 5-10 • , 10-15 • and 15-50 • , and these were given ratings of 10, 8, 6 and 4 respectively. According to the literature, less than 5 • is considered the best slope for landfill locations to avoid the transportation of waste leachate to adjacent areas ( Figure 3F).
In this study, buffer zones should be more than 30 m from a landfill site to power lines on both sides, and they were given a score of 10, with distances of less than 30 m being given a grading of 0. The reason why these values are chosen to prevent the high voltage energy resulting from these lines ( Figure 4A). The "archaeological sites" criterion was categorized into three sub-criteria, a buffer zone of less than 1 km has been classified as zero. Buffer zones of 1-3 km and those of more than 3 km were rated as 5 and 10 respectively ( Figure 4B).
Buffer zones of less than 5 km and 1 km were provided zero ratings in the "oil and gas" and "village" layers. Buffer zones higher than these values were given 10 grading values. (Figure 4C,D).
Sub-criteria for "Geological formations" consisted of seven groups. G1, D1, E2, F1, C4, B4 and A3 were given ratings of 0, 2, 3, 4, 6, 8 and 10, respectively ( Figure 4E). The suitability index in these groups was classified according to the lithology and permeability of the sediments due to grain size distribution. The lithology of the geological formations are described as below [42]: The 2017 Landsat 8 (OLI) satellite image was downloaded for land use and land cover classification using remote sensing software ENVI 5.3 for image processing. Geometric correction was used through pre-processing calibration. The methods for image processing include statistical algorithms that change the visual appearance of the image, and fast line-of-sight atmospheric analysis of hypercubes (FLAASH) converted the digital images to surface reflection to remove atmospheric effects and generate a reflection of the surface image. Using the pre-processed classified images from the training data (region of interest, ROI) and considering the maximum likelihood of estimating probabilistic model parameters in a supervised classification technique based on Bayesian theory [59], the layer of "land use" was classified into ten sub-criteria: Airport, Build up, Water body, Factories and Industrial area, Agriculture and fertile land, non-fertile and agriculture land. These were given a rate of zero due to unsuitability, while Forest, Pasture, Rock and Unused land were given a rate of 10, 8, 6 and 4, respectively ( Figure 4F).
To generate the groundwater layer map, the Kriging method in geostatistical analysis tools was used in GIS to obtain the best interpolation between groundwater depths in 235 drilled wells throughout the study area. The groundwater level data measurements were carried out from 2000 to 2015 and were obtained from the (Groundwater Directorate in Sulaimaniyah). Generally, landfill site selection is influenced by groundwater depths due to the extent of groundwater contamination through leaching. In this study, groundwater depths were classified into five intervals, from 19-

Analytical Hierarchy Process (AHP)
The Analytic Hierarchy Process (AHP) is a multi-criteria decision-making strategy where all criteria are assigned to different levels using pairwise comparison in the matrix to derive the relative weights for criteria [58]. It provides the judgment matrix by comparing the degree of significance of the relative element [58,60], and each index's weight is relative to the general objective, thereby simplifying the decision-making process. The process judgment of pairwise comparisons is done through the priority importance of intensity using a numerical scale of 9 degrees, according to the sequence values from 1 (Equal importance) to 9 (Extreme importance) as presented in Table 4. AHP output evaluates the general preference for each of the options, which ultimately allows decision makers to choose the best strategy [53,58]. Thus, every alternative can be evaluated in terms of the decision criteria and each criterion can be estimated by its weight in order to implement a pair comparison in the matrix (relative scale of

Analytical Hierarchy Process (AHP)
The Analytic Hierarchy Process (AHP) is a multi-criteria decision-making strategy where all criteria are assigned to different levels using pairwise comparison in the matrix to derive the relative weights for criteria [58]. It provides the judgment matrix by comparing the degree of significance of the relative element [58,60], and each index's weight is relative to the general objective, thereby simplifying the decision-making process. The process judgment of pairwise comparisons is done through the priority importance of intensity using a numerical scale of 9 degrees, according to the sequence values from 1 (Equal importance) to 9 (Extreme importance) as presented in Table 4. AHP output evaluates the general preference for each of the options, which ultimately allows decision makers to choose the best strategy [53,58]. Thus, every alternative can be evaluated in terms of the decision criteria and each criterion can be estimated by its weight in order to implement a pair comparison in the matrix (relative scale of importance). The values of (a ij ) when (i i = 1, 2, 3, . . . , m) and (j = 1, 2, 3, . . . , n) are used to indicate the performance values in the rows and columns of the matrix. The upper diagonal triangle of the matrix is filled with the values of comparison criteria, while the lower triangle of the matrix represents the reciprocal values of the upper diagonal, using Equation (1) [10,26,62].
where, (a ij ) is the element of the row (i) and column (j) of the matrix. The typical comparison matrix for the relative importance of the criteria can be represented as follows in a decision matrix: The eigenvectors for each row are calculated using geometric principles, multiplying the value for each criterion in each column in the same row of the original pair-wise comparison matrix, and then applying this to each row. This is represented as follows [58]: where, Eg i = eigenvalue for the row (i); n = number of elements in row (i) The priority vector or AHP weight is determined by normalizing the eigenvalue to 1 (divided by their sum) [58], as follows: The maximum eigenvalue or maximum lambda (λ max ) is obtained from the summation of products between each component of the priority vector and the sum of columns of the reciprocal matrix, as shown in the following formula [58]: where, a ij is the criteria in each column in the matrix, and W i is the value of weight for each criterion which corresponds to the priority vector in the decision matrix.
The consistency index represents the equivalent of the mean deviation of each comparative element and the standard deviation of the evaluation error from the true ones [63], which will often turns out to be larger than the value describing a fully consistent matrix in order to provide a measure of severity of this deviation [58,64].
where, CI is the consistency index, n is the size or order of the matrix. The consistency ratio (CR) is obtained by dividing the consistency index value (CI) by the random index value, here (RI = 1.56), and n is the size of the matrix (n = 13). Table 5 shows the mean random index value RI for a matrix with different sizes [58].
In this study, RI 13 = 1.56 and CR = 0.027. CR should be less than 0.1 because, if CR is less than 0.1, the ratio indicates a reasonable of consistency level in the pairwise comparison. The importance and weight of each parameter were compared with each parameter in the weight assessment procedure. This was done by taking account the views of the expert that have worked in this field. Each parameter specifies a value for the weight that it merits by assuming the method of straightforward additive weighting. These weights are then used to organize the comparison matrix to revise the correct weight for each parameter. The pairwise comparison matrix for the 13 criteria is presented in (Table 6).  In addition, multi-criteria decision-making (MCDM) methods were used to obtain the comparative weights for each criterion using different methods including: Analytical Hierarchy Process (AHP), Simple Additive Weighting (SAW), Ratio Scale Weighting (RSW) and Straight Rank Sum (SRS).
To obtain the suitability index value for the potential areas, the WLC (Weighted Linear Combination) method was used based on the following formula where, A i : is the suitability index, W j : is the relative importance weight of the criterion, C ij : is the grading value (i) under criterion (j), and n: is the total number of criteria [28,29]. The WLC method was applied on all criteria using the special analysis tools "map algebra" in GIS to all criteria using the special analysis tools "map algebra" in GIS to assess the suitability index map. This is achieved by summing up the products by multiplying the sub-criteria rating values for each criterion (based on expert opinion in this sector) by the corresponding relative importance weight.

Simple Additive Weighting (SAW) Method
Simple additive weighting (SAW) is a basic ranking method and is specified as a weighted linear combination or scoring method using a multi-attribute decision approach [56]. The approach is based on weighted averages, with the weights of comparative significance directly allocated by the decision-maker to assess the weight for each criterion and to select the significance and weight of each criterion relative to the other criterion. The normalization of the relative weight of importance is done by dividing the weight of the criterion by their summation. The simple additive weighting method evaluates each alternative using the following formula [44,66].
where, Wi is the normalized weight of each criterion by its sum; A i is the weight of each criterion in the row (i); A j is the weight of each criterion in the column (j); n is the number of criteria.

Straight Rank Sum (SRS) Method
The straight rank sum method (SRS) is considered to be a kind of ranking method. It is a straightforward technique used to determine the weights of the criteria by arranging the criteria according to comparative significance in decreasing order from the most to the least important depending on the preferences of the decision-makers using the formula (n − r i + 1). Then, the weights of the criteria are normalized by dividing each weight of criterion by its sum Σ (n − r g + 1) according to Equation (9) [18,67].
where, W i is the relative importance of normalized weight for (i) criterion; n is the number of criteria in the research under consideration (g = 1, 2, . . . , n), and r i is the ranking position for the criterion.

Ratio Scale Weighting (RSW) Method
In the Ratio Scale Weighting (RSW) technique, the decision process is based on the allocation of an appropriate ratio score value for each criterion, where the value of 100 is provided to the most important criterion as the basis for the other criteria values. The values below 100 are assigned proportionally to criteria which are lower in order based on the importance of each criterion compared to the other criteria [44]. In order to estimate the current criteria weights (NW i ) using the Ratio Scale Weighting (RSW) technique, the proportional weight value of each criterion is divided by the lowest in importance proportional weight value of the criterion. Then, the normalized weights of the criteria, RSW method, are estimated as follows [67]: where, W i is the normalized weight of each criterion divided by the sum of new weight of each criterion; NW i is the new weight of each criterion of the area (i) under criterion (j); n is the number of criteria.

Results and Discussion
From constructing a matrix of comparison, the maximum eigenvalue or maximum lambda (λ max ) in this study is equal to 13.51 and consistency index (CI) equals 0.04, criterion weight results defined by all MCDM methods (including AHP, SAW, RSW and SRS methods) are show in Table 7. The suitability index map for landfill sites in the Sulaimaniyah Governorate adapted using the WLC methods for AHP, SAW, SRS and RSW, where criteria branched with higher sub-criteria take higher weighting in the AHP method than the others, as presented in Figure 6. Therefore, unsuitable areas take a higher portion of the area in the AHP method compared to other methods, as presented in Figure 7.
The suitability index value in each method is categorized into four classifications of suitability areas, these categories are unsuitable areas, moderately suitable areas, suitable areas and most suitable areas, each method represent a different value in area, (Figure 7A-D) Each category's area is calculated over the total study area which resulted in the final index of suitability for landfill locations as assessed in Figure 8.
The suitability index map shows that each method slightly differs from the other methods in terms of area, but all methods were overlaid to determine the intersecting area within the map index, as seen in Figure 9. The suitability index map shows that each method slightly differs from the other methods in terms of area, but all methods were overlaid to determine the intersecting area within the map index, as seen in Figure 9.
Total matching of suitable areas for all method was estimated as being 80%, while non-matching area represents 20% of the area. Thus, seven candidate sites were suggested for the landfill in the Sulaimaniyah Governorate in the "most suitable index" category according to the AHP, SAW, RSW and SRS methods, with areas of 2.74, 4.48, 2.80, 3.67, 2.67, 3.18 and 4.36 km 2 respectively as shown in Figure 10. Given the expectation of huge annual waste generation of 473,000 tons over a small open dumping disposal site of about 0.45 km 2 in area [58], this value needs to be more than 2.5 km 2 for future planning  Total matching of suitable areas for all method was estimated as being 80%, while non-matching area represents 20% of the area. Thus, seven candidate sites were suggested for the landfill in the Sulaimaniyah Governorate in the "most suitable index" category according to the AHP, SAW, RSW and SRS methods, with areas of 2.74, 4.48, 2.80, 3.67, 2.67, 3.18 and 4.36 km 2 respectively as shown in Figure 10. Given the expectation of huge annual waste generation of 473,000 tons over a small open dumping disposal site of about 0.45 km 2 in area [58], this value needs to be more than 2.5 km 2 for future planning.

Conclusions
The existing waste disposal sites in the Sulaimaniyah Governorate do not fulfill science and environmental requirements. The goal of this study was to pursue proper selection processes while avoiding environmental problems, and to suggest suitable sites for landfill using GIS and Multicriteria decision-making (MCDM) techniques to facilitate shifting the present dumpsite to an alternative location. To decide on the most adequate sites for landfills in the Sulaimaniyah Governorate, thirteen layers of criteria were incorporated into GIS to make a decision analysis process using a pairwise comparison through constructing a matrix of comparison. These layers, according to their importance were urban area, villages, rivers, groundwater, slope, elevation, soil, geology, road, oil and gas, land use, archaeology and power lines.
The AHP, SAW, RSW and SRS methods were used in MCDM to derive the weightings of the criteria, and the WLC method was applied on all methods to produce the output map, including a landfill suitability index for the Sulaimaniyah Governorate. The final map was produced from the intersection of all methods' maps, with total matching for suitable areas for all methods estimated at 80%, while non-matching area presents 20% of the area. Seven candidate locations were suggested on the final map for landfill in the Sulaimaniyah Governorate within the category of "most suitable" area.

Conclusions
The existing waste disposal sites in the Sulaimaniyah Governorate do not fulfill science and environmental requirements. The goal of this study was to pursue proper selection processes while avoiding environmental problems, and to suggest suitable sites for landfill using GIS and Multi-criteria decision-making (MCDM) techniques to facilitate shifting the present dumpsite to an alternative location. To decide on the most adequate sites for landfills in the Sulaimaniyah Governorate, thirteen layers of criteria were incorporated into GIS to make a decision analysis process using a pairwise comparison through constructing a matrix of comparison. These layers, according to their importance were urban area, villages, rivers, groundwater, slope, elevation, soil, geology, road, oil and gas, land use, archaeology and power lines.
The AHP, SAW, RSW and SRS methods were used in MCDM to derive the weightings of the criteria, and the WLC method was applied on all methods to produce the output map, including a landfill suitability index for the Sulaimaniyah Governorate. The final map was produced from the intersection of all methods' maps, with total matching for suitable areas for all methods estimated at 80%, while non-matching area presents 20% of the area. Seven candidate locations were suggested on the final map for landfill in the Sulaimaniyah Governorate within the category of "most suitable" area.