simulation of flood hazard using gis-based cellular automata (case study: chirchir catchment)

introduction flood is an inevitable natural phenomenon occurring from time to time in all rivers and natural drainage systems, which not only damages the lives, natural resources and environment, but also causes the loss of economy and health. so estimation and prediction flood hazard is very important spatially in the watersheds without measurement station. there are many models in the water and environmental studies which have been used for runoff and flood in the watersheds without measurement station. one of the newest is cellular automata that has been combined well with the gis for simulation runoff and flood hazard. cellular automata as a tool for modeling and simulation of processes taking place in the real world are now increasingly used, as evidenced by their use not only as a tool for creating simulations, but also by their use in areas of crisis management. using gis knowledge, it is possible to create cellular automata, which can appropriately and authentically reflect the water flow on the earth’s surface. cellular automata (ca) is a mathematical model that can be used for computation and simulation of the systems. in this method, the basin is defined with a network of the rectangular cells, and the interactions between the cells together with the geographic rules that govern the area result in the runoff modeling. this model relies on the gis and satellite images. cellular automata model uses various information such as digital elevation model (dem), land use, hydrologic soil groups, rainfall, slope and etc for runoff estimation. in the present study, the runoff of the chirchir catchment in east azarbaijan province, iran, has been modeled by means of the gis-based cellular automata. materials and methods in this study, gis-based cellular automata were used to simulate flood in the chirchir catchment in the northwest of iran. ca models use several primary components including the cells arranged in a regular mosaic pattern (raster, grid), transition rules determining the changes in cell properties, neighborhood of the cell, and boundary conditions. these components affect the status of each individual cell in a network in a given time step. in this research, chirchir catchment in east azerbaijan province is modeled using cellular automata. first, scs formula is used to predict the runoff in each cell. map of hydrological soil groups of chirchir catchment is determined by means of soil texture map, then land use and shg maps are prepared for calculating the runoff curve number (cn) map for the normal conditions. since the soil has dry moisture condition and the slope is greater than 5%, adjusted cn is calculated for dry antecedent moisture condition and catchment slope using the common relationship (ponce and hawkins, 1996; hawkins et al., 2009). after reading the rainfall and the cn map for dry antecedent moisture condition, the runoff was calculated using the scs equation. then kinematic wave model is used for flow depth in the cells and runoff production within each cell is simulated by determining the cell state (water surface elevation) that included both the cell altitude and the water depth. the distribution of water flow among cells was determined by applying ca transition rules based on conservation of energy and continuity equations. d8 algorithm is used to simulate flow direction during the calculation of the surface convergence. the procedures for channel network delineation are based on the d8 model for flow over a terrain surface represented by a grid dem. in this model, a single flow direction in the direction of steepest slope towards one of the eight (cardinal and diagonal) grid cells neighboring is used to represent the flow field. also for calculating flood hydrograph, travel time is calculated using flow length and flow velocity. so roughness coefficient and flow depth is used for flow velocity and then travel time map is obtained. finally, python programming language is used to estimating flow hydrograph due to simplicity, powerful and object-oriented programming language and supported with gis. results and discussion results show that the most area of the chirchir catchment has pasture and type d soil hydrological group, so that it has very low permeability which means that a large amount of rainfall convert to runoff. runoff depth is high in east and southeast of the chirchir catchment due to physical characteristics and rainfall of the catchment but among these parameters, the slope was the most important parameter in the runoff generation. also, map of the flood hazard shows that downstream river has high potential in flood hazard due to receiving water from upstream. then, for simulation flood hydrograph, travel time was calculated using ratio between flow lengths to flow velocity. flood hydrograph estimated for two events, 27/3/1388 and 12/3/1386 so that the computational runoff is very adapt to observational runoff. the correlation coefficient for two events (0.82 and 0.70) indicate the good accuracy of the model. low error rates also indicated that the cellular automata model has the high efficiency to predict the flood peak and the time of its occurrence in the chirchir catchment. the results of this study matches with the results of researchers such as aboudagga (2005), rinaldi et al. (2012) and cirbus and podhoranyi (2013). they stated that the use of cellular automata model with compared to the conventional methods that use gis environment, has higher accuracy and capable to estimate flood hydrograph. therefore, the use of cellular automata with gis, not only speed up the calculation of runoff, but also increases the accuracy of the results. conclusion comparison of the results with the observation proves that the results are well accurate. besides to the advantages of this method in simplicity and implementation of the realistic rules, this method is good at gaining the runoff data at any point of the basin besides to the exit point. good agreement between the model output and the empirical measurements revealed that a ca approach can provide realistic results for a complex natural process like flood. keywords: cellular automata, chirchir catchment, flood hazard, gis