Adequacy of methodologies for determining SCS / CN in a watershed with characteristics of the Pampa biome

The Soil Conservation Service Curve Number Model is a conceptual model intended for estimating effective rainfall (ER). This model is grounded in a parameter – referred to as Curve Number (CN), which is determined from information on the characteristics of the watershed. The Standard Method (M1) for determining the CN is based on soil and land-use tables; however, some authors have proposed alternative methodologies for defining the CN value from monitored rainfall-runoff events, such as those described by Hawkins (1993) (M2), Soulis and Valiantzas (2012) (M3), and Soulis and Valiantzas (2013) (M4). The objective of this study was to evaluate the impact of using these methods for determination of the CN parameter on the estimation of ER, taking as reference forty rainfall-runoff events monitored between 2015 and 2018 in the Cadeia River Watershed, which has characteristics of the Pampa biome. The different methods assessed for definition of the CN parameter resulted in contrasting performances with respect to the estimation of ER for CRW, as the following findings: i) M1 gave ER values with little reliability, mainly due to the classification of antecedent moisture content classes; ii) M3 provided the best results in determining ER, followed by M2; and iii) the ER values estimated according to M4 differed from those observed, mainly for events with lower rainfall depths.


INTRODUCTION
Knowledge of the hydrological cycle is essential for the adequate management of water resources. In this sense, direct surface runoff (DSR) is one of the most important components of the hydrological cycle, especially when considering practical issues involving hydrological engineering (Alves et al., 2019).
The generation of DSR is dependent on variables associated with the rainfall event, such as its intensity and duration, as well as characteristics of the watershed, e.g. soil classes, land uses, and topography (Elhakeem and Papanicolaou, 2009;Lal et al., 2016). Thus, according to Araújo Neto et al. (2012) and Soulis and Valiantzas (2013), the use of hydrological models to estimate DSR in watersheds is necessary.
Among the existing models, the Soil Conservation Service Curve Number (SCS/CN) (NRCS, 1986) is the simplest conceptual model for estimation of DSR (Ajmal and Kim, 2014), which is numerically equal to effective rainfall (ER). According to Ponce and Hawkins (1996), the SCS/CN model aims to estimate the DSR volume yielded during a rainfall event based on a parameter known as the Curve Number (CN). Although SCS/CN model has been designed for agricultural watersheds in the United States, the model has been widely applied in locations with different edaphoclimatic characteristics (Ajmal and Kim, 2014;Oliveira et al., 2016;Walega et al., 2015;Endale et al., 2015;D'Asaro et al., 2014;Deshmukh et al., 2013). Given its simplicity, the SCS/CN model and its adaptations are often coupled in more complex hydrological models, e.g., models originated from the Unit Hydrograph (HU) theory, the Soil and Water Assessment Tool (SWAT) (Arnold et al., 1998) and the Lavras Simulation of Hydrology (LASH) (Beskow et al., 2011).
The application of the SCS/CN model is greatly impacted by the relationship established between the initial rainfall abstraction (Ia) and the maximum potential for water infiltration in the soil (S), expressed by the coefficient of initial rainfall abstractions (λ). Ponce and Hawkins (1996) reported that λ depends on the relief, vegetation, and climatological characteristics of the watershed. Due to the direct influence of λ on the estimates of ER, the definition of its value has been discussed by several authors (Alves et al., 2019;Mishra and Singh., 2006;Valle Junior et al., 2019).
It is known that the SCS/CN model requires the CN parameter, which can be calculated for each rainfall-runoff event by using data on ER and Ia, as in the study by Alves et al. (2019). However, when there are limitations or absence of observed data, CN is determined based on the tabulated CN values included in the method documentation for various soil/land use complexes, without considering the characteristics of the rainfall event. Furthermore, the determination of this parameter is subjective when considering the effect of antecedent soil 3 Adequacy of methodologies for determining SCS … Rev. Ambient. Água vol. 16 n. 4, e2715 -Taubaté 2021 moisture content conditions (Soulis and Valiantzas, 2012). Therefore, the practicality of assigning a constant CN parameter for the estimation of DSR, based only on physical characteristics of the watershed, is subjective and susceptible to errors. This is expected to occur mainly when applied to watersheds with different physiographic characteristics from those used for the development of the SCS/CN model, such as the watershed assessed in this study, situated in the Pampa biome in Southern Brazil. According to Lupatini et al. (2013), the Pampa biome covers part of Southern Brazil, has distinct characteristics of vegetation, climate, and soil classes, and is characterized by a great plant and animal diversity.
According to Soulis (2018) in cases of watersheds with different characteristics, it is more appropriate to estimate the CN values based on the analysis of actual hydrological observations, comparing the hydrological response in different conditions. In this context, several methods have been developed and recommended for the determination of the CN parameter using measured rainfall-runoff data, considering variables of influence on the generation of DSD, e.g., Hawkins (1993), Soulis and Valiantzas (2012), and Soulis and Valiantzas (2013). Nevertheless, little is known so far about the coherence of the CN parameter derived from these methods with the reality of each watershed, thereby making it necessary to compare the results of these methods with the values obtained from observed rainfall-runoff data.
The aforementioned methods have been frequently applied for modelling of rainfall events with high intensities and short durations (less than one day), that is, in consonance with the original development of the SCS/CN model (Ajmal et al., 2015;Mishra et al., 2008;Endale et al., 2015;D'Asaro et al., 2014;Deshmukh et al., 2013). Mishra et al. (2008) stress that the generation of DSR in a watershed is also dependent on the rainfall duration, as losses by infiltration and evapotranspiration are directly related to the time the water remains on the surface of the watershed. Consequently, it is necessary to investigate the application of the SCS/CN model for rainfall events with long durations, verifying the variability of the CN parameter over the duration of the event (Žlábek et al., 2012;Mishra et al., 2008).
In view of the above, the objective of this study was to evaluate the impact of using different methodological proposals for determination of the CN parameter on the estimation of ER according the SCS/CN model, taking as reference a watershed with typical characteristics of the Pampa biome and with rainfall durations lasting more than one day.

Characterization of the study area
The Cadeia River Watershed (CRW) was evaluated in this study. CRW is located in the south of Rio Grande do Sul (Figure 1), has a drainage area of 121.3 km² and is one of the main sub-watersheds of the Pelotas River Watershed (PRW). The Research Group on Hydrology and Hydrological Modeling in Watersheds is responsible for a hydrological monitoring network installed in the PRW (Figure 1). In this work, rainfall data were obtained from seven recording rain gauges and streamflow data were acquired from 1 hydrological gauging station ( Figure 1a).
CRW is inserted in the Pampa biome, in the southeast region of Rio Grande do Sul State/Brazil. This region is characterized by a subtropical climate with great influence from polar systems (Kuplich et al., 2018). According to the Köppen classification, the climate of the region is characterized as Cfa, with an average temperature of the warmest month exceeding 22ºC and well-distributed monthly rainfall (Alvares et al., 2014). The average annual rainfall depth in CRW is 1,367 mm (Steinmetz et al., 2019).
The Pampa biome is characterized by typical vegetation, such as natural grassland and sparse formations of shrubs and trees (Overbeck et al., 2007), and by a wide variety of soil types (Lupatini et al., 2013). Information on the soil classes existing in the CRW were extracted from the soil map of the Conselho Regional de Desenvolvimento da Região Sul (COREDE-SUL), prepared by the Laboratório de Planejamento Ambiental of the Empresa Brasileira de Pesquisa Agropecuária (EMBRAPA/Clima Temperado) taking as reference the study conducted by Cunha et al. (2006). The study area is characterized by the presence of Ultisols, according to Soil Taxonomy classification (Soil Survey Staff, 2014). This corresponds to the Yellow Argisol (78.7%), Red Argisol (14.7%) and Grayish Brown Argisol (6.6%) (Figure 1b), according to the Brazilian Soil Classification System -SiBCS (Santos et al., 2018). According to Beskow et al. (2018), the land-use classes identified in CRW are forest (31.0%), grassland (30.0%), bare soil (25.0%), pasture (13.9%), and water body (0.1%) (Figure 1c).

Soil Conservation Service/Curve Number Model (SCS/CN)
The present study was carried out using the SCS/CN Model developed by the Soil Conservation Service (NRCS, 1986) of the United States Department of Agriculture (USDA), which has been widely used for estimating ER in watersheds. According to Cao et al. (2011), this model considers that the ER resulting from a given rainfall event is a function of both rainfall depth (Rtotal) and initial abstraction (Ia). The Ia is associated with the initial losses that occur due to the processes of interception by land cover, infiltration, evaporation, and water retention by depressions of the terrain. SCS/CN is represented by Equations 1 and 2 and assumes there is a relationship between Ia and S (maximum potential retention) (Equation 3) 5 Adequacy of methodologies for determining SCS … Rev. Ambient. Água vol. 16 n. 4, e2715 -Taubaté 2021 expressed by λ (initial abstraction ratio) (NRCS, 1986). On investigating the Ia values in agricultural watersheds in the United States, (NRCS, 1986) stated that the values of Ia corresponded to 20% of S on average, in other words, λ was equal to 0.2.
Where: ER is the effective rainfall (mm), Rtotal refers to the total rainfall depth (mm), Ia is the initial abstraction (mm), S corresponds to the maximum potential for water retention in the soil (mm), and CN is the Curve Number parameter.

Methods for determining the CN parameter
The values of the CN parameter were defined according to four methods analyzed in this study to evaluate SCS/CN model taking as reference rainfall-runoff events observed in the CRW: i) estimating CN values from tabulated values (NRSCS, 1986) based on the physiographic characteristics of the watershed (treating the water as ungauged) (M1), combined with the physiographic characteristics of the watershed under analysis; ii) asymptotic CN (Hawkins, 1993) (M2); iii) Two CN (Soulis and Valiantzas, 2012) (M3); and iv) Heterogeneous CN (Soulis and Valiantzas, 2013) (M4). For comparative purposes, the CN for each rainfallrunoff event was also calculated and considered as observed CN (CNobs).

Standard tables for CN (NRCS, 1986) -(M1)
The soils were classified in accordance with the hydrologic soil groups established by the SCS/CN model, following the classification proposed by Sartori et al. (2005) for Brazilian soils. Information on the potential of a specific soil type for generation of DSR, e.g., texture, depth and transition between horizons, was used for the aforementioned classification (Sartori et al., 2005). Therefore, the Red Argisol was classified as Group B, as it is characterized by a moderate depth and infiltration rate, a moderate to good drainage, and a moderately fine to moderately coarse texture. The Yellow Argisol was classified as Group C due to its low infiltration rate, fine texture, and layers that hinder the vertical movement of water. Finally, the Grayish Brown Argisol was defined as Group D because it has a high potential for DSR generation, with a very low infiltration rate.
After the classification of soil types, they were merged with the land-use classes found in the watershed in order to determine the CN value corresponding to each combination, considering the intermediate conditions of antecedent moisture content (AMCII). The CN values (CNII) used in this study for the different combinations of hydrologic soil groups and land-use classes are presented in Table 1. Afterwards, the CNII values were spatialized in the watershed, weighted by the areas of occurrence of each value, to obtain an average CNII value.
However, according to the method, the actual antecedent moisture content is related to the rainfall depth that occurred in the last five days (R5), so it must be evaluated for each event. Following the framework proposed by ASCE (1996) referred to as Antecedent Moisture Content (AMC), the CN values were grouped into three conditions: AMC I (dry soil and R5 up to 35 mm), AMC II (moist soil corresponding to field capacity and R5 between 35 and 52.5 mm) and AMC III (very moist soil and R5 greater than 52.5 mm). In this sense, the average CNII value was corrected for the AMC class according to the characteristics of each rainfall event assessed in this study, using Equations 4 and 5, as presented by Chow et al. (1988).
Where: CNI and CNIII are the corrected CN values corresponding to the classes AMC I and AMC III, respectively.

Asymptotic Determination Method of runoff Curve Numbers (Hawkins, 1993) -(M2)
The M2 was developed by Hawkins (1993) and allows the determination of CN values from observed rainfall-runoff events, using the frequency matching concept. First, the observed values of Rtotal and ER were independently classified in decreasing order, forming Rtotal:ER pairs with equal return periods (Hawkins, 1993). Subsequently, the CN values were determined for each Rtotal:ER pair found using Equation 6 (Hawkins 1993).
Where: ER is the effective rainfall (mm), Rtotal refers to the rainfall depth (mm), λ is the initial rainfall abstraction ratio, equal to 0.2, and CN corresponds to the Curve Number.
The CN values determined in accordance with Equation 6 were analyzed in order to classify the behavior of CRW with respect to the observed CN values and Rtotal. The following behaviors can be identified (Hawkins 1993 Where: CN(P) is the CN value adjusted for Rtotal; CN∞ corresponds to the constant value of CN when Rtotal tends to infinity; k1 and k2 refer to coefficients adjusted by means of the leastsquares method.

Two-CN System Model Approach Method (Soulis and Valiantzas 2012) -(M3)
Soulis and Valiantzas (2012) proposed the Two-CN System Model Approach method, which aims to divide the watershed into two homogeneous areas characterized by relatively similar land uses. According to Walega et al. (2015), M3 is based on the hypothesis that, during a rainfall event, DSR is initially yielded in areas with less permeability, that is, with higher CN values. Thus, as the event continues, DSR also begins to occur in areas with greater permeability. Therefore, two CN values are attributed to these areas for the application of M3 -CNa and CNb, in which CNa > CNb. If α represents the fraction of the watershed corresponding to CNa, then (1-α) is the area represented by CNb. The procedures recommended by Soulis and Valiantzas (2012) were performed to apply M3 as follows: 1) The values of Rtotal and ER were independently classified in decreasing order, forming pairs of values of Rtotal and ER with equal return periods; 2) The CN values obtained by M1 for CRW were related to their respective areas of occurrence within the watershed and organized in a decreasing order to determine the CNa, CNb, and α parameters. Then, considering the highest values of CN related to their respective areas of occurrence, a weighted average CN was calculated, referring to the initial value of CNa. Likewise, the weighted average CN value calculated for the areas of lower CN values corresponded to the initial value for CNb.
3) For the initial CNa and CNb values, the values of S and Ia were determined from Equations 2 and 3, respectively. 4) Equations 9, 10, and 11 were applied to calculate ER values, considering the observed Rtotal values. Subsequently, the CN values were calculated for each event using Equation 6. 5) Finally, the parameters CNa, CNb, and α were adjusted by using statistical optimization in an attempt to approximate the CN values calculated in Step 1 to those calculated in Step. The objective function used to optimize the parameters in CRW was the Root Mean Square Error (RMSE).

Heterogeneous CN Method (Soulis and Valiantzas, 2013) -(M4)
This method is based on M3; however, M4 extends the number of CN categories to the actual number of combinations existing between soil classes and their land uses, i.e., the combinations obtained in M1, assuming that the DSR is the partial sum of ERi computed from these CN values. For the application of M4, Equations 12, 13 and 14 is rewritten (Soulis and Valiantzas 2013): The following steps were performed (Soulis and Valiantzas, 2013) to obtain the CNi values for CRW: 1) The values of Rtotal and ER were ordered as in M2; 2) From Equation 6, the CN values were calculated for each pair formed by Rtotal and ER; 3) The CN values obtained in M1 and their respective areas of occurrence within the CRW were used as initial values of the parameters CNi and ai; 4) ERi values were calculated for each Rtotal value considering the conditions described in 6) The values of CNi and ai were readjusted through the optimization of Equation 12; seeking to approximate the values of CN calculated in Step 2 to those determined in Step 5. As in M3, RMSE was an objective function used for the optimization of the parameters.

Selection and characterization of rainfall-runoff events
A stage-discharge rating curve developed for the CRW's outlet was used to convert the water level data monitored at its outlet into streamflow values. The mean rainfall was calculated according to the Thiessen Polygons method (Thiessen and Alter, 1911). Thus, the observed hydrographs and the mean rainfall hyetographs in the CRW were determined for a 30-minute interval. This time interval is considered satisfactory for hydrological modeling in the CRW since it has an estimated time of concentration of 7.5 hours (Moura et al., 2021). Forty rainfall-runoff events between 2015 and 2018 were selected, which were then represented by a hyetograph of mean rainfall and a hydrograph. The runoff separation was performed based on the hydrograph of each event by means of a graphical method, as described by Chow et al. (1988). Subsequently, the following characteristics of each event were determined: rainfall depth (Rtotal), 5-day antecedent rainfall depth (R5), maximum rainfall intensity in 30 minutes (Imáx), initial abstraction (Ia), effective rainfall (ER) and observed CN (CNobs).
ER was obtained for each observed event using the SCS/CN model (NRCS, 1986), and Ia was defined as the total rainfall depth that occurred between the beginning of rainfall and the beginning of DSR, following the methodology described in Mishra and Singh (2006). R5 was computed based on the start date and time of each event.

Statistical analysis
The following evaluation criteria were used to verify the accuracy of the methods for definition of CN and, consequently, estimation of ER: relative error of each event (Re) (%) (Equation 15), Pbias coefficient (%) (Equation 16) and Nash-Sutcliffe coefficient (NS) (Nash and Sutcliffe, 1970) (Equation 17).
Where: Obsi represents the observed ER value for an event i, Esti refers to the estimated ER value for an event i and Obs ̅̅̅̅̅ i is the average observed ER value.

Rainfall-runoff events
The main characteristics of the rainfall-runoff events assessed in this study are presented in Table 2. Rainfall duration ranged from 3.5 to 63.0h, with an average value equal to 26.0h. On analyzing Rtotal and ER (Table 2), one can verify a noticeable variation in their values, which implies large variations in the CNobs values.

CN values obtained from M1, M2, M3, and M4 methods
Twelve CNII values, varying between 58 and 100, were determined when using M1 ( Figure  2a), resulting in average values of 74, 55 and 87 for, respectively, CNII, CNI and CNIII. For M3, the values found for the parameters CNa, CNb, and α were 87, 25, and 0.2814, respectively (Figure 2b). In the case of M4, 7 CNi values were derived. The values of CN and the area fractions of the CRW for M1, M3 and M4 are shown in Table 3. The spatial distribution of CN values obtained from M1, M3, and M4 is illustrated in Figure 2. Out of the three possible behaviors for M2, CRW fell within the standard behavior, as the calculated CN(P) values decreased as Rtotal increased, tending to a constant value (Figure 2d).   and c) M4, d) behavior of CRW with respect to M2, according to Hawkins 1993, and e) Asymptotic curve according to Hawkins (1993) in comparison with the CNobs values.
The parameters CN∞ (57) Where: CN(P) is the CN value adjusted according to the Rtotal.
The CN values determined from the methods applied in this study differed from the CNobs values. Some reasons for that must be highlighted: i) M1, M3, and M4 require information on soil classes and land uses to estimate the CN parameter, while CNobs is determined from observed rainfall-runoff data; ii) the characteristics of the watersheds used for the development of the SCS/CN model differ from those existing in the CRW; iii) Rtotal: ER pairs formed for the M2, M3 and M4 methods may not be representative for the reality of the events that occur in CRW; iv) it is difficult to identify the direct surface runoff in a watershed with such a large area, such as CRW, given the spatial and temporal variation of rainfall and other processes.
The low accuracy of M1 for determination of CN parameters in the CRW can be explained by the fact that the CN parameter does not depend only on the soil classes, land uses, and soil moisture, as recommended by M1. According to Kowalik and Walega (2015), the CN parameter is also affected by rainfall, which is not considered in M1. Besides, it is also expected to have lower performance when the watershed is treated as ungauged. Several studies that applied M1 to estimate the CN parameter in watersheds reported that its values obtained from the standard tables differed significantly from those empirically computed from observed data (Kowalik and Walega, 2015). In this sense, authors from different regions of the world have investigated and reported better estimates of CN values for watersheds using methods based on the calibration of observed data (D'Asaro et al., 2014;Ajmal and Kim, 2014;Alves et al., 2019).

Performance analysis
One can analyze in Table 4    It can be inferred from Table 4 that M1 did not generate reliable results for ER in terms of performance statistics. Pbias indicated a tendency to overestimate ER values (-0.4%) and the NS value (-0.8) reinforced the lack of accuracy of M1 for determining ER. The performance statistics for M2 indicated a better adjustment when compared to M1, given that Pbias and NS were equal to 0.0 and 0.6, respectively. Analyzing the NS values of the methods to estimate ER, M3 outperformed the other methods (NS = 0.7); however, the Pbias value indicated an underestimation of ER. The values of statistics in the estimation of ER according to M4 demonstrated that the method was not accurate, but this method resulted in a NS greater than that generated by M1, which presented the worst adjustment. On the other hand, analyzing the Pbias coefficient, M4 presented a high tendency of underestimations in the determination of ER, unlike the methods M1 and M2.
Based on the Rtotal values obtained for each event, there was a relationship between the magnitude of the event and the percentage differences between the observed ER values and those estimated by M1. M1 tended to overestimate ER values for events with higher Rtotal values and underestimate ER with lower Rtotal values. A similar behavior occurred in the study by Cunha et al. (2015), where the authors explained that overestimations in ER values may be associated with the reduction in the CNobs value as Rtotal value increases. On the other hand, underestimations in ER occur because most of the Rtotal of the event is considered as part of Ia values (Cunha et al., 2015). Another factor of influence in the underestimates was the classification of the AMC conditions, as 22 out of the 28 underestimated events were classified as AMC I, five events as AMC II, and only one event as AMC III. AMC I assumes that soil is very dry, resulting in high values of Ia (Shaw and Walter, 2009).
The standard behavior observed in the events from CRW for M2 was similar to that found in the watershed analyzed by Hawkins (1993). In addition, the CN∞ obtained for M2 was considerably different from the average CN obtained for M1. This difference observed for CRW may be due to the use of rainfall-runoff events from different rainfall regimes, that is, both convective rainfall and frontal rainfall. M2 produced results with better performances in estimating ER values when compared to M1. According to Alves et al. (2019), among the advantages of using an asymptotic approach to determine the CN parameter, it should be mentioned the fact that this method takes into account all events in the data set.
The parameters found for M3 in CRW indicate that the behavior of DSR in the watershed is more influenced by the characteristics of the rainfall, mainly by its intensity and duration. In this sense, even if the CRW has a larger area with permeable surfaces, Mishra and Singh (2006) underlines that higher rainfall intensities reduce the time water remains on the soil, thus resulting in a higher generation of DSR. In M3, the parameter α corresponds to CNa representing the portion of the watershed characterized by low permeability capacity ; that is, where DSR is formed faster after the beginning of the rainfall. Therefore, analyzing the low calibrated value for α (0.2814) in M3 for CRW, one can infer that CRW has a larger area with permeable soils. According to Walega et al. (2015), the division of the watershed into two relatively homogeneous areas in M3 allows a very accurate estimate of the amount of ER. Walega et al. (2015) point out that, differently from M1 and M2, M3 considers the impact of spatial variability of the characteristics of the watershed.
The ER values obtained from M4 differed from those observed, mainly for events characterized by smaller Rtotal. This behavior was also verified in the studies by Alves et al. (2019), Walega et al. (2015), and Ajmal and Kim (2014). According to Alves et al. (2019), M4 indicates that the ER generated in events with lower Rtotal is associated with areas of the watershed represented by higher CN values, i.e, in areas with greater difficulty of infiltration, such as urban areas, roads, and compacted soils. On the other hand, M4 did not perform well to estimate ER for events with greater magnitude, even though the CN values determined for M4 have been similar to the CNobs values. This fact can be justified by the value of λ considered (0.2), which implies very high Ia values and less ER generation. In addition, the ER values calculated according to M4 have a direct impact on the use of the CN values resulting from M1 to calibrate its parameters.
Based on the ER values derived from M1, M2, M3 and M4, it can be observed that the results found corroborate those obtained by other authors. Walega et al. (2015) concluded in their study that M3 was more accurate when describing the observed ER values. Ajmal et al. (2015) observed that the CN values estimated from observed rainfall-runoff data provided better results than the values determined from standard tables of the SCS/CN model.
The results found in the ER estimate by the applied methods can be elucidated by the value used for λ. When analyzing different methods for estimating CN, some authors found that the parameter λ varied considerably depending on the watershed analyzed. Beskow et al. (2009) calibrated a value of λ equal to 0.105 for a watershed located in the south of Minas Gerais State/Brazil. In the study by Caldeira et al. (2019), the calibration of the LASH hydrological model for the Fragata River Watershed, located near CRW, resulted in values of λ of 0.147 (distributed modeling) and 0.039 (semi-distributed modeling). In this sense, it appears that it is not advisable to use the fixed value recommended by the SCS/CN Model for λ, since the variations in the value of the parameter λ in the watersheds occur due to differences in climatic conditions, land use, and antecedent moisture content (Beskow et al., 2011).
It should be noted that the SCS/CN model considers rainfall events of short duration and, therefore, without significant variation in the maximum potential for water infiltration in the soil during the event. However, the events that occur in watersheds with predominance of frontal systems, such as CRW, are long-lasting and experience varying intensities over time (Beskow et al., 2018). Thus, the maximum potential for water infiltration in the soil may vary during the event, given the variability of losses by evapotranspiration and infiltration that significantly depend on the time that water remains on the soil surface (Mishra et al., 2008).

CONCLUSIONS
The different methods evaluated in this study for obtaining the CN parameter resulted in contrasting performances with respect to the estimation of ER for CRW. Based on the results obtained, the following conclusions can be drawn: • The CN values estimated derived from the four applied methods produced significant differences compared to the CNobs values.
• M1 resulted in ER values with little reliability for CRW, mainly due to the classification of soil AMC classes, such that, in general, the events with CN values classified under AMC I resulted in underestimations of ER values for CRW.
• The CN values obtained from M2 were not close to the CNobs values since the hydrological reality of CRW does not present a good relationship between Rtotal and ER. In this sense, for watersheds with this behavior, it would be interesting to evaluate the influence of the temporal distribution of the rainfall event on the variation of S for modeling ER.
• M2 and M3 provided the most accurate results in determining ER for CRW, with better performance for M3, which had CNa and CNb parameters with very clear physical meaning for this watershed.