A Typical Meteorological Year Generation Based on NASA Satellite Imagery (GEOS-I) for Sokoto, Nigeria

1 Mechanical Engineering Department, Covenant University, Ota, Ogun State, Nigeria 2 Department of Ecology and Natural Resources Management, The Norwegian University of Life Sciences, Ås, Norway 3Department of Mechanical Engineering, University of Ibadan, Oyo State, Nigeria 4Department of Mechanical Engineering, Obafemi Awolowo University, Ile-Ife, Nigeria 5 Department of Mechanical Engineering, Michael Okpara University of Agriculture, Umudike, Abia State, Nigeria


Introduction
Energy remains the convergence point of most critical economic, environmental, and developmental issues confronting the world at the moment.Clean, efficient, stable, and sustainable energy services are ideal for global prosperity.Energy is paramount to achieving Nigeria's Vision 20 : 2020 needed by the country to be among the top 20 industrialized nations of the world.Lack of energy or its insufficiency in an economy is a potential source of social and economic poverty [1].In general, a larger proportion of energy is found to be consumed in buildings in Nigeria as is the case in many countries.There is thus a growing concern about energy consumption in buildings and the attendant adverse impacts on the environment and with the increasing population and socioeconomic growth (as witnessed in Nigeria), building unit will continue to be a key energy intensive sector [2].With the country's location on the tropics [3] (latitudes 4 ∘ 1  and 13 ∘ 9  North of the Equator and longitudes 2 ∘ 2  and 14 ∘ 30  East of the Greenwich Meridian), a large portion of the available energy is consumed on air conditioning.Energy consumption and electric load variation can be very sensitive to weather changes [1].In Koyuncuolu [4], the unique load variation along with the highly increasing energy consumed for air conditioning called for immediate energy saving programs in which substantial savings can be achieved by decreasing the cooling or heating loads of buildings.
According to Jiang [5], an accurate prediction of building thermal performance relies very much on the weather data such as dry bulb temperature, relative humidity, wind speed, and solar radiation data.Energy simulation in buildings thus offers a valuable tool for engineers and architects to evaluate building energy consumption before the building is constructed; since weather conditions vary significantly from year to year due to the impact of climate change, there is a need to derive a customized weather data set that can well 2 International Journal of Photoenergy represent the long-term averaged weather conditions over a year which may be referred to as a typical year [4,5].Forms of typical years include typical meteorological years (TMY), test reference year (TRY), and the weather year for energy calculations (WYEC).Typical weather year files like TMY are commonly used in building simulations to estimate the annual energy consumptions since analysis using a multiyear dataset is often not feasible and economical for the common design and analysis problems; representative days are too limited and often not accurate enough for a specific design and analysis problem [2].The typical year approach can thus reduce the computational efforts in simulation and weather data handling by using one year instead of multiple years.However, in Skeiker and Abdul Ghani [6], the term "typicality" could be interpreted in many ways as used by several investigators, namely, (i) the selection of weather which appears to be typical of an appropriate portion of the year, (ii) selection of a year which appears to be typical from several years of solar radiation data, and (iii) some investigators having run long periods of observational data in an attempt to simulate typical weather for the calculation.
Several works have been carried out on the development of TMY for various locations worldwide (e.g., [7][8][9][10][11][12][13][14]) whereas in Nigeria, existing work on the generation of weather data may be limited to (i) TRY generation for Ibadan, (ii) TMY generation for Port Harcourt zone, and (iii) generation of a typical meteorological year for northeast, Nigeria [1,15,16].The primary aim of this work is to generate TMY for Sokoto (located in North Nigeria with latitude-13.01∘ N, longitude-05.15∘ E, and altitude-350.8m) using a 23-year (1984-2006) hourly weather data from available satellite data acquired from NASA Geostationary Operational Environmental Satellite-1 (GEOS-1) Multiyear Time-Series Data. Figure 1 shows the studied site (Sokoto) and the other locations for which TMYs have been previously generated.Satellite and model-based products have been shown to be accurate enough to provide reliable solar and meteorological resource data over regions where surface measurements are sparse or nonexistent [17,18].This advantage of wide GEOS-1 data coverage has been utilized in the work of Dike et al., [19] for correlating meteorological estimates for seventeen locations within the southern region of Nigeria while Fadare [20] adopted the same dataset for modeling solar energy potential in Nigeria.
The selected weather data (including global solar radiation, dew point temperature, mean temperature, maximum temperature, minimum temperature, relative humidity, and wind speed) were chosen due to their impact on (and significant contribution to) the solar energy systems performance and heat gain and loss in buildings.For instance, in building applications, the knowledge of solar radiation is vital for accurate determination of cooling load in tropical regions (like Nigeria).It should also be noted that the air ambient temperature influences the thermal response of a building and the amount of heat gain and loss through it.The relative humidity is essential for the determination of latent heat for air-conditioning systems and evaporation levels.For solar energy conversion systems (Concentrated Solar Power and Photovoltaic systems), these meteorological data are also essential for design, selection, and performance evaluation of these systems.

Materials and Method
A TMY consists of 12 typical meteorological months (TMMs) selected from various calendar months in a multiyear weather database [5].Several methods including Festa and Ratto, Danish and Sandia commonly exist for the development of TMY [21][22][23].Sandia method is adopted in this work because of its wide applicability in the generation of weather data in tropical countries when compared with other methods [24].However, Sandia procedure is difficult to implement directly in spite of its wide deployment in generating TMY [25].In view of this, several procedures were described in literatures for generating TMY using Sandia method [13][14][15]25]; this work employed the procedural steps developed and outlined in Sawaqed et al., [25] due to its simplicity, for the development of TMY for Sokoto.
In Sawaqed et al., [25], TMY development begins with the generation of typical meteorological month (TMM).The TMM procedure involves the selection of five candidate years that are closest to the composite of all the years under study for each of the 12 calendar months.The statistical analysis commenced with the generation of daily weather indices from the adopter parameters for each day of record.Each of the seven sets of daily indices is sorted into bins by month and they are used to establish 12 long-term cumulative distribution functions (CDFs) given by [25] CDF  = 1  ,  = 1, 2, . . ., , where  is ranked order number and  is the total number of elements.A nonparametric method known as the Finkelstein-Schafer (FS) statistics [26] is then used to measure the closeness between the corresponding short-term and long-term CDFs as using (2) [26]: where   is the absolute difference between the short-term and the long-term CDFs for day  in the month, n is the number of days in the month, and  is the parameter (index) considered.For mean temperature j = 1, for max temperature, j = 2, and for solar flux j = 7.For each of the candidate months, the seven different FS statistics calculated for the seven indices are grouped into a composite weighted index using (3) [25] together with the weighting factors listed in Table 1.The resulting product gives the weighted sum (WS): where m is the number of indices (parameters/elements) considered,   is the weight for index j, FS  is the FS statistic for index j.
A weighting factor of 0.5 is assigned to global solar radiation while other weather indices were given equal   weights that add up to 0.5.This is because (i) in tropical region with warm weather (like Nigeria), solar heat gain can be more significant in cooling load calculations for building as well as for passive and active heating systems and (ii) other weather parameters are directly or indirectly affected by global solar radiation.
All individual months are ranked in ascending order of the WS values.A typical month is then selected by choosing from among the five months with the lowest WS values by following the procedure developed by Sawaqed et al. [25].The persistence of mean dry bulb temperature and daily global horizontal solar radiation are evaluated by determining the frequency and run length above and below fixed longterm percentiles.For mean daily dry bulb temperature, the frequency and run length above the 67th percentile and below the 33rd percentile are determined.For global horizontal radiation, the frequency and run length below the 33rd percentile are also determined.The persistence data are used to select, from the five candidate months, the month to be used in the TMY.The highest ranked candidate month in ascending order of the WS values that meet the persistence criterion is used in the TMY.
The generated TMY values were compared with the long-term average values for the location under study.The performance of the TMY was investigated using methods of stochastic analysis to calculate the root-mean-square error (RMSE) given in (4) as expressed in Ohunakin et al. [27]:

Results and Discussion
In this work, the mean temperature and global solar radiation are used in selecting 12 typical meteorological months (TMMs) and together with the remaining parameters, TMY for Sokoto is developed.For each month, the FS statistic is also estimated for every year and for all of the six parameters that have been considered.For example, the FS statistics of global solar radiation for Sokoto is presented in Table 2.It can be observed from this table that the FS statistics show variation from one month to the other.Similar variation is also observed among other variables.
With the daily values of the selected parameters determined from the measured data, the values of the WS for each month of the 23-year period were then calculated as shown in Table 3.The five candidate years for each months having the lowest values of WS together with the corresponding year are all arranged in ascending order.The persistence of mean dry bulb temperature and daily global horizontal radiation on the five candidate months were evaluated by determining the frequency and run length above and below fixed percentiles (33rd and 67th) as discussed in [8].The selected 12 TMMs are shown in Table 4.It can be seen that the 12 TMMs spread across the whole 23-year period; three TMMs are found in  1996 and two in 1999, whereas only one TMM exists in 1989, 1992, 1994, 1995, 2000, 2001, and 2002, respectively.The relationship between the long-term mean and the TMY for selected weather parameters are presented in Fig- ures 2-5.
Figure 2 showed that good relationships exist between the long-term mean (23 year) and the TMY values for the global solar radiation, though with very small variations.It can be observed that the TMY and the long-term solar fluxes exceed 24 MJ/m 2 /day from March to June with the peak as 25.1 and 25.8 MJ/m 2 /day in June and April for the TMY and long-term fluxes, respectively.Figure 3 shows the monthly average mean temperature based on TMY and long-term values.It can be seen that Figure 3    related to high weighting factor placed on each parameter being 0.5 for global solar radiation and 2/12 for mean temperatures (Table 1).The root-mean-square error (RMSE), which measures level accuracy between the TMY and longterm, is estimated as 0.62 ∘ C and 0.68 MJ/m 2 /day for the mean temperature and global solar radiation.These low values further confirm good agreement between the TMY and longterm values for these parameters.
The relationship between the relative humidity (Figure 4) and dew point temperature (Figure 5) is reasonable but not as close as that for the global solar radiation and the mean temperature.The RMSE values for the relative humidity and dew point temperature are 3.23% and 2.04 ∘ C, respectively.These values are relatively small, but higher than RMSE  values for the mean temperature and global solar radiation.This may be due to nonuniformity of the weighting factors and may also result from the fact that the selection of the TMMs is based on persistence of mean temperature and daily global horizontal radiation.Hence, the TMMs selected for the TMY can represent the long-term meteorological weather conditions for the Sokoto as considered in this study.

Conclusions
The TMY generated will be very useful for optimal design and evaluation of solar energy conversion systems, HVAC, and other solar energy dependent systems to be located within the vicinity of Sokoto.This work is a part of ongoing research work to develop a comprehensive TMY weather database for hourly building energy simulations in Nigeria and design of solar energy systems.In addition, the adopted satellite data have significantly advanced our ability to observe weather systems by providing frequent interval visible and infrared imagery of the earth surface, atmospheric moisture, and cloud cover that otherwise would have been impossible with ground observation.Findings based on the adoption of Filkenstein-Schafer statistical method utilized for the creation of a TMY for Sokoto show the following.
(i) The developed 12 TMMs are evenly spread across the whole 23-year period with three TMMs found in 1996 and two in 1999, and in 1989, 1992, 1994, 1995, 2000, 2001, and 2002 only one exists.
(ii) The close fit observed in the values suggests a close agreement between the TMY predictions and longterm averages.
(iii) Observed agreement between TMY and long-term values of the respective fluxes may be related to weighting factor assigned to each parameter.
New TMY site Existing TMY site Existing TRY site

Figure 1 :
Figure 1: Map of Nigeria showing the studied site and locations with existing TMY values.

Figure 2 :
Figure 2: Predicted TMY with long-term average values for the global solar radiation.
follows similar trend with Figure2.The peak values of the mean temperature are observed in April as 32.0 and 31.6 ∘ C for the TMY and long-term average values, respectively.The good agreement between TMY and longterm global solar radiation and mean temperature may be Jan Feb Mar Apr May Jun Jul Aug Sep Oct Nov Dec Month

Figure 3 :
Figure 3: Predicted TMY with long-term average values for the mean temperature.

Figure 4 :
Figure 4: Predicted TMY with long-term average values for the relative humidity.

International
Journal of PhotoenergyJan Feb Mar Apr May Jun Jul Aug Sep Oct Nov Dec Month

Figure 5 :
Figure 5: Predicted TMY with long-term average values for the dew point temperature.

Table 1 :
Weighting factors of FS statistics.

Table 2 :
FS statistics for global solar radiation.

Table 3 :
WS of FS statistics for the indices.