# Improving HelioClim-3 estimates of surface solar irradiance using the McClear clear-sky model and recent advances in atmosphere composition

Zhipeng Qu, Benoît Gschwind, Mireille Lefèvre, Lucien Wald

► To cite this version:

Zhipeng Qu, Benoît Gschwind, Mireille Lefèvre, Lucien Wald. Improving HelioClim-3 estimates of surface solar irradiance using the McClear clear-sky model and recent advances in atmosphere composition. Atmospheric Measurement Techniques, 2014, 7, pp.3927 - 3933. 10.5194/amt-7-3927-2014. hal-01087291

HAL Id: hal-01087291

https://minesparis-psl.hal.science/hal-01087291

Submitted on 25 Nov 2014

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
<figure><img src="image_2.png" /></figure>

<figure><img src="image_3.png" /></figure>

# Improving HelioClim-3 estimates of surface solar irradiance using the McClear clear-sky model and recent advances in atmosphere composition

Z. Qu, B. Gschwind, M. Lefevre, and L. Wald

MINES ParisTech, PSL Research University, O.I.E. Centre for Observation, Impacts, Energy, CS 10207, rue Claude Daunesse, 06904 Sophia Antipolis CEDEX, France

Correspondence to: L. Wald (lucien.wald@mines-paristech.fr)

Received: 8 June 2014 – Published in Atmos. Meas. Tech. Discuss.: 29 July 2014

Revised: 7 October 2014 – Accepted: 23 October 2014 – Published: 25 November 2014

**Abstract.** The HelioClim-3 database (HC3v3) provides records of surface solar irradiation every 15 min, estimated by processing images from the geostationary meteorological Meteosat satellites using climatological data sets of the atmospheric Linke turbidity factor. This technical note proposes a method to improve a posteriori HC3v3 by combining it with data records of the irradiation under clear skies from the new McClear clear-sky model, whose inputs are the advanced global aerosol property forecasts and physically consistent total column content in water vapour and ozone produced by the MACC (Monitoring Atmosphere Composition and Climate) projects. The method is validated by comparison with a series of ground measurements for 15 min and 1 h for 6 stations and for daily irradiation for 23 stations. The correlation coefficient is large, greater than respectively 0.92, 0.94, and 0.97, for 15 min, 1 h and daily irradiation. The bias ranges from -4 to 4% of the mean observed irradiation for most sites. The relative root mean square difference (RMSD) varies between 14 and 38% for 15 min, 12 and 33% for 1 h irradiation, and 6 and 20% for daily irradiation. As a rule of thumb, the farther from the nadir of the Meteosat satellite located at latitude 0° and longitude 0°, and the greater the occurrence of fragmented cloud cover, the greater the relative RMSD. The method improves HC3v3 in most cases, and with no degradation in the others. A systematic correction of HC3v3 with McClear is recommended.

## 1 Introduction

The downwelling solar irradiance observed at ground level on horizontal surfaces and integrated over the whole spectrum (total irradiance) is called surface solar irradiance (SSI). It is the sum of the direct irradiance from the direction of the sun and the diffuse from the rest of the sky vault, and is also called the global irradiance. The SSI is an essential climate variable (ECV) as established by the Global Climate Observing System in August 2010 (GCOS, 2013). Knowledge of the SSI and its geographical distribution is of prime importance for numerous domains where SSI plays a major role, as e.g. in weather, climate, biomass, and energy.

The HelioClim project is an ambitious initiative of MINES ParisTech launched in 1997 to increase knowledge about the SSI and to offer SSI values for any site, any instant over a large geographical area, and a long period of time, to a wide audience (Blanc et al., 2011). The project comprises several databases that cover Europe, Africa and the Atlantic Ocean. The HelioClim-1 (HC1) database offers daily means of the global SSI for the period 1985–2005. The HelioClim-3 (HC3) database contains 15 min values of the global SSI. It was created in 2004, and is updated daily from images taken by the Meteosat Second Generation satellites. Its recent improvements have taken place in the framework of the European MACC and MACC-II (Monitoring Atmosphere Composition and Climate) projects funded by the European Commission. The HelioClim-4 database is under creation in these MACC projects. It will contain 15 min values of the global, direct and diffuse components of the SSI with a daily update.
<table><caption>**Table 1.** List of stations, ordered from north to south. Data from 2005 to 2009, except if stated otherwise.</caption><thead><tr><th>Station</th><th>Country</th><th>Latitude</th><th>Longitude</th><th>Elevation (m)</th><th>Available summarization</th></tr></thead><tbody><tr><td>Toravere</td><td>Estonia</td><td>58.250</td><td>26.467</td><td>70</td><td>15 min, 1 h, 1 day</td></tr><tr><td>Rucana</td><td>Latvia</td><td>56.150</td><td>21.167</td><td>18</td><td>1 day</td></tr><tr><td>Hamburg</td><td>Germany</td><td>53.633</td><td>10.000</td><td>16</td><td>1 day</td></tr><tr><td>Valentia</td><td>Ireland</td><td>51.933</td><td>-10.250</td><td>9</td><td>1 day</td></tr><tr><td>Uccle</td><td>Belgium</td><td>50.800</td><td>4.350</td><td>100</td><td>1 day</td></tr><tr><td>Camborne</td><td>UK (2004–2007)</td><td>50.217</td><td>-5.317</td><td>88</td><td>15 min, 1 h, 1 day</td></tr><tr><td>Vienna</td><td>Austria</td><td>48.250</td><td>16.367</td><td>203</td><td>1 day</td></tr><tr><td>Kishinev</td><td>Moldova (2005–2007)</td><td>47.000</td><td>28.817</td><td>205</td><td>1 day</td></tr><tr><td>Payerne</td><td>Switzerland</td><td>46.815</td><td>6.944</td><td>491</td><td>15 min, 1 h, 1 day</td></tr><tr><td>Carpentras</td><td>France</td><td>44.083</td><td>5.059</td><td>100</td><td>15 min, 1 h, 1 day</td></tr><tr><td>Nice</td><td>France</td><td>43.650</td><td>7.200</td><td>4</td><td>1 day</td></tr><tr><td>Thessaloniki</td><td>Greece (2005–2006)</td><td>40.633</td><td>22.967</td><td>60</td><td>1 day</td></tr><tr><td>Casablanca</td><td>Morocco (2005)</td><td>33.567</td><td>-7.667</td><td>62</td><td>1 day</td></tr><tr><td>Mersa Matruh</td><td>Egypt</td><td>31.333</td><td>27.217</td><td>25</td><td>1 day</td></tr><tr><td>El Arish</td><td>Egypt</td><td>31.083</td><td>33.750</td><td>31</td><td>1 day</td></tr><tr><td>Sede Boqer</td><td>Israel</td><td>30.905</td><td>34.782</td><td>500</td><td>15 min, 1 h, 1 day</td></tr><tr><td>Asyut</td><td>Egypt (2005–2007)</td><td>27.200</td><td>31.167</td><td>52</td><td>1 day</td></tr><tr><td>Aswan</td><td>Egypt</td><td>23.967</td><td>32.783</td><td>192</td><td>1 day</td></tr><tr><td>Tamanrasset</td><td>Algeria (2005–2007)</td><td>22.783</td><td>5.517</td><td>1378</td><td>15 min, 1 h, 1 day</td></tr><tr><td>Rochambeau</td><td>French Guiana</td><td>4.822</td><td>-52.365</td><td>4</td><td>1 day</td></tr><tr><td>Brasilia</td><td>Brazil (2006–2007)</td><td>-15.601</td><td>-47.713</td><td>1023</td><td>1 day</td></tr><tr><td>Bulawayo</td><td>Zimbabwe (2005)</td><td>-20.150</td><td>28.620</td><td>1343</td><td>1 day</td></tr><tr><td>Maputo</td><td>Mozambique</td><td>-25.967</td><td>32.600</td><td>70</td><td>1 day</td></tr></tbody></table>

<table><caption>**Table 2.** Comparison of correlation coefficients for 15 min irradiation and clearness index. Best values are in bold.</caption><thead><tr><th>Station</th><th>I<sub>HC3v3</sub></th><th>I<sub>HC3McClear</sub></th><th>KT<sub>HC3v3</sub></th><th>KT<sub>McClear</sub></th></tr></thead><tbody><tr><td>Toravere</td><td>0.921</td><td><strong>0.924</strong></td><td>0.765</td><td><strong>0.773</strong></td></tr><tr><td>Camborne</td><td>0.950</td><td><strong>0.952</strong></td><td>0.829</td><td><strong>0.830</strong></td></tr><tr><td>Payerne</td><td>0.958</td><td><strong>0.959</strong></td><td>0.846</td><td><strong>0.853</strong></td></tr><tr><td>Carpentras</td><td>0.970</td><td><strong>0.972</strong></td><td>0.842</td><td><strong>0.845</strong></td></tr><tr><td>Sede Boqer</td><td>0.973</td><td><strong>0.975</strong></td><td>0.829</td><td><strong>0.846</strong></td></tr><tr><td>Tamanrasset</td><td>0.965</td><td><strong>0.970</strong></td><td>0.824</td><td><strong>0.830</strong></td></tr></tbody></table>

The HelioClim databases are available on the Internet through the SoDa website (http://www.soda-pro.com), and support research and business by providing data of known quality on surface solar irradiance. More than 100 000 requests were made in 2012 to HC1 by users and more than 2 million to HC3, demonstrating the large use of HelioClim databases. Lefèvre et al. (2014) performed a review of the scientific literature citing HelioClim, and found many examples of usages in various domains: oceanography, climate, energy production, life cycle analysis, agriculture, ecology, human health, and air quality.

The HC1 and HC3 databases are derived from images of the Meteosat series of satellites using the Heliosat-2 method (Rigollier et al., 2004). The Heliosat-2 method needs a so-called clear-sky model to predict the SSI that should be observed under a clear sky. The European Solar Radiation Atlas

(ESRA) clear-sky model (Rigollier et al., 2000) modified by Geiger et al. (2002) was selected, with the climatology of the Linke turbidity factor from Remund et al. (2003) as input. The Linke turbidity factor is a convenient approximation for modelling the atmospheric absorption and scattering of the solar radiation under clear skies. The climatology of Remund et al. (2003) comprises 12 maps, one per month, covering the world with cells of 5' of arc angle in size. The use of this climatology is one of the drawbacks of the HC1 and HC3, and especially HC3, whose high temporal resolution (15 min) is in principle well suited for monitoring and reproducing rapid changes in SSI. Aerosols have different scattering and absorbing properties according to their type and the spatial and temporal heterogeneity of their number, size, chemical composition, and shape (Elias and Roujean, 2008; Xu et al., 2011). These properties as well as total column content in water vapour and ozone may vary rapidly within a day or from day to day, thus influencing the SSI under clear skies. Climatology cannot account for such changes, and HC3 estimates are often underestimated in the case of clear skies (Lefèvre et al., 2013). In addition, the Linke turbidity factor has a drawback inherent to its definition. It is a single value that summarises the effects of many variables. Simultaneous changes in these variables induce changes in irradiation under clear-sky conditions that may not be reflected in the Linke turbidity factor and therefore not in the irradiation estimated by the ESRA model.
The MACC and MACC-II projects are preparing the operational provision of global aerosol property analyses and forecasts together with physically consistent total column content in water vapour and ozone available every 3 h (Benedetti et al., 2011; Kaiser et al., 2012; Peuch et al., 2009). Up to now, a multi-annual reanalysis data set has been provided, and is used here (Inness et al., 2013). Such information has not been available so far from any operational numerical weather prediction (NWP) centre. A new clear-sky model called McClear has been developed to exploit this new input data source for estimating the direct and global SSI (Lefèvre et al., 2013). Validation of McClear outputs against beam and global irradiances measured at 1 min by BSRN stations in the world reveals satisfactory results. Good correlation is attained; bias, standard deviation and root mean square error (RMSE) are small (Lefèvre et al., 2013).

How can such advanced data sets on aerosol properties, water vapour and ozone be exploited to bring a significant improvement to the widely used HC3 without re-factoring the Heliosat-2 method and re-processing all Meteosat images since 2004? If this is possible, the dynamics of the aerosol properties, water vapour and ozone would be taken into account in the enhanced HC3, thus possibly yielding better estimates under clear-sky conditions. This technical note investigates the changes brought to HC3 in an a posteriori manner, i.e. by applying post-processing to the HC3 estimates, and assesses the benefit compared to the original HC3.

## 2 Data sets and method

The method is the following. A standard request to HC3v3 (version 3 of HC3) for a given site integrated over a given period, called summarization, e.g. 1 h or 1 day, yields several data, including the global SSI $I_{HC3v3}$, that under clear-sky condition $I_{ESRA}$, and $I_0$ the irradiance received on a horizontal surface at the top of atmosphere. The clear-sky index $K_c$ is computed:

$$K_c = I_{HC3v3} / I_{ESRA}. \qquad (1)$$

The McClear model may be invoked through the SoDa website. It yields the clear-sky value $I_{McClear}$ for the requested summarization and site, and the new version of the SSI $I_{HC3McClear}$ is obtained:

$$I_{HC3McClear} = K_c I_{McClear}. \qquad (2)$$

A series of ground measurements of surface solar irradiation $I_{ground}$ was assembled and serves as a reference in the comparison of $I_{HC3v3}$ and $I_{HC3McClear}$. Comparison was performed for the period 2005–2009. Measurements were collected from 23 stations located in the field of view of the Meteosat satellite.

Measurements of 15 min, hourly, and daily irradiation were collected at six stations of the BSRN (Baseline Surface

Radiation Network). BSRN stations record global irradiation $I_{ground}$ as well as its direct B and diffuse D components every minute. Roesch et al. (2011) recommend keeping only $I_{ground}$ measurements that obey the following constraints:

$$\text{if } \theta_S \le 75^{\circ}, 1.08 \ge (D+B)/I_{\text{ground}} \ge 0.92 \qquad (3)$$

$$\text{if } \theta_S > 75^{\circ}, 1.15 \ge (D+B)/I_{\text{ground}} \ge 0.85, \qquad (4)$$

where $\theta_S$ is the solar zenith angle. Roesch et al. (2011) note a percentage of missing data of 4% for global irradiance and 13% for direct irradiance for all studied BSRN data. The original measurements passing Eq. (3) are summed up to yield 15 min, hourly and daily irradiation, provided a sum is made with at least 90% valid measurements, e.g. 54 valid 1 min data to yield 1 h irradiation.

Measurements of daily irradiation were collected for another set of 17 stations of the meteorological networks archived in the WRDC (World Radiation Data Center) and available through the WRDC website. Table 1 lists the 23 stations that have been selected in order to represent various climates in Europe, Africa and tropical South America.

The clearness index KT, also called the global transmissivity of the atmosphere, or atmospheric transmittance, or atmospheric transmission, is defined as

$$KT = I_{\text{ground}} / I_0. \qquad (5)$$

For clear skies, KT is close to 0.8, and is close to 0 for overcast skies. This index has the advantages of removing most of the effects due to sun position and indicating the type of sky. While irradiation for clear-sky conditions but a large solar zenith angle may be similar to that under cloudy conditions but with a low solar zenith angle, the clearness indices will be different. The clearness index is useful for analysing causes of discrepancies between the data sets. The clearness indices $KT_{HC3v3}$ and $KT_{McClear}$ are computed:

$$KT_{HC3v3} = I_{HC3v3} / I_0 \qquad (6)$$

$$KT_{McClear} = I_{HC3McClear} / I_0. \qquad (7)$$

The irradiation $I_{HC3McClear}$, and hence the clearness index $KT_{McClear}$, are computed for each summarization: 15 min, 1 h, and 1 day:

$$(I_{HC3McClear})_{\text{hour}} = [(I_{HC3v3})_{\text{hour}}(I_{HC3McClear})_{\text{hour}}] / (I_{ESRA})_{\text{hour}} \qquad (8)$$

$$(I_{HC3McClear})_{\text{day}} = [(I_{HC3v3})_{\text{day}}(I_{HC3McClear})_{\text{day}}] / (I_{ESRA})_{\text{day}}, \qquad (9)$$

where the quantities $(I_{HC3McClear})_{\text{hour}}$, $(I_{HC3v3})_{\text{hour}}$, $(I_{HC3McClear})_{\text{hour}}$, $(I_{ESRA})_{\text{hour}}$, $(I_{HC3McClear})_{\text{day}}$, $(I_{HC3v3})_{\text{day}}$, $(I_{HC3McClear})_{\text{day}}$, and $(I_{ESRA})_{\text{day}}$ are directly retrieved from the SoDa website. Another approach could be to compute $I_{HC3McClear}$ every 15 min, and then to perform the summarization for 1 h or 1 day, though this is less practical for the many users of the SoDa website.
<table><caption>**Table 3.** Comparison of differences for 15 min irradiation, in J cm⁻². The mean value is obtained from the measurements. The first value is *I*<sub>HC3v3</sub> and the second is *I*<sub>HC3McCle</sub>, with the best value in bold. Bias and RMSD of *I*<sub>HC3McCle</sub> relative to the mean irradiation are given in brackets.</caption><thead><tr><th>Station</th><th>Mean 15 min irradiation</th><th>Bias</th><th>Standard deviation</th><th>RMSD</th></tr></thead><tbody><tr><td>Toravere</td><td>20.5</td><td>0.3 / <strong>-0.1</strong> (0%)</td><td>7.9 / <strong>7.8</strong></td><td>7.9 / <strong>7.8</strong> (38%)</td></tr><tr><td>Camborne</td><td>23.2</td><td><strong>0.4</strong> / 0.5 (2%)</td><td>6.8 / 6.8</td><td>6.8 / 6.8 (29%)</td></tr><tr><td>Payerne</td><td>25.5</td><td>-1.6 / <strong>-0.3</strong> (-1%)</td><td>6.8 / 6.8</td><td>7.0 / <strong>6.8</strong> (27%)</td></tr><tr><td>Carpentras</td><td>31.9</td><td>0.6 / 0.6 (2%)</td><td>6.3 / <strong>6.1</strong></td><td>6.4 / <strong>6.2</strong> (19%)</td></tr><tr><td>Sede Boqer</td><td>47.6</td><td>-3.0 / <strong>-1.8</strong> (-4%)</td><td>6.5 / <strong>6.4</strong></td><td>7.1 / <strong>6.6</strong> (14%)</td></tr><tr><td>Tamanrasset</td><td>47.6</td><td><strong>1.3</strong> / 1.7 (4%)</td><td>8.2 / <strong>7.7</strong></td><td>8.3 / <strong>7.8</strong> (16%)</td></tr></tbody></table>

<table><caption>**Table 4.** Comparison of differences for hourly irradiation, in J cm⁻². The mean value is obtained from the measurements. The first value is *I*<sub>HC3v3</sub> and the second is *I*<sub>HC3McCle</sub>, with the best value in bold. Bias and RMSD of *I*<sub>HC3McCle</sub> relative to the mean irradiation are given in brackets.</caption><thead><tr><th>Station</th><th>Mean hourly irradiation</th><th>Bias</th><th>Standard deviation</th><th>RMSD</th><th>Correlation coefficient</th></tr></thead><tbody><tr><td>Toravere</td><td>76.8</td><td>1.1 / <strong>-0.4</strong> (-1%)</td><td>25.9 / <strong>25.2</strong></td><td>25.9 / <strong>25.2</strong> (33%)</td><td>0.945 / <strong>0.949</strong></td></tr><tr><td>Camborne</td><td>87.4</td><td><strong>1.6</strong> / 2.0 (2%)</td><td>21.5 / <strong>21.2</strong></td><td>21.5 / <strong>21.3</strong> (24%)</td><td>0.968 / <strong>0.970</strong></td></tr><tr><td>Payerne</td><td>95.0</td><td>-5.9 / <strong>-1.2</strong> (-1%)</td><td>20.9 / <strong>20.8</strong></td><td>21.8 / <strong>20.8</strong> (22%)</td><td>0.974 / <strong>0.976</strong></td></tr><tr><td>Carpentras</td><td>120.4</td><td>2.3 / <strong>2.2</strong> (2%)</td><td>20.6 / <strong>19.7</strong></td><td>20.7 / <strong>19.8</strong> (16%)</td><td>0.980 / <strong>0.982</strong></td></tr><tr><td>Sede Boqer</td><td>184.3</td><td>-11.1 / <strong>-6.8</strong> (-4%)</td><td>20.8 / <strong>20.1</strong></td><td>23.5 / <strong>21.2</strong> (12%)</td><td>0.984 / <strong>0.985</strong></td></tr><tr><td>Tamanrasset</td><td>179.8</td><td><strong>4.9</strong> / 6.5 (4%)</td><td>27.1 / <strong>24.4</strong></td><td>27.5 / <strong>25.3</strong> (14%)</td><td>0.977 / <strong>0.982</strong></td></tr></tbody></table>

For each summarization, the deviations (*I*<sub>HC3v3</sub> − *I*<sub>ground</sub>), (*I*<sub>HC3McCle</sub> − *I*<sub>ground</sub>), (*K*<sub>HC3v3</sub> − *K*<sub>ground</sub>) and (*K*<sub>HC3McCle</sub> − *K*<sub>ground</sub>) are computed and synthesised by the bias, the standard deviation, the root mean square difference (RMSD), and the correlation coefficient. Tables 2 to 5 report the results of the comparison for 15 min, hourly and daily irradiation respectively.

### 3 Results and discussion

The correlation coefficient for 15 min irradiation is reported in Table 2. For both *I*<sub>HC3v3</sub> and *I*<sub>HC3McCle</sub>, the correlation coefficient is very large, greater than 0.95, except for Toravere (0.91), indicating that the 15 min irradiation is well reproduced by both estimates. The correlation is slightly greater for *I*<sub>HC3McCle</sub> than for *I*<sub>HC3v3</sub>, showing that the combination of *I*<sub>HC3v3</sub> with McClear brings a better reproduction of the observed changes in irradiation. This observation is supported by the fact that the correlation coefficient for the clearness index *K*<sub>McCle</sub> is greater than that for *K*<sub>HC3v3</sub> (Table 2).

The bias for *I*<sub>HC3v3</sub> ranges from −3.0 to 1.3 J cm⁻² (Table 3). The bias for *I*<sub>HC3McCle</sub> is similar to or smaller than that for *I*<sub>HC3v3</sub> for all cases. An exception to the overall decrease in bias is Tamanrasset, where the bias is 1.3 J cm⁻² (3% of the mean irradiation) for *I*<sub>HC3v3</sub> and 1.7 J cm⁻² (4%) for *I*<sub>HC3McCle</sub>. A closer examination of the data sets of

irradiation and the clearness index for Tamanrasset reveals that *I*<sub>HC3v3</sub> exhibits a negative bias for clear-sky conditions and a positive bias for cloudy situations. The balance between these negative and positive biases yields an overall bias of 1.3 J cm⁻². The combination of *I*<sub>HC3v3</sub> with McClear yields more accurate results for clear-sky conditions, as expected. The bias in these conditions is now strongly reduced and close to 0. On the contrary, there is almost no change in the results for cloudy situations, which exhibit a positive bias. Contrary to *I*<sub>HC3v3</sub>, this positive bias is not counterbalanced in *I*<sub>HC3McCle</sub> by an equivalent but negative bias for clear skies. The result is that the bias in *I*<sub>HC3McCle</sub> is slightly greater than that of *I*<sub>HC3v3</sub>.

The standard deviation is fairly similar for *I*<sub>HC3v3</sub> and *I*<sub>HC3McCle</sub> for all stations. It ranges from 6.3 to 8.2 J cm⁻² for *I*<sub>HC3v3</sub>. It is smaller for *I*<sub>HC3McCle</sub>, and ranges from 6.1 to 7.8 J cm⁻². The smaller standard deviation may be linked to the better correlation coefficient observed for *I*<sub>HC3McCle</sub>. Similarly, the RMSD is slightly less for *I*<sub>HC3McCle</sub> than for *I*<sub>HC3v3</sub>. Tables 2 and 3 show that the combination of HC3 and McClear brings a benefit for 15 min irradiation for the six studied stations.

Table 4 reports results for hourly irradiation. The correlation coefficient for *I*<sub>HC3McCle</sub> is large, greater than 0.97, except for Toravere (0.95), and is greater than that for *I*<sub>HC3v3</sub>. The bias ranges from −11.1 to 4.9 J cm⁻² for *I*<sub>HC3v3</sub>, and from −6.8 to 6.5 J cm⁻² for *I*<sub>HC3McCle</sub>. The bias for
<table><caption>**Table 5.** Comparison of differences for daily irradiation, in J cm⁻². The mean value is obtained from the measurements. The first value is *I*<sub>HC3v3</sub> and the second is *I*<sub>HC3McCle</sub>, with the best value in bold. Bias and RMSD of *I*<sub>HC3McCle</sub> relative to the mean irradiation are given in brackets.</caption><thead><tr><th>Station</th><th>Mean daily irradiation</th><th>Bias</th><th>Standard deviation</th><th>RMSD</th><th>Correlation coefficient</th></tr></thead><tbody><tr><td>Toravere</td><td>1237</td><td>30 / 8 (1%)</td><td>204 / 184</td><td>206 / 184 (15%)</td><td>0.969 / 0.974</td></tr><tr><td>Rucana</td><td>1336</td><td>91 / -11 (-1%)</td><td>237 / 211</td><td>254 / 211 (16%)</td><td>0.966 / 0.971</td></tr><tr><td>Hamburg</td><td>1112</td><td>-26 / 6 (1%)</td><td>114 / 110</td><td>117 / 111 (10%)</td><td>0.989 / 0.991</td></tr><tr><td>Valentia</td><td>1065</td><td>42 / 49 (5%)</td><td>200 / 188</td><td>205 / 194 (18%)</td><td>0.968 / 0.972</td></tr><tr><td>Uccle</td><td>1113</td><td>-23 / 20 (2%)</td><td>108 / 110</td><td>111 / 112 (10%)</td><td>0.990 / 0.991</td></tr><tr><td>Camborne</td><td>1150</td><td>24 / 38 (3%)</td><td>169 / 156</td><td>171 / 160 (14%)</td><td>0.978 / 0.982</td></tr><tr><td>Vienna</td><td>1237</td><td>-36 / 0 (0%)</td><td>119 / 124</td><td>124 / 124 (10%)</td><td>0.989 / 0.989</td></tr><tr><td>Kishinev</td><td>1348</td><td>37 / 21 (2%)</td><td>171 / 154</td><td>175 / 156 (12%)</td><td>0.980 / 0.984</td></tr><tr><td>Payerne</td><td>1275</td><td>-79 / 5 (0%)</td><td>145 / 138</td><td>165 / 138 (11%)</td><td>0.985 / 0.987</td></tr><tr><td>Carpentras</td><td>1552</td><td>28 / 27 (2%)</td><td>162 / 126</td><td>164 / 129 (8%)</td><td>0.982 / 0.989</td></tr><tr><td>Nice</td><td>1589</td><td>48 / 60 (4%)</td><td>152 / 130</td><td>160 / 143 (9%)</td><td>0.984 / 0.988</td></tr><tr><td>Thessaloniki</td><td>1646</td><td>-74 / 3 (0%)</td><td>128 / 102</td><td>148 / 102 (6%)</td><td>0.988 / 0.992</td></tr><tr><td>Casablanca</td><td>1954</td><td>-45 / -44 (-2%)</td><td>176 / 172</td><td>181 / 178 (9%)</td><td>0.972 / 0.974</td></tr><tr><td>Mersa Matruh</td><td>1853</td><td>-36 / 52 (3%)</td><td>185 / 165</td><td>189 / 173 (9%)</td><td>0.962 / 0.972</td></tr><tr><td>El Arish</td><td>1754</td><td>-16 / 53 (3%)</td><td>205 / 177</td><td>206 / 185 (11%)</td><td>0.956 / 0.964</td></tr><tr><td>Sede Boqer</td><td>2087</td><td>-128 / -85 (-4%)</td><td>151 / 134</td><td>198 / 159 (8%)</td><td>0.978 / 0.984</td></tr><tr><td>Asyut</td><td>2092</td><td>94 / 84 (4%)</td><td>168 / 185</td><td>193 / 204 (10%)</td><td>0.962 / 0.951</td></tr><tr><td>Aswan</td><td>2230</td><td>-26 / 15 (1%)</td><td>175 / 190</td><td>177 / 191 (9%)</td><td>0.933 / 0.920</td></tr><tr><td>Tamanrasset</td><td>2319</td><td>58 / 76 (3%)</td><td>236 / 186</td><td>229 / 170 (7%)</td><td>0.903 / 0.951</td></tr><tr><td>Rochambeau</td><td>1750</td><td>157 / 123 (7%)</td><td>221 / 200</td><td>271 / 234 (13%)</td><td>0.927 / 0.934</td></tr><tr><td>Brasilia</td><td>1964</td><td>205 / 162 (8%)</td><td>257 / 228</td><td>328 / 280 (14%)</td><td>0.815 / 0.857</td></tr><tr><td>Bulawayo</td><td>1948</td><td>241 / 264 (14%)</td><td>300 / 292</td><td>385 / 394 (20%)</td><td>0.842 / 0.852</td></tr><tr><td>Maputo</td><td>1801</td><td>180 / 65 (4%)</td><td>215 / 201</td><td>281 / 211 (12%)</td><td>0.944 / 0.951</td></tr></tbody></table>

$I_{HC3McCle}$ is similar to or smaller than that for $I_{HC3v3}$ for all cases.

The standard deviation ranges from 20.6 to 27.1 J cm⁻² for $I_{HC3v3}$. In all cases, the standard deviation is smaller for $I_{HC3McCle}$ than that for $I_{HC3v3}$, and ranges from 19.7 to 25.2 J cm⁻². Like previously, the smaller standard deviation may be linked to the better correlation coefficient observed for $I_{HC3McCle}$. The RMSD is less for $I_{HC3McCle}$ than for $I_{HC3v3}$. It ranges from 19.8 to 25.3 J cm⁻². Like for 15 min irradiation, Table 4 shows that the combination of HC3 and McClear brings a benefit for 1 h irradiation for the six studied stations.

Table 5 reports the results of the comparison for daily irradiation. The correlation coefficient for $I_{HC3McCle}$ is large, greater than 0.93, except for Aswan (0.92), Brasilia (0.86) and Bulawayo (0.85). For all stations except Asyut and Aswan, the correlation is greater for $I_{HC3McCle}$ than for $I_{HC3v3}$. The day-to-day changes in daily irradiation are well reproduced by $I_{HC3v3}$, and slightly better by $I_{HC3McCle}$.

The bias ranges from -128 to 241 J cm⁻² for $I_{HC3v3}$. The bias for $I_{HC3McCle}$ is similar or smaller for 16 stations out of 23, and ranges from -85 to 264 J cm⁻². Several stations exhibit spectacular decreases, such as Rucana (from 91 down to -11 J cm⁻²), Thessaloniki (from -74 down to 3 J cm⁻²), or Maputo (from 180 down to 65 J cm⁻²). Seven stations exhibit

greater bias for $I_{HC3McCle}$ than for $I_{HC3v3}$: Valentia, Camborne, Nice, Mersa Matruh, El Arish, Tamanrasset, and Bulawayo.

The standard deviation for $I_{HC3McCle}$ ranges from 102 (Uccle) to 292 J cm⁻² (Bulawayo). It is similar to or less than that for $I_{HC3v3}$, except for Asyut and Aswan. The RMSD for $I_{HC3McCle}$ ranges from 102 to 394 J cm⁻², that is, from 6% to 20% of the mean observed value. It is similar to or less than that for $I_{HC3v3}$, with the exception of Asyut, Aswan and Bulawayo. Actually, the difference in standard deviation or RMSD is small for these three sites, and is less than 15 J cm⁻². This is less than the 66% uncertainty required by the World Meteorological Organization for the measurement of the daily irradiation (WMO, 2008), which is 40 J cm⁻² for $I_{ground} < 800$ J cm⁻² and 5% for $I_{ground} > 800$ J cm⁻². Taking this into account, it is found that $I_{HC3McCle}$ exhibits similar or better accuracy than $I_{HC3v3}$ for daily irradiation.

One may observe that the relative RMSD for $I_{HC3McCle}$ is less than 12% in most cases. Exceptions are Toravere (15%), Rucana (16%), Camborne (14%), Valentia (18%), Rochambeau (13%), Brasilia (14%), and Bulawayo (20%). These stations are seen with a large viewing angle by the Meteosat satellite. Schutgens and Roebeling (2009) or Marie-Joseph et al. (2013) argue that such angles induce a shift in the actual locations of clouds. The sensor aboard the Meteosat satellite
does not see exactly what is happening in the atmospheric column right above a measuring station. This contributes to the deviation between HC3 and ground measurements. The effects of the parallax are enhanced in the case of fragmented cloud cover, especially when the pixel size is large, which happens for large viewing angles. Marie-Joseph et al. (2013) mention that cloud fragmentation may contribute to a larger bias for intermediate skies because of the limited spatial resolution of the Meteosat sensor that prevents one from detecting small broken clouds such as cumulus. This patchwork of small clouds may be interpreted by the sensor and furthermore by the Heliosat-2 method as a large thin cloud. This mistake contributes to the deviation. As a rule of thumb, the farther from the nadir of the Meteosat satellite located at latitude 0° and longitude 0°, and the greater the occurrence of fragmented cloud cover, the greater the bias, relative standard deviation and RMSD.

## 4 Conclusions

This technical note proposes a very simple method to improve HC3v3 records by combining them with data records of the irradiation under clear skies from the new McClear clear-sky model. Inputs to McClear are the advanced global aerosol property forecasts and physically consistent total column content in water vapour and ozone produced by the MACC projects. All irradiation data sets may be retrieved on the SoDa website (http://www.soda-pro.com), and therefore the method is easily applicable. The method can be applied at any scale; it is not necessary to correct HC3v3 at 15 min resolution and then to sum up to obtain hourly or daily irradiation. Hourly and daily irradiation can be corrected using the corresponding irradiation from McClear.

The method has been validated against ground measurements made at several summarizations: 15 min, 1 h, and 1 day. The correlation coefficient is large, greater than respectively 0.92, 0.94, and 0.97, for 15 min, 1 h and daily irradiation. The bias ranges from -4 to 4% of the mean observed irradiation for most sites. The relative root mean square difference (RMSD) varies between 14 and 38% for 15 min, 12% and 33% for 1 h irradiation, and 6 and 20% for daily irradiation.

For all studied scales, 15 min, 1 h and 1 day, and almost all stations, the corrected irradiances $I_{HC3McClear}$ are closer to the ground-based measurements than those of $I_{HC3v3}$ obtained with a climatology of the Linke turbidity factor. There are few stations for which $I_{HC3McClear}$ does not show better performances than $I_{HC3v3}$, and in these cases, the difference is not large, and is less than the 66% uncertainty required for daily irradiation by the World Meteorological Organization (WMO, 2008). It is believed that the main cause of the benefit of this combination of HC3 and McClear is due to the fact that the inputs to McClear, aerosol properties and total column content in water vapour and ozone, are estimated

every 3 h. The main advantage of combining HC3v3 and McClear is that the large irradiances are reproduced better. The method brings an improvement in most cases and no degradation in the others, and a systematic correction of HC3v3 with McClear is recommended.

*Acknowledgements.* The research leading to these results has received funding from the European Union's Seventh Framework Programme (FP7/2007-2013) under grant agreement no. 283576 (MACC-II project) and no. 262892 (ENDORSE project). The authors thank all ground station operators of the BSRN network for their valuable measurements and the Alfred Wegener Institute for hosting the BSRN website, and acknowledge the help of the Baseline Surface Radiation Network (BSRN), the World Radiation Data Center (WRDC) and the meteorological networks supplying measurements of irradiation to the WRDC.

Edited by: R. Engelen

## References

Benedetti, A., Kaiser, J. W., and Morcrette, J.-J.: Global Climate, Aerosols in: State of the Climate in 2010, Bull. Am. Meteorol. Soc., 92, S65–S67, 2011.

Blanc, P., Gschwind, B., Lefevre, M., and Wald, L.: The HelioClim project: Surface solar irradiance data for climate applications, Remote Sens., 3, 343–361, doi:10.3390/rs3020343, 2011.

Elias, T. and Roujean, J.-L.: Estimation of the aerosol radiative forcing at ground level, over land, and in cloudless atmosphere, from METEOSAT-7 observation: method and case study, Atmos. Chem. Phys., 8, 625–636, doi:10.5194/acp-8-625-2008, 2008.

GCOS: Global Climate Observing System Essential Climate Variables, available at: www.wmo.int/pages/prog/gcos/index.php?name=EssentialClimateVariables, last access: 20 September 2013.

Geiger, M., Diabate, L., Menard, L., and Wald, L.: A web service for controlling the quality of measurements of global solar irradiation, Sol. Energy, 73, 475–480, doi:10.1016/S0038-092X(02)00121-4, 2002.

Inness, A., Baier, F., Benedetti, A., Bouarar, I., Chabriillat, S., Clark, H., Clerbaux, C., Coheur, P., Engelen, R. J., Errera, Q., Fleming, J., George, M., Granier, C., Hadji-Lazaro, J., Huijnen, V., Hurtmans, D., Jones, L., Kaiser, J. W., Kapsomenakis, J., Lefever, K., Leitão, J., Razinger, M., Richter, A., Schultz, M. G., Simmons, A. J., Suttie, M., Stein, O., Thépaut, J.-N., Thouret, V., Vrekoussis, M., Zerefos, C., and the MACC team: The MACC reanalysis: an 8 yr data set of atmospheric composition, Atmos. Chem. Phys., 13, 4073–4109, doi:10.5194/acp-13-4073-2013, 2013.

Kaiser, J. W., Peuch, V.-H., Benedetti, A., Boucher, O., Engelen, R. J., Holzer-Popp, T., Morcrette, J.-J., Wooster, M. J., and the MACC-II Management Board: The pre-operational GMES Atmospheric Service in MACC-II and its potential usage of Sentinel-3 observations, ESA Special Publication SP-708, Proceedings of the 3rd MERIS/(A)ATSR and OCLI-SLSTR (Sentinel-3) Preparatory Workshop, held in ESA-ESRIN, Frascati, Italy, 15–19 October 2012, 2012.
Lefèvre, M., Oumbe, A., Blanc, P., Espinar, B., Gschwind, B., Qu, Z., Wald, L., Schroedter-Homscheidt, M., Hoyer-Klick, C., Arola, A., Benedetti, A., Kaiser, J. W., and Morcrette, J.-J.: McClear: a new model estimating downwelling solar radiation at ground level in clear-sky conditions, Atmos. Meas. Tech., 6, 2403–2418, doi:10.5194/amt-6-2403-2013, 2013.

Lefevre, M., Blanc, P., Espinar, B., Gschwind, B., Menard, L., Ranchin, T., Wald, L., Saboret, L., Thomas, C., and Wey, E.: The HelioClim-1 database of daily solar radiation at Earth surface: an example of the benefits of GEOSS Data-CORE, IEEE J. Sel. Topics Appl. Earth Obs. Remote Sens., 7, 1745–1753, doi:10.1109/JSTARS.2013.2283791, 2014.

Marie-Joseph, I., Linguet, L., Gobindass, M.-L., and Wald, L.: On the applicability of the Heliosat-2 method to assess surface solar irradiance in the Intertropical Convergence Zone, French Guiana, Int. J. Remote Sens., 34, 3012–3027, doi:10.1080/01431161.2012.756598, 2013.

Peuch, V.-H., Rouil, L., Tarrason, L., and Elbern, H.: Towards European-scale Air Quality operational services for GMES Atmosphere, 9th EMS Annual Meeting, EMS2009-511, 9th European Conference on Applications of Meteorology (ECAM) Abstracts, held 28 September–2 October 2009, Toulouse, France, EMS2009-511, 2009.

Remund, J., Wald, L., Lefevre, M., Ranchin, T., and Page, J.: Worldwide Linke turbidity information, In Proceedings of ISES Solar World Congress, 16–19 June 2003, Göteborg, Sweden, CD-ROM published by International Solar Energy Society, 2003.

Rigollier, C., Bauer, O., and Wald, L.: On the clear sky model of the 4th European Solar Radiation Atlas with respect to the Heliosat method, Sol. Energy, 68, 33–48, doi:10.1016/S0038-092X(99)00055-9, 2000.

Rigollier, C., Lefevre, M., and Wald, L.: The method Heliosat-2 for deriving shortwave solar radiation from satellite images, Sol. Energy, 77, 159-169, doi:10.1016/j.solener.2004.04.017, 2004.

Roesch, A., Wild, M., Ohmura, A., Dutton, E. G., Long, C. N., and Zhang, T.: Assessment of BSRN radiation records for the computation of monthly means, Atmos. Meas. Tech., 4, 339–354, doi:10.5194/amt-4-339-2011, 2011.

Schutgens, N. A. J. and Roebeling, R. A.: Validating the validation: The influence of liquid water distribution in clouds on the intercomparison of satellite and surface observations, J. Atmos. Ocean. Technol., 26, 1457–1474, doi:10.1175/2009JTECHA1226.1, 2009.

WMO: Guide to Meteorological Instruments and Methods of Observation, World Meteorological Organization, WMO-No 8, 7th Edn., Geneva, Switzerland, 2008.

Xu, J., Li, C., Shi, H., He, Q., and Pan, L.: Analysis on the impact of aerosol optical depth on surface solar radiation in the Shanghai megacity, China, Atmos. Chem. Phys., 11, 3281–3289, doi:10.5194/acp-11-3281-2011, 2011.