Muluneh Beyene Industry Site Selection
Muluneh Beyene Industry Site Selection
JUNE, 2019
ADDIS ABABA, ETHIOPIA
GIS AND REMOTE SENSING BASED ANALYSIS ON INDUSTRIAL SITE
SELECTION IN HAWASSA TOWN, ETHIOPIA
June, 2019
Addis Ababa, Ethiopia
STATEMENT OF APPROVAL
Signed by
____________ ___________
Thesis Advisor Signature Date
____________ ___________
____________ ___________
_________________________________________
Chair of Department or Graduate Program Coordinator
i|Page
STATEMENT OF DECLARATION
This thesis is my original work and has not been presented for academic purpose or degree
award in the Addis Ababa University or any other university. Therefore, hereby, I declare that
this thesis Paper contains no section copied in whole or in part from any sources unless explicitly
quoted.
Declared by:
Signature__________________ Date:_________________
ii | P a g e
ACKNOWLEDGMENT
First and foremost, I would thank my Almighty God, for the source of all fulfillments and his
endless support to complete this thesis paper in my life.
I am deeply indebted to my advisor, Teferi Mekonnen(Ph.D) for his time and guidance in doing
the research work.
I would like to thank the following organizations: Goffa Zone, Demba Goffa Woreda Finance
and Economic Development Sector for their cooperation, Hawassa town Trade and Industry
Department, SNNPRS Trade and Industry Bureau and Urban Development and Housing Bureau
for their kindly provision of the required data for this study. I would express a great pleasure
from the bottom of my heart to Addis Ababa University, Collage of Social Sciences for their
financial support to execute this thesis.
I am interested to provide special thanks and heartfelt gratitude for my mother Batale Beyene
and my uncle Belachew Beyene for their moral initiation and financial support in strengthening
my academic level.
I want to thank my classmates who showed me love and support from the first day of class up to
final thesis submission.
I am so grateful and would like to extend my heartfelt thanks to my remaining parents, brothers
and sisters. May the Almighty God Bless Them!
iii | P a g e
TABLE OF CONTENT
STATEMENT OF APPROVAL...................................................................................................... i
ACKNOWLEDGMENT................................................................................................................ iii
1. INTRODUCTION ................................................................................................................... 1
iv | P a g e
2.1.2. Site selection process ................................................................................................ 7
2.2.3. GIS and Multi-Criteria Analysis for Industrial Site Selection ...................................... 9
v|Page
3.1. Description of the Study Area ........................................................................................ 25
vi | P a g e
CHAPTER FOUR ......................................................................................................................... 45
4.5. Distance from the Residential Area and Industrial Location ......................................... 57
4.10. Integration of Criteria Maps and Preparation of Final Suitability Map ..................... 66
vii | P a g e
CHAPTER FIVE .......................................................................................................................... 70
REFERENCES
Appendix 1:
Appendix_2:
Appendix _3:
Appendix_4:
Appendix 5:
Appendix_ 6:
viii | P a g e
ACRONYMS AND ABBREVIATIONS
ix | P a g e
LIST OF FIGURES
Figure 1: Site selection process-----------------------------------------------------------------------------8
Figure 2: Conceptual framework of the study-----------------------------------------------------------24
Figure 3: Location map of the study area---------------------------------------------------------------25
Figure 4: Temperature distribution of the study area--------------------------------------------------27
Figure 5: Rainfall distribution of the study area--------------------------------------------------------27
Figure 6: Geology fault map of the study area ----------------------------------------------------------29
Figure 7: Soil types of the study area-------------------------------------------------------------------30
Figure 8: Road network of the study area----------------------------------------------------------------31
Figure 9: Major industrial distribution in the study area-----------------------------------------------32
Figure 10: Slope map of the study area-------------------------------------------------------------------35
Figure 11: Lake map of the study area--------------------------------------------------------------------36
Figure 12: River map of the study area ------------------------------------------------------------------37
Figure 13: Swamp map of the study area ---------------------------------------------------------------37
Figure 14: Land use land cover map of the study area ------------------------------------------------38
Figure 15: Road network of the study area --------------------------------------------------------------39
Figure 16: Residential area map of the study area-------------------------------------------------------39
Figure 17: Figure 17: Map of distance from geological fault of the study area---------------------40
Figure 18: soil type map of the study area----------------------------------------------------------------41
Figure 19: Over all workflow of the study---------------------------------------------------------------44
Figure 20: Suitability map of slope -----------------------------------------------------------------------46
Figure 21: Lake suitability map of the study area ------------------------------------------------------47
Figure 22: River suitability map of the study area ------------------------------------------------------49
Figure 23: Swamp suitability map of the study area ---------------------------------------------------50
Figure 24: Land use land cover suitability map of the study area ------------------------------------55
Figure 25: Road suitability map of the study area -----------------------------------------------------56
Figure 26: Residential area suitability map -------------------------------------------------------------58
Figure 27: Suitability map of geological fault --------------------------------------------------------59
Figure 28: Soil type suitability map of the study area--------------------------------------------------60
Figure 29: Weighted overlay suitability map of the study area ---------------------------------------68
x|Page
LIST OF TABLES
Table 1: Industrial Types in the Study Area-------------------------------------------------------------32
Table 2: Types of Data Used for this Study--------------------------------------------------------------33
Table 3: Software and Material Applied in this Study--------------------------------------------------33
Table 4: Suitability requirements rating for the selected factors in the study area------------------34
Table 5: The preference scale for pair wise comparison in AHP--------------------------------------43
Table 6: Reclassified Slope Value and Suitability Level-----------------------------------------------46
Table 7: Reclassified Proximity to Lake Suitability Level---------------------------------------------47
Table 8: Reclassified Proximity to River and Suitability Level---------------------------------------48
Table 9: Reclassified Proximity to Swamp and Suitability Level-------------------------------------49
Table 10: Accuracy Assessment Matrix-----------------------------------------------------------------51
Table 11: Error Matrix of the Estimated Area Proportion----------------------------------------------52
Table 12: Estimated Land Use Land Cover and confidence interval---------------------------------53
Table 13: Area estimates of each land use class based on reference data (Km2) -------------------54
Table 14: Land Use Land Cover Classes and Suitability Level---------------------------------------54
Table 15: Road Distance and Suitability Level----------------------------------------------------------56
Table 16: Residential Area Suitability--------------------------------------------------------------------57
Table 17: Geological Fault and its Suitability-----------------------------------------------------------59
Table 18: Soil Type and its Suitability--------------------------------------------------------------------60
Table 19: pair-wise comparison, 9-point weighting scale----------------------------------------------61
Table 20: The Criteria used in industrial site selection-------------------------------------------------63
Table 21: Normalized Criteria Selection for Industrial Site Selection--------------------------------64
Table 22: Classified Factor Values and Weights By Pair-Wise Comparison Method -------------65
Table 23: Statistical Analysis for Overall Weights for the Industrial Sites Selection -------------67
xi | P a g e
ABSTRACT
The significance of GIS and Remote Sensing techniques and technologies currently in every
aspect of human beings day to day activities initiated the researcher to conduct this study.
Industry is important to enhance the economic level of a given country and GIS and remote
sensing technologies would provide a variety of options during site selection process. The
overall objective of the study is to identify suitable industrial site location by using Geographic
Information System and Remote Sensing based analysis in Southern Nations, Nationalities and
Peoples Regional State, Hawassa town. In order to conduct this study, multi criteria evaluation
system to identify factors influence during site selection by pair-wise comparison matrix was
employed. Nine determining factors/parameters identified in this study to identify suitable site
for the industrial establishment. These factors were slope, land use land cover, soil types,
distance from geological faults, proximity to main road, distance from residential area, distance
from water bodies (lake, river and swamp), The findings revealed that land use land cover
factors and distance from geological faults highly determines the selection for industrial sites;
and elevation and slope determines minimum compared with other factors. Majority of the study
area (47.5%) is suitable to establish industries and a small portion (0.2%) are less suitable.
Conducting environmental impact assessment prior to establishing industries in the study area
was not fully covered by the concerned government organization. Application of GIS and Remote
Sensing technologies in the study area is not fully employed. In order to select suitable sites for
the industry, the town administrators should focus on those important parameters by using
Geographic Information System and Remote Sensing technologies.
xii | P a g e
CHAPTER ONE
1. INTRODUCTION
1.1. Background of the Study
Principally in developing countries of big cities, a number of directors took the initial step of
creating special regions for industry, essentially to isolate them from heavily populated or
affected areas. Government agencies, in charge of regulating the industry site selection
commenced considering the inclusion of environmental criteria in the selection process as a
measure to lessen potential environmental impacts to local communities. In the past, site
selection was based almost purely on commercial and technical criteria. Today, a higher degree
of refinement is expected. Selection criteria must also meet a number of social and
environmental elements, which are enforced by law and government regulations (Aleksandar, et
al 2013). GIS is an efficacious tool for the management and analysis of data required for any
land development activity (Manish and Vivekananda, 2013). GIS system has the ability of
display needed in the context of decision making. A set of tools has been used to manage the
proper site for a city improvement facility, including GIS, and Multi-Criteria Decision Making
(MCDM) techniques (Zahra et al, 2014).
Industrial site selection is critical point in the process of starting, expanding or changing the
location of industrial systems of all kinds. One of the main objectives in industrial site selection
is finding the most appropriate site with desired conditions defined by the selection criteria. In a
site selection process, the analyst strives to determine the optimum location that would satisfy
the selection criteria. The selection process attempts to optimize a number of objectives desired
for a specific facility. Such optimization often involves numerous decision factors, which are
frequently contradicting, and the process often involves a number of possible sites each has
advantages and limitations. Decision making is based on numerous data concerning the problem
of selection appropriate site. Decisions about industrial location typically involve the evaluation
of multiple criteria according to several, often conflicting, objectives. While many decisions we
make are prompted by a single objective (Eldin and Sui, 2003).
The principles of liberalization, globalization and privatization and the relative changes at the
global economy have been very important for the industrial zones development. Most of the
countries use this tool for development of industrial field with minimizing its negative impacts
1|Page
on people and the environment. This has been very important as a better environ-mental
management tool which can be used for sustainable development of the industrial sector. The
best location for the industrial zone is the key factor of the accomplishment of its main objectives
which are the requirement of all facilities and the minimum of social and environmental negative
impacts (Thanuja et al, 2015).
Site selection of industrial zones in a given area has become a critical issue and a sensitive
decision making process that may create a range of socioeconomic and environmental problems
over time. Hence, several site selection criteria and appropriate methods for establishing
industries have to be concerned by the decision makers and authorities before locating industrial
zones in particular regions. Though, at present the enormous data volume and complex criteria
regarding this field are available, the suitable site selection process is still problematic (Fernando
et al, 2015).
Site selection is one of the primary decisions in the start-up process and town planning,
development or relocation of industries of all kinds. Building of a new industrial system is
related with long term investment, and in this sense planning the site are a significant point on
the road to progress or failure of the industrial system and that effect on all other services in the
city. One of the main objectives in industrial site selection is finding the most suitable site with
required conditions. A large number of researchers depend on GIS because of its availability and
its wide uses in site selection process. GIS is used in combination with other systems and
methods such as the method for multi-criteria decision making (MCDM). Multi-criteria decision
analysis (MCDA) techniques can be used in such conditions to classify and rank options for
subsequent complete evaluation, or to specify acceptable from unacceptable potentiality for
many sites (Mahmud, 2015).
This research thesis focuses on how to apply GIS based analysis to select industrial sites in
Southern Ethiopia, Hawassa town. Site selection with minimum impacts is very difficult to
achieve considering criteria one by one, because of the big list of criteria including very
important spatial factors. But still in Hawassa town this process is conducted as politicians
need. In this study, the researcher considered the criteria that need to be undertaken before
establishing industrial zones in Hawassa town. That is the one of major significance of this
study. The GIS based analysis was developed in this research that can be used for every
where to get very accurate steps for site selection process. It is very important especially for
the sustainable development of industrial field. Therefore, this paper assumes that, in some
2|Page
extent, such failures can be overcome by applying Geographic Information System (GIS)
and Multi-criteria Decision Making Techniques into the site selection process of industrial
zones. In this study, a GIS and Remote Sensing based analysis was applied to screen most
suitable locations for establishing Industrial Zones in Hawassa town.
The goal in a site selection exercise is to find the best location with desired conditions that
satisfy predetermined selection criteria. Site selection typically involves two main phases: a) site
screening (identification of a small number of candidate sites from a broad geographic area and a
range of selection factors) and b) site evaluation (in-depth examination of each candidate site to
the most suitable (Davis, 2001).
According to Davies (2001), the selection process attempts to optimize a number of objectives in
determining the suitability of a particular site for a defined facility. Such optimization often
involves a multitude of factors, sometime contradicting. A number of tools were used to
determine the proper site for capital improvement facilities. These tools include Expert Systems,
geographic information systems (GIS), and Multi-criteria decision-making (MCDM) techniques.
These tools have played an important role in solving site selection problems.
According to the Hawassa Town Sanitation, Beautification and Parks Development Department
report (2012), there are many of potential sites that could be chosen for the facility. Each of these
sites usually satisfies the selection criteria at varying levels. It is fairly common to find
3|Page
contradiction between the multiple objectives for a sitting problem. Land use designation
problem (residential, industrial, commercial) sites were not as such differentiated. i.e. residential,
commercial and industrial sites located near to near. The town administrators are simply
selecting those areas that are inhabited/bare lands/ for industrial sites without identifying other
preconditions like social integration of communities, proximity to the town, safety of the
community since industrial sites has their own negative impacts on society.
For instance, the objective of keeping minimum capital investment may contradict with the
objective of keeping a long-term safe environment. Many objectives lack means for quantitative
measurement. Examples of those are the aesthetic deterioration of the view of a natural mountain
scene as a result of the installation of transmission towers/lines, the social disruption felt by a
community as a result of the expected rapid influx of workers during construction, and similar
issues.
Different studies undertaken in relation with industrial site selection in Hawassa town mainly
focused on identifying sites without considering parameters and factors based on GIS and remote
sensing techniques. For example, Mamo Mihretu and Gabriela Llobet (2017) conducted a case
Study of PVH’s Commitment in Ethiopia’s Hawassa Industrial Park. In this study, on page 17,
they tried to describe how industrial sites screened/selected. But they ignored about the way how
industrial sites screened or selected by the help of GIS and remote sensing technologies
Another study undertaken by Dejene et al (2018), on industrial park development in Ethiopia mainly
discussed about integrated agro-industrial parks are under development in the selected regions.
The study focused on the importance of industry in relation with job creation, strong linkages to
the agricultural sector, high export potential and capacities to attract private sector investment.
The study did not stated about site selection process since effective industrial development has
direct relation with the site it established.
Therefore, this study considered different factors/parameters that affect the selection of industrial
sites by applying GIS and remote sensing techniques.
4|Page
1.3. Objective of the Study
1.3.1. General objective of the study
The general objective of this study is to identify optimal industrial site location by using GIS and
remote sensing based analysis in SNNPRS, Hawassa town.
To identify different parameters that determines the selection of industrial site location in
the study area.
To identify suitable sites for industrial location for the future by applying GIS based
analysis in the study area.
To produce final suitability output map that shows suitable sites for industrial locations in
the town.
What are different parameters that determine the selection of industrial site location in the
study area?
Which sites are important to locate industries in the future in Hawassa town?
What seems final suitability output map that shows suitable sites for industrial locations
in the town?
5|Page
industrial site selection identified as the research focusing part and from different towns of
Ethiopia, Hawassa town taken as the study area. From this perspective, the researcher interested
to focus on GIS modeling in industrial site selection in SNNPRS, the case of Hawassa town.
6|Page
CHAPTER TWO
2. LITERATURE REVIEW
In this part, that information that was studied before by other authors was taken as stepping
stone to compare gaps. The researcher also tries to put the theoretical frame works that are
relevant for this study and conceptualizing those terms on his account.
7|Page
citizen can be one of the key factors for decision makers. All so risk management is an
indispensable analysis in site selection process. Managing the risks involved in selecting a new
industrial location is one of the most critical factors in determining the ultimate success or failure
of a business. To keep risks at a minimum, investors should first be familiar with the stages of
the site selection process and what are the key risks that need to be considered and managed
during each of these stages (Florida, 2002).
One of the most important and far reaching decisions faced by operations managers is deciding
where to locate new industrial facilities. This is a strategic decision involving irreversible
allocation of the firm’s capital, and often has a crucial impact on key measures of the firm’s
supply chain performance such as lead time, inventory, and responsiveness to demand
variability, flexibility, and quality (Bhatnagar, 2005).
Collection of information allows the generation of a potential industrial sites that can be grouped,
while the use of certain term criteria, through several iterations, gradually narrowing to a choice
(Fig. 1). In such way, the total number of available sites, the customer is aware of a certain
number of them. Of these, only a certain number of locations meet the selection criteria of the
decision maker, so that makes group of sites for consideration. By collecting information on
these sites, it remains just making a group of sites that are included in the shortlist. Out of this
group, based on the criteria used by the decision maker (investor) chooses one location (Kotler,
et al, 1993).
Phase I : Site Screening Phase II : Site Evaluation
Locations those Decision
Total Locations which the Locations that
are included in /selected
number of decision makers is are being
the selection locations
locations aware of consider
Figure 1: Site Selection Process (Source: Kotler, et al, 1993)
The process of site selection includes: Establishing a set of influential factors relevant to site
selection, predicting and evaluating the intensity and direction of their effects in time and given
conditions and evaluation of possible variants of solutions and selection of optimal variants.
8|Page
2.2.3. GIS and Multi-Criteria Analysis for Industrial Site Selection
9|Page
In general, the MCDA in GIS should be viewed as a process of conversion of data to information
that adds extra value to the original data. GIS techniques & procedures have an important role to
play in analyzing decision problems recognized as a decision support system for industrial site
selection especially in site screening phase. In industrial site screening phase role of GIS is to
geo-referenced and analyzes feasible alternatives that will be later consider in evaluation phase.
In evaluation phase roll of GIS is to produce criteria, constraints and suitability maps according
to the results from Multi-criteria decision analysis and value judgments of decision makers
(Malczewski, 2006).
GIS techniques & procedures have an important role to play in analyzing decision problems
recognized as a decision support system for industrial site selection especially in site screening
phase. In industrial site screening phase role of GIS is to geo-referenced and analyzes feasible
alternatives that will be later consider in evaluation phase. In evaluation phase roll of GIS is to
produce criteria, constraints and suitability maps according to the results from Multi-criteria
decision analysis and value judgments of decision makers.
10 | P a g e
alternatives is explicit. Multi-attribute decision problems require that choices be made among
alternatives described by their attributes. This implies that attribute objective relationships are
specified in such a form that attributes can be regarded as both objectives and decision variables.
Attributes are used as both decision variables and decision criteria (John Wiley and Sons, 1999).
The most significant factors that describe decision problems or affect the choice and
implementation of MCDA methods the most significant: Number of decision makers, Number of
objectives, Number of alternatives, existence of constraints and risk tolerance (Bal, 2006).
11 | P a g e
2.2.1. Problem definition
The first step in all decision-making processes is to recognize and define the relevant
characteristics of the decision problem. Depending on the business environment and the
company’s overall objectives, the problem definition specifies the essential inputs for the
decision making process. Managers would like to have everything in one location: low labor
costs and a quality workforce, simultaneously. These objectives are contradictory and very
difficult to achieve. This problem can only be solved through compromise. Therefore, it is
necessary to explain the problem and to find the space for compromise. In this step, we define
the main objectives that should not be in conflict and identify the decision makers (either an
individual or a group) for the location problem under consideration (Zelenovic, 2003).
12 | P a g e
for the subsequent multi-criteria analysis, which requires clearly defined location alternatives for
analysis. Each location alternative should be a feasible location for the industry and should meet
the basic construction requirements. For mining data and generating the location alternatives, we
developed a GIS by using Arc-GIS (McCoy, 2001)
13 | P a g e
2.2.9. Sensitivity analysis
The sensitivity analysis tests how much the sensitivity map and, therefore, our decisions are
sensitive to criteria variations. Thus, in the sensitivity analysis, we test the accuracy and the
robustness of the multi-criteria analysis. This allows decision makers to be confident in their
choices even though the values of the criteria may not be highly accurate.
In this study, the researcher employed weight normalized procedure how the processed data
ensures its influence in prioritizing parameters.
In the past, site selection was based almost purely on economical and technical criteria. But
currently, a higher degree of sophistication is expected. Selection criteria must also satisfy a
number of social and environmental requirements. Selection of industrial site means complex
multi-criteria analysis which includes a complex array of factors involving economic, social,
technical, environmental and political issues that may result in conflicting objectives (Reisi,
2011).
Nowadays, in the post-industrial society and knowledge-based society, people become the
most important resource. Managing the risks involved in selecting a new industrial location is
one of the most critical factors in determining the ultimate success or failure of a business.
To keep risks at a minimum, investors should first be familiar with the stages of the site
selection process and what are the key risks that need to be considered and managed during
each of the stages (Aleksandar et al, 2013).
According to USAID research conducted on Serbia, the following points are important
criteria to select industrial location: labor costs, geographic position, availability of a quality
workforce, transport infrastructure availability of raw materials, licensing and permitting
procedures for the land, telecommunication infrastructure, efficiency of local authorities,
availability and cost of business premises, references from local partners and previous
experience, cost of construction land, availability of construction land, level of political
interference in business, ecology utility costs accommodations (Aleksandar et al, 2016).
As it was stated above, in Serbia, industrial sites were not selected as simply and only vacant
land availability is not enough to establish industry. Different accessibilities have to exist in
the area in which industry located. Besides to physical factors that determining site selection,
government commitment and security issues were also highly determining private investors
14 | P a g e
to invest in the given area. In addition, previous experience of industrial establishment is also
one of criteria to select industrial sites, which is important to take a lesson in order to identify
failures and achievements. This would help the newly establishing industry from being
stagnant and failure. Therefore, a given administration has to desire a specific criterion that
ensures environmental, social wellbeing and security issues of the communities during
industrial site selection phases.
Africa. As a result of our trip, we identified a great opportunity for the industry to invest nearly a
billion dollars to create a vertically integrated apparel supply chain in some regions of Africa. In
other words, create jobs in Africa not just in apparel, but also in the textile and even agricultural
industry through better cotton production techniques (McRaith, 2015).
The apparel industry’s mistakes of the past have often been the result of near-sighted investment
in lawless environments. That model must and will change rapidly countries in Africa where we
invest will be the beneficiaries of a new and more inclusive model of investment and growth in
which we will be able to put in place, right from the beginning, facilities, norms and values that
will guide the work at the factories and the relationships between workers, managers,
associations, civil society groups, governments and any other stakeholders. That is why when we
look at Africa we do not just look for a place to quickly set up a sewing operation.
So, it simple to understand in that establishing industries and inviting private investors on
industries requires a clear process and regulations to be effective and productive. But as stated
above, Africans did not conduct impact identification mechanisms before establishing and
15 | P a g e
participating on industries. Therefore, Africans have to conduct impact mitigation analysis
before establishing industries.
16 | P a g e
local entrepreneurs, diversify the structure of the country’s export, and generate employment.
The IPDC, which has the vision of becoming an “innovative and leading eco-industrial parks
developer and operator in Africa by 2025, is mandated to activate both pre and post investment
servicing, availing land, and pre-built sheds equipped with all-encompassing utilities and
infrastructural facilities with international standards of quality of service, labor security, and
environmental safety.
About 12 industrial parks for export processing, have been identified across the country based on
proximity to market outlets, infrastructure, economic potential, and regional balance in
development. Almost all industrial parks are established along corridors of railway networks that
are either completed or under construction and are on schedule. The construction of two major
industrial parks in Addis Ababa and Hawassa has been completed with nine firms having started
production in Addis Ababa and Hawassa (Kenichi, 2017)
The Ministry of Industry of the Federal Democratic Republic of Ethiopia has prepared three
documents on the industrial strategy, roadmap to industrial development, and institutional setup
for the implementation of the industrial development. These documents are believed to lead the
country’s direction on industrial development for the next fifteen years. The industrial
development objectives, goals and strategies on the GTP II document are therefore the
reflections and highlights of the more elaborated documents prepared by the Ministry of
Industry.
18 | P a g e
and construction of Hawassa will be conceived around energy and water conservation principles
including maximization of natural lightning and natural ventilation, fitting of low consumption
bulbs, recycling of rain water, and solar powered led street lights making it Ethiopia’s most
important eco-friendly site. In Phase 1 cycle 1, 100 hectares of land were developed.
19 | P a g e
ii. The costs of transportation and
iii. The cost of labor
Of these, the relative price differential of raw materials from various sources can be expressed in
terms of transportation costs. For theoretical purposes, it is equivalent to thinking in terms of
cheap deposits situated relatively near to the plant and dear deposits further away. Consequently,
we can work with two general regionally operating location factors, namely the cost of
transportation and cost of labor. Being strictly function of space, regional factors can all be
analyzed from the point view of the individual isolated production process.
All the other factors of location work between industries and therefore are not to be found in any
examination of an isolated production process. They are grouped together under the title
"agglomeration factors" and work to create groupings of industrial in agglomerations of various
sizes. The aggregate of agglomerative factor is treated as a "uniform agglomerating force". This
is the third and last of the general location factors.
From this concept, the researcher understood that the relative price of raw materials and its
transportation cost determines industrial establishment, i.e. it the area is not accessible, private
investors will not have interest to invest on that inaccessible areas. Therefore, proximity to the
accessible road is important for industrial development.
20 | P a g e
Thus, inter- industry linkages and theory of industrial interdependence play a major role in the
growth pole theory. As a matter of fact, this theory together with the Schumpeterian theory of
development can be said to constitute the two cornerstones upon which Perroux bases his theory.
The external economies that become available in the area constituting the growth pole of a
region are basically of the following three types:
Economies internal to the firm: These are the lower average costs of production resulting
from an increased rate of output. These are the economies which any single firm by their
own organization and effort can enjoy.
Economies external to the firm but internal to the industry: These are associated with
localization of industry. On account of close location proximity of linked firms, as
industry expands at a particular location, cost per unit of output to a firm declines.
Economies external to the industry but internal to the urban area: These can be termed
urbanization economies. They include development of urban labor markets, access to a
larger market, and provision of a wider range of services.
The growth pole theory is a useful concept for spatial analysis. The theory is attractive to policy
makers because it presents opportunities for integrating industrial policy, physical planning and
inter-regional and intra-regional economic planning.
In this concept, the existence of one industry may attract the establishment of another industry.
For example, if there is car production/seller in a given area, garages may interest to open at near
there to obtain customers easily. So, industrial site selection also considers the existence of other
industries near or far away from the newly desired industrial establishment.
21 | P a g e
lowest delivered price and the boundary between the market area of two producers will be a line
joining the points at which delivered price is the same from both sources.
As long as the cost of extraction does not vary with output, transport costs are the only variable
affecting price, but Hoover extends his analysis to include the influence of diminishing returns to
scale. He follows Weber, fairly closely at first, pointing that in the absence of production cost
differences the best location will be at the point of minimum transport costs, which may be at a
material source, at the market or at an intermediate point: The least-transport-cost location is
found by drawing isotims around given material and market points, from which lines of equal
total transport cost 1 can be constructed. But Hoover goes further than Weber by showing how
different sections of the market will be served by different producing points.
He also differs with Weber's emphasis on least transport cost points within location triangles.
Even with the assumption of uniform transport costs, the possibility of a separate minimum point
not at one corner of the triangle is much less than might be thought at first sight. It is far more
likely than what Weber suggested that a material or the market will have a pull which will be
greater than that of the other corners. The chance of a location not at one corner is even less
likely, if the fact that transfer costs are actually less than proportional to distance is also
considered. In addition, loading costs and other terminal changes operate against least cost
location inside the triangle. If a separate point away from material sources and market does
occur, Hoover suggests that perhaps this is a sign that industry is not primarily transport oriented
at all and that possibly a low labor cost location enters into the picture.
In this theory, transportation only does not determine industrial location, but also labor cost also
determines the effectiveness of industries. That means, even though there is transport
accessibility in the area, shortage of labor force would impose negatively the effectiveness of
industries. This directly related with the multi-criteria decision support system in that only one or
two factors did not give a clue whether the desired site is optimum or not.
22 | P a g e
2.6. Conceptual Framework of the Study
Different studies have been undertaken in how to apply GIS based models in selecting industrial
site location. Those studies tried to indentify factors and constraints that hinder suitable industrial
location. In our country, updating the workflow of software is still lagging behind when compared
with other countries experience.
Appropriately evaluating physical factors and constraints allows cities to tactically attract
developers with minimum complication. The decision process identifies optimal sites while
taking into account economic benefits and environmental sustainability.
The following factors, assessed by proximity to industrial land, are sewer service areas, fiber-
optic networks, and public transportation. The model inventories and prioritizes industrial land in
a two-part process. First, zoning policies and GIS data identify town districts that allow industrial
activity. After selecting these zones, a constraint analysis removes all undevelopable land in each
area. The purpose is to eliminate any existing development hurdles, according to environmental
regulations and development standards, inhibiting the development process, in order to improve
project feasibility.
The second part of the procedure, also conducted in GIS, applies a factor analysis, which
incorporates specific variables essential for industry. A factor analysis examines the spatial
relationship between developable industrial land and proximity to existing infrastructure
amenities (or factors). The closer the distance of industrial land to each factor, the more suitable
the location becomes. The weighted overlay analysis then measures suitability for development
and ranks locations as either having a high priority or a low priority for protection and
investment. This process applies weighted values to each factor, which vary based on specific
requirements or business needs.
Finally, after suitability is measured, a more detailed analysis identifies highest priority industrial
parcels not adjacent to residential parcels. Large contiguous parcel sizes are compatible for
future investment; they allow a community to readily absorb industrial growth with the least
amount of restrictions. Communities can efficiently use existing infrastructure to take advantage
of investment opportunities for green technology, niche manufacturing, or eco-industrial parks.
These spaces are also critical locations, positioned at the crossroads of major arterials, within
utility service areas, and in walking distance to public transportation systems. With respect to
industrial land, proximity to residential areas is least desirable.
23 | P a g e
Determining Factors to Select Industrial sites
Dependent variable is: suitable site for industry. The characteristics and nature of the
independent variables determines the characteristic of dependent variable. In short, site selection
for industries primarily depends on the above listed independent variable selection mechanisms.
24 | P a g e
CHAPTER THREE
Hawassa city is bounded by Lake Hawassa in the West, Oromia Region in the North,
Wendogenet woreda in the East and Shebedino Woreda in the South. The town administration
has an area of 157.7 km2. The town administration is divided in to 8 sub-cities and 32 Kebeles.
These Eight sub cities are Hayk Dar, Menehariya, Tabor, Misrak, Bahil Adarash, Addis Ketema,
Hawela-Tula and Mehal sub city (Ayele, 2017).
25 | P a g e
3.1.2. Physiography and drainage
Hawassa is a town in Ethiopia, on the shores of Lake Awasa in the Great Rift Valley. Other
prominent volcanic features surrounding the city are Hill Tabor and Hill Alamura and spring
water(http://www.geology.cz/projekt681700/vystupy/Hawassa_subsheet_0738_C4.pdf,2019:4:1
8:24 am).
Hill Tabor: Tabor hill has majestic presence in the western end of the city. This attractive hill is
perfect for a walk in the afternoon or early morning, and offers beautiful panoramic views across
the whole area. Said to resemble Tabor Mountain in Israel, it sweeps dramatically down to the
lakeshore, mirroring the ridge of the mountain on the opposite side of the lake Hawassa.
Hill Alamura: It is a small mountain situated at the south western end of the town. From the top
of the mountain, one can enjoy a partial view of the town and the lake stretched below it. Paying
a visit to this mountain and ascending to its top is worth a day of any one’s time.
Spring Waters: There are three natural heritage sites (spring waters) surrounding Hawassa
town. These are: Burqito (12km from Hawassa city), Wendogenet located 39 km from Hawassa
and Gidabo 41 km south of Hawassa.
26 | P a g e
35
30
25
Temperature in 0C
20
Mean of
15 MaxTem
10 Mean of Min
Tem
5
Mean Temp
0
jan feb mar apr may jun jul aug sep oct nov dec
Mean monthly RF
180
160
140
120
RF in mm
100
80
60
40
20
0
jan feb mar apr may jun jul aug sep oct nov dec
27 | P a g e
3.1.4. Population of the study area
Based on the 2007 Census conducted by the Central Statistical Agency of Ethiopia, the study
area has a total population of 258,808, of whom 133,123 are men and 125,685 women. A total of
61,279 households were counted in this town, which results in an average of 4.22 persons to a
household, and 57,469 housing units. Much of the population growth has been the result of
internal migration and expansion of educational and other facilities, also widening of the city’s
boundaries has caused some of the increase. Hawassa has a young population. Around 65% of
the people are under 25 years of age, and only about 5.5% of the population is over 50 years of
age. (CSA, 2007).
Hawasa basalt: it forms a ridge transecting the Hawasa Caldera in a NNE-SSW direction. The
ridge subdivides the caldera into the western part (with Lake Hawasa) and the eastern part (with
Chaleleka Lake –Shashemene subsheet). Lavas of the Hawasa basaltic ridge were emitted from
several scoria and tuff-cones arranged along the ridge.
Wendo Koshe obsidian IV: In the accessible eastern part of the Wendo Koshe Volcano,
obsidian lava flows were documented beneath the younger obsidian and pumiceous
accumulations of the older Wendo Koshe Volcano.
28 | P a g e
Wendo Koshe pumice fall and minor flow deposits: The deposits are mostly well sorted and
clast-supported suggesting fall deposition. These were produced during the last big explosive
event of the Wendo Koshe Volcano. The pyroclastic (mostly fall) deposits cover a large area of
the northern part of Hawasa subsheet and extend as far as Bura and Shashemene.
Along the Shashemene to Alaba-Kulito road (just north of the subsheet), 0.5 m
of young Wendo Koshe pumice overlays a 6 m thick accumulation of ochre fine ash with
abundant accretionary lapilli of an unclear source. Even along the Hawasa-Shashemene road, the
thickness of the young Wendo Koshe pumice reaches 2 m. The exact age of this violent eruption
has not been determined, as no charcoal was found at the base or within this deposit (JICA
(2012).
29 | P a g e
Organic soil material has organic carbon content (by weight) of 12 to 18 percent, or more,
depending on the clay content of the soil.
Eutric Fluvisol: these soils are genetically young soils. They are zonal soils in alluvial deposits,
which found in alluvial plains, river fans, valleys and tidal marshes on all continents and in all
climate zones. Many Fluvisols under natural conditions are flooded periodically. Most
Fluvisols are wet in all or part of the profile due to stagnating groundwater and/or flood
water from rivers or tides. So they are poorly drained soils.
Chromic Luvisols: is a soil with a clearly developed clay leaching horizon and a high cation
exchange capacity . The parent material of Luvisolen consists of unconsolidated material.
Luvisols mainly occur in landscapes with a flat to slightly hilly relief with a cool
temperate climate or a warm climate with clearly dry and wet seasons.
Eutric Vertisol: they are dark colored cracking and swelling clays. Vertisols are soils with a
high content of clay minerals that shrink and swell as they change water content. The clay
minerals adsorb water and increase in volume (swell) when wet and then shrink as they dry,
forming large, deep cracks. Surface materials fall into these cracks and are incorporated into the
lower horizons when the soil becomes wet again.
30 | P a g e
3.1.7. Infrastructure in the Town
The town finance and economic development department report on socio-economic profile 2015
of the town states that, the economic growth and overall development of the town depends on the
presence and improvement of infrastructure facilities. These facilities are Road, transportation.
Road Network: Hawassa City has transportation service since 1952. During the time
transportation activity was done by Ministry of works. First, gravel roads constructed followed
by Asphalt roads after the Municipality established.
Road infrastructure of the City has made a radical change since 2004. Upgrading of the gravel
roads to an Asphalt road (Arterial Asphalt Road) has began in 2004, and made a significant
change in road infrastructure of the City.
Transportation: Transport systems in Hawassa include the buses, Taxis/Bajaj/, Motor bicycles,
Carts and other types of vehicles. During the past years there has been tremendous growth in the
number of vehicles specially Bajaj Taxis in the city, resulting in immense traffic pressure and
increasing road bottlenecks. So, road network and transportation availabilities also considered as
one of those factors that determining industrial site selection, the researcher will give high
attention on those factors in this research study.
31 | P a g e
3.1.8. Industrial Development in the Study Area
The Ethiopian economy has recorded continuous and double-digit growth since 2003/04. The
service sector outperformed other sectors in this period, and its share of GDP increased from
39.7% in 2003/04 to 46% in 2009/10. Despite 10 per cent annual average growth in the same
period, the agriculture share in GDP declined from 47% to 42%. The industry sector also grew
by about 10% per annum between 2003/04 and 2009/10. However, its share in the economy
remained relatively static (12–14%) over the last decade (WB, 2013).
The following industrial types were currently operating in the town: grain mills, house and office
furniture, textile, soap factory, stone crusher, printing, laundry, bakery, HCB manufacturing,
metal works, oil mills, industrial parks, etc (Hawassa City Trade and Industry Office, 2017).
And also, the data collected by the hand GPS about the distribution of main industrial types
shows that industries condensed almost all in one parts of the town as described in the figure 6
below. In the future, it needs a high attention to focus on those newly incorporated structural plan
of the town. Especially at the fringe of the town and the surroundings, it is possible to establish
industries and in doing that it is had better to avoid industrial noises emission at the central part
of the town.
32 | P a g e
3.2. Research Methods
3.2.1. Types of data used for this study
Both primary and secondary data were used for the industrial site mapping of the study
area, which were obtained from field survey and concerned institutions. The data and materials
used include satellite imagery Landsat8), climatic data (rainfall and temperature), Digital
Elevation Model (DEM) 30 m resolution data, soil map and industrial data of the
study area. Types of data used and their sources are described in Table 2 below.
ERDAS IMAGINE 2015 Land use land cover classification and for accuracy assessment
ENVI 5.3 Land use land cover image classification for accuracy assessment
33 | P a g e
3.2.3. Suitability requirements rating for the selected factors
Table 4: Suitability requirements rating for the selected factors in the study area
Parameters Units Suitability Range
Unsuitable(1) Less- Moderately Suitable(4) Most
suitable(2) suitable(3) suitable(5)
LULC classes type Settlement recreation Business site shrub land Vacant
land
Soil type type Eutric Chromic haplic luvisol Fabric Eutric
vertisol luvisol histosol fluvisol
Distance from meter 0-2000 2001- 4000 4001- 10000 >1000 -
geological fault
Distance from meter 0-200 201 -500 501-1000 1001- -
lake 10000
Distance from meter 0-250 251-500 501- 1000 1001-2000 >2000
river
Distance from meter 0-250 251-500 501- 1000 1001-2000 >2000
swamp
slope % 0-6 7-10 11-15 15-20 >2-0
Distance from meter 0-1000 1001-1500 1501-3000 3001-6000 >6000
residential area
Proximity to road meter 0-300 301-500 501-3000 3001-5000 >5000
Source: Sarath et al, 2018, Ebraim et al, 2015, Amita Johar 2013, FAO, 2009
34 | P a g e
Figure 10: Slope map of the study area
35 | P a g e
Figure 11: Lake map of the study area
Proximity to rivers
Traditionally, industrial land is located along rivers and canals because of its proximity to
inexpensive hydropower (Mark, 2014). Sites overlapping or adjacent to these areas pose greater
insurance costs and risk potential flooding, requiring expensive flood mitigation strategies like
barriers and retaining walls. Therefore, development within river zones is an unattractive
characteristic and considered an industrial location constraint.
The river through Tikur Woha and Wosha River the way to Wondegenet was digitized from
Google earth. The map showing proximity to rivers is generated by using ring buffer distance in
spatial analyst tools of Arc Map from digitized rivers and canals. The reclassification of rivers
map is based on suitability criteria. Areas found 0- 250 m, 250 – 500 m, 500 – 1000 m, 1000 –
2000 m and above 2000 m from rivers. The assigned new values as 1, 2, 3, 4 and 5 and
classified as unsuitable, less suitable, moderately suitable, suitable and most suitable,
respectively.
36 | P a g e
Figure 12: River map of the study area
Proximity to swamps
As it was mentioned above, swamps were digitized from shapefile data obtained from FAO 2006
and the map was generated in ArcMap 10.3.1 environment by using ring buffer distance
calculation. The reclassification was done on the basis of industrial site selection criteria. New
values were assigned for regions located within 0 – 200 meter, 200 - 500m, 500 – 1000 m, 1000
– 2000 m and above 2000m. The assigned new values were 1, 2, 3, 4 and 5.
37 | P a g e
3.2.5. Land use land cover factor
Land use situation with regard to different environmental effects and the condition of industrial
establishment is very important. Forest and farming lands shouldn’t change to industrial land use
and only a negligible damage to the plant cover can be allowed. Undersigned, in boundary, open
space and free lands is preferable (Eldin etal, 2013).
The reclassification was done on the basis of industrial site selection criteria. New values were
assigned for regions located within forest and water bodies are restricted it does not be represented with
value whereas, settlements assigned as new value 1, recreation center as 2, business sites as 3, shrub lands
as 4 and vacant land as 5. Then reclassified as unsuitable, less suitable, moderately suitable, suitable
and most suitable, respectively.
Figure 14: Land use land cover map of the study area
3.2.6. Road network
Roads and other infrastructures that facilitate different cities and regions connection are very
important in industrial site selection. According to the data obtained from Ebrahim et al, (2015),
industries needs transport facility and therefore its establishment should be near to the road. The
reclassification was done on the basis of industrial site selection criteria. New values were assigned
for regions located within 300 meter as 5, 301 – 500 meter as 4, 501 m – 1000 meter as 3, 1001 - 5000
meter as 2 and above 5000 meter as 1.
38 | P a g e
Figure 15: Road network of the study area
3.2.7. Residential area factor
The reclassification was done on the basis of industrial site selection criteria. New values were
assigned for regions located within 1000 meter as 1, 1000 – 1500 meter as 2, 1500 m – 3000 meter as 3,
3000 – 6000 meter as 4 and above 6000 meter as 5. Then, the newly reclassified value represented as
unsuitable, less suitable, moderately suitable, suitable and most suitable respectively.
39 | P a g e
3.2.8. Geological fault factor
The Awassa town has a volcanic origin. It was formed due to intensive tectonic activity. This
tectonic activity is a direct result of the Main Ethiopian Rift valley. The geology of the Awassa
area consist of different volcanic deposits such as alkaline and peralkaline rocks Late Miocene,
Basaltic lava flows from Pleistocene to recent, Acidic volcanics from Pleistocene to recent, and
Volcano clastic lacustrine sediments from Pleistocene to recent (Alemayehu A, 2008).
Distance to geological fault is one of the hindering factors for industrial site selection. The more
far away from this event is the better. The map of fault distribution was constructed from
geological maps. The geological fault zones of the study area were digitized from Google Earth
pro and converted to the shapefile data format in order to undertake further analysis.
The reclassification was done on the basis of industrial site selection criteria. New values were
assigned for regions located within 2,000 meter as 1, 2,001 – 4,000 meter as 2, 4001 m – 10,000 meter as
4, and above 10,000 meter as 5. Then, the newly reclassified value represented as unsuitable, less
suitable, suitable and most suitable respectively.
Figure 17: Map of distance from geological fault of the study area
40 | P a g e
3.2.9. Soil factor
According to Anurag Ohri et al (2010), those soils with clay, alluvial and loam deposits are
suitable for industrial establishment, whereas soils with cracking, swelling and porous are not
suitable for industrial establishment.
The reclassification was done on the basis of industrial site selection criteria. New values were
assigned for regions located within eutric vertisol as 1, chromic luvisol as 2, haplic luvisol as 3, and fabric
histosol as 4 and eutric fluvisol as 5. Then, the newly reclassified value represented as unsuitable, less
suitable, moderately suitable, suitable and most suitable respectively.
Direct and indirect unstructured interviews were conducted with the experts during the field
survey to gather more information. The information derived from this study was used to identify
and develop priority criteria and factors for the selection of industrial site. It was also used to
41 | P a g e
identify problem in the study area and prioritize the potential industrial sites in the study area.
Additional discussions were made particularly with those experts, who are involved in the urban
land preparation of Hawassa town, in order to acquire some appreciation of the study area. In
particular, the following experts working at the town Industrial and Enterprise department,
Housing and Urban Development Department and finally with Regional Bureau Industry experts
discussion were held. Those who were involved in the discussion were interviewed: Ato
Matusala Semaye, Yaregal Debebe, Sileshi Getahun and Walelign Awashe.
This research study was mainly used spatial Multi-Criteria Decision Making (MCDM) approach
integrated with Geographic Information System (GIS) techniques in order to determine suitable
site for industry. Spatial multi-criteria decision making is a process that combines and transforms
geographical data into a decision (Malczewski, 1999).
MCDM, combined with GIS data, is a powerful approach to systematically and comprehensively
analyze a problem. The fundamental components of a multi-criteria problem are human value
judgment and assessments of the importance of criteria. The main purpose of the multi-criteria
evaluation techniques is to investigate a number of alternatives in the light of multiple criteria
and conflicting objectives (Voogd, 1983).
To combine and determine the importance of each factor, multi-criteria decision making method
was applied.
42 | P a g e
criterion is less important than the column criterion. To ensure the credibility of the relative
significance used, AHP also provides measures to determine inconsistency of judgments
mathematically. Based on the properties of reciprocal matrices, the consistency ratio (CR) can be
calculated. CR < 0.1 indicates that level of consistency in the pair wise comparison is acceptable.
Saaty (1980) suggests that if CR is smaller than 0.10, then the degree of consistency is fairly
acceptable. But if it is larger than 0.10, then there are inconsistencies in the evaluation process,
and AHP method may not yield meaningful results.
1/2, 1/3, 1/4., 1/5, 1/6, 1/7, 1/8 and 1/9 Reciprocal values of the previous appreciation
43 | P a g e
For each environmental factors mentioned above, maps were generated and then reclassified
depending on their suitability for industrial sites. In the process of reclassification, new values 1,
2, 3, 4 and 5 were assigned and they represent unsuitable, Less suitable, moderately suitable,
suitable and most suitable respectively. Finally, weighted overlay carried out using ArcMap
10.3.1 to generate the final industrial map of the study area.
The overall work flow used to accomplish the study was presented in the figure (20) below:
Proximity to:
Land use Slope suitability Soil type Water body /lake map, River/streams
land cover map suitability based map, Swampy based map, Road
suitability Map map, and Residential area based map
map and Geological fault suitability map
Weight Computation
Weighted Overlay
44 | P a g e
CHAPTER FOUR
The following factors were deeply discussed in this part. These are: Land use land cover types,
distance from geological faults, distance from main roads, soil types, and distance from
residential areas, distance from water bodies (lake, swamp and river/stream), slope and elevation.
In this study slope factor was generated from the Digital Elevation Model (DEM) using the
ArcGIS 10.3.1 spatial analyst extension of surface module, which enabled to classify the area
according to the steepness and the gentleness of the terrain. Every cell in the output raster had a
slope value. The lower the slope value, the flatter the terrain was and the higher the slope value
the steeper was the terrain. Then the slope raster was reclassified in to five classes of slope
percent by examining the value and the frequency of slope percent in the study area table 6)
The reclassified slope was given a rank value 1 to 5 with the higher value of 5 showing high
influence, i.e. highly suitable, while the lower value of 1 showing low influence, least
suitable. According to Ebrahim et al, 2015, slope with degree ranges from 0-6 considered as level
slope, 6-10 degree as gentle, 10 -15 degree as strong sloppy, 15- 20 degree as moderately steep
slope and > 20 degree slope as considered as steep slope). As indicated in the figure 21 below, most
of (68.6%) the study area has gentle slope and 3.2% has a moderately steep slope
45 | P a g e
Table 6: Reclassified Slope Value and Suitability Level
Slope in Degree Rank Area (km2) Area (%) Suitability Level
>20 1 9.8 6.2 Unsuitable
15 - 20 2 5.1 3.2 Less suitable
10 - 15 3 11.7 7.4 Moderately Suitable
6 – 10 4 23 14.6 Suitable
0-6 5 108.2 68.6 Most Suitable
Total 157.7 100
Source: Extracted from Reclassified Slope, 2019
The less the slope amount, the more suitable will be the area for the industry. Based on the
suitability of the slope for industrial sites, the reclassified slope map shows 108.2 km² (68.6 %)
area is most suitable, 23 km² (14.6%) is suitable, 11.7 km² (7.4%) is moderately suitable, 5.1 km²
(3.2%) is less suitable and 9.8 km² (6.2%) unsuitable for industrial location.
From the total area of the study, majority of the slope initiates public and private investors to
participate on industrial sectors in which Plain areas have especial importance because of lesser
costs during site preparation.
46 | P a g e
4.2.Proximity to Water Bodies and Industrial Locations
Water accessibility is one of the main affecting factors for industrial site selections. Industry
must be located at some safe distance from water bodies. Water bodies and areas within 1000
meter are considered unsuitable for an industrial establishment. Waste water from industry is one
of the major sources of pollution in river distance of at least 1000m (Jain, 2008). According to
Ebrahim et al (2015), land suitability for industrial location within proximity to lake as follows:
The proximity to water bodies map is the result of overlaid proximity to Lake Hawassa and
Chelelka maps. As indicated in the figure below, 6.8 km2 of the total area is within 200 meter of
47 | P a g e
the water bodies, 8.1km2 is found between 200 meter up to 500 meter, 21.3km2 is found
between 500 meter up to 1000 meter, and the remaining 121.2km2 found beyond 10000 meter.
This implies that most parts (76.8%) of the study area are suitable, 13.5% is moderately suitable,
5.2% is less suitable, and 4.3 % is unsuitable for industrial location.
48 | P a g e
Figure 22: River suitability map of the study area
According to Amita Johar (2013), land suitability for industrial location within proximity to
swamps as follows:
Table 9: Reclassified Proximity to Swamp and Suitability Level
Distance from Swamp (m) Rank Suitability Level
0 - 250 1 Unsuitable
250 - 500 2 Less suitable
500 - 1000 3 Moderately Suitable
1000 - 2000 4 Suitable
>2000 5 Most Suitable
Source: Extracted from Reclassified Swamp Value, 2019
49 | P a g e
Based on the proximity to swamps, the majority of study area is most suitable for industrial
location. As described in the table below, 94.5 km2 (60%) of the total area is most suitable and
34.7 km2 (22%) is unsuitable to locate industry in the study area. Besides that, 3.7km2 (2.3%)
less suitable, 7.7km2 (4.88%) is moderately suitable and 16.6km2 (10.5%) is suitable to locate
industries in the study area. Figure (24) below portrays proximity to swamps based industrial site
selection.
Based on the supervised classification done for the satellite image, the study area has seven
major land use land cover types, such as recreational center, shrub lands, water body, vacant
land, forest coverage, commercial or business centers and residential/settlements. Mixed type of
land use is dominant, especially residential development mixed with commercial and other
activities which are common in the central parts of the study area. In general, the majority (65.7
50 | P a g e
km2) of land use land cover map of the study area was vacant land since newly desired master
plan of the town incorporates the surroundings of Sidama zone and 32.2km2 of the total area is
residential/settlements and 30.6 km2 of the total area is shrub lands and 17.7km2 is about forest
coverage. On the other hand, a small portion of the study area was covered by water bodies,
recreational centers and business activities.
Business
%
Vacant
shrub
water
forest
Total
Settlement 18 1 0 2 3 0 0 24 0.2
Business 2 13 0 1 0 1 0 17 0.03
Recreation 0 0 7 0 1 0 1 9 0.03
Vacant 1 0 0 115 2 0 0 118 0.4
shrub 0 0 1 0 20 1 0 22 0.2
forest 4 0 0 1 2 40 0 47 0.11
water 0 0 0 0 0 1 2 3 0.02
Total 25 14 8 119 28 43 3 240 1
Source: Compiled by the Researcher, 2019
51 | P a g e
4.3.2. Error matrix of the estimated area proportion
The proportion of mapped area (wi) is calculated as the ratio of mapped area to the total area of
the study area. Eg: - The proportion of forest mapped area: 17.75/157.7 = 0.11
An error matrix can be reported in terms of sample counts. However, error matrix is usually
expressed in terms of (the unbiased estimator) estimated area proportion rather than in terms
of sample counts. After this the area of each land use/cover class is estimated based on the
reference classification, which can be calculated using the equation proposed by Olofsson et al.
(2013) as follows:
Where, wi: (the proportion of area mapped as land use class i that means mi/mt);
mi: is the mapped area of land use class i, mt: is the total area of the map, nij: is the cell value at
row i and columun j; ni.: is the total sample count of land use class i(total number of samples for
the mapped category but not for the reference category).
For example: e.g: for Forest. Wi=17.7(Mapped area for forest/157.7(total area of the study area)
= 0.11; nij=40; ni.=47 pij= 0.13*40/47= 0.1 as stated in the table (11) below.
LULC
Recreation
Settlement
Producer
Business
Mapped
Vacant
shrub
water
forest
Total
Area
User
Settlement 0.15 0.008 0 0.002 0.025 0 0 0.17 32.14 0.937 0.88
Business 0.0035 0.024 0 0.002 0 0.002 0 0.03 3.6 0.75 0.8
Recreation 0 0 0.022 0 0.003 0 0.003 0.026 4.06 0.733 0.846
Vacant 0.003 0 0 0.4 0.007 0 0 0.42 65.5 0.983 0.95
shrub 0 0 0.008 0 0.18 0.008 0 0.19 31.2 0.818 0.947
forest 0.1 0 0 0.03 0.005 0.11 0 0.13 17.7 0.846 0.846
water 0 0 0 0 0 0.01 0.026 0.034 3.5 0.896 0.764
Total 0.166 0.032 0.03 0.407 0.22 0.13 0.029 157.7
Overall Accuracy: 0.856
Overall Kappa Statistics = 0.7941
Source: Compiled by the Researcher, 2019
52 | P a g e
Estimated Land Use Land Cover and confidence interval
Confidence interval for the estimated area of each land use class should be
estimated after the standard error of each land use class is calculated. The estimated standard
error (SE) of the estimated area proportion (the column total) of each land use class is
estimated as Olofsson et al. (2013):
pj refers to land use class j (based on reference class). The standard error for the error
adjusted estimated area of each land use class is calculated as the product of the proportion
area of a given class and the total map area.
Table 12: Estimated Land Use Land Cover and confidence interval
53 | P a g e
Table 13: Area estimates of each land use class based on reference data (Km2)
Land use land cover class Estimated Mapped
Forest 28.18 17.75
Vacant/open land 4.5 65.74
recreational Center 4.7 4.06
Residential area/settlement 64.8 32.14
commercial site/business 34.7 3.59
water body 20.3 3.2
shrub land 4.6 30.6
Total 157.7 157.7
Source: Compiled by the Researcher, 2019
Figure (25) below shows land use land cover based industrial site suitability. According Amit
Johar 2013, water bodies, swamps, forest coverage and residential areas are not suitable for
industrial establishment. Whereas areas like open spaces, vacant lands, range lands area suitable
for industrial site location.
Table 14: Land Use Land Cover Classes and Suitability Level
54 | P a g e
Figure 24: Land use land cover suitability map of the study area
As indicated in the table (15) below, road distance in meter and its suitability level used for this research
study.
55 | P a g e
Table 15: Road Distance and Suitability Level
Access to road (in meter) Rank Suitability level Area in km2 Area in %
0-300 5 Most suitable 26.6 16.8
301-500 4 Suitable 15.2 9.6
501-1000 3 Moderately suitable 32.3 20.4
1001 - 5000 2 Less Suitable 76.4 48.4
>5000 1 Unsuitable 4.8 3.04
Total 157.7 100
Source: Compiled by the Researcher, 2019
From the total area, about 26.6 km2 (16.8%) is most suitable for industrial sites and only 4.8km2 (3.04%)
is unsuitable for the industrial establishment. About 76.4km2 (48.4%) of the study area is less suitable to
settle industry, which is related transport accessibility. As newly proposed structural plan of the town
stretched to the fringe of Sidama Zone, it requires road construction in the future. The figure (26) below:
road suitability map of the study area.
56 | P a g e
4.5.Distance from the Residential Area and Industrial Location
A logical distance should be selected between residential areas and industrial sites. A buffer
should be defined as well. These areas are barriers for industries development. The more distance
from industrial areas the higher suitability and consistency will be.
According to the data obtained from Ebrahim et al (2015), industries needs to be established far from
the residential areas. These areas are barriers for industries development. The more distance from
industrial areas, the higher suitability and consistency will be. As indicated in the table (16) below,
residential areas proximity to industrial sites was described for this research study.
Table 16: Residential Area Suitability
Distance from residential areas Rank Suitability Level Area in Krm2 Area in %
0-1000 1 Unsuitable 70.03 44.4
1000-1500 2 Less Suitable 14.5 9.2
1500-3000 3 Moderately Suitable 40.5 25.7
3000 - 6000 4 Suitable 30.4 19.3
>6000 5 Most Suitable 2.07 1.4
Total 157.7 100
Source: Compiled by the Researcher, 2019
The buffer distance zones have been categorized into five levels based on the level of distance
from residential areas. Accordingly, the low buffer distance ranked as unsuitable whereas the
longer buffer distances ranked as most suitable. Thus, the rank value of 5 was given for highly
suitable road buffers and the rank value 1 was given for unsuitable road buffers.
As indicated in the table above, from the total area 70.03km2 (44.4%) is unsuitable as the area
occupied by the residents of the town and only 2.07km2 (1.4%) is most suitable for industrial site
selection. And also, 40.5km2 (25.7%) is moderately, 40.2km2 (19.3%) is suitable and 14.5km2
(9.2%) is less suitable for the site selection of the industry. Figure (28) below portrays residential
area suitability for the industrial site selection.
57 | P a g e
Figure 26: Residential area suitability map
Faults are one of the geological phenomena that cause problem for industrial buildings and
infrastructure construction. Industries should be 2000 m away from faults. The area that are
within 2000 meter are unsuitable, 2001 – 4000 meter are less suitable, 4001 – 10000 meter are
suitable, and area >10000 meter are most suitable to establish industries (Marzieh et al, 2011).
58 | P a g e
Table 17: Geological Fault Suitability
Distance from residential areas Rank Suitability Level Area in Krm2 Area in %
0 - 2000 1 Unsuitable 32.4 20.5
2000-4000 2 Less Suitable 50.1 31.7
4000-10000 4 Suitable 74.9 47.5
>10000 5 Most Suitable 0.3 0.2
Total 157.7 100
Source: Compiled by the Researcher, 2019
59 | P a g e
Luvisols: Soils with subsurface accumulation, low activity clay, high basic saturation, are porous
and crumb characteristics.
Vertisols: they are dark colored cracking and swelling clays.
Histols: Soils which are composed of organic materials.
According to Anurag Ohri, et al, 2010 report, those soils with clay, alluvial and loam deposits are
suitable for industrial establishment, whereas soils with cracking, swelling and porous are not
suitable for industrial establishment. Based on this information, the researcher classified soil
types as described in the table (18) below:
Table 18: Soil Type and its Suitability
Soil Types Rank Suitability Level Area in Krm2 Area in %
Eutric Vertisol 1 Unsuitable 1.4 0.9
Chromic Luvisol 2 Less Suitable 11.8 7.5
Haplic Luvisol 3 Moderately Suitable 18.2 11.5
Fabric Histosol 4 Suitable 31.6 20.03
Eutric Fluvisol 5 Most Suitable 92.5 58.7
Total 157.7 100
Source: Compiled by the Researcher, 2019
60 | P a g e
As indicated in the table (18) above, from the total area 92.5km2 (58.7%) is most suitable and
only 1.4km2 (0.9%) is unsuitable for industrial site selection as the area filled by the Eutric
vertisol. And also, 18.2km2 (11.5%) is moderately, 31.6km2 (20.03%) is suitable and 11.8km2
(7.5%) is less suitable for the site selection for the industry. Figure (29), above portrays soil type
suitability for the industrial site selection
Before combining them, the following procedures were taken place; first, rasterization was done
for the vector data layers in order to produce similar data layers to perform GIS analysis and
secondly standardization of each data set to a common scale of 1 to 5 was done in ArcGIS
Software. Prior to combining the factors, weights have to be given based on Satty's Analytic
Hierarchy Process (AHP), where a pair-wise comparison matrix will be prepared for each map
using a nine point importance scale (Table 19 ). Weighting is used to express the relative
importance of each factor relative to other factor. The larger the weight, the more important is
the factor in overall utility.
Moderately
Extremely
Extremely
Important
Strongly
Strongly
Strongly
strongly
Equally
Very
Very
Less More
Important Important
Source: Adapted from Saaty, 1980 as cited in Malczewsk, 1999
In the process of AHP, the prime task of calculation is the eigenvector corresponding to the
largest eigen value of the matrix. Each element in the eigenvector indicates the relative priority
61 | P a g e
of corresponding factor, i.e. if a factor is preferred to another; its eigenvector component is larger
than that of the other. A sum/product method is used to obtain the eigen value and the subsequent
eigenvector. The weights finally derived by AHP are used for developing the HIS model. To
examine the rationality of AHP, it is necessary to determine the degree of consistency that has
been used in developing the judgments. In AHP, an index of consistency, known as the
consistency ratio (CR), is used to indicate the probability that the matrix judgments were
randomly generated.
CR = CI
RI
Where RI is the average of the resulting consistency index depending on the order of the matrix
given by Saaty, 1980 and consistency index (CI) is defined as:
Where:
is the principal Eigen value of the matrix, n is the order of the matrix
The AHP also provides measures to determine inconsistency of judgments mathematically. The
CR, which is a comparison between Consistency Index (CI) and Random Consistency Index
(RI), can be calculated using the following formula: CR= CI/RI
It is recommended that the consistency ratio presents values below 0.1. CR* was also calculated
and found to be 0.06, which is acceptable to be used in the site suitability analysis. The computed
Eigen vector is used as a coefficient for the respective factor maps to be combined in weighted
overlay in ArcGIS Map.
62 | P a g e
Table 20: The Criteria Used in Industrial Site Selection
Soil
LULC geology type residential road lake swamp river slope
LULC 1 3 5 2 3 3 3 4 5
geology 1/3 1 2 5 5 5 3 5 3
Soil
type 1/5 1/2 1 5 3 5 5 5 3
reside 1/2 1/5 1/5 1 2 2 2 3 5
road 1/3 1/5 1/3 1/2 1 3 3 2 3
lake 1/3 1/5 1/5 1/2 1/3 1 1 2 5
swamp 1/3 1/3 1/5 1/2 1/3 1 1 2 5
river 1/4 1/5 1/5 1/3 1/2 1/2 1/2 1 3
slope 1/5 1/3 1/3 1/5 1/3 1/5 1/5 1/3 1
total 3.4 5.9 9.4 14.9 15.4 19.8 18.7 24.3 30.3
Consistency ratio = 0.06
Consistency is acceptable.
Source: Compiled by the Researcher, 2019
63 | P a g e
4.9.Calculation of the Weight Normalized
The normalization of criteria values allows for the direct comparison of criteria when calculating
suitability and reflects the conceptual relationship between the criteria value and a suitability
score (Jiang & Eastman, 2000). Priority vector is also called normalized principal Eigen vector.
To normalize the values, divided the cell value by its column total and calculated the priority
vector or weight by determining the mean value of the rows (Saaty, 1980). In order to check
whether the set criteria is correct or not, the researcher computed normalization.
Influenc
geology
swamp
weight
LULC
e (%)
value
slope
river
road
lake
LULC 0.3 0.51 0.5 0.12 0.2 0.15 0.15 0.16 0.15 2.24 0.25 25
geology 0.09 0.18 0.2 0.33 0.3 0.24 0.15 0.2 0.09 1.78 0.2 20
Soil
type 0.06 0.09 0.11 0.33 0.2 0.24 0.25 0.2 0.09 1.57 0.17 17
reside 0.15 0.03 0.02 0.07 0.14 0.1 0.1 0.12 0.15 0.88 0.1 10
road 0.09 0.03 0.04 0.04 0.07 0.14 0.25 0.09 0.09 0.84 0.09 9
lake 0.09 0.03 0.03 0.04 0.02 0.05 0.04 0.09 0.16 0.55 0.06 6
swamp 0.09 0.05 0.03 0.04 0.02 0.05 0.03 0.09 0.15 0.55 0.06 6
river 0.07 0.03 0.03 0.02 0.03 0.02 0.02 0.04 0.09 0.35 0.04 4
slope 0.06 0.05 0.04 0.01 0.02 0.01 0.01 0.01 0.03 0.24 0.03 3
total 1 1 1 1 1 1 1 1 1 9 1 100
Source: Compiled by the Researcher, 2019
64 | P a g e
Table 22: Classified Factor Values and Weights By Pair-Wise Comparison Method
65 | P a g e
Factor Class Value Level of Suitability Influence (%)
Lake 0 – 200 1 Unsuitable 6
200 – 500 2 Less Suitable
500 – 1000 3 Moderately Suitable
1000 - 10000 4 Suitable
Swamp 0 – 250 1 Unsuitable 6
250 – 500 2 Less Suitable
500 – 1000 3 Moderately Suitable
1000 - 2000 4 Suitable
>2000 5 Most Suitable
River 0 - 250 1 Unsuitable 4
250 – 500 2 Less Suitable
500 – 1000 3 Moderately Suitable
1000 - 2000 4 Suitable
>2000 5 Most Suitable
Slope 0-6 5 Most Suitable 3
6 - 10 4 Suitable
10 - 15 3 Moderately Suitable
15 - 20 2 Less Suitable
>20 1 Unsuitable
Overall Parameters Influence 100
Source: Compiled by the Researcher, 2019
66 | P a g e
4.11. Industrial site suitability analysis in the study area
Table 21 below shows that only 0.3 km2 (0.2%) of the study area are less suitable, 54.8km2
(34.7%) is moderately suitable, 74.9 km2 (47.5%) is suitable, 4.9km2 (3.1%) is most suitable and
finally 18.5 km2 (11.7%) is restricted areas.
Table 23: Statistical Analysis for Overall Weights for the Industrial Sites Selection
Suitability class Area (km2 ) Area (%)
Restricted area 18.5 11.7
Less suitable 0.3 0.2
Moderately Suitable 54.8 34.7
Suitable 74.9 47.5
Most Suitable 4.9 3.2
Total 157.7 100
Source: Extracted from the Industrial Site Suitability Map, 2019
As can be seen from the suitability map (Fig 30), about half (50.7%) of the study area is suitable
for industrial location. The suitable areas are mainly located in the Southern, South Eastern,
South Western and central part approaching to the corner of the town. These areas were mainly
vacant, bare lands and away from water bodies. In the Eastern part of the study area, on the way
to Wondogenet, the new structural plan incorporated vacant lands to the town which favors good
condition for the industrial site selection; and areas are free from urban settlements; and are
occupied by sparse vegetation. For this reason, those parts are suitable for industrial
implantation.
On the other hand, most parts of Huwella Tula sub-city were moderately suitable for industrial
sites. Even if the areas were sparsely populated and free from water bodies, because of optimum
transport inaccessibility up to the fringe the sub-city; it was a little bit important for the industrial
sites. And also, those areas with shrub trees and bare lands are suitable to plant industry for the
future. Obviously, those parts in the stated sub-city, the existing weather roads were considered
as hindrance for the transportation of raw materials. This would impose private investors not to
invest on that transport inaccessible areas. So in this area, small proportions of the study area
were leveled as moderately suitable.
67 | P a g e
Swamps, lakes, streams and forest coverage through Lake Hawassa up to Tikur Wuha area were
considered as restricted. For this reason, the researcher masked out those areas and leveled as
restricted parts. Therefore, about 11.7% of the total study areas were not permitted to suitability
criterion on this study.
In the study area, about 35% of the study areas were less suitable and moderately suitable for
industrial site selection. This is because of, central parts were occupied by the recreational
centers, settlements and built up areas, business sites and around Hayk Dar areas were
considered to be constraints for industrial establishment.
In the study area, industries were administered by the town Industry and Enterprise Development
Department, except Hawassa Industrial Park which administered by Federal Government of
Ethiopia (Industrial Park Development Corporation). As the information obtained from the
departments manufacturing and industrial controlling expert, the office has the following
68 | P a g e
stakeholders in executing its objectives, vision and mission. Those stakeholders were:
Construction department, the town Municipality Office, Environmental Protection and Forest
Conservation Authority and Housing and Urban Development Department.
In order to have erect industry, the industry department submits the desired action plan for those
stakeholders in order to implements the project properly. By doing this, the housing and Urban
Development and Municipality offices provides a vacant land for the proposed project,
Environmental Protection and Forest Conservation Authority undertakes environmental impact
assessment cases and construction department issues and evaluates construction process.
As the information obtained from urban and housing development department, vacant land
identified for industries based on the following factors of importance: transport facilities,
availability of construction land, and availability of labor costs, references from local partners
and previous experience, accommodation, projects public importance’s, compensation paying
capacity, noisy issues, capital and availability of local inputs.
Besides the above, the expert raised about the intentions of private sectors to invest in the town.
As he told, the private sector tends to invest more in the service sector shying away from the
manufacturing sector. The agglomeration of industries at central part of the town is also
considered as one of the problem in the study area.
During site selection, identifying vacant land for different purposes should be based on documented
parameters and criteria. And now a days, different software’s were employed to identify land for the
respective purposes like Land Information System (LIS), Geographic Information System ( GIS), Remote
Sensing techniques, AutoCAD, etc. but in the study area, land identified for the desired purposes merely
by the experts only by AutoCAD.
69 | P a g e
CHAPTER FIVE
GIS and Remote Sensing as analysis tools are valuable tools that can support the decision makers
to find best possible industrial sites. The GIS analysis requires collecting data from different
sources with different formats to create a complete uniform database. Thus, the GIS data should
be updated regularly in order to reflect the current situation of an area under investigation.
Remote sensing data can assist to have updated information of the study area. Also, it can
support the decision makers to monitor the investigated area using different dates of satellite
images to extract the urban land use/cover class for example. The findings have shown the ability
of GIS and remote sensing as a genuine tool for analyzing the criteria for decision support.
From those land use land cover factors, water bodies and forest coverage were masked out since
these factors were incorporated in suitability level. Vacant/open lands and shrub land covers high
percentage from the total area of the study and also they are suitable for industrial location. On
the other hand, settlements, business sites and recreational center are not suitable for industrial
location since interrupts them.
The analysis has taken land use/cover, slope, surface water, proximity to main roads and
streams/river, swamps, distance from residential areas, distance from geological faults and soil
types. From those parameters, land use land cover, geological faults and soil types were highly
70 | P a g e
determining the selection of industrial sites in the study area. Since the study area is relatively
flat slope, elevation and slope parameters play small portion in determining site selection.
After weight over lay computed, about 50% of the study area is suitable for industrial location
and only 0.2 % is less suitable to locate industry.
The application and utilization of GIS and Remote Sensing in identifying sites for different
purposes were very limited. Currently, there are a lot of technological instruments that are
important to administer land and land related properties. But in the study area, the town Urban
Development and Housing Department utilizes AutoCAD highly to administer land. Besides
this, there is also shortage of experts who are GIS and Remote Sensing backgrounds. This
resulted in ignorance of appreciating those best technologies like GIS and Remote Sensing.
GIS model is limited to the available data, which in this study; nine parameters were considered.
Therefore, any additional information such as wind direction, land price, raw material, security
issues, and government loan for the private investors and other social and economic factors can
enhance the outputs of the GIS model, and provide more realistic results. The planners and the
decision makers can get useful information about the possible locations of industry sites using
this methodology. Specially, the site ranking process allows for easily readjustment of the criteria
weights in case a sensitivity analysis is required. Nevertheless, defining detailed and standard
criteria by the Environmental Agency that comply with the local conditions can enhance the
outputs of GIS models used for the purpose of finding a suitable industry site. However, getting
public agreement on any candidate site is a must, and cannot be avoided. Therefore, the local
community should participate in the selection process of industry site to avoid any opposition in
the future.
5.2. Recommendations
Based on the findings, the following recommendations are forwarded by the researcher:
The GIS based multi criteria evaluation technique is simple and flexible which can be
used to analyze the potential sites for urban development and encourage public
participation in the urban decision making process. Thus, planners and authorities in
order to formulate suitable plan for sustained development of the town, they have to
undertake the application of GIS and Remote Sensing technologies.
71 | P a g e
In the study area, industries were accumulated in the central and main parts of the town.
This may result in different social and economic problems (interruption) on the residents.
Therefore, the town administration has focus in order to plant industries far away from
residential and other economic areas of the town.
Data handling system in the industry and enterprise department and Urban Development
and Housing Department was very poor. Therefore, they have to give recognition for
those currently accepted technologies like land information system, GIS and Remote
Sensing to address information flow for the public.
The study area selects industrial sites if there was a vacant land. Existence of a vacant
land itself is not enough to establish industry. Those listed factors in this study and
additional parameters have to be taken in to account for further industrial development.
This study has been intended to serve for solving the location problems associated with
industries. For private investors and enterprises, which are expected to participate in the
country in the near future, some additional factors will probably be required to be
considered, depending on the nature and complexity of industries.
The town Industry and Enterprise, Housing and Urban Development, Forest and
Environmental Protection Departments have to provide short term and long term GIS and
remote sensing training for the experts.
72 | P a g e
REFERENCES
Aleksander Rikalovic, Dordo Cosic, Srdan Popov and Dorde Lazarevic. 2013. Spatial Multi-
Criteria Decision Analysis for Industrial Site Selection: The state of the art, XI Balkan
Conference on Operational Research-Balcor, Belgrade.
Aleksandar Rikalovic, Ilija Cosic, Ruggero Donida Labati, Vincenzo Piuri. 2016. A
Comprehensive Method for Industrial Site Selection: The Macro-Location Analysis.
Alemayehu Abiye, T. 2008. Environmental resources and recent impacts in the Awassa
collapsed caldera, Main Ethiopian Rift. Quaternary International, 189, 152-162.
Amita Johar. 2013. Land suitability analysis for industrial development using GIS, Indian
Institute of Technology Roorkee.
Anurag Ohri, Prabhat Kumar Sing and Priyanka Kumari Singh. 2010. Spatial Multi Criteria
Analysis for Siting Industries, Banaras Hindu University, Varanasi.
Ayele Abebe. 2017. Land Use and Land Cover Change Detection Analysis using Remote
Sensing Techniques : The Case of Hawassa Town, Southern Ethiopia.
Ayenew, Tassew. 2004. Environmental implications of changes in the levels of lakes in the
Ethiopian Rift since 1970. Regional Environmental Change, 4, 192-204.
BASALFEW Z. 2012. Geological map 1: 250 000, sheet Hoseina with explanatory Note.
Geological Survey of Ethiopia.
Bhatnagar, R and Sohal, S. 2005. Supply Chain Competitiveness: Measuring the Impact of
Location Factors, Uncertainty and Manufacturing Practices, Technovation, 25(5), pp. 4434.
Church, R. L. 2002. Geographical Information Systems and Location Science. Computers &
Operations Research, 29(6), Pp. 541–562.
i|Page
Claus Renner. 2013. Multi-Criteria Decision Analysis for ArcMap. An Open-Source Multi-
Criteria Decision Analysis and Geovisualization Tool for ArcGIS 10,” Feature Article,
Cartouche, Newsletter of the Canadian Cartographic Association,
Cunningham. 2009. More than just the Kappa coefficient: A program to fully characterize inter-
rater reliability between two raters.SAS Global Forum, paper 242.University of Pittsburgh,
Pittsburgh, PA.
Czogala and J. Leski. 2000. Fuzzy and Neuro-Fuzzy Intelligent Systems, ser. Studies in
Fuzziness and Soft Computing. Springer, vol. 47.
Davis. 2001, Applications of GIS in Community Based Forest Management in Australia and
Nepal.
Dereje Tenaw,Enyew Abera, Xiaodi Zhang, Zhen Wan, Jie Zhao, Ciyong Zou and Jaidev
Dhavle. 2018. Industrial Park Development in Ethiopia Case Study Report, United Nations
Industrial Development Organization, Vienna.
Dessie, G. And Kleman, J. 2007. Pattern and Magnitude of Deforestation in the South Central
Rift Valley Region of Ethiopia. Mountain Research and Development, 27, 162–168.
Drobne, Samo, Lisec and Anka. 2009. Multi-attribute Decision Analysis in GIS: Weighted
Linear Combination and Ordered Weighted Averaging, Informatica 33, pp. 459–474.
Dragutin Zelenovic. 2003. Location of production systems, in the Design of Production Systems,
2nd ed. Novi Sad: Faculty of Technical Sciences, pp. 373-394.
Eastman, J.R., P.A.K. Kyem., J. Toledano, and Willian. 1993. GIS and Decision Making, Vol. 4,
Explorations in Geographic Information System Technology, UNITAR, Geneva.
Ebrahim Fataei, Mehdi Aalipour erdi, Habib Farhadi and Ali Mohammadian. (2015). Industrial
State Site Selection Using MCDM Method and GIS in Germi, Ardabil, Iran Journal of
Industrial and Intelligent Information Vol. 3, No. 4, December.
Eldin and Sui. 2013. A COM-Based Spatial Decision Support System for Industrial Site
Selection, Geography Department Texas A&M University College Station, TX 77843,
USA.
ii | P a g e
Eldin, K. Eldrandaly, M. Shouman and G. Na wara4. 2013. An Industrial Site Selection system
Use of COM in Integrating an Intelligent GIS-AHP Application. The International Arab
Journal of Information Technology, Vol. 10, No. 3.
FAO. 2006. World reference base for soil resources. A framework for international
classification, correlation and communication. Rome: Italy
Fernando Garcia-Menendez, F., E. Monier, and N. E. Selin. 2015. The Role of Natural
Variability in Projections of Climate Change Impacts on U.S. Ozone Pollution,
Geophysical Research Letters, 44, Doi:10.1002/2016GL071565
Florida, R. 2002. The Rise of the Creative Class, the Washington Monthly, pp. 15-25.
Francois Perroux. 1955. "Economic space: Theory and Applications" Quarterly Journal of
Economics, vol.64, No.1.
Greenhut. 1956. Plant Location in Theory and in Practice, University of North Carolina Press
Chapel Hill, N.C., p. 17.
Grim Haw D J. 1996. Bringing Geographical Information Systems into Business (2nd edition).
Cambridge, Geo-information International.
Hawassa Town Finance and Economic Development Department. 2007. Socio-Economic Profile
of Hawassa Town, Published, Hawassa, Ethiopia.
Hawassa Town Sanitation, Beautification and Parks Development Department (SBPDD) Report,
2012. Annual Planned and Executed Activities Report on Solid Waste Management
Practice in Hawassa Town, Ethiopia. un published material.
Hawassa Town Trade and Industry Development Department. 2017. Profile of Major Industrial
Distribution, Job creation and Industrial Distribution, unpublished, Hawassa, Ethiopia.
Internet.(2019).(http://www.geology.cz/projekt681700/vystupy/Hawassa_subsheet_0738_C4.pdf
_2019:4:18:24 am).
iii | P a g e
JICA. 2012. The Study on Groundwater Resources Assessment in the Rift Valley Lakes Basin in
the Federal Democratic Republic of Ethiopia, Japan International Cooperation Agency
(JICA), Kokusai Kogyo Co., Ltd., Ministry of Water and Energy (MoWE), The Federal
Democratic Republic of Ethiopia. Final Report (Main Report).
John Wiley and Son. 1994. A GIS-Based Approach to Mapping Probabilities of River Bank
Erosion: Regulated River Tummel, Scotland.
Kenichi Ohno. 2017. Hawassa Industrial Park, Impression from a Japanese perspective Tokyo,
February 17, 2017
Kotler, P.J., Haider, D.H., Rein, I. 1993. Marketing places: Attracting investment, industry and
tourism to cities, states and nations, New York, itd: Free Press.
Malczewski J. 1999. GIS and Multi criteria Decision Analysis, John Wiley and Sons, Toronto.
Malczewski J. 2006. GIS Based Multi-Criteria Decision Analysis: A Survey of the Literature.
International Journal of Geographical Information Science, 20(7), pp. 703–726.
Mamo Mihretu and Gabriela Llobet. 2017. Looking Beyond the Horizon: A Case Study of
PVH’s Commitment in Ethiopia’s Hawassa Industrial Park, Hawassa, Ethiopia.
Manish and Vivekanda. 2013. Identification of Potential Sites for Urban Development Using
GIS Based Multi Criteria Evaluation Technique. A Case Study of Shimla Municipal Area,
Shimla District, Himachal Pradesh, India.
Mark Berube. 2014. A GIS Multi-criteria Evaluation for Identifying Priority Industrial Land in
Five Connecticut Cities, Master's of Regional Planning University of Massachusetts,
Amherst.
Marzieh Reisi1, Lu Aye1 and Alireza Soffianian. 2011. Industrial Site Selection by GIS in
Isfahan, Iran, Department of Natural Resources, Isfahan University of Technology.
Masood A. Badri. 2007. Dimensions of Industrial Location Factors: Review and Exploration,
Journal of Business and public affairs Volume 1, Issue 2.
iv | P a g e
McCoy, K. Johnston. 2001. Using ArcGIS Spatial Analys, Environmental Systems Research
Institute, Redland.
McRaith and William. 2015. African Growth and Opportunity Act. Written Testimony to the
Senate Committee on Foreign Relations, Africa and Global Health Policy Subcommittee.
National Meteorological Agency of Ethiopia. 2019. Precipitation and Temperature Data of the
Selected Station from 2012 up to 2018, Addis Ababa, Ethiopia.
Pontus Olfson, Martin Herold and Foody and Stephen. 2013. Good Practices for Assessing
Accuracy and Estimating Area of Land Change Boston, USA.
Reisi, L. Aye, A. 2011. Industrial site selection by GIS in Isfahan, 19th Iran International
Conference on Geo-informatics.
Rizaludin Mahumud. 2015. Multiple-Criteria Analysis to Perform the Tasks of Site Selection.
Study and Analysis of Chennai Flood 2015 Using GIS and Multi-Criteria Technique, Civil
Engineering, Excel Engineering College, Namakkal, India.
Saaty. 1997. A Scaling Method for Priorities in Hierarchical Structures. J. Math. Psychol. 15.
Saaty, T.L. 1980. A scaling method for priorities in hierarchical structures: Journal of
Mathematical Psychology, Vol, 15, No, 3,Pp. 234 – 281.
Sarath.M, Sameer Saran and K.V.Ramana. 2018. Site Suitability Analysis for Industries Using
GIS and Multi Criteria Decision Making, ISPRS Annals of the Photogrammetry, Remote
Sensing and Spatial Information Sciences, Volume Iv-5, 2018 Isprs Tc V Mid-Term
Symposium “Geospatial Technology Pixel To People”, 20–23 November 2018, Dehradun,
India.
Thanuja, Leo. P, Pichaya.R, and Mohamed. 2015. Suitable Site Selection for Establishing SWM
Centers by Considering All Environmental, Social and Economic Criteria
Voogd H. 1983. Multi-criteria evaluation for urban and regional planning. London.
v|Page
Williams, E.A., and Massa, A.K. 1983. Siting of Major Facilities: A Practical Approach. New
York, USA: McGraw-Hill Inc.
World Bank. 2013. Ethiopian Economic Update II: Laying the Foundation for Achieving Middle
Income Status. Washington, DC: World Bank Group.
World Bank. 2010. World Development Indicators. Washington, DC: World Bank, available at:
<http://databank.worldbank.org/ddp/home.do>, accessed 10 July 2013.
Yemane Gebreziabher. 2004. Assessment of Water Balance at Lake Hawassa Catchment Area,
International Institute for Geo-Information Science and Earth Observation, Enschede,
Netherland, published.
Zahra, Aarushi Singh, Sanjay Kumar and Malik. 2014. Major MCDM Techniques and their
application-A Review. M. Tech. Student, Hindu College of Engineering, Sonepat Dept of
CSE, Hindu College of Engineering.
vi | P a g e
Appendix 1:
Questions Used During Discussions with Experts
vii | P a g e
Appendix_2: GPS Reading Taken in Order to Locate Industrial Distribution in the
Study Area
Name of Industries Longitude Latitude
Hawassa Industrial Park 38.4947 7.0677
Tabor Ceramic Factory 38.491 7.0795
National Tobacco Monopoly 38.502 7.016
Textile Factoey 38.501 7.034
Markan Soap Factory 38.503 7.032
Admasu Demeke Garage 38.48 7.045
Hawassa Fish Market 38.46 7.04
Hawassa Gloval Automotive 38.47 7.057
Nyala Motors 38.477 7.066
Monaco Automotive 38.479 7.069
A/T/Steal Factory 38.491 7.071
Sky Printing Press 38.492 7.056
Reis Engineering 38.487 7.081
viii | P a g e
Appendix _3: GPS Readings Used for Land Use Land Cover Map Accuracy Assessment
Longitude latitude LULC class Longitude latitude LULC class
38.483422° 7.084835° Forest 38.474160° 7.041648° Recreation site
38.482354° 7.083824° Forest 38.459349° 7.046367° Recreation site
38.479899° 7.071127° Forest 38.471404° 7.028855° Recreation site
38.470564° 7.057837° Forest 38.476784° 7.041539° Recreation site
38.470127° 7.057040° Forest 38.495594° 7.040151° Recreation site
38.469960° 7.055591° Forest 38.464736° 7.008926° Recreation site
38.466613° 7.051973° Forest 38.462797° 7.006421° Recreation site
38.457730° 7.044649° Forest 38.458121° 6.999136° Recreation site
38.460304° 7.048092° Forest 38.451511° 6.993929° Recreation site
38.540977° 7.001802° Forest 38.458402° 6.999233° Recreation site
38.542740° 6.999905° Forest 38.461641° 7.040042° Recreation site
38.529693° 7.002961° Forest 38.464040° 7.035738° Recreation site
38.527465° 6.995638° Forest 38.478600° 7.077299° Recreation site
38.529051° 6.990151° Forest 38.475286° 7.064631° Recreation site
38.538778° 6.990472° Forest 38.496536° 7.064692° Recreation site
38.526979° 6.973656° Forest 38.542639° 7.048126° Recreation site
38.483979° 7.094545° Forest 38.544945° 7.043989° Recreation site
38.485228° 7.098134° Forest 38.501390° 7.072204° Business site
38.487147° 7.097087° Forest 38.500790° 7.072178° Business site
38.487804° 7.097241° Forest 38.500297° 7.067948° Business site
38.487566° 7.097319° Forest 38.498353° 7.064447° Business site
38.484575° 7.094759° Forest 38.495175° 7.065935° Business site
38.483628° 7.093887° Forest 38.494774° 7.068068° Business site
38.482951° 7.089986° Forest 38.492494° 7.070574° Business site
38.482890° 7.090377° Forest 38.493114° 7.064889° Business site
38.484911° 7.087895° Forest 38.491921° 7.080075° Business site
ix | P a g e
38.489440° 7.083738° Business site 38.523574° 6.988737° Vacant land
38.501267° 7.045141° Business site 38.498822° 6.954458° Vacant land
38.502040° 7.029669° Business site 38.502947° 6.965053° Vacant land
38.484415° 6.961632° Settlement 38.509205° 6.964049° Vacant land
38.489884° 6.962129° Settlement 38.522706° 6.970680° Vacant land
38.479214° 6.960006° Settlement 38.497240° 7.075006° Vacant land
38.476758° 6.958158° Settlement 38.506624° 7.057388° Vacant land
38.482228° 6.954142° Settlement 38.506559° 7.056211° Vacant land
38.491119° 7.039858° Settlement 38.515764° 6.922166° Vacant land
38.482034° 7.040622° Settlement 38.520660° 6.923297° Vacant land
38.470807° 7.048367° Settlement 38.476779° 7.075826° Water body
38.486432° 7.060467° Settlement 38.463552° 7.020793° Water body
38.490518° 7.090974° Settlement 38.464995° 7.013753° Water body
38.498486° 7.087838° Settlement 38.449940° 6.990588° Water body
38.496174° 7.082727° Settlement 38.446600° 6.986550° Water body
38.486918° 7.076236° Settlement 38.515113° 6.922768° Water body
38.496522° 7.078295° Settlement 38.512536° 6.927671° Water body
38.513905° 7.046099° Settlement 38.508389° 6.959080° shrub land
38.514648° 7.042855° Settlement 38.501551° 6.934354° shrub land
38.511539° 7.035725° Settlement 38.527652° 6.995193° shrub land
38.491288° 7.015194° Settlement 38.540675° 6.922011° shrub land
38.526689° 7.011506° Vacant land 38.540372° 6.918625° shrub land
38.524142° 7.007826° Vacant land 38.528368° 6.934358° shrub land
38.517398° 7.007124° Vacant land 38.523181° 6.939710° shrub land
38.517317° 7.006146° Vacant land 38.516597° 6.944223° shrub land
38.515958° 7.003167° Vacant land 38.518440° 6.952380° shrub land
38.521660° 6.997780° Vacant land 38.540433° 6.997507° shrub land
x|Page
Appendix_4: APH Weight Derivation Procedure
xi | P a g e
Appendix 5:
Sample for Geological Fault Sites of Study Area
xii | P a g e
Appendix_ 6: Factors Suitability Map
A B C
D E F
G H I