Geological – Geomorphological Features of River Catchments in Flood Susceptibility Assessment (on the Example of Middle Struma Valley, Bulgaria)

. The geological-geomorphological features of the river catchments control flood propagation by influencing the surface runoff and its redistribution to surface and underground flow. The paper analyses the geological-geomorpho-logical susceptibility to floods in a part of the river Struma valley – Municipalities Kresna and Strumyani (Bulgaria, EU). The following parameters are considered: lithology (presented by water permeability, capacity and type of rock reservoirs), distance from streams/river, slopes and fluvial landforms. Raster layers for each one of these parameters are created in geographic information system (GIS) environment and susceptibility rates are assigned. The flood susceptibility map is elaborated by spatial weighted sum overlay of the raster layers. The results show that the susceptibility is the highest in the low part of the valley, south of Kresna. The susceptibility of flooding is not equal over the whole area of the floodplain but depends on the retention capacity of the alluvial deposits and the hydro-geo-morphological state of the floodplain. The river bank height above the water level and the river bed type are also taken into consideration by application of fuzzy logic. Using the fuzzy logic minimizes some imperfections of the initial data and gives good results in data-scarce areas. The fuzzy logic model is designed as a two-level hierarchical system with three inputs and one output. The output gives the complex assessment for geological-geomorphological flood susceptibility of study area regarding the interaction of the three inputs. The results of the current research aims to direct the attention of planning experts and decision-makers to the importance of geological-geomorphological parameters in flood management and can be considered as a first step of flood hazard assessment.


First Section
Floods are ones of the most important natural hazards in Europe in terms of economic loss [1,2] and when they concern human life this is a risk.According to the Directive 2007/60/EC [3] the flood risk means the combination of the probability of a flood event and of the potential adverse consequences for human health, the environment, cultural heritage and economic activity associated with a flood event.In this regard the risk is considered as a function of hazard (the probability of flood event).The flood hazard is determined by the hydro-climatological features of the area which behavior is controlled by the other environmental conditions (lithological composition, physical-mechanical properties of soils, character of the topography, vegetation) and also on the land use type.The review of the literature and practice in flood management shows different approaches.The traditional strong engineering approach is based on hydrological and hydraulic analyses and modelling [4][5][6].This approach is often limited by the availability of the data which requires the fusion of hydrologic observations, historical information about past flood events, estimates of surface roughness, boundary conditions, and topographical features [7].In this regard many publications consider the importance of sedimentary and geomorphological approaches in providing base-line flood data for validation hydrologic and hydraulic models and emphasize on the maintenance of natural geomorphological functioning of the river channels and floodplains for effective flood management [8][9][10][11].A combined application of geomorphological and hydrological approaches in flood hazard investigation is considered in publications Patton and Baker, Lastra, Dimitriu, etc. [10,12,13].
The complex nature of the flood event requires an integrated approach in flood hazard / flood risk assessment and management.An integrated approach employing remote sensing data coupled with digital elevation models (DEM), hydrological study and field observations is recommended to understand the causative factors of flooding [14].The need of processing a large volume of spatial data motivates the geoinformation approach and use of geographic information systems (GIS), and remote sensing in flood hazard assessment [14][15][16][17][18].
Taking into account the complex nature of flood events, the aim of the current research is to assess flood susceptibility of geological-geomorphological features in the Middle Struma valley by application of GIS technology and fuzzy logic.The susceptibility is considered as a sensitivity degree of the area and determining the flood susceptibility is the first step in flood risk assessment.Although the intensive precipitation have triggering role for flood propagation the differences in the distribution of floods in the river basin show the control of geomorphological and geological factors.The geological-geomorphological features of the area influence the flow of surface runoff and the redistribution of the runoff into the underground and surface component.This influences the water capacity of the area and possible occurrence of floods.The composition of the rocks, their grain size and physical properties determine the water permeability of the rocks as well as their retention capacity.In this relation, in flood susceptibility assessment the water permeability of the rocks should be considered in a complex with other environmental characteristics such as water retention capacity, groundwater level, location of the site relative to the water level / river bed, landforms, etc.
In this research an integrated assessment of these factors is elaborated in GIS environment and by fuzzy logic.The following geological-geomorphological parameters of the investigated area are evaluated: lithological composition, slopes, fluvial landforms and the distance from the rivers.A flood susceptibility map is created in GIS environment by application of weighted sum overlay.A similar approach combined with application of analytic hierarchy process is used many publications about different areas and shows good result, even for data-scarce areas [17,[19][20][21].In difference, the present study focuses the geological and geomorphological conditions of flooding by evaluating the water permeability, retention capacity and type of the rocks reservoirs.For the purpose of modelling a fuzzy logic is applied additionally to the results of the spatial GIS overlay.

Data and Methods of the Research
The flood susceptibility modelling in the current research is elaborated for a part of the Middle Struma valley (Bulgaria), the areas lower than 600 m.Special attention is given on the fluvial landformsfloodplains and low river terraces, as well as alluvial / proluvial fans.Field investigation and literature review, including analyses of geological (sheets Blagoevgrad, Razlog and Petrich, 1:100 000) and geomorphological maps [22,23] are done for data collection.The analysis of geological and geomorphological data is proceeded in AcrGIS environment (ESRI Inc.).For this purpose a GIS database was built which includes:  Municipalities boundariespolygon layer (source: GIS, Ministry of regional development, Bulgaria, EU);  Stream / river networkgenerated on the base of DEM;  Slopesgenerated on the base of DEM;  Geology data: rock permeability; capacity and type of the reservoirpolygons layers created on the base of the geological map of Bulgaria, 1: 100 000;  River terraces, alluvial / proluvial fanspolygons determined on the base of publications [23,24], topographic map of Bulgaria 1: 50 000 and field investigations;  Settlementspolygons (CORINE Land Cover 2012 project [25]);  Digital elevation model (DEM), SRTM, 30 m NASA JPL [26].
In order to assess flood susceptibility a spatial overlay analysis is done on the base of geological and geomorphological features of the area.The complex influence of the lithological composition, slopes, fluvial landforms and the distance from the rivers were evaluated in determining the susceptibility to floods.The evaluation is done using weighted sum overlay in ArcGIS.For this purpose, the geological-geomorphological factors are entered in raster layers, and depending on their characteristics, each one is classified from 1 to 5 (1 -the lowest flood susceptibility, 5 -the highest).In the overlay analysis a weight of 30% is given to each one of the layers of river distance and fluvial relief, 20% -for lithological composition and 20% for slopes.As a result of the overlay analysis a map of geological-geomorphological susceptibility to flooding is created.
Taking into account the limitation of the data we applied fuzzy logic in order to further develop the assessment.The following additional parameters are considered: river bank height above the water level and the river bed type.The first one is considering the aptitude of river water to flood the adjacent areas.About the river bed we define 3 types: sandy, rocky and mixed.Sandy bed has higher water permeability which is a prerequisite for better infiltration and less surface stream flow, while rocky bed determines less infiltration and higher flood susceptibility.On the other side because of the high water permeability of sands they are usually related with high level of ground water and this can lead to high flood susceptibility if the river bed is saturated.In this relation all available geological and geomorphological conditions of the area should be considered in interaction in flood susceptibility modelling.
It is necessary to note, that the fuzzy logic approach is an appropriate tool for the risk assessment from natural disasters [27][28][29].This approach provides adequate processing the expert knowledge and uncertain quantitative data [30,31].In particular, the flood risk assessment is performed under the subjective and uncertain conditions (for example, there are short records or only few measured points) [10,20,32].For these reasons the fuzzy logic are used for the flood risk assessment [21,33,34] In this research, a fuzzy logic model for complex assessment of flood susceptibility based on geologicalgeomorphological features of river catchments is proposed.The fuzzy logic model is designed as a two-level hierarchical system with three inputs and one output.Each level of the system is consisted of one fuzzy logical subsystem with two inputs.
The fuzzy logic system output gives the complex assessment for geological-geomorphological flood susceptibility of study area regarding the interaction of the three inputs.A scheme of this two level hierarchical fuzzy system is presented on Fig. 1.

Fig. 1. Two-level hierarchical fuzzy system with three inputs
The fuzzy system inputs are defined on the basis of expert knowledge and current research of study area, as follow:  Input 1 "Flood susceptibility according to the river bank height above the water level";  Input 2 "Flood susceptibility according to the river bed type";  Input 3 "Flood susceptibility according to the overlay analysis of geological-geomorphological parameters (map results)".
An intermediate variable "Flood susceptibility according to the river bank height and bed type" is defined.
Here, the inputs of the first fuzzy logic subsystem are Input 1 "Flood susceptibility according to the river bank height above the water level" and Input 2 "Flood susceptibility according to the river bed type", and the output variable is the Intermediate variable "Flood susceptibility according to the river bank height and bed type".
The inputs of the second fuzzy logic subsystem are Intermediate variable "Flood susceptibility according to the river bank height and bed type" and Input 3 "Flood susceptibility rate according to the overlay analysis of geological-geomorphological parameters (map results)".The output of second fuzzy logic subsystem is the variable "Complex geological-geomorphological flood susceptibility", which is actually the output of the proposed fuzzy logic system.
The value of the variable "Complex geological-geomorphological flood susceptibility" can be considered as a useful assessment in informative decision making about the flood susceptibility regarding geologicalgeomorphological features of river catchments in the study area.The higher value corresponds to the higher flood susceptibility rate.
In the proposed fuzzy logic model, the linguistic variables, corresponding of the three inputs, intermediate variable and system output are defined by different trapezoid membership functions.The input variables (Input 1, Input 2, Input 3 and Intermediate variable) are assessed in the interval [0, 5] and the system output -in the interval [0, 100].
The numeral values of trapezoidal membership functions and the corresponding linguistic variables of the Input 1, Input   IF "Flood susceptibility according to the river bank height and bed type" is "High" and "Flood susceptibility according to the overlay analysis of geological-geomorphological parameters (map results)" is "Very high" THEN "Complex geological-geomorphological flood susceptibility" is "High" IF "Flood susceptibility according to the river bank height and bed type" is "Very high" and "Flood susceptibility according to the overlay analysis of geological-geomorphological parameters (map results)" is "High" THEN "Complex geological-geomorphological flood susceptibility" is "Very high".
The logic matrix with inference rules of one fuzzy logic system with two linguistic input variable can be expressed also in the form of table.The logic matrices corresponding to the proposed two fuzzy logic subsystem are represented in Tables 2 and 3. Here, the fuzzy logic hierarchical system is designed in MATLAB computer environment using Fuzzy Logic Toolbox [35].The two fuzzy logic subsystem are based on Mamdani's inference machines, max/min operations and center of gravity defuzzification [36].
The inference surfaces in 3D for the two fuzzy logic subsystems are shown on Figures 5 and 6, respectively as follow: (Input 1, Input 2, Intermediate variable) and (Intermediate variable, Input 3, Fuzzy logic system output).The rock varieties in the study area can be divided according to their water retention properties, the water permeability and the distribution of surface and underground runoff.Porous, cracked, cavernous and mixed type (porous and cracked) reservoirs are distinguished in the area.Porous reservoir properties have quaternary alluvial and deluvial deposits, and Neogene and Paleogene clastic sediments.The weathering crust, which is a few meters thick on the granite and the metamorphic rocks can be presented as a porous reservoir, as well as the deluvial deposits.Granite and the metamorphic rocks have a cracked reservoir.Cavernous reservoirs are found only in highly soluble rocks and most often they are carbonates -limestone, dolomite or marble.
The landforms, which are most vulnerable to floods, and on the other hand, impact on the floods propagation are floodplain.Two floodplains (low and high) are developed in the river Struma valley, south of the Kresna gorge.There are also dynamic river deposits -islands that are 1 to 2 m above the average water level.Due to the active hydrodynamics of the river, the height of the floodplains is not constant.Often the connectivity of the low floodplain is cut and this increase the flood susceptibility and hazard.

4
Flood Susceptibility Assessment

Flood Susceptibility According to the Lithology
The lithological composition is assessed as a complex indicator, calculated as a weighted sum of the susceptibility rates of water permeability, capacity and type of rock reservoirs.In this sum 40% weight is given to the indicator type of the reservoir, 35% -to the water permeability and 25% -to the reservoir capacity.For this purpose a GIS polygon layer is created on the base of geological map [22] and field investigations.The different rocks are rated from 1 (very low) to 5 (very high) according to their importance for surface runoff retention and the probability for high water (Table 4).The flood susceptibility rates are determined by expert assessment and taking into account the specific geological-geomorphological conditions of the environment.For example, sands and gravels have the highest water permeability, they can collect a significant part of slope runoff, convert it into an underground flow, and thus minimize the risk of flooding, but already saturated, these can help to form larger surface runoff.Regarding that sands and gravels are usually related to high level of ground water they have to be rated with higher flood susceptibility in the complex assessment.About the reservoir type, the karst and pore reservoirs are rated with lower flood susceptibility considering their surface water-retaining properties relative to surface runoff, and cracked reservoirs have a higher coefficient because they determine a greater surface runoff.
Based on the susceptibility rates for the above indicators the resultant indicator of flood susceptibility is determined by weighted sum overlay in ArcGIS environment (Fig. 8).

Flood susceptibility according to the slopes
The slopes of the topographic surface affect the flow rate of slope runoff and the water retention, and in this relation the possibility of flooding.The slopes are calculated by ArcGIS (Spatial Analyst Tools > Surface) based on the digital elevation model (SRTM) [26] and are rated to susceptibility as follow (5very high, 1very low), (Table 5).
The suggested rates concern river floods and are not applicable for torrential floods or debris floods.In interpreting and assessing the slopes, ridges and watersheds lines are excluded.The spatial distribution of the slope flood susceptibility is presented on Figure 9.The areas steeper than 20⁰ are also excluded from the assessment.Very high and high flood susceptibility is observed south of the town Kresna in the low parts of the river Struma valley.This is determined by small slopes of the topographic surface which is a condition for retention of the surface runoff and also higher level of ground water contribute to high flood susceptibility.

Flood Susceptibility According to the Fluvial Relief
Floodplains and fluvial fans (alluvial and proluvial) are evaluated for flood susceptibility of fluvial relief.Having regards the closest relation of floodplains and fans to hydrological processes, floodplains are rated as most susceptible to floods (5) and fluvial fans -as highly susceptible (4) (Figure 10).

Flood Susceptibility According to the Distance from Streams
The distance from rivers and streams is assessed in relation to the possibility of high water to flood adjacent areas.Five distance zones are generated by ArcGIS Analysis Tools -Proximity -Multiple Ring Buffer.They are rated according to the possibility to be flooded from 1 (very low) to 5 (very high).For the aim of spatial overlay analysis the polygon layer is converted in raster format.The result of the assessment of flood susceptibility by distance from streams is presented on Figure 11.Areas higher than 600 m are excluded of the assessment.

Geological-geomorphological Flood Susceptibility
The assessment of the geological-geomorphological flood susceptibility is done by the weighted sum overlay in ArcGIS environment of the 4 generated raster layers, presenting the flood susceptibility rates for lithology, slopes, fluvial relief and distance from the streams.In order to determine the weighting rates of these factors, an analytical hierarchy process is used.As a result of the analysis, the following weight coefficients were applied:  distance from streams -45%;  fluvial relief -35%;  lithology -13% and  slopes -7%.
The output raster of the weighted sum overlay is shown in Figure 12.Taking into account the environmental characteristics of the investigated area, the analysis excludes areas with mountainous relief, higher than 600 m above sea level.If it considers the floodplain as a most susceptible to flooding, the results of the overlay analysis show that not always the whole floodplain is very high susceptible.The susceptibility rates depends on the other environmental parameters, for instance distance from river and retention capacity of the alluvial layer.

Simulation results from fuzzy logic model
Regarding to the above results, additional 2 parameters are considered for more detailed assessment of geological-geomorphological flood susceptibility.These are "River bank height above the water level, [m]"and "River bed type".Each one of them is rated according the possibility for flood propagation on the base of literature review and expert assessment (Table 1 and Table 6).Taking into account the discontinuous nature of available data about river bank height and river bed type, a fuzzy logic model is suggested for determining the complex geological-geomorphological flood susceptibility.
Here, several simulations with proposed fuzzy logic model were done using particular input values in order to calculate the rate of the "Complex geological-geomorphological flood susceptibility" (Fuzzy logic system output) of the study area.Some simulation results are show in Table 7.
The obtained output values (crisp values, established numerical value) from the fuzzy logic model after the defuzyfication belong to one of the five sub-intervals of the interval [0, 100].
These five sub-intervals are used to determine the rate of the "Complex geological-geomorphological flood susceptibility" of study the region as shown in  The simulation results show the importance of the river bank height in flood susceptibility, following the observed cases where floodplains are not equally inundated during flood events (Table 7).
For example if the geological-geomorphological flood susceptibility is Very low to Low (Input 3 = 0.8), the river bed is partly sandy and partly rocky (Moderate flood susceptibility, Input 2 = 2.5) and the height of the river bank is 3 m (Input 1 = 1.5, Low flood susceptibility), then the Complex geological-geomorphological flood susceptibility is Low (Output = 30).
If the initial assessment of flood susceptibility of the overlay analysis of geologicalgeomorphological indicators is moderate to High (Input 3 = 2.8) and the height of the river bank above the water level is 0.5 -1m (Moderate to High flood susceptibility, Input 1 = 3.3), and the flood susceptibility according to the river bed type is High to Moderate (Input 2 = 3.7), then the Complex geological-geomorphological flood susceptibility is High (Output = 70).
On the other side if the river bank is higher than 3.5 m (Input 1 = 0.6, Very low flood susceptibility) and the flood susceptibility resulted by the overlay analysis of geological-geomorphological indicators is Moderate to High (Input 3 = 3), although the high value of Input 2 = 4.3, then the Complex geological-geomorphological flood susceptibility is Low (Output = 39).
The observed changes in the flood susceptibility rates after applying suggested logic frame show increasing of the of the susceptibility rates and in some cases could lead to overestimating the susceptibility but following the precautionary principle in disaster risk reduction we consider that the simulation results can be used at the first stage of flood management process and are helpful for informative decision making.

Conclusions
The applied weighted sum overlay analysis of flood susceptibility rates according to the lithology, slope, fluvial relief and distance from streams for the investigated areas of the river Struma valley (municipalities Kresna and Strumyani) shows high and very high flood susceptibility in the lowest parts of the river Struma valley, south of the town Kresna to the Bulgaria -Greece border.The geological-geomorphological features of the area influence the redistribution of the runoff of the surface and the underground component as well as the time for which the slope flow reaches the riverbed.Regarding the assessment of geological prerequisites for floods special attention have to be given on the indicators related to rock resistance of weathering and erosion, porosity and cracks that influence the water permeability and flood retention capacity.About the geomorphological prerequisites the most sensitive are floodplains, alluvial and proluvial fans.The result of the analysis provides information about flood sensitive areas and can be used as the basis for organizing monitoring to minimize the risk of flooding.The data obtained should be interpreted in relation to the hydro-climatic conditions and, particularly, taking into account the quantity and duration of precipitation.Having regard the limit of the available data application of fuzzy logic gives good result.The suggested flood susceptibility model can be used as a basis for further development of a complex flood susceptibility assessment by adding other indicators as precipitation, stream flow, soils, vegetation etc. and can be considered as a first step of flood hazard assessment.
The fuzzy logic model is designed as a two-level hierarchical system with three inputs and one output.Each level of the system is consisted of one fuzzy logical subsystem with two inputs.The fuzzy logic system output gives the complex assessment for geological-geomorphological flood susceptibility of study area regarding the interaction of the three inputs.
The fuzzy logic model for flood susceptibility rate can be used as a basis for further development of a complex flood susceptibility assessment by adding other indicators as precipitation, stream flow, soils, vegetation etc. and can be considered as a first step of flood hazard assessment.

Table 1
2, Input 3, Intermediate variable and Output are presented in

Table 1 .
The numeral values of trapezoidal membership functions and the corresponding linguistic variables of the Input 1, Input 2, Input 3, Intermediate variable and Output

Table 2 .
The logic matrix of the first fuzzy logic subsystem

Table 3 .
The logic matrix of the second fuzzy logic subsystem

Table 4 .
Flood susceptibility according to the lithology

Table 5 .
Flood susceptibility according the slope of the topographic surface

Table 7 .
The sub-intervals are defined as follows: IF the value of "Complex geological-geomorphological flood susceptibility" ϵ [0, 20) THEN the "Flood susceptibility rate" is Very low; IF the value of "Complex geological-geomorphological flood susceptibility" ϵ [20, 40) THEN the "Flood susceptibility rate" is Low; IF the value of "Complex geological-geomorphological flood susceptibility" ϵ [40, 60) THEN the "Flood susceptibility rate" is Moderate; IF the value of "Complex geological-geomorphological flood susceptibility" ϵ [60, 80) THEN the "Flood susceptibility rate" is High; IF the value of "Complex geological-geomorphological flood susceptibility" ϵ [80, 100] THEN the "Flood susceptibility rate" is Very high.

Table 7 .
Simulation results for the rate of the Complex geological-geomorphological flood susceptibility from the fuzzy logic model