基于DEM的数字流域特征提取研究进展

数字高程模型(DEM)是地表形态高程属性的数字化表达,利用流域DEM数据构建数字水系模型并提取流域水文特征,是分布式水文过程模拟的重要基础。DEM中闭合洼地和平坦区域的处理、流向的确定以及DEM分辨率的大小都直接影响着流域水系特征提取的质量与效率。针对当前基于DEM提取河网与流域特征的诸多问题,阐述了DEM数据提取流域水系特征的原理,回顾了数字水系模型与流域特征提取方法,评述了洼地和平地的处理方法及水流方向的确定方法的研究进展,介绍了当前基于不同DEM数据类型的提取方法研究,探讨了DEM尺度和分辨率对提取流域特征的影响,总结了平缓区域数字水系和河网提取的研究进展。     

地形指数的物理意义分析

地形指数模型 ( TOPMODEL)用地形指数在流域中的空间格局来确定流域饱和缺水量的空间分布和产流区的空间位置与范围 ,物理意义明确。本文介绍了地形指数的物理基础及其与土壤水分的关系。分析了地形指数空间变化与汇流面积 a及局地坡度空间变化的关系 ,a范围内有效汇流面积的变化 ,提出了一种确定 a上限值的方法。     

Investigating the surface process response to fault interaction and linkage using a numerical modelling approach

In order to better understand the evolution of rift‐related topography and sedimentation, we present the results of a numerical modelling study in which elevation changes generated by extensional fault propagation, interaction and linkage are used to drive a landscape evolution model. Drainage network development, landsliding and sediment accumulation in response to faulting are calculated using CASCADE, a numerical model developed by Braun and Sambridge, and the results are compared with field examples. We first show theoretically how the ‘fluvial length scale’, L_f, in the fluvial incision algorithm can be related to the erodibility of the substrate and can be varied to mimic a range of river behaviour between detachment‐limited (DL) and transport‐limited (TL) end‐member models for river incision. We also present new hydraulic geometry data from an extensional setting which show that channel width does not scale with drainage area where a channel incises through an area of active footwall uplift. We include this information in the coupled model, initially for a single value of L_f, and use it to demonstrate how fault interaction controls the location of the main drainage divide and thus the size of the footwall catchments that develop along an evolving basin‐bounding normal fault. We show how erosion by landsliding and fluvial incision varies as the footwall area grows and quantify the volume, source area, and timing of sediment input to the hanging‐wall basin through time. We also demonstrate how fault growth imposes a geometrical control on the scaling of river discharge with downstream distance within the footwall catchments, thus influencing the incision rate of rivers that drain into the hanging‐wall basin. Whether these rivers continue to flow into the basin after the basin‐bounding fault becomes fully linked strongly depends on the value of L_f. We show that such rivers are more likely to maintain their course if they are close to the TL end member (small L_f); as a river becomes progressively more under supplied, i.e. the DL end member (large L_f), it is more likely to be deflected or dammed by the growing fault. These model results are compared quantitatively with real drainage networks from mainland Greece, the Italian Apennines and eastern California. Finally, we infer the calibre of sediments entering the hanging‐wall basin by integrating measurements of erosion rate across the growing footwall with the variation in surface processes in space and time. Combining this information with the observed structural control of sediment entry points into individual hanging‐wall depocentres we develop a greater understanding of facies changes associated with the rift‐initiation to rift‐climax transition previously recognised in syn‐rift stratigraphy.     

Catchment-scale (dis)connectivity in sediment flux in the upper Hunter catchment, New South Wales, Australia

(Dis)connectivity within and between landscape compartments affects the extent and rate of transfer of energy and matter through catchments. Various landforms may impede sediment conveyance in a river system, whether laterally to the channel (termed buffers) or longitudinally along the channel itself (termed barriers). A generic approach to analysis of landscape (dis)connectivity using slope threshold analysis in GIS, tied to air photograph interpretation and field mapping of buffers and barriers, is tested in the upper Hunter catchment, Australia. Under simulated conditions, effective catchment area, which is a measure of the proportion of a catchment that has the potential to contribute sediment to the channel network, varies from 73% to just 3% of the total catchment area for differing subcatchments in the upper Hunter. This variability can be explained by the spatial distribution and assemblage of buffers and barriers in each subcatchment. Multiple forms of disconnectivity are evident in some subcatchments, such that when one buffer or barrier is breached, other features still impede sediment transfer within the system. The importance of the position of buffers and barriers within any given subcatchment is emphasised. Spatial variability in valley width exerts a critical control on catchment connectivity, with more efficient sediment conveyance in narrow valleys relative to wider valleys characterised by piedmonts, terraces, fans and extensive floodplains in which conveyance is impeded. This variability reflects the landscape history and geological setting of each subcatchment. The framework developed in this paper can be used to assess the impact of natural or human-induced buffers and barriers on catchment-scale sediment flux in any landscape setting, providing a physical template atop which other biogeochemical fluxes could be examined.     

大尺度分布式水文模型数字流域提取方法研究

构建大尺度分布式水文模型是当前大气水文模型耦合研究的一项重要内容。本文介绍一 种根据1km DEM 生成更大网格尺度DEM 数据, 同时可以保持流域河网信息并减缓高程、坡度 等地貌参数信息量衰减速度的有效方法———ZB算法。利用该方法和常规的网格平均法生成黄河 唐乃亥以上流域的5km、10km、15km 和20km 两套DEM 数据, 分别提取高程、坡度、地形指数、河 网密度、主河道长度、流域面积等流域特征参数, 并与1km DEM 提取的上述参数进行比较, 对两 种方法作出评价。结果显示, 随着网格尺度的增大, ZB 算法获得的DEM 数据可以保持河网的连 续性, 提取出合理的流域范围, 减缓地形信息量的衰减速度。该方法满足构建大尺度分布式水文 模型提取数字流域的需要。     

The Fractal Characteristics of Drainage Networks and Erosion Evolution Stages of Ten Kongduis in the Upper Reaches of the Yellow River,China

The fractal characteristics of drainage in the ten kongduis of the upper Yellow River were obtained using the box counting dimension, and the evolution stages of the watershed topography were defined by different ranges of the fractal dimensions of river networks (Dg). The results show that the fractal scaleless range of the Maobula River is 20–370 m based on a combination of artificial judgment, correlation coefficient test and fitting error. Other kongduis show good fractal characteristics in this fractal scaleless range as well. The box counting dimension can be used as a quantitative index of watershed topography fractal characteristics. The fractal dimension of stream networks is independent of the threshold contributing area used for extracting the drainage networks from the DEM. The values of Dg in the upper ten kongduis are in the range of 1.08?1.14. Both the runoff yield and the sediment yield are positively and linearly related with Dg. The positive relation between the sediment yield and Dg reflects the effect of landform features on sediment yield in the young and/or mature stages of landform evolution of the study area. By revising the critical value of Dg, the value of Dg of the basin in the young evolution stage is less than 1.06, while it is more than 1.06 for the basin in mature or old evolution stage. The upper ten kongduis are in the mature stage of landform evolution.     

基于栅格数字高程模型提取特征地貌技术研究

本文对近年来基于栅格数字高程模型提取特征地貌技术进行了详细的研究,认为该技术的关键在于两个方面：一是如何定义地貌形态结构,二是提取算法的设计。本文提出了基于地貌学角度来定义地貌形态结构的方法,利用有限个数的形态要素的空间组合和对比分析来获取特征地貌,并可以对各种特征地貌形态进行符合物理意义的改进。     

Adjustment of a drainage network to capture induced base-level change: an example from the Sorbas Basin, SE Spain

Quaternary catchments in the south of the Sorbas Basin, SE Spain have been affected by two regionally significant river captures. The river captures were triggered by changes in regional gradients associated with sustained Quaternary uplift in the region of 160 m Ma⁻¹. The first capture occurred in the early Pleistocene and re-routed 15% of the original Sorbas Basin drainage into the Carboneras Basin to the south. The second occurred in the late Pleistocene and re-routed 73% of the original Sorbas Basin drainage to the east. This latter capture had dramatic consequences for base-level in the Sorbas Basin master drainage. Local base-level was lowered by 90 m at the capture site, 50 m at 7 km upstream and 25 m at 13 km upstream of the site. The base-level change instigated a complex re-organisation of the drainage networks in systems tributary to the master drainage over the ensuing period (some 100 ka). After the capture, drainage systems closer to the capture site experienced a tenfold increase in incision rates over most of their network. Those located some 13 km upstream of the capture site experienced a fivefold increase in incision, although in this instance, the changes do not appear to have propagated to the headwater regions of the drainage nets. The sensitivity of individual catchments was largely governed by geological controls (structure and lithology). The detailed network evolution in the most sensitive areas can be traced by reconstructing former drainage pathways using abandoned drainage cols and the alignment and degree of incision of the drainage networks. Three main stages of evolution can be identified which record the progressive spread of base-level changes from the master drainage. These are Stage 1 (pre-capture): original south-to-north consequent drainage; Stage 2 (early stage, post capture): aggressive subsequent southwest-to-northeast and east–west drainage developed along structural lineaments first in the east of the area (Stage 2a), and later in the west of the area (Stage 2b); and Stage 3 (late stage, post capture): obsequent drainage developing on the topography of the Stage 2 drainage. All stages of the network evolution are associated with drainage re-routing as a function of river capture at a variety of scales. The results highlight the complex response of the fluvial system, and the very different geomorphological histories of adjacent catchments, emphasising the need for regional approaches for examining long-term changes in fluvial systems.     

基于不同DEM数据源的数字河网提取对比分析——以韩江流域为例

基于HYDRO1K、SRTM3和ASTER GDEM三种DEM数据,利用BTOPMC地形子模型提取韩江流域河网,并作对比分析。结果表明:①SRTM3提取的河网精度最高,HYDRO1K相对最低。②DEM的垂直精度对提取的河网精度起控制作用。ASTER GDEM的水平分辨率较高,但垂直精度不如SRTM3,因而提取的河网精度不如SRTM3。③HYDRO1K提取大尺度流域河网具有一定的精度,但在地势平坦区域的效果较差,HY-DRO1K不宜用来提取小尺度流域河网。④由DEM提取的数字河网精度与当地的地面坡度以及处理DEM的填洼算法有关。     

不同时空尺度下沟壑对流域侵蚀产沙的影响——以黄土丘陵沟壑区岔巴沟流域为例

基于DEM,分别以不同的沟壑发育所需的汇水面积为阈值,提取岔巴沟流域等6个不同空间尺度流域的沟壑.并求出不同级别沟壑的差值.结合黄河子洲观测站多年的水沙观测资料,分析不同级别的沟壑及沟壑的不同部位在不同时间尺度下对流域侵蚀产沙的影响.分析结果表明:沟壑密度并不是简单地随着流域面积的增大而减小.沟壑密度的空间尺度效应与所提取沟壑的级别有关;黄土丘陵沟壑区梁峁坡、沟坡的次暴雨产沙量占流域产沙总量的比重与降雨特性有关;重力侵蚀多发生在毛沟的中、上部及沟头处,而支沟、干沟主干道的沟坡比较稳定,且其发生的重力侵蚀对流域产沙的影响多为"无效"部分;沟壑发育汇水面积的阈值为1500、1800、2000个栅格提取的沟壑对流域侵蚀产沙的作用较大,同时沟壑对流域侵蚀产沙的影响与时间尺度有关,在多年平均尺度下,高级别沟壑只起到泥沙输送通道的作用.     

Recent fold growth and drainage development: The Janauri and Chandigarh anticlines in the Siwalik foothills, northwest India

The active growth of a fault-and-thrust belt in frontal zones of Himalaya is a prominent topographical feature, extending 2500 km from Assam to Pakistan. In this paper, kinematical analysis of frontal anticlines and spatial mapping of active faults based on geomorphological features such as drainage pattern development, fault scarps and uplifted Quaternary alluvial fans are presented. We analyse the geomorphic and hydrographic expressions of the Chandigarh and the Janauri active anticlines in the NW India Siwaliks. To investigate the morphological scenario during the folding process, we used spatial imagery, geomorphometric parameters extracted from digital elevation models and fieldwork. Folding between the Beas and Sutlej Rivers gives clear geomorphological evidence of recent fold growth, presumably driven by movements of blind thrust faults. Structural style within the Janauri and Chandigarh anticlines is highly variable (fault-propagation folds, pop-up structures and transfer faults). The approach presented here involves analysis of topography and drainage incision of selected landforms to detect growth of active anticlines and transfer faults. Landforms that indicate active folding above a southwest-dipping frontal thrust and a northeast-dipping back-thrust are described. Along-strike differences in ridge morphology are measured to describe the interaction of river channel patterns with folds and thrust faults and to define history of anticline growth. The evolution of the apparently continuous Janauri ridge has occurred by the coalescence of independent segments growing towards each other. By contrast, systematic drainage basin asymmetry shows that the Chandigarh anticline ridge has propagated laterally from NW to SE.     

Characteristics, controlling factors and importance of deep gullies under cropland on loess-derived soils

Several studies describe the formation and importance of shallow ephemeral gullies under cropland (depth<0.8 m). Some of these gullies may develop into channels with a depth of more than 0.8 m (up to 4 m). Despite their spectacular nature, these deep gullies have not yet been studied in detail. Therefore, the objectives of this study are to analyze the characteristics and the controlling factors of these deep gullies, as well as their importance in terms of sediment production.Comparison of a dataset with 28 deep gullies, formed in the period 1985–2003 and 123 shallow ephemeral gullies formed in the period 1994–1999 in central Belgium indicates that the deep gullies have clearly different morphological characteristics compared to shallow ephemeral gullies. Several factors were analyzed to understand the formation of deep gullies. Plotting runoff contributing area versus slope of the soil surface at the gully head indicates that the topographical threshold for deep gully formation is significantly larger compared to ephemeral gully formation. Deep gullies form on short, steep valley sides and their position is strongly affected by the presence of linear landscape elements. All deep gullies incised on landscape positions with a very erodible soil horizon at shallow depth. Analysis of causative rainfall showed no significant differences in incision thresholds between rills and shallow ephemeral gullies on the one hand and deep gullies on the other hand.The relation between area-specific sediment yield in central Belgium and drainage area, indicates that the development of deep gullies contributes significantly to the sediment output of small rural catchments and causes peaks in the mean area-specific sediment yield that are up to a factor of three higher compared with catchments where shallow ephemeral gullying occurs.     

Role of landscape parameters in riverine run-off,and sediment and organic matter yield on Disko Island,West Greenland

Measurements of discharge and suspended sediment and organic matter yield from nine different drainage basins on the island Disko in central West Greenland were carried out in the period 19–30 July 1997. A series of landscape properties (glacier cover, altitude, slope, aspect and vegetation cover) were measured for each drainage basin from a digital elevation model and a satellite image. Principal Component Analysis and regression statistics have been used to examine which landscape parameters dominate the measured discharges and yields. It is concluded that differences in suspended sediment and organic matter yield between the drainage basins can be explained by the measured morphometric properties while differences in water run-off can not. Glacier cover percentage and percent area above 800 m a.s.l. were the most important parameters influencing suspended sediment and organic matter yield in the drainage basins on Disko during the study period.     

Relationships between morphological and sedimentological parameters in source‐to‐sink systems: a basis for predicting semi‐quantitative characteristics in subsurface systems

The study of source‐to‐sink systems relates long‐term variations in sediment flux to morphogenic evolution of erosional–depositional systems. These variations are caused by an intricate combination of autogenic and allogenic forcing mechanisms that operate on multiple time scales – from individual transport events to large‐scale filling of basins. In order to achieve a better understanding of how these mechanisms influence morphological characteristics on different scales, 29 submodern source‐to‐sink systems have been investigated. The study is based on measurements of morphological parameters from catchments, shelves and slopes derived from a ∼1 km global digital elevation model dataset, in combination with data on basin floor fans, sediment supply, water discharge and deposition rates derived from published literature. By comparing various morphological and sedimentological parameters within and between individual systems, a number of relationships governing system evolution and behaviour are identified. The results suggest that the amount of low‐gradient floodplain area and river channel gradient are good indicators for catchment storage potential. Catchment area and river channel length is also related to shelf area and shelf width, respectively. Similarly to the floodplain area, these parameters are important for long‐term storage of sediment on the shelf platform. Additionally, the basin floor fan area is correlative to the long‐term deposition rate and the slope length. The slope length thus proves to be a useful parameter linking proximal and distal segments in source‐to‐sink systems. The relationships observed in this study provide insight into segment scale development of source‐to‐sink systems, and an understanding of these relationships in modern systems may result in improved knowledge on internal and external development of source‐to‐sink systems over geological time scales. They also allow for the development of a set of semi‐quantitative guidelines that can be used to predict similar relationships in other systems where data from individual system segments are missing or lacking.     

Scaling regimes of local slope versus contributing area in digital elevation models

Four scaling regimes are observed in Digital Elevation Models (DEMs) when the average local slope is calculated for pixels grouped according to the values of the contributing area. Threshold criteria, proposed by various researchers to identify the extent of the channel network, are examined relative to the slope-area scaling diagram. The scaling response observed in DEMs is reproduced with a landscape evolution and channel network growth model originally developed by Willgoose et al. (1989). The threshold criterion proposed in this model is a useful tool to separate two different slope-area scaling regimes observed in DEM data.     

Automatically derived catchment boundaries and channel networks and their hydrological applications

A method of generating drainage networks and catchment boundaries from gridded ground elevation data (a digital terrain model) is presented. The method includes an algorithm for coping with localised low points in the data. The paper suggests two objective ways of comparing the conventionally mapped with the automatically generated drainage network. Emphasis is placed on the way that this approach is capable of benefiting a variety of hydrological applications, in particular the automatic estimation of flow statistics from catchment characteristics.     

Channel network morphology and sediment dynamics under alternating periglacial and temperate regimes: a numerical simulation study

The occurrence of permafrost in a highly permeable catchment has a profound effect on runoff generation. The presence of permafrost effectively makes the subsoil impermeable. Therefore, overland flow can be the dominant runoff-generating process during periglacial conditions. The absence of permafrost will promote subsurface drainage and, therefore, saturation excess overland flow can become the dominant runoff-generating process during temperate conditions. In this paper, we present a numerical modelling study in which the effect of alternating climate-related phases of permafrost and nonpermafrost on catchment hydrology and geomorphology is investigated. Special attention is given to the characteristics of the channel network being formed, and the sediment yield from these catchments. We find that channel networks expand under permafrost conditions and contract under nonpermafrost conditions. A change from permafrost to nonpermafrost conditions is characterised by a decrease in sediment yield, while a change towards permafrost conditions is marked by a peak in sediment yield. This peak is explained by the build-up of a reservoir of erodible sediment during the nonpermafrost phase. The driving force behind this reservoir build-up may be local base-level change due to tectonic uplift or eustacy. We present a number of experiments, which show the details of this process. The results are in line with existing reconstructions of climate and fluvial dynamics during the Pleistocene in Europe and offer a new explanation to these observations.     

Geomorphic constraints on the active tectonics of southern Taiwan

Taiwan is a region of rapid active tectonics, yet the study of the tectonic processes that shape the interior of the island is difficult due to the high rates of erosion and dense vegetation. We use digital topography to look for indications of active deformation preserved in the local geomorphology. In particular, anomalies in the regional pattern of drainage are used to infer zones of enhanced tectonic activity. The apparent anticlockwise rotation of major river systems in plan view indicates the presence of a diffuse zone of left-lateral shear running down the southeastern side of Taiwan. Asymmetries in the catchments of individual drainage basins show the influence of varying rates of uplift across southern Taiwan, with the most rapid uplift close to Taitung at the indentation point of the Luzon arc with the Chinese continental margin. Our interpretations, though based predominantly on remote-sensing observations, are consistent with the available field evidence. This study demonstrates the usefulness of drainage basins as tectonic markers in the quantification of regional strain and uplift, which may have wider applicability in other deforming parts of the world.     

基于DEM的复杂地形流域特征提取

利用流域数字高程模型(DEM)构建数字水系模型并提取流域水系特征是分布式水文过程模拟的重要前提。提出了面向分布式水文过程模拟和流域特征提取的数字水系模型,并针对现有方法对复杂地形DEM中含有的平地、洼地及其嵌套情形的处理不足,提出了栅格水流分类、填洼分类与归并及有效填平处理、河谷平地的出流代价法构建栅格流向和流序等新的处理方法,并在开发的软件系统得到实现。使用该方法创建的黄土岭流域数字水系模型和提取的水系等流域特征结果表明:本文方法可有效应对复杂地形流域的处理,提取的流域水系特征与实际自然水系比较吻合,能够有效地消除现有方法在地形平坦区域容易产生的平行河道、奇异河道、河道变形等不足。     

Characterising the origin,nature and fate of sediment exported from catchments perturbed by active tectonics

Changes to the tectonic boundary conditions governing erosional dynamics in upland catchments have a significant effect on the nature and magnitude of sediment supply to neighbouring basins. While these links have been explored in detail by numerical models of landscape evolution, there has been relatively little work to quantify the timing, characteristics and locus of sediment release from upland catchments in response to changing tectonic boundary conditions that are well‐constrained independently. We address this challenge by quantifying the volume and granulometric characteristics of sediment exported from modern rivers draining across active normal faults in the Central Apennines in Italy. We demonstrate that catchments undergoing a transient response to tectonics are associated with significant volumetric export of material derived primarily from the zone upstream of the fault, producing bi‐modal grain‐size distributions with elevated D₈₄ values within the transient reach. This is in direct contrast to the headwaters, where the fluvial capacity to transport sediment is low and the grain‐size distribution of material in transit is fine and uni‐modal. The grain‐size response is driven by landslides feeding coarse material directly into the channel, and we show the amplitude of the signal is modulated by the degree of tectonic perturbation, once the threshold for bedrock landsliding is exceeded. Additionally, we evaluate the length‐scale over which this transient grain‐size signal propagates downstream into the basin. We show that the coarse‐fraction sediment released is retained in the proximal hanging‐wall if rates of tectonic subsidence are high and if the axial river system is small or far from the fault‐bounded mountain front. Our results therefore provide some of the first quantitative data to evaluate how transient landscape responses affect the locus, magnitude and calibre of sediment supply to basins.