Developing Intensity-Duration-Frequency (IDF) Curves Based on Rainfall Cumulative Distribution Frequency (CDF) for Can Tho City, Vietnam

: Information on the relationship between rainfall intensity, duration and accumulation frequency or return period (IDF) is commonly utilized in the design and management of urban drainage systems. Can Tho City, located in the Vietnamese Mekong Delta, is a city which has recently invested heavily in upgrading its stormwater drainage systems in the hope of preventing reoccurring ﬂood events. Yet, much of these works were designed based on obsolete and outdated IDF rainfall curves. This paper presents an updated IDF curve for design rainfall for Can Tho City. For each duration and designated return period, a cumulative distribution function (CDF) was developed using the Pearson III, Log-Pearson III, and Log-Normal distribution functions. In order to choose the best IDF rainfall curve for Can Tho City, the CDF rainfall curve and empirical formulas used in Vietnam and Asia (Vietnamese standard 7957:2008, Department of Hydrology, Ministry of Transportation, Talbot, Kimijima, and Bermard) were compared. The goodness of ﬁt between the IDF relationship generated by the frequency analysis (CDF curve), and that predicted by the IDF empirical formulas was assessed using the efﬁciency index ( EI ), and the root mean squared error ( RMSE ). The IDF built from Vietnam’s standard TCVN 7957:2008 with new parameters ( A = 9594, C = 0.5, b = 26, n = 0.96) showed the best performance, with the highest values of EI (0.84 ≤ EI ≤ 0.93) and the lowest values of RMSE (2.5 ≤ RMSE ≤ 3.2), when compared to the other remnants.


Introduction
Determining the three key variables of rainfall intensity (I), frequency (F), and duration (D) are important prerequisites in the design of urban drainage systems and stormwater management practices [1,2]. The relationship between the three variables is usually represented as a graph, in which rainfall duration is placed on the horizontal axis, rainfall intensity is placed on the vertical axis, and the curves corresponds to the period of reappearance or frequency. The relationship has been a commonly applied input criteria for the design of drainage systems worldwide since the 1920s. Urbanization and floodplain management were of interest to Meyer [3]. Sherman [4] used the basic formula to determine rainfall intensity for different areas throughout each period, and Bernard (1932) [5] further developed the IDF curve so that it could be applied to larger study regions. Thereafter, investigating varying return periods of 2, 5, 20, 50, 100, and 200 years and varying rainfall intensity durations, Hershfield (1961) [6] created numerous IDF rainfall curves. Whereas, Earth 2022, 3 867 Bell created an IDF rainfall curve in 1969 using a scale zoom calculation formula for the then Soviet Union, and Chen (1983) [7] presented a straightforward and generic method to create an IDF curve that was applicable to all study areas.
In a more recent study, Di Baldassarre et al. (2006) analyzed seven individual IDF curves to determine how well they represented the varying relationships between rainfall intensity, duration, and frequency. They found that rainfall intensity could only be accurately determined for short rainfall events (typically less than 1 day). However, if the rainfall monitoring data at a single station was sufficiently large, cumulative distribution frequency (CDF) could be used successfully to generate an IDF rainfall curve [8,9].
Similar research applied the IDF curve for precipitation data in Korea using the cumulative distribution function [10]. Here, the authors estimated the IDF curve parameters from the CDF using a variety of techniques, including single-and multi-objective algorithms, ROBUST, OneStep, and linear and non-linear regression (MOGA). They then roughly divided the parameters into short-term and long-term precipitation using the SEPDUR approach to produce a more precise calculation. The effectiveness of each method was compared using both the RMSE and RRMSE between the estimated value and the findings from an in-situ frequency study. The results demonstrated that, the IDF curve based on the CDF showed the least relative standard error in the IDF (RRMSE) and relative mean square error (RMSE), respectively. These results found that the most accurate method is a MOGA's IDF curve, which links CDF and SEPDUR.
The work investigating the design rainfall for Qatar by [11], employed a linear moments (L-monents) regional index to develop input data for an IDF curve. With an average relative error of only 5.5%, the study also highlighted that employing the Pearson III distribution function to construct the IDF rainfall curve was optimal. An IDF curve based on the Gumbel and Person III functions was also developed in a study focusing on Saudi Arabia [12]. Here, these methods both provided outcomes that were consistent with previous research. The parameters of the IDF equation were calculated for different return periods using non-linear multivariable regression with correlation numerical analysis for 2, 5, 10, 25, 50, and 100 years. The obtained results showed that every return period exhibited a strong correlation and were suitable for use in the computation of the study area's IDF curve. The work also highlighted that the IDF curve's construction based on the Gumbel distribution function displayed a high level of accuracy.
It is now well-known and acknowledged that climate change exerts a direct impact on the hydrological cycle, influencing precipitation intensity, duration, and frequency [13]. Rainfall is frequently taken into account for planning and designing stormwater collection systems, run-off, flood estimation, calculating urban drainage capacities, and computing the effects of different climate change scenarios in order to estimate the future risk and the efficiency of different adaptation options in urban areas.
However, because each location exhibits a unique I, D, F relationship, it is still challenging to evaluate the link between rainfall intensity for specific periods (I), duration of rain events (D), and the frequency of precipitation (F), or repeatability (N). Currently in Asia, IDF rainfall curves have been developed, for monsoon-affected regions with high rainfall, such as China (Sherman), Indonesia (Talbot, Sherman), and Japan (Kimijima, Talbot, Bermard), in addition to the application of some methods for calculating rain intensity according to the limit intensity of the Soviet Union (former) D.F. Gorbachev (1920) or British, American, German, and Polish formulas. In Vietnam, the Department of Hydrology, Ministry of Transport, and Vietnamese standards are used to develop experimental formulations (TCVN 7957:2008). The deterministic character of the hydrological system (geography) provides the foundation for the aforementioned approaches, which ignores the randomness of hydrological occurrences. Using the case of Can Tho City in the Vietnamese Mekong delta, this study applies experimental methods and probability distribution, which have been used successfully in many previous studies [14][15][16] to construct a locally relevant IDF curve.

Materials and Methods
Some definitions of rainfall intensity, design rain time step, and repetition period (repeatability) are provided in this study.

1.
Rainfall intensity (i, q), is the amount of rain falling per unit area in a unit of time.
Rainfall intensity can be in the form of accumulative rainfall depth (i.e., mm/min, cm/h) or rainfall volume (q, L/s.ha) and is shown to be related by the formula q = 166.7i.

2.
Design rainfall time step (period, t), is the amount of rainfall in the calculation period (15 min or 30 min) assuming that the intensity does not change.

3.
Rainfall return period or frequency (N, years): is the time in years for the reappearance of rainfall of the equal amount (or greater than that of design intensity). The larger the rainfall, the less frequent it will occur (in case P ≤ 50%).
Figures 1 and 2 depict the location of Can Tho, a city in the VMD, and outline of the research methodology used in this study respectively. The study collected rainfall data from the official Can Tho station, Can Tho City over a period of 38 years . For each year, the maximum precipitation value was selected for the calculation of time steps (15 min, 30 min, 60 min, 90 min, and 180 min respectively). distribution, which have been used successfully in many previous studies [14][15][16] to construct a locally relevant IDF curve.

Materials and Methods
Some definitions of rainfall intensity, design rain time step, and repetition period (repeatability) are provided in this study.
1. Rainfall intensity (i, q), is the amount of rain falling per unit area in a unit of time.
Rainfall intensity can be in the form of accumulative rainfall depth (i.e., mm/min, cm/h) or rainfall volume (q, L/s.ha) and is shown to be related by the formula q = 166.7i. 2. Design rainfall time step (period, t), is the amount of rainfall in the calculation period (15 min or 30 min) assuming that the intensity does not change. 3. Rainfall return period or frequency (N, years): is the time in years for the reappearance of rainfall of the equal amount (or greater than that of design intensity). The larger the rainfall, the less frequent it will occur (in case P ≤ 50%).     [17]. According to the research results of Tuan [18], who compared the frequency distribution functions for rainfall in the VMD, the distribution functions of Pearson III (P.III) and Log-Pearson III, Log -Normal are the most suitable. Furthermore the research of Anh [19] shows that the Log-Normal, Pearson III and Log Pearson III all provided equivalent calculation results and large biased result compared to other curves in the region with rare return periods for curve for Hoa Binh and Son Tay stations in Vietnam. Resultantly, this study assessed the development a cumulative rain frequency curve for each period according to the distribution of Pearson III (P.III), Log-Pearson III, Log-Normal and re-adjusted the coefficient of variation (Cv) suitable for rainfall in each period.
In which: α, β, γ are the 3 parameters which indicate the scale, shape and position of the data set, respectively; while Γ(β) is the Gamma function.
Determined rainfall intensity return period (N, year): Determine the rainfall intensity return period (N, years) for the rainfall intensity for the selected event 2 years, 10 years, 20 years, 50 years, and 100 years, respectively. Based on the event's occurrence frequency, determined the return period: N = 100/P (years) if P ≤ 50%.
Develop the IDF rainfall curve from the CDF rainfall curve: The rainfall intensity according to the return period can be determined from the curve of cumulative duration frequency (CDF) for each duration of rainfall. Figure 3 offers an illustration of how to Some European and Southeast Asian countries have typically used 3-parameter Log-Normal distribution or External I and External III. While, Russia, China and Vietnam have typically relied Pearson III distribution [17]. According to the research results of Tuan [18], who compared the frequency distribution functions for rainfall in the VMD, the distribution functions of Pearson III (P.III) and Log-Pearson III, Log-Normal are the most suitable. Furthermore the research of Anh [19] shows that the Log-Normal, Pearson III and Log Pearson III all provided equivalent calculation results and large biased result compared to other curves in the region with rare return periods for curve for Hoa Binh and Son Tay stations in Vietnam. Resultantly, this study assessed the development a cumulative rain frequency curve for each period according to the distribution of Pearson III (P.III), Log-Pearson III, Log-Normal and re-adjusted the coefficient of variation (C v ) suitable for rainfall in each period.
Distribution of Log-Pearson III : Distribution of Log-Normal : In which: α, β, γ are the 3 parameters which indicate the scale, shape and position of the data set, respectively; while Γ(β) is the Gamma function.
Determined rainfall intensity return period (N, year): Determine the rainfall intensity return period (N, years) for the rainfall intensity for the selected event 2 years, 10 years, 20 years, 50 years, and 100 years, respectively. Based on the event's occurrence frequency, determined the return period: N = 100/P (years) if P ≤ 50%.
Develop the IDF rainfall curve from the CDF rainfall curve: The rainfall intensity according to the return period can be determined from the curve of cumulative duration Earth 2022, 3 870 frequency (CDF) for each duration of rainfall. Figure 3 offers an illustration of how to calculate the intensity of rain that will fall over a period of 60 min, corresponding to a return period of 100 years and 2 years, respectively. Using this method, the duration and return periods for other rainfall events were determined.
Earth 2022, 3, FOR PEER REVIEW 5 calculate the intensity of rain that will fall over a period of 60 min, corresponding to a return period of 100 years and 2 years, respectively. Using this method, the duration and return periods for other rainfall events were determined.

Development of the IDF Rainfall Curve Using Empirical Formula
The study developed the IDF rainfall curve from empirical formulas according to (i) the Vietnamese standard (TCVN 7957: 2008) [20], (ii) the according to the Department of Hydrology, and (iii) the Ministry of Transport.
To develop IDF rainfall curve according to the Vietnamese standard (TCVN 7957:2008) To develop IDF rainfall curve according to the Department of Hydrology To develop IDF rainfall curve according to the Ministry of Transport where: q-rainfall intensification ( . ℎ ); t, Td-rainfall duration (minutes); P, f (N)-rainfall event return period (year); A, C, b, n, q20-These parameters were determined according to local conditions ((Vietnamese standard TCVN 7975:2008; for areas not listed in the table, adjacent or the closest areas were used) [20]. In addition, the study was also expanded to include the development of some empirical IDF rainfall curves commonly used in Japan (Kimijima, Bermard, Talbot).
Bermard: = Talbot: = where: i-Rainfall intensity; d-Rainfall duration; a, b, c-These available parameters are listed in tables.

Development of the IDF Rainfall Curve Using Empirical Formula
The study developed the IDF rainfall curve from empirical formulas according to (i) the Vietnamese standard (TCVN 7957: 2008) [20], (ii) the according to the Department of Hydrology, and (iii) the Ministry of Transport.
To develop IDF rainfall curve according to the Vietnamese standard (TCVN 7957:2008) To develop IDF rainfall curve according to the Department of Hydrology To develop IDF rainfall curve according to the Ministry of Transport where: q-rainfall intensification ( L s .ha); t, T d -rainfall duration (minutes); P, f (N)-rainfall event return period (year); A, C, b, n, q 20 -These parameters were determined according to local conditions ((Vietnamese standard TCVN 7975:2008; for areas not listed in the table, adjacent or the closest areas were used) [20].
In addition, the study was also expanded to include the development of some empirical IDF rainfall curves commonly used in Japan (Kimijima, Bermard, Talbot). where: i-Rainfall intensity; d-Rainfall duration; a, b, c-These available parameters are listed in tables.

Selecting the Most Suitable IDF Rainfall Curve for Can Tho City
To select the most suitable IDF rainfall curve for Can Tho City, the IDF curves built from the empirical formulas were compared with the IDF curve built from the CDF rainfall curve. The study used the EI index, and the root mean standard error (RMSE) for comparison and selection.

Various Rainfall Intensity Duration Assessments
Rainfall measurements collected over 38 years included the long-term rainfall duration (yearly data), and the maximum rainfall intensify for short periods (180 min, 30 min, 60 min, 90 min, 120 min, and 180 min, respectively) as shown in Figure 4.

Selecting the Most Suitable IDF Rainfall Curve for Can Tho City
To select the most suitable IDF rainfall curve for Can Tho City, the IDF curves built from the empirical formulas were compared with the IDF curve built from the CDF rainfall curve. The study used the EI index, and the root mean standard error (RMSE) for comparison and selection.

Various Rainfall Intensity Duration Assessments
Rainfall measurements collected over 38 years included the long-term rainfall duration (yearly data), and the maximum rainfall intensify for short periods (180 min, 30 min, 60 min, 90 min, 120 min, and 180 min, respectively) as shown in Figure 4.  Figure 4 shows that the rainfall in the rainy season and the dry season both increased and decreased alternately over the years. Moreover, the coefficient of variation (Cv) for all periods was quite high ( ≥ 0.2), yet the fluctuations in the periods were found to be not significantly different from each other (Cv ranges from 0.2 to 0.3) ( Table  1). The high Cv showed that the level of precipitation varied from period to period over the years. Although Cv was quite high, the data series of 38 years for Can Tho station is quite long, so ̅ ≤ 6% in all time periods. According to the Code of Calculation of Design Hydrological Characteristics (1997), the 38-year data series is reliable enough to represent the Can Tho station for calculating the cumulative rainfall frequency distribution. On the one hand, Figure 5 shows rainfall in the rainy season has tended to decrease, while on the other hand, the rainfall in the dry season has tended to increase. Rainfall in each period (15 min, 30 min, 60 min, 90 min, 120 min, and 180 min) has also tended to increase. Furthermore, rainfall over shorter periods has tended to increase, more than rainfall over longer periods. However, these uptrends showed low R-square values.   Figure 4 shows that the rainfall in the rainy season and the dry season both increased and decreased alternately over the years. Moreover, the coefficient of variation (C v ) for all periods was quite high (C v ≥ 0.2), yet the fluctuations in the periods were found to be not significantly different from each other (C v ranges from 0.2 to 0.3) ( Table 1). The high C v showed that the level of precipitation varied from period to period over the years. Although C v was quite high, the data series of 38 years for Can Tho station is quite long, so σ x ≤ 6% in all time periods. According to the Code of Calculation of Design Hydrological Characteristics (1997), the 38-year data series is reliable enough to represent the Can Tho station for calculating the cumulative rainfall frequency distribution. On the one hand, Figure 5 shows rainfall in the rainy season has tended to decrease, while on the other hand, the rainfall in the dry season has tended to increase. Rainfall in each period (15 min, 30 min, 60 min, 90 min, 120 min, and 180 min) has also tended to increase. Furthermore, rainfall over shorter periods has tended to increase, more than rainfall over longer periods. However, these uptrends showed low R-square values.

Rainfall Duration (Min)
C

Developing an IDF Rainfall Curve Based on CDF Rainfall Curve
The CDF was built from the three distribution curves, P.III, Log-P.III, and Log-Normal 3 parameters in a short period of time (<1 day). The Log-P.III distribution pattern showed that it was only suitable for the rainfall measured in the 60-min duration and 180-min duration, while the remaining durations were deemed not suitable with the measured data obtained from Can Tho station. While P.III, Log-Normal III were suitable

Developing an IDF Rainfall Curve Based on CDF Rainfall Curve
The CDF was built from the three distribution curves, P.III, Log-P.III, and Log-Normal 3 parameters in a short period of time (<1 day). The Log-P.III distribution pattern showed that it was only suitable for the rainfall measured in the 60-min duration and 180-min duration, while the remaining durations were deemed not suitable with the measured data obtained from Can Tho station. While P.III, Log-Normal III were suitable for the rain durations (15 min, 30 min, 60 min, 90 min, and 180 min). The study calculated and compared the RMSE between P.III and Log-N-III, with the results suggesting that the P.III curve is more consistent with the actual measured points at each calculation period than the with Log-N.III three parameters. Therefore, the study selected P.III to build the CDF curve for each period in this study. Figure A1 shows the distribution of P.III and Log-N-III for the rain monitoring points of each period at Can Tho station.
Based on each calculation period and repeatability, the study found the corresponding rainfall intensity (mm/h) shown in Figure 6.
Earth 2022, 3, FOR PEER REVIEW 8 the CDF curve for each period in this study. Figure A1 shows the distribution of P.III and Log-N-III for the rain monitoring points of each period at Can Tho station. Based on each calculation period and repeatability, the study found the corresponding rainfall intensity (mm/h) shown in Figure 6.

Development of IDF Rainfall Curve from Empirical Formulas According to the Vietnamese Standard (TCVN 7957:2008), the Department of Hydrology, and the Ministry of Transport
According to the Vietnamese design standards (TCVN 7957:2008) on drainage-external networks and structures, to determine the intensity of rainfall in each return period, it is possible to use a curve or other empirical calculation formulas for comparison. Therefore, we attempted this with two experimental formulas proposed by TCVN 7957:2008. Climatic parameters were not shown for the study area (Can Tho City) in the table according to TCVN 7957:2008. Therefore, the study compared the data of two neighboring stations (stations in Ho Chi Minh City and Ca Mau). The results show that the parameters in Ca Mau were best suited to the study area, as shown in Figure 7.

Development of IDF Rainfall Curve from Empirical Formulas According to the Vietnamese Standard (TCVN 7957:2008), the Department of Hydrology, and the Ministry of Transport
According to the Vietnamese design standards (TCVN 7957:2008) on drainage-external networks and structures, to determine the intensity of rainfall in each return period, it is possible to use a curve or other empirical calculation formulas for comparison. Therefore, we attempted this with two experimental formulas proposed by TCVN 7957:2008. Climatic parameters were not shown for the study area (Can Tho City) in the table according to TCVN 7957:2008. Therefore, the study compared the data of two neighboring stations (stations in Ho Chi Minh City and Ca Mau). The results show that the parameters in Ca Mau were best suited to the study area, as shown in Figure 7.
Earth 2022, 3, FOR PEER REVIEW 8 the CDF curve for each period in this study. Figure A1 shows the distribution of P.III and Log-N-III for the rain monitoring points of each period at Can Tho station. Based on each calculation period and repeatability, the study found the corresponding rainfall intensity (mm/h) shown in Figure 6.

Development of IDF Rainfall Curve from Empirical Formulas According to the Vietnamese Standard (TCVN 7957:2008), the Department of Hydrology, and the Ministry of Transport
According to the Vietnamese design standards (TCVN 7957:2008) on drainage-external networks and structures, to determine the intensity of rainfall in each return period, it is possible to use a curve or other empirical calculation formulas for comparison. Therefore, we attempted this with two experimental formulas proposed by TCVN 7957:2008. Climatic parameters were not shown for the study area (Can Tho City) in the table according to TCVN 7957:2008. Therefore, the study compared the data of two neighboring stations (stations in Ho Chi Minh City and Ca Mau). The results show that the parameters in Ca Mau were best suited to the study area, as shown in Figure 7.   Similarly, the IDF rainfall curve built on the empirical formula of the Department of Hydrology and the Ministry of Transport is shown in Figure 8.
Earth 2022, 3, FOR PEER REVIEW 9 Similarly, the IDF rainfall curve built on the empirical formula of the Department of Hydrology and the Ministry of Transport is shown in Figure 8.

The IDF Rainfall Built from Empirical Formulas That Have Been Used in Some Asian Countries
Furthermore, the results of IDF curves according to Talbot, Kimijima and Bermard are shown in Figure 9.

The IDF Rainfall Built from Empirical Formulas That Have Been Used in Some Asian Countries
Furthermore, the results of IDF curves according to Talbot, Kimijima and Bermard are shown in Figure 9.

Compareration the Selection of IDF Curves for the Study Area
The study compared the IDF rainfall curve built from CDF with empirical formulas. The results shown in Figure 10 are efficiency factor (EI) comparisons.

Compareration the Selection of IDF Curves for the Study Area
The study compared the IDF rainfall curve built from CDF with empirical formulas. The results shown in Figure 10 are efficiency factor (EI) comparisons.
(a) (b) Figure 10. Comparation of the EI index (a) and RMSE index (b) among empirical formulas. Figure 9 shows the three empirical formulas, including those from the Department of Hydrology, the Vietnamese national Standard (TCVN 7957:2008), and the Ministry of Transport which have EI ≥ 0.6. The two datasets calculated according to the Department of Hydrology and TCVN 7957:2008 have EI ≥ 0.9, although they show slight decrease (EI ≈ 0.8) at N = 2 years. In general, the RMSE is within the acceptable range (RMSE ≈ 10). Meanwhile, at N = 2 years, N = 5 years, and N = 10 years, the practical formula according to the Ministry of Transport is the best (EI ≥ 0.9). However, as N increases, the EI decreases and the RMSE increases (the highest RMSE = 35 at N = 100 years). The formula according to TCVN 7957:2008 has good EI values at all return periods (EI ≥ 0.8) and RMSE ≈ 12 is quite high. The study used the trial and error method to find the suitable coefficient for the study area (A = 9594, C = 0.5, b = 26, n = 0.96). The results of comparing EI and RMSE between the data using the neighbor approach and adjustment coefficients for the study area are shown in Figure 11.   Figure 9 shows the three empirical formulas, including those from the Department of Hydrology, the Vietnamese national Standard (TCVN 7957:2008), and the Ministry of Transport which have EI ≥ 0.6. The two datasets calculated according to the Department of Hydrology and TCVN 7957:2008 have EI ≥ 0.9, although they show slight decrease (EI ≈ 0.8) at N = 2 years. In general, the RMSE is within the acceptable range (RMSE ≈ 10). Meanwhile, at N = 2 years, N = 5 years, and N = 10 years, the practical formula according to the Ministry of Transport is the best (EI ≥ 0.9). However, as N increases, the EI decreases and the RMSE increases (the highest RMSE = 35 at N = 100 years). The formula according to TCVN 7957:2008 has good EI values at all return periods (EI ≥ 0.8) and RMSE ≈ 12 is quite high. The study used the trial and error method to find the suitable coefficient for the study area (A = 9594, C = 0.5, b = 26, n = 0.96). The results of comparing EI and RMSE between the data using the neighbor approach and adjustment coefficients for the study area are shown in Figure 11.

Compareration the Selection of IDF Curves for the Study Area
The study compared the IDF rainfall curve built from CDF with empirical formulas. The results shown in Figure 10 are efficiency factor (EI) comparisons.
(a) (b) Figure 10. Comparation of the EI index (a) and RMSE index (b) among empirical formulas. Figure 9 shows the three empirical formulas, including those from the Department of Hydrology, the Vietnamese national Standard (TCVN 7957:2008), and the Ministry of Transport which have EI ≥ 0.6. The two datasets calculated according to the Department of Hydrology and TCVN 7957:2008 have EI ≥ 0.9, although they show slight decrease (EI ≈ 0.8) at N = 2 years. In general, the RMSE is within the acceptable range (RMSE ≈ 10). Meanwhile, at N = 2 years, N = 5 years, and N = 10 years, the practical formula according to the Ministry of Transport is the best (EI ≥ 0.9). However, as N increases, the EI decreases and the RMSE increases (the highest RMSE = 35 at N = 100 years). The formula according to TCVN 7957:2008 has good EI values at all return periods (EI ≥ 0.8) and RMSE ≈ 12 is quite high. The study used the trial and error method to find the suitable coefficient for the study area (A = 9594, C = 0.5, b = 26, n = 0.96). The results of comparing EI and RMSE between the data using the neighbor approach and adjustment coefficients for the study area are shown in Figure 11.   The comparison of the EI and RMSE according to the TCVN between the neighboring station and with the adjustment dataset for the study area show that, at each return period, EI was ≥0.97 in the adjustment dataset, while at the neighboring station, EI is only in the range of 0.84 ≤ EI ≤ 0.93. When comparing the RMSE according to the Vietnamese standard between the neighboring station and with the adjustment dataset, the adjusted dataset that RMSE is in the range of 2.5 ≤ RMSE ≤ 3.2, while the neighboring station is comparatively high, 11 ≤ RMSE ≤17. A similar study also comparing EI and RMSE for Talbot, Kimijima and Bermard are shown in Table 2.  Table 2 shows that the empirical formula according to Talbot has a high EI and lower RMSE, than the other two formulas. For the return period of 25 years or more, the EI at three stations gave good results, but the RMSE at all periods were high (RMSR ≥ 19.4). Therefore, all three stations Talbot, Kimijima, and Bermard were deemed not suitable for the study area.
From the results of the comparison of efficiency coefficients, and mean standard errors, along with comparison of IDF curves built according to empirical formulas in Vietnam and Asia, we developed an IDF rainfall curve set for Can Tho City as shown in Figure 12.
Earth 2022, 3, FOR PEER REVIEW 12 The comparison of the EI and RMSE according to the TCVN between the neighboring station and with the adjustment dataset for the study area show that, at each return period, EI was ≥0.97 in the adjustment dataset, while at the neighboring station, EI is only in the range of 0.84 ≤ EI ≤ 0.93. When comparing the RMSE according to the Vietnamese standard between the neighboring station and with the adjustment dataset, the adjusted dataset that RMSE is in the range of 2.5 ≤ RMSE ≤ 3.2, while the neighboring station is comparatively high, 11 ≤ RMSE ≤17. A similar study also comparing EI and RMSE for Talbot, Kimijima and Bermard are shown in Table 2.  Table 2 shows that the empirical formula according to Talbot has a high EI and lower RMSE, than the other two formulas. For the return period of 25 years or more, the EI at three stations gave good results, but the RMSE at all periods were high (RMSR ≥ 19.4). Therefore, all three stations Talbot, Kimijima, and Bermard were deemed not suitable for the study area.
From the results of the comparison of efficiency coefficients, and mean standard errors, along with comparison of IDF curves built according to empirical formulas in Vietnam and Asia, we developed an IDF rainfall curve set for Can Tho City as shown in Figure 12.   Figure 12 shows the proposed set of IDF curves for Can Tho City, in which at a period of less than 60 min, the monitoring points are approximately lower than the values of the proposed IDF curve. However, for durations of more 60 min, the values at the monitoring points are approximately greater than the values of the proposed IDF curve at all repetition periods. Therefore, although the results of the proposed IDF curve fit very well with the monitoring points, further studies are needed to adjust the parameters (A, C, b, n).

Conclusions
This study calculated and selected the most suitable IDF rainfall curve for Can Tho City according to the empirical formulas in Vietnam (Department of Hydrology, TCVN 7957:2008, and the Ministry of Transport), as well as, via the formulas used previously for several other Asian countries (Talbot, Kimijima, and Bermard). In order to compare and select a suitable IDF rainfall curve, we constructed an IDF rainfall curve from the CDF according to rainfall intensity for each duration (15 min, 30 min, 60 min, 90, 120 min and 180 min), and also using EI and RMSE as the basis for selection. Research results for the three experimental formulas (Talbot, Kimijima, and Bermard) were found not suitable for the study area. When comparing the empirical formulas in Vietnam, the EI was relatively high, and the RMSE was also quite large. Here, the experimental formula according to TCVN 7957:2008 gave the best results at all calculation periods, compared to the other two formulas. As there is no set coefficients for Can Tho City in TCVN 7957:2008, this study used the closest location (Ca Mau) for construction. Therefore, the study used (A = 9594, C = 0.5, b = 26, n = 0.96) according to TCVN 7957:2008. With the new set of coefficients, the developed and adjusted IDF rainfall curve fitted very well with IDF built from CDF with 0.84 ≤ EI ≤ 0.93 and 2.5 ≤ RMSE ≤ 3.2. The suitable IDF rainfall curve for Can Tho City was selected according to the empirical formula of TCVN 7957:2007, with a new developed correction factor for the city. Although the proposed IDF rainfall curve results for Can Tho City fitted well with the monitoring points, it is necessary to continue to adjust the parameters (A, C, b, n) in further studies.