Spatially distributed modelling of soil erosion and sediment yield at regional scales in Spain

Initiated by the need to quantify erosion rates and the impacts of global changes on erosion, several attempts have been made to apply erosion models at regional scales. However, these models have often been directed towards on-site soil erosion estimates, emphasising sheet and rill erosion processes, and disregarding gully erosion, channel erosion and sediment transport. These models are therefore of limited use for the assessment of sediment yield, off-site impacts of erosion, and for the development of environmental management to control these impacts at regional scale. This study analyses and compares three spatially distributed models for the prediction of soil erosion and/or sediment yield at regional scales: the WATEM-SEDEM model that is based on the empirical Revised Universal Soil Loss Equation (RUSLE) in combination with a sediment transport equation, the physics-based Pan European Soil Erosion Risk Assessment model (PESERA), and a newly developed Spatially Distributed Scoring model (SPADS). The three models were applied to 61 Spanish drainage basins and model predictions were evaluated against data on measured reservoir sedimentation rates. Global data sets on land use, climate, elevation and soil characteristics were used as model input for WATEM-SEDEM and SPADS, whereas published soil erosion estimates of PESERA at 1 km² resolution were used directly. SPADS and WATEM-SEDEM provided best results after separate calibration for basins with a Sediment Delivery Ratio (SDR) higher than 5% and those with an SDR lower than 5%. In this way, SPADS explained 67% of variation in sediment yield, while WATEM-SEDEM explained 48% of the variation. PESERA represents a promising alternative to the use of empirical models at the regional scale as it can be applied to very diverse environments with little calibration. However, PESERA provides soil erosion rates and not sediment yield estimates. For most basins PESERA soil erosion rates vary between fifty and close to zero percent of total sediment yield. Two major factors may explain this discrepancy between modelled soil erosion rates and measured sediment yield. First, it may be that PESERA underestimates soil erosion under Mediterranean conditions, although PESERA soil erosion rates are of the same order of magnitude as erosion rates measured in erosion plot studies. Second, gully-, river channel erosion and sediment transport processes may be much more important than sheet- and rill erosion for regional scale sediment yield in these environments. These issues therefore require further attention in future model development. Although spatially lumped models provide better predictions of sediment yield at the basin scale, and while validation of the predicted spatial patterns of sources and sinks of sediment requires further research, spatially distributed models are expected to be of value to support management decisions regarding the assessment of on-site and off-site impacts of erosion at the regional scale.     

黄土高原沟谷地貌发育演化研究进展与展望

黄土沟谷是黄土地貌中最有活力、最具变化、最富特色的对象单元,黄土高原千沟万壑的地貌形态以及触目惊心的侵蚀状态也让区域内沟谷地貌的形成、发育及演化问题成为研究中焦点及前沿性科学问题。近年来,诸多学者采用地学测年法、特征表达法、监测模拟法力图实现对黄土沟谷发育演化进程中“过去-现代-未来”的科学认知。这些研究在相当程度上丰富了黄土沟谷发育过程的认知。本文梳理了黄土高原沟谷地貌演化相关研究的现状,并从黄土高原地貌演化、黄土沟谷发育、基于DEM的沟谷信息提取与表达等研究进行了系统的回顾、梳理与分析。此外,本文提出“黄土沟道剖面群组”概念与方法,试图从新的视角审视黄土沟谷地貌发育演化过程。沟道剖面在黄土沟谷发育演化进程中传递物质能量和累积地形动力,并通过径流节点的串联实现剖面群的连接与组合,形成独特的剖面“群组”模式;该沟道剖面群组是集黄土沟谷地貌特征与过程于一体的综合信息集成体,其三维空间结构是对黄土沟谷地貌发育演化的高度抽象与映射,并可望进一步丰富黄土高原数字地形分析理论与方法体系,为黄土高原黄土地貌成因机理与空间分异格局带来创新的认识。     

DEM空间插值方法对土壤侵蚀模拟的影响研究——以USPED分析干热河谷典型冲沟为例

为探索不同空间插值方法得到的DEM如何影响土壤侵蚀模拟效果,本文选择金沙江干热河谷区典型冲沟为研究对象,利用野外测量高程数据,采用反距离加权(IDW)、析取克里格(DK)、局部多项式(LPI)和张力样条函数(ST)4种方法构建高精度DEM。基于USPED模型模拟冲沟的土壤侵蚀,对比不同空间插值方法的精度、土壤侵蚀的空间分布,采用相对差系数对比不同插值在土壤侵蚀研究中的相似性。结果表明：DEM空间插值的精度排序为ST相似文献   

Threshold criterion for debris flow initiation in seasonal gullies

A series of flume experiments were done to investigate the effect of grain composition on the critical gradient and discharge of debris flows initiated in seasonal gullies. The results indicated that the critical gradient and discharge for debris flow initiation decrease initially, and then increase as the mass content of fine particles (<2 mm) increases. As the mass content of fine particles increases, the angle of repose, permeability of widely graded gravel soils, and the incipient motion conditions of the coarse grains in non-uniform sediments decrease at first, and then increase. The mass content of fine particles of all inflection points is the same. The theoretical model based on the combination of hydrodynamic force and shear stress is more applicable to the prediction of the critical gradient for debris flow initiation. The critical discharge model considering the effect of non-homogeneity of the soil and the size of coarse and fine grains provides a more accurate prediction of debris flow initiation than other models based on the mean diameter.     

Soil erosion in the Anthropocene: Research needs

Soil erosion is a geomorphological and, at the same time, a land degradation process that may cause environmental and property damage, loss of livelihoods and services as well as social and economic disruption. Erosion not only lowers soil quality on‐site, but causes also significant sediment‐related problems off‐site. Given the large number of research papers on this topic, one might therefore conclude that we know now almost everything about soil erosion and its control so that little new knowledge can be added. This conclusion can be refuted by pointing to some major research gaps. There is a need for more research attention to (1) improved understanding of both natural and anthropogenic soil erosion processes and their interactions, (2) scaling up soil erosion processes and rates in space and time, and (3) innovative techniques and strategies to prevent soil erosion or reduce erosion rates. This is illustrated with various case studies from around the world. If future research addresses these research gaps, we will (1) better understand processes and their interactions operating at a range of spatial and temporal scales, predict their rates as well as their on‐site and off‐site impacts, which is academically spoken rewarding but also crucial for better targeting erosion control measures, and (2) we will be in a better position to select the most appropriate and effective soil erosion control techniques and strategies which are highly necessary for a sustainable use of soils in the Anthropocene. Copyright © 2017 John Wiley & Sons, Ltd.     

Assessment of soil factors controlling ephemeral gully erosion on agricultural fields

The soil factor is crucial in controlling and properly modeling the initiation and development of ephemeral gullies (EGs). Usually, EG initiation has been related to various soil properties (i.e. sealing, critical shear stress, moisture, texture, etc.); meanwhile, the total growth of each EG (erosion rate) has been linked with proper soil erodibility. But, despite the studies to determine the influence of soil erodibility on (ephemeral) gully erosion, a universal approach is still lacking. This is due to the complex relationship and interactions between soil properties and the erosive process. A feasible soil characterization of EG erosion prediction on a large scale should be based on simple, quick and inexpensive tests to perform. The objective of this study was to identify and assess the soil properties – easily and quickly to determine – which best reflect soil erodibility on EG erosion. Forty‐nine different physical–chemical soil properties that may participate in establishing soil erodibility were determined on agricultural soils affected by the formation of EGs in Spain and Italy. Experiments were conducted in the laboratory and in the field (in the vicinity of the erosion paths). Because of its importance in controlling EG erosion, five variables related to antecedent moisture prior to the event that generated the gullies and two properties related to landscape topography were obtained for each situation. The most relevant variables were detected using multivariate analysis. The results defined 13 key variables: water content before the initiation of EGs, organic matter content, cation exchange capacity, relative sealing index, two granulometric and organic matter indices, seal permeability, aggregates stability (three index), crust penetration resistance, shear strength and an erodibility index obtained from the Jet Test erosion apparatus. The latter is proposed as a useful technique to evaluate and predict soil loss caused by EG erosion. Copyright © 2018 John Wiley & Sons, Ltd.     

Hydraulic properties of concentrated flow of a bank gully in the dry‐hot valley region of southwest China

To quantify spatiotemporal variation in hydraulic properties of bank gully concentrated flow, a series of scour experiments were run under water discharge rates ranging from 30 to 120 l min^?1. Concentrated flows were found to be turbulent and supercritical in the upstream catchment area and downstream gully beds. As discharge increased, values of the soil erosion rate, Reynolds number (Re), shear stress, stream power, and flow energy consumption (ΔE) increased while values of the Froude number (Fr) and the Darcy–Weisbach friction factor (resistance f ) did not. With the exception of gully headcut collapse under discharge rates of 60, 90, and 120 l min^?1, a declining power function trend (P < 0.05) in the soil erosion rate developed in the upstream catchment area, headcuts, and downstream gully beds. However, increasing trends were observed in temporal variations of hydraulic properties for downstream gully beds and the upstream catchment area. Despite significant differences in temporal variation between the soil erosion rate and hydraulic property values, relative steady state conditions of the soil erosion rate and ΔE were attained following an initial period of adjustment in the upstream catchment area, headcuts, and downstream gully beds under different discharge rates. A logarithmic growth of flow energy consumption per unit soil loss (ΔE_u) was observed in bank gullies and the upstream catchment area as the experiment progressed, further illustrating the actual reason behind the discrepancy in temporal variation between soil erosion rates and ΔE. Results demonstrate that ΔE can be used to estimate headcut erosion soil loss, but further quantitative studies are required to quantify coupling effects between hydraulic properties and vertical variation in soil mechanical properties on temporal variation for bank gully soil erosion rates. Copyright © 2015 John Wiley & Sons, Ltd.     

Uncertainty assessment of ephemeral gully identification,characteristics and topographic threshold when using aerial photographs in agricultural settings

Manual digitizing on aerial photographs is still commonly used for characterizing gully erosion over large areas. Even when automated detection procedures are implemented, manual digitizing is frequently being resorted to in order to constitute reference datasets used for training and validation. In both cases, manual digitizing entails some subjective decisions on behalf of the operator, which introduces uncertainty into the resulting datasets. To assess the magnitude of this uncertainty, 11 experienced operators were asked to digitize and classify ephemeral gullies (EGs) on cropland following a standardized methodology. The resulting 11 datasets were compared in terms of number, type and location of EGs. Furthermore, for EGs located on a well‐defined runoff flow concentration axis, the slope versus contributing area topographic thresholds required for initiating gully channels were assessed using four thresholding methods, and compared across the 11 datasets. The operators identified 259 different EGs. However, the number (52–139) and sum total length (8.9–23.7 km) of EGs varied widely across operators. Only 34% of the EGs were digitized by more than half of the operators, and 7% were identified by all. Identification of EGs located on a well‐defined flow concentration axis proved least subjective. The longer the EG and the more fields the EG crossed, the larger the number of operators that were able to identify it. EGs were also most easily identified when located in sugar beet fields as compared to other crops. EG classification and topographic threshold lines were also found to be strongly operator‐dependent. Quantile regression appeared to be one of the most robust thresholding methods. Operator subjectivity when digitizing EGs on orthophotographs introduces uncertainty that should be taken into account in future remote sensing‐based studies of EG erosion whenever they rely, in part or in full, on manual photograph interpretation. Copyright © 2014 John Wiley & Sons, Ltd.     

安宁河断裂带勒帕沟泥石流发育特征及防治建议

2011年6月17日,勒帕沟发生近50年来最大的一次泥石流,经计算泥石流断面平均流速为2.03 m/s,峰值流量为200.11 m3/s,为大型泥石流,对帕勒村村民生命财产造成了巨大危害。勒帕沟上游岸坡前缘不稳定斜坡体及勒帕沟沟道内松散堆积的卵砾石为勒帕沟泥石流的主要物源。上游支沟中碎屑物块度较小,易被水流冲走,不会形成堰塞湖和溃决型泥石流。勒帕沟流域现今仍有再次暴发泥石流的可能,建议于流域中上游进行人工植树造林,在下游清理河道,进行监测预警,以实现泥石流综合防治。      

黄土高原丘陵沟壑区切沟侵蚀与地形关系分析——以纸坊沟流域为例

沟蚀是土壤侵蚀研究的主要内容之一,地形地貌是沟蚀的一个重要影响因子。本文以安塞纸坊沟流域作为研究区域,选择土地利用方式、土壤类型、坡度坡长因子、平面曲率、坡向和地形湿度指数6个因子,通过因子内切沟所占比重/整个研究区切沟所占比重计算各个因子的权重值,通过空间叠加分析土壤侵蚀敏感性,并通过重分类的方法把土壤侵蚀敏感性分为基本无侵蚀、轻度侵蚀、中度侵蚀、强度侵蚀、剧烈侵蚀5个等级,来研究切沟侵蚀与地形的关系。结果表明：切沟多发生在坡度坡长较大、地表湿度较高的林草地区域和更容易发生侵蚀的黄绵土区域,并且多分布在阴坡的凹面;对比分析切沟侵蚀和土壤侵蚀敏感性,切沟大多分布在中等侵蚀敏感性以上的区域,约占总切沟的90%;实验权重值对验证区冲沟的响应精度为82.43%（中度侵蚀及其以后阶段）,与实际值90.53%相差不大,说明此种方法对黄土丘陵沟壑区具有一定适用性,对黄土丘陵沟壑区水土保持工作有重要意义。