全球土壤可蚀性因子砾石效应的影响因素研究

为了全面认识土壤中的砾石对土壤可蚀性因子(K)的影响,即砾石效应(stoniness effect on erodibility,SE_K),利用全球范围的砾石覆盖影响(S_t),剖面砾石影响(K_cf)、砾石覆盖和剖面砾石综合影响(K_f-cs)以及高程、坡度、坡向、植被覆盖、土地利用类型、多年平均降雨量、多年平均气温等环境要素数据,通过相关性分析、格局分析以及随机森林回归分析,探索了SE_K的空间格局以及影响SE_K的主控因素。结果表明:(1)高程和坡度与SE_K呈正相关; 植被覆盖度、降雨量和温度与S_t和K_f-cs呈负相关; 降雨量与K_cf呈负相关; 温度与K_cf呈正相关; 阳坡的砾石效应大于阴坡; 耕地和林地的砾石效应较小,草地、荒漠和裸地的砾石效应较大。(2)S_t和K_f-cs的主控因子为高程和坡度,降雨量和温度是K_cf的主控因子。因此,在大区域土壤侵蚀评价中,对于砾石含量较大、高程较高、坡度较大的区域,均应给予特别关注,如果不加考虑,则有可能使土壤侵蚀评价结果出现偏差。     

土壤可蚀性K值的计算和K值图的制作方法研究——以南京市方便水库小流域为例

本文介绍了土壤可蚀性K值的计算和K值图的制作方法。涉及利用公式法计算该区的土壤可蚀性K值的参数设置方法,和运用地统计学的协同克里格空间插值的方法进行K值图的编制方法,并对K值图在水土保持、水土流失定量监测、生态与环境等方面的应用作了简要说明。     

工程堆积体标准小区界定与可蚀性因子改进

将USLE模型应用于工程堆积体侵蚀预报时结果偏差较大,这是因为其标准小区与土壤可蚀性因子不适合直接应用于工程堆积体侵蚀预报。针对这一问题,结合对我国6大水蚀类型区工程堆积体参数的调查、统计、分析,界定工程堆积体标准小区的坡度、坡长;结合国内外已有研究成果,分析论证堆积体可蚀性因子的改进办法,并引用相关数据验证其可行性。结果表明,6大区域工程堆积体坡度、坡长及坡面物质组成均与USLE中规定的标准小区相差较大;为了尽量消除误差,建议工程堆积体标准小区坡度采取众数35°,坡长采取平均值5m;堆积体可蚀性因子更名为土石质因子Tr,并将单位体积土石混合体中石砾总表面积Cs,作为石砾因素指标纳入堆积体土石质因子Tr中,与堆积体物质中土壤可蚀性因子K共同构建工程堆积体土石质因子函数。验证结果表明,堆积体可蚀性因子改进办法是可行的。工程堆积体坡度、坡长的界定及可蚀性因子的改进,对提高基于USLE模型的堆积体侵蚀预报精度有重要意义。     

A rational method for estimating erodibility and critical shear stress of an eroding rill

Soil erodibility and critical shear stress are two of the most important parameters for physically-based soil erosion modeling. To aid in future soil erosion modeling, a rational method for determining the soil erodibility and critical shear stress of rill erosion under concentrated flow is advanced in this paper. The method suggests that a well-defined rill be used for shear stress estimation while infinite short rill lengths be used for determination of detachment capacity. The derivative of the functional relationship between sediment yield and rill length at the inlet of rill flow, as opposed to average detachment rate of a long rill, was used for the determination of detachment capacity. Soil erodibility and critical shear stress were then regressively estimated with detachment capacity data under different flow regimes. Laboratory data of rill erosion under well defined rill channels from a loess soil was used to estimate the soil erodibility and critical shear stress. The results showed that no significant change in soil erodibility (K_r) was observed for different slope gradients ranging from 5 to 25 while critical shear stress increased slightly with the slope gradient. Soil erodibility of the loess soil was 0.3211 ± 0.001 s m^− 1. The soil erodibility and critical shear stress calculations were then compared with data from other resources to verify the feasibility of the method. Data comparison showed that the method advanced is a physically logical and feasible method to calculate the soil erodibility and critical shear stress for physically-based soil erosion models.     

石砾参数对土壤水流和溶质运移影响研究进展

土壤水流和溶质运移一直是土壤学研究的热点,溶质运移理论主要应用于地下水污染、污染物运移、土壤重金属污染研究等方面。溶质主要通过优先流和基质流进行运移,影响溶质运移因素很多,主要包括土壤结构、质地、水力传导率、体积质量、初始含水量、根系、石砾等。石砾作为土壤质地中的一个分级单位,与溶质运移关系较为复杂。本文综合介绍了石砾基本内涵以及石砾对土壤水流和溶质运移影响研究进展;系统阐述了石砾内部参数(石砾覆盖度、含量、粒径、空间异质性等)和外部参数(根石结构、干湿冻融、耕作等),指出目前研究主要量化土壤表面及土壤表层石砾参数对水文效应、土壤侵蚀、入渗以及径流的影响,然而石砾参数对溶质运移影响研究不够系统,石砾参数与溶质运移关系研究尚处于初步阶段,对土壤深层石砾研究缺乏;归纳了石砾参数研究技术手段及模型;探讨了目前石砾参数对土壤水流和溶质运移影响研究存在的问题以及今后研究趋势。     

四川自然土壤和旱耕地土壤可蚀性特征研究

应用美国通用土壤流失方程 (USLE)和土壤侵蚀预报模型 (WEPP)中的土壤可蚀性K值 ,对四川各类自然土壤和旱耕地土壤可蚀性特征进行了研究。结果表明 :土壤可蚀性K值与土壤理化性质直接相关 ,自然土壤和旱耕地土壤可蚀性K值在 0 2 68～ 0 3 44之间 ,紫色土的分布面积和K值较大 ,是易遭受侵蚀的土壤。应采取增施有机肥、实行坡改梯等措施 ,加强对耕地、高可蚀性土壤侵蚀的综合防治     

Soil erosion risk assessment using GIS and CORINE model: a case study from western Shiraz,Iran

This study employs the Coordination of Information on the Environment (CORINE) model with geographic information system to assess soil erosion risk for restoring and protecting areas within the Bonrod Zangane watershed, western Shiraz, Iran. Actual soil erosion risk was determined by combining two main parameters including potential soil erosion risk and vegetation cover. The potential soil erosion risk was generated by integrating soil erodibility, erosivity and slope parameters. Soil texture, depth and stoniness layers were overlaid to form a soil erodibility map. Modified Fournier index and Bagnouls–Gaussen aridity index were integrated to generate the erosivity layer. The slope classes also were generated from digital elevation model. In order to estimate vegetative land cover, the normalized difference vegetation index (NDVI) was used. The raster-based layers were then integrated to produce erosion risk map. The results showed that 34.7% of the study area has high and only 31.4% of the study area has low soil erosion risk. It is concluded that CORINE model can be used to delineate the soil erosion risk and also to discriminate the potential soil erosion risk areas.     

Dynamic study of infiltration rate for soils with varying degrees of degradation by water erosion

Ultisols, widely distributed in tropical and subtropical areas of south China, are suffering from serious water erosion, however, slope hydrological process for Ultisols under different erosional degradation levels in field condition has been scarcely investigated. Field rainfall simulation at two rainfall intensities (120 and 60 mm/h) were performed on pre-wetted Ultisols with four erosion degrees (non, moderate, severe and very-severe), and the hydrological processes of these soils were determined. The variation of soil infiltration was contributed by the interaction of erosion degree and rainfall intensity (p < 0.05). In most cases, time to incipient runoff, the decay coefficient, steady state infiltration rate, and their variability were larger at the high rainfall intensity, accelerating by the increasing erosion severity. Despite rainfall intensity, the infiltration process of Ultisols was also significantly influenced by mean weight diameter of aggregates at the field moisture content, soil organic carbon and particle size distribution (R² > 30%, p < 0.05). The temporal erodibility of surface soil and soil detachment rate were significantly and negatively correlated with infiltration rate (r < -0.32, p < 0.05), but less significant correlation was observed between sediment concentration and infiltration rate for most soils, especially at the high rainfall intensity. The variation of surface texture and soil compactness generated by erosion degradation was the intrinsic predominant factors for the change of infiltration process of Ultisols. The obtained results will facilitate the understanding of hydrological process for degraded lands, and provide useful knowledge in managing crop irrigation and soil erosion.     

晋北侵蚀区土壤侵蚀潜在危险度评价研究

选取人口环境容量失衡度、年降水量、植被覆盖度、土层厚度、地形坡度、土壤可蚀性、岩性和坡耕地占坡地面积比例等八个因素,作为评价晋北侵蚀区土壤侵蚀潜在危险度的主要因素,利用特尔斐法(Delphi)和层次分析法(AHP)确定8个因子的权重,运用数学模型计算,对晋北侵蚀区土壤潜在危险度进行了评价。     

广东省土壤可蚀性现状及影响因素分析

土壤可蚀性是土壤侵蚀预报和土地利用规划的重要参数,本文采用EPIC(Erosion Productivity Im-pact Calculator)模型中土壤可蚀性因子K值为指标,利用第二次土壤普查资料,探讨广东省土壤可蚀性K值及分布特征,并绘制了广东省土壤可蚀性K值图,结果表明:广东土壤可蚀性K值为0.116~0.415,加权平均K值为0.25,主要分布在较低-中高可侵蚀性范围;以铁铝土为例,成土母质对土壤侵蚀影响是多因素的,由于母质的特性差异,母质所发育土壤可蚀性K值并不能完全代表其侵蚀危害性,从总体上看,土壤经过多年耕种,抗侵蚀能力明显下降。     

Removal of rock fragments and its effect on soil loss and crop yield, Tigray, Ethiopia

Abstract. When the farmers of the Highlands of Tigray (northern Ethiopia) consider rock fragment cover in their fields to be excessive, they remove some of them. In addition, large amounts of rock fragments of all sizes are removed from fields for building stone bunds. Semi-structured interviews indicate that the farmers are often reluctant to take away the smaller rock fragments (i.e. < 5 cm across) from their fields, since they believe these benefit soil moisture conservation and protect topsoil from erosion. A field experiment was carried out on a Vertic Cambisol (average slope: 0.125 m m^–1), 2 km east of Hagere Selam (subhumid climate). Rock fragments were totally, partially or not removed from the 12 runoff plots (5 m × 6 m) before the beginning of the 1999 cropping season, during which a local mixture of wheat varieties ( Triticum spp.) was sown. After harvest, erosion rates were assessed by measuring deposited sediment volume in trenches at the lower side of each subplot, and grain and straw yields were assessed. We found a significant negative relationship between rock fragment cover and soil loss by water erosion. However, the resulting positive relationship between rock fragment cover and grain and straw yield was weak. This might be explained by the fact that the plot did not suffer from drought due to soil and climatic conditions. Detailed analysis showed that cover by medium and large rock fragments (> 2 cm diameter) showed an optimum percentage cover above which crop yields decrease. A recommendation resulting from this study is to rely on the farmers' experience: smaller rock fragments should never be removed from the surface of fields during soil and water conservation works; instead rock fragment rich soil can be used to top the stone bunds.     

Filling the European blank spot—Swiss soil erodibility assessment with topsoil samples

Soil erodibility, commonly expressed as the K‐factor in USLE‐type erosion models, is a crucial parameter for determining soil loss rates. However, a national soil erodibility map based on measured soil properties did so far not exist for Switzerland. As an EU non‐member state, Switzerland was not included in previous soil mapping programs such as the Land Use/Cover Area frame Survey (LUCAS). However, in 2015 Switzerland joined the LUCAS soil sampling program and extended the topsoil sampling to mountainous regions higher 1500 m asl for the first time in Europe. Based on this soil property dataset we developed a K‐factor map for Switzerland to close the gap in soil erodibility mapping in Central Europe. The K‐factor calculation is based on a nomograph that relates soil erodibility to data of soil texture, organic matter content, soil structure, and permeability. We used 160 Swiss LUCAS topsoil samples below 1500 m asl and added in an additional campaign 39 samples above 1500 m asl. In order to allow for a smooth interpolation in context of the neighboring regions, additional 1638 LUCAS samples of adjacent countries were considered. Point calculations of K‐factors were spatially interpolated by Cubist Regression and Multilevel B‐Splines. Environmental features (vegetation index, reflectance data, terrain, and location features) that explain the spatial distribution of soil erodibility were included as covariates. The Cubist Regression approach performed well with an RMSE of 0.0048 t ha h ha^?1 MJ^?1 mm^?1. Mean soil erodibility for Switzerland was calculated as 0.0327 t ha h ha^?1 MJ^?1 mm^?1 with a standard deviation of 0.0044 t ha h ha^?1 MJ^?1 mm^?1. The incorporation of stone cover reduces soil erodibility by 8.2%. The proposed Swiss erodibility map based on measured soil data including mountain soils was compared to an extrapolated map without measured soil data, the latter overestimating erodibility in mountain regions (by 6.3%) and underestimating in valleys (by 2.5%). The K‐factor map is of high relevance not only for the soil erosion risk of Switzerland with a particular emphasis on the mountainous regions but also has an intrinsic value of its own for specific land use decisions, soil and land suitability and soil protection.     

MINErosion 4: Using measurements from a tilting flume-rainfall simulator facility to predict erosion rates from post-mining catchments/landscapes in Central Queensland,Australia

The use of draglines to remove overburden in Queensland opencut mines, results in landscapes that consist of long parallel tertiary overburden spoil-piles that are generally highly saline, dispersive, and highly erodible. The height of these spoil-piles may exceed 50–60 m above the original landscapes and the slopes are at the angle of repose of around 75% or 37°. Legislation and public opinion require that these highly disturbed open-cut post-mining landscapes should be satisfactorily rehabilitated into an approved post-mining land use with acceptable erosion rates. Therefore, these slopes must be reduced before the landscape can be rehabilitated. The most expensive component of the rehabilitation process is the re-shaping and preparation of the overburden to create a suitable landscape for vegetation growth. As soils and overburden varies greatly in their erodibilities, the extent and cost of earthworks can be minimized, and rehabilitation failures avoided, if soil erosion from designed landscapes can be predicted using laboratory-based parameters prior to construction of these landscapes. This paper describes the development of a model for that purpose.A catchment or landscape erosion model MINErosion 4 was developed by upscaling the existing hillslope model MINErosion 3 (So, et al., 2018) and integrate it with both ESRI ArcGIS 10.3 or QGIS 3.16 (freeware), to predict event based and mean annual erosion rate from a postmining catchment or landscape. MINErosion 3 is a model that can be used to predict event and annual erosion rates from field scale hillslopes using laboratory measured erodibility parameters or routinely measured soil physical and chemical properties, and to derive suitable landscape design parameters (slope gradient, slope length and vegetation cover) that will result in acceptable erosion rates. But it cannot be used to predict the sediment delivery from catchments or landscapes. MINErosion 4 was validated against data collected on three instrumented catchments (up to 0.91 ha in size) on the Curragh mine site in Central Queensland. The agreement between predicted (Y) and measured (X) values were very good with the regression equation of Y = 0.92X and an R² value of 0.81 for individual storm events, and Y = 1.47X and an R² value of 0.73 for the average annual soil loss. This is probably the first time that a catchment scale erosion is successfully predicted from laboratory measured erodibility parameters.     

Factors responsible for soil erosion hazard in submontane Punjab, India

Abstract. Various factors of soil erosion (erodibility of the soil, slope of the land and nature of the plant cover) were studied in a representative area of submontane Punjab. In the north-eastern part of the area, slope steepness, slope length, convexo-concave and concave slope patterns and sparse vegetation were the dominant factors governing erosion hazard. In the central and lower south-western part of the area, high soil credibility and sparse vegetation were the most important factors.
Rill erosion was positively correlated ( r = 0.87) with slope steepness. The erodibility of the soil was strongly correlated ( r = 0.98) with amounts of silt + very fine sand.
In the area with steeper slopes, the lower segments of transects were prone to more erosion because of high credibility. Slope shapes such as convexo-concave and concave could result in sloughing and greater erosion hazard if not properly protected. Slopes facing southwest were more prone to erosion than others because of more solar energy, greater aridity, less vegetation and concave shape.     

中国土壤可蚀性值及其估算

土壤可蚀性是评价土壤对侵蚀敏感程度的重要指标，也是进行土壤侵蚀预报的重要参数。本文运用野外观测资料，研究了我国不同水土流失区的土壤可蚀性值问题。根据实测资料，计算给出了一组土壤可蚀性实测值。并利用这组实测值。对我国土壤可蚀性估算问题进行了探讨。结果表明，国外现有的可蚀性估算模型不能直接应用于我国土壤的可蚀性计算，估算值明显大于实测值。但估算值与实测值之间存在有良好的线形关系。最后提出了我国不同地区及不同资料占有情况下的土壤可蚀性估算方法。本文研究结果可以直接用于我国土壤侵蚀预报中土壤可蚀性计算。     

Determination of rill erodibility and critical shear stress of saturated purple soil slopes

The hydrological conditions near the soil surface influence the soil erosion process, as determined by the soil erodibility and critical shear stress. The soil erodibility and critical shear stress of saturated purple soil slopes were computed and compared with those of unsaturated purple soil slopes. The detachment capacities computed through the numerical method (NM), modified numerical method (MNM) and analytical method (AM), from rill erosion experiments on saturated purple soil slopes at different flow rates (2, 4, and 8 L min^?1) and slope gradients (5, 10, 15, and 20°), were used to comparatively compute the soil erodibility and critical shear stress. The computed soil erodibilities and critical shear stresses were also compared with those of unsaturated purple soil slopes. At the different slope gradients ranging from 5° to 20°, there were no significant differences in the soil erodibilities of the saturated purple soil and also in those of the unsaturated purple soil. The critical shear stresses slightly varied with the slope gradients. The saturated purple soil was relatively significantly more susceptible to erosion. The NM overestimated the soil erodibility of both saturated and unsaturated soils by 31% and underestimated the critical shear stress. The MNM yielded the same soil erodibility and critical shear stress values as the AM. The results of this study supply parameters for modeling rill erosion of saturated purple soil slope.     

三峡库区植物篱模式对土壤理化性质和可蚀性的影响

为探究不同植物篱模式对土壤理化性质及土壤可蚀性空间上的影响,以三峡库区秭归县张家冲水土保持试验站为研究区,选取经济林地小区(H1)和农耕地小区(H3)分别为"植物篱+经济林地小区"(H2)和"植物篱+农耕地小区"(H4)对照小区,对其3个坡位(上、中、下坡)和2个土层(0—20,20—40 cm)进行土壤理化性质对比分析,探讨2种植物篱配置下对土壤理化性质和可蚀性的影响。结果表明:(1)同一小区,土壤机械组成和土壤养分含量空间分异较明显,上坡与中坡、下坡差异显著(p0.05),其中土壤机械组成以砂粒为主(59.01%～63.51%),且分布于上坡,而土壤细颗粒主要分布在中、下坡,0—40 cm土层土壤养分均表现为下坡中坡上坡,全钾含量均值较其对照小区分别高3.75%～19.61%。(2)不同小区,土壤细颗粒和土壤养分含量表现为植物篱小区高于无植物篱小区,其中土壤细颗粒占比表现为H4小区(39.94%)H2小区(38.92%)H1小区(38.34%)H3小区(37.84%);土壤可蚀性K值大小与土壤细颗粒含量呈反比,即土壤细颗粒占比越大,K值越小,越不易被侵蚀,"农耕地+植物篱"(H4)较"经济林+植物篱配置"(H2)更不易被侵蚀。(3)随着土壤有机质、黏粒、全氮和全钾含量的增加,可以有效增强土壤抗侵蚀能力。土壤可蚀性K值与粉粒、黏粒分别呈极显著正相关和负相关关系(相关系数分别为1.000,-0.708),而与砂粒、有机质、全氮和全钾均无显著相关关系(p0.05)。     

Interrill erosion from disturbed and undisturbed samples in relation to topsoil aggregate stability in red soils from subtropical China

The assessment of soil erodibility to water erosion in the field is often expensive and time-consuming. This study was designed to reveal the effects of aggregate breakdown mechanisms on interrill erosion dynamics and develop an improved model for assessing interrill soil loss, which incorporated the soil aggregate stability tests as a substitute for the interrill erodibility parameter, from both disturbed and undisturbed samples for red soils in subtropical China. Six cultivated areas of sloping land with red soils were selected, and topsoil aggregate stability was analyzed using the Le Bissonnais method to determine the different disaggregation forces. Laboratory rainfall simulations were designed to distinguish the effects of slaking (at different wetting rates) and mechanical breakdown (with and without screening) on soil erosion characteristics. Field rainstorm simulations with medium and high rainfall intensities were conducted on runoff plots (2 m 1 m) with slope gradients varying from 10% to 20% for each soil type. A new instability index, Ka, which considers aggregate breakdown mechanisms in interrill erosion processes, was proposed based on the disturbed sample results. Ka showed a close relationship with erosion rates in both disturbed and undisturbed samples. Following from the results of undisturbed sample experiments, Ka was used as a substitute for the erodibility factor, and introduced into the WEPP model, establishing a new erosion predication formula for red soils which had a good correlation coefficient (R² = 0.89**). This research made a good attempt at estimating the interrill erosion rate on the basis of aggregate stability from simple laboratory determinations. These results extend the validity of soil aggregation characterization as an appropriate indicator of soil susceptibility to interrill erosion in red soils from subtropical China. The formula based on the instability index, Ka, has the potential to improve the methodology used for assessing interrill erosion rates.     

Effects of land use on soil erosion in a tropical dry forest ecosystem,Chamela watershed,Mexico

During the past few decades, Mexico has been converting tropical dry forest (TDF) into cropland and pasture, with land degradation expressed as soil erosion being the main environmental consequence. The factors and processes influencing soil erosion are related to scale. At a microscale, the stability of soil aggregates has a significant impact on soil erodibility and strongly influences other soil properties. However, at plot and watershed scales, these relationships are less well known. The relationships between the distribution of soil aggregate size, soil properties and soil erosion were examined for two soil geomorphological units (hillslopes over granite and hillslopes over tuffs) and three land uses (TDF, unburned pasture and burned pasture) within the Chamela watershed of west–central Mexico. To evaluate soil aggregation as a parameter for upscaling soil erosion, the researchers measured microtopographic features at plot scales and interpreted 1:35,000 panchromatic aerial photographs at a watershed scale. Analysis of variance indicated significant differences in soil organic carbon (P < 0.05) and soil moisture (P < 0.01) contents between the two soil geomorphological units, and field tests showed differences in soil texture and structure.     

Erodibility of agricultural soils on the Loess Plateau of China

Soil erodibility is thought of as the ease with which soil is detached by splash during rainfall or by surface flow. Soil erodibility is an important factor in determining the rate of soil loss. In the universal soil loss equation (USLE) and the revised universal soil loss equation (RUSLE), soil erodibility is represented by an erodibility factor (K). The K factor was defined as the mean rate of soil loss per unit rainfall erosivity index from unit runoff plots. Although high rate of soil loss from the Loess Plateau in China is well known and widely documented, it is remarkable that there is little systematic attempt to develop and validate an erodibility index for soils on the Loess Plateu for erosion prediction. Field experimental data from four sites on the Loess Plateau were analyzed to determine the K factor for USLE/RUSLE and to compare with another erodibility index based on soil loss and runoff commonly used for the region. The data set consists of event erosivity index, runoff, and soil loss for 17 runoff plots with slope ranging from 8.7 to 60.1%. Results indicate that the K factor for USLE/RULSE is more appropriate for agricultural soils on the Loess Plateau than the erodibility index developed locally. Values of the K factor for loessial soils range from 0.0096 to 0.0269 t h/(MJ mm). The spatial distribution of the K value in the study area follows a simple pattern showing high values in areas with low clay content. For the four sites investigated, the K factor was significantly related to the clay content, (K=0.031−0.0013 Cl, r²=0.75), where Cl is the clay content in percent. The measured values of the K factor are systematically lower than the nomograph-based estimates by a factor of 3.3–8.4. This implies that use of the nomograph method to estimate soil erodibility would considerably over-predict the rate of soil loss, and local relationship between soil property and the K factor is required for soil erosion prediction for the region.