Hidrologia

Post on 04-Jan-2016

28 views 3 download

description

Hidrologia. Evapotranspiração. Carlos Ruberto Fragoso Jr. http://www.ctec.ufal.br/professor/crfj/ Marllus Gustavo Ferreira Passos das Neves http://www.ctec.ufal.br/professor/mgn/ Ctec - Ufal. Evapotranspiração. Conceito Geral Fatores que afetam a evapotranspiração Medição da evaporação - PowerPoint PPT Presentation

Transcript of Hidrologia

Hidrologia

Evapotranspiração

Carlos Ruberto Fragoso Jr.http://www.ctec.ufal.br/professor/crfj/

Marllus Gustavo Ferreira Passos das Neveshttp://www.ctec.ufal.br/professor/mgn/

Ctec - Ufal

• Conceito Geral

• Fatores que afetam a evapotranspiração

• Medição da evaporação

• Evaporação em lagos e reservatórios

• Estimativa da evapotranspiração– Medição– Cálculo

Evapotranspiração

Evaporação (E) – Taxa de conversão da água líquida em vapor, água esta presente nos oceanos, lagos, rios e solo

Transpiração (T) – parte do total evaporado para a atmosfera proveniente do solo, através das plantas

Conceito Geral - Evapotranspiração

DefiniçõesVapor d’água moléculas 10 vezes mais distantes umas das outras que na fase líquida Calor latente de vaporização da água Energia requerida para separá-las

T0,0023612,501 λ

em MJ.kg-1

T temperatura na superfície da água em ºC

Definições

E diferença entre duas taxas: taxa de vaporização (função da temp.) e taxa de condensação (determinada pela pressão de vapor)Diferença nula ar saturado (equilíbrio)

Evaporação natural trocas de moléculas de água entre o ar e a superfície livre

T237,3T17,27

exp0,6108eS

T é a temperatura em ºC

em kPa

DefiniçõesTão importante quando es é o seu gradiente = des/dT

2s

T237,3

e4098

Δ em kPa.C-

1

Calor sensível parte da energia radiante que atinge a superfície da terra que não é utilizada para E aquece a atmosfera em contato com o solo movimento ascendenteH fluxo de calor sensível (fluxo de calor por unidade de área)

Evapotranspiração (ET) Processo simultâneo de transferência de água para a atmosfera através da evaporação (E) e da transpiração (T).

TEET

Potencial (ETP)

Real (ETR)

Transpiração no Sistema Solo Planta Atmosfera

Transpiração das raízes às folhas, pelo sistema condutor, pelo estabelecimento de umgradiente de potencial desde o solo até o ar

Transpiração no sistema solo planta atmosfera

Local de maior resistência ao fluxo

O gradiente de tensão de vapor de água também favorece o fluxoQuanto mais seco estiver o ar (menor umidade relativa), maior será esse gradiente

proporcional à resistência ao fluxo da água na planta

DefiniçõesETP Quant. de água transferida para a atmosfera por E e T, em uma unidade de tempo, de uma superfície extensa, completamente coberta de vegetação de porte baixo e bem suprida de água (Penman,1956)

ETR Quant. de água transferida para a atmosfera por E e T, nas condições reais (existentes) de fatores atmosféricos e umidade do solo. A ETR é igual ou menor que a evapotranspiração potencial (Gangopadhyaya et al, 1968)

Fatores que afetam

• Radiação solar•Umidade do ar• Temperatura do ar• Velocidade do vento•Tipo de solo• Vegetação (transpiração)

19%

26%

4%

51% + 19% = 70%

evaporação

ar quente

Radiação Solar

30% 70%

Radiação Solar

Radiação no topo da atmosfera (Stop)

St (Ssup) Radiação incidente

albedo

Sto

p

Ssu

p

RL

.Ssu

p

)(1SR SUPL αTOPssSUP SNn

baS

Radiação Solar

MJ.m-2.dia-1 TOPssSUP SNn

baS

N insolação máxima possível em horasn isolação medida em horas a fração de atinge a superfície em dias encobertos (quando n=0)b fração de atinge a superfície em dias sem nuvens (quando n = N)

Radiação Solar

• Quando não existem dados locais medidos que permitam estimativas mais precisas, são recomendados os valores de 0,25 e 0,50, respectivamente, para os parâmetros as e bs;

• Quando a estação meteorológica dispõe de dados de insolação, a equação acima é utilizada com n medido e N estimado pela equação. Quando a estação dispõe de dados de fração de cobertura, utiliza-se o valor de n/N diretamente

Radiação Solar

Radiação Solar

para

a

su

perfí

cie

Ln radiação líquida de ondas longas que deixa a superfície terrestre 4

n 273,2)(TfL σε

Radiação Solar

f fator de correção devido à cobertura de nuvensT [ºC] temperatura média do ar a 2 m do solo emissividade da superfícies constante (σ = 4,903.10-9 MJ.m-2.ºK-4.dia-1)

4n 273,2)(TfL σε MJ.m-2.dia-1

Balanço por unidade de áreaRL Radiação líquida (ondas

curtas)

Calor transferido para o solo por condução

Fluxo de calor sensível Energia usada

na evaporação

Arm. temp. Energia

absorvida por processos bioquímicoso

didd AAA

Perda associada com o mov. do ar

Balanço por unidade de áreaRL

E + H = RL – G – S – P - Ad

A = E+HA = RL – G – S – P - Ad

A energia disponível para calor latente e para o calor sensível

MJ.m-2.dia-1

Geralmente 2% de RL

P

Geralmente negligenciados

S Ad

Por simplicidade, G pode ser considerado nulo

3dd TT0,38G

• Td é a temperatura do solo no dia que se deseja calcular a ET• T3d é a temperatura do solo 3 dias antes

Balanço por unidade de área

Condução principal mecanismo de transferência de

calor para o solo

Para flutuações diárias de temperatuta (prof. efetiva

de um solo típico igual a 0,18 m)

• Quanto maior a temperatura, maior a pressão de saturação do vapor de água no ar, isto é, maior a capacidade do ar de receber vapor

Temperatura (oC)

0 10 20 30

P0 (atm) 0,0062 0,0125

0,0238 0,0431

Temperatura

• Para cada 10oC, P0 é duplicada

Umidade relativa (UR) medida do conteúdo de vapor de água do ar em relação ao conteúdo de vapor que o ar teria se estivesse saturado

sww

100UR(%)

Ondew massa de vapor por massa de ar e ws é a massa de vapor por massa de ar no ponto de saturaçãoe pressão parcial de vapor no ar e es é pressão de saturação

Umidade do Ar

see

100UR(%)

• Renova o ar em contato com a superfície que está evaporando (superfície da água; superfície do solo; superfície da folha da planta)

• Vento forte turbulência maior transferência para regiões mais altas da atmosfera mais rápida umidade próxima à superfície menor aumenta taxa de E

Vento

Vento remove ar úmido da superfície onde ocorre ET menos umidade mais ET

Vegetação e solo•A vegetação:

o Controla a transpiraçãoo Pode agir fechando os estômatoso Busca a umidade de camadas profundas do solo

•Umidade do solo uma das variáveis mais importantes na transpiração

o Solo úmido plantas transpiram livremente taxa de transpiração controlada pelas variáveis atmosféricas

o Solo começa a secar fluxo de transpiração começa a diminuir

o Condições ideais de umidade do solo ETPo Condições reais de umidade do solo ETR

•Solos arenosos úmidos tem evaporação maior do que solos argilosos úmidos

Determinação da evaporação e da ET

Relação entre a evaporação e a pressão de vapor, com a introdução do efeito do vento

Leva em conta a radiação solar: efetiva de ondas curtas, efetiva de ondas longas, a energia de evaporação, calor sensível por condução, características aerodinâmicas método de Penman

Ajuste por regressão das variáveis envolvidas

Medida direta tanque classe A, ...

Baseia-se na equação da continuidade do lago ou reservatório

Evaporação

Lisímetros e umidade do solo

ETP Método de thornthwaite, método de Blaney-Criddle. Para determinar ET ET = ETP .kc, onde kc coeficiente de cultura (determinado em lisímetros)

Baseados na variável meteorológica radiação. Equação de Jesen e Haise, ...

Chamada de equação de Penman adaptar o cálculo da evaporação de superfícies livres para a superfície de interesse ETP

Para intervalos de tempo superiores a 1 semana

Evapotranspiração

• Tanque classe A• Evaporímetro de Piché

Evaporímetros medição direta

Tanque Classe A

• O mais usado forma circular com um diâmetro de 121 cm e profundidade de 25,5 cm• Construído em aço ou ferro galvanizado• Pintado na cor alumínio• Instalado numa plataforma de madeira a 15 cm da superfície do solo• Permanecer com água variando entre 5,0 e 7,5 cm da borda superior

Tanque classe A

Tanque "Classe A" – US Weather Bureau

• O fator que relaciona a evaporação de um reservatório e do tanque classe A oscila entre 0,6 e 0,8, sendo 0,7 o valor mais utilizado

Tanque classe A

Fonte : Sabesp

Tanque classe A

Tanque classe A

Tanque classe A

• manutenção da água entre as profundidades recomendadas evita erros de até 15%• a água deve ser renovada turbidez evita erros de até 5%• as paredes sofrem com a influência da radiação e da transferência de calor sensível superestimação da evaporação• próximos a cultivos de elevada estatura subestimação da evaporação

Constituído por um tubo cilíndrico, de vidro (≈ 30 cm) de comprimento e um 1 cm de diâmetro, fechado na parte superior e aberto na inferior

A extremidade inferior tapada, depois do tubo estar cheio com água destilada, com um disco de papel de feltro, de 3 cm de diâmetro, que deve ser previamente molhado com água

Este disco é fixo depois com uma mola. A seguir, o tubo é preso por intermédio de uma argola a um gancho situado no interior do abrigo

Evaporímetro de Piché

Evaporímetro de Piché

• Piché é pouco confiável

• Medição (mais complicada)• Cálculo

Estimativa da evapotranspiração

• Lisímetro– Depósitos enterrados, abertos na

parte superior, preenchidos com solo e vegetação característica

– Controle das variáveis:• Peso• Medir chuva• Coletar água percolada• Coletar água escoada• Superfície homogênea

Lisímetros medição direta

Lisímetros medição direta

Precipitação no solo drenagem para o fundo do aparelho água é coletada e medida

O depósito é pesado diariamente, assim como a chuva e os volumes escoados de forma superficial e que saem por orifícios no fundo

ET calculada por balanço hídrico entre 2 dias subseqüentes

ET = P - Qs – Qb – ΔV

E evapotranspiraçãoP chuva (medida num pluviômetro)Qs escoamento superficial (medido)Qb é o escoamento subterrâneo (medido no fundo do tanque)ΔV variação de volume de água (medida pelo peso)

Lisímetros medição direta

Lisímetros medição direta

Lisímetros medição direta

http://jararaca.ufsm.br/websites/matasul-ufsm/1ca53f95af2a6c15feea202899377cc9.htm

• Usando a temperatura e a umidade do ar

• Usando a temperatura e a radiação solar

• Equações de Penmann (insolação, temperatura, umidade relativa, velocidade do vento)

Mais Equações de cálculo da ET

• Jensen Haise, Turc, Grassi, Stephens – Stewart, Makkink, Blaney-Morin, Hamon, Hargreaves, Papadakis

• Estimação a partir da temperatura somente é recomendada se não houver outra opção equação de Hargreaves provoca erros da ordem de 10 a 15% Blaney-Criddle diminui estes erros

Métodos baseados na temperatura e radiação

Cálculo da ETP baseado na temperatura

Thornthwaite: empírica, caracterizada por um único fator, a temperatura média. Foi desenvolvida para climas temperados (inverno úmido e verão seco)

E = c.Ta

t = temperatura de cada mês ºCT = temperatura média ºC

Blaney-Criddle: também utiliza a temperatura média e horas do dia com insolação, para regiões semi-áridas

ETP=(0,457.T + 8,13).pET = ETP.kc

p = % luz diáriakc = é o coeficiente de cultura.

Cálculo da ETP baseado na temperatura

a

IT10

16ET

12

1j

1,514

j

5

TI

Para estimar evapotranspiração potencial mensalT = temperatura média do mês (oC)a = parâmetro que depende da regiãoI = índice de temperatura

0,49239I101,792I107,71I106,75a 22537

Thornthwaite

j cada um dos 12 meses do anoTj temperatura média de cada um dos 12 meses

ExemploMês Temperatura

Janeiro 24,6

Fevereiro 24,8

Março 23,0

Abril 20,0

Maio 16,8

Junho 14,4

Julho 14,6

Agosto 15,3

Setembro 16,5

Outubro 17,5

Novembro 21,4

Dezembro 25,5

Calcule a evapotranspiração potencial mensal para o mês de Agosto de 2006 em Porto Alegre onde as temperaturas médias mensais são dadas na figura abaixo.

Suponha que a temperatura média de agosto de 2006 tenha sido de 15,3°C

Exemplo

O primeiro é o cálculo do coeficiente I a partir das temperaturas médias obtidas da tabela. O valor de I é 96. A partir de I é possível obter a= 2,1. Com estes coeficientes, a evapotranspiração potencial é:

mm/mês 53,19616,510

16E2,1

Portanto, a evapotranspiração potencial estimada para o mês de agosto de 2006 é de 53,1 mm/mês.

Equações combinadas

•Penman evaporação•Christiansen•Van Bavel•Penman - Monteith ampliação de Penman para ETR de uma superfície vegetada

• Combina – poder evaporante do ar (temperatura,

umidade, velocidade do vento)– poder evaporante da radiação

W

a

s

a

dspAL

rr

ree

cGR

E

1

1

Penman

Penman - MonteithPenman + introdução de um fatores de resistência que leva em consideração o stress de umidade da vegetação e do solo

Analogia com a resistência elétrica

Penman - MonteithToda a energia disponível para E é acessível pela planta e o vapor d´água se difunde contra

Resistência superficial (estômatos)

Resistência aerodinâmica

Penman - Monteith

Resistência estomática das folhas resistência ao fluxo de vapor por difusão molecular (solo-caule-galho-folhas (estômatos)

Resistência superficial da vegetação (rs)

Estômatos abrem ou fecham segundo a demanda por UR (atmosfera) e umidade do solo)

Diferente para os diversos tipos de plantas e depende de variáveis ambientais(umidade do solo, temperatura do ar e radiaçãorecebida pela planta)

Penman - Monteith

Valores de referência (boas condições de umidade)– Grama: rs = 69 s/m (ETP)– Florestas superficiais: rs = 100 s/m

2ma

s U0,34r

r

Resistência estomática das folhas depende da disponibilidade de água no solo em condições favoráveis valores de resistência estomática e, em consequência, os de rs mínimos

Resistência superficial da vegetação (rs)

Penman - Monteith

Controla a taxa de transferência de vapor d’água do solo (difusão turbulenta) inversamente proporcional ao vento e varia com a altura da vegetação que cobre o solo

Resistência aerodinâmica (ra)

Penman - Monteith

m10hparau

94r

m10hparaz

10ln

u

6,25r

m,10a

2

0m,10a

Velocidade do vento a 10 m de altura (m/s)

z0 é a rugosidade da superfície

h é a altura da vegetação em m

Resistência aerodinâmica (ra)

Penman - Monteith

Inversamente proporcional à altura dos obstáculos enfrentados pelo vento, pois estes que geram a turbulência

Resistência aerodinâmica (ra)

ra

Penman - Monteith

Velocidade do vento a 10 m de altura

0

0m,2m,10

z2

ln

z10

ln

uu

• um,2 é a velocidade do vento a 2 m de altura em m/s• z0 é a rugosidade da superfície

Estações climatológicas normalmente dispõe de dados de velocidade do vento medidas a 2 m de altura. Para converter estes dado a uma altura de referência de 10m é utilizada a equação ao lado

W

a

s

a

dspAL

rr

ree

cGR

E

1

1

Penman - Monteith

água; da específica massa ][kg.m

ar; do específica massa ][kg.m

solo; o para energia de fluxo ]s.[MJ.m

;superfície na líquida radiação ]s.[MJ.m

vapor;do saturação de pressão da variaçãode taxa]C[kPa.

ão; vaporizaçde latentecalor ][MJ.kg

água; da evaporação de taxa][m.s

3-

3-

1-2-

1-2-

1-

1-

-1

W

A

L

G

R

E

)(1SR SUPL α

3dd TT0,38G

T0,0023612,501λ

RL

A = E+HA = RL – G

A energia disponível para calor latente e para o calor sensível

MJ.m-2.dia-1

Penman - Monteith

W

a

s

a

dspA

r

r

r

eecA

E

1

1

• Vegetação molhada (durante ou imediatamente após uma chuva ou depois da irrigação) fonte do vapor de água não está dentro das folhas os estômatos não são mais efetivos em restringir a evaporação rs = 0

• Evaporação potencial superfície extensa de água livre rs = 0

Penman - Monteith

W

a

dspA r

eecA

E

1

ETR

ETP

Umidade do solo Smx

ETR = evapotranspiração depende da umidade do solo

Relações

Relações

Períodos de estiagem mais longos ET retira umidade do solo ET diminui

A redução da ET não ocorre imediatamente

Para valores de umidade do solo entre a capacidadede campo e um limite ET não é afetada pelaumidade do solo

A partir deste limite ET diminuída mínimo(normalmente zero) no ponto de murcha permanente

Neste ponto rS atinge valores altíssimos

Penman - Monteith

ca;aerodinâmo aresistênci ]s.m[ r

vegetação; da al superficiaresistênci ]s.m[ r

0,66);( icapsicrométr constante ]C[kPa.

vapor; do pressão ][kPa e

vapor; do saturaçãode pressão ][kPa e

);C. MJ.kg10.013,1(C úmido ar do específico calor ]C.[MJ.kg C

1-a

1-s

1-

s

s

113p

-1-1p

m2a

s U34,0r

r

Penman - Monteith

•Massa específica do ar T275P

3,486 A

PA é a pressão atmosférica em kPaT é a temperatura do ar a 2m da superfície em ºC

• Massa específica do água

594634

23W

2W

WWW

W

T106,536332T101,120083T101,001685

T109,09529T106,793952999,845259ρ

TW é a temperatura da água em ºC

Penman - Monteith

• Pressão de saturação do vapor (es)

T237,3

T17,27exp0,6108eS

T é a temperatura do ar a 2m da superfície em ºC• Pressão real de vapor de água no ar (ed)

100

URee Sd

UR é a umidade relativa do ar em %

Penman - Monteith

• Calor latente de vaporização (λ) T0,0023612,501λ em MJ.kg-1

T é a temperatura do ar a 2m da superfície em ºC

• Constante psicrométrica (γ) )(constante0,66γ

• Taxa de variação da pressão de saturação do vapor com a temperatura do ar (Δ)

2s

T237,3

e4098Δ

T é a temperatura do ar a 2m da superfície em ºC

Penman – Monteith passos

1. Obter o dia Juliano (J) para a data que se deseja calcular a ET

2. Obter a latitude (), em graus, do local que se deseja calcular a ET

3. Calcular a declinação solar em radianos

4. Calcular a distância relativa da terra ao sol (dr)

1,405J3652π

sen0,4093

J

3652

cos0,0331dr

5. Calcular o ângulo ao nascer do sol em radianos (ωs) δω tantanarccos s

é a latitude do local em radianos é declinação solar em radianos

Penman – Monteith passos

6. Calcular a insolação máxima (N) para a localização desejada

s

24N ω

π

7. Calcular a radiação solar que atinge o topo da atmosfera (STOP), em MJ.m-2.dia-1

)( ss sencoscossensend

100015,392S r

WTOP

é a latitude do local em radianos é declinação solar em radianosωs é o ângulo do sol ao nascer em radianos

Penman – Monteith passos

Penman – Monteith passos• Radiação no topo da atmosfera (Stop) função da latitude, distância sol-terra e época do ano

8. Calcular a radiação solar que atinge o topo da superfície (SSUP), em MJ.m-2.dia-1

TOPssSUP SNn

baS

N insolação máxima possível em horasn isolação medida em horas a fração de atinge a superfície em dias encobertos (quando n=0)b fração de atinge a superfície em dias sem nuvens (quando n = N)

Penman – Monteith passos

• Estime a evapotranspiração média, em mm/dia, através da equação de Penman-Monteith para a cidade de Maceió (posto Inmet Ufal), no sábado, dia 23/04/2011

Exercício

Abrir planilha Acessar http://www.inmet.gov.br/ observações estações automáticas ou convencionais localizar a estação da Ufal no mapa colocar a latitude na planilha clicar em dados escolher data e baixar os dados podem ser colocados no Excel

Exercício

Dados meteorológicos

http://meteo.infospace.ru/main.htm

Dados meteorológicos

• Método de estimativa simples com base nos dados precipitação e vazão de uma bacia

• A equação da continuidade

S(t+1)=S(t) + (P –E - Q).t

• Desprezando a diferença entre S(t+1) – S(t)

Q= P- E

• Simplificação aceita para t longos como o um ano

ou sequência de anos

Balanço hídrico

• Exemplo:

Uma bacia (Rio Passo Fundo) com Precipitação média

1.941 mm e Vazão de 803 mm (valores médios de 10

anos). A evaporação real é E= 1941 – 803 = 1137 mm

O coeficiente de escoamento é a relação entre Q/P

C = 803/1941 = 0,41

ou 41% da precipitação gera escoamento.

Balanço hídrico

mm/ano m3/s

A = Área da baciaQ = vazão

1000365243600

)A(kmQ(mm/ano)/s)Q(m

23

Conversão de unidades

Reservatórios são criados para regularizar a vazão dos rios, aumentando a disponibilidade de água e de energia nos períodos de escassez

Entretanto, cria-se uma vasta superfície líquida que disponibiliza água para evaporação, o que pode ser considerado uma perda de água e de energia

Evaporação em reservatórios

A evaporação da água em reservatórios estimada a partir de medições de tanques de Classe A

Entretanto é necessário aplicar um coeficiente de redução em relação às medições de tanque a água do reservatório normalmente está mais fria do que a água do tanque, que tem um volume pequeno e está completamente exposta à radiação solar

Evaporação em reservatórios

Assim, para estimar a evaporação em reservatórios e lagos costuma-se considerar que esta tem um valor de aproximadamente 60 a 80% da evaporação medida em Tanque Classe A na mesma região, isto é:

Onde Ft tem valores entre 0,6 e 0,8.

ttanquelago FEE

Evaporação em reservatórios

Evaporação em lagos e reservatórios

Reservatório de Sobradinho

constituindo-se no maior lago artificial do mundo, está numa das regiões mais secas do Brasil

área superficial de 4.214 km2

Evaporação direta deste reservatório é estimada em 200 m3.s-1 10% da vazão regularizada do rio São Francisco Esta perda é superior à vazão prevista para o projeto de transpiração do rio São Francisco

• Um rio cuja vazão média é de 34 m3/s foi represado por uma barragem para geração de energia elétrica. A área superficial do lago criado é de 5.000 hectares.

Medições de evaporação de um tanque classe A correspondem a 1.500 mm por ano.

Qual é a nova vazão média a jusante da barragem após a formação do lago?

Exercício

1000365243600

)A(kmE(mm/ano)/s)E(m

23

E = 1.500 x 0,7 mm/ano

E = 1,66 m3/s

Q = 34 – 1,66 = 32,34 m3/s

Redução de 4,9 % da vazão

Solução

• Deseja-se construir um reservatório em um rio, cuja bacia possui uma área de 50 km2. A área de inundação do reservatório é de 10 km2. Estime qual deve ser a redução de vazão média disponível na bacia. Considere que a evaporação potencial da superfície da água é de 1.400 mm por ano. A evaporação estimada por balanço hídrico antes da construção do reservatório foi de 1.137 mm por ano. Nestas mesmas condições, a vazão média era de 1,41 m3/s e a precipitação de 1.941 mm por ano.

Exercício

• ET após a construção

ET = (0,7.1.400.10 + 1.137.40)/50 = 1.105,6 mm/ano

• Q após a construção

Q = 1.941 - 1.105,6 = 835,4 mm/ano

• Redução de Q

Qantes = 1,41 m3/s

Qdepois = 835,4 mm/ano = 1,325 m3/s

Redução de 6,45%

Exercício