Landfill Site Selection Using GIS and Multi-Criteria Decision-Making AHP and SAW Methods: A Case Study in Sulaimaniyah Governorate, Iraq ()
1. Introduction
Solid waste management is considered to be a significant issue in developing countries. Population growth, economic recovery and industrial growth are all reasons to increase the generation of solid waste in developing countries [1].
Despite the use of many efficient processes such as reuse and recycling, appropriate landfill disposal is still the most prevalent method of minimizing adverse effects on the environment and waste management [2]. The common issue facing all developing countries is the disposal of solid waste and the availability of land, considering its significant adverse impacts on the environment [3] [4].
Currently, there is no landfill site in the study area that fulfils scientific and environmental requirements to resolve the waste dump site issue, the process of selecting a landfill site is considered complicated task, the combination of multi-criteria decision-making (MCDM) approaches and the Geographic Information System (GIS) generates a powerful spatial decision support that provides the opportunity to effectively create land suitability maps for waste disposal sites [5], multi-criteria decision-making (MCDM) and GIS have been used commonly in different fields and implementations, including the integrated eco-environment Assessment of Soil [6], land evaluation for peri-urban agriculture [7], possibility of groundwater pollution [8]. GIS becomes an important tool for smart decisions on landfill site selection [9].
In order to achieve the research objective, thirteen important criteria that have influenced the environment and waste management have been considered as a data set for the decision model, two techniques of multi-criteria decision-making (MCDM) have been used in this research, which are analytic hierarchy process (AHP) and simple additive weighting (SAW) methods to assess the weights of prospective variables for selecting landfill sites [10] [11] [12]. Geographic Information System (GIS) with AHP and SAW methods are the most common techniques in MCDM with a high capacity to manage complicated problems with large data during the decision-making process [13] [14] [15].
This study aimed to compare both methods and establish an appraisal blueprint to find the best candidate landfill sites that realize the environmental and scientific criteria.
2. Study Area
Sulaimaniyah is among the major cities in the Kurdistan region. The city is situated northwest of Iraq between latitude 35˚45'0''N, 36˚0'0''N and longitude 44˚45'0''E, 45˚45'0''E approximately 370 km north east of Baghdad, Iraq’s capital. The city is bounded in the north-east and south-west by the Mountains and is situated in a low-lying land covering an area of approximately 2400 km2. The study area is characterized by a separate Mediterranean-type continental interior climate with average annual precipitation ranging from (500 to 700 mm). The Sulaimaniyah governorate had a population of approximately 856,990 in 2017 [16]. Figure 1 shows the administrative boundary unit of the Sulaimaniyah governorate [17]. All types of waste dumped without treatment in an open area overlooking the Tanjaro River.
3. Materials and Methods
3.1. Dataset Criteria Map
In order to proceed suitable landfill site map, thirteen criteria as layer maps were prepared using GIS spatial analysis tools over the study area, these layers were urban area, villages, rivers, groundwater depth, slope, elevation, soil types, geological formations, roads, oil and gas field, land use classification, archaeological site and power lines, in accordance with environmental standards, natural and artificial factors for landfill site selection.
The source of the data obtained from official government authorities and International organization data base, official government authorities data includes
![]()
Figure 1. Location map of the study area.
urban area, archaeological sites, power lines, oil and gas fields, villages, soil and geological formations. The river, road, and elevation data downloaded from the United States Geological Survey USGS Earth Explorer, Spatial analysis tools used in GIS to convert the slope map from a digital elevation model.
Water level depth data were obtained from the Sulaymaniyah groundwater authority and GIS was applied to water level using the “Kriging” method in special analysis tools to create a groundwater table map of the study area. The land use classification was prepared using satellite data and processed by remote sensing software (ENVI 5.4).
3.2. Criteria Restriction
Determining the allowable distance from landfill sites requires consideration of government regulations, prospective environmental risks, public health and economical evaluation for each criterion [18] [19].
Specific geographical features established using buffer zones by spatial analysis of GIS software around each criterion, buffer zones were created based on previous literature studies to determine the distance from each feature to the specified criteria. A buffer zone is an area that can be divided by grade to reduce or eliminate the impact of land use activities on vulnerable regions or natural features, restricted criteria and suggested buffer values for the study area as shown in Table 1.
3.3. Sub-Criteria Rating Values
Each criterion was classified into sub-criteria and assigned a suitability rating
![]()
Table 1. Restricted criteria and suggested buffer values for the study area.
value from zero to ten [36] [37]. The criteria rating and importance of its priority were specified based on restrictions on category priorities for the field of study and on the basis of literature and research experts in the field of selecting solid waste sites.
The ranking value for each criterion and sub-criteria was determined following several steps, including in a series (Buffer, Clip, Extract, Overlay, Proximity, Convert, Reclassify and Map Algebra) using GIS spatial analysis tools. Sub-criteria buffer zone and rating values for the input layer are shown in Table 2.
![]()
Table 2. Layers buffer zone with sub-criteria ratings.
(a.m.s.l.): Above Mean Sea Level.
In this revise, sub criteria rating value of 0 is corresponding to the nearest restricted area from the landfill, and a rating value of 10 was provided best area, for example the sub-criteria “Geological Formations” consisted of seven groups G1, D1, E2, F1, C4, B4 and A3 respectively (Figure 3(H)) were given. ratings of 0, 2, 3, 4, 6, 8 and 10 respectively, The suitability index for these groups was graded according to the lithology and permeability of the sediments due to the distribution of grain size [38]. Buffer zones and suitability index maps as shown in Figures 2-4.
3.4. Multi-Criteria Decision-Making Methods
Pairwise comparison implemented in the matrix for all criteria through the priority of the importance intensity of one activity over another using a numerical scale of 9 points [39].
![]()
![]()
Figure 2. Buffer zones and suitability index maps: (A) Urban area; (B) Villages; (C) Rivers; (D) Groundwater depth; (E) Slope; (F) Elevation.
![]()
![]()
Figure 3. Buffer zones and suitability index maps: (G) Soil types; (H) Geological formations; (I) Roads; (J) Oil and gas field; (K) Land use classification; (L) Archaeological site.
![]()
Figure 4. Buffer zone and suitability index map: (A) Power lines.
The upper triangular matrix is filled with the comparative criteria values and the lower triangular matrix is completed with the upper reciprocal values [10] [11]. The eigenvalue is calculated by multiplying the value for each criterion in each column in the same row in the matrix of the pairwise comparison. The priority vector (Pri) is determined by normalizing the eigenvalue to 1 [39] as follows:
(1)
where, Egi = eigenvalue for the row (i)
; n = number of elements in matrix row (i).
The consistency index calculated according to [39].
The maximum lambda (λmax) is obtained from the summation of products between each element of priority vector and the sum of columns of the reciprocal matrix as shown in the following formula:
(2)
where, Wj is the value of weight for each criterion which corresponds to the priority vector in the decision matrix and aij is the criteria in each column in the matrix.
(3)
where, CI consistency index and n is size or order of the matrix, (λmax) which is equivalent to the priority vector in the matrix of decision [39].
The consistency ratio (CR) depends on the size of the matrix (n = 13) thus, random index value (RI = 1.56) [39]. Table 3 shows the Random inconsistency value RI in different sizes for a matrix [24] [40].
(4)
Simple additive weighting (SAW) is a ranking method and defined as a weighted linear combination or scoring method [39].
(5)
where, Wi is the normalized weight of each criterion which was, Ai is the weight of each criterion of area (i) under criterion (j); n is criteria number.
4. Results and Discussion
The matrix of pairwise comparisons with SAW and AHP weighs as presented in Table 4. The maximum lambda (λmax) =13.51, CI = 0.04 and CR = 0.027, If CR is less than 0.1 the ratio indicates a reasonable consistency level in the pairwise comparison [41]. The final map shows the suitability index for landfill sites in Sulaimaniyah Governorate which was divided into four categories of suitable areas, including: unsuitable, moderately suitable, suitable and most suitable areas [12], suitability index with areas for all categories of the SAW and AHP methods as shown in Figure 5.
![]()
Table 3. Random inconsistency indices for different values of (n) [24] [40].
![]()
Table 4. Pairwise comparison matrix with (AHP) and (SAW) methods.
Pairwise comparison matrix (A): Urban area; (B): Villages; (C): Rivers; (D): Groundwater depth; (E): Slope; (F): Elevation; (G): Soil types; (H): Geological formations; (I): Roads; (J): Oil and gas field; (K): Land use classification (L): Archaeological site; (M): Power lines. NW = Normalized weight.
![]()
Figure 5. Suitability index area for landfill site using SAW and AHP methods.
5. Conclusion
This research used the MCDM techniques with the GIS method to evaluate the suitable selection of landfill sites in the study region. The result shows the index values that have been categorized into 4 areas with calculated area using the pixel calculation in GIS. The results indicate that the most suitable area covered the area of 16.37% and 24.35% or 392.92 and 591.71 km2 respectively in SAW and AHP methods, the compatibility of the most suitable area in both methods is 91.71 percent, while the compatibility of all zone areas in both methods is 99.8, 94.7 and 96.85 percent, respectively, for unsuitable, moderately suitable and suitable respectively.