基于DEM的黄土地貌逼近度因子构建及特征分析

沟谷源点作为沟沿线上最为活跃的部位,其分别到上游分水线、下游沟谷线的流线空间比对关系,是表征三线空间结构的重要突破口。它在空间上向流域分水线逼近的程度,是量化黄土流域地貌系统发育程度的重要切入点。为探究黄土地貌区流域沟谷源点向流域分水线逼近的程度,揭示黄土流域地貌发育进程及该进程所表现的主要侵蚀方式,本文从水平和垂直2个维度,构建量化三线空间结构关系的核心因子--逼近度(PI),其中包括水平逼近度(HPI)和垂直逼近度(VPI),基于5 m分辨率数字高程模型,在陕北黄土高原遴选了包含16种地貌类型的42个样区,利用数字地形分析方法,探讨其平均值(MHPI, MVPI)的空间分异规律。另外,选择其中南北序列分别代表黄土塬、残塬、梁状丘陵沟壑和峁状丘陵沟壑区的淳化、宜君、甘泉和绥德4个地区,完备包含一至五级别沟谷的典型流域为重点实验样区,探讨陕北黄土高原流域尺度平均逼近度变异指数(MPIV)序列性特征。实验结果表明:①陕北黄土高原平均逼近度(MPI)存在强烈的空间自相关,MHPI在南北序列上先增大后减小,东西方向上逐渐减小,在黄河沿岸峡谷丘陵区MHPI达到最大;MVPI由西南-东北先减小后增加,由西北-东南逐渐减小,在渭北黄土台塬区达到最小;②在流域尺度上,MPIV值的正负与黄土塬区、丘陵沟壑区敏感相关;③4个重点样区MHPI、MVPI与其它地形因子存在南北序列上的一致性。104个外部汇流区平均水平逼近度与平均坡度相关性较好(P=0.43, a<0.001),平均垂直逼近度与面积高程积分强烈相关(P=0.75, a<0.001)。平均逼近度指标综合考量了黄土高原地区最典型的三条具有结构控制意义特征线的空间关系,对黄土地貌的发育程度有明显的指示性意义。     

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

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

Hierarchy structure characteristics analysis for the China Loess watersheds based on gully node calibration

A land surface region can be decomposed into a series of watershed units with a hierarchical organizational structure. For loess landform, the watershed is a basic spatial–structural unit that can express natural landforms, surface morphology characteristics, spatial organization and developmental evolution. In this research we adopted the concept of node calibration in the watershed structure unit, selected six complete watersheds on China Loess Plateau as the research areas to study the quantitative characteristics of the hierarchical structure in terms of watershed geomorphology based on digital elevation model(DEM) data, and then built a watershed hierarchical structure model that relies on gully structure feature points. We calculated the quantitative indices, such as elevation, flow accumulation and hypsometric integral and found there are remarkably closer linear correlation between flow accumulation and elevation with increasing gully order, and the same variation tendency of hypsometric integral also presented. The results showed that the characteristics of spatial structure become more stable, and the intensity of spatial aggregation gradually enhances with increasing gully order. In summary, from the view of gully node calibration, the China Loess watershed structure shows more significantly complex, and the developmental situation variation of the loess landforms also exhibited a fairly stable status with gully order increasing. So, the loess watershed structure and its changes constructed the complex system of the loess landform, and it has the great significance for studying the spatial pattern and evolution law of the watershed geomorphology.     

复杂网络理论在黄土高原沟谷地貌特征研究中的应用

沟谷侵蚀是塑造黄土地表侵蚀形态的主要动力,沟谷的发育过程深刻地影响着黄土地貌的发育及演化。本文在黄土高原选择6个典型地貌样区,以样区的数字高程模型为基本数据源提取沟谷系统。将沟谷系统中的沟谷节点、沟谷源点和流域出水口点作为网络节点,网络节点之间的空间拓扑关系为边,高程差为权值,构建黄土高原沟谷加权复杂网络模型。对黄土沟谷地貌的节点特征和空间结构进行定量刻画和分析,得到黄土高原不同地貌类型网络特征的空间格局及其变化,并进一步映射地貌的发育过程及演化机理。研究结果表明：① 黄土高原沟谷加权网络的点强度累积概率分布呈指数分布,相关性系数皆达到0.80以上,该网络正处于向无标度网络转化的过渡期;② 样区从南到北,沟谷特征点的点强度值呈现逐渐减小的态势,且点强度的分布具有不对称性,沟谷右侧侵蚀强度较高,点强度分布较多;③ 平均路径长度和网络结构熵值在绥德一带最大,分别为30.94、6.31,并向南北两侧呈减少的趋势,网络密度值的变化与之相反;3个指标分别从网络结构的连通性、稳定性和紧密性反映了不同沟谷地貌类型的侵蚀程度以及地貌系统的演化机理;④ 网络指标与传统地貌指标的相关性系数均超过0.85,其可以科学、准确地表达地貌形态的复杂性及地貌的发育阶段,有望作为沟谷地貌地学特征研究的参数。该方法考虑了沟谷地貌的空间拓扑关系以及系统的整体性,为复杂表面形态的地貌研究提出了一种新的思路和方法。     

Clustering gully profiles for investigating the spatial variation in landform formation on the Chinese Loess Plateau

The gully is the most dynamic and changeable landform unit on the Loess Plateau, and the characteristics of gully landforms are key indicators of gully evolution. Different gully profiles are connected and combined through runoff nodes. Thus, it is necessary to cluster gully profiles into a gully profile combination(GPC) to reveal the spatial variation in gully landforms throughout the Loess Plateau. First, the gradient and gully evolution index(GEI) of two sample areas in Changwu and Suide in Shaanxi Province, China are calculated and analysed based on GPC. Then, the gradient and GEI are calculated by using 90-m-resolution digital elevation model(DEM) data for the severe soil erosion area with the basin as the research unit. On this basis, the spatial variation in the development degree is analysed with Getis-Ord Gi*. The results show that the degree of gully undercutting decreases from southeast to northwest under the influence of rainfall. Due to the soil properties, the loess in the northwest is more prone to collapse, resulting in the decrease of GEI from northwest to southeast. The development degree of gullies is closely related to rivers. The strong erosive capacity of rivers leads to greater differences in gullies within the basin. At the same time, the skewness and kurtosis of the gully index in the basin are correlated; when the distribution of the gully index in the basin is less normal, the distribution of the gully index is more concentrated. These results reveal the spatial variation characteristics of the Loess Plateau based on GPC.     

Extraction and analysis of gully head of Loess Plateau in China based on digital elevation model

In China′s Loess Plateau area, gully head is the most active zone of a drainage system in gully areas. The differentiation of loess gully head follows geospatial patterns and reflects the process of the loess landform development and evolution of its drainage system to some extent. In this study, the geomorphic meaning, basic characteristics, morphological structure and the basic types of loess gully heads were systematically analysed. Then, the loess gully head′s conceptual model was established, and an extraction method based on Digital Elevation Model(DEM) for loess gully head features and elements was proposed. Through analysing the achieved statistics of loess gully head features, loess gully heads have apparently similar and different characteristics depending on the different loess landforms where they are found. The loess head characteristics reflect their growth period and evolution tendency to a certain degree, and they indirectly represent evolutionary mechanisms. In addition, the loess gully developmental stages and the evolutionary processes can be deduced by using loess gully head characteristics. This study is of great significance for development and improvement of the theoretical system for describing loess gully landforms.     

Quantifying spatial scale of positive and negative terrains pattern at watershed-scale: Case in soil and water conservation region on Loess Plateau

The positive and negative terrains (P-N terrains) widely distributed across China’s Loess Plateau constitute the dual structure characteristic of loess landforms. Analysis of loess P-N terrains at the watershed scale can serve to elucidate the structural characteristics and spatial patterns of P-N terrains, which benefits a better understanding of watershed evolution and suitable scales for loess landform research. The Two-Term Local Quadrat Variance Analysis (TTLQV) is calculated as the average of the square of the difference between the block totals of all possible adjacent pairs of block size, which can be used to detect both the scale and the intensity of landscape patches (e.g., plant/animal communities and gully networks). In this study, we determined the latitudinal and longitudinal spatial scale of P-N terrain patterns within 104 uniformly distributed watersheds in our target soil and water conservation region. The results showed that TTLQV is very effective for examining the scale of P-N terrain patterns. There were apparently three types of P-N terrain pattern in latitudinal direction (i.e., Loess Tableland type, Loess Hill type, and Transitional Form between Sand and Loess type), whereas there were both lower and higher values for P-N terrain pattern scales in all loess landforms in the longitudinal direction. The P-N terrain pattern also clearly presented anisotropy, suggesting that gully networks in the main direction were well-developed while others were relatively undeveloped. In addition, the relationships between the first scales and controlling factors (i.e., gully density, nibble degree, watershed area, mean watershed slope, NDVI, precipitation, loess thickness, and loess landforms) revealed that the first scales are primarily controlled by watershed area and loess landforms. This may indicate that the current spatial pattern of P-N terrains is characterized by internal force. In selecting suitable study areas in China’ Loess Plateau, it is crucial to understand four control variables: the spatial scale of the P-N terrain pattern, the watershed area, the main direction of the watershed, and the loess landforms.     

DEM分辨率对黄土侵蚀沟形态特征表达的不确定性分析

黄土侵蚀沟的地形表达是开展黄土沟谷侵蚀研究的基础工作,利用数字高程模型(DEM)定量描述侵蚀沟特征有助于研究侵蚀沟的形态变化和发育过程。基于DEM数据计算多种指标对黄土侵蚀沟特征进行描述是目前侵蚀沟研究中最为常用的方法。但是,受到格网DEM数据结构的限制,其计算结果会存在一定的不确定性。在侵蚀沟地形表达时,对形态特征的表达会受到DEM数据分辨率的影响,进而造成表达结果的不确定性。尤其在黄土高原地区,地形特征更为破碎,地形要素更为复杂,其表达结果受DEM分辨率的影响更为明显。本文以黄土高原典型样区为例,基于点云数据建立不同分辨率的DEM数据集,通过不同地形因子对侵蚀沟特征进行表达,分析DEM分辨率在黄土侵蚀沟形态特征表达时的不确定性。结果显示,分辨率的降低对主沟支沟比和纵比降等侵蚀沟形态特征因子产生了较大影响,且指标与分辨率多呈现线性变化关系。但是,随着侵蚀沟的横向扩张,DEM分辨率对其特征表达的影响逐渐被削弱。此外,在使用固定分析窗口进行侵蚀沟特征计算时,由于分辨率的降低,格网尺寸增大,其实际分析半径随之增大,使得计算范围内地表形态变化增加,导致沟谷切割深度随着分辨率的降低反而增加。同时,侵蚀沟主沟道区域受分辨率影响较小,沟头区域指标与分辨率的关系较弱。     

利用数字摄影测量方法估算半干旱区小流域沟谷侵蚀产沙

土壤侵蚀不仅表现黄土高原强烈的现代地貌过程,也是其生态环境恶化的象征,其地貌演化和土壤侵蚀互为因果,现代研究表明黄土高原半干旱小流域的沟谷侵蚀产沙占其水土流失总量的重要部分,近年有关沟谷侵蚀的定量研究成果,加深了人们对沟道侵蚀危害的认识,为土地资源的合理利用、沟谷危害的评价等提供了科学依据。传统的沟谷侵蚀定量研究是采用实地量测的方法,现代的数字摄影测量技术的发展为估算沟谷侵蚀研究提供了一个经济高效方法。本研究采用的髙分辨率数字高程模型是由3期历史航空照片根据数字摄影测量方法制作而成,利用它对黄土高原半干旱区一个小流域的1959～1981年及1981～1999年两个时段内沟谷侵蚀产沙量进行估算,将流域侵蚀产沙的实测值与估算结果进行评价。并指出此方法在沟谷侵蚀产沙量估算方面的统计学意义及在未来应用时要注意的问题;根据流域侵蚀产沙结果分析历史时期人类活动的影响:人类活动在一定程度上改善了环境,绝大部分泥沙被拦截在流域里,同时沟道侵蚀加剧,如何采取措施减少流域坡面和沟道的侵蚀,是目前人们所面对的问题;此外,结合流域水沙统计数据的DEM的侵蚀空间分析,结果表明坡面流、重力侵蚀及沟谷的下切是导致流域土壤侵蚀的主要原因。     

坡谱信息熵尺度效应及空间分异

选取陕北黄土高原48个不同地貌类型区,以其对应的1∶10000及1∶50000数字高程模型(DEM)为数据源,运用比较分析与数理统计的方法,分析了坡谱信息熵在陕北黄土高原的空间分异特征及尺度效应。结果显示,坡谱的信息熵在一定程度上可以反映地表的复杂度,坡谱信息熵和沟壑密度间有较好的幂函数关系,随着坡谱信息熵的增大,样区沟壑密度也增大。坡谱信息熵的空间分异和陕北黄土高原的黄土地貌形态在空间上的变异是相关的,可将其作为地貌类型划分的判别因子之一。     

Morphology and controlling factors of the longitudinal profile of gullies in the Yuanmou dry-hot valley

The morphology of the gully longitudinal profile (GLP) is an important topographic index of the gully bottom associated with the evolution of the gullies. This index can be used to predict the development trend and evaluate the eroded volumes and soil losses by gullying. To depict the morphology of GLP and understand its controlling factors, the Global Positioning System Real-time Kinematic (GPS RTK) and the total station were used to measure the detail points along the gully bottom of 122 gullies at six sites of the Yuanmou dry-hot Valley. Then, nine parameters including length (Lt), horizontal distance (Dh), height (H), vertical erosional area (A), vertical curvature (Cv), concavity (Ca), average gradient (Ga), gully length-gradient index (GL), normalized gully length-gradient index (Ngl), were calculated and mapped using CASS, Excel and SPSS. The results showed that this study area is dominated by slightly concave and medium gradient GLPs, and the lithology of most gullies is sandstone and siltstone. Although different types of GLPs appear at different sites, all parameters present a positively skewed distribution. There are relatively strong correlations between several parameters: namely Lt and H, Dh and H, Lt and A, Dh and A, H and GL. Most GLPs, except three, have a best fit of exponential functions with quasistraight shapes. Soil properties, vegetation coverage, piping erosion and topography are important factors to affect the GLP morphology. This study provides useful insight into the knowledge of GLP morphology and its influential factors that are of critical importance to prevent and control gully erosion.     

Experimental Study of the Interaction between Building Clusters and Flash Floods

In recent years, flash flood disasters have occurred frequently in southwest China due to the increased frequency of extreme climate events. To solve this problem, great efforts have been made in studying the process of flash flood. However, little attention was paid on bearing body of hazard, the clusters of buildings. Thus the real disaster mechanism of flash flood remains unclear.Accordingly, based on the experiments of artificial flash floods in a conceptual solid model, this paper focuses on the flood-impacted inundation characteristics of the building clusters at different locations of the gully model, in order to obtain a better understanding of the disaster process and the interaction between the flash floods and building clusters. The results showed that, in a typical smallscale flash flood gully with hot and dry climate, 1)clusters of buildings on an alluvial fan could reduce about 35% of the flooding area by blocking the diffusion of the flood to the depression areas, and could also promote the deposition in lower reaches of the river channel by blocking the overbank flow from going back into the channel, making the width-depth ratio of the channel larger. 2) The flash flood rates of disaster and hazard on the alluvial fan are generally higher than that of the inner gully. For the inner gully,buildings located on the beaches along the lower river and the transitional areas of the straight channel and channel bends can easily be affected because of their lower elevations. For the alluvial fan, buildings nearby the meanders suffer the greatest impacts because of bank collapsing and flooding. 3) The safe vertical distance from a building to the river channel is 13 m for the buildings in the inner gully under extreme floods. Below this threshold, the smaller the vertical distance is, the greater the risk exposure is. For the buildings on the alluvial fan, especially for the buildings near the concave bank of the top rush point,the horizontal distance is more important, and the safe value is 80 m under extreme floods.     

Snake Model for the Extraction of Loess Shoulder-line from DEMs

Shoulder lines are the most important landform demarcations for geographical analysis,soil erosion modeling and land use planning in the Loess Plateau area of China.This paper proposes an automatic,effective and accurate method of determining loess shoulder line from DEMs by integrating a hydrological D8 algorithm and a snake model.The watershed boundary line is adopted as the initial contour which evolves to identify the exact position of loess shoulder-line by the guidance of an external force of snake model from DEMs.Experiments show that the method overcomes the difficulties in both threshold selection for edge detection and the disconnecting issues in former extraction approaches.The accuracy evaluation of shoulder-line maps from the two test sites of the loess plateau area show obvious improvements in the extraction.The average contour matching distance of the new method is 12.0 m on 5 m resolution DEM,and shows improvement in the accuracy and continuity.The comparisons of accuracy evaluations of the two test sites show that the snake model method performs better in the loess plain area than in the area with high gully density.     

Multitemporal UAV-based photogrammetry for landslide detection and monitoring in a large area:a case study in the Heifangtai terrace in the Loess Plateau of China

With high spatial resolution,on-demand-flying ability,and the capacity for obtaining threedimensional measurements,unmanned aerial vehicle(UAV) photogrammetry is widely used for detailed investigations of single landslides,but its effectiveness for landslide detection and monitoring in a large area needs to be investigated.The Heifangtai terrace in the Loess Plateau of China is a loess terrace that is extremely susceptible to irrigation-induced loess landslides.This paper used UAV-based photogrammetry for a series of highresolution images spanning over 30 months for landslide detection and monitoring of the terrace with an area of 32 km~2.Dense and evenly distributed ground control points were established and measured to ensure the high accuracy of the photogrammetry results.The structure-from-motion(Sf M) technique was used to convert overlapping images into orthographic images,3D point clouds,digital surface models(DSMs) and mesh models.Using multitemporal differential mesh models,landslide vertical movements and potential landslides were detected and monitored.The results indicate that a combination of UAV-based orthophotos and differential mesh models can be used for flexible and accurate detection and monitoring of potential loess landslides in a large area.     

基于DEM的黄土崾岘提取及其地形特征分析

崾岘是将要被切穿的鞍部,是正负地形矛盾斗争的结果,也是重要的地形控制点。典型的崾岘多位于黄土高原黄土地貌区,又称黄土崾岘,其对识别沟间地与沟谷的斗争程度有一定的指示作用。本文以黄土高原样区为例,基于1:1万DEM（5 m分辨率）和影像分辨率为0.95 m的遥感影像,利用流域边界算法和缓冲区标定,分析窗口选择5×5,实现了崾岘点位的半自动化提取。并对各崾岘点位求取坡度等地形因子,总结崾岘的空间格局和地形特征。结果显示,崾岘多分布在主流域边界和垂直于主沟道的最宽部分,地形控制作用明显。崾岘的坡度、起伏度、切割深度等值均大于鞍部值,同时,高级流域区的崾岘值大于低级流域区的崾岘值,反映出崾岘具有侵蚀程度强、表层完整性低、地表破碎度高的特点。总体而言,崾岘受沟道蚕食度高,从侧面反映了黄土地貌的发育阶段,是黄土地貌发育到中期的标志性产物。     

黄土高原丘坡信息DEM提取算法的应用

选择代表黄土地貌类型(丘陵沟壑区、梁峁区、高原丘陵过渡区)的75个样区(每个样区约4km²),以1∶1万水平分辨率为5m的DEM为研究对象,研究不同算法对提取地面坡度精度的影响。采用6种不同算法分别提取坡度。运用方差分析,对比分析,排序分析等方法,通过比较不同算法所提取坡度矩阵的平均值、最大值、标准差,以及坡度中误差等指标,认为在黄土丘陵沟壑区3类地貌类型区域如果整体考虑坡度的以上指标时选用三阶反距离平方权差分算法较为合理。在黄土高原沟壑区与丘陵沟壑区交错过渡地带也可以选用算法4提取坡度。同时,根据水土保持工作实际生产部门需要,提出并实现了对坡度信息评价时须分级别进行的思路。研究认为,在黄土丘陵沟壑区域相同的地貌类型区,算法对所获得的相同的坡度级别面积影响是不同的,并且规律不相同。在前述三种地貌类型区域,如果需要获得坡度分级面积信息,使用三阶反距离平方权差分算法和三阶反距离权差分最为理想。由于简单差分算法小于其他算法提取坡度所涉及的区域单元大小,不适宜在本区域使用。     

Gully erosion regionalization of black soil area in northeastern China

Gully erosion is the frequent and main form of soil erosion in the black soil area of the northeastern China, which is one of the most important commodity grain production bases in China. It is encroaching upon the fertile farmland there. Regionalization of gully erosion can reveal the spatial distribution and regularity of the development of gully erosion. Based on the eco-geographical regional background features of the black soil area, this study combined the regionalization with influencing factors of the development of gully erosion. GIS spatial analysis, geostatistical analysis, spatial statistics, reclassification, debris polygon processing and map algebra methods were employed. As a result, the black soil area was divided into 12 subregions. The field survey data on type, length, volume and other characteristics indicators of gully erosion were used to calibrate the results. Then the features of every subregion, such as where the gully erosion is, how serious it is, and why it happens and develops, were expounded. The result is not only an essential prerequisite for gully erosion surveys and monitoring, but also an important basis for gully erosion prevention.     

Digital extraction of altitudinal belt spectra in the West Kunlun Mountains using SPOT-VGT NDVI and SRTM DEM

In this paper, a digital identification method for the extraction of altitudinal belt spectra of montane natural belts is presented. Acquiring the sequential spectra of digital altitudinal belts in mountains at an acceptable temporal frequency and over a large area requires extensive time and work if traditional methods of field investigation are to be used. Such being the case, often the altitudinal belts of a whole mountain or the belts at a regional scale are represented by single points. However, single points obviously cannot accurately reflect the spatial variety of altitudinal belts. In this context, a digital method was developed to extract the spectra of altitudinal belts from remote sensing data and SRTM DEM in the West Kunlun Mountains. By means of the 1km resolution SPOT-4 vegetation 10-day composite NDVI, the horizontal distribution of altitudinal belts were extracted through supervised classification, with a total classification accuracy of 72.23%. Then, a way of twice-scan was used to realize the automatic transition of horizontal maps to vertical belts. The classification results of remote-sensing data could thus be transformed automatically to sequential spectra of digital altitudinal belts. The upper and lower lines of the altitudinal belts were then extracted by vertical scanning of the belts. Relationships between the altitudinal belts based on the montane natural zones concerning vegetation types and the geomorphological altitudinal belts were also discussed. As a tentative method, the digital extraction method presented here is effective at digitally identifying altitudinal belts, and could be helpful in rapid information extraction over large-scale areas.     

IMPACT EVALUATION OF HAIZUKA DAM ON ITS UP STREAM： A CASE STUDY IN HIROSHIMA PREFECTURE, JAPAN

Japan ranks fifth in the world for the number of large dams. Environmental impacts of large dams are known, such as enormous losses of water or disruption of fish spawning, however, impacts of the dams on their up streams are functions of topography of the up stream. Haizuka Dam is located in Hiroshima Prefecture, Japan and its implementation will start in 2006. This large dam influences its up stream through dam making activities, which occurs in order and with different spatial presences that were categorized into chronological and spatial impacts. In this case study, spatial impacts were further divided into horizontal and vertical ones. The horizontal impacts were identified as new roads, diversion tunnel, dam lake, and submerged cultivated land, while vertical impacts were recognized as submerged historical monuments, slope protections, dam body, and deforested area in the reservoir. There were convergences of spatial and temporal impacts, however, the extent of the impacts was limited to the lake boundary.     

Geomorphic characteristics of hillslope and channelized debris flows: A case study in the Shitou area of central Taiwan

The data on the hillslope and channelized debris flows in the Shitou area of central Taiwan occurred during Typhoons Toraji and Nali in 2001 were applied in this paper. The geomorphic parameters, including the flow length, gully gradient, drainage area and form factor of the debris flows were determined by spatial analysis using a Geographic Information System （GIS） based on the data derived from field investigation, aerial photographs, and topographical maps. According to such determined geomorphic parameters, the threshold conditions and empirical equations, such as the relationship between the gully gradient and drainage area and that between gully length and drainage area and topographic parameter, are presented and used to distinguish the geomorphic characteristics between the channelized and hillslope debris flows.