Analysis of runoff in ungauged mountain watersheds in Sichuan, China using kinematic-wave-based GIUH model

Floods are one of the most common natural hazards occurring all around the world. However, the knowledge of the origins of a food and its possible magnitude in a given region remains unclear yet. This lack of understanding is particularly acute in mountainous regions with large degrees in Sichuan Province, China, where runoff is seldom measured. The nature of streamflow in a region is related to the time and spatial distribution of rainfall quantity and watershed geomorphology. The geomorphologic characteristics are the channel network and surrounding landscape which transform the rainfall input into an output hydrograph at the outlet of the watershed. With the given geomorphologic properties of the watershed, theoretically the hydrological response function can be determined hydraulically without using any recorded data of past rainfall or runoff events. In this study, a kinematic-wave-based geomorphologic instantaneous unit hydrograph (KW-GIUH) model was adopted and verified to estimate runoff in ungauged areas. Two mountain watersheds, the Yingjing River watershed and Tianquan River watershed in Sichuan were selected as study sites. The geomorphologic factors of the two watersheds were obtained by using a digital elevation model (DEM) based on the topographic database obtained from the Shuttle Radar Topography Mission of US’s NASA. The tests of the model on the two watersheds were performed both at gauged and ungauged sites. Comparison between the simulated and observed hydrographs for a number of rainstorms at the gauged sites indicated the potential of the KW-GIUH model as a useful tool for runoff analysis in these regions. Moreover, to simulate possible concentrated rainstorms that could result in serious flooding in these areas, synthetic rainfall hyetographs were adopted as input to the KW-GIUH model to obtain the flow hydrographs at two ungauged sites for different return period conditions. Hydroeconomic analysis can be performed in the future to select the optimum design return period for determining the flood control work.     

Application of SCS model in estimation of runoff from small watershed in Loess Plateau of China

Soil Conservation Service （SCS） model, developed by U. S. Soil Conservation Service in 1972, has been widely applied in the estimation of runoff from an small watershed. In this paper, based on the remote sensing geo-information data of land use and soil classification all obtained from Landsat images in 1996 and 1997 and conventional data of hydrology and meteorology, the SCS model was investigated for simulating the surface runoff for single rainstorm in Wangdonggou watershed, a typical small watershed in the Loess Plateau, located in Changwu County of Shaanxi Province of China. Wangdonggou watershed was compartmentalized into 28 sub-units according to natural draining division, and the table of curve number （CN） values fitting for Wangdonggou watershed was also presented. During the flood period from 1996 to 1997, the hydrograph of calculated runoff process using the SCS model and the hydrograph of observed runoff process coincided very well in height as well as shape, and the model was of high precision above 75%. It is indicated that the SCS model is legitimate and can be successfully used to simulate the runoff generation and the runoff process of typical small watershed based on the remote sensing geo-information in the Loess Plateau.     

Development of a Cell-based model to derive direct runoff hydrographs for ungauged mountainous basins

A model to derive direct runoff hydrograph for an ungauged basin using the physical properties of the basin is presented. The basin is divided into grid cells and canal elements. Overland flow is generated from each grid cell of the basin by application of continuous effective rainfall of 1 mm/hr to the basin. The flow generated is routed through downstream grid cells and the canal elements using the kinematic wave approach. The travel time for direct runoff from each grid cell to the basin outlet is calculated and the S-curve is derived for the basin. The S-curve is used to derive the unit hydrograph of a given duration for the basin. The model, referred as Cell-basin model was applied to the Upper Kotmale Basin in Sri Lanka and the model predictions of direct runoff hydrographs for rainfall events agreed with the observations to a reasonable accuracy. Comparison of the unit hydrographs obtained from the model and from the conventional Snyder’s synthetic unit hydrograph using regionalized parameters assuming the basin as an ungauged basin, with the unit hydrograph derived from the observations showed that the model predicted unit hydrograph was more suitable than that obtained by Snyder’s method for Sri Lankan up country basins. Thus, the present model is a useful tool to obtain direct runoff hydrograph for ungauged basins.     

Suspended sediment transport analysis in two Italian instrumented catchments

Suspended sediment transport in streams is an effective indicator of soil erosion at the watershed scale. This process was studied using the data obtained from two continuous monitoring stations installed in Italian watersheds, the Rio Cordon and Torrente Carapelle. The catchments have substantially different climatic, morphological and land use characteristics. The Rio Cordon, a small Alpine watershed, has been monitored since 1986, while in the Torrente Carapelle, a medium-size Mediterranean watershed, the monitoring station has provided data since 2007. Several flood events with suspended sediment transport were isolated in the two catchments, excluding those determined by snowmelt in the Rio Cordon as this process does not affect the Carapelle watershed. Analysis of the events showed different behavior in terms of hysteresis loop trends between water discharge (m³ s^-1) and suspended sediment concentration (g l^-1) values, as the Rio Cordon confirms clockwise relationships most frequently, while counter-clockwise trends represent the majority of cases for the Carapelle Torrent. The different behavior of the two watersheds was further analyzed by evaluating relationships between the main hydrological parameters related to suspended sediment transport. Runoff controls the total sediment budget in both catchments. In contrast, it was noted that the runoff process does not interfere in the magnitude of the suspended sediment transport in the Rio Cordon catchment, while variations due to the larger size of the catchment area characterize the events in the Carapelle watershed. Lastly, a flow peak threshold that determines the advance or delay of the hydrograph peak with respect to the sedigraph peak was established for the Carapelle Torrent, while it was not reported in the Rio Cordon.     

Hydrological Response to Environment Change in Himalayan Watersheds: Assessment from Integrated Modeling Approach

Land use changes such as deforestation,increase in cropping or grazing areas and built-up land, likely modify the water balance and land surface behavior in the Himalayan watersheds.An integrated approach of hydrological and hydraulic modeling was adopted for comparative analysis of hydrological pattern in three Himalayan watersheds i.e.Khanpur,Rawal and Simly situated in the Northern territory of Pakistan.The rainfall-runoff model SWAT- Soil and water assessment tool and Hydro CAD were calibrated for the selected watersheds.The correlation analysis of the precipitation data of two climate stations i.e.Murree and Islamabad, with the discharge data of three rivers was utilized to select best suitable input precipitation data for Hydro CAD rainfall-runoff modeling.The peak flood hydrograph were generated using Hydro CAD runoff to optimize the basin parameters like CN, runoff volume, peak flows of the three watersheds.The hydrological response of the Rawal watershed was studied as a case study to different scenarios of land use change using SWAT model.The scenario of high deforestation indicated a decline of about 6.3% in the groundwater recharge tostream while increase of 7.1% in the surface runoff has been observed under the scenario of growth in urbanization in the recent decades.The integrated modeling approach proved helpful in investigating the hydrological behavior under changing environment at watershed level in the Himalayan region.     

Snowmelt modeling using two melt-rate models in the Urumqi River watershed,Xinjiang Uyghur Autonomous Region,China

In this paper,the performance of the classic snowmelt runoff model(SRM)is evaluated in a daily discharge simulation with two different melt models,the empirical temperature-index melt model and the energy-based radiation melt model,through a case study from the data-sparse mountainous watershed of the Urumqi River basin in Xinjiang Uyghur Autonomous Region of China.The classic SRM,which uses the empirical temperature-index method,and a radiation-based SRM,incorporating shortwave solar radiation and snow albedo,were developed to simulate daily runoff for the spring and summer snowmelt seasons from 2005 to 2012,respectively.Daily meteorological and hydrological data were collected from three stations located in the watershed.Snow cover area(SCA)was extracted from satellite images.Solar radiation inputs were estimated based on a digital elevation model(DEM).The results showed that the overall accuracy of the classic SRM and radiation-based SRM for simulating snowmeltdischarge was relatively high.The classic SRM outperformed the radiation-based SRM due to the robust performance of the temperature-index model in the watershed snowmelt computation.No significant improvement was achieved by employing solar radiation and snow albedo in the snowmelt runoff simulation due to the inclusion of solar radiation as a temperature-dependent energy source and the local pattern of snowmelt behavior throughout the melting season.Our results suggest that the classic SRM simulates daily runoff with favorable accuracy and that the performance of the radiation-based SRM needs to be further improved by more ground-measured data for snowmelt energy input.     

流域空间离散化及其对径流过程模拟的影响研究

本文分析比较了流域空间离散化的两种途径 :子流域法和网格法 ,并以此对流域径流过程模拟的影响进行了探讨。结果表明空间尺度和时间尺度之间应该是相互关联、相互匹配的 ;时空尺度具体大小的选定要根据研究目的、研究区域的面积、流域下扩建面信息源和模型输入资料的可获性综合考虑。     

Using uncertainty and sensitivity analysis for finding the best rainfall-runoff model in mountainous watersheds(Case study: the Navrood watershed in Iran)

In watersheds that have not sufficient meteorological and hydrometric data for simulating rainfall-runoff events, using geomorphologic and geomorphoclimatic characteristics of watershed is a conventional method for the simulation. A number of rainfall-runoff models utilize these characteristics such as Nash-IUH, Clark-IUH, Geomorphologic Instantaneous Unit Hydrograph(GIUH), Geomorphoclimatic Instantaneous Unit Hydrograph(GcIUH), GIUH-based Nash(GIUH-Nash) and GcIUH-based Clark(GcIUH-Clark). But all these models are not appropriate for mountainous watersheds. Therefore, the objective of this study is to select the best of them for the simulation. The procedure of this study is: a) selecting appropriate rainfall-runoff events for calibration and validation of six hybrid models, b) distinguishing the best model based on different performance criteria(Percentage Error in Volume(PEV); Percentage Error in Peak(PEP); Percentage Error in Time to Peak(PETP); Root Mean Square Error(RMSE) and Nash-Sutcliffe model efficiency coefficient(ENS)), c) Sensitivity analysis for determination of the most effective parameter at each model, d) Uncertainty determination of different parameters in each model and confirmation of the obtained results by application of the performance criteria. For application of this procedure, the Navrood watershed in the north of Iran as a mountainous watershed has been considered. The results showed that the ClarkIUH and GcIUH-Clark are suitable models for simulation of flood hydrographs, while other models cannot simulate flood hydrographs appropriately. The sensitivity analysis shows that the most sensitive parameters are the infiltration constant rate and time of concentration in the Clark-IUH model. Also, the most sensitive parameters include the infiltration constant rate and storage coefficient in the GcIUHClark model. The Clark-IUH and GcIUH-Clark models are more sensitive to their parameters. The Latin Hypercube Sampling(LHS) on Monte Carlo(MC) simulation method was used for evaluation of uncertainty of data in rainfall-runoff models. In this method 500 sets of data values are produced and then the peak discharge of flood hydrographs for each produced data set is simulated with rainfall-runoff models. The uncertainty of data changes the value of simulated peak discharge of flood hydrograph. The uncertainty analysis shows that the observed peak discharges of different rainfall-runoff events are within the range of values of simulated by the six hybrid rainfall-runoff models and IUH that inputs of these models were the produced data sets. The range of the produced peak discharge of flood hydrographs by the Clark-IUH and GcIUH-Clark models is wider than those of other models.     

黄土地貌区径流节点提取及其特征信息指数的DEM实验分析

径流节点是一种真实存在于流域单元中的地学实体对象,携带了大量的水文与地貌特征信息。针对当前研究中缺乏定量地度量径流节点所承载的地学特征信息量的相关探讨,本文结合径流节点的已有理论,总结了其地学特征信息,借鉴数学粗集理论提出了径流节点的特征信息指数（FIIRN）,并分析了相关的定义与计算方法。本文以黄土高原的7个典型黄土地貌样区5 m分辨率DEM数据为实验材料,提取了各个实验样区的径流节点,并计算了每个样区的FIIRN指数值。实验结果表明：（1）FIIRN指数可有效地映射流域单元的径流分级、发育状况等地学特征信息;（2）FIIRN指数可定量地表达流域单元的地学特征信息量;（3）FIIRN指数可作为挖掘流域地学特征信息的一个重要参数。     

Influence of vegetation on runoff and sediment in wind-water erosion crisscross region in the upper Yellow River of China

All characteristics of vegetation,runoff and sediment from 1960 to 2010 in the Xiliu Gully Watershed,which is a representative watershed in wind-water erosion crisscross region in the upper reaches of the Yellow River of China,have been analyzed in this study.Based on the remote sensing image data,and used multi-spectral interpretation method,the characteristics of vegetation variation in the Xiliu Gully Watershed have been analyzed.And the rules of precipitation,runoff and sediment's changes have been illuminated by using mathematical statistics method.What′s more,the influence mechanism of vegetation on runoff and sediment has been discussed by using the data obtained from artificial rainfall simulation test.The results showed that the main vegetation type was given priority to low coverage,and the area of the low vegetation coverage type was reducing year by year.On the country,the area of the high vegetation coverage type was gradually increasing.In a word,vegetation conditions had got better improved since 2000 when the watershed management project started.The average annual precipitation of the river basin also got slightly increase in 2000–2010.The average annual runoff reduced by 37.5%,and the average annual sediment reduced by 73.9% in the same period.The results of artificial rainfall simulation tests showed that the improvement of vegetation coverage could increase not only soil infiltration but also vegetation evapotranspiration,and then made the rainfall-induced runoff production decrease.Vegetation root system could increases the resistance ability of soil to erosion,and vegetation aboveground part could reduce raindrop kinetic energy and splash soil erosion.Therefore,with the increase of vegetation coverage,the rainfall-induced sediment could decrease.     

Influence of Land Use/Cover Change on Storm Runoff— A Case Study of Xitiaoxi River Basin in Upstream of Taihu Lake Watershed

Land use/cover change (LUCC) is one of the main boundary conditions which influence many hydrologic processes. In view of the importance of Taihu Lake Watershed in China and the urgency of discovering the impacts of LUCC on storm runoff, two flood events under five land cover scenarios in the Xitiaoxi River Basin of the upstream of Taihu Lake watershed were simulated by distributed hydrologic modeling system HEC-HMS. The influences of each land cover on storm runoff were discussed. It was concluded that under the same rainstorm the ascending order of runoff coefficient and peak flow produced by the 5 different land covers were woodland, shrub, grassland, arable land, and built-up land; the descending order of swelling time were woodland, shrub, grassland, arable land, and built-up land. Scenario of built-up land was the first to reach peak flow, then arable land, grassland, shrub, and woodland. There were close relationships between the runoff coefficients produced by the 5 different land covers. The degrees of impacts on runoff coefficient of land cover change modes were sorted by descending: woodland to built-up land, shrub to built-up land, grassland to built-up land, arable land to built-up land, woodland to arable land, shrub to arable land, arable land to grassland, shrub to grassland, grassland to arable land, and woodland to shrub. Urbanization will contribute to flood disaster, while forestation will mitigate flood disaster.     

黄河中游黄土高原区河川基流特点及变化分析

由于气候变化和人类活动的影响。黄土高原河川基流呈显著下降趋势。对黄河径流量减少产生一定影响。采用日流量过程线切割法计算了河龙区间黄土高原主要支流的基流量。结合地形地貌、水文地质条件，分析了不同地貌类型支流基流的特点，初步分析了影响基流变化的各种因素。认为河龙区间基流变化对黄河水资源影响较大。     

Non-point-source nitrogen and phosphorus loadings from a small watershed in the Three Gorges Reservoir area

Non-point-source pollution has become a major threat to the water quality of the Three Gorges Reservoir(TGR);however,nutrient loadings from terrestrial sources are unclear due to a lack of in situ monitoring.A representative small watershed in the central part of the TGR area was selected to monitor the loss of nitrogen(N) and phosphorous(P) continuously along with the runoff from 2007 through 2009 to understand the exact sources and loadings.Results show that the non-point-source nitrogen and phosphorus comes mainly from the storm runoff from residential areas,citrus orchards and sloping croplands,which contributes up to 76% of the loadings in this watershed.Thus,a crucial measure for controlling non-point-source pollution is to intercept storm runoff from the three sources.Paddy fields provide a sink for non-point-source N and P by intercepting the runoff and sediment along with their different forms of nitrogen and phosphorus.The N and P removal efficiency by paddy fields from residential areas is within the range of 56% to 98%.Paddy fields are an important land-use option for reducing the non-point-source loading of N and P in the TGR area.     

Searching for an Optimized Single-objective Function Matching Multiple Objectives with Automatic Calibration of Hydrological Models

In the calibration of hydrological models, evaluation criteria are explicitly and quantitatively defined as single-or multi-objective functions when utilizing automatic calibration approaches.In most previous studies, there is a general opinion that no single-objective function can represent all important characteristics of even one specific hydrological variable(e.g., streamflow).Thus hydrologists must turn to multi-objective calibration.In this study, we demonstrated that an optimized single-objective function can compromise multi-response modes(i.e., multi-objective functions) of the hydrograph, which is defined as summation of a power function of the absolute error between observed and simulated streamflow with the exponent of power function optimized for specific watersheds.The new objective function was applied to 196 model parameter estimation experiment(MOPEX) watersheds across the eastern United States using the semi-distributed Xinanjiang hydrological model.The optimized exponent value for each watershed was obtained by targeting four popular objective functions focusing on peak flows, low flows, water balance, and flashiness, respectively.Results showed that the optimized single-objective function can achieve a better hydrograph simulation compared to the traditional single-objective function Nash-Sutcliffe efficiency coefficient for most watersheds, and balance high flow part and low flow part of the hydrograph without substantial differences compared to multi-objective calibration.The proposed optimal single-objective function can be practically adopted in the hydrological modeling if the optimal exponent value could be determined a priori according to hydrological/climatic/landscape characteristics in a specific watershed.     

Impacts of precipitation variation and soil and water conservation measures on runoff and sediment yield in the Loess Plateau Gully Region,China

The Loess Plateau of China has experienced a lengthy drought and severe soil erosion.Changes in precipitation and land use largely determine the dynamics of runoff and sediment yield in this region. Trend and mutation analyses were performed on hydrological data(1981–2012) from the Yanwachuan watershed in the Loess Plateau Gully Region to study the evolution characteristics of runoff and sediment yield. A time-series contrasting method also was used to evaluate the effects of precipitation and soil and water conservation(SWC) on runoff and sediment yield. Annual sediment yield declined markedly from 1981 to 2012 although there was no significant change in annual precipitation and annual runoff. Change points of annual runoff and annual sediment yield occurred in 1996 and 1997,respectively. Compared with that in the baseline period(1981–1996), annual runoff and annual sediment yield in the change period(1997–2012)decreased by 17.0% and 76.0%, respectively, but annual precipitation increased by 6.3%. Runoff decreased in the flood season and normal season, but increased in the dry season, while sediment yield significantly declined in the whole study period. The SWC measures contributed significantly to the reduction of annual runoff(137.9%) and annual sediment yield(135%) and were more important than precipitation. Biological measures(forestland and grassland) accounted for 61.04% of total runoff reduction, while engineering measures(terraces and dams) accounted for 102.84% of total sediment yield reduction. Furthermore, SWC measures had positive ecological effects. This study provides a scientific basis for soil erosion control on the Loess Plateau.     

基于DEM数字河网提取时集水面积阈值与河源密度关系的研究

以研究区3030个流域水文数据为基础,利用数字高程模型（DEM）并运用非线性拟合法分析中小流域集水面积阈值与河源密度的相关关系,分析集水面积阈值与河源密度的幂函数关系可得,幂指数a近似为-1,拟合公式常系数k具有区域分布特征,反映了不同地区的河网发育程度,k值越大,河源密度越大,河网越密集,河网发育程度越高。对河源密度随集水面积阈值的变化趋势进行分析,发现2个具有统计意义的阈值临界点,分别为河源从坡面到沟道及从沟道到河道的集水面积分界点,在研究区内选择不同k值区间的贡曲流域、辰清河流域和藤条江流域作为样本流域,利用样本流域数字正射影像（DOM）手工提取其沟谷河网及主河道河网的图形信息和水文信息,并与分界点处阈值提取的河网作对比,将分析得出的河网密度相对误差率作为检验标准,对临界点的物理意义进行检验。     

广西龙江流域水环境状况及地下水资源潜力评价

龙江流域是西南岩溶区一个典型的水文盆地。由于降水丰沛 ,地下水、地表水资源十分丰富 ,但受工矿企业废水的污染 ,龙江河水在枯水季节已基本失去了作为人畜饮用水水源的功能 ,沿江许多城镇人畜饮用水困难。在水资源十分丰富的西南岩溶区 ,类似的许多沿江城市出现水质性缺水 ,给水文地质工作提出了新的要求。文章认为以一定规模的水文盆地为基本单元 ,以流域地表水地下水资源的可利用性的综合评价为出发点 ,有针对性开展缺水城镇及岩溶干旱片的供水论证及地下水生态环境调查 ,是新一轮水文地质工作的重点 ,也是基础性公益性水工环工作的首要任务。     

Effects of slope gradient on runoff from bare-fallow purple soil in China under natural rainfall conditions

Purple soil is highly susceptible for overland flow and surface erosion, therefore understanding surface runoff and soil erosion processes in the purple soil region are important to mitigate flooding and erosion hazards. Slope angle is an important parameter that affects the magnitude of runoff and thus surface erosion in hilly landscapes or bare land area. However, the effect of slope on runoff generation remains unclear in many different soils including Chinese purple soil. The aim of this study was to investigate the relationship between different slope gradients and surface runoff for bare-fallow purple soil, using 5 m × 1.5 m experimental plots under natural rainfall conditions. Four experimental plots(10°, 16°, 20° and 26°) were established in theYanting Agro-ecological Experimental Station of Chinese Academy of Science in central Sichuan Basin. The plot was equipped with water storage tank to monitor water level change. Field monitoring from July 1 to October 31, 2012 observed 42 rainfall events which produced surface runoff from the experimental plots. These water level changes were converted to runoff. The representative eight rainfall events were selected for further analysis, the relationship between slope and runoff coefficient were determined using ANOVA, F-test, and z-score analysis. The results indicated a strong correlation between rainfall and runoff in cumulative amount basis. The mean value of the measured runoff coefficient for four experimental plots was around 0.1. However, no statistically significant relationship was found between slope and runoff coefficient. We reviewed the relationship between slope and runoff in many previous studiesand calculated z-score to compare with our experimental results. The results of z-score analysis indicated that both positive and negative effects of slope on runoff coefficient were obtained, however a moderate gradient(16°-20° in this study) could be a threshold of runoff generation for many different soils including the Chinese purple soil.     

Effect of storm pattern on soil erosion in damaged rangeland; field rainfall simulation approach

Raindrop size,rainfall intensity and runoff discharge affect the detachment and transportation of soil particles.Among these three factors,the rainfall intensity seems to be more important because it can change other two factors.Storm patterns can be determined by changing the rainfall intensity during the storm.Therefore,the objective of this research is to test the influence of storm pattern on runoff,soil erosion and sediment concentration on a rangeland soil slope under field rainfall simulation.Four storm rainfall intensity patterns were selected for examining the effects of variations in storm event characteristics on soil erosion processes.The selected storm patterns were:I(45,55 and 70 mm h-1);II(45,70 and 55 mm h-1);III:(70,55 and 45 mm h-1);and IV(55,45 and 70 mm h1).The last pattern is a new one instead of the uniform pattern which has been sufficiently studied in previous researches.The experiments were conducted in field plots(in Kojour watershed,Mazandaran Province,Iran)with an area of one square meter and an constant slope gradient of 18%,surrounded by galvanised sheets.Following the nonuniform prioritization of the storm patterns for the studied variables,time to runoff(I>II>IV>III),runoff volume(III>IV>II>I),sediment concentration(IV>III>I>II)and soil erosion(III>IV>II>I)),it can be generally inferred that each pattern has specific effect on soil erosion processes during a storm.The results of the general linear model(GLM)test indicated that the effects of storm pattern on time to runoff,total runoff volume,runoff coefficient and soil erosion were significant at a level of 99%.The Duncan test showed that the storm patterns can be divided into three groups of III,IV;II;I(for time to runoff),I,II;IV,III(for runoff coefficient),and I;II;IV,III(for runoff volume and soil erosion).     

Rainfall and Tillage Impacts on Soil Erosion of Sloping Cropland with Subtropical Monsoon Climate- A Case Study in Hilly Purple Soil area,China

Under global warming, storm events tend to intensify, particularly in monsoon-affected regions.As an important agricultural area in China, the purple soil region in the Sichuan Basin, where it has a prevailing monsoon climate, is threatened by serious soil erosion. Tillage operations alter runoff and soil erosion processes on croplands by changing the physical properties of the soil surface. To clarify the relationship between tillage and soil erosion in the purple soil region, three different tillage practices in this region were investigated at the plot scale over 4 years: bare land with minimum tillage(BL),conventional tillage(CT) and seasonal no-tillage ridges(SNTR) which was initially designed to prevent soil erosion by contoured ridges and no-tillage techniques. The results showed that although there were no significant differences in the surface runoff and soil erosion among the three practices, BL causedrelatively high surface runoff and soil erosion,followed by CT and SNTR. Classification and comparison of the rainfall events based on cluster analysis(CA) verified that the surface runoff was not significantly different between most intensive event and long intensive events but was significantly different between most intensive and short and medium-duration events. Only the rainfall events with the highest rainfall intensity could trigger serious soil erosion, up to 1000 kg ha-1 in the region. Further detailed investigations on the effects of tillage operations on the soil erosion in a subtropical region with a monsoon climate are needed to provide a basis for modeling catchments and designing better management practices.