Print version ISSN 0102-7786
Rev. bras. meteorol. vol.26 no.3 São Paulo Sept. 2011
Downward longwave radiation estimates for clear-sky conditions over northeast Brazil
Estimativas da radiação de onda longa incidente para condições de céu claro no nordeste do Brasil
Carlos Antonio Costa dos SantosI; Bernardo Barbosa da SilvaI; Tantravahi Venkata Ramana RaoI; Prakki SatyamurtyII; Antonio Ocimar ManziIII
IUniversidade Federal de Campina Grande, Unidade Acadêmica de Ciências Atmosféricas (UFCG/UACA), Campina Grande, PB, Brasil
IIUniversidade do Estado do Amazona (CESTU), Manaus, AM, Brasil
IIIInstituto Nacional de Pesquisas da Amazônia (INPA), Manaus, AM, Brasil
firstname.lastname@example.org, email@example.com, firstname.lastname@example.org, email@example.com, firstname.lastname@example.org
The main objective of this paper is to assess the performance of nine downward longwave radiation equations for clear-sky condition and develop a locally adjusted equation using the observed vapor pressure and air temperature data. The radiation and atmospheric parameters were measured during the months of October 2005 to June 2006 at a micrometeorological tower installed at the experimental site in a banana orchard in the semiarid region of Northeast Brazil. The comparative statistics for the performance of the downward longwave radiation calculation models during daytime and nighttime compared to measured data have shown that the parameterizations with more physical foundations have the best results. The locally adjusted equation and Sugita and Brutsaert model developed in 1993 showed errors less than 1.0% in comparison with measured values. Downward longwave radiation is one of the most expensive and difficult component of the radiation budget to be monitored in micrometeorological studies. Hence, the locally adjusted equation can be used to estimate downward longwave energy, needed as input to some agricultural and hydrological models, in semi-arid regions of the Northeast Brazil, where this component is not monitored.
Keywords: Radiation balance, banana crop, semiarid, evapotranspiration
O principal objetivo desse artigo é avaliar o desempenho de nove equações, para obter radiação de onda longa incidente para condições de céu claro e desenvolver uma equação localmente ajustada utilizando dados observados de pressão de vapor e temperatura do ar. Os parâmetros atmosféricos e de radiação foram medidos durante os meses de outubro/2005 e junho/2006, em uma torre micrometeorológica instalada no sítio experimental, em um pomar de bananeiras na região semi-árida do Nordeste do Brasil. A estatística comparativa para o desempenho dos modelos, para o cálculo da radiação de onda longa durante os períodos diurnos e noturnos, comparados aos dados medidos, tem mostrado que as parametrizações com mais fundamentações físicas têm melhores resultados. A equação ajustada localmente e o modelo de Sugita e Brutsaert, desenvolvido em 1993, mostraram erros menores que 1.0%, em comparação com os valores medidos. A radiação de onda longa incidente é um dos componentes do balanço de radiação mais caro e difícil de ser monitorado em estudos micrometeorológicos. Logo, a equação localmente ajustada pode ser usada para estimar a energia de onda longa incidente, necessária como dado de entrada em alguns modelos agrícolas e hidrológicos, nas regiões semi-áridas do Nordeste do Brasil, onde esta componente não é monitorada.
Palavras-Chave: Balanço de Radiação, banana, semi-árido, evapotranspiração
Solar radiation () is absorbed and reflected by the Earth´s surface and atmosphere (including clouds), and part of the absorbed shortwave radiation can be emitted back as longwave radiation, additionally, the emitted longwave radiation can as well be reflected and absorbed by Earth´s surface and atmosphere (Oke, 1996). Many biological and physical processes are controlled by the amounts of solar energy received, reflected and emitted by the Earth and the atmosphere.
A measure of the quantity of energy available for processes such as evapotranspiration is often needed to be determined in agricultural, biological and engineering sciences, and in many other studies, such as climate modeling and forecasting, water resources planning and plant growth modeling (Sridhar and Elliott, 2002; Kjaersgaard et al., 2007). The methods for calculating the energy available at an evaporating surface are, in general, derived from the energy balance and for which the net radiation has to be known. Net radiation is sometimes measured. The accuracy of the net radiation term obtained from the radiation balance depends on the quality of the measurements or estimations of other terms used as input (Kjaersgaard et al., 2007).
Considerable efforts have been made to develop simplified methods to calculate atmospheric radiation. In spite of recent advances in measurement technology, there is considerable interest in determining the downward longwave radiation (Ld), which is emitted by the atmosphere and clouds, from the available meteorological measurements such as air temperature and vapor pressure near the ground surface (Sugita and Brutsaert, 1993). The studies of surface energy balance are often complicated due to the various uncertainties in the estimation of downward longwave radiation at the surface (Crawford and Duchon, 1999). Accurate estimates of surface radiation balance, which modulates the magnitude of the surface energy balance components (including the latent heat flux), depend on good estimates of downward longwave radiation.
Instrument that measures directly the downward longwave radiation is the pyrgeometer, which is expensive and sensitive in comparison with the pyranometer which is used to measure the incoming shortwave radiation (Duarte et al., 2006). Downward longwave radiation (Ld) is the surface radiation balance component rarely available at meteorological stations (Kruk et al., 2010). Thus, in recent decades many techniques have been developed to estimate Ld where there are no pyrgeometers, based only on surface observations. These methods have presented varied accuracy, and new techniques continue to be developed on a regional scale (Crawford and Duchon, 1999; Kruk et al., 2010).
Based on a perceived similarity between heat conduction and radiative transfer, Brunt (1932) first theorized that Ld was related to the square root of surface vapor pressure (ea) and developed the first empirical relationship between these quantities, resulting in a correlation coefficient of 0.97 between observations and estimates. Swinbank (1963) concluded that Ld was not related to ea at all, but to the square of the temperature (T), resulting in a correlation coefficient of 0.99. The equation developed by Idso and Jackson (1969), also dependent on T2, produced a correlation coefficient of 0.99 between measured and estimated Ld. However, Brutsaert (1975) was the first to develop a more physical parameterization of atmospheric emissivity (εa), based on the equation of Schwarzschild (Liou, 1980) and assumptions of standard atmospheric lapse rates of T and ea.
In Brazil, some studies have been made to analyze those proposed equations and to adjust the equations for clear-sky conditions with observed data. For example, Duarte et al. (2006) assessed the performance of Ld equations at an experimental site near the city of Ponta Grossa in Paraná State, Southern Brazil. Kruk et al. (2010) adjusted the equations for Ld using data observed in Sertãozinho, São Paulo State, Southeast Brazil. However, no study has been developed to assess the equations and to calibrate locally the parameters for the models in the Northeast Brazil. Thus, the objective of this study is to assess the performance of downward longwave radiation equations for clear-sky conditions and to develop an equation using the observed data at an experimental site in a banana orchard for the semiarid region of Northeast Brazil.
2. MATERIALS AND METHODS
The radiation parameters and atmospheric parameters were measured at every 10 min at the micrometeorological tower installed at the Frutacor farm (5º08'44"S; 38º05'53"W; altitude 147 m), within an approximate area of 250 ha of banana crop (Musa sp.), during the months of October 2005 to June 2006. This farm is located within the irrigation district of Quixeré, in the Lower Jaguaribe basin, in Ceará state, Northeast Brazil. The area has a semiarid climate, type BSw'h', according to the Köeppen classification. The mean air temperature is 28.5 ºC, with maximum and minimum values of 36 ºC and 22 ºC, respectively, and mean relative humidity of 62%. The mean total annual precipitation is 772 mm.
The radiation balance equation at the surface is given by
where Ŝd is the incident solar radiation flux, Sref is the reflected solar radiation flux which is related to albedo, Ld is the atmospheric longwave radiation flux or downward longwave radiation and Lsurface is the surface longwave radiation flux which is related to the Earth's surface temperature and albedo as well (Tanaka et al., 2008). The net radiation (Rn) was computed using a four component net radiometer (model CNR1, Kipp & Zonnen, Delf, The Netherlands), placed at an elevation of 7.0 m above the ground level and 2.0 m above a row of plants. Downward longwave radiation was measured with a CG3 pyrgeometer contained in the CNR1 instrument. Microclimatic data of air temperature and relative humidity were measured using a Vaisala HMP45C probe. All sensors signals have been sampled at 1 minute intervals and averaged at 10 min intervals and were recorded using a Campbell Scientific CR23X datalloger, where the data were stored for processing later. This study was developed with a new net radiation sensor and used calibration constants provided by the manufacturer. The experimental station was continuously powered by four batteries and four solar panels.
The equation for the clear-sky Ld has the general form:
where T is the air temperature near-surface, εa is the emissivity of the clear-sky atmosphere and σ = 5.67051 x 10-8 W m-2 K-4 is the Stefan-Boltzmann constant. According to Duarte et al. (2006) and Kruk et al. (2010) Equation (2) is inspired by Stefan-Boltzmann's law for the radiation emitted by a grey body at uniform temperature. All the equations analyzed in this study make the assumption that εa is a function of temperature and/or vapor pressure near the ground.
Several authors have proposed models for calculating εa. Six popular εa formulations and those developed by Duarte et al. (2006) and Kruk et al. (2010) inserted into Equation (2) were tested against Ld data collected on the experimental field. To analyze the consistency of the clear-sky models, proposed by different authors are presented in Table 1. A sample of 12 clear-sky days was selected. These days presented fractional cloud coverage (c), calculated from the Equation 3 (Crawford and Duchon, 1999), less than or equal to 0.05.
where Sd is the observed downward solar (shortwave) radiation at the surface and is the theoretical downward clear-sky solar radiation for the period obtained as in Crawford and Duchon (1999).
Statistical criteria for the evaluation of the models are given in Table 2. For each parameterization, the estimated value of Ld, which is denoted by e, were compared with the measured values, which is denoted by m. The statistical analysis includes bias (BIAS), root mean square error (RMSE), mean absolute error (MAE) and percent mean relative error (PMRE).
3. RESULTS AND DISCUSSION
The data needed for testing values of Ld (downward longwave radiation) using Equations 4 - 11 given in Table 1, included near-surface observations of downward shortwave radiation (Ŝd), vapor pressure and temperature. Initially, an equation (Equation 12) has been adjusted for clear-sky conditions (c < 0.05) using the measured vapor pressure and temperature data from the experimental site in Quixeré, Ceará. The locally calibrated equation is:
The least square linear regression obtained between estimated downward longwave radiation fluxes and observed radiation fluxes in clear-sky days are presented in Figure 1. The Equations 4 and 5, which utilize air temperature only do not employ the vapor pressure as a measure of the humidity effect on atmospheric path length, presented slope values of 1.64 and 1.70, respectively. These results indicate overestimation as showed in Figures 1a-1b. They are in agreement with those presented in Duarte et al. (2006). The other parameterizations (Equations 6 - 11) presented underestimation. The highest slope was obtained with the parameterization proposed by Sugita and Brutsaert (1993) with a = 0.95.
In Figure 2 the comparisons between downward longwave radiation fluxes estimated by Equation 12, i.e. locally calibrated equation, and the observed downward longwave radiation fluxes on clear-sky days. The least square linear regression presented slope value of 0.92 indicating underestimation and showing similar results when compared with those presented by Brutsaert (1975), Sugita and Brutsaert (1993) and Prata (1996) models (see Figure 1).
Table 3 shows the comparative statistics for the performance of the nine downward longwave clear-sky calculation models during daytime compared to measured data from the experimental site in Northeast Brazil. A 10 min input dataset was used with a total of 460 observations presenting mean value of 389.3 W m-2. Overall, those parameterizations, which use vapor pressure and air temperature to obtain the atmospheric emissivity, have the best scores during daytime. The best results, which resulted in the smallest BIAS, RMSE, MAE and PMRE were obtained by the Equation 12 followed by the Sugita and Brutsaert (1993) model. The results with the highest BIAS, RMSE, MAE and PMRE were presented by the Idso and Jackson (1969) and Swinbank (1963), respectively. These results presented by Equation 12 are in agreement with those obtained by Duarte et al. (2006) and Kruk et al. (2010).
Equations proposed by Sugita and Brutsaert (1993), Prata (1996), Duarte et al. (2006) and Equation 12 showed the highest correlation coefficient (r = 0.94), showing the strong correlation between the calculated and measured values of Ld. All models presented mean relative percentage error less than 4.0% indicating good accuracy of the models. However, the adjusted Equation 12 and Sugita and Brutsaert (1993) model showed errors less than 1.0% in comparison with measured values. These results are better than those found by Duarte et al. (2006), Kjaersgaard et al. (2007) and Kruk et al. (2010).
The comparative statistics of the nine downward longwave clear-sky calculation models during nighttime compared to measured data is shown in Table 4. The 10 min input dataset with a total of 596 observations and mean value of 365.2 W m-2 is used for this purpose. During the nighttime, the best results were obtained by the Equation 12, resulting in the smallest RMSE, MAE and PMRE followed by the Sugita and Brutsaert (1993) model. Idso (1981) and Kruk et al. (2010) models showed results with highe BIAS of 17.4 W m-2 and 7.2 W m-2, respectively. The higher r value equal to 0.94 has been presented by the Equation 12, Prata (1996) and Duarte et al. (2006) showing a strong correlation between the calculated and measured values of downward longwave radiation. PMRE values are less than 5.0% showing a good accuracy of the parameterizations and the adjusted Equation 12 and Sugita and Brutsaert (1993) model showed errors less than 1.0% in comparison with those measured. Again, the results have been better than those found by Duarte et al. (2006), Kjaersgaard et al. (2007) and Kruk et al. (2010). Surprisingly, the parameterizations with less physical foundations (Swinbank, 1963 and Idso and Jackson, 1969) had better scores during nighttime in comparison with those presented in the daytime analysis. Stull (1988) affirms that if there is strong convection in the daytime, both lapse rates above the surface layer should be close to constant through the convective boundary layer (CBL), while in the nighttime under stable conditions this idea might not be valid especially for temperature lapse rate.
Figure 3 shows the mean diurnal cycle of Ld for clear-sky days calculated, using Equations 2 and 12, and field measurements. The behavior of calculated values of Ld, between midnight and early morning, is smaller than the measured values. In the early hours after the sunrise the model overestimated the values. However, between 10:00 AM and 15:00 PM the calculated values are underestimated and after that time they are overestimated until the sunset. Meanwhile, it can be observed that the average of Ld calculated for clear-sky days from Equation 12 closely follows the average of observed values (Figure 3), showing a good agreement between them.
The results found in this research will contribute to a better understanding of the surface energy budget in the semi-arid of Northeast Brazil and it can results in improved characterization of evapotranspiration for many applications in agriculture, biology, ecology, climatology and hydrology (Choi et al., 2008). The Ld estimates are vitally important for numerous applications requiring surface radiation and energy balance, including prediction evapotranspiration, surface temperature and frost occurrence (Flerchinger et al., 2009). Solar radiation is fundamental for crop growth, but there are many other factors that influence plant development. For example, water availability for the plants is one of the most important factors for their development, and under some optimal conditions of water availability, the radiation use remains almost constant throughout most of the plant growth cycle, and shows no variations caused by local atmospheric conditions.
These models can be used as research tools to help improve current understanding of the basic physiology of crop growth and development, and as decision support tools to help optimize crop and soil management strategies. Most crop simulation models, however, require daily records of temperature, precipitation, and radiation. Records of daily maximum and minimum temperatures, and precipitation, are available for many sites, but accurate solar radiation measurements are available for only a few sites, a situation which probably reflects not only the cost of the required instrumentation, but also the problems of maintenance and calibration of radiation sensors.
The equations which utilize air temperature only, i.e. Swinbank (1963) and Idso and Jackson (1969), overestimate the atmospheric longwave radiation while other parameterizations present underestimation. The comparative statistics for the performance of the nine downward longwave clear-sky calculation models during daytime and nighttime compared to measured data have shown that the parameterizations with more physical foundations have the best scores. The results obtained by the locally adjusted equation presented the smallest BIAS, RMSE, MAE and PMRE. The adjusted Equation 12 and Sugita and Brutsaert (1993) model showed errors less than 1.0% in comparison with measured values during the daytime and nighttime.
Since downward longwave radiation is one of the most difficult component of the radiation budget to be monitored in micrometeorological studies, it is important to evaluate a locally adjusted equation to estimate it. Hence, the locally adjusted equation (Equation 12) can be used to estimate downward longwave energy, needed as input to some agricultural and hydrological models, in semi-arid regions of the Northeast Brazil, where this component is not monitored.
The authors are grateful for the National Council for Scientific and Technological Development (CNPq) and thank the anonymous reviewers whose suggestions were helpful for improving the quality of this paper.
BRUNT, D. Notes on radiation in the atmosphere. Quarterly Journal of the Royal Meteorological Society. v. 58, p. 389-418, 1932. [ Links ]
CHOI, M., JACOBS, J. M., KUSTAS, W. P. Assessment of clear and cloudy sky parameterizations for daily downwelling longwave radiation over different land surfaces in Florida, USA. Geophysical Research Letters, v. 35, L20402, doi: 10.1029/2008GL035731, 2008 [ Links ]
CRAWFORD, T. M., DUCHON, C. E. An improved parameterization for estimating effective atmospheric emissivity for use in calculating daytime downwelling longwave radiation. Journal of Applied Meteorology, v. 38, p. 474-480, 1999 [ Links ]
CULF, A. D., GASH, J. H. C. Longwave radiation from clear skies in Niger: a comparison of observations with simple formulas. Journal of Applied Meteorology, v. 32, p. 539-547, 1993 [ Links ]
DUARTE, H. F., DIAS, N. L., MAGGIOTTO, S. R. Assessing daytime downward longwave radiation estimates for clear and cloudy skies in Southern Brazil. Agricultural and Forest Meteorology, v. 139, p. 171-181, 2006 [ Links ]
FLERCHINGER, G. N., XAIO, W., MARKS, D., SAUER, T. J., YU, Q. Comparison of algorithms for incoming atmospheric long-wave radiation. Water Resources Research, v. 45, W03423, doi:10.1029/2008WR007394, 2009 [ Links ]
IDSO, S. B. A set of equations for full spectrum and 8 to 14 μm and 10.5 to 12.5 μm thermal radiation from cloudless skies. Water Resources Research, v. 17, p. 295-304, 1981 [ Links ]
IDSO, S. B., JACKSON, R. D. Thermal radiation from the atmosphere. Journal of Geophysical Research, v. 74, p. 5397-5403, 1969. [ Links ]
KJAERSGAARD, J. H., PLAUBORG, F. L., HANSEN, S. Comparison of models for calculating daytime long-wave irradiance using long term data set. Agricultural and Forest Meteorology, v. 143, p. 49 - 63, 2007 [ Links ]
KRUK, N. S., VENDRAME, I. F., ROCHA, H. R, CHOU, S. C., CABRAL, O. Downward longwave radiation estimates for clear and all-sky conditions in the Sertãozinho region of São Paulo, Brazil. Theoretical and Applied Climatology, v. 99, v. 115-123, 2010 [ Links ]
LIOU, K. N. An Introduction to Atmospheric Radiation. Academic Press, 392 pp., 1980 [ Links ]
OKE, T. R. Boundary Layer Climates, 2nd ed.Routledge, 435 pp., 1996 [ Links ]
PRATA, A. J. A new long-wave formula for estimating downward clear-sky radiation at the surface. Quarterly Journal of the Royal Meteorological Society, v. 122, p. 1127-1151, 1996 [ Links ]
SRIDHAR, V., ELLIOTT, R. L. On the development of a simple downwelling longwave radiation scheme. Agricultural and Forest Meteorology, v. 112, p. 237-243, 2002 [ Links ]
STULL , R. B. An Introduction to Boundary Layer Meteorology, Kluwer Academic Publishers, 666 pp.,1988 [ Links ]
SUGITA, M., BRUTSAERT, W. H. Cloud effect in the estimation of instantaneous downward longwave radiation. Water Resources Research, v. 29, p. 599-605, 1993 [ Links ]
SWINBANK, W. C. Long-wave radiation from clear skies. Quarterly Journal of the Royal Meteorological Society, v. 89, p. 339-348, 1963 [ Links ]
Received May 2010
Acepted December 2010