International Journal of Photoenergy

Volume 2014 (2014), Article ID 680596, 9 pages

http://dx.doi.org/10.1155/2014/680596

Research Article

Estimation of Monthly Sunshine Duration in Turkey Using Artificial Neural Networks

H. M. Kandirmaz,¹ K. Kaba,¹ and M. Avci²

¹Department of Physics, Cukurova University, 01330 Adana, Turkey

²Biomedical Engineering Department, Cukurova University, 01330 Adana, Turkey

Received 22 May 2014; Accepted 18 July 2014; Published 10 August 2014

Academic Editor: Shahed Khan

Copyright © 2014 H. M. Kandirmaz et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

This paper introduces an artificial neural network (ANN) approach for estimating monthly mean daily values of global sunshine duration (SD) for Turkey. Three different ANN models, namely, GRNN, MLP, and RBF, were used in the estimation processes. A climatic variable (cloud cover) and two geographical variables (day length and month) were used as input parameters in order to obtain monthly mean SD as output. The datasets of 34 stations which spread across Turkey were split into two parts. First part covering 21 years (1980–2000) was used for training and second part covering last six years (2001–2006) was used for testing. Statistical indicators have shown that, GRNN and MLP models produced better results than the RBF model and can be used safely for the estimation of monthly mean SD.

1. Introduction

The importance of solar energy is increasing because demands for energy are increasing from day to day and the world’s energy resources are limited. Due to this fact, a large number of studies are carried out and will continue to be so in the future in order to benefit from the solar energy which is priceless and clean. On the other side, the significant variations of solar radiation over long years cause important climatic changes which affect the life in the world from many aspects. For example, a decrease in the amount of solar radiation between 1950s and early 1980s was observed, which is known as global dimming, and this phenomenon was reversed from second half of 1980s known as global brightening [1–6]. It is then crucial to know and analyze the variations of solar radiation in time and space. Global solar radiation measurements are generally made with the actinographs which are often not reliable due to the need of routine calibration of thermal sensitivity of the mechanical components of their sensors. More accurate measurements can be done by constructing networks with calibrated modern pyranometers but this is not the case for many countries because this type of instruments is expensive. It has been proved that solar radiation is highly correlated to the sunshine duration (SD) [7–10]. This obviously means that if one has information about the SD over an area then one also has information about the solar radiation over that area or vice versa. Hence, long-term accurate observations of global SD become very important and necessary for climatologic and some other applications [1, 6, 7, 11–13]. In fact, SD measurements have been achieved at many locations over the world more accurately than the solar radiation with cheaper instruments for long time periods. For instance, it has been measured in Ankara station in Turkey since 1928 [14]. SD measurements have been done in The Netherlands since 1901 [9] and at eight meteorological stations in Taiwan since 1910 or earlier [15]. Records of sunshine hours have been available from four stations of Ireland since the late 19th century [16]. Measurements of duration of sunshine in the United States began in 1891 at 20 U.S. Weather Bureau (USWB) stations [17]. Although SD can be measured almost in all meteorological stations, the networks of meteorological stations are still insufficient and/or limited for some parts of the world due to geographical and sometimes financial problems. This means that some countries have limited or sometimes unreliable SD database and maps. In order to overcome this problem, estimation techniques are improved for the areas where SD measurements are not available or reliable, especially for remote and inaccessible regions. Up to now, many researchers have concentrated themselves on the estimation of global solar radiation so that numerous studies can be found in the literature. In contrast to this, limited numbers of studies have been done for the estimation of SD, its relation with geographic and climatologic variables, and its variation throughout time and space.

Reddy developed an empirical model to calculate sunshine from total cloud amount [18]. The relationship between point cloudiness and sunshine derived cloud cover was investigated using data recorded at 34 stations in India by Raju and Kumar [19]. A cubic regression equation between monthly mean values of fraction of the sky and duration of bright sunshine was proposed by Rangarajan et al. and using this relationship they computed SD from cloud cover data [20]. Harrison and Coombes derived a second order regression equation which defines the statistical relationship between long-term averages of monthly cloud shade and point cloudiness using data of 43 Canadian weather stations [21]. Yin used monthly data of 729 worldwide stations for finding a generic algorithm that captures global variation of bright SD data in relation to temperature, precipitation, and geographic location [22]. El-Metwally proposed a nonlinear relationship based on maximum and minimum air temperatures and cloud cover fraction for estimating relative SD [23]. Matzarakis and Katsoulis estimated mean annual and seasonal duration of bright sunshine from an empirical formula which depends on distance of each station from the nearest cost, height of above sea level for each station location, percentage of land cover around each station, latitude of each station, and longitude of each station [24]. Kandirmaz proposed a simple model for the estimation of daily global SD and constructed spatially continuous map of SD over Turkey using meteorological geostationary satellite data [25]. Robaa established a simple model using observed cloud data and estimated SD at any region in Egypt [26]. Kandirmaz and Kaba showed that the statistical relation between daily satellite-derived cloud cover index and measured SD is quadratic rather than linear [27].

Recently, artificial neural networks (ANNs) have been extensively used for the solution of many problems in many areas, such as optimization, prediction, and modeling in engineering, climate science, pattern classification, and economy. ANNs are effective tools for modeling nonlinear systems [28–31]. Many previous studies have shown that usage of ANN techniques is an alternative and strong key for prediction of global solar radiation as compared to classical regression models. In ANN based solar radiation studies, meteorological, geographical, and climatological variables such as SD, temperature, cloud cover, relative humidity, wind speed, vapor pressure, precipitation, elevation, latitude, longitude, month, and satellite recorded or derived variables were used as input variables for obtaining the solar radiation as an output [32–45]. Although many ANN based studies have been in existence for the estimation of solar radiation in the literature, there have been only two studies for the estimation of SD. Jervase et al. generated contour maps for sunshine hours and sunshine ratios for Oman using a radial basis function (RBF) neural network model [46]. In their study latitude, longitude, altitude, and month of the year were used as input parameters. Mohandes and Rehman estimated the SD over Saudi Arabia using two neural network algorithms in which maximum possible day length, extraterrestrial solar radiation at that particular location, latitude, longitude, altitude, and month number were used as input parameters [47]. On the other side, only a few studies have been devoted to the estimation, determination, and distribution of SD for Turkey. Aksoy analyzed the changes in SD for the Ankara station in Turkey [48]. Kandirmaz [25] and Kandirmaz and Kaba [27] used satellite data for predicting the SD in Turkey. Şahin developed a simple model and estimated SD for some stations in Turkey [49]. Trends of the measured SD of 36 stations were analyzed by Yildirim et al. [50]. The present work is the first study where an ANN approach is proposed for the estimation of SD for Turkey. In order to get the best possible results, we have constructed three different ANN models in which cloudcover, day length, and month have been used as input parameters.

2. Data Sources and Methodology

This research was conducted for Turkey which is located between latitudes 36° and 42°N and longitudes 26° and 45°E and has seven different climatic zones, namely, Black Sea, Marmara, Eastern Anatolian, Southeastern Anatolian, Mediterranean, Aegean, and Central Anatolian [51]. Its average annual solar radiation was determined as 1311 kWh/m²-year (3.6 kWh/m²-day) and the annual average total SD as 2640 hours (7.2 hours/day) for the period from 1966 to 1982 [52]. Ground-measured daily SD and cloud fraction data are collected from 34 selected stations of the Turkish State Meteorological Service (TSMS). These stations cover almost the whole country and hereby reflect all climatic properties of the country. The geographical positions of these stations can be seen in Figure 1.

Figure 1: The geographical distributions of 34 meteorological stations in Turkey.

Previous studies have shown that SD can be interrelated with the cloud cover, air temperature, precipitation, relative humidity, wind speed, and geographical variables or combination of some of variables given above [53–58]. Clouds, consisting of liquid water droplets or ice particles, decrease incoming solar radiation in many ways before reaching earth surface. Physical explanation of the interactions between clouds and solar rays is rather difficult because these interactions depend on the size and shape of droplet or particles and total mass of water and spatial distribution. Other meteorological variables may affect the incoming solar radiation but the main and the greatest effect comes from the clouds [10, 27].

Day length, which is a function of latitude and solar declination of the site interested, gives the maximum possible duration of sunshine in a day. Since the cloud cover and day length are the most responsible elements for determining the daily bright sunshine hours, they were used as input data in our ANN models. Daily cloud coverage and bright sunshine hours used in this study were made available by Turkish State Meteorological Service (TSMS). Cloud coverage observations have been made three times in a day by trained meteorologists and their average value is assigned to a day. Bright sunshine hours are recorded by Campbell-Stokes type sunshine recorder. In this instrument the solar energy is concentrated by a special lens and it detects the sunshine if the energy of the beam is high enough to burn a special dark paper card. A mean value of 120 Wm⁻² of direct solar irradiance is accepted as the threshold value by World Meteorological Organization (WMO). SD should be measured with an uncertainty of 0.1 h and a resolution of 0.1 h. Adjustments and calibrations of the sunshine recorders are routinely made by TSMS. Day length is calculated for the station of interest by using the following formula [59]:

where is the latitude of location in the range

and is the solar declination

and is the number of day of the year starting from first of January.

Validation of the simulated values of the models is achieved using the five statistical indices: root mean square error (RMSE), percentage root mean square error (%RMSE), mean bias error (MBE), mean absolute error (MAE), and percentage mean absolute error (%MAE) defined by the following mathematical equations:

where is the estimated value and is the measured value and is the number of data.

MBE is used to test whether the proposed model tends to overestimate or underestimate the measured value and generally it provides good information for long-term observations. On the other side, RMSE generally provides valuable information for short-term applications and it explains the measure of differences between measured and estimated values. As the values of RMSE, %RMSE, MBE, MAE, and %MAE get closer to zero, performances of the models get better.

2.1. Artificial Neural Network

An ANN consists of biological neuron like operation units (nodes) linked together according to a specific architecture. ANNs have generally one input layer, one output layer, and some hidden layers. Hidden and output layers contain activation functions. Neurons between adjacent layers are interconnected. Connection weights are multiplied by inputs to obtain product terms. Sum of products and biases then applied to a transfer function through the output layer. The result of the output layer contains total effect of all the neurons in the network [34, 60].

2.1.1. Generalized Regression Neural Network

Nadaraya-Watson kernel based general regression neural network (GRNN) is introduced by Specht in 1991 [61]. It is used for many applications ever since. Some of these are function approximation, prediction, and control, medical diagnosis, engineering, speech recognition, and 3D modeling [62–69]. In many studies GRNN performed better function approximation feedforward networks and other statistical neural networks on some datasets [62, 64–66, 68–70]. Although GRNN is proposed for function approximation purposes [62], in some works, it is applied to classification problems with small modifications [65, 66]. Main advantages of GRNN are fast learning, consistency, and optimal regression with large number of samples [67]. GRNN has four layers: input, pattern, summation, and output as shown in Figure 2.

Figure 2: Structure of the GRNN.

Input layer holds corresponding input without any operation and supplies input vector to pattern layer. Pattern layer consists of neurons for each training datum. In this layer, weighted squared Euclidean distance is calculated as shown in (5). Any test input applied to network is first subtracted from pattern layer neuron values and then according to the distance function either squares or absolute values of subtracts are summed and applied to activation function. Exponential function is the most popular activation function. Results are transferred to summation layer. Summation layer neurons add dot product of pattern layer outputs and weights. In Figure 2 weights are shown by and , their values are determined by values of training data stored at pattern layer, and denotes weighted outputs of pattern layer where is a Parzen window associated constant. denotes multiplication of pattern layer outputs and training data output values. At output layer, is divided by to estimate desired , given in (6), (7) [61, 71, 72]:

GRNN is also known as normalized RBF-NN. RBF units in GRNN are probability density functions shown in (7). In GRNN structure only smoothing parameter (), also known as bandwidth, is updated during training phase [66, 73]. values are important for network performance; smaller limits the number of effective samples; on the other hand the larger one extends radius of effective neighbors [62–68].

2.1.2. Multilayer Perceptron (MLP)

Unlike single perceptron, multilayer perceptron (MLP) contains input, output, and one or more hidden layers with computation nodes [73]. Figure 3 is an example of architecture of MLP.

Figure 3: Structure of the MLP.

MLP provides a solution to overcome the limitations of the perceptron. Input layer holds input to first hidden layer without any operation. Main computational tasks are done in hidden layers. Last hidden layer transmits results to output layer. Model of each neuron includes a differentiable nonlinear activation function. Connectivity is provided by synaptic weights of the network. Default training method of MLP is error backpropagation algorithm containing a special kind of LMS algorithm. Gradient descent is one of the most popular ones. Training includes forward and backward steps. In the forward step, synaptic weights of the network are fixed. Operation results of inputs with these fixed weights and activation functions determine the changes. In this step, each hidden and output neuron compute continuous nonlinear function of the input and weights associated with that neuron. Typically, activation functions used in hidden neurons are tangent hyperbolic sigmoid and logarithmic sigmoid functions given in (8) and (9), respectively. Forward step calculations are summarized as shown in (10). Each unit in layer receives from the previous layer of processing units and sends activations to the next layer of units:

In the backward step, error is calculated by subtracting network output with desired output. Then, square error is tried to be minimized while propagating through the network in the backward direction. During propagation, weights are updated to minimize the error. In backward step, each hidden and output neuron computes the gradient vector for learning. Update operation of weight is given in (11) where denotes the output layer neuron for data, denotes learning rate, is error function, denotes target value of data for output layer neuron, and denotes the weight value between layer neuron and layer neuron:

In our study, optimum results were obtained with three inputs, by 13 neurons in one hidden layer and a single neuron in the output layer for the MLP. Tansig and purelin were used in hidden layer and output layer as activation function, respectively.

2.1.3. Radial Basis Function Network

Radial Basis Function Neural Network (RBF-NN) is three layered feed-forward network type applicable to various regression and classification problems. Structure of the RBF is given in Figure 4.

Figure 4: Structure of the RBF.

Input layer has neurons for dimensional inputs. Hidden layer contains the same number of neurons as the size of the training sample, . In this layer, activation function used by neurons is Radial Basis Function given in (12), where data point denotes the center point of RBF, is input, and is width of Gaussian function:

Output layer consists of neurons to calculate sum of weighted hidden layer outputs, given in (13). Although there is no limitation on the size of output layer, typically the size of output layer is smaller than hidden layer:

3. Results and Discussion

In this study three different ANN models, namely, GRNN, MLP, and RBF were employed in order to estimate the monthly mean daily SD for 34 stations in Turkey. The inputs of the networks were monthly mean values of cloud coverage and day length, and the output was monthly mean daily SD. Data belonging to selected 34 stations were subdivided into two separate datasets. First part covering first 21 years (1980–2000) was used in the training process and the second part covering last six years (2001–2006) was used in the testing process. The models were constructed and optimized by varying the number of neurons in hidden layer. Matlab ANN toolbox was used for each modeling.

Monthly mean daily SD values were estimated for each station by the proposed ANNs and their averages were compared with the values recorded at meteorological stations from 2001 to 2006. Figure 5(a) shows the comparison of observed values with the estimated values for each station for the considered period. It has been generally observed that simulated values were very close to those observed at meteorological stations. In order to further clarify these results, MAE and RMSE were calculated for each model and station and the graphical representations were given in Figures 5(b) and 5(c), respectively.

(a)

(b)

(c)

Figure 5: (a) Measured and estimated values of average monthly mean SD for each station from 2001 to 2006. (b) MAE for each station. (c) RMSE for each station.

It can be concluded from Figure 5(b) that the GRNN model produced much closer results to the observed values at 14 stations (Rize, Trabzon, Bursa, Erzincan, Erzurum, Sivas, Ankara, Kirikkale, Denizli, Adana, Elazig, Mersin, Izmir, and Antalya). MLP model produced closer results at 11 stations (Kocaeli, Samsun, Kastamonu, Zonguldak, Tokat, Afyon, Istanbul, Yozgat, Iskenderun, Adiyaman, and Van) than at the others. On the other hand, RBF model predicted better results for the remaining nine stations (Bolu, Sakarya, Kutahya, Kahramanmaras, Gaziantep, Kayseri, Konya, Diyarbakir, and Sanliurfa). Lower MAE values (less than 0.2 h) were found for Bolu, Zonguldak, Tokat, Ankara, Konya Izmir, Van, and Antalya stations while higher MAE values (greater than 0.4 h) were found for Bursa, Kahramanmaras, and Sanliurfa stations by all three models. On the other side MLP model produced lower RMSE results for 17 stations (Bolu, Kocaeli, Samsun, Kutahya, Kastamonu, Afyon, Istanbul, Ankara, Iskenderun, Kirikkale, Denizli, Yozgat, Adana, Diyarbakir, Adiyaman, Van, and Antalya) than the GRNN (Rize, Trabzon, Sakarya, Tokat, Bursa, Erzincan, Erzurum, Kahramanmaras, Gaziantep, Sivas, Konya, Sanliurfa, Elazig, Mersin, and Izmir,) and RBF (Zonguldak and Kayseri) models. The RMSE values were found to be less than 1.0 h for all stations for GRNN and MLP models. RBF model produced high RMSE values (greater than 1.0 h) for stations Bolu, Sakarya, Samsun, Kastamonu, Erzincan, Istanbul, Gaziantep, Iskenderun, Kirikkale, and Yozgat. The average of observed monthly mean daily sums of SD was estimated as 6.7723 h, 6.7952 h and 0.6715 h by GRNN, MLP and RBF models with MBE 0.0488 h, 0.0719 h and −0,0544 h, MAPE 9.88%, 9.36% and 10.75%, RMSE 0.7971, 0.6256 and 1.1279, %RMSE 11.75%, 9.30% and 16.78% and linear regression of observed versus predicted values () was found to be 0.9530, 0.9563, and 0.8902, respectively.

The range of monthly mean daily SD values for the all stations is around 4.0 h (minimum value is around 4.4 h and maximum value is around 8.4 h) which really indicates that the country has geographically and climatologically different zones. Generally, lower SD values have been observed for the stations especially located inside the Black Sea Region (Rize, Trabzon, Bolu, Samsun, Kastamonu, and Zonguldak) and Marmara Region (Sakarya, Istanbul, Kocaeli, and Bursa), which is highly affected by Black Sea climate, north part of East Anatolia (Erzincan and Erzurum), and Central Anatolia (Kutahya and Tokat). This was expected because Black Sea Coasts have much more cloudy days than the other regions and receive the greatest amount of rainfall. On the other hand, higher SD values observed for the stations that stay inside Southeastern Anatolia (Sanliurfa, Diyarbakir, Adiyaman, and Gaziantep), Mediterranean (Antalya, Mersin, Adana, Iskenderun, and Kahramanmaras), and Aegean (Izmir, Denizli) and Central Anatolia (Ankara, Kirikkale, Konya, Kayseri, and Yozgat) and south part of East Anatolia (Van having latitude as 38,47) regions. It is certain that the cloudiness increases as one goes from lower latitude (south of country) to the higher latitude (north of country).

Results of the current study are comparable with previous studies in which neural network approaches and other methodologies were used for other geographical locations. Jervase et al. used RBF neural network model and found that the estimated values deviate from the measured values in between 0 and 1 hour for the stations inside Oman [46]. Mohandes and Rehman obtained the minimum mean absolute percent errors (MAPE) for Al-Madina station as 2.3% and 2.7% and maximum values for Al-Numas station in Saudi Arabia as 22.9% and 16.7% using PSO and SVM methods [47]. The duration of sunshine varies between 7.4 h and 9.4 h per day and its average daily value is approximately 8.89 h for Saudi Arabia [47] and 9.5 h for Oman [74]. Due to their locations both of these countries have less cloudy days and have much solar energy and SD, in comparison to Turkey. Probably this was the reason why Mohandes and Rehman [47] and Jervase et al. [46] did not use the cloudiness as one of the input parameter in their studies. That is, the effect of cloudiness on SD is more impressive on Turkey and it has been used as an input parameter in the present study.

El-Metwally estimated relative SD for six sites in Egypt. MBE% and RMSE% values varied from −0.2% to −13.3% and 2.3% to 14.5%, respectively [23]. Temporal and spatial distribution of bright sunshine hours over Greece were estimated by Matzarakis and Katsoulis and the correlation coefficient () and RMSE were found to be 0.87 and 9.90 h, 0.58 and 6.15 h, 0.89 and 4.69 h, 0.86 and 6.22 h, and 0.84 and 5.33 h for annual sunshine, winter, spring, summer, and autumn, respectively [24]. The relative SD was estimated from the cloud data using three empirical formulae for Egypt in the study of Robaa. It was found that relative percentage error (), mean percentage error (MPE), MBE, and RMSE were changed from −7.2698% to +3.7908%, −0.6240% to +0.8069%, −0.0053 to +0.0070, and 0.0046 to 0.0160, respectively [26]. A simple model was set up by Stanghellini [75] to evaluate monthly SD for various sites in Italy on the basis of the mean daily cloudiness. It is found that monthly MBE was in the range of −17.3 h to 14.9 h.

4. Conclusions

This paper presents a study on the monthly mean estimation of daily SD using three ANN methods, GRNN, MLP, and RBF, which were applied to 34 stations in Turkey. Month, day length, and cloud coverage data were selected and used as input parameters of the constructed models. Since day length can be calculated using astronomical factors and cloud cover obtained visually, monthly mean SD can be determined very accurately for any region by choosing an appropriate ANN model and without using any measuring instrument if enough historical SD and cloud cover database exist. The statistical indicators have shown that GRNN and MLP models work better than RBF model. Results obtained here seem to be good enough because Turkey has geographically and climatologically diverse zones, meaning that the range of distribution of SD over the country is not homogeneous. Also, it should be noted that the visual determination of cloud coverage is a subjective work which may also cause some error and affects the accuracy of the models used here.

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

References

G. M. Abakumova, E. M. Feigelson, V. Russak, and V. V. Stadnik, “Evaluation of long-term changes in radiation, cloudiness, and surface temperature on the territory of the former Soviet Union,” Journal of Climate, vol. 9, no. 6, pp. 1319–1327, 1996.

G. Stanhill and S. Cohen, “Global dimming: a review of the evidence for a widespread and significant reduction in global radiation with discussion of its probable causes and possible agricultural consequences,” Agricultural and Forest Meteorology, vol. 107, no. 4, pp. 255–278, 2001.

B. Liepert and I. Tegen, “Multidecadal solar radiation trends in the United States and Germany and direct tropospheric aerosol forcing,” Journal of Geophysical Research D, vol. 107, no. 12, pp. 7–15, 2002.

M. Wild, H. Gilgen, A. Roesch et al., “From dimming to brightening: decadal changes in solar radiation at earth's surface,” Science, vol. 308, no. 5723, pp. 847–850, 2005.

R. T. Pinker, B. Zhang, and E. G. Dutton, “Do satellites detect trends in surface solar radiation?” Science, vol. 308, no. 5723, pp. 850–854, 2005.

A. Sanchez-Lorenzo, J. Calbó, M. Brunetti, and C. Deser, “Dimming/brightening over the Iberian Peninsula: trends in sunshine duration and cloud cover and their relations with atmospheric circulation,” Journal of Geophysical Research, vol. 114, no. 8, Article ID D00D09, 2009.

A. Angström, “Solar and terrestrial radiation,” Quarterly Journal of the Royal Meteorological Society, vol. 50, pp. 121–125, 1924.

H. Suehrcke, R. S. Bowden, and K. G. T. Hollands, “Relationship between sunshine duration and solar radiation,” Solar Energy, vol. 92, pp. 160–171, 2013.

Y. B. Hinssen and W. H. Knap, “Comparison of pyranometric and pyrheliometric methods for the determination of sunshine duration,” Journal of Atmospheric and Oceanic Technology, vol. 24, no. 5, pp. 835–846, 2007.

H. Suehrcke, “On the relationship between duration of sunshine and solar radiation on the earth's surface: angström's equation revisited,” Solar Energy, vol. 68, no. 5, pp. 417–425, 2000.

P. T. Nastos and A. Matzarakis, “Weather impacts on respiratory infections in Athens, Greece,” International Journal of Biometeorology, vol. 50, no. 6, pp. 358–369, 2006.

B. East, “Mean annual hours of sunshine and the incidence dental caries,” American Journal of Public Health, vol. 29, pp. 777–780, 1939.

J. McGrath, J. Selten, and D. Chant, “Long-term trends in sunshine duration and its association with schizophrenia birth rates and age at first registration—data from Australia and the Netherlands,” Schizophrenia Research, vol. 54, no. 3, pp. 199–212, 2002.

B. Aksoy, “Variations and trends in global solar radiation for Turkey,” Theoretical and Applied Climatology, vol. 58, no. 1-2, pp. 71–77, 1997.

S. C. Liu, C. H. Wang, C. J. Shiu, H. Chang, C. Hsiao, and S. Liaw, “Reduction in sunshine duration over Taiwan: causes and implications,” Terrestrial, Atmospheric and Oceanic Sciences, vol. 13, no. 4, pp. 523–545, 2002.

E. Pallé and C. J. Butler, “Sunshine records from Ireland: cloud factors and possible links to solar activity and cosmic rays,” International Journal of Climatology, vol. 21, no. 6, pp. 709–729, 2001.

G. Stanhill and S. Cohen, “Solar radiation changes in the United States during the twentieth century: evidence from sunshine duration measurements,” Journal of Climate, vol. 18, no. 10, pp. 1503–1512, 2005.

S. J. Reddy, “An empirical method for estimating sunshine from total cloud amount,” Solar Energy, vol. 15, no. 4, pp. 281–285, 1974.

A. S. N. Raju and K. K. Kumar, “Comparison of point cloudiness and sunshine derived cloud cover in India,” Pure and Applied Geophysics, vol. 120, no. 3, pp. 495–502, 1982.

S. Rangarajan, M. S. Swaminathan, and A. Mani, “Computation of solar radiation from observations of cloud cover,” Solar Energy, vol. 32, no. 4, pp. 553–556, 1984.

A. W. Harrison and C. A. Coombes, “Empirical relationship of cloud shade to point cloudiness (Canada),” Solar Energy, vol. 37, no. 6, pp. 417–421, 1986.

X. Yin, “Bright sunshine duration in relation to precipitation, air temperature and geographic location,” Theoretical and Applied Climatology, vol. 64, no. 1-2, pp. 61–68, 1999.

M. El-Metwally, “Sunshine and global solar radiation estimation at different sites in Egypt,” Journal of Atmospheric and Solar-Terrestrial Physics, vol. 67, no. 14, pp. 1331–1342, 2005.

A. P. Matzarakis and V. D. Katsoulis, “Sunshine duration hours over the Greek region,” Theoretical and Applied Climatology, vol. 83, no. 1–4, pp. 107–120, 2006.

H. M. Kandirmaz, “A model for the estimation of the daily global sunshine duration from meteorological geostationary satellite data,” International Journal of Remote Sensing, vol. 27, no. 22, pp. 5061–5071, 2006.

S. M. Robaa, “Evaluation of sunshine duration from cloud data in Egypt,” Energy, vol. 33, no. 5, pp. 785–795, 2008.

H. M. Kandirmaz and K. Kaba, “Estimation of daily sunshine duration from Terra and Aqua Modis data,” Advences in Meteorology, vol. 2014, Article ID 613267, 9 pages, 2014.

S. Khoshjavan, M. Mazlumi, B. Rezai, and M. Rezai, “Estimation of hardgrove grindability index (HGl) based on the coal chemical properties using artifical neural networks,” Oriental Journal of Chemistry, vol. 26, no. 4, pp. 1271–1280, 2010.

Z. K. Yip and M. K. Yau, “Application of artificial neural networks on North Atlantic tropical cyclogenesis potential index in climate change,” Journal of Atmospheric and Oceanic Technology, vol. 29, no. 9, pp. 1202–1220, 2012.

B. M. Ozyildirim and M. Avci, “Generalized classifier neural network,” Neural Networks, vol. 39, pp. 18–26, 2013.

F. Özkan, “Comparing the forecasting performance of neural network and purchasing power parity: the case of Turkey,” Economic Modelling, vol. 31, no. 1, pp. 752–758, 2013.

S. M. Al-Alawi and H. A. Al-Hinai, “An ANN-based approach for predicting global radiation in locations with no direct measurement instrumentation,” Renewable Energy, vol. 14, no. 1–4, pp. 199–204, 1998.

M. Mohandes, S. Rehman, and T. O. Halawani, “Estimation of global solar radiation using artificial neural networks,” Renewable Energy, vol. 14, no. 1–4, pp. 179–184, 1998.

A. Sözen, E. Arcaklioglu, and M. Özalp, “Estimation of solar potential in Turkey by artificial neural networks using meteorological and geographical data,” Energy Conversion and Management, vol. 45, no. 18-19, pp. 3033–3052, 2004.

F. S. Tymvios, C. P. Jacovides, S. C. Michaelides, and C. Scouteli, “Comparative study of Ångström's and artificial neural networks' methodologies in estimating global solar radiation,” Solar Energy, vol. 78, no. 6, pp. 752–762, 2005.

J. Mubiru and E. J. K. B. Banda, “Estimation of monthly average daily global solar irradiation using artificial neural networks,” Solar Energy, vol. 82, no. 2, pp. 181–187, 2008.

J. L. Bosch, G. López, and F. J. Batlles, “Daily solar irradiation estimation over a mountainous area using artificial neural networks,” Renewable Energy, vol. 33, no. 7, pp. 1622–1628, 2008.

D. A. Fadare, “Modelling of solar energy potential in Nigeria using an artificial neural network model,” Applied Energy, vol. 86, no. 9, pp. 1410–1422, 2009.

N. Lu, J. Qin, K. Yang, and J. Sun, “A simple and efficient algorithm to estimate daily global solar radiation from geostationary satellite data,” Energy, vol. 36, no. 5, pp. 3179–3188, 2011.

A. Moreno, M. A. Gilabert, and B. Martínez, “Mapping daily global solar irradiation over Spain: a comparative study of selected approaches,” Solar Energy, vol. 85, no. 9, pp. 2072–2084, 2011.

A. Ouammi, D. Zejli, H. Dagdougui, and R. Benchrifa, “Artificial neural network analysis of Moroccan solar potential,” Renewable & Sustainable Energy Reviews, vol. 16, no. 7, pp. 4876–4889, 2012.

M. Rumbayan, A. Abudureyimu, and K. Nagasaka, “Mapping of solar energy potential in Indonesia using artificial neural network and geographical information system,” Renewable and Sustainable Energy Reviews, vol. 16, no. 3, pp. 1437–1449, 2012.

A. Rahimikhoob, S. M. R. Behbahani, and M. E. Banihabib, “Comparative study of statistical and artificial neural network's methodologies for deriving global solar radiation from NOAA satellite images,” International Journal of Climatology, vol. 33, no. 2, pp. 480–486, 2013.

Y. Eissa, P. R. Marpu, I. Gherboudj, H. Ghedira, T. B. M. J. Ouarda, and M. Chiesa, “Artificial neural network based model for retrieval of the direct normal, diffuse horizontal and global horizontal irradiances using SEVIRI images,” Solar Energy, vol. 89, pp. 1–16, 2013.

T. Khatib, A. Mohamed, K. Sopian, and M. Mahmoud, “Assessment of artificial neural networks for hourly solar radiation prediction,” International Journal of Photoenergy, vol. 2012, Article ID 946890, 7 pages, 2012.

J. A. Jervase, A. Al-Lawati, and A. S. S. Dorvlo, “Contour maps for sunshine ratio for Oman using radial basis function generated data,” Renewable Energy, vol. 28, no. 3, pp. 487–497, 2003.

M. A. Mohandes and S. Rehman, “Estimation of sunshine duration in Saudi Arabia,” Journal of Renewable and Sustainable Energy, vol. 5, no. 3, Article ID 033128, 2013.

B. Aksoy, “Analysis of changes in sunshine duration data for Ankara, Turkey,” Theoretical and Applied Climatology, vol. 64, no. 3-4, pp. 229–237, 1999.

A. D. Şahin, “A new formulation for solar irradiation and sunshine duration estimation,” International Journal of Energy Research, vol. 31, no. 2, pp. 109–118, 2007.

U. Yildirim, I. O. Yilmaz, and B. G. Akinoğlu, “Trend analysis of 41 years of sunshine duration data for Turkey,” Turkish Journal of Engineering and Environmental Sciences, vol. 37, pp. 286–305, 2013.

S. Erinc, Klimatoloji ve Metotları, Gur-ay Pres, Istanbul, Turkey, 3rd edition, 1984.

General Directorate of Electrical Power Resources Survey and Development Administration (EIE), Turkey, 2014, http://www.eie.gov.tr/eie-web/turkce/YEK/gunes/tgunes.html

D. V. Hoyt, “Percent of possible sunshine and the total cloud cover,” Monthly Weather Review, vol. 105, pp. 648–652, 1977.

D. Matuszko, “Influence of cloudiness on sunshine duration,” International Journal of Climatology, vol. 32, no. 10, pp. 1527–1536, 2012.

G. R. Weber, “On the seasonal variation of local relationships between temperature, temperature range, sunshine and cloudiness,” Theoretical and Applied Climatology, vol. 50, no. 1-2, pp. 15–22, 1994.

M. Rebetez and M. Beniston, “Changes in sunshine duration are correlated with changes in daily temperature range this century: an analysis of Swiss climatological data,” Geophysical Research Letters, vol. 25, no. 19, pp. 3611–3613, 1998.

Y. H. Yang, N. Zhao, X. H. Hao, and C. Q. Li, “Decreasing trend of sunshine hours and related driving forces in North China,” Theoretical and Applied Climatology, vol. 97, no. 1-2, pp. 91–98, 2009.

M. Dolinar, “Spatial interpolation of sunshine duration in Slovenia,” Meteorological Applications, vol. 13, no. 4, pp. 375–384, 2006.

J. A. Duffie and W. A. Beckman, Solar Engineering of Thermal Processes, John Wiley & Sons, New York, NY, USA, 1991.

A. Hasni, A. Sehli, B. Draoui, A. Bassou, and B. Ameur, “Estimating global solar radiation using artificial neural network and climate data in the south-western region of Algeria,” Energy Procedia, vol. 18, pp. 531–537, 2012.

D. F. Specht, “A general regression neural network,” IEEE Transactions on Neural Networks, vol. 2, no. 6, pp. 568–576, 1991.

A. Amrouche and J. M. Rouvaen, “Efficient system for speech recognition using general regression neural network,” International Journal of Computer Systems Science and Engineering, vol. 1, pp. 183–189, 2006.

B. Asad, D. Zhijiang, S. Lining, K. Reza, and M. A. Fereidoun, “Fast 3D reconstruction of ultrasonic images based on generalized regression neural network,” IMFBE Proceedings, vol. 14, no. 5, pp. 3117–3120, 2007.

M. Firat and M. Gungor, “Generalized regression neural networks and feed forward neural networks for prediction of scour depth around bridge piers,” Advances in Engineering Software, vol. 40, no. 8, pp. 731–737, 2009.

K. Kayaer and T. Yildirim, “Medical diagnosis on pima indian diabetes using general regression neural networks,” in Proceedings of the Artificial Neural Networks and Neural Information Processing (ICANN/ICONIP '03), pp. 181–184, Turkey, Istanbul, Turkey, 2003.

T. Kiyan and T. Yildirim, “Breast cancer diagnosis using statistical neural networks,” Journal of Electrical & Electronics Engineering, vol. 4, no. 2, pp. 1149–1153, 2004.

I. Popescu, I. A. Kanatas, A. P. Constantinou, and I. Nafornita, “Application of general regression neural networks for path loss prediction,” in Proceedings of International Workshop Trends and Recent Achievements in Information Technology, 2002.

S. Ren, D. Yang, F. Ji, and X. Tian, “Application of generalized regression neural network in prediction of cement properties,” in Proceedings of the International Conference on Computer Design and Applications (ICCDA '10), pp. V2-385–V2-388, Qinhuangdao, China, June 2010.

Z. Wang and H. Sheng, “Rainfall prediction using generalized regression neural network: case study Zhengzhou,” in Proceedings of the International Conference on Computational and Information Sciences (ICCIS '10), vol. 2, pp. 1265–1268, Chengdu, China, December 2010.

T. Yildirim and H. K. Cigizoğlu, “Comparison of generalized regression neural network and MLP performances on hydrologic data forecasting,” in Proceedings of the 9th International Conference on Neural Information Processing (ICONIP '02), vol. 5, pp. 2488–2491, 2002.

B. Erkmen and T. Yildirim, “Improving classification performance of sonar targets by applying general regression neural network with PCA,” Expert Systems with Applications, vol. 35, no. 1-2, pp. 472–475, 2008.

E. Al-Daoud, “A comparison between three neural network models for classification problems,” Journal of Artificial Intelligence, vol. 2, no. 2, pp. 56–64, 2009.

S. Haykin, Neural Networks: A Comprehensive Foundation, Prentice Hall, Upper Saddle River, NJ, USA, 2nd edition, 1999.

A. S. S. Dorvlo and D. B. Ampratwum, “Summary climatic data for solar technology development in Oman,” Renewable Energy, vol. 14, no. 1–4, pp. 255–262, 1998.

C. Stanghellini, “A simple method for evaluating sunshine duration by cloudiness observations,” Journal of Applied Meteorology, vol. 20, no. 3, pp. 320–323, 1981.