Comparative Analysis of Developed Rainfall Intensity–Duration–Frequency Curves for Erbil with Other Iraqi Urban Areas

: Rainfall Intensity–Duration–Frequency (IDF) relationships are widely used in water infrastructure design and construction. IDF curves represent the relationship between rainfall intensity, duration, and frequency, and are obtained by analyzing observed data. These relationships are critical for the safe design of ﬂood protection structures, storm sewers, culverts, bridges, etc. In this study, the IDF curves and empirical IDF formulas for the city of Erbil were developed for the ﬁrst time by employing the annual maximum rainfall data for a period of 39 years (1980–2018), which is the only available recorded data. Statistical techniques such as Gumbel and Log-Pearson Type III (LPT III) were utilized to determine the IDF curves and empirical equations from daily rainfall data for several standard durations and return periods. The correlation between the rainfall intensities obtained from IDF curves and the empirical formula presented a reliable match, with a coefﬁcient of determination of (R 2 = 1). The results were compared to previously developed IDF curves and empirical formulas in Iraqi cities to show their reliability. Moreover, the results can be an initial step for authorities to establish required guidelines in the studied area, and in the design process of the storm water infrastructure of urban basins in the future.


Introduction
Rain is an important source of fresh water on Earth, provided by the long hydrological cycle; however, it can also be a source of natural disasters, such as urban flash floods. Al-Amri and Subyani [1] claimed that water project design is difficult in arid regions because rainfall is unpredictable and erratic in both time and space. This is also true for semi-arid regions. The precipitation varies significantly between the years, which makes it a challenge for water resources management [2]. Rainfall intensity-duration-frequency (IDF) relationships are graphical representations of the amount of water that falls within a given period of time in catchment areas [3]. Thus, they are critical for preventing flooding and, as a result, minimizing the loss of life, property, and the cost of water damage insurance and weather-related risk assessments. IDF curves have been estimated in developed countries since 1932, but they are still not available in many developing countries due to a lack of data. As Koutsoyiannis et al. [4] pointed out, the availability of accurate and long-term rainfall data in developed countries and the lack of such information in developing countries is the major factor in the construction of IDF curves. Many researchers in different parts of the world have developed IDF curves and empirical formulas [5][6][7][8][9][10]. Lima et al. [11] came up with a new way to figure out the parameters of a Generalized Extreme Value (GEV) distribution at the sub-daily level. They used the scaling-invariant property of rainfall duration versus intensity. When they looked at historical data, they found that the model was good at predicting IDF curves for both gauged and partially gauged sites. Global Satellite Mapping of Precipitation (GSMaP) and the Bartlett-Lewis rectangular pulses (BLRP) model, which were used by Sun et al. [12], disaggregate the daily in situ rainfall, which was then used to derive more accurate IDF curves. Kristvik et al. [13] used the concept of scale invariance to downscale the distribution of extreme precipitation from daily to sub-daily timescales in three Norwegian cities. From that, they derived the local-scale intensity-duration-frequency (IDF) curves for future precipitation. Some studies looked into how satellite-based precipitation could be used to make IDF curves [14,15]. Schardong et al. [16] conducted their research with the help of a web-based tool for the development of IDF curves under changing climate conditions at gauged and ungauged locations. Using the controlled random search method, Weinerowska-Bords [17] analyzed the practical issues associated with developing IDF curves and formulas (such as sample homogeneity and independence, and storm episode selection) for the development of the Gdańsk IDF-formula. Krvavica and Rubinić [18] investigated the effects of design storms and critical rainfall durations on flood prediction in partially urbanized catchments. In Saudi Arabia, which is considered an arid region, and the neighboring country of Iraq, in the southwest, researchers developed IDF relationships with different return periods in different parts of the country [1,[19][20][21][22]. In Turkey, Deger and Yuce [23] analyzed the rainfall intensity-frequency of Gaziantep city, which is located in the northwest of the studied area. They expected their results to be useful for the local authorities.Şen [24] proposed a new method for constructing IDF curves based on the Annual Daily Maximum Rainfall (ADMR) data. The suggested methodology was applied to four ADMR records from meteorology stations in Turkey, in the Ceylanpinar region next to the Syrian border. In the west of the study area, using the Gumbel method, Hadadin [25] found a correlation between rainfall data and the intensity and duration of rain in the Mujib Basin in Jordan. In Iraq, starting from the north, Hamaamin [26] developed rainfall intensity-duration-frequency curves for Sulaymaniyah city, which is situated to the east of Erbil city. The study's findings revealed that the IDF curves computed with the coefficient of determination (R 2 = 1) fit perfectly with the predicted values of rainfall intensities. Rainfall data for a number of stations in Northern Iraq was analyzed by Hasan and Saeed [27]. Mahdi and Mohamedmeki [28] found that the maximum rainfall intensity of 118.052 mm/hr occurred at the duration of 0.25 h of the return period of 100 years when they analyzed IDF curves of the capital city Baghdad in the middle of Iraq. In the south of Iraq, IDF curves were estimated for the cities of Basrah, Najaf, and Nasiriyah by [29][30][31].
There have been more urban flash floods since 2010 in Erbil, and this has caused many sewer pipes to overflow, allowing sewage to escape and mix with storm water. As a result of the flooding, many small and local power generators lost their ability to generate electricity, resulting in more outages [32]. Mustafa et al. [33] presented extreme rainfall events as a cause of urban flash floods in Erbil. Moreover, according to Mustafa and Szydłowski [34,35], the city's urbanization is causing an increase in runoff peak and volume. As a result, the possibility of flash floods increases, posing a risk to life quality in the studied area, particularly in the city center. In fact, the hydrological monitoring system is important for the modeling of the flood risk in urban areas [36]. However, researchers have studied the effect of urbanization on floods not only in Erbil, but also in other cities of Iraq. Recently, Shareef and Abdulrazzaq [37] investigated a river flood to mitigate flood risk in the Fallujah district in Anbar Province, which is 69 km to the west of Baghdad, the capital city of Iraq. Some Iraqi cities have already developed their IDF curves and empirical formulas, but Erbil is the only one that has not yet done so.
The main objective of this study was to develop IDF curves for Erbil city center using daily rainfall depth measurement data from Erbil meteorological station. Additionally, empirical formulas were generated to estimate rainfall intensity for various return periods and rainfall durations. The obtained results were compared to previously developed IDF curves and empirical formulas in urban areas of Iraq, in order to show the reliance of the findings. While we recognize that the methods and equations used in this study are standard and have been used in other studies, the novelty of our work is that we developed IDF curves and formulas for a location where such work is not available. In particular, the findings of this study could help to improve water resource management in local urban areas, as well as the design and operation of storm water management systems. Moreover, the results of the study aim to assist authorities in Erbil municipality and the General Directorate of Water and Sewerage in implementing strategies and making good decisions.

Study Area and Data Collection
Erbil is the capital city of the Kurdistan Region of Iraq (KRI) and the second biggest city in the north of Iraq after Mosul. It lies approximately 350 km (220 miles) from Baghdad, the capital of Iraq. The total area of Erbil Province is 14,873.68 km 2 ; however, the area considered as the studied area in this work is the center of the city ( Figure 1). It lays on the coordinates: 36.11 28 latitude and 44.00 33 longitude, with an average elevation of 406 m above sea level [38]. There are long, hot summers and mild winters in Erbil, which has a Mediterranean climate. Aridity prevails during the summer months, with little to no precipitation between June and September. Wet and humid weather dominates the winter, with January being the wettest month of the year. The average annual temperature is around 21 • C, with temperatures ranging from 7 to 47 • C. Approximately 400 mm of rain falls on average each year, ranging from 200 mm in dry years to 650 mm in wet years, with the majority of the rainfall starting in October, increasing through December to January, and continuing until the end of May. The meteorological dataset for this study has been obtained from the Climate Forecast System Reanalysis (CFSR) of the National Center for Environmental Prediction and General Directorate of Meteorology and Seismology in Kurdistan region/Iraq (GDMS). From 1980 to 1991, the dataset was obtained from the Climate Forecast System Reanalysis The meteorological dataset for this study has been obtained from the Climate Forecast System Reanalysis (CFSR) of the National Center for Environmental Prediction and General Directorate of Meteorology and Seismology in Kurdistan region/Iraq (GDMS). From 1980 to 1991, the dataset was obtained from the Climate Forecast System Reanalysis (CFSR) of the National Centre for Environmental Prediction. The General Directorate of Meteorology and Seismology in Kurdistan region/Iraq (GDMS) provided the dataset from 1992 to 2018. A classic rain gauge inside Erbil was used to record the rainfall data consistently and daily. The dataset consists of annual maximum daily (24 h) rainfall (mm/h); over 39 years from 1980 to 2018 were taken from [33] (see Table 1). The study area is big in size, a series of meteorological stations is required, but there is only one meteorological station available.

Precipitation Duration Reduction Formula
The only available rainfall data for Erbil city is the annual maximum 24 h measurements for a period of 39 years. Indeed, there are two ways to develop IDF curves for standard periods: either by obtaining detailed rainfall data from meteorological stations on short timescales (such kinds of data can be obtained only from automatic rain gauges) or by downscaling annual maximum rainfall from classic rain gauges to shorter durations. In order to obtain IDF curves for the study area, annual maximum rainfall data for standard periods of 10, 20, 30, 60, 120, 180, 360, and 720 min is required. Due to the lack of these standard times, the Indian Meteorological Department (IMD) utilized Equation (1) to determine the rainfall data for periods shorter than 24 h [39]. The IMD method is simply the application of an empirical formula for downscaling the amount of precipitation falling on a daily or 24 h basis into shorter-duration precipitations. Rainfall analysis using the IMD method can be performed in stations where there are no automated rain gauges that can record both rainfall amount and duration data at the same time, but where there are rain gauges that can record data of cumulated 24 h rainfall amounts: where P i is the required precipitation depth for a duration less than 24 h in mm, P 24 is the daily precipitation depth in mm, and t is the required duration time in hours.

Frequency Distribution and the Development of IDF Curves
An analysis of the dataset for the past 39-year period has been conducted utilizing existing rainfall data and statistical tools to derive the intensity-duration-frequency relationships. Estimates of rainfall intensities over various timescales and recurrence intervals have been used to construct the IDF relationships. For example, the Generalized Extreme Value (GEV), normal, lognormal, Gumbel, and Log-Pearson Type III distributions are widely used in many parts of the world [3,9,21,22,25,[40][41][42]. The relationship between rainfall intensity, storm duration, and return periods from rainfall data for the study area was developed using the following two common frequency analysis techniques.

Gumbel Distribution
The Gumbel method was established by the German mathematician Emil Gumbel, and has been widely used for modeling extreme events in hydrology and other fields [43,44]. Kotz and A Nadarajah [45] published a book that details this distribution and lists over 50 applications, ranging from accelerated life testing to earthquakes, floods, horse racing, rainfall, supermarket queues, sea currents, wind speeds, and track race records. This method was chosen to perform the flood probability analysis because it has become the most widely used distribution for IDF [46][47][48][49]. As Elsebaie [22] pointed out, the Gumbel method is a simple and straightforward approach that can be applied in extreme events (maximum value or peak rainfalls). The frequency of precipitation P T (in mm) for each duration with a specified return period T (in years) is given by the following equation: K is the Gumbel frequency factor, which can be found by: where P ave is the average of the maximum precipitation corresponding to a specific duration. The arithmetic average in Equation (1) is used when using the Gumbel distribution.
The individual extreme value of rainfall is P i , and the number of events or years of record is n. The following equation is used to calculate the standard deviation: where S is the standard deviation of P T data. The frequency factor (K), which is a function of the return period and sample size, when multiplied by the standard deviation gives the departure from the average rainfall of a desired return period. The rainfall intensity I T (mm/h) for the return period T d is then calculated as follows: where T d is the duration in hours.

Log-Pearson Type III
As previously stated in this section, there are several methods for determining the intensity-duration relationship, with Log-Pearson Type III (LPT III) being one of them. LPT III contains a logarithm of the measured values. The mean and the standard deviation are determined using the logarithmically transformed data. The precipitation frequency is calculated using the LPT III method in the same way as the Gumbel method. The following is a simplified expression for the latter distribution: P * T , P * ave , and S * are the same as before, but based on logarithmically transformed P i values; i.e., P * of Equation (7). The Pearson frequency factor K T is determined by the return period (T) and the skewness coefficient (C s ).
The frequency factor for this distribution must be computed using the skewness coefficient C s . Equation (11) calculates the skewness coefficient: K T values can be obtained from tables in many sources of hydrologic data; such as [50]. By knowing the skewness coefficient and the recurrence interval, the frequency factor K T for the LPT III distribution can be extracted. The antilog of the solution in Equation (8) represents the estimated extreme value for the given return period.

Derivation of IDF Empirical Formula
The relationship between the rainfall intensity (I), rainfall duration (t), and return period (T) is described by the IDF empirical formula. Both the Gumbel distribution method and the Log-Pearson Type III method can use the formula. To derive an equation for calculating the rainfall intensity (I) for the study area, there are several steps that must be completed in order to establish an equation that is suitable for the calculation of rainfall intensity for a specific recurrence interval and specific rainfall period, which is dependent primarily on the results obtained from the IDF curves. The steps below are used to derive the equation: where I is the intensity of rainfall (mm/h), t is the duration of rainfall (minutes), T is the return period (years), and the constants (C, m, and a) are empirical parameters that are dependent on precipitation data, and the shape, size, and location of the study area, and can be obtained by using the logarithmic transformation of Equation (12). After performing a logarithmic transformation on Equation (12), we can obtain Equation (13): Assuming (CT m = K), Equation (13) can be rewritten as in Equation (14): Plotting the logarithm of precipitation intensity (logI) value against the logarithm of time (logt) value will present a linear relationship. To find the linear relationship's slope (the average value of constant a for all return periods), solve Equation (13), where represents the constant a, and log(K) in Equation (14) represents the intercept from each return period's plot. The linear relationship's slope will be displayed.
Finding C and m is as simple as plotting the logarithmic intercept (logCT m ) against the logarithmic return period (logT) in a new graph in Equation (16).
Assuming that: taking the log of both sides of Equation (15) will result in the following Equation (16): Plotting logK and logT in Equation (16), a linear equation of the plot can be obtained, then the m value can be found, which represents the slope of the linear relationship. The value of the anti-log of the intercept from the plotted curve represents the C coefficient for Equation (12).

Goodness of Fit
With this test, we can determine how well the observed frequency of occurrence in a sample corresponds to the expected frequencies obtained from the hypothesized distributions in a sample. It is possible to perform a goodness-of-fit test between observed and expected frequencies using the chi-square quantity, which is expressed as: where χ 2 is a random variable with a chi-square sampling distribution. The observed and expected frequencies of class intervals are shown in the histogram by O i and E i . The number of class intervals is k. The χ 2 value is small if the observed frequencies match the expected frequencies; otherwise, it is large. A good fit accepts the null hypothesis, while a bad fit reject it. The critical region will thus be in the chi-square distribution's right tail. The critical value is found using readily available chi-square tables, and the critical region is χ 2 > [20,22,41].

Results and Discussion
The IDF curves were developed using 39 years of rainfall intensity records for the city center of Erbil in the KRI. The IDF curves are utilized to assist in the design of drainage structures for any type of engineering project. The curves enable the engineers from local authorities to design flood control measures that are both safe and economical. The two techniques (Gumbel and LPT III) were used to estimate rainfall in mm and their intensities in millimeters per hour (mm/h) for a variety of return periods and different durations. The results of the Gumbel method are listed in Table 2. It has been demonstrated in all of them that the rainfall intensity decreased with an increase in the rainfall duration at the same return period, and that the rainfall intensity increased with an increase in the return period. Figures 2 and 3 show the IDF curves using the Gumbel method on an ordinary scale and on the log-log scale (for T = 2, 5, 10, 25, 50, 100 years). The results of this study are compared to the findings of previous work. The results revealed that the intensity of rainfall determined using the Gumbel method for different return periods in Erbil is lower than that determined in Sulaymaniyah and higher than that calculated in Baghdad and Basrah [26,29,30]. As an example, for a two-year return period and a duration of 60 min, the rainfall intensity in Erbil was 15.583 mm/h, while in Sulaymaniyah, which is located in the east of the study area and higher in elevation, the rainfall intensity was 20.21 mm/h, but in Baghdad, in the middle of Iraq, the amount decreased to about 11 mm/h and in Basrah decreased even further.    A possible explanation for these results is because of the different climates of the mentioned locations, and their elevation above sea level. In fact, the average annual rainfall in Iraq is 154 mm, but it varies widely across the country, from less than 100 mm in the south to more than 1200 mm in the northeast. The rainy season is limited to the months of October to April. Salman et al. [51] divide Iraq into four zones based on topography and climate ( Figure 4): (1) Zone-I encompasses the mountainous region in the north and northeast (21% of the country) with a Mediterranean climate; (2) Zone-II encompasses the undulating land in the south and west (9.6%) with a steppe climate; (3) Zone-III encompasses the western Plateau in the west, which is an extension of the Syrian and Arabian deserts (39.2%) having a continental desert climate; and (4) Zone-IV covers the great Mesopotamian plain in the center and the south (30.2%) having a subtropical semi-arid climate. Moreover, Al-Ansari [52] indicated that Iraq is shaped like a basin on the topographic map, containing the great Mesopotamian plains of the Tigris and Euphrates rivers. The climate is predominantly continental, subtropical, and semi-arid, with a Mediterranean climate in the north and northeastern mountainous regions. As is obvious, Iraq's climate divisions have an effect on the annual rainfall in each province. For example, the annual rainfall in Erbil, the study area is approximately 419.2 mm, while in Sulaymaniyah, the average annual precipitation varies between 328 mm in dry years and 848 mm in wet years [26]; both provinces are located in Zone I, which is considered to have a Mediterranean climate. Baghdad, in the middle of Iraq, gains an annual rainfall depth of around 280 mm [28]. The average annual precipitation in Nasiriyah Province's central district from 1965 to 2015 was 117.64 [51]. The mean annual precipitation in the central district of Basrah Province from 1965 to 2015 was 131.65 mm. Those respectively mentioned provinces are located in Zone IV, which is considered to have a subtropical, semi-arid climate [51]. Al-Ansari [53] classified Iraq into four regions based on the land's topography (Mountain Region, Plateau and Hills Regions, the Mesopotamian plain, and Jazera and Western Plateau). Erbil city center could be defined as being in the Plateau and Hills Regions; however, the northern and north eastern parts of Erbil Province are considered as the Mountain Region according to the Al-Ansari [53] classification. Figure 5 shows the elevation above sea level (ASL) and the administrative boundaries of the Iraqi Provinces. A possible explanation for these results is because of the different climates of the mentioned locations, and their elevation above sea level. In fact, the average annual rainfall in Iraq is 154 mm, but it varies widely across the country, from less than 100 mm in the south to more than 1200 mm in the northeast. The rainy season is limited to the months of October to April. Salman et al. [51] divide Iraq into four zones based on topography and climate ( Figure 4): (1) Zone-I encompasses the mountainous region in the north and northeast (21% of the country) with a Mediterranean climate; (2) Zone-II encompasses the undulating land in the south and west (9.6%) with a steppe climate; (3) Zone-III encompasses the western Plateau in the west, which is an extension of the Syrian and Arabian deserts (39.2%) having a continental desert climate; and (4) Zone-IV covers the great Mesopotamian plain in the center and the south (30.2%) having a subtropical semi-arid climate. Moreover, Al-Ansari [52] indicated that Iraq is shaped like a basin on the topographic map, containing the great Mesopotamian plains of the Tigris and Euphrates rivers. The climate is predominantly continental, subtropical, and semi-arid, with a Mediterranean climate in the north and northeastern mountainous regions. As is obvious, Iraq's climate divisions have an effect on the annual rainfall in each province. For example, the annual rainfall in Erbil, the study area is approximately 419.2 mm, while in Sulaymaniyah, the average annual precipitation varies between 328 mm in dry years and 848 mm in wet years [26]; both provinces are located in Zone I, which is considered to have a Mediterranean climate. Baghdad, in the middle of Iraq, gains an annual rainfall depth of around 280 mm [28]. The average annual precipitation in Nasiriyah Province's central district from 1965 to 2015 was 117.64 [51]. The mean annual precipitation in the central district of Basrah Province from 1965 to 2015 was 131.65 mm. Those respectively mentioned provinces are located in Zone IV, which is considered to have a subtropical, semi-arid climate [51]. Al-Ansari [53] classified Iraq into four regions based on the land's topography (Mountain Region, Plateau and Hills Regions, the Mesopotamian plain, and Jazera and Western Plateau). Erbil city center could be defined as being in the Plateau and Hills Regions; however, the northern and north eastern parts of Erbil Province are considered as the Mountain Region according to the Al-Ansari [53] classification. Figure 5 shows the elevation above sea level (ASL) and the administrative boundaries of the Iraqi Provinces.
The results of the LPT III method are listed in Table 3. Additionally, the results of this method, which were similar to those obtained from the Gumbel method, revealed that the rainfall intensity decreased with an increase in the rainfall duration in the same return period, and that the rainfall intensity increased with an increase in the return period. It was demonstrated that there were minor discrepancies in the results obtained using the two methods. The Gumbel method produces slightly higher results than LPT III for short return periods such as two, five, and ten years. While the situation becomes contradictory when the return period is increased to 25, 50, or 100 years, as the LPT III method produces slightly higher results than the Gumbel method. However, the LPT III method is believed to be more accurate because it takes into account the skewness of the data distribution.   The results of the LPT III method are listed in Table 3. Additionally, the results of this method, which were similar to those obtained from the Gumbel method, revealed that the rainfall intensity decreased with an increase in the rainfall duration in the same return period, and that the rainfall intensity increased with an increase in the return period. It was demonstrated that there were minor discrepancies in the results obtained using the two methods. The Gumbel method produces slightly higher results than LPT III for short return periods such as two, five, and ten years. While the situation becomes contradictory when the return period is increased to 25, 50, or 100 years, as the LPT III method produces Al Hassoun [21] recommended that when accuracy is required in storm design studies, the results of the LPT III method be used. Figures 6 and 7 show the intensity-durationfrequency (IDF) curves using the Gumbel method on an ordinary scale and on the log-log scale (for T = 2, 5, 10, 25, 50, 100 years). The findings of that study are consistent with the findings of previous studies conducted in some parts of the study area and neighboring countries, according to the researchers.      The parameter values obtained from analyzing the IDF data using the two methods are shown in Table 4, and these values are used to derive formulas for the study area. In addition, the IDF equation parameters (C, m and a) obtained from the two distributions can be considered to be relatively close, with the exception of parameter C from the Gumbel Theory distribution, which was slightly higher than in the other distribution (see Table 4). Moreover, the IDF equation parameters (C, m and a) of other Iraqi Provinces have been shown, revealing that the (m and a) parameters are slightly different between the developed IDF equations of all the mentioned governorates, while the parameter (C) totally varies. Sulaymaniyah and Erbil have the highest (C) values of 288.4 and 234.020, respectively, using the Gumbel method. While in Baghdad, Nasiriyah, and Basrah, the value of the (C) parameter decreased to 9.074, 8.283, and 7.95, respectively. This variation in parameter values in different Iraqi cities is evidence of the different climate conditions in each city, as explained earlier in this section. During the development of the IDF empirical formulas, a logarithmic conversion was applied to the precipitation parameters for both methods, which allowed the relationship to be converted into a linear equation, which in turn allowed the equation to be used to calculate all of the parameters associated with it. The parameters of the empirical equations (Equations (18) and (19)), which were determined for the Gumbel and LPT III distribution methods, respectively, are presented in Table 4.
Furthermore, goodness-of-fit tests were applied to determine which statistical distribution was the best fit among the various techniques. Table 5 presents the results of the chi-square goodness-of-fit test performed on an annual series of rainfall data. Table 5 shows that there is no statistically significant difference between the rainfall intensity values for all return periods, based on the results of the chi-square test performed on the dataset. This finding validates the usefulness of local IDF formulas instead of one general equation for the whole country because, as mentioned previously, the rainfall is unpredictable and erratic in both time and space.

Conclusions
In this study, IDF curves and IDF empirical formulas were derived for rainfall durations (10, 20, 30, 60, 120, 180, 360, 720, and 1440 min) for return periods (2, 5, 10, 25, 50, and 100 years) for Erbil city in the north of Iraq. The IDF curves were obtained using two different methods (Gumbel and LPT III). Although, the suitability of these distribution types has not been determined in this analysis, in order to maintain consistency with previous studies, we decided to use the same methodology as other researchers. The rainfall IDF curves developed in this study are immensely useful in the estimation of rainfall intensity for the design of flood control structures and water resources.
In this research, it was discovered that there is no significant difference between the results of the Gumbel and LPT III methods. Results in this study show that the IDF curves for Erbil city are fully fit to predict the values of rainfall intensities after subjecting the results to calibration and checking the empirical equation with the coefficient of determination (R 2 = 1) and the accepted chi-square values. In addition, the derivation of IDF formulas is useful to make an approximate prediction more easily in comparison to IDF curves. The findings of the study are consistent with the findings of previous studies conducted in some parts of the study area. The findings confirmed that each province must have its own set of IDF curves because Iraqi cities differ in terms of topography and climate. However, despite the promising results, there are still some issues. The study's validity could not be checked because there were not any recorded data on shorter-duration rainfall events. Analysis only examines the statistical relationship between the variables and does not reveal their physical mechanism.
The findings of this study can be used effectively as a guideline for urban water resource management, particularly in the design of storm water systems. As far as we are aware, no such analysis exists in the studied area. In fact, it would be beneficial to validate the results of this study using some new techniques, such as satellite-based precipitation, in order to confirm their accuracy. Additionally, the study's findings and developed IDF curves will be useful in the field of research on flash flood mitigation in the studied area, as well as in the preparation of different flood risk maps and scenarios.
Author Contributions: D.A.K., Conceptualization, methodology, writing-original draft preparation; A.R.M.A., conceptualization, methodology, formal analysis; A.M., data curation, investigation, writing-reviewing and editing; M.I.Y., supervision; M.S., writing-review and editing, validation, supervision. All authors have read and agreed to the published version of the manuscript.