Using a topographic index to distribute variable source area runoff predicted with the SCS curve‐number equation

Because the traditional Soil Conservation Service curve‐number (SCS‐CN) approach continues to be used ubiquitously in water quality models, new application methods are needed that are consistent with variable source area (VSA) hydrological processes in the landscape. We developed and tested a distributed approach for applying the traditional SCS‐CN equation to watersheds where VSA hydrology is a dominant process. Predicting the location of source areas is important for watershed planning because restricting potentially polluting activities from runoff source areas is fundamental to controlling non‐point‐source pollution. The method presented here used the traditional SCS‐CN approach to predict runoff volume and spatial extent of saturated areas and a topographic index, like that used in TOPMODEL, to distribute runoff source areas through watersheds. The resulting distributed CN–VSA method was applied to two subwatersheds of the Delaware basin in the Catskill Mountains region of New York State and one watershed in south‐eastern Australia to produce runoff‐probability maps. Observed saturated area locations in the watersheds agreed with the distributed CN–VSA method. Results showed good agreement with those obtained from the previously validated soil moisture routing (SMR) model. When compared with the traditional SCS‐CN method, the distributed CN–VSA method predicted a similar total volume of runoff, but vastly different locations of runoff generation. Thus, the distributed CN–VSA approach provides a physically based method that is simple enough to be incorporated into water quality models, and other tools that currently use the traditional SCS–CN method, while still adhering to the principles of VSA hydrology. Copyright © 2004 John Wiley & Sons, Ltd.     

A combined rainfall infiltration model based on Green–Ampt and SCS‐curve number

The Green–Ampt infiltration equation is an incomplete governing equation for rainfall infiltration due to the absence of an inertia term. The estimation of the capillary pressure head at the wetting front is difficult to determine. Thus, a major limitation of the Green–Ampt model is the constant, non‐zero surface ponding depth. This paper proposes an integrated rainfall infiltration model based on the Green–Ampt model and the SCS‐CN model. It achieves a complete governing equation for rainfall infiltration by momentum balance and the water budget based on the Green–Ampt assumption, and uses the curve number from the SCS‐CN method to calculate the initial abstraction, which is used as a basic parameter for the governing equation of the intensity of rainfall loss during the runoff period. The integrated rainfall infiltration model resolves the dilemma for capillary pressure head estimation, overcomes the limitation of constant, non‐zero surface ponding depth, and facilitates the calculation of runoff for individual flood simulations. Copyright © 2014 John Wiley & Sons, Ltd.     

Green‐Ampt Curve‐Number mixed procedure as an empirical tool for rainfall–runoff modelling in small and ungauged basins

The Soil Conservation Service Curve Number (SCS‐CN) method is a popular rainfall–runoff model that is widely used to estimate direct runoff from small and ungauged basins. The SCS‐CN is a simple and valuable approach to quantify the total streamflow volume generated by storm rainfall, but its use is not appropriate for estimating the sub‐daily incremental rainfall excess. To overcome this drawback, we propose to include the Green‐Ampt (GA) infiltration model into a mixed procedure, which is referred to as Curve Number for Green‐Ampt (CN4GA), aiming to distribute in time the information provided by the SCS‐CN method. For a given storm, the computed SCS‐CN total net rainfall amount is employed to calibrate the soil hydraulic conductivity parameter of the GA model. The proposed procedure is evaluated by analysing 100 rainfall–runoff events that were observed in four small catchments of varying size. CN4GA appears to provide encouraging results for predicting the net rainfall peak and duration values and has shown, at least for the test cases considered in this study, better agreement with the observed hydrographs than the classic SCS‐CN method. Copyright © 2012 John Wiley & Sons, Ltd.     

LAND-USE IMPACT ON WATERSHED RESPONSE: THE INTEGRATION OF TWO-DIMENSIONAL HYDROLOGICAL MODELLING AND GEOGRAPHICAL INFORMATION SYSTEMS

The integration of a two-dimensional, raster-based rainfall–runoff model, CASC2D, with a raster geographical information system (GIS), GRASS, offers enhanced capabilities for analysing the hydrological impact under a variety of land management scenarios. The spatially varied components of the watershed, such as slope, soil texture, surface roughness and land-use disturbance, were characterized in GRASS at a user-specified grid cell resolution for input into the CASC2D model. CASC2D is a raster-based, single-event rainfall–runoff model that divides the watershed into grid cell elements and simulates the hydrological processes of infiltration, overland flow and channel flow in response to distributed rainfall precipitation. The five-step integration of CASC2D and GRASS demonstrates the potential for analysing spatially and temporally varied hydrological processes within a 50 square mile semi-arid watershed. By defining possible land-use disturbance scenarios for the watershed, a variety of rainfall–runoff events were simulated to determine the changes in watershed response under varying disturbance and rainfall conditions. Additionally, spatially distributed infiltration outputs derived from the simulations were analysed in GRASS to determine the variability of hydrological change within the watershed. Grid cell computational capabilities in GRASS allow the user to combine the scenario simulation outputs with other distributed watershed parameters to develop complex maps depicting potential areas of hydrological sensitivity. This GIS–hydrological model integration provides valuable spatial information to researchers and managers concerned with the study and effects of land-use on hydrological response.     

Mapping runoff generating areas using AGNPS-VSA model

ABSTRACT

In humid regions, surface runoff is often generated by saturation-excess runoff mechanisms from relatively small variable source areas (VSAs). However, the majority of the current hydrologic models are based on infiltration-excess mechanisms. In this study, the AGricultural Non-Point Source Pollution (AGNPS) model was used to integrate the VSA concept using topographic wetness index (TWI). Both the original and AGNPS-VSA models were evaluated for a small agricultural field in Ontario, Canada. The results indicate that the AGNPS-VSA model performed better than original model. The AGNPS-VSA model predicted that only the saturated portion of the field with higher TWI values produced runoff, whereas the original AGNPS model showed uniform hydrologic response from the entire field. The results of this study are important for accurately mapping the locations of VSAs. This new model could be a powerful tool in identifying critical source areas for applying targeted best management practices to minimize pollutant loads to receiving waters.     

Development and application of a physically based landscape water balance in the SWAT model

Watershed scale hydrological and biogeochemical models rely on the correct spatial‐temporal prediction of processes governing water and contaminant movement. The Soil and Water Assessment Tool (SWAT) model, one of the most commonly used watershed scale models, uses the popular curve number (CN) method to determine the respective amounts of infiltration and surface runoff. Although appropriate for flood forecasting in temperate climates, the CN method has been shown to be less than ideal in many situations (e.g. monsoonal climates and areas dominated by variable source area hydrology). The CN model is based on the assumption that there is a unique relationship between the average moisture content and the CN for all hydrologic response units (HRUs), and that the moisture content distribution is similar for each runoff event, which is not the case in many regions. Presented here is a physically based water balance that was coded in the SWAT model to replace the CN method of runoff generation. To compare this new water balance SWAT (SWAT‐WB) to the original CN‐based SWAT (SWAT‐CN), two watersheds were initialized; one in the headwaters of the Blue Nile in Ethiopia and one in the Catskill Mountains of New York. In the Ethiopian watershed, streamflow predictions were better using SWAT‐WB than SWAT‐CN [Nash–Sutcliffe efficiencies (NSE) of 0·79 and 0·67, respectively]. In the temperate Catskills, SWAT‐WB and SWAT‐CN predictions were approximately equivalent (NSE > 0·70). The spatial distribution of runoff‐generating areas differed greatly between the two models, with SWAT‐WB reflecting the topographical controls imposed on the model. Results show that a water balance provides results equal to or better than the CN, but with a more physically based approach. Copyright © 2010 John Wiley & Sons, Ltd.     

Curve‐Number/Green–Ampt mixed procedure for streamflow predictions in ungauged basins: Parameter sensitivity analysis

A procedure combining the Soil Conservation Service‐Curve Number (SCS‐CN) method and the Green–Ampt (GA) infiltration equation was recently developed to overcome some of the drawbacks of the classic SCS‐CN approach when estimating the volume of surface runoff at a sub‐daily time resolution. The rationale of this mixed procedure, named Curve Number for Green–Ampt (CN4GA), is to use the GA infiltration model to distribute the total volume of the net hyetograph (rainfall excess) provided by the SCS‐CN method over time. The initial abstraction and the total volume of rainfall given by the SCS‐CN method are used to identify the ponding time and to quantify the hydraulic conductivity parameter of the GA equation. In this paper, a sensitivity analysis of the mixed CN4GA parameters is presented with the aim to identify conditions where the mixed procedure can be effectively used within the Prediction in Ungauged Basin perspective. The effects exerted by changes in selected input parameters on the outputs are evaluated using rectangular and triangular synthetic hyetographs as well as 100 maximum annual storms selected from synthetic rainfall time series. When applied to extreme precipitation events, which are characterized by predominant peaks of rainfall, the CN4GA appears to be rather insensitive to the input hydraulic parameters of the soil, which is an interesting feature of the CN4GA approach and makes it an ideal candidate for the rainfall excess estimation at sub‐daily temporal resolution at ungauged sites. Copyright © 2013 John Wiley & Sons, Ltd.     

Two‐dimensional continuous simulation of spatiotemporally varied hydrological processes using the time‐area method

Distributed, continuous hydrologic models promote better understanding of hydrology and enable integrated hydrologic analyses by providing a more detailed picture of water transport processes across the varying landscape. However, such models are not widely used in routine modelling practices, due in part to the extensive data input requirements, computational demands, and complexity of routing algorithms. We developed a two‐dimensional continuous hydrologic model, HYSTAR, using a time‐area method within a grid‐based spatial data model with the goal of providing an alternative way to simulate spatiotemporally varied watershed‐scale hydrologic processes. The model calculates the direct runoff hydrograph by coupling a time‐area routing scheme with a dynamic rainfall excess sub‐model implemented here using a modified curve number method with an hourly time step, explicitly considering downstream ‘reinfiltration’ of routed surface runoff. Soil moisture content is determined at each time interval based on a water balance equation, and overland and channel runoff is routed on time‐area maps, representing spatial variation in hydraulic characteristics for each time interval in a storm event. Simulating runoff hydrographs does not depend on unit hydrograph theory or on solution of the Saint Venant equation, yet retains the simplicity of a unit hydrograph approach and the capability of explicitly simulating two‐dimensional flow routing. The model provided acceptable performance in predicting daily and monthly runoff for a 6‐year period for a watershed in Virginia (USA) using readily available geographic information about the watershed landscape. Spatial and temporal variability in simulated effective runoff depth and time area maps dynamically show the areas of the watershed contributing to the direct runoff hydrograph at the outlet over time, consistent with the variable source area overland flow generation mechanism. The model offers a way to simulate watershed processes and runoff hydrographs using the time‐area method, providing a simple, efficient, and sound framework that explicitly represents mechanisms of spatially and temporally varied hydrologic processes. Copyright © 2015 John Wiley & Sons, Ltd.     

Rainfall–runoff simulation using a normalized antecedent precipitation index

Abstract

The normalized antecedent precipitation index (NAPI) model by Heggen for the prediction of runoff yield is analytically derived from the water balance equation. Heggen's model has been simplified further to a rational form and its performance verified with the Soil Conservation Service Curve Number (SCS-CN) model. The simplified model has three coefficients specific to a watershed, and requires two inputs: rainfall and the derived parameter, NAPI. The characteristic behaviour of the NAPI has resonance with the curve number (CN) of the SCS model. The proposed NAPI model was applied to three watersheds in the semi-arid region of India to simulate runoff yield. The model showed improved correlation between the observed and predicted runoff data compared to the SCS-CN model. The F test and paired t test also confirmed the reliability of the model with significance levels of 0.01 and 0.001%, respectively. The proposed model could be used successfully for rainfall–runoff modelling in a watershed.

Citation Ali, S., Ghosh, N. C. & Singh, R. (2010) Rainfall–runoff simulation using a normalized antecedent precipitation index. Hydrol. Sci. J. 55(2), 266–274.     

Effect of the inter-annual variability of rainfall statistics on stochastically generated rainfall time series: part 2. Impact on watershed response variables

This study analyzes how the stochastically generated rainfall time series accounting for the inter-annual variability of rainfall statistics can improve the prediction of watershed response variables such as peak flow and runoff depth. The modified Bartlett–Lewis rectangular pulse (MBLRP) rainfall generation model was improved such that it can account for the inter-annual variability of the observed rainfall statistics. Then, the synthetic rainfall time series was generated using the MBLRP model, which was used as input rainfall data for SCS hydrologic models to produce runoff depth and peak flow in a virtual watershed. These values were compared to the ones derived from the synthetic rainfall time series that is generated from the traditional MBLRP rainfall modeling. The result of the comparison indicates that the rainfall time series reflecting the inter-annual variability of rainfall statistics reduces the biasness residing in the predicted peak flow values derived from the synthetic rainfall time series generated using the traditional MBLRP approach by 26–47 %. In addition, it was observed that the overall variability of the peak flow and run off depth distribution was better represented when the inter-annual variability of rainfall statistics are considered.     

FLOOD MODELLING WITH GIS-DERIVED DISTRIBUTED UNIT HYDROGRAPHS

The concept of a spatially distributed unit hydrograph is based on the fact that the unit hydrograph can be derived from the time–area curve of a watershed by the S-curve method. The time–area diagram is a graph of cumulative drainage area contributing to discharge at the watershed outlet within a specified time of travel. Accurate determination of the time–area diagram is made possible by using a GIS. The GIS is used to describe the connectivity of the links in the watershed flow network and to calculate distances and travel times to the watershed outlet for various points within the watershed. Overland flow travel times are calculated by the kinematic wave equation for time to equilibrium; channel flow times are based on the Manning and continuity equations. To account for channel storage, travel times for channel reaches are increased by a percentage depending on the channel reach length and geometry. With GIS capability for rainfall mapping, the assumption of a uniform spatial rainfall distribution is no longer necessary; hence the term, spatially distributed unit hydrograph. An example of the application for the Waiparous Creek in the Alberta Foothills is given. IDRISI is used to develop a simple digital elevation model of the 229 km² watershed, using 1 km × 1 km grid cells. A grid of flow directions is developed and used to create an equivalent channel network. Excess rainfall for each 1 km × 1 km cell is individually computed by the Soil Conservation Service (SCS) runoff curve method and routed through the equivalent channel network to obtain the time–area curve. The derived unit hydrograph gave excellent results in simulating an observed flood hydrograph. The distributed unit hydrograph is no longer a lumped model, since it accounts for internal distribution of rainfall and runoff. It is derived for a watershed without the need for observed rainfall and discharge data, because it is essentially a geomorphoclimatic approach. As such, it allows the derivation of watershed responses (hydrographs) to inputs of various magnitudes, thus eliminating the assumption of proportionality of input and output if needed. The superposition of outputs is retained in simulating flood hydrographs by convolution, since it has been shown that some non-linear systems satisfy the principle of superposition. The distributed unit hydrograph appears to be a very promising rainfall runoff model based on GIS technology.     

Reliability assessment of water structures subject to data scarcity using the SCS-CN model

ABSTRACT

When discharge measurements are not available, design of water structures relies on using frequency analysis of rainfall data and applying a rainfall–runoff model to estimate a hydrograph. The Soil Conservation Service (SCS) method estimates the design hydrograph first through a rainfall–runoff transformation and next by propagating runoff to the basin outlet via the SCS unit hydrograph (UH) method. The method uses two parameters, the Curve Number (CN) and the time of concentration (T_c). However, in data-scarce areas, the calibration of CN and T_c from nearby gauged watersheds is limited and subject to high uncertainties. Therefore, the inherent uncertainty/variability of the SCS parameters may have considerable ramifications on the safety of design. In this research, a reliability approach is used to evaluate the impact of incorporating the uncertainty of CN and T_c in flood design. The sensitivity of the probabilistic outcome against the uncertainty of input parameters is calculated using the First Order Reliability Method (FORM). The results of FORM are compared with the conventional SCS results, taking solely the uncertainty of the rainfall event. The relative importance of the uncertainty of the SCS parameters is also estimated. It is found that the conventional approach, used by many practitioners, might grossly underestimate the risk of failure of water structures, due to neglecting the probabilistic nature of the SCS parameters and especially the Curve Number. The most predominant factors against which the SCS-CN method is highly uncertain are when the average rainfall value is low (less than 20 mm) or its coefficient of variation is not significant (less than 0.5), i.e. when the resulting rainfall at the design return period is low. A case study is presented for Egypt using rainfall data and CN values driven from satellite information, to determine the regions of acceptance of the SCS-CN method.

EDITOR D. Koutsoyiannis; ASSOCIATE EDITOR A. Efstratiadis     

SWAT application to estimate design runoff curve number for South Korean conditions

Runoff estimations based on the standard USDA–NRCS curve number (CN) table without calibration have a tendency to give inaccurate results when the CN values are applied in South Korea which has many high slope watersheds and that has a continental monsoon climate. Particularly for the design flood estimation, accurately calibrated CN values are required because the estimated peak flow is very sensitive to the selection of CN. However, the lack of flood data makes it difficult to calibrate and assign runoff CNs to Korean watersheds. Even if sufficient data are available to estimate CN values, it is also difficult to obtain the direct flows by separating base flows from total runoff hydrographs due to the temporal irregularity of rainfall events and the resulting complex pattern of runoff. Therefore, an alternative method for estimating CNs needs to be developed to overcome these issues. The purpose of this study is to present a method for estimating runoff CNs using the soil and water assessment tool (SWAT) model which can take into account watershed heterogeneities such as climate conditions, land use and soil types. The proposed CN estimation method uses the simulated flow data by SWAT instead of using measured flow data. This method has advantages in estimating CN values spatially for each subbasin division considering watershed characteristics. The use of daily data can reduce the sensitivity to the abnormality that is commonly involved in flow data with a small time scale. The SWAT‐based CN estimation method, combined with the asymptotic CN method, was applied to the Chungju dam watershed in South Korea. A regression equation was then developed from this approach, which was used to estimate CN values that decrease exponentially as rainfall amounts increase and that converge to 60·6 and 79·4 without and with considering subsurface lateral flow, respectively. Furthermore, the CN values for the antecedent moisture conditions were determined using the probabilistic approach. The CN associated with the 50% probability for the Chungju dam watershed is 87·8 which can be taken to be representative of antecedent moisture condition (AMC) II. The CN_I and CN_III associated with 90% and 10% probabilities are 78·9 and 94·1, respectively. The estimated CN_II = 87·8 differs markedly from the geographic information system (GIS)‐based CN 65·0, which implies that the standard USDA–NRCS CN method should be calibrated to the studied area of interest. Copyright © 2010 John Wiley & Sons, Ltd.     

Dynamics of stormflow generation—A hillslope‐scale field study in east‐central Pennsylvania,USA

A 40 m × 20 m mowed, grass hillslope adjacent to a headwater stream within a 26‐ha watershed in east‐central Pennsylvania, USA, was instrumented to identify and map the extent and dynamics of surface saturation (areas with the water table at the surface) and surface runoff source areas. Rainfall, stream flow and surface runoff from the hillslope were recorded at 5‐min intervals from 11 August to 22 November 1998, and 13 April to 12 November 1999. The dynamics of the water table (0 to 45 cm depth from the soil surface) and the occurrence of surface runoff source areas across the hillslope were recorded using specially designed subsurface saturation and surface runoff sensors, respectively. Detailed data analyses for two rainfall events that occurred in August (57·7 mm in 150 min) and September (83·6 mm in 1265 min) 1999, illustrated the spatial and temporal dynamics of surface saturation and surface runoff source areas. Temporal data analyses showed the necessity to measure the hillslope dynamics at time intervals comparable to that of rainfall measurements. Both infiltration excess surface runoff (runoff caused when rainfall intensity exceeds soil infiltration capacity) and saturation excess surface runoff (runoff caused when soil moisture storage capacity is exceeded) source areas were recorded during these rainfall events. The August rainfall event was primarily an infiltration excess surface runoff event, whereas the September rainfall event produced both infiltration excess and saturation excess surface runoff. Occurrence and disappearance of infiltration excess surface runoff source areas during the rainfall events appeared scattered across the hillslope. Analysis of surface saturation and surface runoff data showed that not all surface saturation areas produced surface runoff that reached the stream. Emergence of subsurface flow to the surface during the post‐rainfall periods appeared to be a major flow process dominating the hillslope after the August rainfall event. Copyright © 2002 John Wiley & Sons, Ltd.     

Estimation of SCS Curve Number variation following forest fires

The Soil Conservation Service Curve Number (CN) method is routinely used to estimate the effects of forest fires on hydrological response. However, despite recent efforts, CN values are still not well known for burned conditions. A major forest fire in Attica, Greece, which affected the Lykorrema stream experimental watershed, provided an opportunity for the estimation of post-fire CN variation using detailed pre-fire and post-fire rainfall–runoff datasets. The CN values for both periods were estimated and compared using a wide range of available methods. Methods considering the spatial variability of soil-cover complexes were also used to investigate the effect of spatial heterogeneity. The post-fire watershed response changed from complacent to standard. Direct runoff depths and peak flows increased by a factor of more than 7.7 and 11.8, respectively. On average, the estimated post-fire CN values for the studied soil-cover complexes increased by about 25 units. This study may assist the improvement of existing post-fire hydrological assessment tools.     

A relook at NEH‐4 curve number data and antecedent moisture condition criteria

This paper investigates the variation of the popular curve number (CN) values given in the National Engineering Hand Book–Section 4 (NEH‐4) of the Soil Conservation Service (SCS) with antecedent moisture condition (AMC) and soil type. Using the volumetric concept, involving soil, water, and air, a significant condensation of the NEH‐4 tables is achieved. This leads to a procedure for determination of CN for gauged as well as ungauged watersheds. The rainfall‐runoff events derived from daily data of four Indian watersheds exhibited a power relation between the potential maximum retention or CN and the 5‐day antecedent rainfall amount. Including this power relation, the SCS‐CN method was modified. This modification also eliminates the problem of sudden jumps from one AMC level to the other. The runoff values predicted using the modified method and the existing method utilizing the NEH‐4 AMC criteria yielded similar results. Copyright © 2006 John Wiley & Sons, Ltd.     

A spatially distributed Clark’s unit hydrograph based hybrid hydrologic model (Distributed-Clark)

ABSTRACT

A hybrid hydrologic model (Distributed-Clark), which is a lumped conceptual and distributed feature model, was developed based on the combined concept of Clark’s unit hydrograph and its spatial decomposition methods, incorporating refined spatially variable flow dynamics to implement hydrological simulation for spatially distributed rainfall–runoff flow. In Distributed-Clark, the Soil Conservation Service (SCS) curve number method is utilized to estimate spatially distributed runoff depth and a set of separated unit hydrographs is used for runoff routing to obtain a direct runoff flow hydrograph. Case studies (four watersheds in the central part of the USA) using spatially distributed (Thiessen polygon-based) rainfall data of storm events were used to evaluate the model performance. Results demonstrate relatively good fit to observed streamflow, with a Nash-Sutcliffe efficiency (E_NS) of 0.84 and coefficient of determination (R²) of 0.86, as well as a better fit in comparison with outputs of spatially averaged rainfall data simulations for two models including HEC-HMS.     

Estimating effective impervious area in urban watersheds using land cover,soil character and asymptotic curve number

Knowledge of the effective impervious area (EIA) or the degree to which impervious surfaces are hydraulically connected to the drainage system is useful for improving hydrological and environmental models and assessing the effectiveness of green stormwater infrastructure in urban watersheds. The goal of this research is to develop a method to estimate EIA fraction in urban watersheds using readily available data. Since EIA is dependent on rainfall–runoff response and cannot be solely determined based on the physical characteristics of a watershed, the EIA is linked with the asymptotic curve number (CN), a watershed index that represents runoff characteristics. In order for the method to be applicable to ungauged watersheds, the asymptotic CN is predicted using land cover and soil data from 35 urban catchments in Minnesota and Texas, USA. Similar data from 11 other urban catchments in Wisconsin and Texas, USA, are used to validate the results. A set of runoff depth versus EIA fraction curves is also developed to assess the impact of EIA reduction on discharge from an urban watershed in land-use planning studies.     

Sensitivity and fuzzy uncertainty analyses in the determination of SCS-CN parameters from rainfall–runoff data

Spatial and seasonal variations of curve number (CN) and initial abstraction ratio (λ) in a watershed can result in inaccurate runoff volume estimations when using the US Natural Resources Conservation Service (SCS-CN) method with constant values for these parameters. In this paper, parameters of CN and λ are considered as calibration parameters and the sensitivity of estimated runoff to these parameters using the SCS-CN method is scrutinized. To incorporate the uncertainty associated with CN and λ, fuzzy linear regression (FLR) is applied to derive the relationships of CN and λ with rainfall depth (P) by employing a large dataset of storm events from four watersheds in Iran. Results indicate that the proposed approach provides more accuracy in estimation of runoff volume compared to the SCS method with constant values of CN and λ, and gives a straightforward technique for evaluating the hydrological effects of CN, λ, and P on runoff volume.     

A modification to the Soil Conservation Service curve number method for steep slopes in the Loess Plateau of China

The Soil Conservation Service curve number (CN) method is widely used for predicting direct runoff from rainfall. However, despite the extent of cultivation on hillslope areas, very few attempts have been made to incorporate a slope factor into the CN method. The objectives of this study were (1) to evaluate existing approaches integrating slope in the CN method, and (2) to develop an equation incorporating a slope factor into the CN method for application in the steep slope areas of the Loess Plateau of China. The dataset consisted of 11 years of rainfall and runoff measurements from two experimental sites with slopes ranging from 14 to 140%. The results indicated that the standard CN method underestimated large runoff events and overestimated small events. For our experimental conditions, the optimized and non‐optimized forms of the slope‐modified CN method of the Erosion Productivity Impact Calculator model improved runoff prediction for steep slopes, but large runoff events were still underestimated and small ones overpredicted. Based on relationships between slope and the observed and theoretical CN values, an equation was developed that better predicted runoff depths with an R² of 0·822 and a linear regression slope of 0·807. This slope‐adjusted CN equation appears to be the most appropriate for runoff prediction in the steep areas of the Loess Plateau of China. Copyright © 2005 John Wiley & Sons, Ltd.