Proc.of the gthWSEAS Int. Conf. on Mathematicaland ComputationalMethods in Scienceand Engineering,Trinidad and Tobago, November 5-7, 2007

The influence of hydrological parameters on the stream floods of a mountainous area in ,

P. LOKKASI, S. KOTSOPOULOS2, J. ALEXIOU2, G. GRAVANISI, V. VASSILOGLOUI. S. MAGALIOS3. V. KASSOS4 1: Department of Civil Engineering Infrastructure Works, Technological Educational Institute (TEI) of , GREECE 2: National Agricultural ResearchFoundation, Institute of Soil Classification and Mapping of Larissa, GREECE 3: MScE Surveyor Eng, Larissa, GREECE 4: Electro-MechanicalEng, Larissa, GREECE p. lokkas @tei lar.gr http ://users.teilar. grI -p. lokkas

Abstract: - The sloping area on the North and North-Eastof ,opposite , in the region of Thessaly,Greece, receives considerable rainfalls. However, despiteits extendednatural vegetation, severeproblems have arisendue to floods. The existenceof theseproblems is mainly basedon two reasons:a) on the peculiarity of rainfall eventscombined with the local conditions(hydrologic, geologic, land use, etc.) and b) on the existing infrastructureworks. In order to analyzeand interpretthe phenomenaof local floods, in 'oArchimedes" the frame of the researchprogram which is materializedby the TEI of Larissaunder the title "Spatial mapping and estimationof hydrological risk, emphasizingthe floods and draughtsin urban and non urban areas of Thessaly and their environmentalimpacts" it is essentialto record, evaluateand validate various hydrologic parametersof the correspondingstudy area.The GeographicInformation Systems(GIS) is utilized to analyzeand validate the information abouttopographic, hydrologic, geologic and soil data, along with land cover and infrastructure.Finally, the GIS are also used to produce and map complex hydrologic parameters,which, in tum consistinput datafor hydrologic modelsas well as presentationof results.

Key-Words:- Hydrology, geology,vegetation, stream, flood, infrastructure.

1 Introduction . the particularity of the rainfall events in this The region of the North andNorth East slopesof the region, mount Ossa,in Thessaly,Greece, is coveredby rich r the local conditions, i.e. hydrologic, geologic, and densenatural vegetation, forest trees, shrubs etc, land use.etc and especially chest nut trees, which, locally, is o the lack of the necessaryinfrastructure works. considered a valuable means for agricultural The environmental consequencesof the above development.The level of rainfalls in the area is floods are highly crucial for both the cultivatedareas considerably high; almost three times higher and the local delta of Pinios river. comparedto those taking place on the main part of For recordingand managementof the necessaryfor EasternThessaly. This situation,combined with the this study local data, the Geographic Information local hydrologic conditions and the existing SystemGIS-ArcView has been utilized. The use of infrastructure, has resulted in serious problems GIS [1] is advisable for any presentationof the concemingdamages due to excessiveflood events. spatial variables and for this reason it has been In the frame of the research program treated as the main tool for a large number of "Archimedes - EPEAEK II", co-funded by the applications,like soil-water resourcesmanagement European Social Fund & National Resourcesand L2l,l3l, [4], [5] etc. realizedby the Dept. of Civil Engineeringof TEI of As far as estimation of flood discharge is Larissa, Greece,the above region has been selected concemed,it has to be noted that there are spatial as a study area. The work is aiming at a fluent parameters,that after processing,create input data evaluation of hydrological parameters,which are for flood simulation [6] throughhydrological models relatedto the estimationof flood discharges. [7], [8], like soil data,data ofvegetation cover, land Theseflood discharges,which actually constitute use,relief, streamnetwork characteristicsetc. the main problem,take placeoccasionally during the The processingof the above data,may directly lead winter and are due to the large amount of water in to hydrological parameters, like the area of the streams.The problemtherefore, is originatingto: Proc.of the gthWSEAS Int. Conf. on Mathematicaland ComputationalMethods in Scienceand Engineering,Trinidad and Tobago, November 5-7 ,2007 168

hydrologic basins,the slopesofthe basinsand their position of each work, through GPS, were water streams,runoff coefficientsetc. saved in corresponding files. There were totally kept 46 positions of stream bridges, ranked from very small to very large. Each 2 Procedure one of the correspondingstreams was mainly For recording and plotting both the existing connectedto one or more watersheds. information and the resulted conclusions, the o Topo€traphicdata for plotting the slopes, Geographic Information System GIS-ArcView has along with the correspondingcross - sectional been utilized, so that the spatial variability of areasof the main streamsthat belong to the different parametersmay easily be presented.The study region. GIS is a powerful technology approachableand All the abovecartographic information, related to accessibleto many users.This technologyallows for the study area of mount Ossa is shown in Fig. l, automationin many procedures,therefore working where, the watershedswere drawn with a bold in a complementary way for many analytical contour. The information was then digitized and methods, which are utilized for aquatic resources convertedto GIS-ArcView format file, where, the management[]. inter-correlationsamong the geographic elements, For analysis,estimation and evaluation of flood i.e. their topologicalrelations [1] along with their eventson the study area,the collection and use ofa databaseconnected to them were developed. variety of data is necessary.Some of this data may On the top of them, in GIS files there were concern hydrological, soil and geological createdsecondary information, like: characteristics,use of land, vegetation cover, and o the borders of watersheds of the most terrain, along with the existing infrastructureworks. important water streams, along with their area, The last one may comprisethe road network, the averageslope and database, cross section of natural water streams. technical o the databaseof water streamsaccompanied by works and crucial cross-sections that affect the flow their length, average slope and the corresponding of water. locations of the selectedcross sections,along with For this reason the following tools have been their geometricalcharacteristics, utilized: r the databaseof soil samples,taken out of . Maps of Hellenic Army GeographicalService selectedlocations ofthe study area. (HAGS),scaled at l:5000and l:50,000, All these data were then combined and worked . Geological maps of the Institute of Geology out to give results,which in tum becamenew data and Mineral Exploration (IGME), scaled at necessaryto simulatethe naturalsystem. 1:50,000, The elaborationof elementscomprises: . Soil maps of the Ministry of Agricultural . Topological processing of hydrological and Development and Foods (MADF), scaled at topographicdata, i.e. hydrographicalnetwork, l:50,000, watershedareas, slopes of earth and water . A silviculturalmap scaled at 1:20,000and streams,as well as length and cross sections o A mapof mountOssa's soil cover,i.e. borders of water streams. of forest vegetation species, scaled at . Topological processing of pedologic maps 1:350,000. along with their correspondingdata. Apart from the above official data it was . Topological processingof silvicultural maps necessary to collect data from local held and mapsof vegetationcover. measurements; for this reason the following The averageslope of the watershedshas been activitieswere realized: calculated through the equation of Gregory & r Soil samplestaken from different locationsall Walling,as described by Tsakiris[10]: over the study area.The depthof samplesmay \,L.D ^S^-- . (r, vary in the rangesof 0-30, 30-60 and 60-90 .E cm. Every location of samplingwas recorded whereS" is the averageslope of the watershed,)Z is through a Geographic Positioning System the total length ofthe contour lines, D is the contour (GPS). interval of lines andE is the areaof the watershed. o Elementsof technicalworks, having the form The outflow of the watershedcould be estimated of stream bridges, along the axis-road through the Soil Conservation Service (SCS) Omolion - Stomion- Karitsa, were recorded, method, where, an estimation of the runoff curve i.e. their crosssectional details alons with the number,Cl{ is required[9], [10],[1 1], [12], [13]. proc.of the gthWSEAS Int. Conf. on Mathematicaland ComputationalMethods in Scienceand Engineering,Trinidad and Tobago, November 5-7,2007 169

AN

Fig. 1: Topographicmap of the local studyarea with the contoursof watersheds.Scale l:165,000

This coefficient is given in tables and varies The elaborationof cartographicdata, along with according to the land use, to the managementor the data and types of earthy plotting, result in the practise,to the hydrological conditions and to the estimationof: typeofearth [9], [10], [1 1], [12], [13]. o the area along with the averageslope of the In the casethat there are more than one land uses hydrologicwatersheds, within the watershed,then, the so-calledcomposite o the slopes,the lengths and the cross sections curve numberoCN", must be taken into account.This of water streams. number is an areal weighted average of the Figures 2 and 3 representgraphically the above watershedsub-regions and is calculatedthrough the results. equation An analysisof the soil samplestaken from depths I cl'r,'r n up to 90 cm, wherever feasible, leads to the CN,=#=lCN,.e,, (2) conclusion that, according to the final infiltration , L: H i=l rate estimates,the majority of the hydrologic soil where Ci[ representsthe runoff curve number for groupsbelong to typeA [9], [0], [11]. the sub-regioni, E; is the watershedarea for sub- In relation to the hydrological conditions, their region i, and e; is the relative area of sub-regionl, characterizationas normal, moderate or poor, it over the total watershedarea (Iar: 1). resultsfrom the combinationof the earth slopes,the The composite curve number, CN", for every depths, the earth profile and the density of watershed has been calculated through the GIS- vegetationcover. ArcView program, provided, all the partial runoff Therefore,for each sub-region,i.e. region which curve numbersfor every soil cover categorywithin is coveredby different kind ofvegetation, the value into account. eachwatershed, have been taken of runoff curve number, CN of SCS is defined [9], [10], [117,ll2], [3], and then,through a topologic processing,the composite value of CN, CN", of 3 Results every watershedis obtained(Fig. 4). In the abovestudy areathere are various species of forestvegetation and pasturesas shownin Fig. 2. Proc.of the gthWSEAS Int. Conf. on Mathematicaland ComputationalMethods in Scienceand Engineering,Trinidad and Tobago, November 5-7 , 2007 170

- BOUNDARYOF BASIN ,1 SAMPLINGPOINT -BoUNDARYoFBASIN - , 0.22 AVERAGESLOPEOFWATERSHED - WATERCOURSE WATERCOURSE N/ CONTOURLINES 1OO /\/ CONTOURLINES 1OO F'N FIR FORESTDENSE N:] OAK FORESTDENSE 25,0 COMPOSITERUNOFF rii:!r?FIR FORESTSPARSE 'm CHESTNUTFQREST DENSE CURVENUMBER rT- BEECHFOREST DENSE [TIq-,:6HESTNUT F6REST SPARSE r----r BEECHFOREST SPARSE F.,j.:rPASTURE ,@ EVERGREENFOREST ffi CULTIVATEDOR BARE LAND (SCRUBBY) Fig. 4: Map of compositerunoff curve numbercCNc per watershed.Scale l:80,000 Fig. 2: Part of the study areawith categoriesof soil cover and soil samplinglocations. Scalel: 80.000 4 Conclusion The GeographicInformation Systemshave been proven a very useful tool for the abovehydrological application in the region of the study area, located on the North and North Eastslopes of mount Ossain Thessaly,Greece. It has been especiallyuseful for tracing a mass of cartographic data, i.e. terrain, hydrographicalnetwork, vegetationcover - or land usefrom the geologicalor pedologicpoint of view. It is also useful for process,management and valuation of hydrologicalparameters, along with the area,the slopesofhydrologic watersheds,the slopes, lengths and cross sectionsof water streamsas well asthe runoff curve numberCNof SCS. More and probably more useful conclusionson the aboveresearch are expectedto be obtainedafter a full completion of the program,when, the spatial variation of flood risk in the study area,is going to - BOUNDARYOFBASIN be refereed. - WATERCOURSE /\./ CONTOUR LINES 100 0-0.25 VALUES OF AVEMGE SURFACE SLOPES r.-r s0.25 tr',':r0.41-0.50 References: r- 0.2$0.40 w >o50 [1] Bunough P. A., Principles of geographical information systems for land resources assessment,Clarendon Press, 1986 Fig. 3: Map of averagesurface slopes per category KalfountzosD., I. Alexiou, S.S.Magalios and P. ofsoil cover.Scale l:80,000 [2] Vyrlas, Integratedwater resourcesmanagement Proc.of the gthWSEAS Int. Conf. on Mathematicaland ComputationalMethods in Scienceand Engineering,Trinidad and Tobago, November 5-7,2007 171

using GIS: A case study in TOEB Pinios, Proc. Kotsopoulos, P. Georgiou and G. Terzidis Int. Conf. Protection and Restoration of the "Rainfall - runoff relationshipsof the Almopia Environment IV, Halkidiki, Macedonia,Greece, watershedusing HEC-I ,2'o National Conference 1998,pp. 83-90. of the Hellenic Committeefor Water Resources on Integrated Interstentions [3] KalfountzosD., I. Alexiou, S.S. Magalios,P. Management for Vyrlas and G. Tsitsipa, An expert system of Flood Risk Decrease, National Technical integrated management on inigation water: University of Athens,Greece, 1995 (in Greek). Specialapplication of the Local Organisationfor [7] Hydrologic EngineeringCenter, HEC-I, Flood Land Reclamation of Pinios, Greece, Hydrograph Package, U,S Ar*y Corps of Proceedingsof the 4'nNational Conferenceof the Engineers,Davis, CA, USA, 1980. Engineering Center, Hydrologic Hellenic Committee for Water Resources t8l Hydrologic Managementon ll'ater ResourcesManagement Modeling SystemHEC-HMS, US Army Corps of Davis,CA, USA,2002. for Sensitive Regions of Hellenic Tenitory, Engineers, Volos, 1999,pp. 20-27(in Greek). [9] Soil ConservationService, National Engineering Washington DC, Section 4, [4] KalfountzosD., D. Pateras,L Alexiou and S.S. Handbook, Magalios,Soil Surveyof XeromeroDistrict using Hydrology,1972. GIS for the Supportof IntegratedSoil and Water [10] Tsakiris G. Water Resources,I. Technical Resources Management, International Hydrologt, Symmetria Publications, Athens, Conference on the Agriculture, Fisheries and 1995(in Greek). Agro-Industry in the Mediterraneanwith Special [11] KotsopoulosS. Hydrologt, Ion Publications, Referenceto Islands,Valetta, Malta,1997. Athens,2006 (in Greek). V., D. Maidmentand L. Mays, Applied [5] KalfountzosD, P. Vyrlas, S. Kotsopoulos,I. [2] Chow Alexiou and A. Panagopoulos, Inigation Hydrologt, McGraw-Hill International Book managementusing GIS in Aisonia District, Proc. Company,NY, 1988. Intl. Conf Protection and Restoration of the tl3l Wanielista M., R. Kerten and R. Eaglin EnvironmentVI, Skiathos,Greeceo Vol. III, 2002, Hydrologt: Water Quantity and Quality Control, pp. 1689-1696. J. Wiley & Sons,NY, 1997. t6l Papamichail D., C. Babajimopoulos, S.