Next Article in Journal
Evaluation of Water Quality Parameters in Shatt AL-Arab, Southern Iraq, Using Spatial Analysis
Next Article in Special Issue
Intra-Storm Pattern Recognition through Fuzzy Clustering
Previous Article in Journal
Analytical Protocol to Estimate the Relative Importance of Environmental and Anthropogenic Factors in Influencing Runoff Quality in the Bumbu Watershed, Papua New Guinea
Previous Article in Special Issue
Sensitivity Analysis of the Rainfall–Runoff Modeling Parameters in Data-Scarce Urban Catchment
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

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

by
Ena Gámez-Balmaceda
1,2,
Alvaro López-Ramos
3,
Luisa Martínez-Acosta
2,3,
Juan Pablo Medrano-Barboza
3,
John Freddy Remolina López
4,
Georges Seingier
1,
Luis Walter Daesslé
1 and
Alvaro Alberto López-Lambraño
1,2,5,6,*
1
Instituto de Investigaciones Oceanólogicas, Universidad Autónoma de Baja California, Baja California 22860, Mexico
2
Faculty of Engineering, Architecture and Design, Universidad Autónoma de Baja California, Baja California 22860, Mexico
3
GICA Group, Faculty of Civil Engineering, Universidad Pontificia Bolivariana campus Montería, Montería 230002, Córdoba, Colombia
4
ITEM Group, Faculty of Electronic Engineering, Universidad Pontificia Bolivariana campus Montería, Montería 230002, Córdoba, Colombia
5
Hidrus S.A. de C.V., Ensenada 22760, Mexico
6
Grupo Hidrus S.A.S., Montería 230002, Colombia
*
Author to whom correspondence should be addressed.
Hydrology 2020, 7(4), 78; https://doi.org/10.3390/hydrology7040078
Submission received: 9 September 2020 / Revised: 8 October 2020 / Accepted: 12 October 2020 / Published: 15 October 2020
(This article belongs to the Special Issue Advances in Modelling of Rainfall Fields)

Abstract

:
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 ( P 1 2 ) 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.

1. 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 semi-arid 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 (2000–2020), 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 ( P 1 2 ) of two years. ( P 1 2 ) 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.

2. 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).
F ( X ) = e e X μ α   ( < X )
where 0 < α < 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:
X T = μ + α Y T
Y T = ln ln T T 1
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:
P t T r = a P 1 10 log 10 2 X T r X 1 t 60 ( t + b ) c
where, P t T r is the intensity of precipitation in mm/h, P 1 10 = R ( P 24 10 ) , P 1 10 is the rain in mm, generated in one hour for a return period of 10 years, (X) is the ratio of the rain-return period X = P t 100 P t 10 , P t 100 and P t 10 is the rain of 24 hours and return period of 100 and 10 years respectively. (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 P 1 10 , 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:
R = P 1 2 P 24 2
where, ( P 1 2 ) is the rain in mm, generated in one hour for a return period of 2 years. ( P 24 2 ) is the rain in mm, generated in 24 h for a return period of 2 years. This value is easily estimated by applying Gumbel distribution for each standard rain gauge. ( P 1 2 ) can be estimated from pluviograph records applying the Gumbel distribution for each station. However, the issue to apply formula 5 is finding the value of the rain of one hour and a return period of two years ( P 1 2 ) , when 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 ( P 1 2 ) .
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-h maximum precipitation; and Bell, based on Hershfield method through the following equation:
P 1 2 = 0.17 M N 0.33
Equation (6) could be applied if 0 < M ≤ 2.0, 1 < N ≤ 80, in which ( P 1 2 ) = 2-yr, 1-h rainfall in inches: M = mean of maximum annual observational—precipitation day in inches: and N = mean annual number of thunderstorm days.
The purpose was to find the R-value to 10 rain gauge stations applying Equation (5), where ( P 24 2 ) was easily estimated using daily records, by applying Gumbel distribution (Equations (2) and (3)). In the case of ( P 1 2 ) there are three stations with pluviograph records and seven stations with daily records, the issue was to find ( P 1 2 ) for the seven stations with daily records. Consequently, it was decided to evaluate the methods available for estimating the value of ( P 1 2 ) and choosing the best way of finding the accurate ( P 1 2 ) to calculate the (R) ratio.
The process started by knowing the approximated value of the ( P 1 2 ) along the regional area to compare them with the estimated values through the mentioned methods, started the process. Approximated values of ( P 1 2 ) in the regional area were examined through a literature review of regional studies. This review, with the same rainfall characteristics, includes California in the US (CA), near to the international border, and Baja California in Mexico, as shown in Figure 1.
The first method evaluated to calculate ( P 1 2 ) was the Gumbel distribution using pluviograph records from stations: CIC, 2072, 2036, and 2035. Then, a weighted average of ( P 1 2 ) was estimated 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 ( P 1 2 ) 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) to have different values of (R) like (R2) and (R3). Three different (R) values were estimated. Nevertheless, from the different values of ( P 1 2 ) , it was considered to search for a relationship to provide values of ( P 1 2 ) for the rain gauges stations close to the values of ( P 1 2 ) estimated from pluviograph records. The Hershfield and Bell methods provide the same results. This way, a new method to find ( P 1 2 ) 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:
P 1 2 = 0.12 M N 0.33
The ( P 1 2 ) 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.
R ¯ = i = 1 n R n
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 (r2). 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)).
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:
  • IDF calculated from pluviograph records using Gumbel distribution (Equation (1)); this included only three stations (2072, 2036, and 2035).
  • IDF calculated from daily records of ten stations (standard rain gauge), using the Chen equation, where (R) was calculated through the weighted average of the ( P 1 2 ) from automatic stations, applying Gumbel distribution, implementing (R1).
  • IDF obtained by the isohyetal map of the (SCT) only for stations 2072, 2036, 2035, 2001, 2005, and 2045.
  • IDF calculated from daily records of ten stations (standard rain gauge), using the Chen equation, where (R) was calculated through the ( P 1 2 ) by the modification of Bell equation [28], implementing (R4).
  • IDF calculated from daily records of ten stations (standard rain gauge), using the Chen equation, where (R) was calculated through the ( P 1 2 ) by the world isopluvial map of Reich [27], implementing (R3).
  • IDF calculated from daily records of ten stations (standard rain gauge), using the Chen equation, where (R) was calculated through the ( P 1 2 ) by Hershfield and Bell [26,28], implementing (R2), both methods provide the same results of ( P 1 2 ) .
  • IDF obtained by the Norm of the State of Baja California only for station 2072.
  • IDF calculated by the average of IDF developed through options A, B, C, D, E, F, and G.
  • 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.

3. 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 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 ( P 1 2 ) following the process shown in Figure 3. The result of the regional review of ( P 1 2 ) is shown in Table 4. The regional review reveals that values of ( P 1 2 ) from Southern California to Northern Baja California are in the range between 10 and 16.61 mm, indicating that values of ( P 1 2 ) for the study area should be estimated in this range. For this reason, the evaluation of the different methods to calculate ( P 1 2 ) and estimate (R) followed by the comparison of their results was carefully analyzed to determine the accurate (R) ratio. The results of ( P 1 2 ) are in the range of the regional review, when ( P 1 2 ) 1 is calculated by the average of pluviograph and daily records applying Gumbel, ( P 1 2 ) 2 is calculated by Hershfield [26], ( P 1 2 ) 3 is calculated by Reich [27], and ( P 1 2 ) 4 is calculated by the adjusting and modification of Bell [28]. The results of (R1...4) ratios estimated for each station through the different methods of ( P 1 2 ) were averaged to define and support the proposed (R) for each station. These values are shown in 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 (R2=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.
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 Figure 6, Figure 7 and Figure 8 for the stations 2072, 2036, and 2035.
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,b.
Figure 6, Figure 7 and Figure 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 Figure 6, Figure 7 and Figure 8 can be seen in Tables S1–S3 in Supplementary Materials. 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.
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 ( P 1 2 ) . The results of ( P 1 2 ) 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 ( P 1 2 ) 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.

4. 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 ( P 1 2 ) 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 construction—greater 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 intensity–duration–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 https://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: Established comparison of the IDF relationship estimated through different methods for the station 2035 for different return periods.

Author Contributions

Conceptualization, E.G.-B.; A.A.L.-L.; methodology, E.G.-B., A.A.L.-L., J.P.M.-B., A.L.-R. and J.F.R.L; formal analysis, E.G.-B., G.S., L.W.D., A.A.L.-L., L.M.-A., J.P.M.-B., and A.L.-R.; writing—original draft preparation, E.G.-B., A.A.L.-L., L.M.-A., J.P.M.-B., A.L.-R. and J.F.R.L.; writing—review and editing, E.G.-B., A.A.L.-L., L.M.-A., A.L.-R. and J.F.R.L.; visualization, E.G.-B. and L.M.-A.; supervision, A.A.L.-L.; project administration, E.G.-B., A.A.L.-L., L.M.-A., A.L.-R. and J.F.R.L.; funding acquisition, J.F.R.L. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by Universidad Autónoma de Baja California and HIDRUS S.A de C.V, Grupo HIDRUS S.A.S and Universidad Pontificia Bolivariana Campus Montería, the APC was funded by Universidad Pontificia Bolivariana Campus Montería.

Acknowledgments

We would like to thank the Centro de Investigación Científica y de Educación Superior de Ensenada (CICESE), and especially Santiago Higareda from the meteorology laboratory, who provided the precipitation data of automatic station.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. De Paola, F.; Giugni, M.; Topa, M.E.; Bucchignani, E. Intensity-Duration-Frequency (IDF) rainfall curves, for data series and climate projection in African cities. Springer 2014, 3, 1–18. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  2. Raghunath, H.M. Hydrology Principles-Analysis-Design; New Age International: New Delhi, India, 2006; ISBN 9788122423327. [Google Scholar]
  3. Sherif, M.; Chowdhury, R.; Shetty, A. Rainfall and Intensity-Duration-Frequency (IDF) Curves in United Arab Emirates. In Proceedings of the World Environmental and Water Resources Congress, Portland, OR, USA, 1–5 June 2014. [Google Scholar]
  4. Gericke, O.J.; Du Plessis, J.A. Evaluation of the standard design flood method in selected basins in South Africa. J. S. Afr. Inst. Civ. Eng. 2012, 54, 2–14. [Google Scholar]
  5. Lin, X. Flash Floods in Arid and Semi-Arid Zones; UNESCO: Paris, France, 1999. [Google Scholar]
  6. Cruz Aguirre, R.U. Además de “Rosa”, Ocho Huracanes más han Pasado a Menos de 200 km de Ensenada. Available online: https://todos.cicese.mx/sitio/noticia.php?n=1215#.XdhI15NKjIV (accessed on 31 August 2020).
  7. Pacheco, B. Tiene Baja California Historial de Desastres—El Vigía. Available online: https://www.elvigia.net/general/2015/12/10/tiene-baja-california-historial-desastres-220040.html (accessed on 31 August 2020).
  8. Rodríguez Esteves, J.M. La conformación de los “desastres”—Construcción social del riesgo y variabilidad climática en Tijuana, B.C. Front. Norte 2007, 19, 83–112. [Google Scholar]
  9. Arnell, V. Analysis of Rainfall Data for Use in Design of Storm Sewer Systems; International Conference on Urban Storm Drainage: Southampton, UK, 1978; pp. 1–15. [Google Scholar]
  10. Campos-Aranda, D.F. Relación y estimación de predicciones de lluvia horaria-diaria en dos zonas geográficas de México. Tecnol. Cienc. Agua 2012, 3, 141–152. [Google Scholar]
  11. López-Lambraño, A.; Fuentes, C.; González-Sosa, E.; López-Ramos, A. Pérdidas por intercepción de la vegetación y su efecto en la relación intensidad, duración y frecuencia (IDF) de la lluvia en una cuenca semiárida. Tecnol. Cienc. Agua 2017, 8, 37–56. [Google Scholar]
  12. DeGaetano, A.T.; Castellano, C.M. Future projections of extreme precipitation intensity-duration-frequency curves for climate adaptation planning in New York State. Clim. Serv. 2017, 5, 23–35. [Google Scholar] [CrossRef]
  13. Groisman, P.Y.; Knight, R.W.; Karl, T.R. Changes in intense precipitation over the Central United States. J. Hydrometeorol. 2012, 13, 47–66. [Google Scholar] [CrossRef]
  14. WHO; UNEP. Managing the Risks of Extreme Events and Disasters to Advance Climate Change Adaptation; Cambridge University Press: New York, NY, USA, 2012; p. 594. [Google Scholar]
  15. Jaklič, A.; Šajn, L.; Derganc, G.; Peer, P. Automatic digitization of pluviograph strip charts. Meteorol. Appl. 2016, 23, 57–64. [Google Scholar] [CrossRef] [Green Version]
  16. Llasat, M.-C. An Objective Classification of Rainfall Events on the Basis of their Convective Features. Application to Rainfall Intensity in the North-East of Spain. Int. J. Climatol. 2001, 21, 1385–1400. [Google Scholar] [CrossRef]
  17. Organización Meteorológica Mundial. Guía del Sistema Mundial de Observación; Organización meteorológica mundial: Geneva, Switzerland, 2010; ISBN 9789263304889. [Google Scholar]
  18. Plummer, N.; Allsopp, T.; Lopez, J.A. Guidelines on Climate Observation Networks and Systems; Llansó, P., Ed.; World Meteorological Organization: Geneva, Switzerland, 2003. [Google Scholar]
  19. Secretaría de Infraestructura y Desarrollo Urbano. Normas Técnicas para Elaboración de Proyectos de Alcantarillado Pluvial en el Estado de Baja California. 2013. Available online: http://www.ordenjuridico.gob.mx/Documentos/Estatal/Baja%20California/wo83549.pdf (accessed on 31 August 2020).
  20. Secretaría Comunicaciones y Transporte. Isoyetas de Intensidad-Duración-Periodo de Retorno para la República Mexicana—Baja California; Secretaria Comunicaciones y Transporte: Mexicali, Mexico, 2014; Available online: http://www.sct.gob.mx/carreteras/direccion-general-de-servicios-tecnicos/isoyetas/ (accessed on 31 August 2020).
  21. SEDATU. Términos de Referencia para Elaborar Atlas de Riesgos; Secretaría de Desarrollo Agrario, Territorial y Urbano: Ciudad de México, Mexico, 2016; p. 121. [Google Scholar]
  22. Secretaría de Gobernación y Cenapred. Guía Básica para la Elaboración de Atlas Estatales y Municpales de Peligros y Riesgos; Violeta, R., Ed.; Secretaría de Gobernación y Cenapred: Ciudad de México, Mexico, 2014; ISBN 9706289054. [Google Scholar]
  23. Chen, C. Rainfall intensity-duration-frequency formulas. J. Hydraul. Eng. 1983, 109, 1603–1621. [Google Scholar] [CrossRef]
  24. Campos-Aranda, D. Intensidades máximas de lluvia para diseño hidrológico urbano en la república mexicana Rainfall Maximum Intensities for Urban Hydrological Design in Mexican Republic. Ing. Investig. Tecnol. 2010, 11, 179–188. [Google Scholar]
  25. Lorente Castelló, J.; Casas Castillo, M.C.; Rodríguez Sola, R.; Redaño Xipell, Á. Intensidades extremas y precipitación máxima probable. Fenóm. Meteorol. Advers. Esp. 2013, Volumen 5, 142–155. [Google Scholar]
  26. Hodges, L.H.; Hershfield, D.M.; Washington, D.C.; Reichelderfer, F.W. Rainfall Frequency Atlas of the United States for Durations from 30 Minutes to 24 Hours and Return Periods from 1 to 100 Years; Weather Bureau: Washington, DC, USA, 1961; p. 65. [Google Scholar]
  27. Reich, B.M. Short-Duration Rainfall-Intensity Estimates and Other Design Aids for Regions of Sparse Data. J. Hydrol. 1963, 1, 3–28. [Google Scholar] [CrossRef]
  28. Bell, F.C. Generalized Rainfall-Duration-Frequency Relationships. J. Hydraul. Div. 1969, 95, 311–327. [Google Scholar]
  29. Smith, S.V.; Bullock, S.H.; Hinojosa-Corona, A.; Franco-Vizcaíno, E.; Escoto-Rodríguez, M.; Kretzschmar, T.G.; Farfán, L.M.; Salazar-Ceseña, J.M. Soil erosion and significance for carbon fluxes in a mountainous Mediterranean-climate watershed. Ecol. Appl. 2007, 17, 1379–1387. [Google Scholar] [CrossRef]
  30. López-Lambraño, A. Elaboración de las Curvas de Intensidad-Duración-Frecuencia de la Lluvia en las Estaciones Pluviográficas Aeropuerto los Garzones, Universidad de Córdoba y Turipana, Localizadas en la Cuenca Media del Río Sinú, Montería; Universidad Pontificia Bolivariana: Montería, CO, USA, 2001. [Google Scholar]
  31. Maldonado, M.; Martínez, G.; Matajira, F.F. Elaboración curvas IDF estaciones: Cinera Villa Olga y Santa Isabel—Municipio de Cúcuta—Colombia. Rev. Ambient. Agua Aire Suelo 2007, 2, 80–94. [Google Scholar]
  32. Said, S.E.; Dickey, D.A. Testing for Unit Roots in Autoregressive-Moving Average Models of Unknown Order. Biometrika 1984, 71, 599–607. [Google Scholar] [CrossRef]
  33. Courty, L.G.; Wilby, R.L.; Hillier, J.K.; Slater, L.J. Intensity-duration-frequency curves at the global scale. Environ. Res. Lett. 2019, 14, 1–11. [Google Scholar] [CrossRef]
  34. Hersfield, D.M.; Wilson, W.T. Generalizing of Rainfall-Intensity-Frequency Data. AIHS. Gen. Ass. Toronto 1957, 1, 499–506. [Google Scholar]
  35. López-Lambraño, A.A.; Fuentes, C.; López-Ramos, A.A.; Mata-Ramírez, J.; López-Lambraño, M. Spatial and temporal Hurst exponent variability of rainfall series based on the climatological distribution in a semiarid region in Mexico. Atmosfera 2018, 31, 199–219. [Google Scholar] [CrossRef]
  36. Poveda, G.; Álvarez, D.M. El colapso de la hipótesis de estacionariedad por cambio y variabilidad climática: Implicaciones para el diseño hidrológico en ingeniería. Rev. Ing. 2012, 36, 65–76. [Google Scholar]
  37. Frevert, R.; Schwab, G.; Edminster, T.; Barnes, K. Soil and Water Conservation Engineering; John Wiley & Sons Ltd.: New York, NY, USA, 1963. [Google Scholar]
  38. Robie, R.B. Rainfall Analysis for Drainage Design; Department of Water Resources: Sacramento, CA, USA, 1976. [Google Scholar]
  39. NOAA’s, N.W.S. NOAA Atlas 14 Point Precipitation Frequency Estimates: CA. Available online: https://hdsc.nws.noaa.gov/hdsc/pfds/pfds_map_cont.html?bkmrk=ca (accessed on 31 August 2020).
  40. Nhat, M.L.; Tachikawa, Y.; Takara, K. Establishment of Intensity-Duration-Frequency Curves for Precipitation in the Monsoon Area of Vietnam. Annu. Disas. Prev. Res. Inst. 2006, 49B, 93–103. [Google Scholar]
  41. Tfwala, C.M.; van Rensburg, L.D.; Schall, R.; Mosia, S.M.; Dlamini, P. Precipitation intensity-duration-frequency curves and their uncertainties for Ghaap plateau. Clim. Risk Manag. 2017, 1–9. [Google Scholar] [CrossRef]
  42. Koutsoyiannis, D.; Kozonis, D.; Manetas, A. A mathematical framework for studying rainfall intensity-duration-frequency relationships. J. Hydrol. 1998, 206, 118–135. [Google Scholar] [CrossRef]
  43. Sun, Y.; Wendi, D.; Kim, D.E.; Liong, S.Y. Deriving intensity–duration–frequency (IDF) curves using downscaled in situ rainfall assimilated with remote sensing data. Geosci. Lett. 2019, 6, 1–12. [Google Scholar] [CrossRef] [Green Version]
  44. Domínguez, R.; Carrizosa, E.; Fuentes, G.E.; Arganis, M.L.; Osnaya, J.; Galván-Torres, A.E. Análisis regional para estimar precipitaciones de diseño en la república mexicana. Tecnol. Cienc. Agua 2018, 9, 5–29. [Google Scholar] [CrossRef]
Figure 1. Study area location showing sub-basins and distribution of weather stations.
Figure 1. Study area location showing sub-basins and distribution of weather stations.
Hydrology 07 00078 g001
Figure 2. General diagram of the methodology.
Figure 2. General diagram of the methodology.
Hydrology 07 00078 g002
Figure 3. Procedure to determine the (R) ratio proposed. R ¯ is the average ratio and n is the number of (R1…4) calculated for each ( P 1 2 ) .
Figure 3. Procedure to determine the (R) ratio proposed. R ¯ is the average ratio and n is the number of (R1…4) calculated for each ( P 1 2 ) .
Hydrology 07 00078 g003
Figure 4. Trend analysis for station 2072 using automatic and daily records. (a) Annual maximum daily rainfall depth. (b) Annual Peak rainfall events.
Figure 4. Trend analysis for station 2072 using automatic and daily records. (a) Annual maximum daily rainfall depth. (b) Annual Peak rainfall events.
Hydrology 07 00078 g004
Figure 5. Distribution of the Depth-Duration Ratio (R) for the city of Ensenada.
Figure 5. Distribution of the Depth-Duration Ratio (R) for the city of Ensenada.
Hydrology 07 00078 g005
Figure 6. Established comparison of the IDF relationship estimated through different methods for station 2072 for different return periods. (a) shows the IDF curves from method (A), estimated from pluviograph records using Gumbel (Equation (1)). (b) shows the IDF curves from method (C) estimated from isohyetal map of the SCT. (c) 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 (d) shows the IDF curves from method (I) estimated from Chen equation using the (R) ratio proposed (Equation (8)).
Figure 6. Established comparison of the IDF relationship estimated through different methods for station 2072 for different return periods. (a) shows the IDF curves from method (A), estimated from pluviograph records using Gumbel (Equation (1)). (b) shows the IDF curves from method (C) estimated from isohyetal map of the SCT. (c) 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 (d) shows the IDF curves from method (I) estimated from Chen equation using the (R) ratio proposed (Equation (8)).
Hydrology 07 00078 g006
Figure 7. Established comparison of the IDF relationship estimated through different methods for station 2036 for different return periods. (a) shows the IDF curves from method (A), estimated from pluviograph records using Gumbel (Equation (1)). (b) shows the IDF curves from method (C) estimated from isohyetal map of the SCT. (c) shows the IDF curves from method (H) estimated from the average of IDF developed by methods A, B, C, D, E, and F, and (d) shows the IDF curves from method (I) estimated from Chen equation using the (R) ratio proposed (Equation (8)).
Figure 7. Established comparison of the IDF relationship estimated through different methods for station 2036 for different return periods. (a) shows the IDF curves from method (A), estimated from pluviograph records using Gumbel (Equation (1)). (b) shows the IDF curves from method (C) estimated from isohyetal map of the SCT. (c) shows the IDF curves from method (H) estimated from the average of IDF developed by methods A, B, C, D, E, and F, and (d) shows the IDF curves from method (I) estimated from Chen equation using the (R) ratio proposed (Equation (8)).
Hydrology 07 00078 g007
Figure 8. Established comparison of the IDF relationship estimated through different methods for station 2035 for different return periods. (a) shows the IDF curves from method (A), estimated from pluviograph records using Gumbel (Equation (1)). (b) shows the IDF curves from method (C) estimated from isohyetal map of the SCT. (c) shows the IDF curves from method (H) estimated from the average of IDF developed by methods A, B, C, D, E, and F, and (d) shows the IDF curves from method (I) estimated from Chen equation using the (R) ratio proposed (Equation (8)).
Figure 8. Established comparison of the IDF relationship estimated through different methods for station 2035 for different return periods. (a) shows the IDF curves from method (A), estimated from pluviograph records using Gumbel (Equation (1)). (b) shows the IDF curves from method (C) estimated from isohyetal map of the SCT. (c) shows the IDF curves from method (H) estimated from the average of IDF developed by methods A, B, C, D, E, and F, and (d) shows the IDF curves from method (I) estimated from Chen equation using the (R) ratio proposed (Equation (8)).
Hydrology 07 00078 g008
Table 1. Rain gauge network from the study zone.
Table 1. Rain gauge network from the study zone.
Station NameStation IDElevation above Sea Level (M)Type of RecordsLatLongRecord LengthSource of Data
Emilio Lopez Zamora20724310 min
24 h
31.89−116.601999–2019
1940–2019
CONAGUA
CONAGUA/CICESE
CICESECIC605 min31.86−16.662007–2019CICESE
Guadalupe20363615 min
24 h
32.02−116.612009–2019
1954–2019
CICESE
CONAGUA
Ojos Negros20356805 min
24 h
31.91−116.232009–2019
1948–2019
CICESE
CONAGUA
Agua Caliente200140024 h32.10−116.451969–2011CONAGUA
Boquilla de Santa Rosa200525024 h32.02−116.771969–2011CONAGUA
San Carlos204516424 h31.78−116.461969–2011CONAGUA
Santo Tomas206518024 h31.79−116.401969–2011CONAGUA
El Alamar207971024 h31.83−116.201969–2011CONAGUA
Maneadero21065024 h31.69−116.571969–2011CONAGUA
Punta Banda21081524 h31.71−116.661969–2010CONAGUA
5-min: automatic stations, 24 h (daily): rain gauge stations. Lat: latitude, Long: longitude.
Table 2. Peak rainfall events for ELZ station (2072).
Table 2. Peak rainfall events for ELZ station (2072).
Peak Rainfall Events (Mm) at Different Durations (ELZ -2072 )
Years10 min20 min30 min60 min120 min180 min
19993.305.085.585.587.359.12
20006.358.138.1210.4014.2320.60
20015.596.858.129.1410.9013.40
20023.054.324.577.1110.1611.70
20032.543.555.096.6011.1812.40
20044.837.119.1410.9018.7924.60
20053.565.346.6111.7013.7213.70
20064.327.379.1412.4018.2918.30
20074.328.1312.1918.0324.3735.30
20083.566.356.6011.4018.2721.60
20094.836.869.1411.4017.0018.80
20106.868.899.6512.2019.5623.40
20113.565.337.1111.4016.0020.10
20124.835.346.869.1516.5216.50
20133.056.106.617.6112.9519.60
20146.6010.9212.1912.7013.1913.40
20156.809.9112.4513.2021.3422.90
20164.327.628.389.9014.2117.80
20175.848.899.1413.0020.0722.40
20183.004.576.007.409.8913.20
20194.405.607.4014.6017.4027.60
Table 3. Augmented Dickey–Fuller hypothesis test for precipitation data from stations used in the study area.
Table 3. Augmented Dickey–Fuller hypothesis test for precipitation data from stations used in the study area.
DataAlternative Hypothesis (Significance 0.05)ResultP-Value
2072 (Rainfall Data)10 minStationaryStationary0.0382
20 minStationaryNon stationary0.2505
30 minStationaryNon stationary0.5175
60 minStationaryNon stationary0.5984
120 minStationaryNon stationary0.7600
180 minStationaryNon stationary0.6485
2001StationaryNon stationary0.2235
2005StationaryNon stationary0.2138
2035StationaryNon stationary0.5035
2036StationaryNon stationary0.1968
2045StationaryNon stationary0.0870
2065StationaryNon stationary0.4364
2072StationaryNon stationary0.3331
2079StationaryStationary0.0380
2106StationaryNon stationary0.0603
2108StationaryNon stationary0.0535
Table 4. Literature review of ( P 1 2 ) for California and Baja California.
Table 4. Literature review of ( P 1 2 ) for California and Baja California.
Location of Weather StationsRain ( P 1 2 )
(mm)
Reference
Southern California12.7Frevert et al., 1963 [37]
Southern California12.7Reich, 1963 [27]
Ensenada B.C14Reich, 1963 [27]
San Diego14.22Dedrick et al., 1976 [27]
Coronado San Diego, CA13.13
Ranges (11.48–16.61)
Hodges et al, 1961 [26]
NOAA, 2020 [38]
Chula Vista San Diego, CA11.63NOAA, 2020 [38]
Imperial Beach, CA10.84NOAA, 2020 [38]
San Ysidro, CARanges (9.39–13.63)NOAA, 2020 [38]
Northern of Baja California10CENAPRED, 2016 [22]
Table 5. Depth-duration ratio (R) estimated through the average ratios R1…4 estimated from ( P 1 2 ) and ( P 24 2 ) . Where ( P 1 2 ) 1 was averaged of pluviograph and daily records applying Gumbel, ( P 1 2 ) 2 calculated by Hershfield [26] and Bell [28], ( P 1 2 ) 3 calculated by Reich [27], and ( P 1 2 ) 4 calculated by the adjusting and modification of Bell [28] (Equation (7)).
Table 5. Depth-duration ratio (R) estimated through the average ratios R1…4 estimated from ( P 1 2 ) and ( P 24 2 ) . Where ( P 1 2 ) 1 was averaged of pluviograph and daily records applying Gumbel, ( P 1 2 ) 2 calculated by Hershfield [26] and Bell [28], ( P 1 2 ) 3 calculated by Reich [27], and ( P 1 2 ) 4 calculated by the adjusting and modification of Bell [28] (Equation (7)).
Station ( P 1 2 ) 1
(mm)
( P 1 2 ) 2
(mm)
( P 1 2 ) 3
(mm)
( P 1 2 ) 4
(mm)
( P 24 2 )
(mm)
R1R2R3R4R
207210.5617.3514.2012.2932.900.320.530.430.370.41
203611.0021.1313.8014.9238.400.290.550.360.390.40
203510.9617.1314.5012.0931.200.350.550.470.390.44
200110.2019.4013.9011.6234.700.290.560.400.330.40
200510.2023.3313.7012.8641.000.250.570.330.310.37
204510.2022.3814.5013.0239.900.260.560.360.330.38
206510.2021.8714.5011.0540.700.250.540.360.270.35
207910.2020.3314.7012.9635.100.290.580.420.370.41
210610.2018.3414.5010.8034.000.300.540.430.320.40
210810.2021.2214.510.1639.700.260.530.360.260.35
Table 6. First comparison of IDF relationship, estimated through different methods for station 2072 for different return periods.
Table 6. First comparison of IDF relationship, estimated through different methods for station 2072 for different return periods.
Return PeriodIDF (mm/hr) Station 2072
Method10 (min)20 (min)30 (min)60 (min)120 (min)
10 years(A)38.328.122.314.610.5
(B)45.232.225.817.211.2
(C)48.038.026.016.013.0
(D)52.837.830.119.712.5
(E)59.642.733.921.813.5
(F)70.650.840.125.114.9
(G)99.259.043.525.915.4
Average59.141.231.720.013.0
25 years(A)44.332.425.616.712.0
(B)53.338.030.520.313.2
(C)57.040.031.020.015.0
(D)62.444.635.523.214.7
(E)70.350.440.125.715.9
(F)83.460.047.329.717.6
(G)119.070.052.031.018.0
Average70.047.937.423.815.2
50 years(A)48.735.528.018.213.1
(B)59.542.534.022.614.7
(C)64.045.035.022.018.0
(D)69.649.739.625.916.4
(E)78.556.344.728.717.7
(F)93.066.952.833.119.7
(G)134.780.159.135.120.9
Average78.353.741.925.817.0
100 years(A)53.138.730.419.714.2
(B)65.746.937.525.016.3
(C)71.050.038.024.020.0
(D)76.854.943.828.618.1
(E)86.662.149.331.719.6
(F)102.773.958.336.621.7
(G)150.089.265.839.123.3
Average86.659.446.229.219.0
(A) Estimated from pluviograph records using Gumbel (Equation (1)). (B) Estimated from the Chen equation using R1 from Table 5. (C) Obtained from the isohyetal map of the SCT. (D) Estimated from the Chen equation using R4 from Table 5 (Equation (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.
Table 7. IDF relationship selected for the rain gauge stations estimated through the Chen Method with the (R) ratio proposed in this study.
Table 7. IDF relationship selected for the rain gauge stations estimated through the Chen Method with the (R) ratio proposed in this study.
Return PeriodStationsIntensity Duration Frequency (mm/h) for Different Durations
10 (min)20 (min)30 (min)60 (min)120 (min)
10 years207257.7141.3532.8621.2213.21
203676.7054.9243.6828.3017.69
203563.5845.6336.1923.1714.22
200168.6849.1839.1125.3415.84
200576.6154.7643.6628.5518.12
204566.5947.6237.9324.7315.61
206568.5848.9739.1125.7616.51
207971.3251.1040.6126.2316.32
210661.6044.1135.0822.7314.21
210874.8453.4442.6828.1118.02
25 years207268.1548.8338.8025.0615.59
203693.9667.2853.5134.6621.67
203576.8255.1443.7328.0017.18
200183.9560.1247.8130.9719.36
200594.0067.1953.5635.0322.23
204579.0056.5145.0129.3418.53
206583.2859.4747.5031.2820.05
207987.2862.5449.7032.1019.97
210673.7452.8141.9927.2117.01
210893.0966.4853.0934.9722.42
50 years207276.0454.4843.3027.9617.40
2036107.0176.6360.9439.4824.68
203586.8562.3349.4331.6519.43
200195.5068.3954.3835.2322.03
2005107.1576.5961.0639.9425.34
204588.3963.2250.3632.8320.73
206594.4167.4153.8435.4622.73
207999.3671.1956.5736.5422.74
210682.9259.3847.2230.5919.13
2108106.9076.3460.9740.1525.74
100 years207283.9360.1447.7930.8719.21
2036120.0785.9868.3844.3027.70
203596.8769.5255.1335.3021.67
2001107.0576.6660.9639.4924.69
2005120.3085.9968.5544.8428.45
204597.7969.9455.7136.3222.93
2065105.5375.3560.1839.6325.41
2079111.4379.8463.4540.9825.50
210692.1165.9652.4533.9821.25
2108120.7186.2068.8445.3429.07
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Gámez-Balmaceda, E.; López-Ramos, A.; Martínez-Acosta, L.; Medrano-Barboza, J.P.; Remolina López, J.F.; Seingier, G.; Daesslé, L.W.; López-Lambraño, A.A. Rainfall Intensity-Duration-Frequency Relationship. Case Study: Depth-Duration Ratio in a Semi-Arid Zone in Mexico. Hydrology 2020, 7, 78. https://doi.org/10.3390/hydrology7040078

AMA Style

Gámez-Balmaceda E, López-Ramos A, Martínez-Acosta L, Medrano-Barboza JP, Remolina López JF, Seingier G, Daesslé LW, López-Lambraño AA. Rainfall Intensity-Duration-Frequency Relationship. Case Study: Depth-Duration Ratio in a Semi-Arid Zone in Mexico. Hydrology. 2020; 7(4):78. https://doi.org/10.3390/hydrology7040078

Chicago/Turabian Style

Gámez-Balmaceda, Ena, Alvaro López-Ramos, Luisa Martínez-Acosta, Juan Pablo Medrano-Barboza, John Freddy Remolina López, Georges Seingier, Luis Walter Daesslé, and Alvaro Alberto López-Lambraño. 2020. "Rainfall Intensity-Duration-Frequency Relationship. Case Study: Depth-Duration Ratio in a Semi-Arid Zone in Mexico" Hydrology 7, no. 4: 78. https://doi.org/10.3390/hydrology7040078

APA Style

Gámez-Balmaceda, E., López-Ramos, A., Martínez-Acosta, L., Medrano-Barboza, J. P., Remolina López, J. F., Seingier, G., Daesslé, L. W., & López-Lambraño, A. A. (2020). Rainfall Intensity-Duration-Frequency Relationship. Case Study: Depth-Duration Ratio in a Semi-Arid Zone in Mexico. Hydrology, 7(4), 78. https://doi.org/10.3390/hydrology7040078

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop