Parametrização das equações de Hargreaves & Samani e Angströn-Prescott para estimativa da radiação solar na região de Telêmaco Borba, estado do Paraná

Parameterization of Hargreaves & Samani and Angstrom-Prescott equations to estimate solar radiation in the region of Telemaco Borba, Parana state, Brazil

Parametrização das equações de Hargreaves & Samani e Angströn-Prescott para estimativa da radiação solar na região de Telêmaco Borba, estado do Paraná

Parameterization of Hargreaves & Samani and Angstrom-Prescott equations to estimate solar radiation in the region of Telemaco Borba, Parana state, Brazil

Daniela Jerszurki^I; Jorge Luiz Moretti de Souza^II,1

^IPós-graduação em Ciência do Solo, Universidade Federal do Paraná (UFPR), Curitiba, PR, Brasil

^IIDepartamento de Solos e Engenharia Agrícola, UFPR, 80035-050, Curitiba, PR, Brasil. E-mail: [e-mail]

RESUMO

A estimativa da radiação solar incidente (Rs) por equações matemáticas é útil quando não se dispõe de estações meteorológicas que a registrem. O objetivo do trabalho foi estimar, para os períodos anual, sazonal e mensal, o coeficiente empírico Krs da equação de Hargreaves & Samani e os coeficientes a e b da equação de Angströn-Prescott para estimar a Rs, confrontando valores estimados com valores observados para o município de Telêmaco Borba, Estado do Paraná. Utilizaram-se dados históricos diários (janeiro de 1981 a janeiro de 2011) de temperatura máxima, mínima e média do ar, radiação solar global acumulada, insolação diária e radiação solar incidente no topo da atmosfera (Ra), calculada com metodologia recomendada pela Food and Agriculture Organization of the United Nations - FAO. Ocorre maior associação entre Rs estimada e observada com valores de Krs anual e sazonal. Os coeficientes a e b apresentaram índices "c" "ótimo" ("c" > 0,85) para estimar Rs nos períodos analisados. Os resultados sugerem que os coeficientes encontrados podem ser utilizados para estimar valores de Rs na região.

Palavras-chave: coeficiente empírico, estimativa, radiação solar.

ABSTRACT

The estimation of solar radiation (Rs) by mathematical equations is useful when there are no weather stations that record. The objective of this study was to estimate, for the periods annual, seasonal and monthly, the empirical coefficient Krs of Hargreaves & Samani equation and the coefficients a and b of Angstrom-Prescott equation for estimating Rs, comparing estimated values with observed values for Telemaco Borba, Paraná. Used historical data daily (January 1981 to January 2011) maximum, minimum and average air temperature, cumulative global solar radiation, daily insolation and extraterrestrial radiation (Ra), calculated using a methodology recommended by the Food and Agriculture Organization of the United Nations- (FAO). There is greater association between estimated and observed Rs with values of annual and seasonal Krs. The coefficients a and b show indices "c" "great" ("c">0.85) to estimate Rs in the periods analyzed. The results suggest that the coefficients found are used to estimate values of the region.

Key words: empirical coefficient, estimate, solar radiation.

INTRODUÇÃO

A radiação solar tem grande influência sobre os processos atmosféricos e terrestres, alterando a temperatura e a circulação atmosférica (SOUZA et al., 2008). Segundo CARGNELUTTI FILHO et al. (2007), sua medida em um determinado local e período do ano é imprescindível para estudos ambientais, possibilitando a modelagem meteorológica e climática. Além disso, é utilizada em estudos de necessidade hídrica de culturas irrigadas, modelagem de crescimento, mudanças climáticas, dentre outras (BORGES et al., 2010). Apesar disso, existem poucas estações meteorológicas que registram a radiação solar, fazendo com que a disponibilidade de dados medidos seja reduzida (POLO et al., 2011). No Brasil, o Instituto Nacional de Meteorologia (INMET) dispõe de 321 estações meteorológicas convencionais. Dessas, 232 medem a insolação e somente 58 (18,1%) a radiação solar incidente (Rs), por meio de actinógrafos (LIMA, 2005). Além da falta de estações que estimem esses parâmetros, as que existem estão centralizadas nas regiões mais desenvolvidas do país, enquanto há grandes áreas agrícolas que não dispõem de dados climáticos completos.

Na tentativa de eliminar o problema, vários modelos matemáticos foram desenvolvidos para a estimativa da Rs. No entanto, frequentemente, necessitam de dados climáticos, que nem sempre estão disponíveis nas estações meteorológicas, que medem principalmente temperatura e precipitação, e necessitam também de calibração quando utilizados em condições diferentes das que foram desenvolvidos (BORGES et al., 2010). Dentre os diversos modelos matemáticos que estimam a radiação solar incidente (Rs), destaca-se a equação desenvolvida por HARGREAVES & SAMANI (1982), baseada na amplitude térmica diária, radiação solar no topo da atmosfera (Ra) e no coeficiente de ajuste (Krs). Na equação, a amplitude térmica diária pode ser usada como um indicador da fração da Ra que chega à superfície do solo, ou seja, a Rs.

Para localidades do interior, onde dominam as grandes extensões de terra e as massas de ar não são fortemente influenciadas por um grande corpo de água, o valor do Krs recomendado é de 0,16. No entanto, para localidades costeiras, onde as massas de ar estão influenciadas por uma massa de água nas proximidades, o Krs recomendado é de 0,19 (ALLEN et al., 1998). No entanto, o Krs pode ser melhor estimado para as condições de uma determinada localidade, refletindo melhores estimativas da Rs. Outra equação, amplamente utilizada é a de Angströn-Prescott (PRESCOTT, 1940), que estima a Rs a partir de medidas do número efetivo de horas de brilho solar. ANGSTRÖN (1924) apresentou uma equação em que o quociente entre Rs e Ra foi linearmente correlacionado com o quociente do número de horas de brilho solar (n) pela insolação máxima possível (N). PRESCOTT (1940) simplificou a equação para que os coeficientes linear (a) e angular (b) pudessem ser obtidos a partir de ajustes estatísticos, sendo que o método foi denominado de Angströn-Prescott. A equação de Angströn-Prescott pode ser utilizada tanto para estimar a Rs, quanto para obter n em estações automáticas que normalmente não medem esse valor. A Food and Agriculture Organization of the United Nations (FAO) (ALLEN et al., 1998) apresenta os valores dos coeficientes a (0,25) e b (0,50) considerados padrão, no entanto, esses coeficientes podem ser estimados regionalmente, refletindo-se em valores mais confiáveis de Rs.

O objetivo do presente trabalho foi estimar o coeficiente empírico Krs da equação de Hargreaves & Samani, bem como os coeficientes a e b da equação de Angströn-Prescott para estimar a Rs, confrontando os valores estimados com os valores observados para o município de Telêmaco Borba, Estado do Paraná.

MATERIAL E MÉTODOS

As análises do presente trabalho foram realizadas para o município de Telêmaco Borba, Estado do Paraná. A região está localizada no Segundo Planalto Paranaense e, segundo a classificação de Köeppen, apresenta clima tipo Cfa/Cfb, subtropical úmido transicional para temperado propriamente dito, com temperatura média no mês mais frio inferior a 16°C, com ocorrência de geadas e temperatura média no mês mais quente acima de 22°C, com verões quentes (IAPAR, 2000). Foram utilizados dados históricos diários de temperatura máxima, mínima e média do ar (°C), radiação solar global acumulada (MJ m^-2 dia^-1) e insolação diária (h dia^-1), de janeiro de 1981 a janeiro de 2011, provenientes de uma estação climatológica instalada na Fazenda Monte Alegre, pertencente à empresa Klabin Florestal, localizada na região de Telêmaco Borba, Estado do Paraná, com altitude média de 885 m, nas coordenadas 24°13' de latitude Sul e 50°32' de longitude Oeste. A radiação solar incidente no topo da atmosfera (Ra) foi calculada de acordo com metodologia apresentada por ALLEN et al. (1998), a qual é baseada na latitude do local e no dia juliano. As estimativas do coeficiente Krs foram realizadas empregando-se a equação de HARGREAVES & SAMANI (1982).

Sendo: Rs - radiação solar incidente (MJ m^-2 dia^-1); Krs - coeficiente empírico (°C^-0,5); T_max - temperatura máxima do ar (°C); Tmin - temperatura mínima do ar (°C); Ra - radiação solar incidente no topo da atmosfera (MJ m^-2 dia^-1). Para as estimativas dos coeficientes a e b, foi utilizada a equação de Angströn-Prescott (PRESCOTT, 1940).

Sendo: Rs - radiação solar incidente (MJ m^-2 dia^-1); Ra - radiação solar incidente no topo da atmosfera (MJ m^-2 dia^-1); n - insolação diária (h dia^-1); N - duração máxima teórica do dia (h dia^-1); a - coeficiente linear; b - coeficiente angular.

Os coeficientes a e b da equação Angströn-Prescott, para os períodos anual, sazonal e mensal, foram obtidos a partir de análise estatística de regressão linear, estimando Rs a partir de n. Para verificar a correlação entre os resultados de Rs observados e estimados, foram realizadas análises de regressão linear com a obtenção dos seus respectivos coeficientes (R). Para avaliar o grau de concordância entre as estimativas e os valores observados de Rs, foi utilizado o índice "d" de WILLMOTT et al. (1985):

Sendo: d - índice de concordância de WILLMOTT et al. (1985); Rs_i - radiação solar incidente estimada no i-ésimo dia (MJ m^-2 dia^-1); Rs_oi - radiação solar incidente observada no i-ésimo dia (MJ m^-2 dia^-1); - média da radiação solar incidente observada no período considerado (MJ m^-2 dia^-1).

Para comparação dos valores observados e estimados, foi utilizada uma adaptação do índice de desempenho "c" adotado por CAMARGO & SENTELHAS (1997):

Sendo: c - índice "c" de desempenho proposto por CAMARGO & SENTELHAS (1997) (adimensional); R - coeficiente de correlação (adimensional) obtido na regressão linear (adimensional); d - índice "d" de concordância proposto por WILLMOTT et al. (1985) (adimensional).

O critério de interpretação do índice "c" tem os seguintes desempenhos: ótimo ("c">0,85); muito bom (0,75<"c" ≤0,85); bom (0,65<"c" ≤0,75); mediano (0,60<"c" ≤0,65); sofrível (0,50<"c" ≤0,60); mau (0,40<"c" ≤0,50); e péssimo ("c" ≤0,40).

RESULTADOS E DISCUSSÃO

Observou-se que os valores do coeficiente empírico Krs ficaram entre 0,11 e 0,12 ao longo das estações e meses do ano (Tabela 1). Os resultados discordam do valor de 0,16 recomendado por ALLEN et al. (1998) para regiões do interior. O índice "d" de WILLMOTT et al. (1985), considerando a estimativa da Rs com o coeficiente Krs, foi de 0,91 para o período anual, variando de 0,79 a 0,87 ao longo do ano. Comparando-se os valores observados e estimados de Rs, as análises indicaram índices de desempenho "c" "mediano" (0,60<"c" ≤0,65) a "muito bom" (0,75<"c" ≤0,85) para os períodos anual, sazonal e mensal (Tabela 2). Pode-se inferir que houve maior associação entre a Rs estimada e observada quando se utilizaram valores de Krs anual e sazonal (Tabela 2). Os piores resultados ocorreram para o mês de fevereiro, provavelmente, devido à dispersão que ocorre em períodos chuvosos, quando os valores de Rs medidos são menores do que os calculados. Entretanto, a estimativa de Rs pode apresentar melhor desempenho quando são empregadas metodologias que utilizem dados de temperatura e precipitação pluvial (PODESTÁ et al., 2004; WEISS & HAYS, 2004).

As análises estatísticas evidenciaram que, em média, o parâmetro "a" foi igual a 0,19 para o período anual, apresentando variações ao longo das estações e meses do ano (Tabela 3). Os maiores valores de "a" para os períodos verão e outono evidenciaram que a transmissividade da atmosfera é maior nos meses de inverno e primavera, devido à menor quantidade de poeira no ar e à frequente entrada de ar polar que apresentam menor turbulência em relação às massas tropicais (BLANCO & SENTELHAS, 2002). Os maiores valores de b ocorreram no inverno e os menores no verão (Tabela 3), devido também ao efeito da transmissividade atmosférica. BLANCO & SENTELHAS (2002), estimando os coeficientes a e b da equação de Angströn-Prescott para o município de Piracicaba - SP, em um período de quatro anos (junho de 1996 a maio de 2000), utilizando dados de radiação solar e brilho solar, na escala anual, mensal e sazonal, encontraram valores anuais para a=0,23 e b=0,50, semelhantes ao encontrado no presente trabalho para o coeficiente a, que variou de 0,17 a 0,23. No entanto, diferiram do encontrado para o coeficiente b que variou de 0,35 a 0,45. Em estudo semelhante, considerando uma série de um ano (setembro de 1966 a agosto de 1967), OMETTO (1968) encontrou valores para a e b iguais a 0,26 e 0,51, respectivamente. DANTAS et al. (2003), estimando os coeficientes a e b para o estado de Minas Gerais, também encontraram valores para os coeficientes a=0,23 e b=0,49 semelhantes ao encontrado para o coeficiente a, mas diferindo dos valores encontrados para o coeficiente b, para o presente trabalho.

O valor médio anual do coeficiente a determinado no presente trabalho foi 11% menor que o encontrado por OMETTO (1968), enquanto que o valor médio do coeficiente b foi 18% e 20% menor do que o apresentado, respectivamente, por BLANCO & SENTELHAS (2002) e OMETTO (1968). Os valores de a e b obtidos para Telêmaco Borba são menores, e os resultados podem estar relacionados ao maior período de anos (30 anos) utilizados e à diferença entre os sensores que medem a Rs. Os resultados obtidos evidenciam as considerações de CARGNELUTTI FILHO et al. (2004), quando afirmam que há uma relação inversamente proporcional entre o tamanho da amostra utilizada (número de anos) com os coeficientes a e b e, consequentemente, com a Rs.

Considerando os valores de R (Tabela 4), pode-se inferir que existe grande associação entre a Rs estimada com os coeficientes a e b (Tabela 3) e a Rs observada para o período estudado. O índice "d" de WILLMOTT et al. (1985), considerando a estimativa da Rs com os coeficientes a e b, foi de 0,95 para o período anual, variando de 0,92 a 0,95 ao longo do ano. Comparando-se os valores observados e estimados da Rs, as análises evidenciaram índices de desempenho "c" "ótimos" ("c">0,85) para os períodos anual, sazonal e mensal (Tabela 4), indicando precisão de estimativa dos valores de Rs. Aplicando-se os valores de a e b, obtidos nos períodos anual, sazonal e mensal na equação de Angströn-Prescott para a estimativa de Rs, pode-se afirmar estatisticamente que não houve diferença entre valores estimados e observados (PEREIRA et al., 2002; DORNELAS et al., 2006), sugerindo que os coeficientes a e b podem ser utilizados na estimativa da Rs, para o município de Telêmaco Borba.

CONCLUSÃO

O valor médio de Krs (0,12) está abaixo do valor proposto pela bibliografia, provavelmente devido às diferenças entre os sensores que registram a Rs nas estações meteorológicas. Ocorre maior associação entre a Rs estimada e observada quando são utilizados valores de Krs anual e sazonal, resultando em índices "c", classificados de muito bom (0,75<"c" ≤0,85) a bom (0,65<"c" ≤0,75). Os coeficientes a e b apresentaram índices de desempenho "c", classificados como "ótimo" ("c">0,85) para estimar a Rs de todos os períodos analisados. Os coeficientes encontrados para as equações de Angströn-Prescott e Hargreaves-Samani podem ser utilizados com confiança para estimar valores de Rs na região de Telêmaco Borba-PR.

Recebido para publicação 03.05.12

Aprovado em 16.09.12

Devolvido pelo autor 12.11.12

CR-2012-0108.R1

• ALLEN R.G. et al. Crop evapotranspiration: guidelines for computing crop water requirements. Rome: FAO, 1998. 356p. (FAO: Irrigation and Drainage Paper, 56).
• ANGSTRÖN, A. Solar and terrestrial radiation. Quarterly Journal of the Royal Meteorological Society, v.50, p.121-126, 1924. Disponível em: <http://onlinelibrary.wiley.com/doi/10.1002/qj.49705021008/abstract>. Acesso em: 26 mar. 2012. doi: 10.1002/qj.49705021008.
• BLANCO, F.F.; SENTELHAS, P.C. Coeficientes da equação de Angströn-Prescott para estimativa da insolação para Piracicaba. Revista Brasileira de Agrometeorologia, v.10, p.295-300, 2002. Disponível em: <http://dx.doi.org/10.1590/S1415-43662011000800011>. Acesso em: 26 mar. 2012. doi: 10.1590/S1415-43662011000800011.
• BORGES V.P. et al. Avaliação de modelos de estimativa da radiação solar incidente em Cruz das Almas, Bahia. Revista Brasileira de Engenharia Agrícola e Ambiental, v.14, p.74-80, 2010. Disponível em: <http://www.scielo.br/pdf/rbeaa/v14n1/v14n01a10.pdf>. Acesso em: 26 mar. 2012. doi: 10.1590/S1415-43662010000100010.
• CAMARGO, A.P. de; SENTELHAS, P.C. Avaliação do desempenho de diferentes métodos de estimativa da evapotranspiração potencial no Estado de São Paulo, Brasil. Revista Brasileira de Agrometeorologia, v.5, p.89-97, 1997.
• CARGNELUTTI FILHO, A. et al. Tamanho de amostra para a estimativa das médias decendiais de radiação solar global no estado do Rio Grande do Sul. Ciência e Agrotecnologia, v.31, p.1402-1410, 2007. Disponível em: <http://dx.doi.org/10.1590/S1413-70542007000500020>. Acesso em: 26 mar. 2012. doi: 10.1590/S1413-70542007000500020.
• DANTAS, A. et al. Estimativa da radiação solar global para a região de Lavras, MG. Ciência e Agrotecnologia, v.27, p.1260-1263, 2003. Disponível em: <http://dx.doi.org/10.1590/S1413-70542003000600008>. Acesso em: 26 mar. 2012. doi: 10.1590/S1413-70542003000600008.
• DORNELAS, K.D.S. et al. Coeficientes médios da equação de Angström-Prescott, radiação solar e evapotranspiração de referência em Brasília. Pesquisa Agropecuária Brasileira, v.41, p.1213-1219, 2006. Disponível em: <http://dx.doi.org/10.1590/S0100-204X2006000800001>. Acesso em: 26 mar. 2012. doi: 10.1590/S0100-204X2006000800001.
• HARGREAVES, G.H.; SAMANI, Z.A. Estimating potential evapotranspiration. Journal of Irrigation and Drainage Engineering, v.108, p.225-230, 1982. Disponível em: <http://cedb.asce.org/cgi/WWWdisplay.cgi?35047>. Acesso em: 26 mar. 2012.
• INSTITUTO AGRONÔMICO DO PARANÁ. Cartas climáticas do Estado do Paraná Londrina: IAPAR, 2000. 1v. Disponível em: <http://www.iapar.br/modules/conteudo/conteudo.php?conteudo=677>. Acesso em: 03 set. 2012.
• LIMA, E.P. Evapotranspiração de referência de Penman-Monteith, padrão FAO (1998), a partir de dados de temperaturas máxima e mínima de Minas Gerais 2005. 67f. Dissertação (Mestrado em Meteorologia Agrícola) - Curso de Pós-graduação em Meteorologia Agrícola, Universidade Federal de Viçosa, MG.
• OMETTO, A. Estudo das relações entre radiação solar global, radiação líquida e insolação 1968. 64f. Tese (Doutorado em Agronomia) - Curso de Pós-graduação em Agronomia, Universidade de São Paulo, Escola Superior de Agricultura "Luiz de Queiroz", SP.
• PEREIRA, A.B. et al. Estimativa da radiação solar global diária em função do potencial de energia solar na superfície do solo. Scientia Agricola, v.59, n.2, p.211-216, 2002. Disponível em: <http://dx.doi.org/10.1590/S0103-90162002000200002>. Acesso em: 03 set. 2012.
• PODESTÁ, G.P. et al. Estimating daily solar radiation in the Argentine Pampas. Agricultural and Forest Meteorology, v.123, p.41-53, 2004. Disponível em: <http://dx.doi.org/10.1016/j.agrformet.2003.11.002>. Acesso em: 26 mar. 2012. doi: 10.1016/j.agrformet.2003.11.002.
• POLO, J. et al. A simple approach to the synthetic generation of solar irradiance time series with high temporal resolution. Solar Energy, v.85, p.1164-1170, 2011. Disponível em: <http://www.sciencedirect.com>. Acesso em: 30 out. 2012. doi:10.1016/j.solener.2011.03.011.
• PRESCOTT, J.A. Evaporation from a water surface in relation to solar radiation. Transactions of the Royal Society Science Australian, v.64, p.114-118, 1940.
• SOUZA, J.D. et al. Estimativa da radiação solar global à superfície usando um modelo estocástico: caso sem nuvens. Revista Brasileira de Geofísica, v.2, p.1-36, 2008. Disponível em: <sid.inpe.br/mtc-m17@80/2008/04.14.19.06v22008-04-16>. Acesso em: 26 mar. 2012.
• WEISS, A.; HAYS, C.J. Simulation of daily solar irradiance. Agricultural and Forest Meteorology, v.123, p.187-199, 2004. Disponível em: <http://dx.doi.org/10.1016/j.agrformet.2003.12.002>. Acesso em 26 mar. 2012. doi: 10.1016/j.agrformet.2003.12.002.
• WILLMOTT, C.J. et al. Climatology of terrestrial seasonal water cicle. Internation Journal of Climatology, v.5, p.589-606, 1985. Disponível em: <http://climate.geog.udel.edu/~climate/publication_html/Pdf/WRM_JClimatology_1985.pdf>. Acesso em: 26 mar. 2012.

1

Autor para correspondência.

Datas de Publicação

• Publicação nesta coleção
  01 Fev 2013
• Data do Fascículo
  Mar 2013