Empirical Models for the orrelation of Global Solar Radiation with Meteorological Data for Enugu, Nigeria.. Augustine and M.N. Nnabuchi Department of Industrial Physics, Ebonyi State University, Abakaliki, Nigeria. E-mail: emmyaustine3@yahoo.com ABSTRAT A number of multilinear regression equations were developed to predict the relationship between global solar radiation with one or more combinations of the following meteorological parameters: fraction of sunshine, maximum temperature, cloudiness index and relative humidity for Enugu for seventeen years (199 7). Using the Angstrom model as the base, other regression equations were developed by modifying the Angstrom equation. The values of the global solar radiation estimated by the models and the measured solar radiation were tested using the mean bias error (MBE), root mean square error (RMSE) and mean percentage error (MPE) statistical techniques. The values of the correlatiooefficient (R) and coefficient of correlation (R ) were also determined for each equation. The equation with the highest values of R, R and least values of MBE, RMSE and MPE is given as: R =.7 +.5.3 +.171 Τ +. The developed model can be used for estimating global solar radiation in Enugu and other locations with similar climatological characteristics. (Keywords: sunshine, maximum temperature, cloudiness index, relative humidity, solar radiation) INTRODUTION Knowledge of the global solar radiation is of fundamental importance for all solar energy conversion systems. Information on global solar radiation received at any site (preferably gained over a long period) should be useful not only to the locality where the radiation data is collected but also for the wider world community (Massaquoi,19). A global study of the world distribution of global solar radiation requires knowledge of the radiation data in various countries and for the purpose of world wide marketing, the designers and manufacturers of solar equipment will need to know the mean global solar radiation available in different and specific regions (Ibrahim, 195). Obviously, the best solar radiation information is that obtained from experimental measurements of the global and diffuse components of the solar insolation at the location in question. Unfortunately, there are few meteorological stations conducting such measurements. Therefore, it is rather important to develop models to estimate the solar radiation using whatever weather parameters available. Several empirical models have been developed to calculate the solar radiation using various parameters. Such models include that of Sambo(5), Awachie and Okeke (199), Akpabio et al. (), Badamus and Momoh (5), Hussaini et al. (5), Iheonu (1), Arinze and Obi (193), Falayi and Rabiu (5), El- sabaii et al. (5), Babatunde et al. (199), and Burari et al. (1) to mention but a few. This work tried to investigate the suite of empirical models used in the estimation of the global solar radiation from the meteorological parameters in Enugu and its environs. METHODOLOGY The monthly mean daily sunshine duration, maximum temperature, cloudliness index and relative humidity, were obtained from the Archives of Nigerian Meteorological Agency, Oshodi, Lagos. The monthly mean global solar radiations were collected from the Nigerian Meteorological Agency, Akanu Biam Airport, Enugu. The data obtained, covered a period of seventeen years The Pacific Journal of Science and Technology 93 Volume. Number 1. May 9 (Spring)
(199 7) for Enugu Located at latitude 7.55 N and longitude.5 E and it has an altitude of 1.5m above sea level. The monthly mean daily data processed in preparation for the correlations are presented in Table1. To develop the model, the global solar radiation measured using Gun-Bellani distillate were converted to useful form (MJ/m /day) using a conversion factor of 1.13 proposed by Sambo (195). The liner regression model used iorrelating the Η measured global solar radiation data Η with Ο fraction of sunshine n is given after Angstrom (19) and later modified by Prescott (19): Η n = a + b (1) Where a and b are regressioonstants, Η is the measured monthly mean daily global radiation, Η Ο is the monthly mean extraterrestrial solar radiation on horizontal surface, given by Iqbal (193) as follows: π Η Ο = Ι sceο Wssinφsinδ + cosφcosδ sin ws π 1 () Where Ι sc is the solar constant, E Ο is the eccentricity correction factor, φ is the latitude, δ is the solar declination and W s is the hour angle. The expressions for Ι sc, E Ο, δ and W s are given by Iqbal (193) : 137X 3 Ι sc = (MJm - day -1 ) (3) 3N Ε Ο = 1+.33OS 35 () ( N + ) 3 δ = 3.5Sin (5) 35 1 W = os tanφ tanδ () s ( ) Where N is the day number ranging from N = 1 on 1 st January to N = 35 on 31 st December. Table 1: Meteorological Data and Global Solar Radiation for Enugu. MONTH n Τ ( o ) c R/ Η (MJ/m /day) Η Ο (MJ/m /day) Κ Τ = JAN.5 35.1.1..5 35..397 FEB.5579 37.9..9 15.5 37.1.9 MAR.513 37.3.3.1.77 37.5.393 APR.59 35.57.3.71.7 3..91 MAY.57 3.7.39.77.5 3.1.31 JUN.9 3.5.. 13.1 33.15. JUL.3113 31... 11.5 3.5.333 AUG.3 31...3. 35.7.3 SEPT.37 31.9.1.. 3.95.331 OT.53 3.91.39.7 15.1 37.73.3 NOV.59 3.93.3.5 1.51 35.3. DE.5 35.9.3.5 15. 35..37 Η / Η Ο The Pacific Journal of Science and Technology 9 Volume. Number 1. May 9 (Spring)
The various meteorological data are related to global solar radiation. Multiple linear regression and correlation analysis of four parameters ( n, Τ, c and R ) were employed to estimate the global solar radiation. Where n is the fraction of sunshine duration, Τ is the c maximum temperature, is the loudiness index and R is the relative humidity (Table 1). SPSS omputer Software Program was used in evaluating model parameters. The values of the regressiooefficients obtained for Enugu are found to be different from values for Maiduguri, Bauchi, Ilorin, Nsukka, Onne, and Sokoto obtained by Hussain (5), Bala (1), Burari et al. (1), Awachie and Okeke (199), Akpabio et al. (), and Badamus et al. (5), respectively. These differences suggests that regressiooefficients associated with sunshine hours and other metrological data changes with latitude and atmospheric condition. The accuracy of the estimated values were tested by calculating the Mean Bias Error (MBE), Root Mean Square Error (RMSE), and Mean Percentage Error (MPE). The expressions for the MBE (MJm - day -1 ), RMSE (MJm - day -1 ), and MPE (%) is stated by El Sebaii et al. (5) as follows: ical imeas) n (7) ΡΕ = ( Η, Η, / 1 ical, imeas, / n RSΕ= ( Η Η ) () Ηi meas Η ΡΕ = Ηimeas,, i, cal Χ / n (9) Where Η ical, and Η imeas, is the ith calculated (predicted) and measured values and n is the total number of observations. Iqbal (193), Halouani (1993), Almorox (5) and he et al. (7) have recommended that a zero value for MBE is ideal and a low RMSE is desirable. The RMSE test provides information on the short- term performance of the studied model as it allows a term by term comparism of the actual deviation between the calculated values and the measured values. The MPE test gives long term performance of the examined regression equations, a positive MPE values provide the averages amount of overestimation in the calculated values, while the negative values gives underestimation. A low value of MPE is desirable by Akpabio et al. (). RESULTS AND DISUSSION The values for the fraction of sunshine duration, maximum temperature, cloudiness index, relative humidity, and clearness index are presented in Table 1. Table contains summaries of various linear regression analysis, obtained from the application of Equation (1) to the monthly mean values for the four variables under study. Figure 1 shows the comparison between the measured and predicted values of correlation equations. It is clear that the correlation coefficient R, coefficient of determination R, MBE (MJ/m /day), RMSE (MJ/m /day) and MPE (%) vary from one variable to another variable. ONE VARIABLE ORRELATION The correlatiooefficient of.57 exists between the clearness index and monthly mean daily fraction of sunshine and coefficient of determination of.735 implies that 73.5% of clearness index can be accounted using fraction of sunshine..191.33 n = + () The correlatiooefficient of. exists between the clearness index and monthly mean daily maximum temperature and coefficient of determination of. implies that % of clearness index can be accounted using maximum temperature. =.15 +.1 Τ (11) The Pacific Journal of Science and Technology 95 Volume. Number 1. May 9 (Spring)
Table : Regression Equations and Statistical Indicators. Eqauations R R MBE RMSE MPE.191.33 n.57.735.3.9 -. = +... 1. 1.5 =.15 +.1 Τ.93. c.3..3 1.3-3. = R.53.5.39 1.39-3.9 =.5.1 n.57.735.33. -1.75 =. +.. Τ.9.1.35.79 -.95 =. +.571 +.75 n R.91.79.117.3-1.577 =.37 +.51 +..99..3.79 -.539 =.953 +.5 +.3 +.39 Τ n R.915.37 -.3.9339 -.19 =.39 +.9 +.53 +.11 Τ n R c.9..7.9179 -.79 =.39 +.51.5 +.375 R.91...111. =.7 +.5.3 +.171 Τ +. The correlatiooefficient of.3 exists between the clearness index and monthly mean daily cloudiness index and coefficient of determination of. implies that.% of clearness index can be accounted using cloudiness index. clearness index can be accounted using relative humidity. R =.5.1 (13).93. c = () TWO VARIABLES ORRELATIONS The correlatiooefficient of.53 exists between the clearness index and monthly mean daily relative humidity and coefficient of determination of.5 implies that.5% of The correlatiooefficient of.57 exists fraction of sunshine duration and maximum temperature, also coefficient of determination of.735 implies that 73.5% of clearness index can The Pacific Journal of Science and Technology 9 Volume. Number 1. May 9 (Spring)
be accounted using fraction of sunshine and maximum temperature. n =. +.. Τ () The correlatiooefficient of.9 exists fraction of sunshine duration and cloudiness index, also coefficient of determination of.1 implies that 1.% of clearness index can be accounted using fraction of sunshine and cloudliness index. =. +.571 +.75 1 1 1 1 Equation 1 3 5 7 9 11 Equation 11 1 3 5 7 9 11 (15) Figure 1 (Equation -11): omparison Between the Measured and Predicted Values of orrelation Equation. Equation 1 1 1 3 5 7 9 11 Equation 13 1 1 1 3 5 7 9 11 1 1 Equation 1 3 5 7 9 11 Figure 1 (Equation ): omparison Between the Measured and Predicted Values of orrelation Equation. The correlatiooefficient of.91 exists fraction of sunshine duration and relative humidity, also coefficient of determination of.79 implies that 79.% of clearness index can be accounted using fraction of sunshine and relative humidity. The Pacific Journal of Science and Technology 97 Volume. Number 1. May 9 (Spring)
1 1 Equation 15 1 3 5 7 9 11 1 1 Equation 1 1 3 5 7 9 11 1 1 Equation 17 1 3 5 7 9 11 1 1 Equation 1 1 3 5 7 9 11 1 1 Equation 19 1 3 5 7 9 11 1 1 Equation 1 3 5 7 9 11 Figure 1 (Equation 15 - ): omparison Between the Measured and Predicted Values of orrelation Equation. n R =.37 +.51 +. (1) THREE VARIABLES ORRELATIONS The Equations () (1) can be modified by incorporating extra parameters to the set of correlation equations for two variables. The correlatiooefficient of.99 exists The Pacific Journal of Science and Technology 9 Volume. Number 1. May 9 (Spring)
fraction of sunshine duration, cloudiness index and maximum temperature, also coefficient of determination of. implies that.% of clearness index can be accounted using fraction of sunshine, cloudiness index and maximum temperature. =.953 +.5 +.3 +.39 Τ (17) The correlatiooefficient of.915 exists fraction of sunshine duration, relative humidity and maximum temperature, also coefficient of determination of.37 implies that 3.7% of clearness index can be accounted using fraction of sunshine, relative humidity and maximum temperature. n R =.39 +.9 +.53 +.11 Τ (1) The correlatiooefficient of.9 exists fraction of sunshine duration, relative humidity and cloudiness index, also coefficient of determination of. implies that % relative humidity and cloudliness index. n R c =.39 +.51.5 +.375 (19) FOUR VARIABLES ORRELATION Equations (17) (19) can be modified by incorporating extra parameters to the set of correlation equations for three variables. The correlatiooefficient of.91 exists fraction of sunshine duration, cloudiness index, maximum temperature and relative humidity, also coefficient of determination of. implies % of clearness index can be accounted using fraction of sunshine, cloudiness index, maximum temperature and relative humidity. R =.7 +.5.3 +.171 Τ +. () For better analysis the developed correlation will be considered that has high values of correlation coefficient R and coefficient of correlation R : Equations (), (11), (), (13), (), (15), (1), (17), (1), (19), and (). Equations (), (), (15), (1), (17), (1), (19), and () have the highest values of correlatiooefficient and coefficient of determination, while the others have low values of correlatiooefficient and coefficient of determination. Furthermore, there is a remarkable agreement between the observed and the predicted values of global solar radiation for seventeen years from our correlation (Table ). From Table, based on the highest values of correlatiooefficient, coefficient of determination and least values of MBE, RMSE and MPE, Equation () is the best regression equation suitable for estimating solar radiation in Enugu and its environs. ONLUSION The monthly mean daily global solar radiation, fraction of sunshine duration, maximum temperature, cloudiness index and relative humidity have been employed in this study to develop several correlation equations. Four variables have been developed with different types of equations obtained. It was observed that Equation () has the highest value of correlation coefficient and coefficient of determination, which gives good results wheonsidering statistical indicators that is, MBE, RMSE, and MPE. The equatioould be employed in the estimation of global solar radiation in Enugu and other similar locations. AKNOWLEDGEMENTS The authors are grateful to the management and of Nigerian Meteorological Agency, Oshodi, Lagos State for making the Meteorological data available and the management and staff of Nigerian Meteorological Agency, Akanu Biam Airport, Enugu for making the data of global solar radiation available. The Pacific Journal of Science and Technology 99 Volume. Number 1. May 9 (Spring)
REFERENES 1. Awachie, I.R.N. and Okeke,.E. 199. New Empirical Solar Model and its use in Predicting Solar Irradiation. Nigerian Journey of Solar Energy. 9: 155.. Angstrom, A.S. 19. Solar and Terrestrial Radiation. Meteorological Society. 5:1-. 3. Akpabio, L.E. and Etuk, S.E.. Relationship Between Solar Radiation and Sunshine Duration for Onne. Nigeria.,Turkish J. Physics. 7: 11 17.. Arinze, E.A. and Obi, S.E. 193. Solar Energy Availability and Prediction in Northern Nigeria. Nig. J. Solar Energy. 3: 3. 5. Almorox, J., Benito, M., and Hontoria,. 5. Estimating of Monthly Angstrom Prescott Equation oefficients from Measured Daily Data in Toledo, Spain. Renewable Energy Journey.3: 931 93. Babatunde, E.B. and Aro, T.O. 199. haracteristic Variation of Total (Global) Solar Radiation at IIorin, Nigeria. Nigerian Journal of Solar Energy. 9:157 173 7. Badmus, I.B. and Momoh, M. 5. omparison of Models for Estimating Monthly Average Daily Insolation on a Horizontal Surface. 1st Science Association of Nigeria onference, Sokoto, Nigeria. 5-9 April, 5.. Burari, F.W. and Sambo, A.S. 1. Model for the Prediction of Global Solar for Bauchi using Meteorological Data. Nig. J. Renewable Energy, 91: 3-33. 9. El-Sebaii, A.A. and Trabea, A.A. 5. Estimation of Global Solar Radiation on Horizontal Surfaces Over Egypt. Egypt. J. Solids. : 1.. Halouani M, Nguyen. T., and Vo Ngoc, D. 1993. alculation of Monthly Average Global Solar Radiation on Horizontal Surfaces using Daily Hours of Bright Sunshine. Solar Energy. 5: 7 55. 13. Igbal, M. 193. An Introduction to Solar Radiation. Academy Press: New York, NY.. Iheonu, E.E. 1. Model for the Prediction of Average Monthly Global Solar Radiation on a Horizontal Surface for some Locations in the Tropics using Temperature Data. Nigerian Journey of Solar Energy. 9:-15. 15. Ibrahim, S.M.A. 195. Predicted and Measured Global Solar Radiation in Egypt. Solar Energy. 35: 15 1. 1. Massaquoi, J.G.M. 19. Global Solar Radiation in Sierra Leone (West Africa). Solar Wind Technol. 5: 1 3 17. Prescott, J.A. 19. Evaporation from a Water Surface in Relation to Solar Radiation. Tran. R. Soc. S. Austr. :1 11. 1. Sambo. A.S. 195. Solar Radiation in Kano. A orrelation with Meteorological Data. Nig. J. Solar Energy. : 59-. SUGGESTED ITATION Augustine,. and M.N. Nnabuchi. 9. Empirical Models for the orrelation of Global Solar Radiation with Meteorological Data for Enugu, Nigeria. Pacific Journal of Science and Technology. (1):93-7. Pacific Journal of Science and Technology. Falayi, E.O. and Rabiu, A.B. 5. Modelling Global Solar Radiation Using Sunshine Duration Data. Nig. J. Physics. 17:11-1. 11. Hussaini, A.M., Maina, M. and Onyewuenyi, E.. 5. orrelation of Solar Radiation with some Meteorological Parameters for Maiduguri, Borno State, Nigeria. Nigerian Journey of Solar Energy. 15: 19. The Pacific Journal of Science and Technology 7 Volume. Number 1. May 9 (Spring)