NOTA TÃCNICA - UFV

More documents

Recommendations

Info

£ ¢ ¡ ¥ ¥ 9 7 8 A função densidade de probabilidade gama foi selecionada para representação da duração da precipitação pois, segundo Sediyama et al. (1978), ajusta se muito bem a variáveis aleatórias contínuas em climatologia, sendo sua função acumulativa de probabilidade expressa por x 1 x - F ( x) ¤ -1 = x e dx ( ) 0 (2) em que, F(x) é a probabilidade de ocorrência do evento x; x é a variável aleatória contínua; α é o parâmetro de posição de x; β é o parâmetro de escala de x; e Γ (α ) é a função gama do parâmetro α. Os parâmetros (¦ §©¨¨¨ ! " #$ de posição ção gama foram obtidos para cada mês e para cada faixa de total precipitado. Nos casos em que haviam poucos eventos por faixa, como é o caso de junho a agosto na região Sudeste, optou-se por agrupar os eventos de todas as estações localizadas em uma mesma região homogênea e, a partir daí, determinar os * parâmetros %'&)( Para a geração dos valores de duração para os dias chuvosos, foi utilizada a metodologia proposta por Sediyama et al. (1978), a qual consiste em, primeiramente, gerar um conjunto de números aleatórios, variando entre zero e um, com extensão igual ao número de dias chuvosos da série sintética gerada. Substitui-se, então, a probabilidade de ocorrência da precipitação (Equação 2) pelo valor do número aleatório obtido para o dia, resolvendo-se a equação, sendo x o valor obtido para a duração da precipitação para o evento considerado. A resolução desta equação foi efetuada pelo método das somas de Riemann (Zanetti, 2003). Intensidade máxima instantânea de precipitação A intensidade máxima instantânea possui alta correlação com o total precipitado e baixa correlação com a duração da precipitação. Com isso, a intensidade máxima instantânea foi gerada de forma semelhante à duração da precipitação, com base em uma distribuição de probabilidade gama, sendo os +-,/. 0©12!3 405 parâmetros para cada mês e para cada faixa de total precipitado, considerando-se todos os eventos das estações de uma mesma região homogênea. Tempo padronizado de ocorrência da intensidade máxima instantânea De forma semelhante à forma como foram gerados os dados de duração efetiva e da intensidade máxima instantânea, procedeu-se à geração de dados de tempo padronizado de ocorrência da intensidade máxima instantânea, obtido pela razão entre o tempo de ocorrência da intensidade máxima instantânea e a duração efetiva do evento, considerando-se a correlação com o total precipitado e dividindo-se os eventos por faixa de total precipitado, para a obtenção dos parâmetros de entrada do modelo. Entretanto, para a geração desta variável, utilizou-se distribuição de probabilidade exponencial, que foi a distribuição com melhor ajuste, tendo em vista os resultados obtidos pelos testes de Kolmogorov-Smirnov (20% de significância) e Qui-quadrado (5% de significância). A função acumulativa de probabilidade da distribuição exponencial é expressa por 1 7 6 ; : F(x) = þ f(x) = þ dt = 1- e (3) 0 0 - em que, < é o único parâmetro da distribuição, sendo expresso pelo inverso da média. Temperaturas máxima e mínima e radiação solar A geração de séries sintéticas de temperaturas máxima e mínima e radiação solar foi realizada, com base em processo estocástico multivariado contínuo, conforme metodologia proposta por Matalas (1967) e aplicada para estas três variáveis por Richardson (1981). As séries destas variáveis são geradas, de forma conjunta, - Engenharia na Agricultura, Viçosa, MG, v.13, n.3, 210-220, Jul./Set., 2005 213
= > @ @ @ considerando-se o alto grau de associação existente entre elas, com base nas equações T + máx = Tmáx máx T' máx (4) mín = T mín mín T' mín (5) T + RAD = RAD + σrad RAD' (6) em que, T e T são os valores de máx mín temperaturas máxima e mínima geradas para cada dia da série sintética, ºC; RAD é o valor de radiação solar para cada dia da série sintética; J m -2 ; T máx e T mín são as médias das séries observadas de temperaturas máxima e mínima, ºC; RAD é a média da série observada de radiação solar, J m -2 ? ; e são os desvios máx ? padrão das séries observadas de temperatura ? máxima, mínima, ºC; rad é o desvio padrão da série observada de radiação solar, J m -2 ' ' ; e T máx , T mín e RAD são os resíduos gerados para cada variável, adimensionais. Os resíduos foram gerados para cada dia da série sintética, utilizando-se a equação ⎡ ' Tmáx ⎢ ' ⎢ Tmín ⎢RAD ⎣ ' ⎤ ⎥ ⎥ = A ⎥ ⎦ ⎡ ' Tmáx ⎢ ⎢ ' Tmín ⎢ ⎢RAD ⎣ -1 -1 ' -1 mín ⎤ ⎡ ⎥ ⎢ ⎥ + B⎢ ⎥ ⎢ ⎥ ⎦ ⎣ T máx T mín RAD ⎤ ⎥ ⎥ ⎥ ⎦ (7) em que, T' máx , T' mín e RAD' - 1 -1 -1 são os resíduos para o dia anterior; ∈ ∈ ∈ Tmáx Tmín RAD são os números aleatórios independentes e normalmente distribuídos, com média igual a zero e desvio padrão igual a um, gerados para cada dia da série sintética, adimensionais; e A e B são as matrizes de ordem três, que consideram o grau de correlação entre as variáveis, sendo expressas pelas equações A −1 = ⋅ (8) B⋅ B M 1 M 0 T = M 0 − M 1 ⋅ M −1 0 ⋅ M T 1 (9) em que, M 0 é a matriz de correlação entre as variáveis sem defasagem e M 1 é a matriz de correlação, considerando-se defasagem de um dia entre as variáveis consideradas. Os resíduos de cada dia da série observada foram calculados, subtraindo-se a média e dividindo o resto pelo desvio-padrão da série de dados observados, sendo que, para o primeiro dia da série, a média dos valores dos resíduos de cada variável foi considerada no lugar do resíduo do dia anterior. O emprego da Equação 7 implica que os resíduos das variáveis são, normalmente, distribuídos e que a correlação serial de cada variável pode ser descrita por um modelo de auto-regressão linear de primeira ordem (Matalas, 1967). A matriz B pode ser qualquer matriz que seja a solução desta equação, podendo ser obtida pelo método da decomposição triangular inferior de Cholensky (Bras & Rodriguez-Iturbe, 1985; Press et al., 1986). Velocidade do vento A velocidade média diária do vento foi estimada a partir de uma distribuição de probabilidade Pearson tipo III normalizada, cuja função acumulativa de probabilidade é expressa pela Equação 1. Os valores diários são gerados da mesma forma que os valores de total precipitado diário, usando-se a média, o desvio-padrão e o coeficiente de assimetria dos valores observados de velocidade do vento. Umidade relativa do ar A distribuição de probabilidade beta foi utilizada para a geração da umidade relativa do ar no modelo para geração de dados climáticos. A sua função de distribuição acumulada de probabilidade pode ser expressa pela equação 214 Engenharia na Agricultura, Viçosa, MG, v.13, n.3, 210-220, Jul./Set., 2005
Page 1 and 2: NOTA TÉCNICA PROGRAMA COMPUTACIONA
Page 3: obtidas do Instituto Nacional de Me
Page 7 and 8: Caso o usuário prefira a utilizaç
Page 9 and 10: (a) (b) Figura 5. Telas do CLIMABR
Page 11: SEDIYAMA, G.C. et al. Simulação d

NOTA TÃCNICA - UFV

Create successful ePaper yourself

Delete template?

Save as template?

NOTA TÃCNICA - UFV