Category Archives: EuroSun2008-13

Data transformation

In this study the generation of the variable which is going to be used to forecast future values of half daily solar global radiation is done taking into account gaussian and stationary properties needed in the time series to be used by predictive methods. This approach is based on “lost solar component” (Fig.

1.) defined previously.

Lost component time series presents higher variation in central months of the year due to higher solar radiation levels and bigger influence of cloudiness reflection and absorption on solar irradiance.

Elimination of the trend component is done differencing between successive time steps of lost component time series.

Подпись:W/m2 half day


In the other side, synoptic predictions of future states of sky conditions by statistical post-processing models presents less uncertainty than direct output of global solar irradiance in numerical weather prediction models. The classification of sky conditions is done by national weather services in six levels Fig. 2. As sky condition variable wasn’t available directly from AEMet it has been simulated dividing measured solar radiation in six levels. Future variation of sky conditions is combined with the difference of lost component time series as input to predictive methods. Besides, pure statistical models have the advantage of limiting the upper error of prediction and improving predictions errors of persistence by making relations of patterns of past observations and future values. In the next sections, results of using and not using sky condition as input to neural network is shown.

Numerical weather prediction models

The NWP models are able to provide the solution for seven atmospheric parameters, by solving the momentum, mass and energy conservation equations related to the motion of air and water vapor in the atmosphere. These models are also able to estimate the cloud cover, and incoming solar radiation [11]. By the model GFS (Global Forecasting System) [12], numerical modeling are performed in a (0.5 x 0.5)° earth surface grid with sampling interval of 3h. In order to improve the performance of local forecasts, the data of the global model are assimilated by regional NWP models. By using the hydrostatic model ETA with grid resolution of (0.4 x 0.4)°, Guarnieri [13] obtained a RMSE of 43.9 % and 43.6 % for the daily total of incoming solar radiation, for two different sites in Brazil.

MM5 Setup

The spatial configuration of the MM5 used in this work consisted in five nested domains and twenty four unevenly spaced sigma levels. The domains has an horizontal resolution of 71, 27,9, 3 and 1 km, respectively. The results of the last domain, of 1km resolution, were finally evaluated. Two-ways nesting was used to feed the information between domains. Atmospheric initial and boundary conditions were extracted from the analysis produce by the NCEP. The physical parameterizations used in the simulation were: the GRELL scheme for the cumulus parameterization, the MRF for the planetary boundary layer, the MIXED PASHED for the explicit moisture, the FIVE-LAYER SOIL for the soil model and the RRTM for the radiation scheme. This configuration is maintained for all the integrations. Finally, a 24 hours spin-up period was used in each 72-hours integration.

The simulations were carried out for a set of days selected along the year 2005. Particularly, four sets (one for each season of the year) of three consecutive days with clear-sky conditions were selected, one for each season of the year. As highlighted earlier, the aim of this work is to evaluate the ability of the MM5 to estimate the solar radiation in a complex topography area. To this end, two simulations were carried out for each season of the year. In one of the simulations (called hereinafter T, Topography), the solar radiation were computed using the MM5 subroutines OROSHAW and LEVSLP. These subroutines allows taking into account the effect of the slope, angle and shadow cast caused by the topography on the solar radiation estimates at the earth surface. In the second simulation, these subroutines were not used and, therefore, the MM5 solar radiation estimates do not account for these topographic effects. We will call these simulations as Not Topography (NT).

The MM5 estimates were evaluated in terms of the Mean Error (ME) and the Root-Mean-Square Error (RMSE). The ME quantified the overall bias and detected if the model is producing overestimation or underestimation, while the RMSE accounts for the spread of the error distribution. Al error estimates are computed using hourly values along the whole simulated period.

3. Results and conclusions

Table 2 to 5 shows the evaluation results for, respectively, winter, spring, summer and autumn,. The results are just presented for three representative stations: stations 4, 5 and 11. These stations have a important slope and represent different aspects (stations 4 aspect east, stations 5 aspect west and station 11 aspect south) (Table 1). Evaluation results present the comparison of the MM5 estimates using (T) and not using (NT) the topographic parameterization against the measured horizontal solar radiation (H). For instance, the NT-H notation in Tables 2 to 5 stands for the evaluation of the MM5 estimates not using topographic parameterization against the global horizontal radiation measured values.

Overall, the MM5 shows considerable skills in estimating the solar radiation under clear-sky conditions along the whole year, even for this complex topography area under study. Particularly, the lowest RMSEs are found in summer (~20%) and the highest during winter (more than 30%). Additionally, from the analysis of the ME values, it could be concluded the existence of a general tendency to overestimation in winter, spring and summer and to underestimation in autumn.

Regarding the topographic parameterization, tables 2 to 5 shows and overall improvement in the estimates, although this improvement strongly depends on the season of the year and the topographic characteristics of the location under study. Regarding the season of the year and as can be expected, the most important improvement takes place in winter (Table 2). The relatively low sun elevation angles during this seasons makes the topographic influence on the solar radiation measured ant the surface more important. The use of the topographic parameterization improves the MM5 estimates, in terms of RMSE, ranging from 10% of improvement in station 11 to less than 5% in the stations 4 and 5. This difference can be explained by the fact of that the station 11 has south aspect, while the station 4 has a west aspect and station 5 an east aspect. Therefore the station 11 receives more radiation than the other two and the potential improvement is higher. For the rest of the seasons, the improvement in the estimates provides by the MM5 topographic parameterization are lower than for the winter. Particularly, in spring and autumn, Tables 3 and 5, the improvement in the estimates in terms of the RMSE ranges from 1% to 7%. Particularly, the most important improvement (7%) is found during autumn and for station 5. Similar results are found in terms of the ME. During summer (Table 4), and as can be expected, the improvements in the estimates provide by the topographic parameterization are modest. The only important improvement (4% in terms of RMSE) is found for station 5.


RMSE (Wm-2) (%) ME (Wm-2) (%)











-61.1 (-11.9)





















Table 4. As in Table 2 but for summer.


RMSE (Wm-2) (%) ME (Wm-2) (%)
































Table 5. As in Table 2 but for autumn.



[1] G. A. Grell, J. Dudhia and D. R: Stauffer, (1994). A description of the fifth-generation Penn State/NCAR Mesoscale Model (MM5). Tech. Rep. NCAR/TN-398+STR, National Center for Atmospheric Research.

Analysis of the incoming solar radiation and other significant. parameters to estimate the solar resource at eight sites in Galicia (NW Spain)

Alberto Pettazzi, Jose Antonio Souto Gonzalez*

Department of Chemical Engineering, University of Santiago de Compostela, Lope Gomez de Marzoa St.,

Campus Sur. 15782 Santiago de Compostela, Spain
Corresponding author, ia. souto@,usc. es


The evaluation of the solar resource in a region needs the characterization of different meteorological parameters. In this work, an analysis of global irradiation, clearness index and sunshine hours has been carried out in Galicia (NW Spain). Mean values of irradiation are compared with those derived from satellite images. The analysis is completed with the evaluation of other related parameters: precipitations, temperature and relative humidity. A complete characterization of the solar resource has been achieved in eight areas of Galicia. Results show high variability in the solar irradiation values, describing a more complex pattern than that derived from satellite observations. Evidences of the presence of local atmospheric phenomena, that overlap the diurnal evolution of solar radiation, drive to the conclusion that the distribution of the solar resource in this region cannot be explained without the consideration of climatologic and topographical features.

Keywords: solar resource, ground-satellite comparison

1. Introduction

The development of technologies for renewable energies is a priority target in the policy from different fields. Between the renewable energies, solar thermal and solar photovoltaic energies have a great potential in Spain and in Galicia too [1]. In Spain, the only official reference of the solar radiation climatology is the Atlas of Solar Radiation [2] that collects monthly means of sunshine hours and global irradiation between 1951 and 1983. Unfortunately, the low resolution of its results and their oldness may not represent a valid reference for estimating the solar resource in the present. Other works [3] applies satellite images to obtain solar maps for Spain, but they lack of homogeneous validation with ground measurements.

In Galicia, only an attempt was done, by Vazquez et al. [4]. Using the method Heliosat-2 [5], they collected in an atlas the results of monthly means of solar radiation from observations by satellite Meteosat-6. Nevertheless, the analysis covered a short period of time and the results obtained were only compared with the measurements proceeding from one station. Finally, no information about other relevant parameters was provided.

This work is focused on the evaluation of the solar resource in Galicia from the ground. Over a 5- years period, the surface meteorological parameters that affect the performance of thermal and photovoltaic systems are analysed, keeping in mind the peculiarities of the Galician region. To accomplish this task, global irradiation, sunshine hours and clearness index have been evaluated at

eight locations, using daily measurements. Mean values of global irradiation proceeding from the stations have been compared with those derived from satellite observations. The evaluation is completed with the analysis of the time series of complementary parameters, such as precipitation, temperature and relative humidity, in order to obtain the most complete information about solar resource at these locations.

. Locally linear prediction

The easiest nonlinear method of local prediction was developed by Lorenz [4] and called method of analogies. Let be xt a point within the d-dimensional phase space, the predicted value a time T later, xt+T, will belong to some kind of interpolation between known points xt1+T, …, xtr+T, where x;7, …, xtr are the nearest r points to xt. Farmer and Sidorowich [12] propose an enhanced method by making a least-squared fit of a local linear map of xt into xt+T. Doing T = 1, we can obtain the predicted value as x’t+1 = atxt + bt by minimising:

W = SI lx+i- atxt — bt |2 (1)


in regard to at and bt. Nt is the s-neighborhood of xt.

Подпись: E2 Подпись: < E2 > Подпись: <[xt - xt]2 > <[xt -< xt >]2 > Подпись: (2)

Casdagli [10] suggests to use these models as a test for nonlinearity prediction, obtaining the average forecast error given by:

where <•> is the arithmetic average. If E = 0 then the prediction is perfect, and if E = 1, the prediction is no better than a constant predictor in the time series average. If E is calculated for different s-values, occurring the optimum for large neighbourhood sizes, this will indicate that the data are better described by a linear stochastic process in this embedding space. On the other hand, an optimum at rather small sizes will point out the existence of a nonlinear equation of motion.

The average forecast error of local linear models is plotted as a function of the neighbourhood size s in Fig. 5. Note the high value of E and how the minimum value of E occurs at higher neighborhood sizes. This result is associated with stochastic processes.


Fig. 5. Average forecast error of the local linear modelas a function of the neighborhood size fot the

differenced time series.

Solar Climate of Azores:. results of monitoring at Faial and Terceira islands

Ricardo Aguiar *, Antonio Rocha e Silva and Ricardo Coelho

INETI, Department of Renewable Energies
Lisbon, Portugal

* Tel: +351 210 924 602, E-Mail: ricardo. aguiar@ineti. pt

1. Introduction

The North Atlantic archipelago of Azores (Portugal) consists of nine islands, located about 1,500 km from mainland Europe. Currently the Azores has about 243,000 inhabitants, which depend heavily on imported fossil fuel for their energy supply. This is a concern for the regional Government and its Agency for Energy and Environment, ARENA, which naturally support energy efficiency and use of renewable energies. The Azores are rich in several renewable resources, and have even pioneered Portuguese exploitation of wind, wave, and high enthalpy geothermal energy for electricity generation. However, so far solar energy has not been a priority in the renewable energy panorama, probably because the Azores climate has been considered too cloudy. Nevertheless, this is true only in comparison with the Portuguese mainland, as the radiation levels are probably similar or even better than those of Northern Europe.

In this context ARENA is developing efforts for increasing the use of solar systems. A collaboration, partly financed by an INTERREG III B Project of the European Community, has been established with INETI, the Portuguese Public Laboratory for the area of Energy, to improve the solar climatology of the Azores. This climatology is indeed quite incomplete. It consists mainly on daily measurements for two sites only at the S. Miguel and Terceira islands. Some sunshine records are available from Campbell-Stokes heliographs. Satellite data exists but their quality is uncertain due to the small size of the islands and the shallow view angle.

2. Methods

A measuring campaign was established to gather data on solar radiation, having in mind the data needs for sizing solar thermal and photovoltaic systems in the first place, but also for thermal simulation of buildings. Four meteorological stations were equipped to measure horizontal global solar radiation, ambient temperature and relative humidity. Unfortunately the financial and human resources available did not allow for diffuse radiation sensors.

The first phase of the work addressed the two islands of the eastern group, Sta. Maria and S. Miguel. The current work describes preliminary results for the second phase of the work, which concerned two of the five islands of the central group, Faial and Terceira, see Figure 1 and Table 1. Other efforts are now addressing a third island of this central group, Graciosa.


Figure 1 — Location of monitoring stations.

Station sites were selected according to criteria that included island coverage, representation of microclimatic conditions, safety and reliance of the stations, low horizon obstructions, and placement near populated areas, i. e. those with potential solar energy users. Typical vistas can be appreciated in Fig. 2.

The equipment installed includes Rotronic Hygroclip S3 thermohygrometers, Huksefulx LP 02 pyranometers and Campbell Scientific CR510 dataloggers. The pyranometers were calibrated by early 2005 against a secondary standard (Kipp & Zonnen CM 11) held by the Laboratory for Solar Collector Testing at the INETI Campus, Lisbon.

Tabela 1 — Station sites.






Latitude (N)





Basic School Padre Jeronimo Emilio Andrade

27° 12′ 52.0»

38° 39′ 21.9»




Basic School of Biscoitos

27° 15′ 58.0»

38° 47′ 35.4»




Secondary School Doutor Manuel de Arriaga

28° 37′ 39.6»

38° 31′ 50.4»




Basic School

Prof. Constantino Magno Amaral Jdnior

28° 41′ 34.2»

38° 38′ 9.4»



Fig. 2 — Typical mounting of instruments and vistas around the stations.

3. Results

Measurements started by November 2005 and still continue, but the work reports only on the results obtained up to October 2007 (thus about two years).

Figure 3 shows recorded data at Angra, for the year of 2006, as an example. For solar radiation, the extraterrestrial values and clear-sky estimates are also shown. Clear-sky estimates were obtained with the ESRA model (ESRA, 2000) assuming a value of 2 for the Linke turbidity for airmass 2, TLK(2). Figure 4. shows monthly data for the 4 stations.

It is remarked that sequences of clear sky days — and even single clear sky days — are rare situations. On the contrary there are numerous overcast days. The seasonality of radiation data is low compared with continental sites. This feature is especially visible when the clearness index data is analysed: most of the daily values remain in a narrow range of 0.4 to 0.5, complemented with lower values during the Winter and higher values during the Summer, so at the seasonal level, the oscillation is from 0.2 to 0.55 (Fig. 4).

Temperature series show the well known characteristics for relatively small islands, of low thermal amplitude at daily (see Fig. 5) and seasonal levels. Also relative humidity levels are, as expected, quite high and show little variability.




Fig. 3 — Daily records of solar radiation (top, kWh/m2), ambient temperature (middle, °C), and

relative humidity (bottom, %)

In Figure 5 it can be appreciated that the daily mean profiles of hourly solar radiation are almost symmetrical in respect to solar noon. This is somewhat different from the pattern for continental stations. For the later, it is usual that the afternoon values are a bit smaller than the corresponding morning values, due to larger cloud sizes brought about by the accelerated temperature related convection.


image023 Подпись: (2005-2007)

Monthly Global Horizontal Solar Radiation (kWh/m2)

Jan Feb Mar Apr May Jun Jul Aug Sep Oct Nov Dec

image025 image026 Подпись: 2005-2007

Monthly Clearness Index

Jan Feb Mar Apr May Jun Jul Aug Sep Oct Nov Dec





Fig. 5 — Typical mean daily profiles for 2 stations (May 2006-2007).

4. Climatology

Long term records of sunshine hours are available for some classic meteorological stations operated by the Portuguese Weather Services. It is possible to estimate long term solar radiation data using the well known relationships of relative sunshine and clearness index at the monthly scale.

As simultaneous radiation and sunshine data was not available throughout the Azores region, the usual Angstrom relationship is difficult to apply. The approach taken was instead the one proposed by Suehrcke (e. g. Suehrcke, 2000; Tiba et al., 2002)

Подпись: [1](<Gm> /<Gcm>)2 = <Sm> / <S()m>

where <> denotes long term averages, index m a monthly mean, G is global horizontal solar radiation, Gc the corresponding clear sky value, S is bright sunshine hours and S0 the daylength.

Подпись: Fig. 6 - Adjusting the Linke turbidity for the Angra site.

A clear sky must be implemented to use the method. Figure 6 illustrates the approach used to obtain a Linke turbidity to feed in the ERA model. Maxima for each Julian day where computed from the data and these where fitted by adjusting the Linke turbidity. For many days and the short monitoring period, only cloudy values where available, so the method is not fully objective in this conditions. Nevertheless an estimate of TLK(2)=4 was obtained. This is higher than for continental unpolluted stations, which typically display values around 3.

Using these turbidity estimates and the recorded sunshine data for 1951-80 (IM, 2003), the climatology obtained is the one presented in Table 2.

The differences of the estimated climatology to the short-term monitoring mean data are surprisingly small. The Suehrcke approach points to long term data a bit above the short-term monitoring 2005-2007. The main difference appears for June, and it seems quite possible that this was because of atypical weather at the region, in particular for the Summer 2005.

Table 2 — Solar radiation climatology at Angra do Heroismo estimated from sunhsine records.






































































































5. Conclusions

Preliminary results of monitoring for 2005-2007 reveal that the climate at the four stations — and therefore at the two islands — is very similar. Daily thermal amplitude is low and relative humidity levels are high, features already expected taking into account the nearby ocean. Atmospheric turbidity apparently suffers from high levels of humidity and marine aerosols.

Some morning-afternoon asymmetry exists but it is small. The findings at the Azores islands suggest that most of the nebulosity there is related to macroscale eddies (the Azores Anticyclone in particular). This would also explain why the data at the northern and southern sites of the islands, and at various islands, are so similar. Indeed average yearly global horizontal irradiation at Terceira island was about 3.65 kWh/m2 per day, slightly less than for Faial island, 3.75 kWh/m2 per day.

Daily radiation sequences do not display long runs of overcast or sunny conditions, but frequent oscillations, reflected in low autocorrelation values of the detrended sequences. However, the most striking aspect is the small seasonality of the various parameters. This was expected for temperature and humidity due to the presence of the ocean and the small size of the islands, but is also true for solar radiation, especially when viewed from a clearness index perspective.

As a result of the monitoring, a clear sly model and a methodology to obtain long term climatic values could be obtained for these islands.

6. Further work

The meteorological records are now being processed to yield Reference Meteorological Years, with the objective of providing better meteorological data for the Portuguese governmental schemes supporting solar energy as well as for the Portuguese building regulation codes (subsidiary to the EU Directive on Energy Efficiency in Buildings). These require that both the summer thermal energy balance of buildings and a mandatory contribution of solar energy be computed through detailed numerical simulation with hourly data — currently obtained from stochastic models not yet tuned for these mid-Atlantic climates.

Acknowledgements — This work was partly supported by Project ERAMAC — Maximizagao da Penetragao das Energias Renovaveis e Utilizagao Racional da Energia nas Ilhas da Macaronesia, integrated in Programme INTERREG III-B, also supported by the Regional Government of Azores, and implemented by the Azores Regional Energy and Environment Agency, ARENA.


ESRA (2000). The European Solar Radiation Atlas. Scharmer, K. and J. Greif, Eds. ISBN 2-911762-21­5, 110 pp. Presses de l’Ecole des Mines de Paris, Chapter 3: Basics of Solar Radiation. J. Page and R. Aguiar.

IM (2003). Normais Climatologicas para 1961-90 publicadas pelo Institute de Meteorologia, obtidas via Projecto CLIMAAT, INTERREG_3B — Mac 2.3/A3.

Suehrcke, H. (2000) On the relationship betweeen duration of sunshine and solar radiation on the earth’s surface: Angstrom’s equation revisited, Solar Energy, 68, 5, 417-425.

Tiba, C., N. Fraidenraich e R. Aguiar (2002). Valor Acrescentado do Modelo de Suerhcke para a

Radiagao Global Mensal Avaliado com Dados Brasileiros. Actas do XI Congresso Iberico e VI Ibero Americano de Energia Solar, 29 Set. a 2 Out. 2002, Vilamoura. Edigao SPES, Lisboa.

Yearly Sum and Frequency distribution of DNI

The accuracy of the proposed method to derive DNI was evaluated for the year 2005 for six stations of the Spanish Meteorological Service INM, see table 1. The evaluation was performed for all hourly values with the sun above horizon.

Table 1: INM stations with measurements of direct normal irradiance DNI.














A Coruna















Special focus of the evaluation was on the accuracy of yearly sums and on the investigation of frequency distributions of the time series, as these are relevant for yield estimates. Beside this, the relative bias rBIAS and the relative root mean square error rRMSE are used for quantitative comparisons.

For hourly values an rRMSE of 14.7% for global irradiance and of 31% for direct normal irradiance was found for the Spanish stations. The rBIAS is 1.5% for global irradiance and 1.1% for direct irradiance respectively.

Figure 4 displays rBIAS and rRMSE of DNI for the Spanish stations. The rBIAS for one year and one station is equivalent to the relative deviation of the annual sum for one station. In addition the rBIAS for clear sky situations rBIASclearsky is provided. Hours are assigned as clear sky situations, if two criteria match: The clear sky index lies within 0.9<=k*< =1.1 and the variability of the cloud index in a small region of 3×5 pixel is small.

The comparison of rBIAS and rBIASciearsky illustrates the strong influence of the quality of the clear sky model and the atmospheric input parameters on the quality of the annual sums. The deviation of annual satellite derived irradiance sums from the respective ground measured sums between -2% and 8.5% for DNI.

In table 2 the accuracy information is given on different time scales which are relevant for solar energy applications. The deviation of annual satellite derived irradiance sums from the respective ground measured sums is between -2% and 4.5% for GHI and between -2% and 8.5% for DNI.

Table 2: Accuracy of satellite derived global horizontal irradiance GHI and direct normal irradiance DNI for

Подпись: Fig. 5: Frequency distribution of calculated and measured direct normal irradiance (left) and of global irradiance (right) for the 6 spanish stations.

different time scales.



hourly mean

335 W/m2

366 W/m:

rRMSE hourly



rRMSE daily



rRMSE monthly






For energy conversion systems with non-linear response to the irradiance input a correct representation of the frequency distribution is of special importance. Figure 5 displays the frequency distribution of calculated and measured DNI and GHI for the Spanish stations; a fairly good agreement is achieved.

3. Conclusion

We developed a method to derive direct normal irradiance from MSG data. The influence of clouds on the direct normal irradiance is derived from MSG data with high quality. The quality for clear skies is determined by the accuracy of the aerosol climatology.

4. Acknowledgements

Thanks are due to the Spanish Meteorological Service INM for the supply of ground data.


[1] Hammer A., Heinemann D., Hoyer C., Kuhlemann R., Lorenz E., Muller R., Beyer H. G. (2003): Solar energy assessment using remote sensing technologies. Remote Sensing of Environment, 86, 423-432.

[2] Fontoynont, M., Dumortier, D., Heinemann, D., Hammer, A., Olseth, J. A., Skartveit, A., Ineichen, P., Reise, C., Page, J., Roche, L., Beyer, H. G. and Wald, L.: 1998, ‘Satellight: A WWW Server Which Provides High Quality Daylight and Solar Radiation Data for Western and Central Europe’, Proceeding of the 9th Conference on Satellite Meteorology and Oceanography, Paris, France, 25.-29.05.1998, 434-437

[3] Muller, R. W., Dagestad, K. F., Ineichen, P., Schroedter, M., Cros, S., Dumortier, D., Kuhlemann, R., Olseth, J. A., Piernavieja, C., Reise, C., Wald, L., Heinemann, D. (2004): Rethinking satellite based solar irradiance modelling — The SOLIS clear-sky module. Remote Sensing of the Environment, Volume 91, Issue 2, 160-174

[4] Mayer, B. and Kylling, A. (2005): Technical Note: The libRadtran software package for radiative transfer calculations: Description and examples of use. Atmos. Chem. Phys., 5:1855-1877, 2005

[5] Kinne et al. (2005): An AeroCom Initial Assessment — Optical Properties in Aerosol Component Modules of Global Models. In: Atmospheric Chemistry and Physics Discussions, 5, 8285-8330.

[6] Kalnay, E. et al. (1996): The NMC/NCAR 40-year reanalysis project. Bull. Am. Meteorol. Soc. 77 (3), 437­472.

[7] Kemper A. (2007): Bestimmung der Diffusstrahlung unter Wolken aus Daten des Satelliten MSG. Diploma thesis, Oldenburg Unversity, Germany.

Comparison of meteorological data from different sources. for Bishkek city, Kyrgyzstan

Ruslan Botpaev[2]*, Alaibek Obozov1,

Janybek Orozaliev2, Christian Budig2, Klaus Vajen2,

1 Kyrgyz State Technical University, Department for Renewable Energies, Bishkek (Kyrgyzstan)
2 Kassel University, Institute of Thermal Energy Engineering, Kassel (Germany)
Corresponding Author, botpaev@inbox. ru


Meteorological data from different sources (a local weather station, Meteonorm 5.1 and 6.0 as well as own measurements) have been compared with each other for Bishkek, Kyrgyzstan. It was identified that Bishkek city in Meteonorm 5.1 and 6.0 has wrong altitude (2111 m instead of 760 m), which leads to significantly lower ambient temperatures. Therefore, it is necessary to define Bishkek manually as a new site with the correct coordinates. In this case the ambient temperatures are in agreement with those measured by the weather station Frunze. Annual global solar irradiation from own measurements («1500 kWh/m2a) and weather station Frunze (1572 kWh/m2a) are in a good agreement, while its values from Meteonorm 5.1 and 6.0 are approx. 20% lower. Furthermore, monthly sums of global solar irradiation generated with Meteonorm 5.1 have an untypical trend in summer, having two local maximum points. This is inconsistent with irradiation data from other sources and with the sunshine duration in the relevant period. The data from Meteonorm should be proved on plausibility, particularly for sites or stations in developing countries.

1. Introduction

For designing and dimensioning of solar plants and prediction of solar gains it is important to have precise meteorological data, i. e. data of solar radiation, ambient temperatures and wind velocity. Meteorological data for different places can be taken from own measurements, reference books, weather stations or generated with special software (e. g. Meteonorm1).

For Bishkek, the capital of the Kyrgyz Republic in Central Asia, basically, three sources of meteorological data are available: a local weather station Frunze, a software Meteonorm (only versions 5.1 and 6.0 are considered in this study) and own measurement data (only radiation). Geographical coordinates of Bishkek are 74.5°E longitude, 42.8°N latitude (comparable with Rome) and 760 m altitude. The climate is strong continental (hot summer, cold winter, high number of sunny days).

The objective of this investigation is to compare these sources of meteorological data with each other regarding average annual ambient temperatures and solar radiation data.

Artificial Neural Networks

Predictive methods can be divided fundamentally in Numerical Weather Prediction (NWP) models and statistical techniques. Time series prediction belongs to statistical techniques group, which is based on obtaining future values of the variable to predict in function of past observations.

Artificial Neural Networks (ANNs) paradigm was originally an abstract simulation of neural biological systems which tried to synthetic its abilities [3].

Supervised training networks have been the neural network models more developed from the beginning of its design. Training data is constructed by various input and output patterns. Knowing the output is a fact from which the training benefits as a master which supervises and modify its parameters.

Neural network used is backpropagation type; this class of supervised network has a learning algorithm based on error correction method. Learning in backpropagation networks is done in the following way: an input pattern is selected from the set of input patterns, it determines the activations of input nodes, latter, neurons are activated from hidden layers and finally output units are computed. The activation pattern which links the output layer is then compared with output pattern associated to the input pattern to calculate the error [4].

3. Results

The methods described in the previous section are trained with 75% of the data set. Assessment of the models proposed is performed over the 25% of the clearness sky separated sequentially for this purpose.

Neural Network (NN) is essayed with different layer sizes and neurons in each layer with the following architecture:

• One layer and one neuron.

• Two layers with three neurons in input layer and one neuron in output neuron.

• Three layers with five neurons in input layer, three neurons in hidden layer and one neuron in output neuron.

• Four layers with eight neurons in input layer, five and three neurons in hidden layers and one neuron in output neuron.

• Five layers with fifteen neurons in input layer, eight, five and three neurons in hidden layers and one neuron in output neuron.

Training algorithm is Leverage-Marquardt algorithm [3].

image034 Подпись: / RMSD Подпись: (1)

Accuracy of each one of the models presented is studied in terms of relative Mean Bias Deviation (MBD) and Root Mean Squared Deviation (RMSD):

being N the population size, Kt clearness sky index observed and Kt half daily clearness index predicted.

Each one of the models essayed is compared with a basic model, persistence model (PER), which is based on assumption that clearness index for next temporal step is the same as last temporal observation. Comparisons are done in terms of relative RMSD. Improvements over persistence for deviation index (i-error) (in terms of RMSD) defined previously, is obtained from the following expression:

The neural network model, NN(z), has been tested with different input vector size which represents the temporal value of the difference of lost components for half daily lags z=1..10. Fig 4. shows results for percentage of error in term of RMSD for NN model with just one neuron and one layer and with different input vector size. Improvement of the model with lower RMSD (NN(10)) against persistence is shown in Fig. 5. Difference of lost component time series is a signal with high variability which can be seen in the results in the fact that from second time step of prediction the error gets an upper limit of approximately 26% of RMSD.


Fig. 4. RMSE Error: model NN with [1] layer and only lost component as input.


Fig. 5. Improvement NN(10) over Persistence: NN with [1] layer and only lost component as input.

The second model essayed is based on using as input vector the synoptic situation difference for the time step which is being predicted and the difference of lost component for half daily lags z=1..10. Fig 6. shows results for percentage of error in term of RMSD for NN model with four layers and eight,

five, three and one neuron in consecutive layers. Improvement of model with lower RMSD (NN(10)) against persistence is shown in Fig. 7.



Fig. 6. RMSE Error: model NN with [8 5 3 1] layer and lost component and AEMet prediction as input.



Fig. 7. Improvement NN(10) over Persistence: NN with [8 5 3 1] layer and only lost component and AEMet

prediction as input.


3. Conclusion

Properties of ground global solar radiation has been presented based on requirements of statistical prediction techniques. The solar radiation time series has been transformed taking into account gaussian and stationary properties needed by neural network model used to predict future values.

A new model to predict time series of ground global solar radiation based on previous transformation of solar radiation time series has been presented. The error of the model is limited by an upper level which is due to deterministic nonlinear behaviour of the signal that can’t be followed correctly by neural network model. The second model presented is based on using the difference of lost component and the difference of the synoptic situation for the time step of prediction. The error has a lower level of nine percent. The neural network used has a structure more complex than the network used in previous works.


[1] Veziroglu TN and dot a. 21st Century’s energy: Hydrogen energy system. Energy Conversion and Management 2008; 49:1820-1831.

[2] Patlitzianas KD, Ntotas K, Doukas H, and Psarras J. Assessing the renewable energy producers’ environment in EU accession member states. Energy Conversion and Management 2007; 48:890-897.

[3] Haykin S. Neural networks. A comprhensive foundation. New York (USA): 1994.

[4] Palit AK and Popovic D. Computational Intelligence in Time Series Forecasting: Theory and Engineering Applications. 2005.

Statistical methods for correcting forecast of NWP models

Due to its high uncertainties, the performances of the NWP models have to be improved. As stated by Wilks [14], the predictors are the input variables of a statistic black box model, which has the function to predict an output variable, named as predictand that is used for correction. The MOS can be categorized in two distinct methods that differ in its predictors. While the first method, which is applied in [11] and [13], utilizes as predictors the simulated weather variables of the NWP model, except the solar radiation, the second method utilizes as predictors a time series of the measured minus the forecasted variable as presented by Libonati [15]. The second method is thus

enabled for phase and amplitude corrections of a local weather forecast. The Multiple Linear Regression (MLR) as shown in [11] and [13] or an ANN [13] is used as statistic model for the first method and a Kalman Filter for the second method [12]. Correcting the non-hydrostatic model MM5, which has a resolution of (3 x 3) km, the author of [11] obtained for two different years an RMSE of 28 % and 30 % for the prediction of the total daily horizontal solar energy on a site in Germany. Correcting the hydrostatic model ETA, the author of [13] obtained an RMSE of

25.5 % and 25.6 % for the forecast of the same variable for two different sites in Brazil. At a particular day the radiation was predicted with 17 MJ/m2 (107% of the mean value), whereas the measured energy was only 2 MJ/m2. The author selected the used output variables of the ETA model by the application of a significance test based on the MLR. Substituting the MLR with an ANN model the author didn’t observe considerable improvements of the MOS using the selected variables as predictors. The author applied the obtained model to other sites to test the generality of model performance. For two different cases, he obtained 35.6% and 38.9%. The second MOS method, which utilizes the Kalman Filter model [15], wasn’t already applied for the forecast correction of the solar radiation. The Kalman Filter has the disadvantage that it cannot, in its standard version, handle nonlinear problems [16]. Even applied to strictly linear systems, this model has higher uncertainties compared to an ANN, as shown in [17].