Nowadays, with significant climate change, the trend of environmental hazards is increasing. In the meantime, floods have shown a growing trend than other hazards. Haraz watershed in northern Iran is prone to floods due to the heavy rainfall with irregular pattern. Therefore, combining different methods and examining new approaches is an essential step in the development of methods in this field. In the present study, Analytical Network Process, Analytic Hierarchy Process and Fuzzy Analytic Hierarchy Process models were combined with Ordered Weighted Average, Weighted Linear Combination, Local Weighted Linear Combination models to prepare a flood risk map. The performance of two new models, Weighted Multi-Criteria Analysis and Geo-Technique for Order of Preference by Similarity to Ideal Solution, was also evaluated in this field. The results of the models showed that in general the basin is in a moderate risk situation. Meanwhile, the south-eastern parts of the basin show a high flood risk situation. Also, by comparing the models, it was found that the combination of multi-criteria models and the use of Weighted Multi-Criteria Analysis and Geo-Technique for Order of Preference by Similarity to Ideal Solution models are very effective and efficient for preparing flood risk maps.

Iran is one of the most flood prone areas in the world. The spring flood of 2019 was recorded one of the most devastating flood events in northern region of Iran. In this study, Sentinel-1, Sentinel-2, Sentinel-3 and Landsat-8 images were used to extract the flood map. Then, flood maps of these areas were prepared using Random Forest (RF) algorithm for Sentinel images and Support Vector Machine (SVM) algorithm for Landsat-8 images. In addition, flooding in these areas was assessed using the Fuzzy Best Worse Model - Weighted Multi-Criteria Analysis (FBWM-WMCA). The results of FBWM model showed that the criteria of precipitation, slope, height, land use, drainage density and distance from channel were the highest and the criteria of Curvature, Geology, Topographic Wetness Index (TWI), Stream Transport Index (STI), Stream Power Index (SPI) and The Topographic Ruggedness Index (TRI) played the lowest role in flooding in these areas. According to the FBWM-WMCA model, 38% of the Gorgan watershed in the northern, northwestern, western and southwestern parts and 45% of the Atrak watershed in the eastern, northeastern, northern and western parts are in high flood risk. The overall accuracy of the 2019 flood maps in Gorgan watershed for Sentinel-1, Sentinel-2, Sentinel-3 and Landsat-8 images is 89, 87, 80 and 85% and for Atrak is 91, 88, 82 and 86 percentages respectively. In general, based on the results of this study, FBWM and FBWM-WMCA models are effective and efficient for determining the weight of criteria and preparing flood risk maps, respectively.