Category Archives: EuroSun2008-13

Management and Exploitation of Solar Resource Knowledge

C. Hoyer-Klick1*, H. G. Beyer2, D. Dumortier3, M. Schroedter-Homscheidt4, L. Wald5, M.
Martinoli6, C. Schilings1, B. Gschwind5, L. Menard5, E. Gaboardi6, L. Ramirez-Santigosa7, J.
Polo7, T. Cebecauer8,T. Huld8, M. Suri8, M. de Blas9, E. Lorenz10, R. Pfatischer11, J. Remund12, P.

Ineichen13, A. Tsvetkov14, J. Hofierka15

1 German Aerospace Center (DLR), Institute of Technical Thermodynamics, Pfaffenwaldring 38-40, 70569

Stuttgart, Germany

2 Hochschule Magdeburg-Standal, Germany. 3Ecole Nationale des Travaus Publics de l’Etat (ENTPE), France
4 German Aerospace Center, German Remote Sensing Data Center, Germany 5 Ecole des Mines de Paris, France
6 Icons srl., Italy 7 CIEMAT, Spain 8 European Commission, Joint Research Center, Institute for Energy, Italy
9 Universidad Publica de Navarra, Spain 10 Oldenburg University, Energy — and Semiconductor Research Lab,
Germany 11 meteocontrol GmbH, Germany 12 Meteotest, Switzerland 13 University of Geneva, Switzerland
14 Voeikov Main Geophysical Observatory, World Radiation Data Center, Russia
15 University of Presov, Slovakia

Corresponding Author, carsten. hoyer-klick@dlr. de

Abstract

Knowledge of the solar energy resource is essential for the planning and operation of solar energy systems. In past years there has been substantial European and national funding to develop information systems on solar radiation data, leading to the situations that several data bases exist in parallel, developed by different approaches, various spatial and temporal coverages and resolutions including those exploiting satellite data. The user of these products may end up with different results for the same requested sites. To better guide the users, a benchmarking exercise is under preparation. A set of reference data has been collected and benchmarking measures and rules have been defined. The results of the benchmarking and the feedback from stakeholders will be integrated into a guide of best practices in the application of solar resource knowledge. Access to data has been quite fragmented. Each service has its own way of access to the data and delivery format. A new broker portal based on the experience of the project Soda aims to unify and ease the access to distributed data sources and applications providing solar resource information.

Keywords: Solar resource information, benchmarking, access to data, user guidance

1. Introduction

Knowledge of the solar energy resource is essential for the planning and operation of solar energy systems. In past years there has been substantial funding from the European Commission to develop information systems on solar radiation data, such as the European Solar Radiation Atlas (ESRA), the projects SoDa, Satel-Light, PVGIS, PVSAT, PVSAT-2 or Heliosat-3 and the Envisolar project of the European Space Agency (ESA). In addition national services were set up as Meteonorm by Meteotest in Switzerland and SOLEMI by DLR in Germany. The information on available databases and
integrated systems has been summarised and later updated in [1,2]. From the regional point of view, the projects focused mainly on Europe or its regions, leading to the situation that several different data bases exist in parallel developed by different approaches, various spatial and temporal coverages and different resolutions including those exploiting satellite data. The users comparing information from different data sources for the requested sites may end up with uncertainty that is difficult to deal with.

Large steps forward have been made for the benefit of research, renewable energy industry, policy making and the environment. Nevertheless, these multiple efforts have led to a fragmentation and uncoordinated access: different sources of information and solar radiation products are now available, but uncertainty about their quality remains. At the same time, communities of users lack common understanding how to exploit the developed knowledge.

The project MESoR started in June 2007 and aims at removing the uncertainty and improving the management of the solar energy resource knowledge. The results of past and present large-scale initiatives in Europe, will be integrated, standardised and disseminated in a harmonised way to facilitate their effective exploitation by stakeholders. The project will contribute to preparation of the future roadmap for R&D and strengthening the European position in the international field.

The project includes activities in user guidance (benchmarking of models and data sets; handbook of best practices), unification of access to information (use of advanced information technologies; offering one-stop-access to several databases), connecting to other initiatives (INSPIRE of the EU, POWER of the NASA, SHC and PVPS of the IEA, GMES/GEO) and to related scientific communities (energy, meteorology, geography, medicine, ecology), and information dissemination (stakeholders involvement, future R&D, communication).

First Steps in the Cross-Comparison of Solar Resource. Spatial Products in Europe

M. Suri1*, J. Remund2, T. Cebecauer1, D. Dumortier3, L. Wald4, T. Huld1 and P. Blanc4

1 European Commission, Joint Research Centre, Institute for Energy,

Renewable Energies Unit, via E. Fermi 2749, TP 450, I-21027 Ispra (VA), Italy
2 Meteotest, Fabrikstrasse 14, CH-3012 Bern, Switzerland
3 Ecole Nationale des Travaux Publics d’Etat, Departement Genie Civil et Batiment,

URA CNRS 1652, Rue Maurice Audin, F-69518 Vaulx-en-Velin, Cedex, France
4 MINES ParisTech, Centre Energetique et Procedes, BP 207, 06904 Sophia Antipolis Cedex, France

* Corresponding Author, marcel. suri@jrc. it

Abstract

Yearly sum of global irradiation is compared from six spatial (map) databases: ESRA, PVGIS, Meteonorm, Satel-Light, HelioCliom-2, and NASA SSE. This study does not identify the best database, but in a relative cross-comparison it points out to the areas of higher variability of outputs. Two maps are calculated to show an average of the yearly irradiation for horizontal surface together with the standard deviation that illustrates the combined effect of differences between the databases at the regional level. Differences at the local level are analysed on a set of 37 randomly selected points: global irradiation is calculated from subset of databases for southwards inclined (at 34°) and 2-axis tracking surfaces. Differences at the regional level indicate that within 90% of the study area the uncertainty of yearly global irradiation estimates (expressed by standard deviation) does not exceed 7% for horizontal surface, 8.3% for surface inclined at 34°, and 10% for 2-axis tracking surface. Higher differences in the outputs from the studied databases are found in complex climate conditions of mountains, along some coastal zones and in areas where solar radiation modelling cannot rely on sufficient density and quality of input data.

Keywords: solar radiation database, maps, benchmarking

1. Introduction

Solar energy technologies and energy simulation of buildings need high quality climatic data in the phase of localisation (siting), design, financing, and system operation and management. The choice of the best technological option depends among other things on the geographic region, as the performance of solar energy systems is influenced by solar resource and other climate parameters.

Several spatial databases of solar resource information are now available as a result of European and national projects. They have been developed from various data inputs, covering different time periods, where diverse approaches have been applied. Although quality assessments of the individual databases have been performed, no inter-comparison of the outputs was performed. When comparing various data sources, differences show up which is confusing, especially to users who are not fully aware of the uncertainties and the limits of data application. Therefore, better understanding of the geographic distribution and variability of solar resource in Europe is needed.

In this contribution we open a complex issue of benchmarking the solar radiation databases and underlying models for deriving information relevant to energy technology. We focus on a comparison of six spatial databases and integrated systems that offer solar resource and climate

data and energy-related services for Europe: Meteonorm [1], ESRA [2], Satel-Light [3], NASA SSE/RETScreen [4], HelioClim/SoDa [5], and PVGIS [6]. This list is not exhaustive, and in future also other databases may be considered, including those that cover smaller regions. We compare yearly sum of global irradiation as obtained by querying each database. Map analysis compares horizontal irradiation, while on a set of 37 randomly selected points we compare irradiation received by inclined and 2-axis tracking surfaces.

Meteorological data. Meteorological station “Frunze”

The central administrative board on hydrometeorology of the Ministry of Emergencies of the Kyrgyz Republic is entitled to measure meteorological and hydrological data in the Kyrgyz Republic. There are in total 31 weather stations and 75 hydrological stations. One of the weather stations “Frunze” is situated in the west part of Bishkek. The measurement equipment is remained from the USSR period. An actual value of global and diffuse solar radiation is measured 5 times a day at 6.30, 9.30, 12.30, 15.30 and 18.30. Till 1993 daily solar irradiation on horizontal surface was measured by an integrator. This device is, however, absent since 1993 for technical reasons. Therefore, since 1993 daily solar irradiation is estimated by linear interpolation of solar radiation between 5 measured points taking into account the time of sunrise and sunset (the so-called trapezium method). The central administrative board on hydrometeorology claims the accuracy of this method to be in the range of 10% for monthly sums of solar radiation.

Modeling and Analysis of Chinese Exposure to Solar Radiation Based on the Available Meteorological Data at CMA, China

Yu Qiang1* and Wang Zhifeng2

1,2 Institute of Electrical Engineering, CAS P. O.BOX 2703 Beijing 100190, China Corresponding Author, yuqiang1984@mail. iee. ac. cn

Abstract

The distribution of solar radiation is very important for choosing sites for solar power tower plants. In this paper the observed data which include solar radiation and sunshine duration of eight typical cities of China during the period 1994-2003 are used to establish a correlation equation between monthly average daily values of clearness and relative sunshine. The model is used to estimate the global solar radiation of the whole country. And it is proved to be good results (greater than 94% in most cases). The predictive efficiency of this model is also compared with some other models which are believed to be applicable globally in terms of mean percentage error (MPE), mean bias error (MBE) and root mean square error (RMSE). And the results prove that it is also better than that of those models.

Key Words: global solar radiation, sunshine duration, MBE, MABE, RMSE.

1. Introduction

The design of a solar energy conversion system must always start with a study of solar radiation data at a site. One of the most important requirements in the design is the information on the intensity of solar radiation at a given location [1]. Unfortunately, there are very few meteorological stations that measure global solar radiation. Solar radiation data are still very scarce, especially in developing countries. So we must consider other methods to calculate relative solar data for places where they are not directly measured, many attempts have been made to develop models and empirical correlations that can predict the amount of solar radiation available at a given location from a few input parameters.

While it has been proved that a number of commonly measurable atmospheric and meteorological parameters such as turbidity, relative humidity, degree of cloudiness, temperature and sunshine duration taken severally or jointly, affect the magnitude of the global radiation incident on a given location. And the preponderance of data point to the fact that the greatest influence is exerted by sunshine hours.

There are several correlations [2-7] to have been developed that predict the correlation between the global radiation and the percentage of bright sunshine hours in a simple linear regression form (the Angstrom-Prescott type). And some authors have also developed quadratic correlation [2, 4, 6] model and multiple linear regression. The study of this paper is to establish a linear regression form which uses the data of eight typical cities of China for estimation global solar radiation for the cities where there are no meteorological stations but have similar meteorological conditions.

Materials, methods and models

In the present work the solar radiation is forecasted with the non-hydrostatic model Advanced Regional Prediction System (ARPS). This model is providing its forecast weather variables for a horizontal grid of (0.12 x 0.12)° resolution with a sampling interval of 10 min. The model is simulated at the LEPTEN laboratory (Laboratory of Energy Conversion Process Engineering and Energy Technology), former LABSOLAR, at the Federal University of Santa Catarina. The simulation assimilates the data of the global reanalysis delivered by the National Center for Environmental Prediction (NCEP) [12]. The analysis data characterize the initial condition at every 6 h, necessary to operate ARPS in actual time. The reanalysis data represent improved analysis data of the atmosphere. Both the analysis and reanalysis data are based on atmospheric measurements and their interpolations, as well as on the last forecasts of the GFS, which can accomplish forecasts until a ten days horizon. The operational forecast uncertainty includes both the analysis and the forecast uncertainty. In the present article only the uncertainty based on the reanalysis are verified. Therefore the reanalysis data, based on the GFS model, is assimilated with the regional ARPS model in a 6 h interval. For the uncertainty verification the 24 h mean value of the downward short wave radiation of the ARPS output is compared to the measured mean value of the global radiation. In a second step a statistical correction for the reanalysis uncertainty is accomplished. In the text that follows the daily energy E [Wh/m2] is equal to the daily mean radiation H [W/m2] multiplied by 24 hours.

Ordinary kriging

Kriging refers to a family of least-square linear regression algorithms that attempt to predict values of a variable at locations where data are not available based on the spatial pattern of the available data. The description of kriging theory and its application are given in detail by [5]. A semivariogram, y(h) represents the spatial variability in the data and is defined as Eq. (1):

1 N(h) 2

Y(h) = E [Z(xj + h) — Z(x )]2, (1)

2N(h) 7=1

where N(h) is the number of pairs of points separated by lag distance h, Z(x) and Z(x+h) are random values at locations x and x+h.

In this study, the exponential model have been used to fit the sample semivariograms, this model parameterizes the semivariogram in the following way:

Y(h) = C0 + C1 [1 — exp(-h / a)], (2)

where C0, Cj, a are called nugget, sill, and range respectively.

The objective of ordinary kriging procedure is to estimate data values at unsampled locations x0 using information available elsewhere in the domain (x1, x2,……………………………………………………………………………………………… , xn) . This can be carried out

by expressing Z(x0)as a linear combination of the data Z(x1), Z(x2),…………………… Z(xn) :

Подпись:Z(x0) = EX7Z (x0).

The optimal weights “ A7 ” are calculated assuming that the estimation Z(x0) by Z(x0) is

unbiased, that is, the expected value of the estimates is the same as that of the known data. The

n

condition needed for unbiased estimator is 2 A = 1.

7=1 7

Monthly-averaged daily data

Подпись: 0

Figures 2 and 3 show the annual variation of global irradiation and sunshine hours for the stations considered in this work. Generally, they show a regular annual variation with a maximum between June and July and a minimum in December. The only exception is represented by the station of Pedro Murias, which has a maximum in May and a relative minimum in June. As expected, the highest value of global insolation is found at Ancares, due to its altitude and rural location. The evolution of the sunshine hours shows a similar trend, with maxima occurring in summertime and minima in winter.

image088

Table 2 summarizes these results. Having a look at the mean values of insolation and sunshine hours, Corrubedo has the highest values. It should be noticed that this site, which is located at only 38 km west from Lourizan, has mean sunshine hours and insolation significantly higher than Lourizan, 1 hour and 2 MJ m-2, respectively. The two stations located in an urban environment (Lourizan and Ferrol) are characterized by similar levels of insolation and sunshine hours.

Fig.3: Monthly average daily sunshine hours for the stations analysed

On table 3 are represented the main features of the yearly and monthly mean values of temperature, precipitations and relative humidity. The trend of the temperature is similar in all the locations, with maxima in August and minima in February.

Table 2: Monthly and yearly averaged daily values of Insolation and Sunshine Hours in the stations.

Mean

Max

Min

Mean

Max

Min

Station

Insolation

Insolation

Insolation

Sunshine

Sunshine

Sunshine

(kJ m-2)

(kJ m-2)

(kJ m-2)

Hours

Hours

Hours

Ferrol

12446

20307 Jun

4069 Dec

5.6

8.5 Jul

2.7 Dec

Lourizan

11898

20383 Jun

4194 Dec

5.7

9.1 Jul

2.8 Dec

P. Murias

10511

16628 May

3754 Dec

4.7

6.5 May

2.3 Jan

Alto

Rodicio

12687

20962 Jun

4444 Dec

5.9

8.8 Jul

3.1 Dec

Monte

Aloia

12545

20697 Jun

5013 Dec

6.0

8.5 Jul

3.6 Jan

Corrubedo

13842

22807 Jun

5185 Dec

6.6

9.6 Jun

3.5 Dec

Ancares

13805

23748 Jul

5264 Dec

6.0

9.7 Jul

3.3 Dec

Sergude

10629

18562Jul

3809 Dec

5.3

8.1 Jul

2.3 Dec

As it was already mentioned, mean precipitations are higher than the rest of Spain [10]. The highest values are found in Monte Aloia, due to the westerly air masses with high water content that

channel into the outfall of Mino river and ascend to the slope of the local mountain range, Sierra do Galineiro. In all the stations, maxima occur in October and minima in summertime.

Ferrol

Lourizan

Murias

Rodicio

Aloia

Corrubedo

Ancares

Sergude

14.2

14.8

13.8

10.2

12.7

14.8

8.4

13.5

23.7

Aug

26.3 Aug

22.5

Aug

22.7 Jul

24.1

Aug

23.3 Aug

20.1 Aug

25.9

Aug

7.2 Feb

5.8 Feb

6.4 Feb

2.3 Feb

5.4

Feb

8.6 Feb

0 Feb

4.3 Feb

82

77

88

84

79

82

79

84

85 Oct

83 Dec

92 Jun

91 Nov

87

Oct

88 Oct

87 Oct

90 Jan

79 Apr

73 Jun

84 Apr

77 Aug

75

Aug

72 Jun

72 Jun

78 Aug

1225

1574

918

1108

2097

1032

1406

1409

211Oct

331 Oct

131Oct

220 Oct

465

Oct

150 Oct

230 Oct

261 Oct

45 Aug

46 Jul

39 Jun

41 Jul

62 Jul

30 Jun

44 Aug

24 Aug

Table 3: Monthly and yearly averaged daily values of Temperature, Relative Humidity and accumulated

Precipitations in the stations analysed.

Clearness Index: averaged values and frequency distribution

image089

The monthly averaged daily data of KT are depicted in Fig. 4 for five stations. It can be noticed that the station of Ancares is the one with the highest values of the index all over the year; it can be explained because of its altitude, rural environment and relatively cloud-free atmosphere. These values have a maximum in summertime.

The lowest clearness indexes are found in Pedro Murias and Lourizan. The first station, located in a rural environment, is characterized by high values of relative humidity that produce, together with low values of temperature (Table 3) and persistent fogs, especially in summertime. The station of Lourizan, located in an urban environment, is characterized by the lowest values of relative

humidity (77% yearly average) and the highest value of temperature (table 3). These features suggest that low values of KT are produced by anthropogenic aerosols, generated by local factories and urban pollution.

The station of Ferrol, despite of its suburban location, close to a sea port, has values of KT higher than Lourizan, due to the presence of winds that clean up the atmosphere from aerosols and fogs.

3.1 Ground data vs. satellite observations

Yearly and monthly averaged daily values of global irradiation collected in the meteorological stations were compared with the averaged values derived from satellite images and collected in the Solar Atlas of Galicia [4].

Figure 5 shows, for the stations analyzed, the distribution of the yearly averaged daily values measured by pyranometers at five stations versus the estimated values by Vazquez et al. [4]. Global irradiances estimated from satellite images overestimates the data collected by the meteorological stations in every location, except in Ancares. In three cases (Ancares, Ferrol and Corrubedo) the agreement between ground and satellite is quite good, with relative differences less than 3% (Table 4). In the other sites, differences are higher, ranging from 10.7% (Alto do Rodicio) to 25.3% (Sergude). Relative differences in the monthly averaged daily values are greater, up till 76%. Mainly, overestimations by ground measurements occur during winter.

image090

Fig. 5: Distribution of the yearly mean values of global irradiation for the stations analyzed

The reason for these disagreements may be found in the normalization of the satellite data adopted by Vazquez et al. [4]. To obtain monthly averaged daily values, the authors apply only one coefficient for the entire region, instead of dividing the territory in areas characterized by similar climatologic features and calculating different normalization coefficients for each area. In their work, Vazquez et al. [4], compare their results only with one station, obtaining a very good agreement. For that, they assume the goodness of the results achieved for the entire region that is, as previously stated, characterized by complex topography and high climatologic variability.

Table 4: Relative differences (satellite-pyranometer) between the yearly and monthly averaged daily values

of global irradiation in the location analysed.

Month

Ferrol

Lourizan

Murias

Rodicio

Aloia

Corrubedo

Ancares

Sergude

Jan

-3.9

8.7

16.2

-2.4

3.7

-5.2

-23.6

24.7

Feb

-5.9

-7.7

9.2

-9.0

1.3

-17.7

-23.6

42.2

Mar

13.5

44.0

23.2

22.0

36.6

18.7

10.4

32.2

Apr

-3.06

19.5

12.7

13.5

13.4

1.0

1.8

36.4

May

0.8

19.9

12.6

13.6

11.8

-0.2

-1.4

23.4

Jun

2.8

13.0

24.2

15.1

13.1

5.8

-0.2

29.7

Jul

6.9

34.8

21.2

18.4

18.9

9.5

-4.5

16.4

Ago

-1.5

8.4

18.4

10.3

13.7

3.0

-8.5

12.2

Sep

2.9

-4.9

19.9

4.9

8.5

7.3

-1.8

21.9

Oct

6.4

3.6

75.9

-10.6

15.4

0.6

6.2

22.4

Nov

14.3

26.8

25.5

24.6

5.8

-16.1

1.1

19.3

Dec

6.2

11.6

15.1

1.3

0.5

-6.3

-17.9

18.1

Avg

2.7

15.0

18.2

10.7

13.4

2.7

-2.2

25.3

4. Conclusions

The evaluation of the solar resource in eight sites of Galicia has been carried out for the period 2001-2006. The analysis of global irradiation, sunshine hours, clearness index, together with other meteorological parameters — precipitations, relative humidity and temperature — allowed a characterization of the solar resource in this region.

Monthly averaged daily values of global radiation and sunshine hours point out the complex climatology of the Galician territory. In an area of less than 30.000 km2, distributed over 2° of latitude and longitude, differences in the yearly mean daily values of global irradiation of more than 3 MJm-2 per day were found between the locations considered. This range of variability in the values of global irradiation is comparable to that of Germany (Meteonorm; Bern, Switzerland).

In the same way, yearly averaged daily values of sunshine duration show differences of almost 2 hours per day between the sites.

The analysis of precipitations, temperature and relative humidity, combined with the study of solar radiation, evidence the presence of persistent fogs in particular zones of Galicia. The station of Pedro Murias represents the most remarkable example. In the same way, the analysis of the solar radiation combined with the mean values of precipitations points out that extremely rainy areas, as that represented by the station of Monte Aloia, are not necessarily associated with low insolation.

The complexity of the Galician climate comes also out analysing the mean values of the clearness index: while the highest values are found in high elevation sites, as a consequence of a clearer atmosphere, the lowest values occur in Lourizan and Pedro Murias due, respectively, to the urban pollution and to persistent fog episodes.

The comparison of solar irradiation ground measurements vs. satellite observations gives the opportunity for making several remarks. The solar irradiation maps obtained by Vazquez et al. [4] evidence two different trends: a) the incoming solar radiation increases with the decreasing latitude and b) coastal zones receive more radiation than the inner nearby areas.

The ground measurements of global irradiation recorded in eight sites of Galicia depict a much more complex pattern not so easy to generalize. The latitude effect seems to be inexistent or, at least, is masked by more important local climatologic factors. Thus, Corrubedo and Lourizan, being at the same latitude, show yearly averaged daily differences of 2 MJ m-2. The same

differences are found in Ferrol and Pedro Murias, that have similar geographical features. In these areas, Vazquez et al. found differences of 0.36 MJ m-2 per day.

The stations of Ferrol and Alto do Rodicio, located respectively in the coast and in the interior, have similar mean values of irradiation even they are separated by 1° in latitude. The highest levels of radiation are found in Corrubedo and Ancares, located at the same latitude, but respectively on the western and eastern edges of Galicia. All these facts drive to the conclusion that the distribution of the solar resource in a region such as Galicia cannot be explained without the support of other climatologic and topographic features.

The present work was intended for characterizing the solar resource in eight sites of Galicia over a period of 5 years. In the future, the data record will increase also due to the installation of 17 more first class pyranometers and 30 pyranometers with Photovoltaic sensor. This will allow a more detailed characterization of the solar resource in time and space.

Due to the complex topography and climatology of Galicia, this will not be enough to obtain solar maps of the region since, in these cases, interpolation techniques do not provide sufficiently reliable results [11]. However, data from this solar radiation network will represent a very important tool to validate and calibrate the methods to estimate the solar resource.

Acknowledgements

Meteorological dataset provided by MeteoGalicia (Xunta de Galicia) from its web page is acknowledged This work was partially funded by Galician R&D Programme under project 07REM02CT.

References

[1] Institute Enerxetico de Galicia. Enema solar fotovoltaica na comunidade autonoma de Galicia. Conselleria de Inovacion, Industria e Comercio. Xunta de Galicia; 2003.

[2] Font Tullot I. Atlas de la radiacion solar en Espana. Instituto Nacional de Meteorologia, Ministerio de Transportes, Turismo y Comunicaciones. Madrid, Spain; 1984.

[3] Vera Mella N. Atlas climatico de irradiacion solar a partir de imagenes del satelite NOAA. Aplicacion a la peninsula iberica. PhD thesis. Univ. Politecnica de Catalunya, Barcelona; 2005.

[4] Vazquez Vazquez M., Santos Navarro J. M., Prado Cerqueira M. T., Vazquez Rios D., Rodrigues Fernandes F. M. Atlas de radiacion solar de Galicia. Universidad de Vigo. Vigo, Espana; 2005.

[5] Rigollier, C., Lefevre, M. and Wald, L. The method Heliosat-2 for deriving shortwave solar radiation from satellite images. Solar Energy 2004; 77: 159-169.

[6] Salson S., Souto J. A. Automatic weather stations network of the department of environment of Galicia: data acquisition, validation and quality control, Proceedings of the 3rd international conference on experiences with automatic weather stations, Torremolinos, Spain; 2003.

[7] Pettazzi A., Souto J. A., Salson S. EOAS, a shared joint atmospheric observation site of MeteoGalicia. Proceedings of 4th ICEAWS — International Conference on Experiences with Automatic Weather Stations, Lisbon, Portugal; 2006.

[8] Davies J. A. Validation of models for estimating solar radiation on horizontal surfaces. Report available from the IEA, Downsview, Ontario, Canada; 1988.

[9] Iqbal M. An introduction to solar radiation, Academic Press, San Diego, CA; 1983.

[10] Instituto Nacional de Meteorologia. Guia resumida del clima en Espana 1971-2000. Instituto Nacional de Meteorologia, D25.3, Ministerio de Medio Ambiente. Madrid, Espana; 2001.

[11] Batlles F. J., Martinez-Durban M., Miralles I., Ortega R., Barbero F. J., Tovar-Pescador J., Pozo — Vazquez D., Lopez G. Evaluacion de los recursos energeticos solares en zonas de topografia compleja. XII Congreso Iberico y VII Congreso Ibero Americano de Energia Solar. Vigo, Espana; 2004.

Mapping Solar Radiation over Complex Topography Areas Combining Digital Elevation Models and Satellite Images

J. L. Bosch1*, L. F. Zarzalejo2, F. J. Batlles1 and G. Lopez3

1 Universidad de Almeria, Departmento de Fisica Aplicada, Ctra. Sacramento s/n, 04120-Almeria, Espana
2 CIEMAT, Departamento de Energia, Madrid, Espana
3 EPS-Universidad de Huelva, Departamento de Ingenieria Electrica y Termica, Huelva, Espana

Corresponding Author, jlbosch@ual. es

Abstract

The correlation of solar irradiation data in flat and homogeneous areas is relatively high and classic interpolation methods are very suitable for its estimation. However, in complex topography zones, a simple data interpolation is not adequate. On the other hand, spatial variability of solar irradiation is also affected by site latitude and cloud cover distribution. In this work, a methodology has been implemented consisting in daily solar irradiation estimation for all sky conditions, by means of Meteosat satellite images and additional information from a Digital Terrain Model (DTM) of the studied area. Solar irradiation is calculated following the HELIOSAT-2 methodology; and a method is presented to obtain the horizon of the studied points using the DTM. The effect of the snow covers is also studied. Model performance has been evaluated against data measured in 14 radiometric stations located in a mountainous area, offering good results, with a Root Mean Square Error (RMSE) around 11%. Finally, a daily solar irradiation map has been generated for the complex topography site.

Keywords: Daily Irradiation Mapping, DTM, Meteosat, Complex Terrain

1. Introduction

Incoming solar radiation, through its influence over the energy and water balances of the earth surface, affects processes like air and soil heating, photosynthesis, wind or snow thawing. Therefore, its knowledge is important in diverse fields and necessary for several applications. For most of these applications, global radiation measures are needed over wide regions, for long time periods and with a high spatial resolution.

At local scales the topography is the most important factor in the distribution of the solar radiation on the surface. In plane and homogeneous areas, classic interpolation methods can estimate the solar radiation accurately. However, in zones with a high topographical variability the spatial correlation is difficult to detect, and for distances between 300 and 1000 m is very small or disappear [1]. The use of interpolation in this kind of terrains can lead to large errors and more complex models that include topographical information are needed [2]. In recent years, digital models of terrain (DMT) have been utilized to develop radiation models incorporating topography and consequently, the spatial variability
of terrain mentioned before. In addition, spatial variability of solar radiation is affected by latitude and cloud cover distribution. This issue has been studied recently with the aid of satellite images. The geostationary satellites (METEOSAT) permanently occupy the same zone over the earth surface, acquiring several images per day, for this reason they are suitable for estimating the solar irradiation and evaluating the energy potential of a wide area. The Centre Energetique et Precedes (CEP), the School of Mines of Paris, in cooperation with other European Centers of Investigation, developed a statistical model to estimate the solar irradiation in the terrestrial surface from images of METEOSAT. The mentioned model is known as HELIOSAT [3]. The basic idea of the HELIOSAT is the interrelation between the cloud cover and the global incident irradiation on a point on the earth’s surface. This model was one of earliest used for the evaluation of the global irradiation from images of satellite. It was developed using measurements of French stations and its goal was the estimation of average monthly values of global irradiation. Afterwards different modifications were introduced [4] until a new version named HELIOSAT-2 was implemented.

In this work, a methodology is presented and tested, in which the estimation of solar irradiation is performed by using a modified HELIOSAT-2 [5], together with the information contained in a DTM. Instead of a single irradiation value for a pixel, the horizon of each inner point is used to estimate around 1000 irradiation values for every pixel. Computational cost of the horizon calculation can be a problem when dealing with large areas; that issue has been addressed by developing an algorithm that reduces drastically the time utilized in this process, without loosing much information about the actual horizon. Additionally, the happening of snow covers can lead to a subestimation of the model, because those pixels can be considered covered by very bright clouds instead of snow, this problem has been also addressed in this work with satisfactory results.

The main goal is to perform an irradiation map of daily values from the satellite images, fitting the spatial resolution of pixels (~ 3.5 km) to the resolution of the DTM (100 m). Ground measurements registered at 14 stations located in a complex topography area have been used for validation purposes, observing an error reduction after the consideration of the horizon and snow effects. It is also interesting to note that this procedure can be applied under all kind of sky conditions.

User guidance

1.1 User survey

One of main objectives of the MESoR project is the involvement of key stakeholders throughout the project. Their first task was the participation in a user survey which was conducted as an interview via telephone calls or e-mail. The survey performed a comparative analysis of various solar radiation platforms concerning technical aspects but also addressing usability, integration and pre-commercial information. This analysis based on the evaluations expressed by actual “top-users” and “top- customers” of the platforms. The sample has been composed by current users of the services belonging to various academic, scientific, industrial and business categories, from public and private sectors. The organisations are active in the fields of architecture/building, PV and other solar applications. The initial sample of 53 was selected by each partner according to the criteria of importance, frequency of usage and attitude to scientific cooperation.

The collected answers indicate a very high degree of awareness about the analysed issues. This is witnessed by the fact that most of the respondents use multiple services, they have a deep knowledge of each service, are able to compare the various services and to highlight the related points of strength and weakness.

The survey detected a gap between expectations and the satisfaction as for quality and accuracy of some parameters, reliability of data measurement and calculation, comparability of data across the services and personalisation of services.

The users expect a truly new and integrated service that offers standardised data and protocols.

Specifications of the solar radiation databases and underlying methods

1.1. General specifications

Spatially-distributed (map) solar radiation databases are classified according to several factors:

• Input data from which they have been created: (a) observations from the meteorological stations (global, diffuse and direct irradiances, and other relevant climate data), (b) digital satellite images or (c) combination of both; here also ancillary atmospheric data used in the models are considered, such as water vapour, ozone and aerosols;

• Period of time (typically a number of years) which is represented by the input data;

• Spatial resolution, i. e geographical distribution of the meteorological sites, grid resolution of the satellite data and resulting outputs;

• Time resolution, which characterises periodicity of the measurement of the input data and of the resulting parameters. Thus a primary database may include time series with periodicity of a few minutes up to hourly and daily averages (sums), or it may contain only monthly and long-term averages.

• Methodical approach used for computation of the primary database: typically solar radiation models combined with interpolation methods (e. g. geostatistical methods or splines, in case that ground observations are used) or algorithms for satellite data processing (e. g. Heliosat). Primary database typically consists of global, direct normal or diffuse irradiances (irradiation in case of time-integrated products) and also some auxiliary parameters such as clear-sky index.

• Simulation models used for calculation of derived parameters, such as global irradiance for inclined and sun-tracking surface, spectral products (e. g., illuminances, UV and PV-related irradiances), estimation of terrain effects, derived statistical products (e. g. synthetic time series).

Quality of an individual data set is assessed for a set of locations by comparing them to ground measurements, where the first order statistics is calculated (bias, root mean square deviation, standard deviation, the correlation coefficient) and the frequency distribution is analysed. In this work we focus on the relative map-based cross comparison of several solar radiation products.

Such comparison provides means for improved understanding of regional distribution of the uncertainty by combining all existing resources (calculating the average of all) and quantifying their mutual agreement by the means of standard deviation.