Flash flooding risk impacts can be reduced through the implementation of mitigation strategies plan (MSP) for flood management. This study aims to develop a flash-floods risk mitigation plan, which appears to be beneficial for municipalities, provincial administrators, and authorities to reduce the impact of the flash flood in the Kyrenia region, Northern Cyprus. In this work, rainfall data were collected from the nearest stations for 22 years. The return periods of maximum daily rainfall are also determined by using six formulas. Furthermore, flood inundation and hazard maps were defined by utilizing SAGA, QGIS, ArcGIS, 2D HEC RAS, and HEC -HMS software then determining the degree of risk and identifying strategies based on quantitative risk analysis by developing a risk matrix. As a final result, catastrophic risk areas are distributed significantly downstream. In conclusion, the proposed flash flood mitigation plan includes strategies to reduce flood losses of human life and constructed structures across Kyrenia and proposed hazard and inundation risk maps to assess planners and decision-makers for the potential impact of floods to avoid.   Keywords: DEM; Flash Flood; Hazard map; HEC-RA Mitigation plan; Risk matrix

Kyrenia region is in the northern part of Northern Cyprus that is environmentally fragile and susceptible to natural disasters. Thus, the study of frequency analysis is essential to find the most suitable model that could detect the region’s risk in certain natural phenomena such as rainfall, flood, and so on. The objective of this research is to determine the best fit probability distribution in the case of average daily rainfall and total rainfall using 22 years of data (1995-2016) from the Kyrenia region in Northern Cyprus by using 37 probability distribution models. The best-fit probability distribution in the case of maximum annual daily rainfall is determined using various distribution types. Three goodness-of-fit test statistics were applied. Beta, Dagum, Wakeby, Paretoa, Log-Pearson 3, Gen. Extreme Value, and Gen. Gamma (4P) showed the largest number of best-fit results. The results of this study can be used to develop more accurate models of flooding risk and damage.   Keywords: Distribution function; goodness-of-fit tests;  Northern Cyprus; rainfall.