Parâmetros dos solos voltados para modelagem de fluxos subterrâneos
relacionados aos processos de voçorocamento: bacia do rio Bananal, SP/RJ.
Thiago Aurealino Mulato ([email protected]) Aluno de graduação em Geografia-UFRJ Bolsista: CNPq - PIBIC
Paulo Jorge Vaitsman Leal – Doutorando PPGG André de Souza Avelar – Prof. Adjunto Depto. Geografia
Ana Luiza Coelho Netto – Profa. Titular Depto. Geografia, Pesquisadora 1B-CNPq Financiamento: CNPq (Universal, PRONEX e CT-Hidro), FAPERJ e FUJB/UFRJ
Laboratório de Geo-Hidroecologia – Universidade Federal do Rio de Janeiro, Instituto de Geociências, Departamento de Geografia. Av. Athos da Silveira Ramos, 274 – Edifício do Centro de Ciências Matemáticas e
da Natureza, Cidade Universitária – Ilha do Fundão – CEP 21.941-916 – Rio de Janeiro/RJ – Brasil
RESUMO O presente trabalho está associado a uma tese de doutorado que investiga a rede de fluxos em cabeceira de drenagem de primeira ordem, na bacia do rio Bananal (RJ/SP), aonde o voçorocamento se expandiu entre 1985 e 2001, e hoje se encontra estabilizado (Rocha-Leão, 2005). Tem como objetivo principal mensurar os parâmetros dos solos que contribuem para a modelagem de fluxos subterrâneos e seus processos erosivos correlatos ('piping' e voçorocamento), tal como vem sendo abordado pelo programa de pesquisas do Laboratório de Geo-Hidroecologia da UFRJ. Pretende-se criar subsídios para o entendimento da modelagem de fluxos subterrâneos, utilizando o conhecimento do comportamento hidrológico de solos segundo diferentes metodologias de análise granulométrica, comparando estes resultados com os encontrados nos testes de condutividade hidráulica realizados em campo. Amostras de solo foram coletadas durante sondagens para instalação de aparelhos de monitoramento hidrológico e foram analisadas granulometricamente segundo o método da ABNT (norma técnica NBR 7181/84). Algumas amostras foram analisadas por outra metodologia referente a granulometria floculada, que compreende um método modificado sem uso de defloculante no processo de sedimentação, como descrito por Moraes Silva (2000). Em campo foram mensurados valores de condutividade hidráulica através do "slug test" de injeção (Hvorslev, 1951), através dos piezômetros previamente instalados. Na modelagem de processos hidrológicos, os parâmetros referentes ao solo são determinantes nas análises geradas pelo modelo, portanto é necessário considerar não somente a granulométria dos solos, como também a estrutura que o mesmo possui em campo, assim como os valores de condutividade hidráulica referente aos pontos de monitoramento. Os resultados das análises granulométricas mostram que não há significativa variação entre os solos desta cabeceira de drenagem. Além disso, os resultados dos ensaios de condutividade hidráulica também mostraram-se muito similares, com ordem de grandeza de valores entre 10-2 a 10-4cm/s.
PALAVRAS-CHAVE: Análise Textural, Condutividade Hidráulica e Água Subterrânea.
ABSTRACT This work is associated with a doctoral thesis that investigates the network of headwater streams in the Bananal river basin (RJ / SP), where the erosion process (Gully) expanded between 1985 and 2001, and today is stabilized (Rocha-Leão, 2005). The main goal is to measure the parameters of the soil that contributes to the modeling of groundwater flow and erosion processes related ('piping' and Gully erosion) as has been approached by the Laboratory of Geo-Hydroecology of UFRJ. It is intended to create benefits for the understanding of groundwater flow modeling, using the knowledge of the hydrological behavior of soils by different methods of analysis, comparing these results with those found in tests of hydraulic conductivity made in the field. Soil samples were collected during surveys for installation of equipment for monitoring water and grain size were analyzed by the ABNT (technical standard NBR 7181/84) method. Some samples were analyzed by another method for flocular size, which includes a modified method without using deflocculant in the process of sedimentation, as described by Moraes Silva (2000). In the field were measured values of hydraulic conductivity through the "Slug test" injection (Hvorslev, 1951), through the previously installed piezometer. In the modeling of hydrological processes, the parameters for the soil are decisive in the analysis generated by the model, so it is necessary to consider not only the grain size, but also the structure that it has in the field, as well as the values of hydraulic conductivity related to points of monitoring. The results of granulometric analysis show that there is significant variation between the soils of headwater drainage. Furthermore, the results of tests of hydraulic conductivity also showed to be very similar, with the scale of values between 10-2 to 10-4cm /s. KEYWORDS: Textural Analysis, Hydraulic Conductivity and Groundwater. INTRODUÇÃO
Algumas bacias que drenam a vertente direita (serra da Bocaina) do médio vale do rio Paraíba
do Sul vem expandindo suas respectivas redes de canais, especialmente nos vales de cabeceiras de
drenagem. Isto decorre do crescimento regressivo de canais incisos por ativação de mecanismos
erosivos associados a descargas críticas de água subterrânea em faces de exfiltração (Coelho Netto et
al., 1988). Este tipo de canal erosivo, resultante da ação de água subterrânea, foi denominado por
Pichler (1953) de voçoroca. Na bacia do rio Piracema, principal tributário do rio do Bananal, afluente
à direita do rio Paraíba do Sul, Cambra (1998) verificou que 83% das voçorocas (n=117) estão
topograficamente ajustadas à rede regional de canais. Significa, portanto, que é no domínio das
cabeceiras de drenagem, aonde crescem estes canais de ordem hierárquica inferior, que ocorre a
principal integração dos processos de encostas com os processos fluviais (Coelho Netto, 2003).
A evolução da rede de canais, no entanto, não apresenta um comportamento espacialmente
uniforme, especialmente no dominio de colinas rebaixadas. Ou seja, os canais crescem e bifurcam
irregularmente no espaço e no tempo, de tal modo que hoje se observa uma paisagem aonde às
cabeceiras de drenagem apresentam-se em diferentes estágios evolutivos. Enquanto muitas cabeceiras
ainda estocam espessos pacotes de sedimentos, acumulados no ciclo de instabilidade erosivo-
deposicional entre 10.000 e 8.000 anos atrás (Dietrich et al., 1991; Coelho Netto et al., 1994), outras já
foram erodidas pelo avanço das voçorocas. Uma parte destas voçorocas está em processo ativo de
crescimento vertical e regressivo, outra parte já se encontra estabilizada por ajuste de declive com o
canal coletor adjacente (Gilbert, 1877) e/ou rebaixamento do nível d’água subterrâneo e/ou redução da
descarga do fluxo d’água subterrânea no fundo do canal (Rocha-Leão, 2005).
Avelar e Coelho Netto (1992a e 1992b) verificaram o paralelismo entre a orientação dos eixos
de vales de cabeceiras e a orientação das fraturas locais e, através de uma extensa rede de piezômetros
e poços, observaram a ocorrência de fluxos d´água artesianos ao longo do eixo de um vale de cabeceira
suspenso que acompanha o fraturamento subjacente. Os autores associaram a exfiltração destes fluxos
artesianos à iniciação de canais e origem dos vales de cabeceiras de drenagem. Mais recentemente,
Fonseca et al. (2006) confirmaram a ocorrência de artesianismo no fundo de um vale de cabeceira
ajustado topograficamente ao rio Bananal, e sob voçorocamento ativo. Todos estes autores detectaram
um atraso da ordem de 2 a 4 meses nas resposta das cargas de pressão em relação às chuvas
precipitadas sobre a bacia do rio Bananal, confirmando que estes fluxos artesianos em fraturas não são
regulados por aqüífero local e sim regional.
Alguns trabalhos relevantes discutem a relação entre o crescimento de canais e suas respectivas
áreas de contribuição e gradiente topográfico. Os trabalhos clássicos de Horton (1945) e Dunne (1980)
apontam que os canais crescem regressivamente até uma distância crítica do divisor, a partir de onde a
concentração dos fluxos d´água da área de contribuição não atingiria uma descarga crítica, limitando
trabalho erosivo. Montgomery e Dietrich (1990) ponderam que o tamanho da área de contribuição
varia inversamente com o gradiente topográfico. No entanto, no ambiente de colinas da bacia do rio
Piracema, o desenvolvimento de alguns canais nem sempre é controlado pelo tamanho e/ou gradiente
topográfico das suas respectivas áreas de contribuição.
O presente trabalho focaliza a caracterização textural dos materiais do regolito, e direciona
atenção especial para as amostras mais profundas para relacioná-las com medidas in situ de
condutividade hidráulica, como suporte aos estudos voltados a modelagem da “Rede de fluxos de água
subterrânea entre cabeceiras de vales adjacentes e implicações no crescimento da rede de canais:
médio vale do rio Paraíba do Sul”, tese de doutorado de Paulo V. Leal.
MATERIAL E MÉTODO
Área de estudo: a bacia do rio Piracema
A área de estudo (~4 km2) está localizada na fazenda Bela Vista, a 12km da cidade de Bananal
(SP),e, é formada pelo anfiteatro da Bela Vista (0,25 km2) e pelas sub-bacias com relação de
vizinhança direta. Toda a drenagem da área converge para o rio Piracema, que possui 132 km2 de área
de bacia, e é o principal afluente do rio do Bananal (518km2 de área de bacia).
O Anfiteatro da Bela Vista é formado por uma concavidade principal, a qual é alimentada por
concavidades menores. O fundo da concavidade principal está preenchido por um pacote aluvio-
coluvionar espesso, que vem sendo dissecado pelo avanço remontante da voçoroca da Bela Vista, que
tem aproximadamente 450 metros de comprimento, 40 metros de largura média e profundidades
superiores a 15 metros. A localização regional da área de estudo pode ser vista na figura 1A.
Na figura 2 estão destacados quatro vales de cabeceira de drenagem, vizinhos ao Anfiteatro da
Bela Vista (V1). A cabeceira V2 está topograficamente suspensa em relação ao Anfiteatro da Bela
Vista, possui fundo plano e não é canalizada. As cabeceiras V3, V4 e V5 estão rebaixadas em relação
ao Anfiteatro da Bela Vista, sendo que a cabeceira V3 possui fundo aplainado e divisor topográfico
comum à cabeceira V4 muito rebaixado.
A bacia do rio do Bananal apresenta duas unidades fisiográficas distintas, que podem ser
observadas também na bacia do rio Piracema (figura 1B): uma, o compartimento montanhoso da Serra
da Bocaina, com desnivelamentos topográficos acima dos 600 metros, encostas íngremes e fundos de
vale estreitos; outra, o compartimento de colinas, com fundos de vales largos e ligeiramente inclinados
com desnivelamentos abaixo dos 500 metros (Fernandes, 1990). A área de estudo (figura 2) está
localizada no compartimento colinoso.
A área de estudo drena um substrato geológico de rochas gnáissicas da unidade
metassedimentar São João (Sillimanita-granada-muscovita-biotita gnaisse com intercalações de níveis
ou lentes de rochas calciossilicáticas, gondito, mármores e de sillimanita-muscovita-biotita xisto)
descritas por Almeida et al. (1991 e 1993).
Figura 1 – A- Localização regional da bacia do rio do Bananal (em laranja) e do rio Piracema
Análise Textural e Condutividade Hidráulica
As sondagens foram feitas com um trado mecânico e as amostras foram colhidas nas
profundidades que apresentaram mudança textural. As análises granulométricas foram realizadas no
laboratório de caracterização da COPPE-Geotecnia, através da metodologia da ABNT (NBR 7181/84).
Foram ainda realizadas análises através de uma metodologia adaptada sem o uso de defloculante e
aparelho dispersor (Rodriguez, 2005).
(Fotografias referentes aos ensaios conduzidos em laboratório)
Figura 2 - Mapa da área de estudo, mostrando a localização dos piezômetros (os destacados em amarelo são os que foram usados nos testes de condutividade hidráulica).
Os testes de condutividade hidráulica do tipo “slug test” de injeção, foram conduzidos em
campo, seguindo a metodologia do Hvorslev (1951), em piezômetros previamente instalados. Os testes
consistem na medição do nível de água dentro do aparelho, na injeção de água limpa nos piezômetros e
na observação do tempo que o nível da água dentro do piezômetro leva para descer. Para facilitar a
leitura do nível da água dentro do piezômetro, um tubo de PVC do mesmo diâmetro do piezômetro,
com um tubo transparente acoplado e uma fita métrica foi usado (Figura 3). Para o cálculo da
condutividade hidráulica (K), foi utilizada a seguinte fórmula:
Figura 3 – Ensaios de condutividade hidráulica (slug test – injeção) realizados em campo.
As amostras mais profundas de cada sondagem foram selecionadas para serem analisadas
segundo as duas metodologias (com e sem defloculante e aparelho dispersor), por serem nessas
profundidades que os piezômetros foram instalados. Portanto, quando se faz o “slug test” é a camada
de solo na profundidade da ponta do piezômetro que se analisa. Os resultados foram então comparados
com os valores de condutividade estabelecidos para cada tipo de material (tabela1).
Tabela 1 – Valores de condutividade hidráulica (Freeze & Cherry, 1979)
Sendo:
Qc – a vazão da condição de fluxo constante
F – Fator forma → 2∏L/ln [(L/D)+1+(L/D)2]
Hc – Carga constante
K=Qc/ (F x Hc)
Material K (cm/s) Cascalho 10-1 a 100 Areia Lavada 10-4 a 1 Areia siltosa 10-5 a 10-1 Silte 10-7 a 10-3 Argila 10-10 a 10-6
RESULTADOS E DISCUSSÕES
D1 - domínio de fundo de vale (alúvio-coluvio).
D2 - domínio dos divisores de drenagem
D3 - domínio de encostas - concavidades superiores (colúvios)
De acordo com a caracterização textural do regolito nas diferentes estações piezométricas,
foram estabelecidos diferentes domínios de acordo com as características
Com os resultados da análise granulométrica foram desenvolvidos gráficos em ambiente Excel
para uma melhor análise dos mesmos, que posteriormente foram agrupados em diferentes domínios de
acordo com suas características texturais.
Domínio dos divisores de drenagem: camadas alternadas mais arenosas ou mais sílticas, com
menores teores de argila em E33 e E24 (saprolito no biotita-gnaisse da Unidade São João);
Domínio de fundo de vale (alúvio-coluvio): predomínio de sedimentos finos no metro superior
do solo; em profundidade prevalece a fração areia (alúvio-colúvio, incoeso);
Domínio de encostas - concavidades superiores (colúvios): camadas alternadas no perfil; tende
a apresentar maiores teores de argila nos perfis.
0% 20% 40% 60% 80% 100%
3,7-3,9
6,6-6,8
15,8-16Pro
fun
did
ade
(m)
E 33
P edregulho
Areia
S ilte
Argila
0% 20% 40% 60% 80% 100%
0-0,30,4-0,61,4-1,61,9-2,1
3,8-45,8-6
6,4-6,77,2-7,38,1-8,3
Pro
fun
did
ade
(m
)
E 34
P edregulho
A reia
S ilte
A rgila
0% 20% 40% 60% 80%100%
3.0-3.45
5-5.45
6-6.45
7-7.45
8-8.45
9-9.45
Pro
fun
did
ade
(m)
E 24
P edregulho
A reia
S ilte
A rgila
0% 20% 40% 60% 80% 100%
0,3-0,6
0.6-0.9
1-1.4
2.7-3.0Pro
fun
did
ade
(m)
E 30
P edregulho
Areia
S ilte
Argila
0% 20% 40% 60% 80% 100%
0-0,3
0,3-0,5
1-1,2
5,5-5,7
9,3-9,5
11,5-11,8
Pro
fun
did
ade
(m)
E 35
P edregulho
Areia
S ilte
Argila
0% 20% 40% 60% 80% 100%
0-0,31,3-1,62,4-2,73,1-3,44,4-4,77,4-7,68,4-8,69,5-9,7
11-11,213-13,5
16,4-16,7
Pro
fun
did
ade
(m)
E 27
P edregulho
Areia
S ilte
Argila
0% 20% 40% 60% 80% 100%
0-0,30,3-0,6
0,9-1,202,3-2,5
2,7-33-3,3
4,7-4,96,8-78-8,28,8-9
9,8-1010,6-
Pro
fun
did
ade
(m)
E 29
P edregulho
Areia
S ilte
Argila
0% 20% 40% 60% 80% 100%
0-0,30,6-0,71,6-1,9
3,8-45,9-6,2
7-7,48,4-8,69,4-9,6
12,6-18-18,2
18,8-
Pro
fun
did
ade
(m)
E 26
P edregulho
Areia
S ilte
Argila
0% 20% 40% 60% 80% 100%
0-0,3
1,4-1,6
9,8-10
11,7511,95
Pro
fun
did
ade
(m)
E 36
P edregulho
Areia
S ilte
Argila
0% 20% 40% 60% 80% 100%
0-0,31,3-1,62,4-2,73,1-3,44,4-4,77,4-7,68,4-8,69,5-9,7
11-11,213-13,5
16,4-16,7
Pro
fun
did
ade
(m)
E 27
P edregulho
Areia
S ilte
Argila
0% 20% 40% 60% 80% 100%
0-0,30,3-0,6
0,9-1,202,3-2,5
2,7-33-3,3
4,7-4,96,8-78-8,28,8-9
9,8-1010,6-
Pro
fun
did
ade
(m)
E 29
P edregulho
Areia
S ilte
Argila
0% 20% 40% 60% 80% 100%
0-0,30,6-0,71,6-1,9
3,8-45,9-6,2
7-7,48,4-8,69,4-9,6
12,6-18-18,2
18,8-
Pro
fun
did
ade
(m)
E 26
P edregulho
Areia
S ilte
Argila
0% 20% 40% 60% 80% 100%
0-0,3
1,4-1,6
9,8-10
11,7511,95
Pro
fun
did
ade
(m)
E 36
P edregulho
Areia
S ilte
Argila
Referente às amostras mais profundas, foi feita a classificação destas segundo o triângulo
textural proposto por Trefethen (1950) como pode ser observado abaixo:
Nesta classificação as amostras correspondem às profundidades em que foram realizados os
ensaios de campo para avaliar a condutividade hidráulica saturada. Predominando materiais areno-
siltosos e areno-argilo-siltosos.
Os resultados dos testes de condutividade mostraram valores muito próximos para as estações
E29, E30, E34, E35 e E36 (Ksat=10-4). As estações E24 e E38 apresentaram valores de condutividade
maiores (Ksat=10-2) onde a figura 4 permite uma localização espacial das estações onde foram
conduzidos os testes de condutividade hidráulica.
100 80 60 40 20 0Areia (%)
0
20
40
60
80
100
Silte (%)
0
20
40
60
80
100
Argi
la (%
) Areia
Silte
Argila
Legenda:E30 Cd
E30 Sd
E29 Cd
E29 Sd
E26 Cd
E26 Sd
E28 Cd
E28 Sd
E27 Cd
E27 Sd
E34 Cd
E34 Sd
E36 Cd
E36 Sd
E35 Cd
E35 Sd
E33 Cd
E33 Sd
E2 Cd
E2 Sd
E4 Cd
E4 Sd
E6 Cd
E6 Sd
E11 Cd
E11 Sd
E12 Cd
E12 Sd
Classificação das amostras segundo o triângulo textural proposto por
Trefethen (1950).
100 80 60 40 20 0Areia (%)
0
20
40
60
80
100
Silte (%)
0
20
40
60
80
100
Argi
la (%
) Areia
Silte
Argila
Legenda:E30 Cd
E30 Sd
E29 Cd
E29 Sd
E26 Cd
E26 Sd
E28 Cd
E28 Sd
E27 Cd
E27 Sd
E34 Cd
E34 Sd
E36 Cd
E36 Sd
E35 Cd
E35 Sd
E33 Cd
E33 Sd
E2 Cd
E2 Sd
E4 Cd
E4 Sd
E6 Cd
E6 Sd
E11 Cd
E11 Sd
E12 Cd
E12 Sd
100 80 60 40 20 0Areia (%)
0
20
40
60
80
100
Silte (%)
0
20
40
60
80
100
Argi
la (%
) Areia
Silte
Argila
100 80 60 40 20 0Areia (%)
0
20
40
60
80
100
Silte (%)
0
20
40
60
80
100
Argi
la (%
) Areia
Silte
Argila
Areia
Silte
Argila
Legenda:E30 Cd
E30 Sd
E29 Cd
E29 Sd
E26 Cd
E26 Sd
E28 Cd
E28 Sd
E27 Cd
E27 Sd
E34 Cd
E34 Sd
E36 Cd
E36 Sd
E35 Cd
E35 Sd
E33 Cd
E33 Sd
E2 Cd
E2 Sd
E4 Cd
E4 Sd
E6 Cd
E6 Sd
E11 Cd
E11 Sd
E12 Cd
E12 Sd
Classificação das amostras segundo o triângulo textural proposto por
Trefethen (1950).
Figura 4 – MDT com as estações piezométricas
CONCLUSÕES
Ficou constatado que a porção argilosa desses solos encontra-se totalmente, agregada entre si,
ou com partículas siltosas e arenosas. Deve-se chamar a atenção para o fato de a maioria das amostras
selecionadas serem arenosas, e por esse motivo não mudarem de classe, quando o método de análise
foi alterado. No entanto, pode-se perceber que mesmo quando a amostra se apresentou argilosa (E29),
a metodologia modificada mostrou que o conteúdo de argila estava todo agregado, principalmente à
fração arenosa.
Os valores de condutividade hidráulica encontrados, mostram que a caracterização
granulométrica dos solos feita com a metodologia modificada (sem uso de defloculante e aparelho
disperssor) expressa melhor o comportamento hidrológico dos pacotes sedimentares, já que os valores
encontrados tem condutividade de areias lavadas a areias siltosas.
Futuramente, ensaios de condutividade serão realizados nos outros piezômetros, para entender
o comportamento hidrológico dos solos presentes na área de estudo, pois o entendimento dessas
características é fundamental para a construção de um modelo conceitual que possibilite uma
modelagem hidrológica que represente a realidade de forma satisfatória.
REFERÊNCIAS BIBLIOGRÁFICAS
ALMEIDA, J.C.H.; EIRADO SILVA, L.G.A. & AVELAR, A.S. (1991). Coluna tectono-
estratigráfica de parte do Complexo Paraíba do Sul na Região de Bananal - S.P. Atas do 2o Simp.
Geol. Sudeste, São Paulo. SBG. p.509-513.
ALMEIDA, J.C.H.; EIRADO SILVA, L.G.A. & VALLADARES, C.S. (1993). O Grupo Paraíba do
Sul e as rochas granitóides na região de Bananal (SP) e Rio Claro (RJ): uma proposta de formalização
litoestratigráfica. Atas do 3o Simp. Geol. Sudeste, São Paulo. SBG.
AVELAR, A.S. & COELHO NETTO, A.L. (1992a). Fraturas e desenvolvimento de unidades
geomorfológicas côncavas no médio vale do rio Paraíba do Sul. Rev. Bras. de Geociências.
AVELAR, A.S. & COELHO NETTO, A.L. (1992b). Fluxos d’água subsuperficiais associados a
origem das formas côncavas do relevo; Anais da 1a Conferência Brasileira de Estabilidade de Encostas
/ COBRAE, ABMS e SBGE, Rio de Janeiro; vol. 2:709-719p.
CAMBRA, M.F.E.S. (1998). Movimentos de água na porção superior de solos sob pastagem: o papel
do sistema radicular. Tese de Mestrado, IGEO/UFRJ, Rio de Janeiro, 144p.
COELHO NETTO, A.L, FERNANDES, N.F. & DEUS, C.E. (1988). Gullying in the Southeastern
Brazilian Plateau: Bananal, SP. Sediment Budgets (Proceedings of the Porto Alegre Symposium,
Decemmber 1988). IAHS Publ. nº 174, p.35-42.
DIETRICH. W.E.; MONTGOMERY. D.; COELHO NETTO, A.L. & MOURA, J.R.S. (1991).
Evidence for regional aggradation starting in the Early Holocene in southeastern Brazil and for
degradation due to deforestation. Am. Geoph.Union, Fall Meet.,San
Francisco,USA/EOS,vol.70(43):1124.
DUNNE, T. (1980). Formation and controls of channel networks. Prog. Phis. Geogr.4:221-239.
Fernandes, N.F. (1990). Hidrologia subsuperficial e propridades físico-mecânicas dos Complexos de
Rampa - Bananal (S.P.). Tese de Mestrado, IGEO/UFRJ, Rio de Janeiro 151p.
GILBERT, K.G. (1877). Report on the geology of the Henry Mountains (Utah). U.S.Geogr. Geol.
Survey Rocky Mtn. Region (Powell): 18-98.
FREEZE, R. A. & CHERRY, J. A. (1979). Groundwater. Englewood Cliffs, New Jersey, Prentice
Hall Inc., 604p.
HORTON, R.E. (1945). Erosional development of streams and their drainage basins; hidrophysical
approach to quantitative morphology. Geol. Soc. Am. Bull.v.56, 144-152, 1945.
HVORSLEV, M. J. (1951). Time lag and soil permeability in groundwater observations. Bulletin No.
36, U.S. Army Engineer Waterways Experiment Station, Vicksburg, MS.
LEAL, P.J.V. (2004). Relação entre fluxos subterrâneos de vales vizinhos e a expansão da rede de
drenagem via voçorocamentos na bacia do rio Piracema, Bananal,SP. Dissertação de Mestrado –
Universidade Federal do Rio de Janeiro/PPGG.
MONTGOMERY, D.R. & DIETRICH, W.E. (1989). Source areas, drainage density and channel
initiation. Water Resources Research, 26:1907-1918.
PICHLER, E. (1953). Boçoroca. Boletim de Soc. Brasileira de Geologia, V.2, n.1, pp. 3- 16.
ROCHA LEÃO, O.M. (2005). Evolução regressiva da rede de canais por fluxos de água subterrânea
em cabeceiras de drenagem: bases geo-hidroecológicas para recuperação de áreas degradadas com
controle de erosão. Tese de doutorado – Rio de Janeiro: UFRJ. 247 p.
RODRIGUEZ, T. T. (2005). Proposta de Classificação Geotécnica para Colúvios Brasileiros. Tese de
doutorado. Prog. Pós-Grad. Eng. Civil – COPPE/Geotecnia - UFRJ.
TREFETHEN, J. M. (1950). Classification of sediments. American Journal of Science, 248, 55–62.