Uso de Técnicas Geomorfométricas para Identificação dos Padrões de Relevo na Bacia do Rio Preto, no Planalto Ocidental Paulista

Resumo

A ciência geomorfológica muito se beneficia dos avanços tecnológicos, das ferramentas SIG disponibilizadas por meio de softwares livres e de uma maior acessibilidade às bases de dados planialtimétricas globais ou quase-globais em escala de semidetalhe. Consequentemente, é crescente o uso da geomorfometria aplicada à cartografia geomorfológica e otimização das rotinas de mapeamento utilizando métodos digitais semi-automatizados. O presente trabalho aplica uma classificação geomorfométrica parametrizada a partir de um Modelo Digital de Elevação para compartimentar os padrões de relevo da bacia do Rio Preto (BRP). São três os principais parâmetros morfométricos utilizados na definição dos padrões de formas de relevo, a saber, a densidade de drenagem (DD), o índice de concentração de rugosidade (ICR) e a elevação relativa (ER) - uma aproximação da amplitude altimétrica. As classes do mapa de dissecação global (DDxICR) e a ER foram usadas para criar uma chave de classificação dos modelados e tipos de relevo com base nos conhecimentos teórico e de campo. Apesar da aparente simplicidade e monotonia de um relevo predominantemente colinoso, característico do planalto ocidental paulista, os resultados demonstram riqueza nos padrões de formas, com topos aplainados e colinas suaves (65% da área de estudo), topograficamente favoráveis à expansão agrícola regional; colinas onduladas (9% da área de estudo), associadas a um maior grau de dissecação, onde predominam os processos morfodinâmicos em detrimentos dos processos pedogenéticos; relevos dissecados associados aos fundos de vale (18% da área de estudo), onde predominam processos de dissecação de natureza fluvial; e planícies fluviais (8% da área de estudo), que são feições deposicionais quaternárias localizadas nos fundos dos vales. Para além de uma caracterização morfográfica da área de estudo, este trabalho lança luz sobre alguns aspectos metodológicos do mapeamento e se pretende aplicado à avaliação da vulnerabilidade ambiental em atividades futuras.

Palavras-chave:
Mapeamento Geomorfológico; Geomorfometria; Densidade de Drenagem; Índice de Concentração de Rugosidade; Amplitude Altimétrica

Abstract

Geomorphological science benefits greatly from technological advances, such as GIS tools available through free software and the improved accessibility to global or quasi-global planialtimetric databases at a semi-detailed scale. Consequently, there is growing use of geomorphometric techniques applied to geomorphological cartography and optimization of mapping routines using semi-automated, digital methods. In this sense, this study applies a geomorphometric classification based on a Digital Elevation Model to compartmentalize the relief patterns of the Rio Preto basin (BRP). We used three main morphometric parameters to define seven landform patterns, namely, drainage density (DD), roughness concentration index (RCI), and relative elevation (RE) - an approximation of relative relief. The global dissection map (DDxRCI) and RE classes were used to create a classification key for the models and relief types based on theoretical and field knowledge. Despite the apparent simplicity and monotony of the predominantly hilly topography characteristic of the Western Plateau of the state of São Paulo, the results demonstrate a diversity in landform patterns, with planated hilltops and gentle hills (65 % of the study area), topographically favorable to the expansion of regional agricultural; steep undulating hills (9 % of the study area), associated with a higher degree of dissection, where morphodynamic processes predominate over pedogenetic processes; valley bottoms (18 % of the study area), where fluvial dissection predominates; and fluvial plains (8 % of the study area), which are quaternary depositional features located on the valley bottoms. In addition to a morphographic characterization of the study area, this study sheds light on certainsome methodological aspects of the mapping and it is intended to be applied to environmental vulnerability assessment in future activities.

Keywords:
Geomorphological Mapping; Geomorphometry; Drainage Density; Roughness Concentration; Index Relative Relief

INTRODUÇÃO

A cartografia geomorfológica é a área do conhecimento que integra a geomorfologia e a cartografia visando representar os fatos geomorfológicos segundo a gênese das formas do relevo e suas relações com a estrutura e os processos responsáveis pela sua evolução. Neste sentido, os mapas e as cartas geomorfológicas são importantes fontes de informação das potencialidades e vulnerabilidades do terreno e apoiam o planejamento ambiental e o ordenamento territorial (CHRISTOFOLETTI, 2007; GUERRA e MARÇAL, 2006; ROSS, 2006; ROSS e MOROZ, 2009; MANTOVANI, 2015; RODRIGUES et al., 2016), como nos planos de manejo de áreas protegidas, planos diretores municipais e os variados estudos ambientais solicitados durante o processo de licenciamento ambiental. Entretanto, apesar de sua grande relevância, ainda não há um consenso a respeito dos sistemas metodológicos de representação geomorfológica adotados no Brasil. A dificuldade em adotar uma padronização metodológica que subsidie um sistema brasileiro de classificação do relevo, debate reanimado a partir XVIII Simpósio Brasileiro de Geografia Aplicada (2019) na cidade de Fortaleza, é devida à grande variedade de elementos considerados (morfogênese, morfologia, morfometria, morfodinâmica, morfoestrutura e morfocronologia) e às complexidades intrínsecas aos próprios sistemas geomorfológicos.

Apesar da falta de consenso e pulverização metodológica, no Brasil há uma gravitação em torno das propostas de Jurandir Ross (1992), cuja organização taxonômica encontra suas origens nos trabalhos de Jean Tricart e se expande a partir do projeto RADAMBRASIL (NETO, 2020). O Estado de São Paulo, onde se encontra a área de estudo deste trabalho, recebeu uma das grandes contribuições de Ross, juntamente com Isabel Moroz, a saber, o “Mapa Geomorfológico do Estado de São Paulo”, do ano de 1997, em escala 1:500.000. Este mapeamento se organiza de acordo com a lógica fundamental de Tricart, relacionando as espacialidades e temporalidades segundo ordens de grandeza. Desta forma, realidades geomorfológicas generalizadas correspondem a uma maior temporalidade (táxons superiores) e possuem forte tônica estrutural e morfogenética. À medida que se avança em direção às escalas maiores (táxons inferiores), o mapeamento se aproxima do tempo “atual” e destaca os aspectos morfodinâmicos das realidades geomorfológicas. Os padrões de relevo (3º táxon), segundo a proposição de Ross (1992), agrupam os modelos por tipos genéticos (agradação e dissecação) e morfodinâmicos, deduzidos da morfometria (dimensão interfluvial e grau de entalhe). Neste sentido, em se tratando dos modelados de dissecação, a distância interfluvial representa o grau de dissecação topográfica pela atuação fluvial e/ou controle estrutural e o grau de entalhamento corresponde a profundidade das dissecações.

Os avanços tecnológicos em equipamentos e softwares e as iniciativas de democratização dos dados planialtimétricos globais, ou quase globais - como por exemplo, os Modelos Digitais de Elevação (MDEs) da Missão Topográfica Radar Shuttle / Shuttle Radar Topography Mission / SRTM - favoreceram o uso da geomorfometria aplicada à cartografia geomorfológica e consequente otimização das rotinas de mapeamento. Por conseguinte, pode-se citar as contribuições de Gustavsson (2006); Santos et al. (2006); Evans et al. (2009); Grecu (2009); Pavlopoulos et al. (2009); Silva e Rodrigues (2009); Fonseca et al. (2015); Martins e Rodrigues (2016); Sena Souza et al. (2016); Silveira e Silveira (2016); Diniz et al. (2017); Bortolini et al. (2018); Nazar e Rodrigues (2019), que produziram mapeamentos e compartimentações geomorfológicas com base nos métodos digitais e semiautomatizados em ambiente de Sistema de Informação Geográfico (SIG).

Este trabalho apresenta uma compartimentação do relevo da bacia do Rio Preto, no Planalto Ocidental Paulista, segundo uma adaptação da proposição de Ross (1992), a partir de análises geomorfométricas, utilizando métodos de mapeamento semiautomatizados que conciliam o conhecimento empírico de campo com metodologias modernas. São considerados três aspectos principais na caracterização do relevo: a declividade, a densidade de drenagem (DD) (equivalente à dimensão interfluvial), o Índice de Concentração de Rugosidade (ICR) (equivalente ao grau de entalhe topográfico) e a elevação relativa, esta última medida como uma aproximação da amplitude altimétrica. Em um contexto de grande debate acerca das metodologias utilizadas nos mapeamentos geomorfológicos, para além da questão central da taxonomia e padronização do inventário de formas de relevo, apresentamos procedimentos ainda pouco explorados pela literatura que, devido à sua simplicidade, demonstram grande potencial para uso exploratório nos mapeamentos do 3º táxon e 4º táxon.

ÁREA DE ESTUDO

A bacia do Rio Preto (BRP - 2.866 km²), afluente do Rio Turvo, está localizada na Unidade de Gerenciamento de Recursos Hídricos “Turvo-Grande” (UGRHI 15), à Noroeste do Estado de São Paulo e compreende 16 municípios, entre eles, o município de São José do Rio Preto (SP) (Figura 1). O clima regional, segundo classificação proposta por Köeppen-Geiger, é o “Aw” (tropical com inverno seco) com estação chuvosa no verão, nítida estação seca no inverno e precipitações podendo chegar a 1.800 mm/ano (ALVARES et al., 2014; FOREST-GIS, 2015; EMBRAPA, 2021).

Figura 1
Bacia do Rio Preto, afluente do Rio Turvo, localizada na UGRHI 15, com destaque para o município de São José do Rio Preto, à sudeste.

O mapa hipsométrico representa os diferentes níveis altimétricos da BRP, entre 389m e 582m e registra os principais desníveis controlados pelas variações geológicas (Figuras 2b e 2 c). Segundo Perrota et al. (2005), a área de estudo apresenta uma composição geológica simples onde observamos na porção mais baixa da bacia, próximo ao exutório, o afloramento da Formação Serra Geral, representada neste ponto por rochas básicas intercaladas por arenitos. Sobrepostas à Formação Serra Geral ficam as unidades geológicas predominantes na bacia: a Formação Vale do Rio do Peixe, composta por arenitos de muito finos a finos, bem selecionados e de formação continental desértica/eólica e a Formação São José do Vale do Rio Preto composta por arenitos de muito finos a finos, mal selecionados e de formação continental desértica/fluvial entrelaçado. No eixo da drenagem principal são observados Depósitos Aluvionares quaternários de areia quartzosa, silte e argila com formações locais de turfas.

As informações pedológicas foram obtidas do mapeamento de solos do Estado de São Paulo, revisado e unificado por Rossi (2017), em escala de 1:250.000 (Figura 2a). Na maior parte da BRP são observados os Argissolos Vermelho-Amarelos e nos topos planos os Latossolos Vermelhos. Já no vale do rio Preto se observa a ocorrência de Gleissolos Háplicos e Neossolos Flúvicos.

Figura 2
a) Pedologia; b) Hipsometria; c) Litologia

Os mapeamentos do relevo paulista remontam a Moraes Rego (1932); Monbeig (1949); Ab' Saber (1956); IBGE (1959); Almeida (1964). Já os mapeamentos geomorfológicos mais recentes e com maior detalhamento são do Instituto de Pesquisas Tecnológicas do Estado de São Paulo (IPT), assinado por Ponçano et al. (1981), de Ross e Moroz (1997), de Peixoto (2010) e do IBGE (2017). A área de estudo está localizada na unidade morfoestrutural “Bacia Sedimentar do Paraná” (1º táxon) que, segundo Ross e Moroz (1997), representa uma complexa fossa tectônica de forma elipsoidal, preenchida por sedimentos e cujo embasamento é formado principalmente de rochas cristalinas pré-cambrianas. Avançando mais um táxon, a área de estudo está contida na unidade morfoescultural "Planalto Ocidental Paulista” (2º táxon), unidade que compreende uma área de aproximadamente 100.000km², cerca de 2/5 da área do Estado de São Paulo e possui geologia relativamente simples, quase inteiramente representada pelo Grupo Bauru. O relevo é suavizado e “monótono”, com predomínio de colinas e morrotes levemente ondulados, muito favorável às atividades agrícolas e ao traçado das vias de comunicação (ALMEIDA, 1964). Segundo o autor, o Planalto Ocidental é o menos estudado entre as províncias paulistas, seja por sua aparente simplicidade ou pela distância da Capital.

MATERIAIS E MÉTODOS

A compartimentação do relevo da bacia do Rio Preto foi realizada em duas etapas: i) aquisição de dados e ii) processamento dos dados e compartimentação geomorfológica. A aquisição de dados diz respeito aos levantamentos bibliográficos/cartográficos e aos trabalhos de campo. Foram consultadas as bases de dados cartográficas disponibilizadas pela Companhia de Pesquisas de Recursos Minerais (CPRM) e pelo Instituto Brasileiro de Geografia e Estatística (IBGE), assim como as imagens de alta resolução disponibilizadas pelo Google Maps/Earth através de ferramentas de geoprocessamento que integram essas informações. O Modelo Digital de Elevação (MDE) utilizado nas análises geomorfométricas foi o SRTM (Shuttle Radar Topography Mission) de resolução de 1 segundo de arco (≅ 30m) (NASA, 2009). Os trabalhos de campo foram realizados em janeiro de 2021 e fazem parte da investigação empírica, essencial para a compreensão das características morfográficas e morfodinâmicas da área de estudo, obtenção de pontos de controle utilizando GNSS (Global Navigation Satellite System) e registros fotográficos.

A compartimentação geomorfológica se baseou nos pressupostos de Ross (1992), sendo abordado o táxon que trata dos padrões de formas (3º táxon) e o 4º táxon, que representa as formas individualizadas, neste caso, indicadas no conjunto dos padrões de formas. Os parâmetros morfométricos utilizados para a classificação dos modelados de dissecação foram a declividade e amplitude altimétrica, tendo sido observado os critérios apresentados por CPRM/IPT (2014) (Tabela 1).

Tabela 1
Padrões de Relevo (modelados de dissecação) e parâmetros básicos

TIPO DE DE RELEVO	AMPLITUDE PREDOMINANTE (m)	DECLIVIDADE PREDOMINANTE (%)
Colinas	40 a 70	< 20
Morrotes baixos	40 a 50	> 30
Morrotes	60 a 90	20
Morros baixos	90 a 110	30
Morrotes altos	60 a 90	30 a 40
Morros altos	140 a 200	> 30
Serras	> 300	> 30
Escarpas	100	60

• Fonte: Adaptado de CPRM/IPT (2014).

Neste trabalho calculamos a elevação relativa acima de uma drenagem (ER) matematicamente modelada como uma aproximação da amplitude altimétrica utilizando o MDE em ambiente SIG. Para tanto, utilizamos o módulo Relative Heights and Slope Positions (BOEHNER e SELIGE, 2006), no software livre SAGA GIS (CONRAD et al., 2015). Boehner e Selige (2006) estimaram a elevação acima das áreas de drenagem acumulada (“altitude above drain culmination” / $A{D}_M$) de forma contínua na matriz raster, sem usar dados discretos, como linhas de drenagem ou cristas, segundo a Equação 1. As elevações relativas são designadas como sendo a diferença entre as células da matriz raster de elevação invertida $zi$ e a média ponderada da elevação invertida à montante $z{i}_i$, cada uma ponderada pela raiz quadrada de sua respectiva área de contribuição $C{A}_i^{}$. A área de contribuição foi calculada usando o método de fluxo múltiplo de Freeman (1991). Já as declividades foram calculadas a partir do MDE em ambiente SIG, sendo este procedimento entre os mais elementares e amplamente divulgado.

$A{D}_M=-1.\left[\frac{{{\displaystyle \sum }}_{i=l}^{n}z{i}_i/C{A}_i^{\mathrm{0,5}}}{{{\displaystyle \sum }}_{i=l}^{n}1/C{A}_i^{\mathrm{0,5}}}-z{i}_{}\right]$ (Equação 1)

Seguindo a proposição de Nazar e Rodrigues (2019), foram calculadas a densidade de drenagem (DD) e o índice de concentração de rugosidade (ICR) (SAMPAIO e AUGUSTIN, 2014). A DD é inversamente proporcional à distância interfluvial e, desta forma, quanto maior a DD, menor será a distância interfluvial e maior será a dissecação fluvial. O ICR, ou densidade da declividade, corresponde aos diferentes graus de aprofundamento do vale em 3 dimensões. Ambos os índices foram calculados em ambiente SIG, utilizando um estimador de densidade por kernel (Kernel Density Estimation / KDE) (WANDERLEY, 2013) com raio de 564m (1km²) (HORTON, 1945), como discutido por Sampaio e Augustin (2014) e utilizado por Nazar e Rodrigues (2019) e Nazar (2018). Após aplicar o KDE, os mapas foram reclassificados em uma escala numérica simples (de 1 a 5), sendo (1) muito baixo, (2) baixo, (3) médio, (4) alto e (5), muito alto, respectivamente (NAZAR e RODRIGUES, 2019).

Os índices DD e ICR foram multiplicados por meio de álgebra de mapas, gerando um mapa de dissecação global com classes que variam entre os múltiplos de 1 a 5. As classes do mapa “DDxICR” foram usadas, juntamente com a ER para criar uma chave de classificação com base nos conhecimentos teórico e de campo dos executores do mapeamento (Tabela 2). Para a identificação do modelado de agradação, além da chave de classificação também foram utilizadas as curvas de nível, a rede de drenagem e as imagens de satélite. Para maior detalhamento operacional, consultar Nazar e Rodrigues (2019).

RESULTADOS E DISCUSSÕES

Nesta seção apresentaremos os mapas derivados do MDE, a saber os mapas declividade, ER, DD e ICR e, posteriormente, detalharemos o mapa de padrões de relevo da BRP e suas 5 classes produzidas a partir da chave de classificação geomorfométrica, conforme metodologia previamente descrita.

A declividade foi classificada de acordo com a Tabela 1. A maior parte das células da matriz de declividade se encontram na classe de 0 a 8% (87%), sendo o valor médio de declividade de 4% e a declividade máxima de 33% (Figuras 3a). O mapa de ER, computado pela Equação 1, também foi classificado segundo a Tabela 1 como uma aproximação da distribuição espacial das amplitudes altimétricas na BRP. A maior parte das células da matriz estão entre as amplitudes de 0 e 20m (63%), seguido da classe de 20m e 40m (33%). Na área de estudo, a amplitude máxima registrada é de 75m e o valor médio de amplitude é de 17m (Figura 3b).

Ao confrontarmos os mapas de declividade e ER, observamos que o único padrão de relevo observado é o de colinas (Tabela 1) (CPRM/IPT, 2014), uma vez que estas possuem uma amplitude altimétrica entre 40m e 70m e declividade predominante <20%. Não foram observadas combinações entre elevações superiores a 40m com declividades superiores a 20%, como nos morros baixos e/ou morrotes.

As variáveis DD e ICR refletem os aspectos de dissecação da BRP em cinco classes, que variam de muito alto a muito baixo (Figuras 3c e 3d). Os mapas possuem uma distribuição regular das classes devido ao método de classificação utilizado (Quantil), sugerido por Sampaio e Augustin (2014), utilizado também por Nazar e Rodrigues (2019).

Figura 3
Mapas computados pelo MDE. Em (a), mapa de Declividade; em (b), mapa de Elevação Relativa; em (c), mapa de Densidade de Drenagem; em (d) mapa do Índice de Concentração de Rugosidade.

Padrões de Relevo da Bacia do Rio Preto

São cinco os padrões de relevo da BRP identificados partir da chave de classificação apresentada na Tabela 2: i) relevo aplainado associado a interflúvios (10%); ii) relevo dissecado associado a colinas suaves (55%); iii) relevo dissecado associado a colinas onduladas (9%); iv) relevo dissecado associado aos fundos de vale (18%); v) e planícies fluviais e planícies fluviais isoladas (8%).

Tabela 2
Chave de classificação dos padrões de relevo da Bacia do Rio Preto (SP)

Índice (DDxICR)	ER (m)	Nomenclatura	BRP (%)
1	>7	Relevo aplainado associado a interfúvlios (Int)	10
2,3,4,5,6,8,9,10,12,15,16	>5	Relevo dissecado associado às colinas suaves (Dc)	55
20, 25	>5	Relevo dissecado associado às colinas onduladas (Dco)	9
3,4,5,6,8,9,10,12,15,16,20,25	≤5	Relevo dissecado associado aos fundos de vale (Dv)	18
1,2	≤5	Planícies fluviais e planícies fluviais isoladas (Apl)	8

• Fonte: Os autores (2021).

Apesar da aparente monotonia e simplicidade do relevo colinoso do planalto ocidental paulista, este mapeamento apresenta um ganho de detalhamento em relação aos trabalhos de Ponçano et al. (1981), Ross e Moroz (1997), Peixoto (2010) e IBGE (2017) que se reflete na diversidade dos padrões de relevo na BRP.

De forma geral, a BRP apresenta um relevo constituído por colinas suaves à onduladas, com vertentes convexo-côncavas e topos arredondados e/ou aplainados com morfologia alongada. Já as planícies fluviais e planícies fluviais isoladas representam as áreas de agradação quaternária associadas a ambientes de baixa energia com até 5m de amplitude altimétrica.

A Figura 4 apresenta o mapa das classes dos padrões do relevo, conforme apresentado na Tabela 2, onde foram traçados dois perfis topográficos representativos das classes mapeadas. Na seção a-b (perfil 1), em “a”, observa-se um relevo pouco dissecado associado a uma colina suave convexa (cota 470m), próximo do exutório da BRP. O perfil atravessa a planície fluvial do Rio Preto (cota 420m), seguindo por um compartimento de colinas suaves na formação do Vale do Rio do Peixe, até alcançar as cotas mais elevadas (520m), na formação do Vale do Rio Preto, onde ocorrem colinas onduladas, intercaladas por vale dissecados; finalmente chegando em “b”, em um interflúvio aplainado à 520m de elevação.

Na seção c-d (perfil 2), em “c” há um compartimento mais elevado, na formação do Vale do Rio Preto, onde se observa a ocorrência de um interflúvio aplainado com morfologia alongada na cota 540m, que se estende até o contato transicional entre as formações geológicas, à 1,2km de distância do ponto “c”. No contato observa-se uma mudança do nível de base e ocorrência de colinas onduladas e fundo de vale dissecado. As cotas mais baixas correspondem à uma planície fluvial do Rio Preto (440m) e, seguindo para os limites da bacia à leste, onde o ponto “d” encontra colinas onduladas.

Figura 4
Representação dos padrões de relevo da BRP e perfis topográficos: perfil 1 (seção a-b) e perfil 2 (seção c-d).

Relevo Aplainado Associado a Interflúvios

As formas de relevo aplainadas associadas a interflúvios apresentam morfologia tabular, por vezes alongadas, com baixa rugosidade. Esta classe corresponde a 10% da área da BRP, estando associada a solos profundos e bem drenados, majoritariamente os Latossolos Vermelhos. Importante ressaltar que estudos acerca da dinâmica geossistêmica e fragilidades ambientais do Oeste Paulista (PEREZ FILHO; SEABRA, 2004; SEABRA, 2014; QUARESMA, 2008) têm verificado processos morfogenéticos associados à erosão pluvial e eólica, esta última contribuindo para a instalação de processos de arenização, decorrentes da dinâmica inadequada de uso e ocupação do solo com retirada e substituição da vegetação nativa. A Figura 5 representa o relevo aplainado dos topos, sendo perceptível a presença de processos erosivos eólicos, associada aos arenitos da área.

Figura 5
Interflúvios aplainados no interior da Bacia do Rio Preto.

Relevo Dissecado Associado à Colinas Suaves

Este padrão de relevo ocorre na maior parte da BRP (55%), em encostas pouco dissecadas com curvatura convexa em gradiente suave. É observado principalmente nas bordas dos topos planos, caracterizando uma transição de um relevo plano para um relevo de colinas onduladas ou planícies fluviais. Neste padrão de relevo se observa tanto os processos de pedogênese quanto os processos de morfogênese, com predominância de erosão do tipo laminar sobre os Argissolos Vermelho-Amarelos profundos e bem drenados (Figura 6).

Figura 6
Relevo dissecado associado a colinas suaves no interior da Bacia do Rio Preto.

Os padrões de relevo com topografia suave, a saber os interflúvios aplainados e as colinas suaves, predominantes na área de estudo (65%), facilitam o manejo agrícola na região. Segundo dados obtidos do projeto “Monitoramento da Cobertura e Uso da Terra do Brasil” (IBGE, 2018), a maior parte da BRP é composta por área agrícola e, de acordo com o Plano Estadual de Recursos Hídricos (COBRAPE, 2020), há projeções de expansão dessas áreas. Contudo, a grande pressão agrícola proporcionada pelas condições topográficas favoráveis é objeto de preocupação uma vez que a substituição da vegetação original pode alterar os processos hidrológicos, pedológicos e geomorfológicos e, consequentemente, as trocas de energia/matéria, gerando instabilidades na paisagem (TRICART, 1977; ROSS, 1994), reduzindo a recarga do aquífero e aumentando os processos erosivos.

Relevo Dissecado Associado a Colinas Onduladas

Este padrão de relevo ocorre em 9% da BRP e possui os maiores valores de dissecação global (ICRxDD), sendo frequentemente localizado junto aos vales dissecados. Nesta classe predominam os processos de morfogênese em detrimento dos processos pedogenéticos. O maior grau de dissecação é observado, frequentemente, no terço inferior das vertentes onde também são observados processos difusivos como o rastejo (Figura 7). Nestas áreas se observa uma associação de Argissolos Vermelho-Amarelos e Neossolos Litólicos.

Figura 7
Relevo dissecado associado a colinas onduladas no interior da Bacia do Rio Preto.

Relevos Dissecados Associados aos Fundos de Vale

Nos fundos de vale ocorre o padrão de relevo mais ou menos dissecados (18%), conforme o índice global de dissecação (ICRxDD) da Tabela 2, onde predominam os processos de morfogênese. Este padrão de relevo é observado nas concavidades, entre amplitudes altimétricas menores ou iguais a 5m, podendo apresentar menor grau de dissecação, conforme se pode observar na figura 6, ou incisões mais acentuadas, de natureza fluvial ou por erosão linear, frequentemente associadas aos graus de dissecação ICRxDD = 20 e 25 (Tabela 2), como observado na figura 8.

Figura 8
Relevo dissecado associado aos fundos de vale no interior da Bacia do Rio Preto.

A intensidade da dissecação pode ser considerada um primeiro grande indicador da fragilidade potencial do ambiente natural, conforme Ross (1994), sugerindo proteção e/ou restrição do uso e ocupação destas áreas.

Planícies Fluviais

Trata-se de feições deposicionais quaternárias formadas a partir dos detritos depositados durante as cheias dos rios (8%). São áreas aplainadas nos fundos dos vales e, portanto, associadas a valores baixos de dissecação (Tabela 2). Imagens de satélite do Google Earth revelam feições importantes na delimitação dessas áreas como meandros abandonados e outras paleoformas vinculadas aos processos fluviais produzidas pelos eventos de cheias e vazantes (Figura 9).

Figura 9
Planícies fluviais e planícies fluviais isoladas da bacia do Rio Preto.

As planícies fluviais ou planícies de inundação merecem especial atenção pois são regiões de contato entre biomas, ambientes complexos, frágeis e importantes para o equilíbrio dos sistemas ambientais físico e biótico. As mudanças no uso e ocupação da BRP podem alterar esses sistemas, modificando os padrões de inundações, aumentando a sedimentação, etc.

CONSIDERAÇÕES FINAIS

A cartografia geomorfológica é uma importante ferramenta de planejamento territorial frente a um contexto de mudanças ambientais aceleradas que podem alterar os processos hidrológicos, pedológicos e geomorfológicos gerando instabilidades na paisagem. Na área de estudo, em especial, a expansão agrícola, facilitada pelas características topográficas do terreno, são fonte de preocupação, principalmente no que diz respeito às mudanças de uso e cobertura que podem alterar a recarga do aquífero e gerar aumento dos processos erosivos.

Este trabalho apresenta uma compartimentação do relevo da BRP, segundo uma adaptação da proposição de Ross (1992), a partir de análises geomorfométricas, utilizando métodos de mapeamento semiautomatizados que conciliam o conhecimento empírico de campo com procedimentos metodológicos ainda pouco explorados pela literatura. A técnica permite diminuir o grau de abstração devido ao menor uso de informações de natureza qualitativa. Contudo, os procedimentos foram complementados com trabalhos de campo, essenciais em qualquer pesquisa geomorfológica, e demostram eficiência ao trazer dados mais consistentes e precisos para o mapeamento.

Cinco padrões de relevo (equivalentes ao 3º táxon de Ross) foram identificados a partir da chave de classificação que leva em consideração a gênese, a dissecação global e a ER das formas de relevo. A maior parte da área de estudo possui relevo aplainado e/ou associado às colinas suaves (65%) onde se observa a expansão das atividades agrícolas na BRP. O relevo dissecado associado a colinas onduladas corresponde a 9% da área de estudo sendo naturalmente zonas de maior fragilidade potencial, qualidade inerente ao próprio sistema, resultante de processos morfodinâmicos e morfogenéticos. As áreas dissecadas associadas aos fundos de vale (18%) são áreas importantes para a manutenção da qualidade hídrica regional e, portanto, devem ser pensadas como locais a serem conservados. Já as planícies fluviais e as planícies fluviais isoladas correspondem a 8% da área de estudo e são ambientes muito frágeis e importantes para o equilíbrio dos sistemas ambientais físico e biótico da BRP.

Finalmente, entendemos como necessário o avanço nas pesquisas que buscam otimizar as rotinas de mapeamento geomorfológico através de métodos semiautomatizados em ambiente SIG, sem perder de vista a importância do trabalho e campo. Ainda, recomendamos que diferentes métodos de mapeamento da amplitude altimétrica sejam testados, assim como diferentes técnicas de classificação de mapas. O mapa dos padrões de relevo da BRP será utilizado no mapeamento de fragilidade ambiental emergente (ROSS, 1994), que utiliza também os mapas de pedologia e uso e cobertura.

REFERÊNCIAS

• AB' SABER, A.N. Um conceito de Geomorfologia a serviço das pesquisas sobre o quaternário. Geomorfologia, São Paulo: IGEOG-USP, n.18, p. 1-23, 1969.
• AB' SABER, A. N. A terra paulista. Boletim Paulista de Geografia, v. 23, p. 5-38, 1956.
• ALMEIDA, F. F. M. Fundamentos geológicos do relevo paulista. Boletim do Instituto Geográfico e Geológico, n. 41, 1964.
• ALVARES, C. A. et al. Köppen's climate classification map for Brazil. Meteorologische Zeitschrift, v. 22, n. 6, p. 711-728, 2014. https://doi.org/10.1127/0941-2948/2013/0507
  » https://doi.org/10.1127/0941-2948/2013/0507
• BOEHNER, J.; SELIGE, T. Spatial prediction of soil attributes using terrain analysis and climate regionalisation. In: Boehner, JAnalysis and Modelling Applications. In: BOEHNER, J.;MCCLOY, K. R., et al (Ed.). SAGA - Analysis and Modelling Applications, 2006. p.13-28.
• BORTOLINI, W. et al. Técnicas Geomorfométricas Para a Identificação de Padrões de Relevo: Aplicação nas Cartas Curitiba e Cerro Azul, Estado do Paraná. Revista do Departamento de Geografia (Universidade de São Paulo), v. 36, p. 15-32, 2018. https://doi.org/10.11606/rdg.v36i0.144285
  » https://doi.org/10.11606/rdg.v36i0.144285
• CHRISTOFOLETTI, A. Aplicabilidade do conhecimento geomorfológico nos projetos de planejamento. In: GUERRA, A. J. T.; CUNHA, S. B. Geomorfologia: uma atualização de bases e conceitos. Rio de Janeiro: Bertrand Brasil, 7ª ed. p. 365-39, 2007.
• COBRAPE. Diagnóstico Síntese do Relatório Final do Plano Estadual de Recursos Hídricos (PERH) 2020-2023. São Paulo: Companhia Brasileira de Projetos e Empreendimentos (COBRAPE); Secretaria de Infraestrutura e Meio Ambiente (SIMA), 2020. 469.
• CONRAD, O. et al. System for Automated Geoscientific Analyses (SAGA) v. 2.1.4. Geosci. Model Dev., v. 8, n. 7, p. 1991-2007, 2015. https://doi.org/10.5194/gmd-8-1991-2015
  » https://doi.org/10.5194/gmd-8-1991-2015
• CPRM/IPT. Cartas de suscetibilidade a movimentos gravitacionais de massa e inundações: 1:25.000 (livro eletrônico): nota técnica explicativa. Brasília, DF: CPRM (Serviço Geológico do Brasil) & IPT (Instituto de Pesquisas Tecnológicas do Estado de São Paulo): 42 p. 2014.
• DAEE. Departamento de Águas e Energia Elétrica. Base cartográfica digital, escala 1:50.000 - Projeto GISAT. São Paulo, 2008. Disponível em http://datageo.ambiente.sp.gov.br/ Acesso em: 15 out. 2020
  » http://datageo.ambiente.sp.gov.br/
• DAL RÉ CARNEIRO, C. Os “fundamentos geológicos do relevo paulista” nos dias atuais. Revista do Instituto Geológico, v. 39, n. 3, 2018. https://doi.org/10.33958/revig.v39i3.599
  » https://doi.org/10.33958/revig.v39i3.599
• DINIZ, M. T. M. et al. Mapeamento Geomorfológico do Estado do Rio Grande do Norte. Revista Brasileira de Geomorfologia, v. 18, p. 689-701, 2017. https://doi.org/10.20502/rbg.v18i4.1255
  » https://doi.org/10.20502/rbg.v18i4.1255
• EMBRAPA. Clima. s/d., 2021. Disponível em: https://www.cnpf.embrapa.br/pesquisa/efb/clima.htm Acesso em: 26 març. 2021.
  » https://www.cnpf.embrapa.br/pesquisa/efb/clima.htm
• EVANS, I. S.; HENGL, T.; GORSEVSKI, P. Applications in Geomorphology In: (Ed.). Geomorphometry: concepts, software, applications: Elsevier, v.33, 2009. cap. 22, p.497-525. (Developments in Soil Science). https://doi.org/10.1016/S0166-2481(08)00022-6
  » https://doi.org/10.1016/S0166-2481(08)00022-6
• FONSECA, A. F.; ZÊZERE, J. L.; NEVES, M. Contribuição para o conhecimento da geomorfologia da Cadeia da Arrábida (Portugal). Revista Brasileira de Geomorfologia, v. 16, p. 137-163, 2015. http://dx.doi.org/10.20502/rbg
  » http://dx.doi.org/10.20502/rbg
• FOREST-GIS. Classificação climática de Köppen-Geiger em shapefile. 2015. Disponível em: https://forest-gis.com/2015/10/classificacao-climatica-de-koppen-geiger-em-shapefile.html/. Acesso em: 20 jan. 2020.
  » https://forest-gis.com/2015/10/classificacao-climatica-de-koppen-geiger-em-shapefile.html
• FREEMAN, T. G. Calculating Catchment Area with Divergent Flow Based on a Regular Grid. Computer&Geosciences. 17: 413-422 p. 1991. https://doi.org/10.1016/0098-3004(91)90048-I
  » https://doi.org/10.1016/0098-3004(91)90048-I
• GRECU, F. Geomorphological Map of the Prahova Subcarpathians (Romania). Journal of Maps, v. 5, n. 1, p. 108-116, 2009/01/01 2009. https://doi.org/10.4113/jom.2009.1052
  » https://doi.org/10.4113/jom.2009.1052
• GUERRA, A. J. T.; MARÇAL, M. S., Eds. Geomorfologia Ambiental. Rio de Janeiro: Bertrand Brasil, p.192. 2006.
• GUSTAVSSON, M. Development of a Detailed Geomorphological Mapping System and GIS Geodatabase in Sweden. 2006. 61 Uppasala University
• HORTON, R.E. Erosional development of streams and their drainage basins: hydrophysical approach to quantitative morphology. Bulletin of the Geological Society of America V 56, p. 275-370, 1945. https://doi.org/10.1130/0016-7606(1945)56[275:EDOSAT]2.0.CO;2
  » https://doi.org/10.1130/0016-7606(1945)56[275:EDOSAT]2.0.CO;2
• IBGE. Atlas do Brasil. Rio de Janeiro: Instituto Brasileiro de Geografia e Estatística, 1959.
• IBGE. Unidades de relevo. In: IBGE, Atlas Nacional Digital do Brasil. Rio de Janeiro, 2017. Disponível em: https://www.ibge.gov.br/apps/atlas_nacional/ Acesso em: 26 jan. 2019.
  » https://www.ibge.gov.br/apps/atlas_nacional
• IBGE. Monitoramento da Cobertura e Uso da Terra do Brasil 2014-2016. Rio de Janeiro: Instituto Brasileiro de Geografia e Estatística, 2018. 29.
• ICG. Instituto Geográfico e Cartográfico do Estado de São Paulo. Sedes municipais, escala 1:50.000. São Paulo, 2010. Disponível em http://datageo.ambiente.sp.gov.br/ Acesso em: 15 out. 2020.
  » http://datageo.ambiente.sp.gov.br
• ICG. Instituto Geográfico e Cartográfico do Estado de São Paulo. Limite de UGRHI, escala 1:1.000.000. São Paulo, 2011. Disponível em http://datageo.ambiente.sp.gov.br/ Acesso em: 15 out. 2020
  » http://datageo.ambiente.sp.gov.br
• MANTOVANI, J. D. C. Geomorfologia aplicada aos estudos de impacto ambiental de empreendimentos rodoviários: subsídios aos órgãos ambientais brasileiros. 2015. 203 Dissertação (Mestrado). Departamento de Geografia, Universidade de São Paulo, São Paulo, SP.
• MARTINS, T. I. S.; RODRIGUES, S. C. Compartimentação geomorfológica da folha Piumhi, região do Alto São Francisco, Minas Gerais. Revista Brasileira de Geomorfologia, v. 17, n. 1, p. 145-162, 2016. https://doi.org/10.20502/rbg.v17i1.873
  » https://doi.org/10.20502/rbg.v17i1.873
• MONBEIG, P. A divisão regional do Estado de São Paulo. Anais dos Geógrafos Brasileiros, v. 1, p. 19-36, 1949.
• MORAES REGO, L. F. Notas sobre a geomorfologia de São Paulo e sua gênesis. São Paulo: IAG, 1932. 43p.
• NASA. Shuttle Radar Topography Mission (SRTM), 2009. https://doi.org/10.5066/F7PR7TFT
  » https://doi.org/10.5066/F7PR7TFT
• NAZAR, T. I. S. M. O Chapadão do Diamante na Serra da Canastra/MG, Brasil :caracterização geomorfológica e análise integrada do meio físico a partir de dados multifontes. 2018. 270p. Tese (Doutorado). Programa de Pós-Graduação em Geografia, Universidade Federal de Uberlândia, MG.
• NAZAR, T. I. S. M.; RODRIGUES, S. C. Relevo do Chapadão do Diamante, Serra da Canastra/MG, Brasil: Compartimentação e Análise a partir dos aspectos Geomorfométricos. Revista Brasileira de Geomorfologia, v. 20, n. 1, p. 69-88, 2019. https://doi.org/10.20502/rbg.v20i1.1300
  » https://doi.org/10.20502/rbg.v20i1.1300
• NETO, R. M. Cartografia geomorfológica: revisões, aplicações e proposições. Curitiba CRV, 2020. 174p.
• PAVLOPOULOS, K.; EVELPIDOU, N.; VASSILOPOULOS, A., Eds. Mapping Geomorphological Environments: Springered. 2009. https://doi.org/10.1007/978-3-642-01950-0
  » https://doi.org/10.1007/978-3-642-01950-0
• PEIXOTO, C. A. B., Ed. Geodiversidade do Estado de São Paulo. São Paulo: CPRM Serviço Geológico do Brasil, p.176. 2010.
• PEREZ FILHO, A.; SEABRA, F. B. Análise sistêmica aplicada ao estudo da fragilidade de terras do cerrado paulista. In: V SIMPÓSIO NACIONAL DE GEOMORFOLOGIA I ENCONTRO SUL-AMERICANO DE GEOMORFOLOGIA, UFSM - RS, 02 a 07 de Agosto de 2004, p. 1-12.
• PERROTA, M. M.; SALVADOR, E. D.; LOPES, R.C.; D’AGOSTINO, L. Z.; PERUFFO, N.; GOMES, S. D.; SACHS, L. L. B.; MEIRA, V. T.; GARCIA, M. G. M.; LACERDA FILHO, J. V. Mapa geológico do Estado de São Paulo, escala 1:750.000. Programa de Geologia do Brasil - PGB, CPRM, São Paulo 2005. Disponível em http://datageo.ambiente.sp.gov.br/ Acesso em: 15 out. 2020
  » http://datageo.ambiente.sp.gov.br/ Acesso em: 15 out. 2020
• PONÇANO, W. L. et al. Mapa Geomorfológico do Estado de São Paulo. São Paulo: Instituto de Pesquisas e Tecnologias (IPT), 1981. 94.
• QUARESMA, C. C. Organizações espaciais físico/naturais e fragilidades de terras sob cerrado: abordagem sistêmica aplicada à escala local. 152 f. Dissertação (Mestrado em Geografia) Universidade Estadual de Campinas, Campinas, São Paulo -Brasil, 2008.
• RODRIGUES, J. D. O. N. et al. Licenciamento ambiental: a geomorfologia nos estudos de impactos ambientais e o planejamento ambiental para ocupação de novas áreas. XIV ENEEAmb, II Fórum Latino e I SBEA - Centro-Oeste. Brasília: 634-641 p. 2016. https://doi.org/10.5151/engpro-eneeamb2016-la-003-5031
  » https://doi.org/10.5151/engpro-eneeamb2016-la-003-5031
• ROSS, J. L. S. O registro cartográfico dos fatos geomorfológicos e a questão da taxonomia do relevo. Revista do Departamento de Geografia, v. 6, p. 17-29., 1992. https://doi.org/10.7154/RDG.1992.0006.0002
  » https://doi.org/10.7154/RDG.1992.0006.0002
• ROSS, J. L. Análise empírica da fragilidade dos ambientes naturais e antropizados. Revista do Departamento de Geografia, n. 8, p. 63-74, 1994. https://doi.org/10.7154/RDG.1994.0008.0006
  » https://doi.org/10.7154/RDG.1994.0008.0006
• ROSS, J. L. Geomorfologia aplicada aos EIAs-RIMAs In: GUERRA, A. J. T. e CUNHA, S. B. D. (Ed.). Geomorfologia e meio ambiente. 6ª. Rio de Janeiro, RJ: Bertrabd Brasil, 2006. cap. 6, p.291-336.
• ROSS, J. L. S.; MOROZ, I. C. Mapa geomorfológico do Estado de São Paulo: escal 1:500.000. São Paulo: Departamento de Geografia-FFLCH-USP/IPT/FAPESP., 1997. 41-58 https://doi.org/10.7154/RDG.1996.0010.0004
  » https://doi.org/10.7154/RDG.1996.0010.0004
• ROSS, J. L. Aplicabilidade do conhecimento geomorfológico nos projetos de planejamento. In: (Ed.). Geomorfologia: exercícios, técnicas e aplicações. 3ª ed. Rio de Janeiro, RJ: Bertrand Brasil, 2009. cap. 11, p.311-334.
• ROSSI, M. Mapa pedológico do Estado de São Paulo: revisado e ampliado. São Paulo: Governo do Estado de São Paulo; Instituto Florestal, 2017. 118p.
• 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. https://doi.org/10.20502/rbg.v15i1.376
  » https://doi.org/10.20502/rbg.v15i1.376
• SANTOS, L. J. C. et al. Mapeamento Geomorfológico do Estado do Paraná. Revista Brasileira de Geomorfologia, v. 7, n. 2, p. 03-12, 2006. https://doi.org/10.20502/rbg.v7i2.74
  » https://doi.org/10.20502/rbg.v7i2.74
• SEABRA, F. B. Análise Geossistêmica Aplicada ao Estudo da Fragilidade das Terras em Áreas do Cerrado Paulista. 2014. 119 f. Dissertação (Programa de Pós-Graduação em Geografia - Instituto de Geociências). Universidade Estadual de Campinas, Campinas - São Paulo, Brasil.
• SENA SOUZA, J. P. et al. Comparação Dos métodos de classificação por ângulo espectral e distância euclidiana no mapeamento das formas de terreno. Revista Brasileira de Geomorfologia, v. 17, n. 3, 2016. https://doi.org/10.20502/rbg.v17i3.846
  » https://doi.org/10.20502/rbg.v17i3.846
• SILVA, T. I.; RODRIGUES, S. C. ELABORAÇÃO DE UM TUTORIAL DE CARTOGRAFIA GEOMORFOLÓGICA COMO ALTERNATIVA PARA O ENSINO DE GEOMORFOL. Rev. Geogr. Acadêmica, v. 3, n. 2, p. 85-94, 2009.
• SILVEIRA, R. M. P.; SILVEIRA, C. T. D. Análise digital do relevo aplicada à cartografia geomorfológica da porção central da Serra do Mar Paranaense. Revista Brasileira de Geomorfologia, v. 17, n. 4, p. 615-629, 2016. https://doi.org/10.20502/rbg.v17i4.1063
  » https://doi.org/10.20502/rbg.v17i4.1063
• TRICART, J. Ecodinâmica. Rio de Janeiro: Instituto Brasileiro de Geografia e Estatística (IBGE), 1977. 91.
• WANDERLEY, M. F. B. Estudos em Estimação de Densidade por Kernel: Métodos de Seleção de Características e Estimação do Parâmetro Suavizador. 2013. 96 Tese (Doutorado). Programa de Pós-Graduação em Engenharia Elétrica, Universidade Federal de Minas Gerais, Belo Horizonte (MG).
• WEISS, A. D. Topographic Position and Landforms Analysis. The Nature Conservancy, Northwest Division, 2001, Seattle, WA.

Datas de Publicação

• Publicação nesta coleção
  12 Set 2022
• Data do Fascículo
  2022

Histórico

• Recebido
  24 Fev 2022
• Aceito
  06 Jun 2022
• Publicado
  02 Ago 2022