Using a GIS to Assess the Land Movements Hazard: Application on Berhoum Area, Hodna Basin, Algeria

Modeling land movements hazard by multi-criteria approach is a line of research to provide a methodological framework for risk mapping. This work is intended to establish a model for mapping “land movements hazard” by GIS approach based on multi-criteria analysis. The methodology is to create thematic maps by combining in GIS, determining factors (slope, lithology, water, ...) in triggering landslides phenomena and shrinking-swelling soil, using the method of Weighted Sum Model (WSM). These maps are then combined to provide a hazard map of land movements. The application of this method allows the spatial distribution of different criteria and phenomena in Berhoum area, region of Hodna Basin, eastern Algeria. As a result, the study area has been divided into four different areas: 1) areas with no land movements hazard, 2) areas with a medium land movements hazard, 3) areas with a high land movements hazard, and 4) areas with very high land movements hazard.


Introduction
Natural risks threaten with more or less intensity the human activities.Algeria has been frequently touched by natural disasters; these last ones put in danger life of persons and would make inestimable damages.
The north regions of Algeria, characterized by population density and speed of urban growth, are witnessing difficulties in construction caused by problems related to the land movements.This is due to the weakness of knowledge of phenomenon like: landslides and shrinking-swelling soil.Land movements include all field cases of land masses displacement under the action of gravity.They can be very slow or very fast.
Land movements can be classified into two types: Slope Movements: these are oblique movements, which affect both natural and artificial slopes, collapses, slides, ect., Vertical movements which include subsidence, collapses, settlements, ect., which result from the progressive consolidation of sediments under the influence of their own weight … they can be caused either by local overload, the soil surface as that resulting of construction of a work, or by extraction of subsurface material, as in the exploitation of oil, water [1].
Based on this definition, the land movements are associated with two basic components: Landslides and shrinking-swelling soil.Landslides and shrinking-swelling soil constitute the largest classes of hazards.
The hazard assessment of land movements, including landslides and shrinking-swelling soil has been the subject of several research papers [2]- [11].
Several relationships and mathematical models hazard have been established, based on a single standard or several standards, and specialized in specific phenomenon, such as landslides or shrinking-swelling soil.
Is there any way to collect and evaluate these two phenomena together?
This can be achieved by integrating the Multicriteria Decision Analysis (MCDA) and the analytical capabilities of GIS [12].
There are several methods for solving and evaluating selection problems in MCDM [13] [14].The main technique closely related to this paper is Weighted Sum Model (WSM).WSM can be described as a weighted linear combination or scoring method [15] [16].WSM is a representation of the Decision Maker's preferences in a linear additive function [17].Triantaphyllou and Mann [18] examined the effectiveness of decision-making methods and pointed WSM as one of the simplest methods for solving MCDM problems [19].
In general, suppose that a given MCDA problem is defined on m alternatives and n decision criteria, furthermore, let us assume that all the criteria are benefit criteria, that is, the higher the values are, the better it is.Next suppose that w j denotes the relative weight of importance of the criterion C j and a ij is the performance value of alternative A i when it is evaluated in terms of criterion C j .
Then, the total (i.e., when all the criteria are considered simultaneously) importance of alternative A i , denoted as WSM score i A − , is defined as follows: For the maximization case, the best alternative is the one that yields the maximum total performance value [20].
From (1), the intuitive nature and ease of dealing with Multicriteria Decision-Making problems can be seen, and this simplicity makes WSM one of the most important methods for both exploring possible solutions for a Multicriteria Decision-Making problem and providing solutions in comparison to other methods [13].Journal of Geographic Information System This work aims to use a model of land movements hazard based on multicriteria analysis approach and assessment by GIS.This model, using (WSM) method in the framework of GIS is tested, and the obtained results have been discussed.The results provide useful information about landslides and shrinking-swelling risks mitigation and may serve as guidelines for land use planning in Berhoum area.
The study is part of the contribution to Berhoum area management, which is vulnerable to natural hazards, and also in the context of the production of information and awareness documents.So, the preventive methods remain the best way to overcome these problems or reduce them.One of the means that allows the adoption of effective preventive methods is to understand the phenomenon and predict the places it occurs in.
The choice of site of investigation is focused on Berhoum area, located in the South eastern of Algiers and is part of El Hodna Basin in north-east of Algeria (Figure 1).This choice is dictated by its diversity natural properties and density of equipment.

Methodology
After identification of the parameters involved in land movements, a mathematical model has been well established.
In order to obtain land movements hazard model of the study area the method of Weighted Sum Model (WSM) was applied (Figure 2).
In the analysis, weight values to the layers and rank values to the classes of each layer were assigned.For each class of the layer, rankings were given according to their significance.After the rankings were assigned to the classes of each layer, the weights were assigned to layers according to their importance.
The interaction between the layers was not taken into consideration since the layers were assumed to be independent of each other [15] [16].The weights values, and classification in (WSM) model, are based on opinions of experts (Sounding views was conducted with a group of experts in the framework of the preparation of the doctoral thesis of the author scientific article) [21].
The weight and rank values of the layers and classes of each layer were standardized in order to obtain a common dimensionless unit.Afterwards, the output was created by multiplying the weight value assigned to each layer by the rank value given to the classes of that layer and finally by adding up the products [22].
We apply these steps to get the landslide map, as well as the shrinking-swelling soil map.As a result, the land movements hazard map of Berhoum area was grouped into different zones.
The application of the proposed model, suggested uses the SIG approach.We examine it as follows:

Hypotheses
To achieve this model, assumptions must be made.In this context, two hypotheses arise: • The land movements are a phenomenon that can be evaluated by Landslides and shrinking-swelling soil, phenomena.
• Hazard is the weighted sum of these factors.

Modelling Land Movements Hazard
Land movements hazard is also the weighted sum (WSM) of the determinants of the phenomenon.Land movements hazard is given by the following relationship: where: w i : weight assigned to each determinant, C i : criteria.
According to the definition of land movement's hazard, it is the sum of landslide and shrinking-swelling soil.The combined effect of these two parameters is expressed by the following equation: 0.5 0.5 where: H LM : the land movements Hazard, H LS : the landslides Hazard, H SS : the shrinking-swelling soil Hazard.

Modelling of the Landslide Hazard
Based on the principle of the weighted sum, and in relation to the slope, water and soil, the landslide Hazard is given by the following formula: where: P: slope angle, TS: soil types, PNP: depth of ground water table, EP: rain.

Modelling of Shrinking-Swelling Soil Hazard
On the same basis, and in relation to the types of soil and water, the shrinking-swelling soil hazard is given by the following formula:  3).

Rain (EP):
The rain was divided into four categories (Table 4).Journal of Geographic Information System

Crossing between the Parameters and Classification
To make the crossing between the parameters, it is better to segment them into classes according to their degree of effect.The results are listed in Table 5 & Table 6.
The combination of the weighted results of the different criteria, one by one and taking the common factors in each time, gives the level of each phenomenon, which we have presented in the form of tables.
Intersection matrixes for landslide are given as follows (Table 7).
Intersection matrixes for shrinking-swelling soil are given as follows (Table 8).

Classification for Landslide and Shrinking-Swelling Soil Hazards
The classification for landslides hazard is represented in Table 9.
The classification for shrinking-swelling soil is represented in Table 10.

Modelling Land Movements Hazard
For a combination of sliding and shrinking-swelling, formula (3) is applied for modeling the map of land movement's hazard.The rank of weights and values are follows in Table 11.
Intersection matrixes are given as follows tables (see Table 12 and Table 13).
Table 1.Classifications of slope effect on the hazard (P).

Spatialization by GIS Approach: Application in the Berhoum Area
The approach we have adopted for the spatial distribution map preparation to assess the land movements Hazard is based on the geographic information systems.
The data sets used in this study can be grouped in four main data sources as lithological map (rocks' map), field survey and geotechnical drilling data, topographical map and rainfall map.These data obtained from the topographical map, rainfall map, lithology map or from the field survey and geotechnical drilling, have been implemented in GIS program (Arc Gis 9.3) for the primary maps drawing.
After the building of a database of the geographic reference, organizing and structuring it by the use of programs, a map of the land movements Hazard that shows the best areas for building.

Lithology Map
The lithology map is derived from the geological map of the study area [23].After numeration of the geological map, the rocky facets' database has building.For the study area, nine lithological units were identified (Figure 3).These units include Dolomite, Limestone, Sandstone, Conglomerate, Marls, Clays, Coarse Alluvium Sandy-Clay Alluvium and Fine Alluvium.In the study area, Dolomite, Limestone, Sandstone, and Conglomerate Formations were classified as rock.

Field Survey Data and Geotechnical Borehole
Indeed, three techniques-different but complementary, have been used to complete the database: • The relevant data of 49 study geotechnical surveys were reviewed and evaluated, • Well field investigation based on direct observations, • As well as recordings and field sampling (GPS, images, etc.).
Soil types (TS) and the groundwater table layers are prepared using lithology map and geotechnical borehole (Figure 4 and Figure 5).

The Topographical Map
The topographical map in scale 1/50,000 is used to obtain the Digital Elevation Model (DEM) of the study area.The latter was obtained by digitizing the contour lines of the topographical map.It is the software Arc Gis 9.3 which was used for the interpolation.The slopes map is prepared by degrees, using the DEM of the study area (Figure 6).

The Rainfall Map
The rainfall map in scale 1/50.000 is used to obtain the rain water layer of Berhoum area (Figure 7).

Data Evaluation
These data were consolidated and evaluated to obtain landslide and shrinking-swelling soil hazards maps (Figure 8 and Figure 9).To implement these steps, we use the machine tool calculator networking (Calculator Raster) of spatial analysis toolbar (Suitability Special) in Arc Gis 9.3 [24].

Identification and Classification of Land Movement's Hazard
Combining the results of landslides and shrinking-swelling soil (Formula 3), gives the degree of movements Hazard (H LM ), we have presented in the form of Figure 10.

Validation of Land Movements Hazard Map
Validation of the results was done at two levels.On the one hand, the thematic maps obtained were compared with the field data before their crossing.The perfect match between these field observations (Figure 11) and the medium and high hazard sites revealed by the map made from the GIS allows validating to some extent the final map.On the other hand, during field missions, the survey conducted (whether through technical offices of municipalities, by the Technical Inspection Agency of Constructions (CTC-EST, Agence de M'SILA) or directly in the field) did not record any shows slides in the study area, and therefore did not record any damage to buildings caused by landslides.Also we recorded the presence of some limited damage to buildings (Through Technical Inspection Agency of constructions reports), and are the same observations made by [25] in his research, in each of the municipalities Berhoum and Ouled Addi Gueballa, due to shrinking-swelling soil problems [26].

Discussion
The urban expansion and the increase in population led the urbanization to use the spaces "called at risks".The integration and the treatment of lands movements constitute an important characteristic of the equilibrium established by nature [11].
In order to solve these problems, it was necessary to propose a land movements hazard map, using the geographical information systems (GIS) with the adoption of (WSM) model.These tools enable us to express the perception of space and data processing, and consequently the cartography is carried out in an optimal way.These land movements hazard map have a great part of importance in all levels of a study as information, working paper, alert, and especially is a tool for the decision-making aid, by expressing tendencies and orientations.
The land movements hazard map of the study area (Figure 10) and a database with all the compiled geotechnical information (GEODATA) is a very important contribution to the land planning and management of the study area, and the definition of natural risk zones.It provides an excellent geotechnical data source to be consulted when new projects are being planned and designed for the study area.
An open-database, with information concerning boreholes, water table location, laboratory test data, sample locations and photographs of the samples, should be added to the Lithology, Soil types and the groundwater table maps.Such a database will permit a rapid evaluation of the conditions of the land movements hazard.
The data came from geotechnical reports provided by public institutions and also from our own fieldwork and investigations allowed the data collection to get the landslides and shrinking-swelling soil maps (Figure 8 and Figure 9).These data have been used to calculate the land movements hazard (Figure 10) by the relationship (3).
The proposed model is simple, and easily used by the geotechnical engineer or otherwise.Thread model is very interesting as it is being built on a small number of geotechnical characteristics according to experts' judgment.Similarly, this model can be easily applied to other similar regions in terms of the environmental landscape and climate.This model is not static however; it can be changed sometimes by giving a great attention to one factor and decrease it to another.

Conclusions
This work on cartography of land movements hazard, using geographic information systems (GIS) techniques, with a Weighted Somme Model (WSM), shows the interest caused by these methods, which sounds of the multiple forms, opens the way to a very sensitive improvement of the prevention of the natural risks.These tools enable us to express the perception of space and data processing, and consequently the cartography is carried out in an optimal way.These land movements hazard maps have a great part of importance in all levels of a study as information, working paper, alert, and especially are a tool for the decision-making aid, by expressing tendencies and orientations.

Figure 2 .
Figure 2. Organization chart of a Geographic Information System for Land movements hazard mapping.
water Journal of Geographic Information System Colors reflect the degree of land movements Hazard (H LM ).

Figure 5 .
Figure 5.The ground water table layer.

Figure 10 .
Figure 10.Land movements hazard map of Berhoum area.

Figure 11 .
Figure 11.Field photographs showing the shrinking-swelling soil problems in the study area.

2.3. Criteria Used to Shrinking-Swelling Soil and Landslide Hazards Slope
TS: soil types, PNP: depth of ground water table, EP: rain.2.angle (P): Most large landslides related to the slope angles whose value generally exceeds 45%.The slope angle values have been divided into four categories (Table 1).Soil types (TS): Return to LPC classification, fine soils [Little plastic clay (Ap), Very plastic clay (At), Little plastic silt (Lp) and Very plastic silt (Lt)], are considered the most sensitive to swelling-shrinkage soil phenomena.The granu-Depth of ground water table (PNP): The depth of groundwater table was divided into four categories (Table

Table 2 .
Soil types effect on the hazard (TS).

Table 3 .
Classifications of depths of the groundwater table effect on the hazard (PNP).

Table 4 .
Classifications of rain effect on the Hazard (EP).

Table 5 .
Weights and values of the criteria for landslide Hazard.

Table 6 .
Weights and values of the criteria for shrinking-swelling soil Hazard.

Table 11 .
Weights and values of the criteria.