Soil erosion prediction and sediment yield estimation: the Taiwan experience

Estimating watershed erosion using geographic information systems coupled with the universal soil loss equation (USLE) or agricultural non-point source pollution model (AGNPS) has become a recent trend. However, errors in over-estimation often occur due to the misapplication of parameters in the equation and/or model. Because of poor slope length calculation definitions for entire watersheds, the slope length factor is the parameter most commonly misused in watershed soil loss estimation. This paper develops a WinGrid system that can be used to calculate the slope length factor from each cell for reasonable watershed soil loss and sediment yield estimation.     

草本植被覆盖对坡面降雨径流侵蚀影响的试验研究

通过野外人工模拟降雨试验,研究了草本植被覆盖对坡面降雨径流侵蚀的影响,并从径流侵蚀功率和降雨侵蚀力两个方面对比分析了草本植被对坡面侵蚀动力的调控效果,结果表明草本植被覆盖深刻影响降雨侵蚀动力,并最终对坡面径流侵蚀量产生较大的影响。植被覆盖度为0%～60%时,产流产沙量随植被覆盖度的增加迅速降低,植被覆盖度>80%时,覆盖度的增加不能引起产流、产沙量的大幅度下降,植被水沙调控作用趋于稳定,确定本研究的临界植被覆盖度为60%～80%;以径流深和洪峰流量模数表示的坡面径流侵蚀功率以及降雨侵蚀力等侵蚀动力指标均与侵蚀产沙量呈正相关关系,但径流侵蚀功率与产沙量具有更强的相关性,说明径流侵蚀功率能更好地模拟侵蚀动力;以径流侵蚀功率/侵蚀量表示植被覆盖度对侵蚀结果的影响,反映了临界植被覆盖度的存在,可以作为评价植被侵蚀动力调控效应的一个指标。     

Modelling of snowmelt erosion and sediment yield in a small low-mountain catchment in Germany

Temporal variability and spatial heterogeneity of surface runoff generation triggers the dynamics of source areas of sediment and sediment-associated nutrient transport. Reliable modelling of hydrological special situations i.e. snowmelt is of high importance for the quality of erosion and sediment yield modelling. Data from the research catchment Schäfertal demonstrate the individuality of snowmelt events in terms of runoff coefficient and delivery ratio. This 1.44 km² low mountain catchment is characterised by a high portion of arable land with a winter grain/winter rape crop rotation. The integrated winter erosion and nutrient load model (IWAN) considers these dynamic aspects by coupling a hydrological model with a sediment load model. Cell size of this raster-based approach is 10 × 10 m². Additionally, snowmelt rill erosion is simulated with a newly developed physically based model that is firstly applied on a catchment scale. A sensitivity analysis of this model system component demonstrates the plausibility of the model approach and the overall robustness of the model system IWAN. The results of the long-term hydrological modelling from 1991 to 2003 are reliable and form the basis for the simulation of six snowmelt events which were observed in the Schäfertal catchment. The estimated total runoff volumes for these events match the observations well. The modelled overland runoff coefficients vary from 0.001 to 0.72. The mean values of cell erosion, which were modelled with one set of parameters for all six events range from 0.0006 to 0.96 t ha^− 1. The total modelled erosion for the events with unfrozen soil and low amount of surface runoff is of a factor 50 below those with partly frozen soil. In addition to these distinctions, the major differences are caused by flow accumulation in shallow depressions in variable parts of the catchment. However, the validation of these results on the single event scale is restricted due to limited spatial data. Total simulated sediment yield at the catchment outlet was as high as 13.84 t which underestimates the observed values, with the exception of one event. Oversimplification of the modelled channel processes may be a reason. The temporal variability and spatial heterogeneity of the surface roughness parameter, which was identified to be sensitive, also causes uncertainty in the parameter estimation. Despite these findings, the model system IWAN was applied successfully on the catchment scale and the simulated results are reliable.     

Geomorphological investigation on gully erosion in the Rift Valley and the northern highlands of Ethiopia

Gully erosion phenomena are very common in Ethiopia. They affect large areas with different morphological, pedological and climatic characteristics. The amount of soil loss due to gullying has become a very serious problem in the recent decades as it was associated to remarkable depletion of cultivated land. Field investigations on gully morphology and its genetic processes were carried out in two study areas of Ethiopia, representative of different geo-environmental conditions: the Lakes Region in the Rift portion north of Shashamene and the area surrounding the town of Mekele in Tigray.Two main types of gullies were identified on the basis of their morphological and hydraulic geometry characteristics: (1) discontinuous gullies which generally develop on low gradient slopes (1–5° on average) and the hydraulic radius of which increases from an upstream minimum to a maximum, at approximately their mid length, and decreases again to a relative minimum at their downstream end; and (2) stream gullies, formed by deep erosion processes typically migrating upslope. In order to investigate the main causes originating the different types of gullies, shear stress data were collected in the field from their hydraulic geometry. Hypotheses on the mechanisms responsible for both discontinuous and stream gullies development and for their different characteristics are discussed considering the pattern of shear stress variation in the downstream direction.     

非恒定流对浅沟侵蚀产沙及水动力参数的影响

浅沟侵蚀是黄土高原独特的侵蚀类型,在坡面侵蚀中占有重要的地位。关于浅沟侵蚀的研究,以往多采用恒定流量进行冲刷试验,而在实际中流量是变化的,有关非恒定流的浅沟侵蚀实验尚未见报道。采用野外放水冲刷试验,定量分析10～15、15～20、20～25 L/min 3个非恒定流的放水流量级和10、15、20 L/min 3个恒定流量的放水流量级对14°坡耕地坡面浅沟侵蚀产沙及其水动力参数的影响。结果表明:恒定流量级的流速、雷诺数、弗劳德数、剪切力、水流功率和径流含沙量分别较变流量后的非恒定流小6.64%～13.44%、2.73%～8.26%、6.15%～49.59%、0.27%～10.27%、4.27%～15.77%和2.02%～19.36%;恒定流的阻力系数则较与其相对应的非恒定流增大23.45%～43.41%。     

区域尺度侵蚀产沙估算方法研究

 区域土壤侵蚀模型是大区域土壤侵蚀普查和水土保持宏观决策的支持工具,土壤侵蚀模型的研发是土壤侵蚀学科的前沿领域。基于DEM将区域划分为规则网格,设计产流、产沙过程的单元模型,包括植被截留、入渗、填洼、流速、携沙能力、径流剥蚀量、泥沙沉积等算法。将月降水当作1次降雨事件,并划分若干时段进行迭代计算,利用GIS空间分析功能完成水沙汇集运算,并在ArcGIS支持下进行计算机程序设计,有效地完成了区域侵蚀产沙量的计算。将模型应用于延河流域得到:1995年7月份平均径流深为35.6mm,径流系数为0.237,流域出口径流量为2.72亿m3,流域出口输沙量为0.38亿t,流域平均侵蚀模数为4575t/(km2.月);输出图形空间格局和结构符合实际情况,初步模拟结果令人较满意。     

上方来水来沙对浅沟侵蚀产沙及水动力参数的影响

 浅沟侵蚀是黄土高原重要的侵蚀类型,上方汇水对坡面浅沟侵蚀具有重要的影响。采用野外放水冲刷试验,定量分析26°坡耕地在上方来水量为5、10和15L/min时对坡下方浅沟侵蚀产沙及其水动力参数的影响。结果表明：上方来水的汇入使浅沟水流流速明显增大,雷诺数、弗劳德数、水流功率和剪切力分别增大33%～76%、21%～47%、29%～72%和18%～42%,阻力系数减少11%～13%,导致浅沟侵蚀产沙量明显增大;除流速和弗劳德数外,其余水动力参数随放水时间的延长呈递增趋势;上方来水使浅沟侵蚀产沙量相对增量与水流功率和剪切力相对增量均呈幂函数关系。     

Assessing interrill erosion rate from soil aggregate instability index, rainfall intensity and slope angle on cultivated soils in central Greece

In order to develop a new formula for assessing interrill erosion rate by incorporating the soil aggregate instability index, β, erosion plots at seven sites in central Greece were used to measure interrill erosion rate under natural rainfall conditions during a 39-month period. Soils classified as Alfisols, Inceptisols and Entisols with slopes 7–21%, moderately well to excessively drained, clay to loamy textured, were studied. Runoff and total sediment were collected after each ponding rainfall event. The equation E_i=0.628 β S_t^1.3 e^{0.0967I₃₀} was finally proposed (R²=0.939,P<0.001) to describe interrill erosion rate. The term, S_t represents the tangent of the slope angle, and I₃₀ represents the maximum rainfall intensity in 30 min. The addition of the aggregate instability index to improve existing methodologies provide was considered to provide an easy to determine and reliable measure of soil erodibility. Validation with independent data showed that the model predicted interrill erosion well (R²=0.766, P<0.001). Therefore, the proposed model based on the aggregate instability index, β, has the potential to improved methodology for assessing interrill erosion rate.     

Modeling the impacts of no-till practice on soil erosion and sediment yield with RUSLE, SEDD, and ArcView GIS

The revised universal soil loss equation (RUSLE), the sediment delivery distributed (SEDD) model, and ArcView GIS were used to estimate the impacts of no-till practice on soil erosion and sediment yield in Pataha Creek Watershed, a typical dryland agricultural watershed in southeastern Washington. The results showed that the average cell soil loss decreased from 11.09 to 3.10 t/ha yr for the whole watershed and from 17.67 to 3.89 t/ha yr for the croplands under the no-till scenario. Likewise, the average cell sediment yield decreased from 4.71 to 1.49 t/ha yr for the entire watershed and from 7.11 to 1.55 t/ha yr for the croplands under no-till practices. The major reason why no-till practice can significantly reduce the soil erosion and sediment yield is that it prevents rill generation which through rill erosion ultimately contributes up to 90% of the soil erosion in the Inland Pacific Northwest region.     

Relationships between rainfall characteristics and ephemeral gully erosion in a cultivated catchment in Sicily (Italy)

Recent research has shown a lack of long-term monitoring for detailed analysis of gully erosion response to climate characteristics. Measures carried out from 1995 to 2007 in a wheat-cultivated area in Raddusa (Sicily, Italy), represent one of the longest series of field data on ephemeral gully, EG, erosion. The data set collected in a surface area of almost 80 ha, permits analysis of the influence of rainfall on EG formation and development. Ephemeral gullies formed in the study area were measured on a yearly scale with a Post-Processing Differential GPS for length and with a steel tape for the width and depth of transversal sections. Ephemeral gully formation was observed for 8 years out of 12, which corresponds to a return period of 1.5 years. The measurements show strong temporal variability in EG erosion, in agreement with the rainfall characteristics. The total eroded volumes ranged between 0 and ca. 800 m³ year⁻¹, with a mean of ca. 420 m³ year⁻¹, corresponding to ca. 0.6 kg m⁻² year⁻¹. Ephemeral gully erosion in the study area is directly and mainly controlled by rainfall events. An antecedent rainfall index, the maximum value of 3-days rainfall (H_max3_d), is the rain parameter which best accounts for EG erosion. This index is used here as a simple surrogate for soil water content. An H_max3_d threshold of 51 mm was observed for EG formation. The return period of the H_max3_d threshold is almost the same as the return period for EG formation. Although a mean of seven erosive rain events were recorded in a year, EG formation and development generally occur during a single erosive event, similarly to other semiarid environments. The most critical period is that comprised between October and January, when the soil is wetter and the vegetation cover is scarce. Empirical models for EG eroded volume estimation were obtained using the data set collected at this site. A simple power-type equation is proposed to estimate the eroded volumes using H_max3_d as an independent variable. This equation shows an R² equal to 0.67 and a standard error of estimation of 0.79.     

Runoff and sediment yield in relation to precipitation,temperature and glaciers on the Tibetan Plateau

The riverine sediment is an essential carrier for nutrients and pollutants delivery and is considered as an important indicator of land degradation and environmental changes. With growing interest in environmental changes over the Tibetan Plateau, this study investigated mean annual runoff and sediment yield from eight headwater catchments in relation to annual precipitation, air temperature, and glacier area ratio, etc. Results show that runoff (Q) is positively correlated with both precipitation (P) and temperature (T), i.e., Q = 0.357P+20.3T-6.4, indicating combined water supply from rainfall and meltwater, increase of which may exceed the evapotranspiration water loss caused by temperature raise. Sediment yield (S) shows an inverted parabolic relationship with precipitation and at the same time positive correlation with glacier area ratio (Ag), i.e., S = 0.000609 P²-0.470P+48.5 A_g+202.53, indicating that sediment yield is a minimum at about 500–600 mm of precipitation, increasing sharply on both sides of this minimum in one case owing to decreased vegetation protection and in the other to enhanced erosive power and that erosion rate in the glacierized area is generally higher than non-glacierized area. The variation in sediment yield with precipitation can be explained by the operation of two factors, i.e., rainfall erosive action that increases continuously with increase in precipitation, and vegetation protective action that is unity for zero precipitation and decreases with increases in precipitation. The above results may be useful in visualizing not only variations in rates of erosion among climatic zones on the Tibetan Plateau but also the probable changes of erosion during a climatic change.     

Rill formation and evolution caused by upslope inflow and sediment deposition on freshly tilled loose surfaces

Rill erosion easily occurs on tilled surfaces because the soil shear resistance is less than the runoff shear stress. However, rill erosion formation and evolution on tilled surfaces under upslope inflow conditions remain unclear. The objective of this study was to investigate the rill formation process and rill erosion amount on tilled surfaces via flow scouring under different inflow rates (4, 6, and 8 L min⁻¹) at a 15° slope. Close-range photogrammetric technology was applied to measure the rill morphology during the experiments. The results suggested that the rill formation process due to inflow could be divided into two distinct stages, namely, before and after runoff reached the downslope end, i.e., stages I and II, respectively. At stage I, runoff washed soil particles along the downslope direction. In this process, due to limited transport capacity of runoff, washed soil particles were deposited at the runoff head and formed a soil mound, which blocked the flow path, after which the runoff direction changed within the 6.7°− 50.6° range with an average moving distance of 0.62 m. As a result, a curved rill and a series of soil mounds were left on the surface. The inflow rate affected rill morphology by influencing the runoff direction change angle and runoff moving distance between the mounds on the surface. Herein, the rill formation process is referred to as downslope-trending erosion by inflow (DTEI). At stage II, runoff reached the downslope end, and a rill channel was formed throughout the slope. Thereafter, DTEI was largely reduced, and headward erosion was strengthened. As a result, the rill morphology quickly changed, such as the rill depth and rill width, which gradually increased with ongoing headward erosion. During DTEI, the rill paths were curved due to sediment deposition, and the sediment deposition conditions varied under the different inflow rates. Therefore, the rill curvature (RC) differed (1.039 ± 0.014). The RC decreased with headward erosion progression at stage II. The total sediment yield (TSY) increased with increasing inflow rate. Under inflow rates ranging from 4 to 6 L min⁻¹ and 6–8 L min⁻¹, the TSY increased 2.1–2.4 times. Consequently, DTEI on tilled surfaces significantly affects the initial rill morphology and evolution at the later stage. Hence, on slopes, its role should be considered in rill erosion assessment.     

坡度和雨强对崩岗崩积体侵蚀泥沙颗粒特征的影响

不同侵蚀条件下崩积体的侵蚀产沙特性是阐明崩积体侵蚀机理的关键。采用人工模拟降雨试验,研究不同坡度和雨强条件下崩积体坡面侵蚀泥沙颗粒的变化特征。结果表明:随着雨强和坡度的增大,泥沙粗颗粒含量及粗颗粒的富集率均增加;侵蚀物质随降雨过程逐渐变粗,后趋于稳定,大雨强条件下细沟侵蚀阶段表现为对供试土壤的"整体搬运";侵蚀泥沙颗粒的平均重量直径(Mean weight diameter,MWD)随雨强的增大而增大,1.00 mm min-1和1.33 mm min-1雨强下,细沟间及细沟侵蚀泥沙的MWD随坡度变化均存在临界坡度(30°~35°之间),其他雨强条件下则无此种情况;雨强对侵蚀泥沙MWD的影响大于坡度。     

Modelling soil erosion and sediment yield at a catchment scale: the case of Masinga catchment,Kenya

Development of improved soil erosion and sediment yield prediction technology is required to provide catchment stakeholders with the tools they need to evaluate the impact of various management strategies on soil loss and sediment yield in order to plan for the optimal use of the land. In this paper, a newly developed approach is presented to predict the sources of sediment reaching the stream network within Masinga, a large‐scale rural catchment in Kenya. The study applies the revised universal soil loss equation (RUSLE) and a developed hillslope sediment delivery distributed (HSDD) model embedded in a geographical information system (GIS). The HSDD model estimates the sediment delivery ratio (SDR) on a cell‐by‐cell basis using the concept of runoff travel time as a function of catchment characteristics. The model performance was verified by comparing predicted and measured plot runoff and sediment yield. The results show a fairly good relationship between predicted and measured sediment yield (R²=0·82). The predicted results show that the developed modelling approach can be used as a major tool to estimate spatial soil erosion and sediment yield at a catchment scale. Copyright © 2006 John Wiley & Sons, Ltd.     

Using Cs measurements to validate the application of the AGNPS and ANSWERS erosion and sediment yield models in two small Devon catchments

Distributed erosion and sediment yield models are being increasingly used for predicting soil erosion and sediment yields in agricultural catchments. In most applications, validation of such models has commonly been restricted to comparison of the predicted and measured sediment output from a catchment, because spatially distributed information on rates and patterns of soil redistribution within the catchment has been lacking. However, such spatially distributed data are needed for rigorous model testing, in order to validate the internal functioning of a model and its applicability at different spatial scales. The study reported in this paper uses two approaches to test the performance of the agricultural non-point source pollution (AGNPS) and areal non-point source watershed environmental response simulation (ANSWERS) erosion and sediment yield models in two small catchments in Devon, UK. These involve, firstly, comparison of observed and predicted runoff and sediment output data for individual storm events monitored at the basin outlets and, secondly, information on the spatial pattern of soil redistribution within the catchments derived from ¹³⁷Cs measurements. The results obtained indicate that catchment outputs simulated by both models are reasonably consistent with the recorded values, although the AGNPS model appears to provide closer agreement between observed and predicted values. However, the spatial patterns of soil redistribution and the sediment delivery ratios predicted for the two catchments by the AGNPS and ANSWERS models differ significantly. Comparison of the catchment sediment delivery ratios and the pattern of soil redistribution in individual fields predicted by the models with equivalent information derived from ¹³⁷Cs measurements indicates that the AGNPS model provides more meaningful predictions of erosion and sediment yield under UK conditions than the ANSWERS model and emphasises the importance of using information on both catchment output and sediment redistribution within the catchment for model validation.     

延河流域极端暴雨下侵蚀产沙特征野外观测分析

针对黄土高原目前侵蚀环境明显改善条件下能否经受得住极端暴雨事件的考验及可能出现的问题,该文依据在延河6个典型小流域的植被调查及土壤侵蚀监测结果,分析了2013年极端暴雨条件下延河流域不同尺度的侵蚀产沙特征。结果表明,天然乔木(辽东栎和三角槭)植被的平均侵蚀强度1 000 t/km~2,自然灌木植被在1 118.1~1 161.2 t/km~2之间,演替中后期的草本植被在1 245.2~1 827.8 t/km~2之间,而演替中前期的草本植被在3 087.6~4 408.4 t/km~2之间;人工灌木林(2 119.7~2 183.9 t/km~2)比人工乔木林(2 625.7~5 149.6 t/km~2)具有较强的减蚀能力。坡面良好的植被覆盖能有效抵御暴雨侵蚀,小流域坡面平均侵蚀强度基本2 500 t/km~2;但滑坡侵蚀占主导地位,占各小流域侵蚀总量的49.0%~88.5%,特别是距暴雨中心较近流域的南部2个小流域的滑坡侵蚀可达7 290.3和7 424.9 t/km~2左右。不同小流域内有6.4%~18.3%的侵蚀量淤积在沟道里;延河甘谷驿以上不同河段的淤积量为13.2~145.2万t,平均占总侵蚀产沙量的15%。延河甘谷驿控制区的侵蚀产沙强度变化在2 326.7~7 774.4 t/km~2之间,其中强度为5 000~8 000 t/km~2的面积占59.3%。因此,加强沟间地径流蓄排措施,减轻坡面径流下沟对沟坡重力侵蚀的影响,是土壤侵蚀研究与防治中值得重视的问题。     

Decoupling effects of driving factors on sediment yield in the Chinese Loess Plateau

Investigations of runoff and sediment yield changes and their relationships with potential driving factors provide good insights for understanding the mechanisms of hydrological processes. This study attempted to present a comprehensive investigation on the spatiotemporal variations of sediment yield in the Loess Plateau using continuous observed data at 46 hydrological stations during 1961–2016, and its responses to changes of precipitation, land use/cover and vegetation cover were analyzed by using the Partial Least Squares-Structural Equation Model (PLS-SEM). The results indicated that sediment yield reduced pronouncedly during 1961–2016 in the Loess Plateau, and 77.9% of this variation was explained by the combined effects of precipitation, land-use change, vegetation dynamics and runoff reduction. Indirect effects of precipitation, land-use change, and vegetation cover on sediment yield were 0.242, ?0.528 and ?0.630 (P < 0.05), respectively, and direct effect of runoff on sediment yield was 0.833 (P < 0.05). According to the Pearson Correlation Coefficient, the strongest positive correlation existed between annual sediment yield and runoff (r = 0.88, P < 0.05), followed by vegetation cover (r = ?0.47, P < 0.05) and land-use change (i.e. forest land and grassland) suggesting their significant trapping effects on soil erosion. However, lower correlations were examined between sediment yield and precipitation indices (?0.14<r < 0.34), and a relatively higher relationship was examined between sediment yield and heavy rainfall (P₂₅) (r = 0.34). Overall, changes in runoff and land-use/vegetation cover well explained variations in sediment yield in the Loess Plateau. The findings are expected to provide scientific and technical support for future soil and water conservation planning in the Loess Plateau, and are valuable for sustainable water resources and sediment load management in the Yellow River Basin.     

鲁中山区小流域坡面土壤侵蚀产流、产沙及侵蚀方式变化过程的研究

采用模拟降雨的方法对鲁中山区坡面土壤降雨条件下产流、产沙和侵蚀方式的演变过程进行了研究,结果表明:产流量都是随时间递增最终达到稳定产流,但是各次降雨达到稳定产流的时间均随降雨强度和坡度的增大而递减。相同坡度条件下120mmh-1雨强产流较60mmh-1雨强先到稳态,而相同雨强条件下10°坡面则较5°坡面先达到稳态。产沙量呈递增趋势且逐渐达到稳定产沙,达到稳定产沙的时间随雨强和坡度的增大而递减。相同坡度条件下120mmh-1雨强处理较60mmh-1雨强处理先达到稳态,相同雨强条件下10°坡面比5°坡面产沙过程先达到稳定产沙状态;在降雨过程中坡面首先产生面蚀,后出现细沟侵蚀;面蚀率随降雨历时延长有减小趋势,但变幅较小,细沟侵蚀率随降雨时间先急剧增加后保持稳定。     

Natural and anthropogenic controls on soil erosion in the Internal Betic Cordillera (southeast Spain)

Soil erosion in southeast Spain is a complex process due to strong interactions between biophysical and human components. Significant progress has been achieved in the understanding of soil hydrological behavior, despite the fact that most investigations were focused on the experimental plot scale. Although experimental plots allow exploring the effect of multiple biophysical and anthropogenic factors, they provide limited insights in the combined effect of all factors acting together at the landscape scale. In this study, area-specific sediment yields (SSY) have been estimated based on the volume of sediment trapped behind 36 check dams in the southeast of Spain. Low SSY-values were reported (mean = 1.40 t ha⁻¹ year⁻¹: median = 0.61 t ha⁻¹ year⁻¹). SSY variability could be explained for 67% by catchment characteristics such as drainage area, soil characteristics, land cover, average catchment slope, and annual rainfall. The low SSY values are probably caused by the agricultural abandonment that occurred over the past decades and allowed the recovery of natural vegetation. Furthermore, our results suggest that the soils have eroded in the past to such an extent that nowadays not much sediment is detached by overland flow due to residual enrichment of clay and stones. Also, sediment is to a large extent trapped locally in the catchment, as indicated by the negative relationship between SSY and catchment area.     

非硬化路面侵蚀产沙规律野外模拟试验

为了科学治理非硬化路面水土流失,采用野外放水冲刷试验,对非硬化路面侵蚀产沙规律进行了试验研究。结果表明,同一坡度不同放水流量与同一放水流量不同坡度条件下,非硬化路面产沙随时间的变化形式均有3种,分别为平缓型、多峰型和单峰型。对同一路面不同观测断面的产沙量分析发现,产沙量沿坡面的空间变化形式分别为波动式减小、逐渐式减小和先增大后减小。、得出平均含沙量与放水流量呈对数相关,产沙率与放水流量、平均含沙量与坡度、产沙率与坡度均呈直线相关。以期为非硬化路面的水土流失预测和治理提供科学依据。