Revista Brasileira de Geografia Física, vol 07, n.03, (2014), 476-485.

Revista Brasileira de
Geografia Física
ISSN:1984-2295
Homepage: www.ufpe.br/rbgfe

Potencial de Retenção de Águas Pluviais pelo Método “Curve Number”¹

Marcia Regina Gomes de Jesus Soares², Jorge Antonio Silva Centeno³, Chisato Oka Fiori3,
Claudinei Taborda da Silveira3, Eloy Kaviski3
1
Parte da Tese de Doutorado da primeira autora.
²Doutoranda em Geografia – Universidade Federal do Paraná, UFPR e Bolsista da CAPES. CEP 81531-980, Curitiba, PR. E-mail:
marcia.rgsoares@gmail.com - autor correspondente.
³Professores doutores da Universidade Federal do Paraná – UFPR.

Artigo recebido em 12/05/2014 e aceite em 20/08/2014.

RESUMO
No presente trabalho são apresentados os resultados de um estudo cujo objetivo geral é analisar a evolução temporal do
potencial de retenção de águas pluviais de bacias hidrográficas com base no parâmetro “Curve Number” (CN) do
método de estimativa de vazões do Soil Conservation Service (SCS), visando contribuir para a gestão e manejo dos
recursos hídricos. O método do SCS é um modelo chuva-vazão que é utilizado para estimar vazões em pequenas bacias
com base em dados pedológicos e de evolução de uso e ocupação do solo. A cobertura do solo, derivada de imagens de
satélite Landsat TM 5 foi usada para estimar o parâmetro CN e a evolução temporal deste parâmetro foi analisada para
caracterizar as principais mudanças na área de estudo. O método foi aplicado na sub-bacia hidrográfica do rio Mourão
de 854 km2, localizada no curso médio do rio Ivaí – PR. Comparando os dados para o ano de 1986 e 2011, obtiveram-
se os seguintes resultados: 68% da área permaneceram com o mesmo valor de CN, seguido de 20% de aumento da
impermeabilização e 12% de diminuição da mesma. Observou-se que a principal causa para o aumento da
impermeabilização foi que a urbanização também aumentou em 7% em áreas de solos com restrições ao uso. Estes
resultados explicam o aumento da vazão de 2011 em comparação a 1986, que foi na média de 8%, aliada também ao
aumento da precipitação em 208,7 mm para o mesmo período.

Palavras-chave: Modelo de análise, bacia hidrográfica, impermeabilização do solo.

Potential for retention by method of rainwater “curve number”

ABSTRACT
In this paper we present the results of a study aimed at the analysis of the temporal the evolution of the potential for
stormwater retention in basins using the "Curve Number" (CN) parameter of the Soil Conservation Service (SCS), run-
off model. The SCS model is a run-off model that is used to estimate run-off in small watersheds using soil data and
information about land use and land cover. Land cover was derived from satellite imagery and allowed to compute the
CN parameter for each pixel. Comparing the values obtained for two different years, changes within the basin were
analyzed. The method was tested using data of a 854 km2 sub-basin in the middle course of the Mourão river. The
analysis lead yielded the following results: 68% of the area remained the same value CN, followed by a 20% increase
and 12% of waterproofing decreases. It was observed that the main cause for the increase in urbanization that the
sealing was also increased by 7% in areas with fragile soils. These results explain the increased flow rate of 2011
compared to 1986, which was on average 8%, combined with an increasing precipitation of 208.7 mm for the same
period.

Keywords: Temporal evolution, waterproofing, soil sealing.

Introdução hidrográfica. Estes elementos são importantes

O tipo e a cobertura do solo são duas para o estabelecimento de políticas de gestão
variáveis determinantes na capacidade de de uma bacia, visto que traduzem a evolução
infiltração e na retenção das águas do da ocupação de sua superfície e as atividades
escoamento superficial em uma bacia de preservação.

Soares, M; Centeno, J.; Oka Fiori, C.; Silveira, C.; Kaviski, E 476
 Revista Brasileira de Geografia Física, vol 07, n.03, (2014), 476-485.

Com base no exposto, o presente trabalho tem com os dados de vazão e pluviosidade com o
como objetivo geral analisar a evolução objetivo de conhecer até que ponto a dinâmica
temporal do potencial de retenção de águas multitemporal do uso influencia no
pluviais de bacias hidrográficas através do comportamento hidrológico.
estudo o uso e cobertura do solo derivado de A área de estudo do presente trabalho
imagens de satélite para a gestão e manejo dos compreende a sub-bacia hidrográfica do rio
recursos hídricos. Para isto, foi avaliada a Mourão no curso médio do rio Ivaí entre as
viabilidade do uso do parâmetro CN (curve coordenadas UTM 345.500,2m e
number) do modelo chuva-vazão 7.342.774,9m; 378.837,8m e 7.293.496,2m
desenvolvido pelo Serviço de Conservação do (Figura 1). Seu exutório coincide com o ponto
Solo (SCS), órgão do Departamento de da estação de coleta de vazão e pluviosidade.
Agricultura dos Estados Unidos da América Para atender os objetivos propostos foi
(USDA), que é muito utilizado para estimar construído um banco de dados com
vazões em pequenas bacias com base em informações derivadas de imagens Landsat e
dados que podem ser derivados de dados a respeito dos solos derivados de mapas
sensoriamento remoto e Sistemas de pedológicos na escala de 1:250.000, da
Informação Geográfica. EMBRAPA (1984). A área de estudo abrange
Um método simples e largamente utilizado os municípios de Campo Mourão, Luiziana e
em modelos hidrológicos de bacias Mamborê, no estado do Paraná.
hidrográficas no mundo, o método CN tem Material e Métodos
sido empregado por pesquisadores como A variação temporal da cobertura foi
Hawkins, 1993; Jacobs, 2005; Tedela et al., estimada usando imagens Landsat, de
2008; Santos et al, 2007; Banasik, 2010; Cao diferentes épocas. Através das imagens de
et al., 2011; Araújo Neto, 2012. O método do satélite Landsat TM 5 compreendendo um
SCS baseia-se no conceito de que a lâmina de período de multitemporalidade de vinte e
escoamento superficial produzida em um cinco anos (1986 e 2011), foi verificada a
dado evento é uma função da altura total da evolução do adensamento urbano e seu grau
lâmina precipitada e de perdas que ocorrem, de impacto na hidrologia. Para tanto foi
principalmente, devido à infiltração, à realizada a classificação supervisionada das
interceptação vegetal e à retenção em imagens multiespectrais, utilizando o método
depressões do terreno (TYAGI et al., 2008). da Máxima Verossimilhança Gaussiana,
Neste estudo, são levadas em consideração as dentro do software Spring do INPE (2013).
variáveis físicas pedológicas, O classificador utilizado na elaboração
geomorfológicas e de evolução de uso e da carta de uso e cobertura do solo foi o
ocupação do solo que foram confrontadas “máxima verossimilhança (maxver)” que

Soares, M; Centeno, J.; Oka Fiori, C.; Silveira, C.; Kaviski, E 477
 Revista Brasileira de Geografia Física, vol 07, n.03, (2014), 476-485.

considera a ponderação das distâncias entre obtidos nas cartas de uso e cobertura do solo
médias dos níveis digitais das classes, foram: Agricultura; Água/banhado;
utilizando parâmetros estatísticos. Os temas Vegetação; Área urbana.

Figura 1. Figura de localização geográfica da área de estudo.

O Soil Conservation Service (SCS) do adaptação e simplificação de um de seus

Departamento de Agricultura dos Estados parâmetros.
Unidos propôs um método simplificado para O método consiste em estimar a chuva
estimar o escoamento superficial em pequenas efetiva, a parcela do volume precipitado que
bacias rurais. Esta estimativa se baseia em forma o escoamento superficial, e com isto o
levantamento de dados como uso e cobertura escoamento superficial. A chuva efetiva é
do solo, precipitação e tipo de solo. O método dada por:
se popularizou em estudos hidrológicos
devido à sua fácil aplicação. Posteriormente, Q = ((P - 0.3S)2)/(P + 0.7S) (1)
maneiras de estimar os parâmetros de entrada Onde,
a partir de imagens de satélite foram Q = altura da lâmina do escoamento, em mm
propostos, como o estudo de RAGAN; P = precipitação, em mm
JACKSON (1980), o que introduziu uma

Soares, M; Centeno, J.; Oka Fiori, C.; Silveira, C.; Kaviski, E 478
 Revista Brasileira de Geografia Física, vol 07, n.03, (2014), 476-485.

S = capacidade de recarga máxima da bacia No estudo, para derivar o valor do

após uma precipitação antecedente de 5 dias, parâmetro CN foram analisadas as
em mm. associações entre as classes de cobertura e
A capacidade de recarga máxima está ocupação obtidas das imagens orbitais e o
associada às características físicas da bacia, grupo hidrológico do solo. Para esta
em termos de solos e cobertura vegetal. classificação hidrológica, foram utilizadas as
informações pedológicas na escala 1:250.000
S = (24500/CN) - 254 (2) (Figura 2) fornecidas pela EMBRAPA
(1984). Estes dados permitiram determinar a
Para seu cálculo, é necessário distribuição espacial da classificação
determinar o parâmetro NC, Curve Number, hidrológica dos solos, conforme mostra a
que é um valor tabelado que descreve a Figura 3.
combinação de tipo de solo, a umidade De acordo com a Figura 2 e a Figura
antecedente e o uso e cobertura na bacia. 3, 81% da área de estudo é constituída por
Como uma bacia não apresenta Latossolos (LVd13, LVd9, LVdf10, LVdf4,
características uniformes em termos de solo e LVdf6, LVef2) os quais constituíram o grupo
cobertura vegetal, o valor representativo da hidrológico B, por se tratar de solos arenosos
bacia pode ser obtido ponderando cada valor com pouca argila, proporcionando melhor
presente na bacia pela proporção de área infiltração a escoamento superficial; 5% de
associada a este valor: Argissolo (PVd4), o qual constitui o grupo
hidrológico A, por se tratar de solo arenoso,
profundo e bem drenado; 12% de Neossolo
(3)
Chernossólico (RRe11, RRe12) que
constituíram o grupo hidrológico D, por se
Onde
tratarem de solos com argilas pesadas, muito
CN = número da curva;
impermeáveis e; 0,3% Neossolo típico
A = área da bacia;
(RRe7) que constitui o grupo hidrológico C
CN = número da curva de uma parte da bacia;
por se tratar de solos mais argilosos que o
A = área associada a um número de curva na
grupo B, com baixa permeabilidade.
bacia.

Soares, M; Centeno, J.; Oka Fiori, C.; Silveira, C.; Kaviski, E 479
 Revista Brasileira de Geografia Física, vol 07, n.03, (2014), 476-485.

Figura 2. Mapa de solos. Legenda: Latossolos (LVd13, LVd9, LVdf10, LVdf4, LVdf6, LVef2);
Argissolo (PVd4); Neossolo Chernossólico (RRe11, RRe12); Neossolo típico (RRe7)

A tabela original do método do SCS é capacidade de se mapear a cobertura usando

ampla e apresenta uma grande quantidade de imagens orbitais.
combinações. No estudo de RAGAN; Com base nessa simplificação, e
JACKSON (1980), esta quantidade de considerando as classes presentes na região de
possibilidades foi restrita em função da estudo, uma tabela simplificada foi proposta,
conforme a Tabela 1:
Tabela 1. Número da curva CN do SCS simplificado.
Grupo Hidrológico do Solo
A B C D
Agricultura 62 71 78 81
Uso do solo Água/banhado 100 100 100 100
Vegetação 36 60 70 76
Áreas urbanas 77 85 90 92

Soares, M; Centeno, J.; Oka Fiori, C.; Silveira, C.; Kaviski, E 480
 Revista Brasileira de Geografia Física, vol 07, n.03, (2014), 476-485.

Figura 3. Classificação hidrológica do solo. Legenda: A = alta permeabilidade; B = permeabilidade

alta a moderada; C = permeabilidade moderada a baixa; D = baixa permeabilidade.

A curva de permanência (Figura 4) limite inferior de cada intervalo (ordenada) e

constitui importante fonte de apontamento do a correspondente frequência relativa
potencial de regularização hidrológica, e ela acumulada (abscissa).
foi construída por dispor as vazões diárias
observadas no período considerado em ordem Resultados
decrescente. Com a amplitude da variação das A Figura 4 mostra a Curva de
vazões, foram definidos os intervalos de Permanência da área de estudo, ou a
classe e colocados os intervalos em ordem capacidade de retenção de águas superficiais,
decrescente e verificar o número de eventos mostrando curvatura típica de bacias
ocorridos em cada intervalo à frequência hidrográficas menores (Figura 4), com valores
absoluta. Foi calculada a frequência relativa variando de 3 e 52 m3/s. Notam-se condições
(frequência absoluta / número de dados) para de escoamento superficial de maior equilíbrio
cada intervalo e acumuladas seguindo a na distribuição espacial.
ordem anterior. Plotou-se em um gráfico o

Soares, M; Centeno, J.; Oka Fiori, C.; Silveira, C.; Kaviski, E 481
 Revista Brasileira de Geografia Física, vol 07, n.03, (2014), 476-485.

Figura 4. Curva de Permanência da área de estudo. Valor médio mensal de 1985-2012.

A área total da bacia compreende a podem ser observadas neste período. O

854 km2. Esta área foi classificada no resultado da classificação está apresentado
levantamento do uso do solo em duas datas: nas Figuras 5a e 5b.
1986 e 2011, para verificar as alterações que

(a) (b)

Figura 5. Uso e ocupação do solo: (a) – 1986 e (b) – 2011.

A percentagem de área associada a Utilizando os fatores apresentados na

cada tipo de cobertura foi calculada com base Tabela 1, foi possível determinar os valores
nestes mapas temáticos e é mostrada na tabela do número da curva (CN) para cada pixel da
2. imagem, e com isto obter a distribuição

Soares, M; Centeno, J.; Oka Fiori, C.; Silveira, C.; Kaviski, E 482
 Revista Brasileira de Geografia Física, vol 07, n.03, (2014), 476-485.

espacial deste parâmetro na bacia nas duas enquanto que tonalidades claras
datas em questão. O resultado é mostrado na correspondem a altos valores do CN.
Figura 6. Tonalidades escuras estão
associadas a baixos valores do parâmetro,
Tabela 2. Evolução do uso e ocupação do solo.
Grupo Hidrológico do Solo
Uso do solo Área em 1986 (km2) Área em 2011 (km2)
Agricultura 667 585
Água/banhado 4 18
Vegetação 107 115
Áreas urbanas 76 136

Comparando as figuras das duas datas mais altos, refletindo a expansão da

é notado que na parte norte da bacia valores agricultura (Figura 6).
baixos do CN foram substituídos por valores

Figura 6. Mapa de retenção máxima do solo: (a) – 1986 e (b) – 2011.

A seguir, a evolução do CN foi bacia continua com o mesmo valor de CN. É

analisada comparando pixel a pixel os mapas também observado que 20% dos pixels
temáticos mostrados na figura 6. Esta análise traduzem o aumento da impermeabilização e
permite verificar que grande parte da bacia 12% a diminuição da mesma (Figura 7).
permaneceu inalterada; em torno de 68% da

Soares, M; Centeno, J.; Oka Fiori, C.; Silveira, C.; Kaviski, E 483
 Revista Brasileira de Geografia Física, vol 07, n.03, (2014), 476-485.

Figura 7. Mudança de permeabilidade de 1986 a 2011.

Discussão ocupação do solo através da análise

Observou-se que a principal causa multitemporal. O método permitiu contemplar
para o aumento da impermeabilização foi que variáveis importantes para a análise e gestão
a urbanização também aumentou em 7% em de elementos físicos da paisagem e seus
áreas com solos frágeis ao uso agrícola e processos, ao se trabalhar com os dados de
ocupações urbanas e de textura que favorecem forma conjunta, o que é necessário para se
o escoamento superficial à infiltração. chegar a resultados completos e eficazes.
Embora o uso agrícola tenha diminuído em Desta forma o objetivo da pesquisa
10%, o espaço foi preenchido principalmente proporcionou a integração dos condicionantes
pelas áreas urbanas, resultando no aumento da do meio físico, apontando áreas mais ou
curva de escoamento, ou no agravamento dos menos suscetíveis às formas de modificações
valores do CN. estabelecidas pelas atividades humanas.
Assim se verifica contraposição do uso
do solo às características naturais da bacia Agradecimentos
hidrográfica. Estes resultados explicam o Ao INPE, pelo fornecimento das
aumento da vazão de 2011 em comparação a imagens de satélites para o auxílio na
1986, que foi na média de 8%, aliada também elaboração do mapa de uso e ocupação do
ao aumento da precipitação em 208,7 mm. solo da área de estudo.

Conclusão Referências
O método do CN constitui boa Araújo Neto, J. R.; Palácio, H.; Andrade, E.;
Santos, J.; Pinheiro, E. 2012. Otimização do
ferramenta para identificação de áreas frágeis
número de curva (CN-SCS) para diferentes
e limitantes ao considerar a evolução de uso e manejos na região semiárida, Ceará, Brasil.

Soares, M; Centeno, J.; Oka Fiori, C.; Silveira, C.; Kaviski, E 484
 Revista Brasileira de Geografia Física, vol 07, n.03, (2014), 476-485.

Irriga, Botucatu, Edição Especial, p. 264 – Ragan, R M.; T J Jackson (1980). Runoff
279. Synthesis Using Landsat and SCS Model,
Journal of Hydrology, Divn., ASCE, Vol. 106
Banasik, K.; Woodward, D. 2010. Empirical (HYS5): 667-678.
determination of runoff curve number for a
small agricultural watershed in Poland. 2nd Santos. C. A. G.; Vajapeyam. R. M. S.;
Joint Federal Interagency Conference, Las Srinivasan, S. 2007. Análise das perdas de
Vegas, NV, June 27 - July 1. água e solo em diferentes coberturas
superficiais no semiárido da Paraíba. Revista
Cao, H.; Vervoort, R. W.; Dabney, S. M. OKARA: Geografia em debate, João Pessoa,
2011. Variation in curve numbers derived v. 1, n. 1, p. 1-152.
from plot runoff data for New South Wales
(Australia). Hydrological Processes, v. 25, p. Tedela, N., Mccutcheon, S.; Rasmussen, T.;
3774-3789. Tollner, W. 2008. Evaluation and
Improvements of the Curve Number Method
Embrapa. 1984. Levantamento de of Hydrological Analysis on Selected
reconhecimento dos solos do estado do Forested Watersheds of Georgia. Report
Paraná. SLNCS. submitted to Georgia Water Resources
Institute. Available at
Hawkins, R. H. 1993. Asymptotic http:/water.usgs.gov/wrri/07grants/progress/2
determination of curve numbers from data. 007GA143B .pdf. [Acessado em 01 Agosto
Journal of Irrigation and Drainage Division, 2013].
American Society of Civil Engineers, v. 119,
n. 2, p. 334-345. Tyagi, J. V.; Mishra, S. K.; Singh, R.; Singh,
V. P. 2008. SCS-CN based time-distributed
Inpe – Instituto Nacional de Pesquisas sediment yield model. Journal of Hydrology,
Espaciais. 2013. SPRING Realise 4.2. São v. 352, p. 388–403.
Paulo, conceitos cartográficos.

Jacobs, J. H.; Srinivasan, R. 2005. Effects of

curve number modification on runoff
estimation using WSR-88D rainfall data in
Texas watersheds. Journal of Soil and Water
Conservation, v. 60, n. 5, p. 274-279.

Soares, M; Centeno, J.; Oka Fiori, C.; Silveira, C.; Kaviski, E 485