Parâmetros morfométricos como indicadores de evolução de bacias de drenagem

Autores

Milagres, A. (UFMG) ; Augustin, C. (UFMG) ; Fonseca, B. (UFMG)

Resumo

A análise dos canais fluviais a partir da utilização de parâmetros morfométricos tem se mostrado eficiente para responder a várias questões relativas à evolução de bacias de drenagem. O objetivo deste trabalho é apresentar análise dos perfis longitudinais de rios estudados nas bacias do rio Pardo e Jequitinhonha, no Espinhaço Meridional, MG, como base para a compreensão da evolução e dissecação do relevo. Para isto, foram analisados os resultados do Índice de Hack, da Concentração da Rugosidade, e as Densidades de knickpoints de canais por ordem hierárquica. Foram selecionadas áreas das bacias do Pardo e alto Jequitinhonha localizadas exclusivamente sobre rochas siliciclásticas do Supergrupo Espinhaço, portanto limitando a análise ao mesmo domínio litológico. Os resultados evidenciam padrões de evolução diferenciados dentro do mesmo domínio litológico, possivelmente associados a mudanças nos fatores ligados aos níveis de base e/ou tectônicos na região.

Palavras chaves

Parâmetros morfométricos; perfis longitudinais; Supergrupo Espinhaço

Introdução

Rios possuem grande importância no estudo do modelado da paisagem, uma vez que são os mais ativos agentes na esculturação do relevo (LEOPOLD e MADDOCK, 1953). Nesse sentido o estudo de perfis longitudinais de rios pode ser de grande valia para a compreensão dos processos evolutivos, e das análises morfoestruturais e morfotectônicas. Através da análise de perfis longitudinais de rios, Hack (1973) concluiu que é possível avaliar o estado de equilíbrio de determinado trecho de rio. A partir desses estudos, propôs o Stream Gradient Index, mais conhecido como Índice de Hack (IH), que é utilizado para detectar alterações nos cursos fluviais decorrente de anomalias de drenagem, e avaliar o vigor energético dos mesmos. O valor energético, por sua vez, responderia pelo trabalho executado pelo rio e assim, pelo processo de erosão e transporte envolvidos na elaboração do relevo da bacia. Hack (1975) também considera que em uma topografia ajustada, qualquer mudança na declividade do canal encontra-se relacionada ao relevo total, de tal maneira que a ocorrência de um desequilíbrio em um canal, acarretaria um reajuste de todo o sistema fluvial da bacia. Esse desequilíbrio seria passível de ser identificado nos perfis longitudinais dos níveis de base locais, ou knickpoints. De acordo com Bjornberg (1969, apud GUEDES et al., 2006), estes ocorrem principalmente devido a: 1) confluência com trechos de caudal forte; 2) variação da resistência litológica; 3) variação eustática, o que causa erosão remontante por modificação forte no nível de base; e 4) por deformações tectônicas ou neotectônicas, representando mudanças na rede hidrográfica que são refletidas em toda a morfologia da bacia, na busca de equilíbrio, por isto, dito dinâmico. A análise da posição e distribuição de knickpoints revela-se essencial para a interpretação da incisão da rede fluvial e da evolução das paisagens marcadamente dissecadas pela erosão fluvial (PHILLIPS et al., 2010). Assim, a distribuição desses pontos nas bacias hidrográficas constitui uma tarefa fundamental nos estudos que abordam a evolução da paisagem (SCHUMM et al., 2002; CROSBY et al., 2006). A análise dos perfis longitudinais em uma área relativamente homogênea do ponto de vista litológico pode evidenciar mais claramente condicionantes estruturais, e elucidar a sua relação com a energia dos rios e com o modelado do relevo. O objetivo deste trabalho é o de apresentar os resultados de estudo dos perfis longitudinais de rios de duas bacias, no sentido de auxiliar na compreensão da dissecação e evolução do relevo da área estudada. Parâmetros morfométricos foram utilizados para a análise dos canais de cursos fluviais, como maneira de subsidiar interpretações com relação a energia presente nas bacias, o que é indicativo do grau de equilíbrio dinâmico das mesmas. A pesquisa foi conduzida nas bacias do rio Pardo e do alto Jequitinhonha, em domínios exclusivos das rochas siliciclásticas do Supergrupo Espinhaço, na Serra do Espinhaço Meridional, que se estende abaixo do paralelo 17ºS até as proximidades de Belo Horizonte, MG. De acordo com Fogaça (1997), o Supergrupo Espinhaço encontra-se estruturado por um conjunto de falhas de empurrão/zonas de cisalhamento com direção geral N-S e por amplas dobras de eixo que tendem para N-S, resultado das compressões de direção E-W do Evento Brasiliano, (650-550 Ma). Este fato explica a ocorrência de maior deformação na porção leste da Serra, onde está localizada a bacia do alto Jequitinhonha (FONSECA e AUGUSTIN, 2014), a qual, possui um relevo muito dissecado, característico por seus desníveis acentuados e alta declividade das vertentes, já a bacia do Pardo, mais preservada dessas deformações, possui um relevo menos dissecado. A análise destas bacias do ponto de vista energético dá suporte para compreender os padrões de evolução diferenciados dentro do mesmo domínio litológico do Supergrupo Espinhaço.

Material e métodos

Este trabalho se baseou em dados de radar SRTM, com resolução espacial de 30 metros, obtidos do Projeto TOPODATA (VALERIANO et al., 2009) no qual foram elaborados os Modelos Digitais de Elevação (MDE). A hidrografia foi retirada do Instituto Mineiro de Gestão das Águas – IGAM (2012). Os parâmetros morfométricos aplicados foram: A hierarquização fluvial, obtida pelo método de Strahler (1952). O mapa de Densidade de knickpoints, foi produzido a partir da ferramenta knickpoint Finder (SALAMUNI et al., 2013; QUEIROZ et al., 2015), que trata-se de um software que funciona de maneira acoplada ao ArcGIS®. Segundo Salamuni et al. (2013), esta ferramenta foi criada com base no índice de Hack, e busca, de forma automática, pontos que caracterizam rupturas no terreno a partir de imagens raster que possuam dados altimétricos, de onde a drenagem tridimensional é obtida para a análise do IH. Para isso foi utilizado dados topográficos do MDE. O resultado foi uma malha de pontos que indicam os knickpoints. A esta malha foi aplicado o estimador de densidade de Kernel, considerando um raio de 2000 metros, para a espacialização das zonas com maior adensamento de knickpoints, gerando o mapa. O Índice de Hack (IH), que foi proposto por Hack (1973) trata-se da relação entre a declividade e a extensão do canal fluvial em perfis longitudinais de rios e consiste na seguinte fórmula: SL=ΔH/ΔL•L sendo SL o índice (Declividade-Extensão), ∆H a variação altimétrica do segmento, ∆L a projeção horizontal do trecho a ser estudado e L, o comprimento do ponto mais a jusante do trecho estudado à cabeceira da bacia em linha reta (HACK, 1973). Ele pode também ser calculado pela equação simplificada: SL=ΔH/log_e•ΔL Foram realizados cálculos do IH para toda a extensão de um curso fluvial, considerando-se a amplitude altimétrica total entre a cota da cabeceira e a cota da foz, bem como o logaritmo natural da extensão total do mesmo, conforme sugerido por (ETCHEBEHERE et al., 2004; FONSECA, 2010; FONSECA e AUGUSTIN, 2011; 2014). Os valores calculados para este índice foram plotados na metade do traçado de cada drenagem (ETCHEBEHERE et al., 2004 e FONSECA, 2010), de maneira a possibilitar a aplicação do interpolador IDW (Inverse distance weighted) para espacializar as informações, os dados foram tratados no programa ArcGIS 10.1. O Índice de Concentração da Rugosidade (ICR), proposto por Sampaio (2008) e melhor discutido por Sampaio e Augustin (2014) adota como referencial de análise, a distribuição espacial da declividade, utilizando-se de valores indiretos da rugosidade, calculados a partir da repetição dos valores de declividade por unidade de área (SOUZA e SAMPAIO, 2010). Esta análise espacial também foi realizada utilizando o estimador de densidade de Kernel, em ambiente ArcGis 10.1. Foi utilizada ainda a Densidade de Drenagem (Dd), proposta por Horton (1945), que correlaciona o comprimento total dos canais de escoamento com a área da bacia hidrográfica. O seu cálculo pode ser realizado conforme a seguinte fórmula: Dd=L_t/A onde: Dd é a densidade da drenagem, L_t, o comprimento total dos canais e A é a área da bacia. Foi produzida uma Dd para cada ordem hierárquica da área. Todos os processos foram realizados no arcgis 10.1, e espacializados pelo estimador de densidade de Kernel. Foram confeccionados gráficos no arcgis 10.1, através da ferramenta de Geostatistical analyst, utilizando-se duas variáveis (Densidade de knickpoints e Densidade de Canais de primeira ordem). Essas variáveis foram normalizadas utilizando o shape da área total de estudo. No momento de produção dos gráficos, foi realizada a divisão desta área total em duas bacias (Pardo e alto Jequitinhonha) e os gráficos foram produzidos para cada bacia, de maneira a possibilitar o entendimento das diferenças de densidade entre as duas.

Resultado e discussão

Na Figura 1 são apresentados três mapas que auxiliam na compreensão da evolução geomorfológica da área, baseados em índices espacializados que permitem avaliar as mudanças de energia fluvial dentro de uma área homogênea do ponto de vista litológico. Dessa maneira, é possível dar destaque a outros fatores capazes de “desequilibrar” um sistema fluvial, como exemplo, a mudança de nível de base e reativações tectônicas. As feições herdadas do evento Brasiliano (600 – 550 Ma) influenciam diretamente os processos modeladores das formas de relevo na área estudada, sendo que a bacia do alto Jequitinhonha foi a que sofreu maiores deformações, sendo caracterizada por intensos falhamentos e grandes cavalgamentos (KNAUER, 1990; UHLEIN, 1991; DUSSIN e DUSSIN, 1995; RENGER e KNAUER, 1995, KNAUER, 2007; CHEMALE et al., 2011). Isto é evidenciado pela densidade de knickpoints (Figura 1A), maior nessa bacia, sugerindo que estejam associados a estas zonas de fraqueza herdadas do Brasiliano, com posteriores reativações tectônicas ou não. A bacia do Pardo, quando analisada em toda a extensão, indica a presença de poucos knickpoints, apontando uma baixa ocorrência de rupturas do relevo, indicando tratar-se de área com grandes compartimentos aplainados, circundados por afloramentos, como descrito e mapeados por Augustin et al. (2011) e estudados por Augustin et al. (2014). A maior concentração de knickpoints ocorre na sua borda oeste (Figura 1A), possivelmente por causa da influência do contato litológico dos quartzitos com os carbonatos do Grupo Bambuí, muito bem marcado por grandes escarpas. A ocorrência de maior densidade de knickpoints no alto Jequitinhonha corrobora os resultados de Fonseca e Augustin (2014), no sentido de que nesta porção da bacia há maior energia e, consequentemente, uma intensidade maior de processos tanto de encostas, quanto fluviais, favorecendo maior dissecação do relevo, que é bem distinta entre as duas bacias. O mapa do IH (Figura 1B), que mostra a espacialização dessa energia, permite visualizar a delimitação de grandes compartimentos de valores de alto e baixo índice, por grandes áreas de influência energética. É possível observar que os valores altos correspondentes aos locais de maior energia no sistema e que estes apontam para expressiva diferença dos valores entre a porção leste, ocupada pela bacia do alto Jequitinhonha, e a oeste, pela bacia do Pardo. Na porção oeste, pode-se observar a presença de grandes áreas de baixa energia, classificadas por Augustin et al. (2011) como “áreas erosivas remanescentes”, caracterizadas por extensos platôs pouco recortados por intrusões de rochas metabásicas, o que parece responder por sua evolução mais homogênea, no qual se destacam os afloramentos de quartzito, intercalados por deposições de materiais arenosos. Já a bacia do Jequitinhonha, que apresenta altos valores energéticos, Augustin et al. (2011) classificaram como “áreas modeladas pela incisão vertical da rede de drenagem em planos de falhas e fraturas”, que formam profundas gargantas em quartzitos, como exemplo, o rio Jequitinhonha, o ribeirão do Inferno e córrego Rapadura. Essa característica, apontada pelos autores acima citados, é apoiada também pelos resultados constantes do mapa do ICR (Figura 1C), no qual se observa maior grau de dissecação na bacia do Jequitinhonha É evidente a correspondência espacial entre a densidade de knickpoints e o ICR, o que possibilita confirmar que a densidade de knickpoints possui relevante significado geomorfológico, ao indicar processos de retomada erosiva e aumento na intensidade de dissecação do relevo, sendo de grande importância para os estudos de evolução da paisagem. Outra forma de analisar a maior dissecação na bacia do Jequitinhonha é analisando a densidade da drenagem da área. Na figura 2, é possível verificar uma proporcionalidade relativa da entre as duas bacias, no que se refere ao número de rios de hierarquias superiores por área. Entretanto, quando se trata dos cursos de ordens hierárquicas 2, e principalmente 1 (nascentes), a bacia do Jequitinhonha possui uma proporção bem maior dessas drenagens, principalmente na sua porção mais leste. Isto indica maior energia das bacias, pois valores elevados da Dd são indicativos indiretos do grau também mais elevado de dissecação do relevo pela erosão fluvial, decorrente de fatores como declividade, presença de material subsuperficial pouco permeável, baixa densidade da cobertura vegetal, fatores climáticos, entre outros (FONSECA e AUGUSTIN, 2014) A área leste da bacia do Jequitinhonha com alta Dd apresenta fortes indícios de retomada erosiva recente, o que pode ser interpretado, segundo Fonseca e Augustin (2011) por sua maior proximidade com o mar (nível de base geral), uma vez que a bacia do rio Pardo tem como seu nível de base a bacia do rio São Francisco. Outra possibilidade, de acordo com Augustin (1995) e Fonseca e Augustin (2011; 2014), é que a borda oeste do Supergrupo Espinhaço se comporte como um possível nível de base local para a bacia do Pardo, não permitindo a saída de material que se acumula nas grandes áreas aplainadas dessa bacia. Essas feições são recorrentes em toda essa bacia, áreas aplainadas circundadas por afloramentos. A maior concentração de knickpoints, rugosidade, energia fluvial, e a maior densidade de canais de ordem 1, como pode ser observado nas figuras 1 e 2, são encontrados na bacia do alto Jequitinhonha. Contudo, ao analisar os gráficos de densidade de Knickpoints dos canais de ordem 1, observa se uma proporção maior de knickpoints nos canais de ordem 1 da bacia do Pardo (Figuras 3 e 4). De acordo com Crosby e Whipple (2006), os knickpoints estão associados à dinâmica de ajustamento de canais fluviais, que acompanham alterações climáticas e tectônicas. Os knickpoints são vistos como níveis de base locais, que por quebrarem a energia do fluxo fluvial, funcionam como barreiras para o carreamento dos sedimentos, promovendo a deposição dos mesmos à montante, regulando o entalhe do vale a montante (PHILLIPS et al., 2010) e promovendo uma retomada erosiva. Dessa forma, entende-se que o equilíbrio de uma bacia começa por seus rios principais e segue em cadeia até seus afluentes, em direção às nascentes (erosão à remontante). A resposta do sistema a um impulso de incisão pode ser realizado por etapas. Na primeira, os canais principais ao se ajustarem ao sistema, criam knickpoints, que propagam energia a montante; na segunda etapa, quando este ajustamento fluvial se propaga até chegar aos rios de ordem inferior, eles são incapazes de acompanhar o ritmo da descida e, por isso, desenvolvem knickpoints muito íngremes, que precisam de um tempo bem maior para se ajustarem ao sistema. Esses sucessivos ajustes parecem ter acontecido na bacia do rio Pardo, que já alcançou um maior equilíbrio, ainda que hoje suas nascentes ainda estejam parcialmente se ajustando. Parte desse reajuste ainda está presente em especial nas áreas de maior ruptura de declive, correspondentes aos contatos litológicos entre as rochas do Supergrupo Espinhaço (siliciclásticas) e carbonáticas do grupo Bambuí. Já a bacia do alto Jequitinhonha, em decorrência da sua maior energia, ainda busca um equilíbrio, o que se reflete na maior densidade de knickpoints dos seus rios principais.

Considerações Finais

O estudo de bacias de drenagem com o uso de parâmetros morfométricos como Densidade de knickpoints, Índice de Hack, Densidade de drenagem e Índice de Concentração da Rugosidade mostrou que estes são eficientes indicativos do estado de equilíbrio e ajustamento das mesmas. Sua espacialização foi essencial para a interpretação da ocorrência e distribuição dos fenômenos, o que possibilita compreender melhor as suas influências no modelado do relevo. É possível concluir através dos parâmetros morfométricos, que as bacias estão em processos de equilíbrio distintos, no qual a bacia do alto Jequitinhonha possui maior vigor energético, e a bacia do Pardo está mais ajustada a um perfil de equilíbrio, o que pode estar associado a mudanças maiores no nível de base ou reativações tectônicas na bacia do alto Jequitinhonha. Padrões de comportamento semelhantes entre diversos índices encontrados (Densidade de knickpoints, IH, Dd e o ICR) evidenciam maior taxa de dissecação na bacia do Jequitinhonha do que na bacia do Pardo (Figura 1). Os resultados apontam para a importância do fator litoestrutural na evolução do relevo da área em estudo, evidenciando ainda que as áreas de alta energia sofrem influências de um conjunto de variáveis, entre elas, o clima, a vegetação, e não apenas da energia presente no sistema de drenagem.

Agradecimentos

Agradeço a FAPEMIG pelo apoio financeiro.

Referências

AUGUSTIN, C.H.R.R. Geoökologische Studien im südlichen Espinhaçogebirge bei Gouveia, Minas Gerais, Brasilien unter besonder Berücksichtigung der Landschaftsentwicklung. 294f. Tese de Doutorado em Ciências Naturais – Departamento de Geografia, Johann Wolfgang Goethe Universität, Frankfurt am Main, Alemanha, 1995.

AUGUSTIN, C.H.R.R.; FONSECA, B.M.; ROCHA, L.C. Mapeamento geomorfológico da Serra do Espinhaço Meridional: primeira aproximação. Geonomos vol. n. 19, vol. 2, p. 50-69, 2011.

AUGUSTIN, C.H.R.R.; COE, H.H.G.; CHUENG, K.F. Analysis of geomorphic dynamics in ancient quartzite landscape using phytolith and carbon isotopes, Espinhaço Mountain Range, Minas Gerais, Brazil. Géomorphologie: relief, processus, environment vol. 4, p. 355-376, 2014.

CHEMALE JR. F.; DUSSIN, I.A.; MARTINS, M.; SANTOS, M.N. Nova abordagem tectono-estratigráfica do Supergrupo Espinhaço em sua porção meridional (MG). Geonomos, v. 19, n.2, p.173-41, 2011.
CROSBY, B.T., WHIPPLE, K.X. Knickpoint initiation and distribution within fluvial networks: 236 waterfalls in the Waipaoa River, North Island, New Zealand. Geomorphology v. 82, p. 16–38, 2006.
DUSSIN, I.A. e DUSSIN, T.M. Supergrupo Espinhaço: Modelo de evolução geodinâmica. Geonomos, Belo Horizonte, v.3, n. 1, p.19-26. 1995.
ETCHEBEHERE, M.L.C.; SAAD, A.R.; PERINOTTO, J.A.J.; FULFARO, V.J. Aplicação do Índice “Relação Declividade-Extensão – RDE” na Bacia do Rio do Peixe (SP) para detecção de deformações neotectônicas. Revista de Instituto de Geociências da USP – Séries Científica, v. 4, n.2, p.43-56. 2004.
FOGAÇA, A.C.C. Geologia da Folha Diamantina. In: GROSSI-SAD, J. H.; LOBATO, L.M.; PEDROSA-SOARES, A.C.; SOARES-FILHO, B.S. (coordenadores e editores). Projeto Espinhaço em CD-ROM (textos, mapas e anexos). Belo Horizonte, COMIG - Companhia Mineradora de Minas Gerais. 1997. 1598p.
FONSECA, B.M. O uso do Sistema de Informações Geográficas na análise morfométrica e morfológica de bacias de drenagem na Serra do Espinhaço Meridional – MG. 2010. 93 f. Dissertação (Mestrado em Análise Ambiental), programa de Pós-Graduação em Geografia, Instituto de Geociências, Universidade Federal de Minas Gerais, Belo Horizonte. 2010.
FONSECA, B.M; AUGUSTIN, C.H.R.R. Use of GIS to calculate Hack Index as a basis for comparative geomorphologic analysis between two drainage basins: a case study from SE-Brazil. In: International Geographic Union Regional Geographic Conference - UGI 2011, 2011, Santiago. Abstracts International Geographic Union Regional Geographic Conference – UGI 2011. Santiago: Military Geographic Institute of Chile (IGM), vol. 1, p. 1-12, 2011.

FONSECA, B.M.; AUGUSTIN, C.H.R.R. Análise morfométrica de bacias de drenagem e sua relação com a estrutura geológica, serra do Espinhaço Meridional MG. Revista Brasileira de Geomorfologia, v. 15, n. 2, p. 153-172. 2014.
GUEDES, I.C.; SANTONI, G.C.; ETCHEBEHERE, M.L.C.; STEVAUX, J.C.; MORALES, N.; SAAD, A.R. Análise de perfis longitudinais de drenagens da bacia do rio Santo Anastácio (SP) para detecção de possíveis deformações neotectônicas. Revista da Universidade de Guarulhos, Geociências, v. 5, n. 1, p. 75-102. 2006.
HACK, J.T. Stream-profile analysis and stream-gradient index. Journal of Research of the United States Geological Survey, v.1, n.4, p.421-429. 1973.
HACK J.T. Dynamic equilibrium and landscape evolution. In: MELHORN, W.L.; FLEMAL, R.C. Theories of Landform Development, State University of New York Press, Binghamton, NY. p. 87–102. 1975.
HORTON, R.E. Erosional development of streams their drainage basins: hidrophysical approach to quantitative morphology. Bulletin of the Geological Society of America, Colorado, v. 56, p. 275-370. 1945.
INSTITUTO MINEIRO DE GESTÃO DAS ÁGUAS – IGAM. Disponível em: <http://www.igam.mg.gov.br>. Acesso em: 10. fev. 2016.
KNAUER, L.G. O Supergrupo Espinhaço em Minas Gerais: considerações sobre sua estratigrafia e seu arranjo estrutural. Geonomos, Belo Horizonte, v.15, n.1, p. 81-90. 2007.
KNAUER,L.G. Evolução geológica do Precambriano da porção centro leste da Serra do Espinhaço Meridional e metalogênese associada.1990. 298 f. Dissertação (Mestrado em geologia), Universidade Estadual de Campinas UNICAMP, Campinas, São Paulo. 1990.

LEOPOLD, L.B.; MADDOCK, T. The hydraulic geometry of stream channels and some physiographic implications. Geological Survey Professional Paper. v.252. p. 1-57. 1953.
PHILLIPS, J. D., et al 2010, ‘Origin and interpretation of knickpoints in the Big South Fork River Basin, Kentucky-Tennessee’, Geomorphology, v. 114, n. 3, p. 188-198, 2010.
QUEIROZ, G.L.; SALAMUNI, E.; NASCIMENTO, E.R. Knickpoint finder: A software tool that improves neotectonic analysis. Computers & Geosciences, v. 76, p. 80-87, 2015.
RENGER, F.E. e KNAUER, L.G. Espinhaço – Quo vadis? (Onde está? Aonde vai?) A evolução dos conhecimentos sobre a Cordilheira do Espinhaço Meridional em Minas Gerais entre 1979 e 1995. Geonomos, v.3, n.1, p.31-39. 1995.
SALAMUNI, E.; NASCIMENTO, R.E.; SILVA, P.A.H.; QUEIROZ, G.L.; SILVA, G. Knickpoint Finder: Ferramenta para a busca de geossítios de relevante interesse para o geoturismo. Boletim Paranaense de Geociências, v. 70, p. 200-208, 2013.
SAMPAIO, T.V.M. Parâmetros morfométricos para melhoria da acurácia do mapeamento da rede de drenagem: uma proposta baseada na análise da Bacia Hidrográfica do Rio Benevente - ES. 2008. 147 f. Tese (Doutorado em Geografia) Instituto de Geociências da Universidade Federal de Minas Gerais. Belo Horizonte. 2008.
SAMPAIO T.V.M.; AUGUSTIN C.H.R.R. Índice de Concentração da Rugosidade: uma nova proposta metodológica para o mapeamento e quantificação da dissecação do relevo como subsídio a Cartografia Geomorfológica. Revista Brasileira de Geomorfologia. v. 15, n. 1, p. 47-60. 2014.
SOUZA, L. F.; SAMPAIO T. V. M. (2010). Aplicação do Índice de Concentração da Rugosidade à identificação de classes de dissecação do relevo: uma proposta de quantificação e automatização em ambiente SIG. In: III SIMPÓSIO BRASILEIRO DE CIÊNCIAS GEODÉSICAS E TECNOLOGIAS DA GEOINFORMAÇÃO. 2010, Recife-PE. Anais. P. ?.
STRAHLER, A.N. Dynamic basis Geomorfology. Geological Society of America Bulletin, v. 63, n.9, p. 923-938. 1952.
UHLEIN, A. Transição cráton-faixa dobrada: exemplo do Cráton do São Francisco e da Faixa Araçuaí (Ciclo Brasiliano) no Estado de Minas Gerais. Aspectos estratigráficos e estruturais.1991. 295f. Tese (Doutorado em geologia), Universidade Estadual de São Paulo – USP. 1991.
VALERIANO, M.M.; ROSSETTI, D.F.; ALBUQUERQUE, P.C.G. Topodata: desenvolvimento da primeira versão do banco de dados geomorfométricos locais em cobertura nacional. In: SIMPÓSIO BRASILEIRO DE SENSORIAMENTO REMOTO, 2009, Natal, RN. Anais. São José dos Campos, SP : INPE, 2009. Artigos, p. 1-8. CD-ROM.