Gully erosion has been damaging archaeological sites in Grand Canyon during the last several decades, and there is a need to protect these features through mitigation, monitoring, and better geomorphic understanding. The purpose of this study was to assess the effectiveness of erosion-control structures , determine the accuracy and utility of aerial photogrammetry for monitoring gullies , and understand the geomorphology of the erosion . We performed total-station surveys and other data collection during February and October , 2002, at nine study sites in eastern and western Grand Canyon. Erosion-control structures are more prone to be damaged by flow when they are placed in reaches of very high local gradient. Treatments are generally successful in slowing erosion or causing deposition of sediment, but damaged erosion-control structures were shown to be less effective than intact structures, and actually increase local erosion in cases. Aerial photogrammetry was performed on four eroding archaeological sites in western Grand Canyon in March and October 2002 in order to assess the accuracy and change-detection utility of this tool. Accuracy was assessed on several different levels by comparing photogrammetry data to ground-survey data, and mean absolute vertical error ranged from 6-10 cm. Error of manual photogrammetry digital terrain models (DTMs) increased with topographic ruggedness and decreased with greater photogrammetric point density. Mean error reached a minimum of 5 cm for March and 6.5 cm for October when the ratio of point density to topographic ruggedness was ~40. Ground surveys and repeat photography indicated that two study gullies eroded or aggraded during the study period by 10-20 cm, but these changes were mostly undetected in the photogrammetry DTMs. Repeat ground surveys showed that gullies erode most at knickpoints and in steep reaches , and that new knickpoints tend to form in relatively steep reaches of a given channel. An area-slope erosion threshold was identified for the study sites and applied in a GIS-based model at four sites to show areas that exceed the threshold and are sensitive to gully erosion . Overall results show an upcatchment control of gully erosion and suggest that baselevel changes due to Glen Canyon Dam operation are subordinate controls.

No âmbito da geomorfologia pressupõe-se que existe uma relação entre declividade (S) e área de contribuição (A) que condiciona o escoamento superficial necessário para ocasionar o surgimento das erosões. A relação entre estes parâmetros pode ser expressa pelo índice S = aA-b, que representa o limiar crítico de desencadeamento das erosões. Implementando este índice no SIG, pode-se elaborar produtos cartográficos com a finalidade de reproduzir espacialmente as áreas críticas sujeitas a ocorrência das erosões. O presente trabalho teve por objetivo realizar um estudo sobre o desencadeamento das erosões lineares, utilizando técnicas de geoprocessamento para obtenção dos parâmetros A e S. A área em estudo compreende a Bacia do Córrego do Meio, localizada no município de São Pedro - SP, e apresenta uma área de 48,06 km². Os pontos de feições erosivas foram compilados de trabalhos realizados previamente na área, e também por meio do levantamento em imagens de satélites e ortofotos, totalizando 58 pontos de erosões. O modelo digital de elevação (MDE), mapa de área de contribuição e declividade foram elaborados em ambiente SIG, com diferentes pixels de saída, 2, 10 e 30 metros, e para processamento dos dados hidrológicos foram utilizados dois algoritmos: D8 e D-Infinity. Os parâmetros A e S foram calculados para cada uma das feições, analisadas de maneira integrada, e com base em características especificas como posição no relevo e profundidade dos canais. O tamanho do pixel de saída refletiu significativamente na obtenção dos parâmetros A e S, por isso, foram utilizados os modelos com resoluções maiores, pois corresponderam melhor às características da área. Em relação aos algoritmos de fluxo acumulado o D-Infinity apresentou dados mais condizentes com a topografia, e possibilitou o ajuste matemático do limiar crítico no gráfico S x A. O índice obtido a partir da análise integrada foi S = 0,060A-0,280, o índice para feições localizadas em encosta S = 0,0539A-0,412, e feições em drenagem S = 0,0748A-0,397. Além destes, foram obtidos outros índices específicos, determinados conforme a profundidade das feições (< 0,5; 0,5 a 1; 1 a 1,5 e > 1,5 m). A implementação dos índices no SIG possibilitou a elaboração das Cartas de Zoneamento de Áreas de Ocorrência de Erosão que apresentaram resultados consistentes, quando comparadas com outros produtos cartográficos, e podem orientar a ocupação adequada nestes locais, evitando a deflagração de novos processos. Apesar de ser um modelo de predição de erosão simples, os dados obtidos auxiliaram a compreensão dos mecanismos de desencadeamento, por isso, podem subsidiar as ações de planejamento ambiental, visando à recuperação das áreas comprometidas, com o intuito de minimizar os danos ambientais, e as perdas econômicas e sociais. / Considering the geomorphology aspects, it is assumed that there is a relationship between slope (S) and contributing area (A) which determines the runoff necessary for the onset of erosion. The relationship between the morphometric parameters can be expressed by the index S = aA-b which represents the topographic threshold for the onset of erosion. Implementing this index in the GIS, cartographic products can be performed with the purpose of reproducing spatially the critical areas subject to the erosion occurrence. The aim of this work is the study if the onset of gully erosions using geoprocessing techniques to obtain the parameters A and S. The study area is the watershed of the Córrego do Meio, located in São Pedro, São Paulo state, with an approximated area of 48.06 km². The erosion feature was compiled from previous study in the area and also through the survey on satellite images and orthophotographs, totalizing 58 erosion points. The digital elevation model (DEM), contribution area and slope were performed in GIS, with different output pixels, 2, 10 and 30 meters, and for the hydrological data processing, two algorithms were used: D8 and D-Infinity. The parameters A and S were calculated for each of the erosion feature, which were analyzed as a whole, and also based on the specific characteristics such as position in the landforms and depth of the channels. The size of the output pixel reflected significantly in obtaining the parameters A and S, therefore, models with higher resolutions were used because they better represent the characteristics of the area. In relation to the flow accumulation algorithms, D-Infinity presented more consistent data with the topography, and enabled the mathematical adjustment of the threshold in the S x A chart. The index obtained from the analysis of all erosion was S = 0.060A-0.280, the index for erosions located in hillside S = 0.0539A-0.412, and erosions in drainage S = 0.0748A-0.397. Besides these, other specific indices were determined, according to the depth of the channels (< 0.5, 0.5 to 1, 1 to 1.5 and > 1.5 m). The implementation of the indexes in the GIS made possible to perform the Zoning Map of Probable Erosion Areas, which showed consistent results when compared to other cartographic products and may direct suitable occupancy in these locations, so as to avoid the onset of new processes. Although it is a model for prediction of simple erosion, the obtained data helped to understand the mechanisms of the onset of the gully erosion, and therefore, can subsidize the actions of environmental planning, aiming at the recovery of the committed areas, with the aim of minimizing environmental damage and economic and social losses.

No âmbito da geomorfologia pressupõe-se que existe uma relação entre declividade (S) e área de contribuição (A) que condiciona o escoamento superficial necessário para ocasionar o surgimento das erosões. A relação entre estes parâmetros pode ser expressa pelo índice S = aA-b, que representa o limiar crítico de desencadeamento das erosões. Implementando este índice no SIG, pode-se elaborar produtos cartográficos com a finalidade de reproduzir espacialmente as áreas críticas sujeitas a ocorrência das erosões. O presente trabalho teve por objetivo realizar um estudo sobre o desencadeamento das erosões lineares, utilizando técnicas de geoprocessamento para obtenção dos parâmetros A e S. A área em estudo compreende a Bacia do Córrego do Meio, localizada no município de São Pedro - SP, e apresenta uma área de 48,06 km². Os pontos de feições erosivas foram compilados de trabalhos realizados previamente na área, e também por meio do levantamento em imagens de satélites e ortofotos, totalizando 58 pontos de erosões. O modelo digital de elevação (MDE), mapa de área de contribuição e declividade foram elaborados em ambiente SIG, com diferentes pixels de saída, 2, 10 e 30 metros, e para processamento dos dados hidrológicos foram utilizados dois algoritmos: D8 e D-Infinity. Os parâmetros A e S foram calculados para cada uma das feições, analisadas de maneira integrada, e com base em características especificas como posição no relevo e profundidade dos canais. O tamanho do pixel de saída refletiu significativamente na obtenção dos parâmetros A e S, por isso, foram utilizados os modelos com resoluções maiores, pois corresponderam melhor às características da área. Em relação aos algoritmos de fluxo acumulado o D-Infinity apresentou dados mais condizentes com a topografia, e possibilitou o ajuste matemático do limiar crítico no gráfico S x A. O índice obtido a partir da análise integrada foi S = 0,060A-0,280, o índice para feições localizadas em encosta S = 0,0539A-0,412, e feições em drenagem S = 0,0748A-0,397. Além destes, foram obtidos outros índices específicos, determinados conforme a profundidade das feições (< 0,5; 0,5 a 1; 1 a 1,5 e > 1,5 m). A implementação dos índices no SIG possibilitou a elaboração das Cartas de Zoneamento de Áreas de Ocorrência de Erosão que apresentaram resultados consistentes, quando comparadas com outros produtos cartográficos, e podem orientar a ocupação adequada nestes locais, evitando a deflagração de novos processos. Apesar de ser um modelo de predição de erosão simples, os dados obtidos auxiliaram a compreensão dos mecanismos de desencadeamento, por isso, podem subsidiar as ações de planejamento ambiental, visando à recuperação das áreas comprometidas, com o intuito de minimizar os danos ambientais, e as perdas econômicas e sociais. / Considering the geomorphology aspects, it is assumed that there is a relationship between slope (S) and contributing area (A) which determines the runoff necessary for the onset of erosion. The relationship between the morphometric parameters can be expressed by the index S = aA-b which represents the topographic threshold for the onset of erosion. Implementing this index in the GIS, cartographic products can be performed with the purpose of reproducing spatially the critical areas subject to the erosion occurrence. The aim of this work is the study if the onset of gully erosions using geoprocessing techniques to obtain the parameters A and S. The study area is the watershed of the Córrego do Meio, located in São Pedro, São Paulo state, with an approximated area of 48.06 km². The erosion feature was compiled from previous study in the area and also through the survey on satellite images and orthophotographs, totalizing 58 erosion points. The digital elevation model (DEM), contribution area and slope were performed in GIS, with different output pixels, 2, 10 and 30 meters, and for the hydrological data processing, two algorithms were used: D8 and D-Infinity. The parameters A and S were calculated for each of the erosion feature, which were analyzed as a whole, and also based on the specific characteristics such as position in the landforms and depth of the channels. The size of the output pixel reflected significantly in obtaining the parameters A and S, therefore, models with higher resolutions were used because they better represent the characteristics of the area. In relation to the flow accumulation algorithms, D-Infinity presented more consistent data with the topography, and enabled the mathematical adjustment of the threshold in the S x A chart. The index obtained from the analysis of all erosion was S = 0.060A-0.280, the index for erosions located in hillside S = 0.0539A-0.412, and erosions in drainage S = 0.0748A-0.397. Besides these, other specific indices were determined, according to the depth of the channels (< 0.5, 0.5 to 1, 1 to 1.5 and > 1.5 m). The implementation of the indexes in the GIS made possible to perform the Zoning Map of Probable Erosion Areas, which showed consistent results when compared to other cartographic products and may direct suitable occupancy in these locations, so as to avoid the onset of new processes. Although it is a model for prediction of simple erosion, the obtained data helped to understand the mechanisms of the onset of the gully erosion, and therefore, can subsidize the actions of environmental planning, aiming at the recovery of the committed areas, with the aim of minimizing environmental damage and economic and social losses.

[pt] A ocupação e uso dos solos têm claras implicações nos processos erosivos, porém devemos destacar também um dos atores principais para os processos erosivos, que é a evolução geomorfológica. A bacia do córrego das Palmeiras está inserida em um local com baixas amplitudes de relevo, mas que apresenta um alto potencial de vulnerabilidade a eventos de erosão e movimentos de massa. Destacamos obviamente o grande número de voçorocamentos que atingem determinadas áreas do médio vale do rio Paraíba do Sul, principalmente próximos ao seu eixo principal em um grande alinhamento de direção NE-SW a qual a bacia do córrego das Palmeiras está inserida. O presente estudo tem como objetivo principal descrever a origem e a distribuição dos processos erosivos na bacia do córrego das Palmeiras, na cidade de Barra do Piraí, Rio de Janeiro. Para que o principal objetivo desse trabalho seja alcançado, se torna necessário o cumprimento dos objetivos específicos, os quais são: Descrever a relação da orientação das voçorocas com as estruturas do substrato geológico; Descrever a relação entre processos erosivos e os níveis de base locais e; Descrever a relação entre processos erosivos e uso do solo. Para a realização deste trabalho foram utilizadas técnicas de análise e embasamentos teóricos da ciência geográfica e mais especificamente da área da geomorfologia. Podemos destacar como os principais resultados alcançados as análises referentes à relação entre os condicionantes geológicos/geomorfológicos e a ocorrência e comportamento de processos erosivos, permitindo afirmar que existe uma relação estreita entre a compartimentação geológica, orientação das estruturas, processos de rebaixamento de níveis de base, histórico de uso e ocupação do solo e a presença de processos erosivos. / [en] Land use has clear implications in the erosive processes and their relation with the geomorphological evolution can be seen notoriously. The Palmeiras stream basin is inserted in a place with low amplitudes of relief, but that presents a high potential of vulnerability to events of erosion and mass movements. We highlight obviously the large number of gully erosions that reach certain areas of the middle valley of the Paraíba do Sul river, mainly near its main axis in a great alignment of direction NE-SW to which the basin of the stream of Palmeiras is inserted. This study has as main objective the description of the origin and spatial distribution of the erosive processes in the Córrego dos Palmares basin. In order for the main objective of this work to be achieved, it is necessary to fulfill the specific objectives, which are: Analyze the spatial distribution of erosive processes; analyze the relation of the orientation of the gullies with the lithological structures of the geological substrate; Analyze the relationship between erosive processes and local base levels and; analyze the relationship between erosive processes and land use. For the accomplishment of this work, techniques of analysis and theoretical bases of the geographic science and more specifically of the area of geomorphology were used. We can highlight as the main results of the description the relationship between the geological / geomorphological conditions and the occurrence and behavior of erosive processes, allowing to affirm that there is a close relationship between the geological compartmentation, orientation of the structures. History of land use and occupation and the presence of erosive processes.

Gullies are significant contributors of sediment to streams in the southeastern USA. This study investigated gully erosion in the clay-rich soils of east Tennessee under a humid subtropical climate. The aims of this study were to (1) estimate long-term erosion rates for different gully geomorphic settings, (2) compare patterns of erosion for the different settings, and (3) model the response of gully erosion to freeze-thaw events. Erosion was measured weekly from June 2012 to August 2018 using 105 erosion pins distributed in gully channels, interfluves, and sidewalls. Erosion rates were estimated from average slopes of lines of best fit of pin lengths versus time. Maximum and minimum temperature was calculated daily using an on-site weather station and freeze-thaw events were identified. Gully erosion was modeled using antecedent freeze-thaw activity for the three geomorphic settings. Long-term erosion rates in channels, interfluves, and sidewalls were 2.5 mm/year, 20 mm/year, and 21 mm/year, respectively; however, week-by-week erosion was statistically different between the three settings, indicating different erosive drivers. Models of erosion with lagged freeze-thaw variables explained up to 34.8% of the variability in erosion variables; sidewall erosion was most highly related to freeze-thaw activity. Freeze-thaw in prior weeks was an important variable in all erosion models.

Gully erosion is a global problem that degrades land and reduces its utility for agriculture, development, and water quality. Quantification of sediment yield and control of sediment sources is essential for environmental protection. Five methods to evaluate erosion rates and sediment yield on an east Tennessee, USA, hillslope were compared: (1) physical measurement by removal of accumulated sediment using 10 L buckets; (2) repeated measurement of erosion pins in gully (erosional) and delta (depositional) areas; (3) geometric model using a combination trapezoidal prism-cylinder segment; (4) geometric model using a series of trapezoidal pyramids; and (5) 3D solid computer modeling. The 3D solid model created in SolidWorks was selected as the reference model and all other methods overestimated sediment yield to varying degrees. Erosion pin methods overestimated sediment yield by 368% in deltas and 123% in gullies. Volumetric measurement of sediment using buckets overestimated sediment yield by 160% due to void space in the buckets. The trapezoidal prism-cylinder segment model overestimated sediment yield by 66% and the trapezoidal pyramids method overestimated sediment yield by 5.7%. For estimation of sediment trapped behind an elliptical or circular silt fence dam, use of the trapezoidal pyramid method provides a good approximation comparable to 3D solid computer modeling.

Le bassin supérieur de l’oued Djemaa qui fait partie de la montagne du Djurdjura est affecté par une importante morphogenèse. Les mouvements de terrain présentent une grande variété dans le secteur d’étude et jouent un rôle prépondérant dans l’évolution des versants. L’étude géomophologique détaillée, une enquête auprès des populations locales, des analyses géotechniques et l’analyse sédimentologique des formations déplacées nous ont permis de faire une typologie des mouvements de terrain et de dégager les facteurs de leur développement et de leur localisation. Les fortes précipitations d’hiver et de printemps sont à l’origine de tous les mouvements de terrain et des ravinements, mais l’exposition et la lithologie permettent de comprendre les formes différentes rencontrées sur les adrets et sur les ubacs. Les activités humaines jouent un rôle croissant dans l’accélération de la morphogenèse. Dans cette étude, un système d’Information Géographique a été utilisé pour produire les cartes des risques de dégradation du milieu qui sont des outils pour la planification et la gestion de l’utilisation du sol dans le bassin supérieur de l’oued Djemaa. Il permet de croiser la vulnérabilité avec l’aléa. Le résultat de ces croisements que nous avons codifiés débouche sur une carte globale de niveaux de risque / The upper basin of the oued Djemaa which includes to the mountain of Djurdjura is affected by an important morphogenesis. The ground movements introduces a large variety in the study area and play a preponderant role in the evolution of mountainsides. The geomorphology detailed study an investigation into the local populations, geotechtonic and analyses it sedimentology of displaced trainings allows us to do ground typology movement and to chow theire development factors and their location. The big rainfalls of winter and sprinter are the origin of these ground movements and gully erosions, but exposition and lithology allow to understand different forms on northern and southern slope. Human activities play more and more an important role in morphogenesis acceleration. In this study, a geographical information system has been used to provide degradation maps of risks of which are tools for planning and management of soil using in the upper basin of the oued Djemaa. He allows to cross vulnerability witch vagary. The result of these crossroads which we codified results in a total map of risk levels

Considerando a importância de se definir as condições críticas para a deflagração, desenvolvimento e estabilização das voçorocas e outras feições erosivas lineares, e tendo em conta as poucas abordagens nesse sentido realizadas no Brasil, o objetivo geral desta pesquisa é avaliar os limites críticos topográficos para o desenvolvimento das formas resultantes da erosão linear em Piratininga (SP). Os procedimentos metodológicos envolveram a identificação e caracterização das feições erosivas por meio de fotografias aéreas e in situ; a mensuração da declividade (S) em campo e com diversas resoluções de Modelo Digital de Terreno (MDT) e da Área de Contribuição (A) também a partir do MDT; e, a definição dos limites críticos topográficos na forma da seguinte Função Potência: S = a × . A linha do limite crítico foi ajustada no limite inferior dos dados, paralelamente à inclinação linha ortogonal de regressão de S e A. As feições erosivas foram categorizadas em 4 classes e 13 subclasses de acordo com seus aspectos morfométricos e indicadores dos processos atuantes. A mensuração de S e A por diferentes métodos resultou em algumas variações significativas nos dados gerados, sobretudo na comparação entre os valores obtidos em campo e pelo MDT. Foi obtido o seguinte limite crítico topográfico: S = 0,064 × ,, que tem coeficientes semelhantes a outras pesquisas que analisaram feições erosivas resultantes do escoamento subsuperficial. A distinção entre as feições erosivas onde ainda há indícios de recuo das cabeceiras e feições erosivas que já atingiram o equilíbrio, assim como a análise dos pontos de início dos processos erosivos em relação ao limite crítico, forneceu bases empíricas para a discussão da evolução das feições erosivas lineares em relação ao limite crítico topográfico e também para ao entendimento dos estados estável e instável dos voçorocamentos, permitindo a construção de um modelo teórico de previsão da estabilidade e instabilidade das feições erosivas baseado na relação entre a declividade e a área de contribuição em diferentes perfis de encosta. / Considering the importance of defining the critical conditions for the generation, development and stabilization of gullies and other erosional features, and taking into account the few approaches accordingly held in Brazil, the aim of this research is to evaluate the topographical critical limits for developing of gully erosion in Piratininga (SP). The methodological procedures involved: the identification and characterization of the gullies using aerial photographs and in situ; measurement of the slope (S) in the field and with various resolutions of Digital Terrain Model (DTM) and drainage area (A) also from DTM, and the definition of topographical thresholds in the form of following Power Function: S = a × . The topographical threshold line has been set at the lower limit of data, parallel to the slope of orthogonal regression line of S and A. The erosional features were categorized into 4 classes and 13 subclasses according to their morphometric aspects and indicators of active processes. The measurement of S and A by different methods resulted in some significant changes in the data generated, particularly when comparing the values obtained in the field and the DTM. We obtained the following topographic threshold: S = 0,064 × , which has coefficients similar to other studies that analyzed erosional features resulting from subsurface flow. The distinction between active and stabilized gullies, as well as analysis of the starting points of erosion in relation to threshold line, provide an empirical basis for the discussion of the evolution of gullies in relation to the topographic threshold line and also for the understanding of stable and active states of gullies, allowing the construction of a theoretical model for predicting the stability and instability of the gullies based on the relationship between slope and drainage area in different slope profiles.

Gully erosion has been established as a major source of sediment pollution in the upper Trinity River watershed in north-central Texas. This fact, along with a lack of models appropriate for a large-area gully erosion analysis established a need for a gully erosion study in the upper Trinity basin. This thesis project attempted to address this need by deriving an index indicative of gully erosion risk using remote sensing and geographic information systems (GIS) methodology. In context of previous field studies, the coarse spatial resolution of the input GIS data layers presented a challenge to prediction of gully prone areas. However, the remote sensing/GIS approach was found to provide useful reconnaissance information on gully risk over large areas.

O presente trabalho propõe índices que indicam o limiar topográfico para predição do desencadeamento de canais de erosão linear em uma região de clima tropical. A área escolhida para a realização do estudo localiza-se no município de São Pedro/SP e apresenta uma área de aproximadamente 7 \'KM POT.2\'. Nesta área foi realizado um inventário detalhado das feições erosivas pré-existentes para obtenção dos valores das áreas das bacias de contribuição, considerando as cabeceiras de algumas erosões e as declividades médias dessas áreas, por meio de ferramentas de SIG e também por fotografias aéreas ortorretificadas. Trabalhos de campo, com o auxílio de GPS de precisão, permitiram definir a situação atual das erosões e com o uso de fotografias aéreas de anos anteriores foi possível observar a sua evolução. Com isso, através dos mapas em SIG, como MDT, declividade, direção de fluxo e etc. foram obtidos os valores de área de contribuição (A) e declividade média (S). Os valores de A x S foram apresentados em gráficos que permitiram identificar o ponto mínimo para a ocorrência de um evento erosivo e estabelecer de forma empírica os limiares topográficos que definem as condições para o início do processo erosivo. Para o Índice I foram feitas duas equações (uma para canais menos profundos que 1 metro e outra para canais mais profundos que 1 metro) e para o Índice II foi obtida apenas uma equação, pois em fotografias aéreas nas escalas disponíveis (1:25.000 e 1:30.000) foram possíveis serem observados apenas os canais profundos. Os índices obtidos neste trabalho foram para o Índice I: S = 0,0219A(-0,376) (canais profundos) e S = 0,0128A(-0,443) (canais rasos); enquanto que para o Índice II o valor obtido é: S = 0,0588A(-0,377). Para a validação destes índices propostos, foram realizados os mesmos procedimentos para obtenção das áreas de contribuição (A) e declividade média (S), porém em outra área. A área escolhida para se fazer a validação dos índices foi a Bacia do Ribeirão do Meio, que apresenta 48 \'KM POT.2\' aproximados de área e apresenta características físicas semelhantes da área utilizada neste presente estudo, ou seja, apresenta semelhança em geologia, geomorfologia, clima, uso e ocupação. Foi observado com esta validação, que todos os pontos com problemas de erosão da área usada para validação localizam-se dentro da área prevista pelos índices para ocorrência de problemas de erosão, o que significa dizer que os índices exprimem o limiar topográfico para erosões lineares nas condições encontradas na área de estudo e em áreas semelhantes. / This paper proposes indexes indicating the topographic threshold for predicting the trigger of gully erosion in a tropical climate region. The area chosen for the study is located in São Pedro/SP and presents an area of approximately 7 \'KM POT.2\' . In this area was carried out a detailed inventory of pre-existing erosion to obtain the values of watershed areas of contribution, considering the headwaters of some erosions and the mean slope of these areas, using GIS tools and also orthorectified aerial photographs. Fieldwork, with the assistance of GPS accuracy, allowed defining the current condition of erosion and with the use of aerial photographs from previous years was possible to observe its evolution. Thus, through the maps in GIS, as DEM, Slope, Flow Direct ion among others, were obtained values from contribution area (A) and mean slope (S). The values of A x S were represented in charts al lowing to identify the lower point for the occurrence of an erosion event and establish empirical ly the topographic thresholds that def ine the conditions for initiating the erosion process. For the index I were made two equations (one for channels less than 1 meter deep and one for channels deeperthan 1 meter) and Index II was obtained only one equation, because in aerial photographs available (scales1:25,000 and 1:30,000) were possible only be observed the deep channels. The indexes obtained in this paper were for the Index I: S = 0,0219A(-0,376) (deep channels) and S = 0,0128A(-0,443) (shallow channels), while the index II the value obtained is: S = 0,0588A(-0,377). To validate these proposed indexes, the same procedures were performed to obtain the contribution area (A) and mean slope (S), but in another area. The area chosen to perform the validation of the indexes was the Basin Ribeirão do Meio, which is 48 \'KM POT.2\' extents and shows similarity in geology, geomorphology, climate, use and occupation. The results have been showed that the selected points with erosion problem of the area used for validation are located wi thin the area specified by the indexes for the occurrence of erosion problems, which means that the indexes reflect the topographic threshold for gully erosions in conditions found in the study area and similar areas.