Acessibilidade / Reportar erro

Modelagem do regime térmico de um reservatório tropical de abastecimento público, Juiz de Fora, MG, Brasil

Thermal regime modeling of the public supply tropical reservoir at Juiz de Fora, MG, Brazil

Resumo

As alterações de temperatura da água e a dinâmica de estratificação térmica de um ecossistema aquático influenciam diretamente os processos físicos, químicos e biológicos. Devido às perturbações em que esses ambientes estão expostos, tais como mudanças climáticas e usos múltiplos da água, os estudos por meio de modelagem são importantes para suportar tomadas de decisão emergenciais e planejamentos futuros. Neste contexto, este trabalho teve como objetivo principal avaliar o comportamento térmico do reservatório de São Pedro utilizado para abastecimento de água, localizado no município de Juiz de Fora, MG, por meio da modelagem matemática com o modelo IPH-ECO para avaliar a temperatura superficial da água, com base em dados meteorológicos e hidrológicos, representando as forçantes e as condições de contorno, respectivamente. A modelagem do regime térmico mostrou correspondência significativa entre os valores simulados e observados, caracterizando o reservatório como um ambiente homogêneo. As simulações indicaram que a contribuição dos afluentes determina o padrão de distribuição espaço-temporal da temperatura da água no reservatório, limitando ou regularizando a sua variabilidade da temperatura.

Palavras-chave:
hidrodinâmica; modelagem; temperatura da água.

Abstract

Water temperature changes and the dynamics of thermal stratification of an aquatic ecosystem directly influence physical, chemical and biological processes. Due to disturbances in these environments, such as climate change and multiple water uses, studies which employ modeling are important to support emergency decision-making and future planning. In this context, this study evaluated the thermal behavior of the water supply reservoir of São Pedro, located in the city of Juiz de Fora, MG, using mathematical modeling with the IPH-ECO model to evaluate the surface temperature of water, based on meteorological and hydrological data representing the forcing and boundary conditions, respectively. The modeling of the thermal system showed a significant correlation between predicted and observed values, characterizing the reservoir as a homogeneous environment. The simulations indicated that the tributary contribution determines the pattern of spatio-temporal distribution of water temperature in the reservoir, limiting or regulating the variability of its temperature.

Keywords:
hydrodynamic; modeling; water temperature.

1. INTRODUÇÃO

A construção de reservatórios pode gerar impactos ambientais (i.e. eutrofização, sedimentação, toxicidade e veiculação de doenças), devido à transformação de um ambiente lótico em lêntico. Essas modificações nas condições de escoamento podem afetar, também, o regime térmico do sistema hídrico, além de introduzir novas condições químicas e biológicas (Ferreira e Cunha, 2013FERREIRA, D. M.; CUNHA, C. Simulação numérica do comportamento térmico do reservatório do Rio Verde. Revista Engenharia Sanitária e Ambiental, v. 18, n. 1, jan./mar. 2013. ).

Os reservatórios são construídos em pontos estratégicos de uma bacia hidrográfica com o objetivo de atender diferentes demandas, tais como geração de energia, abastecimento, pesca, irrigação, recreação e aqüicultura (Christofolleti, 1980CHRISTOFOLLETI, A. Geomorfologia. São Paulo: Blucher, 1980.). Além destas também são muito utilizados para a regularização da vazão e amortecimento de cheias, apresentando-se como uma estrutura hidráulica de múltiplas finalidades para a sociedade.

No geral, a variação da temperatura de um corpo hídrico é ditada pelo regime climático natural e, associada à hidrodinâmica, afeta praticamente todos os fluxos ecológicos, como produção primária, nitrificação e respiração (Fragoso Jr. et al., 2009FRAGOSO JR., C. R.; FERREIRA, T. F.; MOTTA MARQUES, D. M. L. Modelagem ecológica em ecossistemas aquáticos. 1. ed. São Paulo: Oficina de Textos, 2009. v. 1. 304p.). Portanto, a avaliação da dinâmica do regime térmico auxilia na compreensão da variabilidade desses processos nos ambientes aquáticos.

Neste contexto, os modelos matemáticos se apresentam como ferramentas que ajudam na obtenção de uma visão dinâmica dos processos naturais, pois permitem a integração de dados dispersos espacialmente, com a possibilidade de prever diferentes situações ao simular cenários futuros (Rosman et al., 2001ROSMAN, P. C. C.; MASCARENHAS, F. C. B.; MIGUEZ, M. G.; CAMPOS, R. O.; EIGER, S. Métodos numéricos em recursos hídricos. 5. ed. Rio de Janeiro: ABRH, 2001.).

Conhecer a variação da temperatura de um corpo de água é particularmente importante por três motivos: (a) as descargas de efluentes em diferentes temperaturas podem causar efeitos negativos no ecossistema aquático; (b) a temperatura influencia as reações químicas e biológicas; e (c) a variação da temperatura afeta a densidade da água, e como consequência, altera os processos hidrodinâmicos e de transporte de massa (COPPE/UFRJ, 2007INSTITUTO ALBERTO LUIZ COIMBRA DE PÓS-GRADUAÇÃO E PESQUISA DE ENGENHARIA - COPPE/UFRJ. Manual do Usuário do SisBaHia. 2007. Disponível em: http://www.sisbahia.coppe.ufrj.br/Arquivos_Gerais. Acesso em: jul. 2013.
http://www.sisbahia.coppe.ufrj.br/Arquiv...
).

As condições de temperatura de um reservatório podem ser diferenciadas ao longo da coluna de água, por exemplo, promovendo um gradiente de densidade que forma uma barreira física que impede que compartimentos verticais se misturem (i.e. superfície e fundo). Essa barreira física cria zonas de circulação na coluna de água, amplificando o efeito da estratificação térmica. A ausência de variação térmica vertical, dentre outros fatores, facilita a homogeneização do corpo de água.

A compreensão do regime térmico em resposta ao clima e às condições hidrológicas é fundamental para o desenvolvimento de estratégias para gestão adaptativa da qualidade da água em reservatórios (Wang et al., 2012WANG, S.; QIAN, X.; HAN, B.; LUO, L.; HAMILTON, D. P. Effects of local climate and hydrological conditions on the thermal regime of a reservoir at Tropic of Cancer, in southern China. Water Research, v. 46, n. 8, p. 2591-2604, 2012. http://dx.doi.org/10.1016/j.watres.2012.02.014
http://dx.doi.org/10.1016/j.watres.2012....
).

A modelagem do regime térmico de um reservatório depende do conhecimento das condições de contorno climáticas, que são as forçantes de temperatura do ar, radiação solar, evaporação, direção e intensidade do vento, precipitação e umidade, além dos dados hidrológicos que são as vazões de entrada e saída do reservatório.

Muitos e variados modelos de ecossistemas aquáticos foram desenvolvidos e publicados durante as últimas quatro décadas, com o objetivo de representar a variabilidade dos processos ecológicos. Dentre os modelos de qualidade da água utilizados em reservatórios podemos citar o IPH-ECO, desenvolvido no Instituto de Pesquisas Hidráulicas da Universidade Federal do Rio Grande do Sul, o SisBaHiA (Sistema Base de Hidrodinâmica Ambiental) desenvolvido na COPPE da Universidade Federal do Rio de Janeiro e o ELCOM (Estuary and Lake Computer Model) acoplado ao CAEDYM (Computacional Aquatic Ecosystem Dinamics Model) ambos desenvolvidos pelo CWR (Center of Water Research) da Universidade da Austrália.

No entanto, apesar de algumas tentativas terem sido feitas, trabalhos que tratam de abordagens integradas no espaço e no tempo em ecossistemas lacustres ainda são escassos (Mooij et al., 2010MOOIJ, W. M.; TROLLE, D.; JEPPESEN, E.; ARHONDITSIS, G.; BELOLIPETSKY, P. V.; CHITAMWEBWA, D. B. R. et al. Challenges and opportunities for integrating lake ecosystem modeling approaches. Aquatic Ecology, v. 44, 2010. http://dx.doi.org/10.1007/s10452-010-9339-3
http://dx.doi.org/10.1007/s10452-010-933...
). Além disso, a complexidade de fatores envolvidos (i.e. poluição, usos múltiplos da água, espécies invasoras, alterações do uso e cobertura do solo na bacia) dificulta a análise e quantificação dos processos.

Assim, este trabalho tem como objetivo investigar o regime térmico do reservatório de São Pedro, localizado no município de Juiz de Fora (MG), por meio da modelagem hidrodinâmica, em duas dimensões na horizontal (2DH), com base nos dados meteorológicos e hidrológicos, utilizando o modelo ecológico desenvolvido no Instituto de Pesquisas Hidráulicas - IPH da Universidade Federal do Rio Grande do Sul, denominado IPH-ECO.

2. MATERIAL E MÉTODOS

2.1. Área de estudo

O Reservatório de São Pedro, também conhecido como Represa dos Ingleses ou Cruzeiro de Santo Antônio, localizada na Mesorregião da Zona da Mata Mineira/Microrregião de Juiz de Fora, está situada nas coordenadas 21° 46' 34,8''S/43° 24' 27,2''W, altitude média de 850 m, possuindo 0,04 km2 de espelho de água próximo a mancha urbana do município de Juiz de Fora/MG, região oeste da cidade. Atualmente, consiste em um manancial urbano inserido em uma área apontada como vetor de crescimento pelo Plano Diretor do município (Juiz de Fora, 1999JUIZ DE FORA. Prefeitura Municipal. Plano diretor de desenvolvimento urbano de Juiz de Fora. Juiz de Fora: FUNALFA, 1999.).

Criado em 1963, o reservatório teve início da sua operação somente em 1967, com o objetivo de abastecer parte da demanda do município. A produção de água da Estação de Tratamento de Água (ETA São Pedro), localizada no bairro São Pedro, é em média de 120 L.s-1, variando de acordo com o nível do reservatório. Atualmente, este reservatório ainda constitui um importante manancial de abastecimento público, operado e mantido pela Companhia de Saneamento Municipal - CESAMA. Este manancial fica a 8 km da densa malha urbana e corresponde a 8% do abastecimento da cidade (CESAMA, 2013COMPANHIA DE SANEAMENTO MUNICIPAL DE JUIZ DE FORA - CESAMA. Website. Disponível em: http://www.cesama.br. Acesso: em abr. 2013.
http://www.cesama.br....
).

Sob o clima Cwb (tropical de altitude) - clima mesotérmico com secas de inverno e chuvas concentradas de verão e verões brandos ou moderadamente quentes, na classificação de Koppen (Machado, 2010MACHADO, P. J. O. (Org.) Diagnóstico físico-ambiental da bacia hidrográfica do Córrego São Pedro: um exercício acadêmico de gestão dos recursos hídricos. Ubá: Geographica, 2010.), a bacia hidrográfica do reservatório de São Pedro (BHRSP) ocupa uma área de 13,04 km2 e seus principais afluentes são o Córrego São Pedro (margem esquerda) e o Córrego Grota do Pinto (margem direita). A bacia do córrego São Pedro tem seu exutório no rio Paraibuna que deságua na margem direita do rio Paraíba do Sul no município de Três Rios - RJ. A localização da bacia hidrográfica do reservatório de São Pedro na área urbana de Juiz de Fora é mostrada na Figura 1.

2.2. O modelo IPH-ECO

O modelo ecológico IPH-ECO é um sistema computacional desenvolvido no Instituto de Pesquisas Hidráulicas (IPH) da Universidade Federal do Rio Grande do Sul (UFRGS), voltado especialmente para o entendimento dos processos físicos, químicos e biológicos de corpos de água rasos e profundos, tais como lagos, reservatórios e estuários (Fragoso Jr., 2009FRAGOSO JR., C. R. Modelagem tridimensional da estrutura trófica em ecossistemas aquáticos continentais. 2009. Tese (Doutorado em Recursos Hídricos e Saneamento Ambiental) - Instituto de Pesquisas Hidráulicas, Universidade Federal do Rio Grande do Sul, Porto Alegre, 2009.). O modelo descreve as mais importantes interações hidrodinâmicas e abióticas, além dos principais processos bióticos, com a finalidade de auxiliar no entendimento de um determinado ecossistema aquático e está disponível para download em https://sites.google.com/site/iphecomodel/.

Figura 1:
Localização da bacia afluente do reservatório de São Pedro.

Para este trabalho foi utilizado o módulo de balanço térmico associado ao módulo hidrodinâmico (bidimensional horizontal) do IPH-ECO para simulação da temperatura da água do reservatório de São Pedro. Uma vez que em corpos de água rasos, geralmente um modelo bidimensional na horizontal (2DH) é aconselhado, já que a dimensão horizontal predomina sobre vertical e nem sempre ocorre estratificação nestes sistemas (Rosman, 1999ROSMAM, P. C. C. Subsídios para modelagem de sistemas estuarinos. In: ROSMAN, P. C. C.; MASCARENHAS, F. C. B.; MIGUEZ, M. G.; CAMPOS, R. O.; EIGER, S. Métodos numéricos em recursos hídricos, Rio de Janeiro: ABRH 1999.).

A escolha do IPH-ECO se baseou na possibilidade de representação do regime térmico para um reservatório raso e acesso à informações fundamentais para utilização do modelo por meio de curso específico, já que existem diversos modelos com variadas adequações e possibilidades de representação. Além disso, é preciso considerar o acesso a informações de forma documental ou até mesmo informal como base para utilização do modelo.

2.3. Módulo hidrodinâmico

As equações governantes para o módulo hidrodinâmico no IPH-ECO são derivadas das equações de Navier-Stokes para um escoamento turbulento sob a suposição de distribuição hidrostática de pressões. Estas equações são baseadas na conservação da massa e quantidade de movimento e podem ser escritas na forma diferencial (Equações 1, 2 e 3).

Equação da Continuidade:

Equação da Quantidade de Movimento:

em que:

t é o tempo;

u e v são as componentes da velocidade nas direções x e y no plano horizontal, respectivamente;

é elevação da superfície da água medida a partir de um nível de referência (Figura 2);

h é a profundidade medida a partir de um nível de referência;

g é a aceleração da gravidade;

f é o parâmetro de Coriolis; e

são as tensões do vento nas direções x e y, respectivamente;

é um operador vetorial no plano x-y; e

Ah é o coeficiente de viscosidade turbulenta horizontal; e

é coeficiente de atrito junto ao fundo.

Figura 2
Elementos integrantes das equações de movimento.

O coeficiente de atrito junto ao fundo ( ) é estimado empiricamente (Equação 4).

em que:

H = h + é a profundidade total que corresponde à soma do nível de referência com o nível da superfície; e

Cz é o coeficiente de atrito de Chèzy.

Como é comum em modelos de escoamentos, a tensão de atrito na superfície livre pelo vento é escrita em termos da velocidade do vento (Rosman, 1999ROSMAM, P. C. C. Subsídios para modelagem de sistemas estuarinos. In: ROSMAN, P. C. C.; MASCARENHAS, F. C. B.; MIGUEZ, M. G.; CAMPOS, R. O.; EIGER, S. Métodos numéricos em recursos hídricos, Rio de Janeiro: ABRH 1999.; Equações 5 e 6).

em que:

CD é o coeficiente de arraste do vento (valores utilizados 1 e -6);

Wx e Wy são as componentes do vetor velocidade do vento nas direções x e y em m.s-1, respectivamente, medidas a 10 metros da superfície livre e

é o módulo do vetor velocidade do vento.

As equações de águas rasas (Equações de 1 a 3) não têm solução analítica genérica. Um esquema Euleriano-Lagragiano semi-implícito de diferenças finitas resolve numericamente as equações governantes para um número finito de pontos no espaço e no tempo (Casulli e Cheng, 1992CASULLI, V.; CHENG, R. Semi-implicit finite difference methods for three-dimensional shallow water flow. International Journal for Numerical Methods in Fluids, v. 15, p. 629-648, 1992. http://dx.doi.org/10.1002/fld.1650150602
http://dx.doi.org/10.1002/fld.1650150602...
). Este método necessita subdividir o domínio de aplicação em uma malha, comum número finito de células. As equações são discretizadas espacialmente em uma grade retangular que consiste em células computacionais quadradas com comprimento e largura .

2.4. Balanço de calor

No modelo IPH-ECO a temperatura da água é estimada a partir do balanço de calor entre a superfície da água e a atmosfera. Esse balanço é aplicação direta da primeira lei da termodinâmica para um fluido incompressível, a qual define que a variação da energia interna é igual à resultante dos fluxos de calor nos contornos do sistema. O balanço do fluxo de calor pode ser representado pela Equação 7.

em que:

/é o fluxo de calor referente à radiação atmosférica de onda curta,

/é o fluxo de calor referente à radiação atmosférica de onda longa,

/é o fluxo de calor referente à radiação da água de onda longa,

/é o fluxo de calor referente à condução para a atmosfera,

/é o fluxo de calor referente ao processo de evaporação e

/é o fluxo de calor resultante das entradas e saídas por advecção tais como escoamento superficial e subterrâneo, fluxo de afluentes, precipitação e fluxos de saída. O calor armazenado durante um intervalo de tempo pode ser convertido em temperatura (Equação 8).

em que:

V é o volume em m3,

é a densidade da água em kg m-3 e

é o calor específico da água em J.kg-1 °C-1 . A variabilidade temporal e espacial da temperatura da água pode ser avaliada para a equação de transporte de massa (Equação 9).

em que:

Kh é o coeficiente de difusão horizontal em m2.s-1,

T é a temperatura da água em graus Celsious (°C) e

St é o termo de perdas e ganhos (dado pela Equação 8). Esta equação é resolvida numericamente por meio de um esquema explícito centrado de diferenças finitas.

2.5. Amostragem e preparação dos dados para modelagem

2.5.1. Monitoramento da temperatura água

A amostragem neste trabalho foi realizada a partir da seleção de três pontos representativos ao longo do espelho de água do reservatório, denominados Ponto 1, Ponto 2 e Ponto 3 (Figura 3).

Figura 3:
Distribuição dos pontos amostrais no reservatório de São Pedro.

A seleção de pontos para o monitoramento deste manancial foi uma etapa importante para se obter um referencial de qualidade do mesmo. Foram selecionados 3 pontos, dois pontos próximos às margens e um ponto na área central do reservatório. Foram realizadas 10 campanhas de amostragem, entre os meses de julho de 2010 a junho de 2011, intencionalmente com frequência mensal, porém nos meses de novembro de 2010 e março de 2011 não foi possível realizar campanhas devido às chuvas e dificuldades de deslocamento para acesso às áreas de coleta. As amostragens in situ foram realizadas com a utilização de uma sonda modelo HQ30d da Hach, com a finalidade de registrar a temperatura da água nos diferentes pontos selecionados no reservatório (Tabela 1).

Tabela 1
Valores de temperatura da água observados in situ.

A temperatura da água foi monitorada próxima a superfície com cerca de 30 a 40 cm de profundidade controlado pela imersão do leitor da sonda utilizada na água. Estes dados foram comparados com os valores de temperatura estimados pelo IPH-ECO a partir dos dados climatológicos e dados das vazões de entrada e de saída do reservatório, os quais serão discutidos mais adiante.

Para a modelagem do regime térmico foi utilizado um intervalo de tempo de 10 segundos durante o período de simulação de 355 dias. Os principais dados de entrada do modelo são descritos a seguir.

2.5.2. Dados meteorológicos

Os dados meteorológicos foram fornecidos pela Estação Climatológica Principal situada na Universidade Federal de Juiz de Fora, operada e mantida pelo Instituto Nacional de Meteorologia na sede de Belo Horizonte/MG (INMET), são eles: temperatura do ar (ºC), evaporação (mm), direção e intensidade do vento (m.s-1), precipitação (mm) e umidade (%) com freqüência 00:00h/12:00h/18:00h. Os dados de radiação solar (w/m²) foram obtidos diretamente da Estação Climatológica Principal do INMET/MG com frequência de 1 hora, em sua página de internet. Esta estação dista cerca de 4,5 km do reservatório, o que garante pouca diferenciação das condições climáticas locais que poderiam afetar na comparação dos dados, como nebulosidade e precipitação.

2.5.3. Levantamento batimétrico

A batimetria do reservatório foi efetuada com a utilização de um Ecobatímetro (Garmin Fishfinder 140), aparelho que descreve a topografia do fundo do lago e os pontos foram marcados com um GPS (Garmin GPSMAP 76CSx) para garantir as amostragens nos mesmos pontos ao longo do monitoramento. Este levantamento foi realizado no dia 11 de julho de 2010, data referente ao primeiro monitoramento.

A batimetria revelou que a profundidade média do reservatório é de 2 metros (Figura 4), caracterizando-o como um reservatório raso, com predomínio da dimensão horizontal sobre a vertical, o que não justifica a aplicação de uma aproximação tridimensional. Desta forma, foi utilizada uma representação bidimensional na horizontal (2DH) do domínio. A discretização espacial consistiu em uma malha de 50x50 m resultando em 450 células computacionais.

Figura 4:
Batimetria do reservatório de São Pedro.

2.5.4. Estimativa das vazões afluentes

Uma estimativa das vazões de entrada do modelo foi obtida por meio do Método Racional Modificado (Euclydes, 1987EUCLYDES, H. P. Saneamento agrícola: atenuação de cheias: metodologia e projeto. Belo Horizonte: Rural Minas, 1987.) que resume todos os processos hidrológicos da bacia em um único coeficiente, o qual determina a parcela da chuva que se transforma em escoamento superficial (Equação 10).

em que:

Q é a vazão na seção de interesse (m3 s-1),

i a intensidade de chuva (mm h-1),

A é a área de drenagem na seção de interesse (km2),

o coeficiente de retardamento (adimensional), e

C é o coeficiente de escoamento superficial ou coeficiente de run off. Este valor é normalmente estimado em função do mapa de uso e ocupação do solo. Para bacia em estudo utilizou-se para os tipos de cobertura do solo os valores 0,20 para Mata, 0,80 para Área Urbana, 0,30 para Pasto/Solo Exposto e 0,0 para Água, conforme Tucci (2001TUCCI, C. E. M. Hidrologia: ciência e aplicação. Porto Alegre: Ed. UFRGS, 2001.).

O cálculo do coeficiente de retardamento é realizado com a Equação 11.

Para tanto é necessário o cálculo das áreas das bacias dos afluentes São Pedro e Grota do Pinto e seus respectivos coeficientes de retardamento.

O coeficiente de retardamento foi estimado com base em dados de vazão levantados in situ fazendo uso de um molinete hidrométrico nos afluentes São Pedro e Grota do Pinto. Os coeficientes de retardamento ajustados para os principais afluentes do reservatório foram: 0,2749944 para o São Pedro e 0,2768576 para o Grota do Pinto, cujas respectivas áreas de drenagem são 8,84 km2 e 3,36 km2.

A vazão de saída foi disponibilizada pela Estação de Tratamento de Água de São Pedro da Companhia de Saneamento Municipal - CESAMA que trata a água do reservatório de São Pedro e faz o controle horário dessa vazão para direcionar as operações dentro da estação.

2.5.5. Saída do modelo

Para a simulação da temperatura da massa de água do reservatório de São Pedro foram selecionados três pontos de saída do modelo (Pontos 1, 2 e 3), como mostra a Figura 5, com a finalidade de comparar os valores simulados com os valores observados em campo.

Figura 5:
Pontos de saída selecionados no modelo IPH-ECO.

2.5.6. Avaliação do desempenho do modelo

O grau de correspondência entre a saída do modelo e os dados observados foi medido com a utilização de dois índices estatísticos: (a) a raiz do erro médio quadrático (REMQ), e (b) o coeficiente de Nash-Sutcliffe (R2). O cálculo da raiz do erro quadrático médio e do coeficiente de determinação de Nash-Sutcliffe são apresentados nas Equações 12 e 13.

em que:

REMQ é raiz do erro médio quadrático,

Yobsé o valor observado,

Ycal é o valor estimado ou calculado pelo modelo,

N é o número total de valores,

R2 é o coeficiente de Nash-Sutcliffe e

/obs é a média dos valores observados.

3. RESULTADOS E DISCUSSÃO

3.1. Simulações da temperatura da água

Ao final da simulação, o modelo IPH-ECO estimou a temperatura da água para os pontos de interesse em um período total de 355 dias, iniciando às 00:00 horas do dia 01/07/2010 e finalizando às 23:00 horas do dia 01/06/2011. Esse período de simulação corresponde ao período de monitoramento realizado no reservatório permitindo observar o comportamento térmico nos períodos seco e chuvoso ao longo do ano em questão. A comparação entre as temperaturas da água simuladas e observadas entre os meses de julho de 2010 a junho de 2011 nos pontos 1, 2 e 3 é mostrada, respectivamente, na Figura 6.

Figura 6
Variação da temperatura da água de Julho de 2010 a Junho de 2011. a) Comparação entre os valores de temperatura da água simulados (linha contínua) e observados (pontos) no Ponto 1 (área central do reservatório). b) Comparação entre os valores de temperatura da água simulados (linha contínua) e observados (pontos) no Ponto 2 (margem esquerda do reservatório). c) Comparação entre os valores de temperatura da água simulados (linha contínua) e observados (pontos) no Ponto 3 (margem direita do reservatório).

Ao se analisar a Figura 6 pode-se perceber que os valores simulados e observados apresentam boa correspondência ao longo do tempo. Os indicadores da eficiência do ajuste apresentaram bons valores (REMQ<2% e R2>0,75) no período de simulação, caracterizando um grau de correspondência significativo para a modelagem da temperatura da água nos pontos de saída selecionados no IPH-ECO. Os dados observados e calculados apresentam correspondência direta com o clima tropical de altitude com temperaturas mais elevadas no período chuvoso (novembro a março) e mais brandas no período seco (abril a outubro).

Praticamente toda energia que controla o metabolismo de um reservatório advém da energia solar. Ferreira e Cunha (2013FERREIRA, D. M.; CUNHA, C. Simulação numérica do comportamento térmico do reservatório do Rio Verde. Revista Engenharia Sanitária e Ambiental, v. 18, n. 1, jan./mar. 2013. ), ao desenvolver e aplicar um modelo unidimensional para simular o comportamento térmico do reservatório do Rio Verde, localizado na Região Metropolitana de Curitiba - PR, observaram que a formação estratigráfica no reservatório depende principalmente da incidência de radiação solar e que os ventos não seriam eficientes em misturar a coluna de água. Porém, para o reservatório de São Pedro, os dados de campo e simulações do comportamento térmico revelam proximidade dos valores de temperatura verificados nos três pontos amostrais caracterizando o reservatório como um ambiente homogêneo (Tabela 2).

Tabela 2
Cálculo da Raiz do Erro Médio Quadrático (REMQ) e coeficiente de determinação R2.

3.2. Representatividade do efeito dos afluentes no regime térmico do reservatório

Com a finalidade de avaliar a influência das vazões afluentes no regime térmico do reservatório de São Pedro, foi realizada uma nova simulação retirando as vazões de entrada e saída. Os resultados desta simulação foram comparados com os dados observados e com a simulação original, considerando as vazões de entrada e saída (Figura 7).

Figura 7
Variação da temperatura do reservatório com e sem os afluentes (Ponto 2).

Nesta circunstância, mostra-se a influência dos afluentes no regime térmico do reservatório, pois a proximidade dos valores observados com os valores calculados mostrou-se relacionada ao fato de se considerar as vazões dos afluentes, ou seja, os afluentes Grota do Pinto e São Pedro controlam a temperatura do reservatório de São Pedro. Fato este também pode ser justificado por conta do tamanho do reservatório e as vazões dos afluentes, onde a mistura e distribuição das temperaturas das águas ocorrem de forma rápida, respeitando a sazonalidade imposta pelo clima tropical de altitude. As alterações climáticas influenciam sobremaneira a dinâmica de temperatura dos afluentes e consequentemente o reservatório.

Destarte, Fantin-Cruz et al. (2010FANTIN-CRUZ, I.; TONDATO, K. K.; MARQUES, D. M.; PEDROLLO, O. Regime térmico em águas correntes e sua importância na estrutura do habitat e na biologia de organismos aquáticos. Revista Caminhos de Geografia, v. 11, n. 36, p. 295-307, 2010.) apontam que o regime térmico dos rios pode ser afetado por outras perturbações provocadas por ações humanas, como efluentes térmicos, redução no fluxo dos rios (irrigação, hidrelétricas) e lançamentos de águas a partir de barragens a montante dos rios. Os lançamentos que atingem os afluentes e o reservatório de São Pedro são basicamente representados por esgotos domésticos e não há barragens a montante. Eles apontam ainda que se devem considerar inclusive alterações na temperatura da água à jusante do reservatório, que são influenciadas por condições ligadas a variáveis operacionais (volume, descarga), limnológicas (tempo de retenção, padrão de estratificação e gradientes térmicos), além da posição da barragem ao longo do perfil longitudinal do rio.

O fluxo dos rios é capaz de influenciar na produtividade biológica de um corpo hídrico, interferindo em processos como eutrofização, ciclagem de nutrientes e atividades fotossintéticas. Para Colischonn e Motta Marques (2003COLISCHONN, W.; MOTTA MARQUES, D. L. Desenvolvimento de modelo bidimensional de transporte de sedimentos e nutrientes para lagos. In: CONGRESSO BRASILEIRO DE LIMNOLOGIA, 9., 2003, Juiz de Fora. Anais... Juiz de Fora: ICB; UFJF, 2003. ), de um modo geral os rios trazem maior parte dos nutrientes e sedimentos para o interior dos lagos e reservatórios, por isso, a circulação originada pelas entradas desses afluentes é particularmente importante.

Os resultados apresentados por Pacheco et al. (2014PACHECO, F. S.; SOARES, M. C. S.; ASSIREU, A. T.; CURTARELLI, M. P.; ROLAND, F.; ABRIL, G. et al. River inflow and retention time affecting spacial heterogeneity of chlorophyll and water-air CO2 fluxes in a tropical hydropower reservoir. Biogeosciences Discuss, v. 11, p. 8531-8568, 2014. http://dx.doi.org/10.5194/bgd-11-8531-2014
http://dx.doi.org/10.5194/bgd-11-8531-20...
), ao investigar o Reservatório do Funil, localizado no município de Resende - RJ, mostram que a variabilidade espacial da concentração de clorofila (Chl) e CO2 está mais relacionada a mudanças no fluxo do rio ao longo do ano do que fatores ambientais, tais como temperatura do ar e radiação solar. Nesse sentido, quando o fluxo do rio atuar como regulador da temperatura de um reservatório, também será capaz de condicionar processos físicos, químicos e biológicos, visto que a temperatura é um fator físico que interfere diretamente na estrutura das comunidades biológicas e processos ecossistêmicos.

A variabilidade das vazões se apresentou conforme esperado para o contexto climatológico, com volumes maiores no período chuvoso e menores no período seco. A maior vazão ocorreu no mês de janeiro de 2011 (0,2294 m3/s) no córrego São Pedro, cuja bacia possui a maior área de contribuição.

4. CONCLUSÃO

Foi possível representar o regime térmico do reservatório de São Pedro por meio da modelagem hidrodinâmica 2DH, com base em dados meteorológicos e hidrológicos, utilizando o modelo IPH-ECO. O contexto climatológico como principal fator influenciador das condições térmicas da água a partir das trocas de calor da atmosfera com a superfície de água foi fundamental para realização deste trabalho, visto que os resultados das simulações mostraram valores de temperatura em consonância direta com as condições do clima tropical de altitude.

Nossos resultados mostraram temperaturas simuladas, em média 22ºC, próximas às temperaturas observadas, em média 22,21°C. A simulação satisfatória do regime térmico é essencial para investigar os processos ecossistêmicos, uma vez que a temperatura é um fator que influencia diretamente nesses processos. Essa representatividade significativa no modelo suportaria testar cenários com mudanças de temperatura e garantir ações previamente necessárias para prevenir eventuais efeitos negativos à diversidade biológica, além da operação e manutenção do sistema de abastecimento.

A contribuição dos afluentes se mostrou muito importante para as simulações hidrodinâmicas, pois a variação das vazões limita ou regulariza as temperaturas do reservatório de São Pedro. Por se tratar de um reservatório pequeno e raso, as temperaturas são mais limitadas às vazões dos afluentes caracterizando, assim, uma baixa capacidade deste reservatório de regularização de suas temperaturas.

5. AGRADECIMENTOS

Os autores agradecem ao Programa de Pós-Graduação em Ecologia da Universidade Federal de Juiz de Fora (PGECOL), ao Departamento de Engenharia Sanitária e Ambiental da Universidade Federal de Juiz de Fora, à Embrapa Gado de Leite (CNPq), à Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES) e ao Corpo de Bombeiros de Juiz de Fora pelo apoio ao trabalho. Os autores agradecem ao CNPq, Processo número 473854/2008-2 e à Embrapa código Infoseg 06.08.06.002.00.00 pelo apoio financeiro ao trabalho.

  • CASULLI, V.; CHENG, R. Semi-implicit finite difference methods for three-dimensional shallow water flow. International Journal for Numerical Methods in Fluids, v. 15, p. 629-648, 1992. http://dx.doi.org/10.1002/fld.1650150602
    » http://dx.doi.org/10.1002/fld.1650150602
  • COMPANHIA DE SANEAMENTO MUNICIPAL DE JUIZ DE FORA - CESAMA. Website. Disponível em: http://www.cesama.br. Acesso: em abr. 2013.
    » http://www.cesama.br.
  • CHRISTOFOLLETI, A. Geomorfologia. São Paulo: Blucher, 1980.
  • COLISCHONN, W.; MOTTA MARQUES, D. L. Desenvolvimento de modelo bidimensional de transporte de sedimentos e nutrientes para lagos. In: CONGRESSO BRASILEIRO DE LIMNOLOGIA, 9., 2003, Juiz de Fora. Anais... Juiz de Fora: ICB; UFJF, 2003.
  • EUCLYDES, H. P. Saneamento agrícola: atenuação de cheias: metodologia e projeto. Belo Horizonte: Rural Minas, 1987.
  • FANTIN-CRUZ, I.; TONDATO, K. K.; MARQUES, D. M.; PEDROLLO, O. Regime térmico em águas correntes e sua importância na estrutura do habitat e na biologia de organismos aquáticos. Revista Caminhos de Geografia, v. 11, n. 36, p. 295-307, 2010.
  • FERREIRA, D. M.; CUNHA, C. Simulação numérica do comportamento térmico do reservatório do Rio Verde. Revista Engenharia Sanitária e Ambiental, v. 18, n. 1, jan./mar. 2013.
  • FRAGOSO JR., C. R.; FERREIRA, T. F.; MOTTA MARQUES, D. M. L. Modelagem ecológica em ecossistemas aquáticos. 1. ed. São Paulo: Oficina de Textos, 2009. v. 1. 304p.
  • FRAGOSO JR., C. R. Modelagem tridimensional da estrutura trófica em ecossistemas aquáticos continentais. 2009. Tese (Doutorado em Recursos Hídricos e Saneamento Ambiental) - Instituto de Pesquisas Hidráulicas, Universidade Federal do Rio Grande do Sul, Porto Alegre, 2009.
  • INSTITUTO ALBERTO LUIZ COIMBRA DE PÓS-GRADUAÇÃO E PESQUISA DE ENGENHARIA - COPPE/UFRJ. Manual do Usuário do SisBaHia. 2007. Disponível em: http://www.sisbahia.coppe.ufrj.br/Arquivos_Gerais Acesso em: jul. 2013.
    » http://www.sisbahia.coppe.ufrj.br/Arquivos_Gerais
  • JUIZ DE FORA. Prefeitura Municipal. Plano diretor de desenvolvimento urbano de Juiz de Fora. Juiz de Fora: FUNALFA, 1999.
  • MACHADO, P. J. O. (Org.) Diagnóstico físico-ambiental da bacia hidrográfica do Córrego São Pedro: um exercício acadêmico de gestão dos recursos hídricos. Ubá: Geographica, 2010.
  • MOOIJ, W. M.; TROLLE, D.; JEPPESEN, E.; ARHONDITSIS, G.; BELOLIPETSKY, P. V.; CHITAMWEBWA, D. B. R. et al. Challenges and opportunities for integrating lake ecosystem modeling approaches. Aquatic Ecology, v. 44, 2010. http://dx.doi.org/10.1007/s10452-010-9339-3
    » http://dx.doi.org/10.1007/s10452-010-9339-3
  • PACHECO, F. S.; SOARES, M. C. S.; ASSIREU, A. T.; CURTARELLI, M. P.; ROLAND, F.; ABRIL, G. et al. River inflow and retention time affecting spacial heterogeneity of chlorophyll and water-air CO2 fluxes in a tropical hydropower reservoir. Biogeosciences Discuss, v. 11, p. 8531-8568, 2014. http://dx.doi.org/10.5194/bgd-11-8531-2014
    » http://dx.doi.org/10.5194/bgd-11-8531-2014
  • ROSMAN, P. C. C.; MASCARENHAS, F. C. B.; MIGUEZ, M. G.; CAMPOS, R. O.; EIGER, S. Métodos numéricos em recursos hídricos. 5. ed. Rio de Janeiro: ABRH, 2001.
  • ROSMAM, P. C. C. Subsídios para modelagem de sistemas estuarinos. In: ROSMAN, P. C. C.; MASCARENHAS, F. C. B.; MIGUEZ, M. G.; CAMPOS, R. O.; EIGER, S. Métodos numéricos em recursos hídricos, Rio de Janeiro: ABRH 1999.
  • TUCCI, C. E. M. Hidrologia: ciência e aplicação. Porto Alegre: Ed. UFRGS, 2001.
  • WANG, S.; QIAN, X.; HAN, B.; LUO, L.; HAMILTON, D. P. Effects of local climate and hydrological conditions on the thermal regime of a reservoir at Tropic of Cancer, in southern China. Water Research, v. 46, n. 8, p. 2591-2604, 2012. http://dx.doi.org/10.1016/j.watres.2012.02.014
    » http://dx.doi.org/10.1016/j.watres.2012.02.014

Datas de Publicação

  • Publicação nesta coleção
    Mar 2016

Histórico

  • Recebido
    20 Ago 2015
  • Aceito
    26 Set 2015
Instituto de Pesquisas Ambientais em Bacias Hidrográficas Instituto de Pesquisas Ambientais em Bacias Hidrográficas (IPABHi), Estrada Mun. Dr. José Luis Cembranelli, 5000, Taubaté, SP, Brasil, CEP 12081-010 - Taubaté - SP - Brazil
E-mail: ambi.agua@gmail.com