Rainfall Intensity-Duration-Frequency Relationship. Case Study: Depth-Duration Ratio in a Semi-Arid Zone in Mexico

: Intensity–Duration–Frequency (IDF) curves describe the relationship between rainfall intensity, rainfall duration, and return period. They are commonly used in the design, planning and operation of hydrologic, hydraulic, and water resource systems. Considering the intense rainfall presence with flooding occurrences, limited data used to develop IDF curves, and importance to improve the IDF design for the Ensenada City in Baja California, this research study aims to investigate the use and combinations of pluviograph and daily records, to assess rain behavior around the city, and select a suitable method that provides the best results of IDF relationship, consequently updating the IDF relationship for the city for return periods of 10, 25, 50, and 100 years. The IDF relationship is determined through frequency analysis of rainfall observations. Also, annual maximum rainfall intensity for several duration and return periods has been analyzed according to the statistical distribution of Gumbel Extreme Value (GEV). Thus, Chen’s method was evaluated based on the depth-duration ratio ( R ) from the zone, and the development of the IDF relationship for the rain gauges stations was focused on estimating the most suitable ( R ) ratio; chosen from testing several methods and analyzing the rain in the region from California and Baja California. The determined values of the rain for one hour and return period of 2 years obtained were compared to the values of some cities in California and Baja California, with a range between 10 and 16.61 mm, and the values of the ( R ) ratio are in a range between 0.35 and 0.44; this range is close to the ( R ) ratio of 0.44 for one station in Tijuana, a city 100 km far from Ensenada. The values found here correspond to the rainfall characteristics of the zone; therefore, the method used in this study can be replicated to other semi-arid zones with the same rain characteristics. Finally, it is suggested that these results of the IDF relationship should be incorporated on the Norm of the State of Baja California as the recurrence update requires it upon recommendation. This study is the starting point to other studies that imply the calculation of a peak flow and evaluation of hydraulic structures as an input to help improve flood resilience in the city of Ensenada .


Introduction
Intensity-Duration-Frequency (IDF) relationship, or IDF curves, is a representation of intense rainfall events that allows for calculation of a peak flow needed to design hydraulic structures (e.g. storm sewers, culverts, drainage systems), to assess and predict flood hazard, and design flood protection structures [1][2][3]. Most of these structures were designed in many developing countries a long time ago without an updated IDF-remarkable, since rain is a variable that changes with space and time. Consequently, any update from the IDF relationship in urban catchments will necessarily imply revision and modification of the local standard structural designs [4]. The updated IDF design would also help to predict flood risk occurrences and map flood hazards of the expected peak flow. This is especially relevant in arid and semiarid zones where rain characteristics indicate that yearly variation in storms is very large, and the intensity of rare storms is always very high for a brief period, and so therefore the flooding is of sudden occurrence and rapid rise [5]. Northwest Mexico is a semiarid region with low annual average rainfall; however, with the presence of rainfall intensities associated with climate variability, that has caused flooding [6][7][8].
The Intensity-Duration-Frequency relationship design is based on the measurement of peak rainfall events, regarding duration and developed for a certain recurrence interval or return period [9], where accuracy depends on the rainfall characteristics, such as magnitude, frequency, and duration [10,11]. The analyzed data is the precipitation time series, modeled for future projections at a regional scale. These projections indicate that the precipitation return period tends to increase if the climate conditions do not change [12,13]. However, due to Climate Change, there are uncertainties of intense rainfall occurrence affecting the return period of the IDF design, which in turn tends to decrease in some global regions [14]. This issue makes it necessary for trend analysis of the extreme rainfall events to design or update the IDF relationship.
The best estimation of the precipitation intensity (unit/time) is directly obtained from the automatic (pluviograph) weather station that measures the sub-hourly rain, and for which the records are automatically transmitted every five or ten minutes [15,16]. There is a low density of this kind of weather station in the Mexican territory, though. The separation between stations should be between 5 and 30 km [17,18] according to the Manual on the Global Observing System proposed by the World Meteorological Organization. However, distance between the stations in Mexico is 70 km on average; thus, there is a lack of information between stations along space and time. In the country, and only since 1999, the available sub-hourly rain records have been registered through automatic weather stations (EMAs), by the National Meteorological Service (SMN).
The Ensenada's city has two IDF designs: one published as an Official Norm by the State of Baja California, with daily rainfall records from 1948 to 2008, just for one station [19]; and the other, by the Federal Communication and Transportation Ministry (SCT) published in 2000, and whose analyzed time series are unknown [20]. These are the only IDF studies found in the literature for the study area, and both serve as official Mexican documents [21,22]. Due to a lack of automatic stations, none of these studies used the pluviograph records from the last two decades , to assess extreme rainfall events necessary to develop the IDF curves. This problem of a lack of data and weather stations is evident in several cities in Mexico; however, to minimize the inconvenience of the periodicity of measurements and the distance between weather stations, most studies in Mexico have used different empirical methods to estimate the IDF relationship based on daily rainfall from standard rain gauges (pluviometric), suggesting that the Chen method [23] is the most appropriate for the IDF estimation [24,25].
Considering flood occurrence caused by intense rains, the insufficient studies on the IDF relationship, lack of pluviographic information, and other important aspects for planning hydraulic structures in Ensenada city and its surroundings, it is necessary to carry out a detailed analysis of the Intensity-Duration-Frequency relationship for the area. Therefore, the purpose of this study is to analyze, estimate, and propose the Depth-Duration-Ratio (R), appropriate for the characteristics of the rains that occur in the study area, to obtain representative IDF curves. For this, it is necessary to estimate the rainfall of one hour, and the return period 2 1 ( ) P of two years. 2 1 ( ) P is traditionally derived from the proposals made by Hershfield [26], Reich [27], and Bell [28] in areas where only pluviometric information is available. However, based on pluviographic information obtained in automatic stations, it was possible to adjust the expression of Bell [28] that represents the characteristics of the rainfall occurring in the study area. Based on this, the Chen method [23] is used to obtain the IDF curves and consequently their updating for the return periods of 10, 25, 50, and 100 years. Additionally, the spatial distribution of the rainfall-duration ratio (R) was obtained, which facilitates obtaining the IDF relationship with the Chen method [23] in places where there are no pluviographic and standard rain gauge stations.
Finally, the adequate estimation of IDF curves will allow and guarantee the planning and optimal design of hydraulic structures. Future flood hazard studies and consequently disaster risk management will also benefit directly from this work.

Materials and Methods
The urban zone of Ensenada and its surroundings has been selected as the study area, located on the Pacific coast of Baja California (31º30'-31º60'N; 116º50'-116º10' W. See Figure 1). This zone, in the northwest of Mexico, has a semiarid climate, with convective type rains characterized for being intense and of short term. The average annual rainfall is 273 mm, where the rainy season usually occurs between November and April. Topographic landscape is variable, with steep slopes, alluvial valleys, and an alluvial coastal plain [29]. This area is divided into four urban and seven rural subbasins that drain toward the Pacific coast. There are 11 weather stations inside the basins located at different elevations, that contain the main data to achieve the purpose of the study (Figure 1). The required data to develop the IDF relationship in the study zone were the maximum rainfall events, annually recorded at different durations; obtained from four automatic stations (pluviograph records). On the other hand, historical rainfall data of ten standard rain gauges (daily rainfall records) were collected and included in the calculations (Table 1). Based on the characteristics of the study area and available precipitation data, the methodology was defined with the steps shown in Figure 2.  Table 1 introduces four automatic stations. The Emilio Lopez Zamora station is the automatic station with more recorded data. It has 20 years of records, with measures from every ten minutes. This station is managed by a governmental agency called Comisión Nacional del Agua (CONAGUA). The other three automatic stations scarcely have 10 years of records, with measures for every five minutes. They are CICESE (CIC), Guadalupe (2036), and Ojos Negros (2035); all three were administered by Centro de Investigación Científica y de Educación Superior de Ensenada (CICESE). From these data, the maximum rainfall events of each year were identified and classified by their magnitude and duration in order to estimate the IDF curves.
Identification of each event was based on the continuous pluviograph records, from the rain start, until it stopped for more than one hour; e.g., if the rain ended at any time, but it suddenly started before an hour passed, this record was considered to be part of the previous event. Nonetheless, it was also considered a separate event during the day if, after one event had occurred, there were two hours without rain, and it started to rain again. Once all events were identified, they were classified by duration and magnitude to select the most intense of every year. Intensity is defined by the rain magnitude in mm, registered for a determined time; e.g., there are events with the same magnitude and different durations in hours, where the most intense is always the one with the shortest duration. Measures of the events recorded by the four automatic stations had different durations. In some automatic stations, the longest event recorded (not the most intense), was one of 180 minutes. Other intense events had a duration of only 10 minutes. However, most of the intense events had a duration of 120 minutes; therefore, this duration was selected to identify the intense events of each year. López-Lambraño developed this method to identify and classify the intense events [30], and it has been satisfactorily applied by Maldonado et al. [31].
Regarding rainfall information, there are 10 stations around the city of Ensenada (Table 1), which have records of maximum rainfall in 24 hours, from 1940 to 2011. In the same location of these stations there are three automatic stations (2072, 2036 and 2035), which have data from 1999 to 2019. Thus, the time series of the pluviometric data of these three stations were completed from the pluviograph records (automatic stations); e.g., for station 2072 a time series of 79 years was obtained. This was achieved through the daily accumulation of rainfall recorded every 10 minutes in that station until 2019. This same procedure was performed for stations 2036 and 2035.
A trend analysis was carried out on the pluviographic and pluviometric data, to determine the trend lines and verify if the registered rainfall tends to increase or decrease in the area. The stationarity in all the stations was also verified through the Augmented Dickey-Fuller test, which is a unit root test that allows accepting or rejecting the stationarity hypothesis in a time series [32]. Subsequently, the Chi square goodness and fit test was performed for the precipitation data analyzed with the probability distribution functions of Gumbel type I, Log-normal, Frechet and double Frechet to verify which of them presented the best fit.
The development of the IDF relationship is a procedure that starts with data availability, type of records, years of records, quality, and coverage, to choose an efficient method and claim the best results. Two methods, for different conditions, were applied in this case to estimate the IDF relationship for the return periods of 10, 25, 50, and 100 years. These periods, proposed by the National Agency for Prevention of Disasters, are suitable for the design of minor urban hydraulic structures and flood hazard assessments [21]. The first method to estimate the IDF relationship is applied to pluviograph data (from automatic stations). These data provide the real intensity that occurred for each rain event through the years and can be extrapolated applying a probability distribution. The Generalized Extreme Value Distribution (GEV) is usually applied to estimate IDF curves, like the extrapolation of the available data to estimate the peak value of the sample. The Gumbel's method, based on the theory of extreme values, is the first method applied here, where the probability of an event of a determined magnitude not being equaled or exceeded, can safely be adopted and have been widely used [30,33]. This is expressed by the Equation (1).
is the scale parameter,      is the location parameter or central value, e is the base of the Naperian logarithms, α, X, and μ; corresponding to the parameters of the statistic moments of the distribution. The distribution derivative provides the probability density function where the values of X, for different return periods (T), are estimated by means of: Confidence intervals are important to estimate the return period and data accuracy. Generally, the data to calculate the IDF curves have a standard error of 10% for short return periods and 20% for periods of 50 and 100 years [28]. By Gumbel's method, applied to the automatic data, the IDF curves can be estimated with a confidence value for a return period of 20 years, considering the 10 and 20 years of records from the three automatic stations. The IDF curves for a return period of 20 years are very useful for minor urban structures design. However, to develop the IDF curves for large return periods, such as 50 or 100 years, there would be uncertainty when applying the Gumbel method. Therefore, the requirement to apply Gumbel's method with a good confidence value, would be to have a long series of precipitation data.
Pluviograph records were useful to have real data of the events to estimate IDF for short return periods but were not enough for extrapolating rainfall intensities for the larger return periods. Therefore, the use of daily rainfall depth, from the standard rain gauges, was necessary to make the calculation of the rainfall intensity for all the return periods proposed. This allows us to counteract the calculations with the automatic data and disseminate better results.
Consequently, a second method to develop the IDF curves recommended for urban hydrological design in the Mexican Republic [24] is applied to daily rainfall records (from the 10 standard rain gauges). This method was developed by Chen as an alternative of the absent pluviograph records [23], by the following equation: Where, T r t P is the intensity of precipitation in mm/h,  (Tr) is the return period in years, (t) is the duration in minutes, (a), (b), and (c) are parameters of regional characterization of the rain defined by the (R) ratio. The depth-duration ratio (R) is the most important parameter to estimate the IDF relationship, from daily records, by using the Chen equation for a specific geographic location, because this ratio is related to rain characteristics of the zone [34]. As shown before, Equation (4) is based on (R) ratio to determine the rain of one hour and a return period of ten years 10 1 P , and therefore, the rain for a given return period. The (R) ratio for any average condition of rainfall over any geographical areas, also, has been proposed by Chen [23], through the following equation: Where, 2 1 ( ) P is the rain in mm, generated in one hour for a return period of 2 years. there are only daily records, not pluviograph records. Therefore, to validate Chen's method in this region, it was considered to make a good estimation of (R) ratio through the value of Several alternatives to apply Equation (5) for any geographic area, based only on daily records, have been given by Hershfield and Wilson through a diagram that relates the mean of maximum annual observations of precipitation days with the mean annual number of thunderstorm days; Reich, through a proposed world isopluvial map of 2-year/1-hour maximum precipitation; and Bell, based on Hershfield method through the following equation: The purpose was to find the R-value to 10 rain gauge stations applying Equation (5), where 2 24 ( ) P was easily estimated using daily records, by applying Gumbel distribution (Equations (2)  from the four automatic stations, through the Thiessen polygon. The averaged value was assigned to the seven rain gauges stations, to replace it in Equation (5) and have the first (R) value in the zone, as (R1).
The next step was to calculate 2 1 ( ) P for the 10 rain gauges stations through the methods previously described methods; Hershfield, Bell, and Reich [26][27][28]. Each value was replaced in Equation (5)  was derived from Bell´s formula (Equation (6)). After several tests of parameters variation, the new relationship resulted from modifying Equation (6) in the following equation: The 2 1 ( ) P values estimated from Equation (7) were replaced in Equation (5) to have new values of (R), such as (R4). Therefore, the evaluated methods provided four different results of (R1…4).
The challenge was to choose the best method of finding (R); thus, to address this challenge, an average of all calculated (R) was made, and this average was proposed as the best option to determine the (R) ratio for the city of Ensenada (Equation (8)). A summary of the process to find (R) ratio for each station in the zone is represented in Figure 3. Once the (R) ratio was determined for each rain gauge station, the spatial distribution of the (R) ratio was drawn on a map by interpolation, using the Inverse Distance Weighting (IDW) method. The error analysis of this spatial interpolation was carried out between actual and interpolated values using three statistical parameters: Mean Absolute Error (MAE), Root Mean Square Error (RMSE) and coefficient of determination (r 2 ). Consequently, the distribution of (R) ratio for the study area can be successfully used to calculate the IDF curves in the zone.
The procedure to develop the IDF relationship for the proposed return periods includes the combination of pluviograph and daily rainfall records, through the methods described above, testing different values of estimated (R), and comparing their results, to choose the most suitable, according to rain characteristics.
The first step was the calculations of the IDF relationship in three automatic stations (2072, 2036, and 2035), using the maximum rainfall events for each year, and applying the Gumbel method (Equation (1)).

P R P
The second step was the estimation of the IDF relationship for the 10 rain gauges stations, with daily records for more than 40 years, applying the Chen method (Equation (4)). There are three pluviographic stations (2072, 2036, and 2035) on the same location from the previous stations. The calculations using Chen equation included the use of the (R) ratios found with the different methods already mentioned, for all stations, followed by the comparison of the different results of IDF developments with the IDF values by the (SCT), and the comparison with the IDF relationship developed in 2011 by the by Norm of the State of Baja California, that only has results for station 2072.
Therefore, for each rain gauge station, with daily records for more than 40 years, there are several IDF calculations to enhance the comparison between each other, and have more options to choose the best results, developed through the following methods: A. IDF calculated from pluviograph records using Gumbel distribution (Equation (1) I. IDF calculated from daily records of ten stations (standard rain gauge), using Chen equation, where (R) is the ratio proposed calculated by Equation (8).
Once the IDF relationships were estimated from the previously established methods, two comparisons were made. The first comparison was made from the use of methods A, B, C, D, E, F, and G. The second comparison was made from the use of methods A, C, H, and I. From the first comparison method A and C were chosen again, considering that method A is very important because it was developed with pluviograph information (2072, 2036, and 2035 stations), representing the continuous measure of rainfall intensity. However, the rainfall time series data at stations do not exceed 20 years of record (Table 1). Thus, the method C was chosen because it corresponds to the official IDF relationship from the country. The chosen method and the corresponding results of the IDF curve for the city of Ensenada and its surroundings are shown in the results and discussion section.

Results and Discussion
Peak rainfall events of each year were identified from the records of the four automatic stations, following the previously described method. Table 2 provides the events extracted from the Emilio Lopez Zamora station (2072). Since this station has the most extensive time series of rainfall data (20 years), it will be taken to illustrate the maximum precipitation events.
It is clearly seen that the highest precipitation depth of each event occurs during the first 10 and 20 minutes. These observations confirm the short periods of rainfall events, as a local rain characteristic. Moreover, to analyze the increase or decrease of extreme events throughout the years, a rainfall trend has been drawn. The line trend of time series was analyzed from data of the station that have pluviograph and daily rainfall record. Figure 4 shows the line trend for data of station (2072). Figure 4a shows the trend of the high daily rainfall depth of every year, and Figure 4b presents the trend of the extreme rainfall events of the years at different durations, not necessarily the most intense (pluviograph records). The series presents a positive trend. In both cases, the trend indicates that maximum daily rainfall depth and peak rainfall events have been increasing with the years and will continue occurring in the future. In the trend analysis shown in Figure 4a, the slope establishes that the average increase in the maximum daily rainfall is 18%, i.e., if the current conditions that govern the occurrence of rainfall in the study area were to be maintained over time, then precipitation would increase at a rate of 1.8 mm per decade. For the case of Figure 4b there would be an increase of approximately 1.2 mm. The above magnitudes are considered significant if we consider that precipitation falls over a given coverage area, which translates into the generation of more direct runoff volume over the city of Ensenada.
These results can also be complemented by a detailed analysis of the rainfall time series for the State of Baja California [35].  The traditional methods used in hydrology to estimate rainfall and extreme flows for different return periods are based on the stationarity hypothesis in the probability distribution function of the series. According to Poveda et al. [36], this hypothesis could not be valid given the effects of climate change, climate variability, changes in land use, and the records of hydrological variables. Therefore, the Augmented Dickey-Fuller test was performed to evaluate the rainfall time series in the city of Ensenada fulfilling the hypothesis of stationarity or non-stationarity (Table 3).  Table 3 shows the results of the Augmented Dickey-Fuller test. The analysis was carried out for a 5% significance level, finding that none of the time series corresponding to weather stations comply with the stationarity assumption because the p-value is greater than 0.05 [32]. Given the previous analysis, the Chi-square goodness and fit test was used for the precipitation data analyzed with the Gumbel type I, Log-normal, Frechet and double Frechet probability distribution functions to verify which of them had the best fit. In the case of the rain gauge stations, eight were better adjusted to the Gumbel type I distribution function, and two stations to the double Frechet. The previous analysis was carried out with a 0.05 confidence level. In the case of the automatic station 2072, both the Gumbel type I distribution function and the double Frechet distribution function presented equivalent results. Given the above, it was decided to use the Gumbel distribution function for the analysis of the Year relationship Intensity-Duration-Frequency of rainfall in the study area. In addition, it has been found that the chosen distribution fits well for precipitation in semi-arid areas [11].
The (R) ratio was estimated based on the 2 1 ( ) P following the process shown in Figure 3. The result of the regional review of  Table 5. The (R) ratio calculated for each station varied from 0.35 to 0.44; this range is close to the (R) ratio of 0.44 for one station in Tijuana reported by Campos-Aranda [24], which indicates that both cities share the same rain characteristics. Thus, once the (R) ratio was estimated, the spatial distribution of this ratio was projected in a map for the study area ( Figure 5). The error analysis was carried out between actual and interpolated values using three statistical parameters: Mean Absolute Error (MAE = 2.341E-14), Root Mean Square Error (RMSE = 1.53011E-07), and coefficient of determination (R 2 =1.0000E+00). There is a significant positive correlation between actual and interpolated values estimated for the depth-duration ratio. Values for MAE and RMSE indicate that the IDW method created a good interpolated surface for the whole area of interest based on the observed data. It can be established that in Ensenada City the average value of R is about 0.39 and this value is proposed as (R) ratio for the area.
It is highly important to highlight that in the absence of pluviograph data, the (R) ratio calculated becomes the key to develop a satisfactory IDF curves for the standard rain gauge stations, by applying Chen equation.  ( ) P calculated by the adjusting and modification of Bell [28] (Equation (7)).

Station
IDF curves are commonly developed using historical annual maximum precipitation, this involves utilization of long-term historical rainfall observations. When sub-daily rainfall records are not available, the characteristics of extreme rainfall intensities, and subsequently, their distribution functions corresponding to the short durations might not be captured [39][40][41][42][43]. This problem is clearly addressed by applying Equation (7) to rain duration ratio (R4) estimated (Table 5), likewise distribution of the depth-duration ratio ( Figure 5). In this way, empirical formulas (e.g. Chen) can be used in areas where there are no pluviographs (high temporal resolution), rain gauges, and/or information available.
The IDF relationship were calculated with pluviograph records, by applying Gumbel distribution and, calculated with daily records by applying Chen methods through different options of (R) ratio. The different results allowed us to make comparisons to choose the most suitable method for calculations of the IDF relationship.
The comparisons of results were focused on the three stations that have pluviograph and daily records (2072, 2036, and 2035). However, station 2072 has been chosen to define the methods of comparisons, because it has the widest time series of pluviograph and daily records. Table 6 presents the results of the first comparison that tested methods A, B, C, D, E, F, and G previously described, and its average likelihood method (H), to assess the IDF relationship. In this table the values of the IDF relationship obtained by the different methods vary considerably, and for this reason it was decided to carry out an average of these methods. However, it is important to highlight that method G presents the highest magnitudes and correspond to the IDF reported in the Baja California State Standards. It should be noted that these magnitudes are far from the average. With these results, it was identified that these IDFs are over-estimated, which implies a greater over determination of dimensions in the designs of hydraulic structures, and so therefore higher construction costs.
From this comparison, another comparison was selected to assess and estimate the IDF relationship for the rest of the rain gauge stations. The comparison defines the selection of the best method of IDF estimation, that includes the IDF estimated using pluviograph data (method A), the IDF reported by the SCT (method C), the IDF estimated with the R ratio proposed in this zone (method I), and the IDF from method H. This comparison is shown in Figures 6 to 8 for the stations 2072, 2036, and 2035.  (7)). (E) Estimated from the Chen equation using R3 from table 5 (Reich, 1963). (F) Estimated from the Chen equation using R2 from table 5 (Equation (6)). (G) Obtained from the of Norm by the State of Baja California.  Figure 6a shows the IDF curves from method (A), estimated from pluviograph records using Gumbel (Equation (1)). Figure 6b shows the IDF curves from method (C) estimated from isohyetal map of the SCT. Figure 6c shows the IDF curves from method (H) estimated from the average of IDF developed by methods A, B, C, D, E, F, and G, and Figure 6d shows the IDF curves from method (I) estimated from Chen equation using the (R) ratio proposed (Equation (8)).
The comparisons of IDF relationship for the 2072 station showed similar values calculated with methods (I) and (H). The method (I) is calculated by using Chen with the proposed (R) ratio, and the method (H) is calculated from the average of IDF developed by methods A, B, C, D, E, F, and G described before. On the other hand, the intensities obtained with method C are lower than those obtained with method I. This may be because the IDFs obtained with method C have not been updated since 2000, therefore, they do not involve the analyzes performed in figure 4a and 4b. Figures 6-8 show that the results of the method (I) are closer to methods (C) and (H) than method (A). Therefore, considering the position of the results of (I) on the range of values of all methods, it could be suggested as the most suitable to calculate the IDF curves. The values to estimate Figures 6-8 can be seen in the supplementary (table S1 to table S3) material section. Method (I) was chosen to calculate IDF for the 10 rain gauge stations like the best option of the methods tested. Table 7 shows the IDF relationship for the 10 stations calculated by the Chen method through the proposed (R) ratio.  Figure 7a shows the IDF curves from method (A), estimated from pluviograph records using Gumbel (Equation (1)). Figure 7b shows the IDF curves from method (C) estimated from isohyetal map of the SCT. Figure 7c shows the IDF curves from method (H) estimated from the average of IDF developed by methods A, B, C, D, E, and F, and Figure 7d shows the IDF curves from method (I) estimated from Chen equation using the (R) ratio proposed (Equation (8)).  Figure 8a shows the IDF curves from method (A), estimated from pluviograph records using Gumbel (Equation (1)). Figure 8b shows the IDF curves from method (C) estimated from isohyetal map of the SCT. Figure 8c shows the IDF curves from method (H) estimated from the average of IDF developed by methods A, B, C, D, E, and F, and Figure 8.d shows the IDF curves from method (I) estimated from Chen equation using the (R) ratio proposed (Equation (8)).
For the stations 2036 and 2035, the results of the IDF relationship were similar by methods (C) and (I). However, when graphing the IDF curves with the C method, it is noted that they do not show the typical fit of said curves. This may be because there are errors in the process of estimating them. A comparison has also been made for the stations 2001, 2005, and 2045, that provided equivalent results by applying methods (C) and (I). From the analysis of the comparison of results, method I is recommended as the most suitable for the IDF design in the evaluated stations. For the IDF relationship calculation from daily records to be used in Mexico, the Chen's method has been mostly recommended by other authors [24,44]. However, considering that rain changes in space and time, the use of this method has been validated through the accurate calculation of (R) ratio, where it depended on the 2 1 ( ) P . The results of 2 1 ( ) P obtained by pluviograph data, Reich and the adjusting from Bell equation, are similar; however, the method widely used is Hershfield that provides the highest values. Therefore, it was considered to calculate (R) with each value of 2 1 ( ) P of each method to make an average of all estimated (R). This average was proposed as the (R) ratio of the zone that can be safe to calculate the IDF curves in other area stations.
The differences between the various methods in this study can be attributed to the length of the time series used. Finally, it has been verified that the quality and length of the pluviographic and rainfall precipitation records are important aspects in the study of the intensity-duration-frequency relationship in an area.

Conclusions
The procedure and methods for developing IDF relationship in Ensenada were carefully tested making effective use of the available rainfall data, including pluviograph and daily records that allowed to achieve the research objectives.
Despite the absence of enough pluviograph data in this study, the IDF relationship was carried out successfully using the combination of Chen's methods, through the average (R) calculated from the average value of the one-hour rainfall and the two-year return period. The values of 2 1 ( ) P obtained were compared to the values of some cities in California and Baja California, with a range between 10 and 16.61 mm, and the values of the (R) ratio in a range between 0.35 and 0.44; this range is close to the (R) ratio of 0.44 for one station in Tijuana, a city 100 km farther North from Ensenada. The values found here correspond to the rainfall characteristics of the zone; therefore, the method used in this study can be replicated to another semi-arid zones with the same rain characteristics. The parameterization of the IDF relation for different durations, allows better understanding and realization of spatio-temporal analysis of the characteristics of rainfall in an area.
Chen's method is applicable to any zone if the (R) ratio is well defined. This value was carefully reviewed in the region and carefully tested with updated pluviograph and daily data. Therefore, the (R) ratio estimated is proposed to develop IDF curves in the absence of pluviograph data.
After analyzing the results of the IDF relationship for the station 2072, it was observed that the IDF relationship published by the Norm of the State of Baja California presents the highest values of all methods. Although they are also safe for designing, they imply the highest cost of constructiongreater sizing-that could be minimized considering the new results. The document review indicates that the IDF relationships were developed with data available up to the year 2011. Therefore, it is suggested that these results from the IDF relationship should be included in the Norm of the State of Baja California, as it requires the recurrence update upon its recommendation.
This study guarantees the following aspects: input to rain-runoff models to improve the information available for an adequate design, planning, estimation of dimensions of civil works, and the integral management of water resources in this semi-arid region. Also, the proposed intensityduration-frequency relationship will facilitate the evaluation of the flood hazard. The city of Ensenada periodically suffers the effects of floods, is out of electricity service, it is not possible to travel, there are drainage problems, health hazards, and work activities are often stopped. Therefore, proper risk management is important to protect not only the lives of citizens, but also the public and private assets, as well as the progress made in the development process. These results constitute a starting point for risk management and flood resilience.
Supplementary Materials: The following are available online at www.mdpi.com/2306-5338/7/4/78/s1, Table S1: Established comparison of the IDF relationship estimated through different methods for station 2072 for different return periods, Table S2: Established comparison of IDF relationship estimated through different methods for the station 2036 for different return periods, Table S3