# On the Assessment of Hourly Means of Solar Irradiance at Ground Level in Clear-Sky Conditions by the ERA5, JRA-3Q, and MERRA-2 Reanalyses

Yves-Marie Saint-Drenan, Lucien Wald

► To cite this version:

Yves-Marie Saint-Drenan, Lucien Wald. On the Assessment of Hourly Means of Solar Irradiance at Ground Level in Clear-Sky Conditions by the ERA5, JRA-3Q, and MERRA-2 Reanalyses. Atmosphere, 2025, 16, pp.949. 10.3390/atmos16080949. [hal-05206474]

HAL Id: hal-05206474

https://minesparis-psl.hal.science/hal-05206474v1

Submitted on 11 Aug 2025

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
Article

# On the Assessment of Hourly Means of Solar Irradiance at Ground Level in Clear-Sky Conditions by the ERA5, JRA-3Q, and MERRA-2 Reanalyses

Yves-Marie Saint-Drenan* and Lucien Wald

MINES Paris, PSL University, Centre Observation Impacts Energie (O.I.E.), 06904 Sophia-Antipolis, France; lucien.wald@minesparis.psl.eu

* Correspondence: yves-marie.saint-drenan@minesparis.psl.eu

## Abstract

Meteorological reanalyses are one of the means to assess the solar irradiance reaching the ground. This paper deals with estimates of the hourly means of irradiance in clear-sky conditions provided by the ERA5, JRA-3Q, and MERRA-2 reanalyses. They are compared to coincident ground-based measurements from 28 BSRN stations located worldwide, selected by a new algorithm for detecting cloud-free instants. Although ERA5 most often underestimates measurements, it is quite reliable over time because it captures the temporal variability of measurements well and provides a constant level of uncertainty. JRA-3Q offers a complex pattern with negative and positive biases depending on station and season. It captures well the temporal variability but, as a whole, is not reliable over time. None of the three reanalyses is reliable in space. Because of its use of the mean solar time instead of the true solar time, MERRA-2 suffers many drawbacks over intraday scales. Its statistical indicators exhibit marked patterns depending on the season and station. Its assimilation of aerosol properties offers advantages when compared to the climatologies used in ERA5 and JRA-3Q. This work exposes the strengths and weaknesses of each reanalysis in clear-sky conditions and formulates suggestions to providers for further improvements.

**Keywords:** solar radiation; clear-sky; meteorological reanalysis; BSRN; aerosol; detection of clear-sky instants

## 1. Introduction

The power of the electromagnetic radiation from the sun that reaches the surface of the Earth is the main driver behind the weather and climate systems on the planet. The density of power received from the sun on a horizontal surface at ground level per unit surface is called the solar irradiance at the surface (SIS) and is expressed in W m⁻². It is hereafter abbreviated as SIS and noted *E*. Equivalent terms may be found in the literature, such as solar flux, surface solar irradiance, downwelling solar irradiance, downwelling shortwave flux, downwelling surface shortwave flux, surface incoming solar radiation, surface incoming shortwave radiation, surface net downward shortwave flux, or global horizontal irradiance. The SIS is usually integrated over the solar spectrum (total SIS).

The Global Climate Observing System (GCOS) program has identified the SIS as an essential climate variable that helps us understand climate evolution and guides adaptation and mitigation efforts [1,2]. The SIS is also one of the essential variables in air quality [3], the terrestrial and marine environment [4,5], and renewable energies [6]. In addition to

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

check for
updates

Academic Editor: Filomena Romano

Received: 24 June 2025

Revised: 11 July 2025

Accepted: 14 July 2025

Published: 7 August 2025

**Citation:** Saint-Drenan, Y.-M.; Wald, L. On the Assessment of Hourly Means of Solar Irradiance at Ground Level in Clear-Sky Conditions by the ERA5, JRA-3Q, and MERRA-2 Reanalyses. *Atmosphere* **2025**, *16*, 949. https://doi.org/10.3390/atmos16080949

**Copyright:** © 2025 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (https://creativecommons.org/licenses/by/4.0/).
these domains, it impacts on human health [7] and many other aspects of our daily lives and our activities (see many examples in [8] or [9] (pp. xxii–xxviii)). As a result, there is a significant and rising demand for precise and trustworthy estimates and measurements of the SIS, anytime, anywhere.

Numerical weather prediction (NWP) models, whether used for forecasting or reanalysis, are one of the methods used to estimate SIS. They can offer a global coverage and span of several years. This study focuses on three widely used atmospheric reanalysis datasets: the ECMWF Reanalysis version 5 (ERA5) produced by the European Centre for Medium-Range Weather Forecasts (ECMWF), the Modern-Era Retrospective analysis for Research and Applications version 2 (MERRA-2) developed by NASA’s Global Modeling and Assimilation Office (GMAO), and the Japanese 55-year Reanalysis with Quarterly updates (JRA-3Q) provided by the Japan Meteorological Agency (JMA). Reanalyses derive from weather forecast models used in a reanalysis mode to replicate observed conditions. A wide range of observations—including weather station measurements, satellite data, and GPS signals—are assimilated to estimate atmospheric state variables such as temperature, humidity, and wind speed. A description of the data assimilated in the three reanalysis is given in [10–13]. In contrast, measurements of the SIS are not assimilated in the three considered reanalyses and the SIS is a diagnostic quantity, derived from all related variables and calculated using a radiative transfer model. It depends on the complete set of radiatively active variables in the atmospheric column above the specific location. Therefore, reanalysis estimates of the SIS should not be confused with direct measurements as they carry the uncertainty inherent in the model simulations such as limited resolution of the model, uncertainties in the parameterizations used, and individual uncertainty of observations.

The SIS $E$ is often written as the product of the SIS received in cloud-free conditions, noted as $E_{clear}$, and a cloud modification factor, also known as the clear-sky index, noted as $Kc$.

$$E = E_{clear} Kc. \qquad (1)$$

A clear-sky index close to 1 would imply a sky without clouds, and a value close to 0 would imply a sky covered by optically thick clouds. Over short time periods, especially in broken-cloud conditions, $Kc$ might be greater than 1, resulting from a so-called cloud enhancement effect. Situations may also arise in which $Kc$ exceeds 1 over extended periods, notably due to the combination of high surface reflectivity and low solar elevation.

Oumbe et al. [14] showed that at first approximation, $Kc$ does not depend on the properties of the cloud-free atmosphere and depends only on the cloud properties and the reflective properties of the ground. This approximation yields errors comparable to the typical uncertainties of pyranometers, except when the ground albedo exceeds 0.7. This result shows that the effects of the cloudless atmosphere are entirely conveyed in the modelling of $E_{clear}$ and highlights the importance of accurate assessment of $E_{clear}$. Clear-sky models aim at estimating $E_{clear}$. They provide realistic upper limits of the SIS and contribute indirectly to quantifying the radiative effects of clouds. $E_{clear}$ varies with the date, time of day, geographic location, and gas concentrations. For total SIS, aerosol, water vapor, and ozone are the most influent variables, while other gases can be assigned standard concentrations. Extinction also depends on local temperature, density, pressure, and volume mixing ratio for the non-ozone, non-water vapor gases, necessitating knowledge of the vertical profiles of these variables, which is crucial for calculating the scattering effects from air molecules. Additionally, $E_{clear}$ is influenced by the elevation and reflective properties of the ground.

The present paper builds upon the result from [14] and focuses on the estimates of the hourly means of $E_{clear}E_{clear}^{h}$, made by each of the three reanalyses considered in the present work. More precisely, it aims at assessing the quality of these estimates by a comparison
against coincident high-quality ground-based measurements of $E_{clear}^h$. It is an attempt to identify potential drawbacks in the modelling of $E_{clear}^h$ and may help the agencies in charge of the reanalyses in further developments as uncertainties in $E_{clear}^h$ affect the estimates of the SIS in all conditions.

There are many scientific articles relating to the comparison between estimates of $E_{clear}$ from reanalyses and ground-based measurements on hourly, daily, monthly, or yearly basis. To the best of our knowledge, based on a comprehensive survey of the literature, none of them was dedicated to clear-sky conditions, except [15], and none of them explored the clear-sky products of these reanalyses except [15,16].

Wang et al. [15] compared the monthly mean of the clear-sky product of MERRA-2 $E_{clear\_merra2}^{month}$ to ground-based measurements of $E_{clear}^{month}$ under cloudless skies in Asia obtained by thresholding the monthly clearness index $KT^{month}$, where the clearness index is the ratio of $E_{clear}$ to the corresponding extraterrestrial irradiance received on a horizontal surface located at the top of the atmosphere, noted as $E_0$. They found that MERRA-2 overestimates $E_{clear}^{month}$ and $KT^{month}$ significantly. Working with hourly measurements in the tropical Atlantic Ocean, Trolliet et al. [16] did the sole study comparing the hourly clear-sky product from MERRA-2, $E_{clear\_merra2}^h$, and the corresponding SIS $E_{McClear}^h$ from the McClear service of the Copernicus Atmosphere Monitoring Service (CAMS) [17,18], whose quality is well documented [19]. Without providing any number, they reported that these datasets are similar and concluded that the MERRA-2 cloud-free SIS are likely accurate. As clear-sky products of ERA5 were unavailable at that time, the authors performed a visual inspection of daily profiles of $E_{era5}^h$ for randomly selected cloud-free days on which daily profiles of $E_{McClear}^h$ were superimposed and found that $E_{McClear}^h$ was underestimated by ERA5 in these cases.

In addition to [15,16], several studies have partially investigated clear-sky conditions. They dealt with ground-based measurements in several places: Scandinavia [20], South Africa [21], West Africa [22], Pakistan [23,24], Indonesia [25], or worldwide [26]. Clear-sky conditions were determined using various approaches. Some works select clear sky instants by thresholding the daily clearness index $KT^{day}$ [22,23,25,26], or modified clearness index [24], and others use the daily clear-sky index when an external clear-sky model is used [20]. Mabasa et al. [21] did not detail how the clear-sky days were determined, but from [19], it appears that they used ERA5 in a first step and then a combination of a modified version of the algorithm of Reno and Hansen [27] and McClear outputs.

This survey of the literature reveals several discrepancies in results. Several of them may be partly explained by discrepancies in the studied timescale, from hour to month. For example, Wang et al. [15] report that MERRA-2 overestimates $E_{clear}^{month}$ and $KT_{clear}^{month}$ significantly while others [23–26] wrote that MERRA-2 slightly underestimates the daily means $E_{clear}^{day}$, and Trolliet et al. [16] concluded that MERRA-2 cloud-free hourly means, $E_{clear\_merra2}^h$, are likely accurate. Other partial explanations of the discrepancies may be due to differences in the studied areas and periods as the quality of the estimates of the cloudless SIS may depend on the climate and changes in aerosol properties and water vapor among others. Discrepancies are also partly due in the manner of determining clear-sky conditions, as discussed above, as it has a great influence on results [19,28].

From this survey of the literature, we observed several discrepancies in results that may be partly explained by discrepancies in the studied timescale, discrepancies in the geographical area, and by discrepancies in determining clear-sky conditions, the latter having a great influence on results [19,28].

Sianturi et al. [25] reported that ERA5 estimates $E_{clear}^{day}$ fairly correctly, in contradiction with other studies [23,26] that found a slight underestimation of $E_{clear}^{day}$ or an overestimation [20]. ERA5 was found to report cloudy conditions while the actual conditions are
cloud-free [16,21,26], and in another study ground-based measurements and ERA5 agree on the clear-sky category 77% of the time [20]. As for JRA-55, the predecessor of JRA-3Q, the sole study reported that JRA-55 tends to underestimate $E_{clear}^{day}$ [23]. As for MERRA-2, it slightly underestimates $E_{clear}^{day}$ [24–26] but Wang et al. found that it overestimates $E_{clear}^{month}$ significantly.

These works give some indications of the ability of ERA5, JRA-3Q, and MERRA-2 to reproduce $E_{clear}$ on hourly, daily, and monthly bases though they provide little quantitative measurement of it, except for [15]. Our study focuses on hourly means of the clear-sky SIS $E_{clear}^h$ and aims at extending the current knowledge as we believe that it would be fruitful to better know the strengths and weaknesses of the reanalyses in clear-sky conditions for further improvements in the diagnostic of E in all conditions. The present paper does not intend to rank reanalyses. There are many reasons for choosing one reanalysis over others; a ranking based solely on comparisons with in situ measurements at a limited number of stations during a limited period is unlikely to cover all needs and constraints.

Figure 1 exhibits an overview of the methodology adopted for this paper, with details given later for each step. The ground-based measurements collected every 1 min in the BSRN network [29], noted as $E_{bsrn}^{min}$, are checked for their plausibility. Detection of clear-sky instants is made on the retained 1 min data, yielding the time series $E_{clear\_bsrn}^{min}$. From this series, hourly means of SIS in cloudless conditions, noted as $E_{ref}^h$, are created and are compared to the coincident estimates of $E_{clear}^h$ given by each reanalysis. The comparisons are also performed on the clearness index $KT_{ref}^h$. As the estimates of the SIS depend on the calculation of the extraterrestrial irradiance $E_0^h$ made by each reanalysis, a comparison is performed between the various $E_0^h$ to check their consistency. This is justified in the case of MERRA-2 as it was reported in [16] that MERRA-2 is using the mean solar time and does not account for the difference between the true solar time and the mean solar time, which is a function of the day in the year as a first approximation, thus bringing additional uncertainty in the estimates of $E^h$.

<figure><img src="image_3.png" /><figcaption>**Figure 1.** Scheme of the methodology. E and KT are, respectively, the irradiance and clearness index. Superscripts “min” and “h” mean, respectively, 1 min and hourly means. Subscripts “0” and “clear” mean, respectively, extraterrestrial and clear-sky. Subscript “ref” means reference values.</figcaption></figure>

The data and models are briefly presented in Section 2. The protocol for validation is described in Section 3. The verification of the extraterrestrial irradiance $E_0^h$ is discussed in Section 4. There are several ways of discussing results when several models are dealt with. Here, the selected approach is to present and discuss the results of each reanalysis separately in Section 5, Section 6, and Section 7, respectively, so that each Section provides
the reader with a full view of the performances of a given reanalysis. Section 8 concludes the paper. How to access the data is described in the “Data availability” section.

## 2. Data and Models

Data used here are measurements of $E_{bsrn}^{min}$ performed in the BSRN network, estimates of $E_{clear}^{h}$ from the reanalyses ERA5, JRA-3Q, and MERRA-2, as well as estimates of several cloudless atmosphere properties derived from CAMS.

### 2.1. Selected Sites and Periods

Ground-based measurements were obtained at 28 sites located worldwide belonging to the Baseline Surface Radiation Network (BSRN) [29] (Figure 2 and Table 1). These stations exist in a variety of climates. In order to perform trimester-based comparisons consistent with climatological convention (DJF, MAM, JJA, and SON), each year is defined as beginning in the preceding December. For instance, DJF 2017 covers the period from December 2016 to February 2017. With this convention, the selected periods are [2015-12-01 to 2016-11-30], [2016-12-01 to 2017-11-30], and [2017-12-01 to 2018-11-30]. For the sake of conciseness, they are designated as 2016, 2017, and 2018. These three years were chosen as they represent the most recent periods with the greatest number of BSRN stations providing high data availability.

<figure><img src="image_4.png" /><figcaption>Figure 2. Map showing the locations of the selected BSRN sites (black diamonds).</figcaption></figure>

<table><caption>Table 1. BSRN stations, from the northernmost station to the southernmost one.</caption><thead><tr><td>Station (Code)</td><td>Code</td><td>Latitude (°)</td><td>Longitude (°)</td><td>Elevation * (m)</td><td>Köppen-Geiger Class **</td><td>Spectral Range (nm) ***</td></tr></thead><tbody><tr><td>Ny-Alesund</td><td>NYA</td><td>78.925</td><td>11.930</td><td>11</td><td>ET</td><td>[210, 3600]</td></tr><tr><td>Toravere</td><td>TOR</td><td>58.254</td><td>26.462</td><td>70</td><td>Dfb</td><td>[285, 2800]</td></tr><tr><td>Lindenberg</td><td>LIN</td><td>52.210</td><td>14.122</td><td>125</td><td>Cfb</td><td>[200, 3600]</td></tr><tr><td>Cabauw</td><td>CAB</td><td>51.971</td><td>4.927</td><td>0</td><td>Cfb</td><td>[200, 3600]</td></tr><tr><td>Palaiseau</td><td>PAL</td><td>48.713</td><td>2.208</td><td>156</td><td>Cfb</td><td>[210, 3600]</td></tr><tr><td>Fort Peck</td><td>FPE</td><td>48.308</td><td>-105.102</td><td>634</td><td>BSk</td><td>[285, 3000]</td></tr><tr><td>Payerne</td><td>PAY</td><td>46.815</td><td>6.944</td><td>491</td><td>Cfb</td><td>[210, 3600]</td></tr><tr><td>Carpentras</td><td>CAR</td><td>44.083</td><td>5.059</td><td>100</td><td>Csa</td><td>[305, 2800]</td></tr></tbody></table>
<table><caption>Table 1. *Cont.*</caption><thead><tr><th>Station (Code)</th><th>Code</th><th>Latitude (°)</th><th>Longitude (°)</th><th>Elevation * (m)</th><th>Köppen-Geiger Class **</th><th>Spectral Range (nm) ***</th></tr></thead><tbody><tr><td>Sapporo</td><td>SAP</td><td>43.060</td><td>141.329</td><td>17</td><td>Dfb</td><td>[285, 2800]</td></tr><tr><td>Cener</td><td>CNR</td><td>42.816</td><td>-1.601</td><td>471</td><td>Cfb</td><td>[210, 3600]</td></tr><tr><td>Boulder</td><td>BOS</td><td>40.125</td><td>-105.237</td><td>1689</td><td>BSk</td><td>[285, 3000]</td></tr><tr><td>Langley Res. Ctr.</td><td>LRC</td><td>37.104</td><td>-76.387</td><td>3</td><td>Cfa</td><td>[200, 3600]</td></tr><tr><td>Desert Rock</td><td>DRA</td><td>36.624</td><td>-116.019</td><td>1007</td><td>BWk</td><td>[285, 3000]</td></tr><tr><td>Billings</td><td>BIL</td><td>36.605</td><td>-97.516</td><td>317</td><td>Cfa</td><td>[285, 2800]</td></tr><tr><td>Southern Great Plains</td><td>E13</td><td>36.604</td><td>-97.485</td><td>314</td><td>Cfa</td><td>[285, 2800]</td></tr><tr><td>Tateno</td><td>TAT</td><td>36.058</td><td>140.126</td><td>25</td><td>Cfa</td><td>[285, 2800]</td></tr><tr><td>Fukuoka</td><td>FUA</td><td>33.582</td><td>130.376</td><td>3</td><td>Cfa</td><td>[285, 2800]</td></tr><tr><td>Minamitorishima</td><td>MNM</td><td>24.288</td><td>153.983</td><td>7</td><td>Ocean</td><td>[285, 2800]</td></tr><tr><td>Petrolina</td><td>PTR</td><td>-9.068</td><td>-40.319</td><td>387</td><td>BSh</td><td>[200, 3600]</td></tr><tr><td>Darwin</td><td>DWN</td><td>-12.424</td><td>130.893</td><td>32</td><td>Aw</td><td>[305, 2800]</td></tr><tr><td>Gobabeb</td><td>GOB</td><td>-23.561</td><td>15.042</td><td>407</td><td>BWh</td><td>[210, 3600]</td></tr><tr><td>Alice Springs</td><td>ASP</td><td>-23.795</td><td>133.889</td><td>546</td><td>BWh</td><td>[285, 2800]</td></tr><tr><td>Florianopolis</td><td>FLO</td><td>-27.602</td><td>-48.518</td><td>31</td><td>Cfa</td><td>[305, 2800]</td></tr><tr><td>De Aar</td><td>DAA</td><td>-30.667</td><td>23.993</td><td>1287</td><td>BSk</td><td>[285, 2800]</td></tr><tr><td>Lauder</td><td>LAU</td><td>-45.045</td><td>169.689</td><td>350</td><td>Cfb</td><td>[305, 2800]</td></tr><tr><td>Syowa</td><td>SYO</td><td>-69.005</td><td>39.589</td><td>18</td><td>EF</td><td>[285, 2800]</td></tr><tr><td>Georg von Neumayer</td><td>GVN</td><td>-70.650</td><td>-8.250</td><td>42</td><td>Ocean</td><td>[210, 3600]</td></tr><tr><td>Concordia Station, Dome C</td><td>DOM</td><td>-75.100</td><td>123.383</td><td>3233</td><td>EF</td><td>[200, 3600]</td></tr></tbody></table>

## 2.2. Reanalyses

This study deals with ERA5 developed at the ECMWF (European Centre for Medium-Range Weather Forecasts) [10], JRA-3Q (Japanese Reanalysis for Three Quarters of a Century) of the Japan Meteorological Agency [12], a follow-on of the JRA-55 reanalysis [11], and MERRA-2 (Modern-Era Retrospective Analysis for Research and Applications, version 2) of the USA National Aeronautics and Space Administration [13]. The sizes of the grid cell of the solar radiation variables are 0.28125° × 0.28125° (31 km in latitude) for ERA5, 0.375° × 0.375° (42 km in latitude) for JRA-3Q, and 0.5° in latitude (56 km) by 0.625° in longitude (70 km at the Equator) for MERRA-2.

Of particular importance for the modelling of the downwelling solar radiation in clear-sky conditions at surface is the modelling of the aerosol properties. ERA5 exploits a prescribed monthly climatology from Tegen et al. [30] for tropospheric aerosols and Tanre et al.’s [31] aerosol climatology as well as GISS dataset [32] for stratospheric aerosols [10]. Similarly, JRA-3Q, like JRA-55, relies on a fixed aerosol climatology derived from the MODTRAN/GINA model [11]. In contrast, MERRA-2 includes a system based on the GOCART model, with assimilation of satellite-derived aerosol optical depth (e.g., MODIS and AVHRR), and a direct coupling to radiative transfer [13,33].

<footnote>* Above mean sea level. ** From the BSRN website—Aw: tropical climate in savannah; BSh: arid and hot steppe; BSk: arid and cold steppe; BWh: arid and hot, and of desert type; BWk: arid and cold, and of desert type; Cfa: temperate, no dry season, and hot summer; Cfb: temperate, no dry season, and warm summer; Csa: temperate, no dry season, and warm summer; Dfb: cold, no dry season, and warm summer; EF: polar, ice cap; and ET: polar tundra. *** Spectral band at 50% sensitivity of the pyranometer measuring $E_{bsrn}^{min}$</footnote>
The solar radiation at the top of the atmosphere and downwelling solar radiation in clear-sky conditions at surface were downloaded from the ECMWF Copernicus Climate Datastore (CDS) for ERA5, Thredds data server operated by UCAR for JRA-3Q, and NASA Earthdata website for MERRA-2 at the four nearest grid cells to each site. All time series are in UT time system. A spatial bilinear interpolation technique with a weighting factor that is inversely proportional to the distance to the site was then applied to yield time series of hourly means of irradiance at the site, noted, respectively, as $E_{0\_era5}^h$, $E_{clear\_era5}^h$, $E_{0\_jra3q}^h$, $E_{clear\_jra3q}^h$, $E_{0\_merra2}^h$, and $E_{clear\_merra2}^h$. Then time series of clearness index are derived as follows:

$$KT_{clear\_era5}^h = \frac{E_{clear\_era5}^h}{E_{0\_era5}^h}, \qquad (2)$$

$$KT_{clear\_jra3q}^h = \frac{E_{clear\_jra3q}^h}{E_{0\_jra3q}^h}, \qquad (3)$$

$$KT_{clear\_merra2}^h = \frac{E_{clear\_merra2}^h}{E_{0\_merra2}^h}. \qquad (4)$$

## 2.3. Astronomical Quantities and Others

Several astronomical quantities were extracted from CAMS in UT time to help analyse the results of the comparisons. For convenience, we used the McClear service in CAMS, which in its verbose mode delivers among other quantities 1 min values of readings from CAMS resampled to the selected location by spatial bilinear interpolation and resampled in time by linear interpolation, namely the solar zenith angle $\theta_S^{min}$ calculated with the SG2 algorithm [34], and the irradiance at the top of the atmosphere $E_0^{min}$ using the nominal total solar irradiance of 1361 W m$^{-2}$, as recommended by [35]. A second call to the McClear service, this time for hourly values, returned the irradiance at the top of the atmosphere, $E_0^h$, as well as the hourly direct components of the clear-sky model McClear on a horizontal surface, $B_{McClear}^h$, and on a plane always normal to the sun's rays, $B_{N\_McClear}^h$. The last two quantities were used to calculate an effective solar zenith angle $\theta_S^{eff\_h}$ for each hour [36]:

$$\theta_S^{eff\_h} = \cos^{-1} \left( \frac{B_{McClear}^h}{B_{N\_McClear}^h} \right). \qquad (5)$$

Because of the high accuracy of the SG2 algorithm, the irradiance $E_0^h$ is considered as a reference to which $E_{0\_era5}^h$, $E_{0\_jra3q}^h$, and $E_{0\_merra2}^h$ are compared.

## 2.4. BSRN Measurements and Their Filtering

Stations in the BSRN network perform high-quality measurements of the SIS and its diffuse and direct components [29]. Each station is equipped with class-A thermopile pyranometers, one of them being equipped with a shadow ball mounted on a tracker, to separately measure $E_{bsrn}^{min}$ and $D_{bsrn}^{min}$, and a pyrheliometer to measure the direct component at normal incidence, $B_{N\_bsrn}^{min}$. The measurements have a temporal resolution of 1 min. Uncertainty requirements are 5 W m$^{-2}$ for $E_{bsrn}^{min}$ and $D_{bsrn}^{min}$ (2% in relative value) and 2 W m$^{-2}$ for $B_{N\_bsrn}^{min}$ (0.5%) [29].

The BSRN-recommended quality check [37] was applied to the measurements and only those flagged as non-suspicious were retained here. For the sake of clarity in writing,
the retained measurements are said to be “valid”. For further use, the clearness index is computed as follows:

$$KT_{bsrn}^{min} = \left( \frac{E_{bsrn}^{min}}{E_{0}^{min}} \right). \qquad (6)$$

## 2.5. Screening of Cloudless Instants

A new algorithm for detecting cloud-free instants has been devised taking into account the comments of Bright et al. [28] on the existing ones. Its principle is that the 1 min SIS in clear-sky conditions exhibit smooth variations, or, said differently, its derivative is continuous. It is also the case for the SIS in overcast conditions, and such cases must be excluded. As it is usually impossible to compute derivative from fragmented series of measurements, the problem is reversed to keep only measurements that are close to an analytical function that fits typical SIS in clear-sky conditions. The proposed algorithm is meant to reliably detect cloud-free instants and not all cloud-free instants. Hence, it may reject actual cloud-free measurements.

The algorithm comprises several cascading dyadic multi-scale filters that apply to the valid 1 min measurements for each day of the time series. It exploits the form of the MLB function [38] for interpolating the SSI as a function of the solar zenith angle $\theta_S^{min}$ in clear-sky conditions:

$$E_{MLB}^{min} = E_{0}^{min} \exp \left[ \frac{b}{\cos(\theta_S^{min})^a} \right], \qquad (7)$$

where the parameters a and b are adjusted by a least-squares technique on all valid $E_{bsrn}^{min}$ within a sliding window of width $\pm n/2$ min, centred on each instant t, expressed in min, in the day. Each interval $[t - n/2, t]$ and $[t, t + n/2]$ must contain at least 30% of valid values, $B_{N_{bsrn}}^{min}$, greater than 10 W m$^{-2}$. Otherwise, the instant t is said to be “likely cloudy” and the next instant, t, is processed. This constraint is meant to discard homogeneous overcast conditions. If the constraint is met, the MLB function is adjusted, and a series of irradiances noted as $E_{MLB}^{min}$ is computed with Equation (7).

Over each window of length 2*n*, the covariance $Cov_n^{min}$ between $E_{MLB}^{min}$ and valid $E_{bsrn}^{min}$ and the variances $Var_{bsrn,n}^{min}$ and $Var_{MLB,n}^{min}$ of the valid $E_{bsrn}^{min}$ and coincident $E_{MLB}^{min}$ are computed for each time step, t, in the window $[t - n/2, t]$ and $[t, t + n/2]$. The instant t is declared as “likely cloud-free” if

$$Cov_n^{min} > \rho_n^{min} \times \sqrt{Var_{MLB,n}^{min}} \times \sqrt{Var_{bsrn,n}^{min} + \sigma^2} \qquad (8)$$

where $\rho_n^{min}$ and $\sigma$ were set empirically to 0.999 and 4 W m$^{-2}$.

This procedure generates a time series of binary values denoting the results of the clear-sky detector over a temporal window of length 2*n*. Four windows were used, with n successively being 120, 60, 30, and 15. Our method finally flags as cloud-free the instants classified as cloud-free over all considered windows. The retained 1 min values are noted as $E_{clear\_bsrn}^{min}$ and $KT_{clear\_bsrn}^{min}$.

The algorithm was first verified with ideal cases of clear conditions simulated by the McClear service, then with simulated cases of clear skies including the slow appearance of thin clouds, and, eventually, by visual inspection of randomly selected days in the BSRN dataset. Its outcomes were also compared to those of the algorithms of Lefèvre et al. [17] and Reno and Hansen [27] for a further check of the proper functioning of the algorithm, though the objective of this paper is not to compare the results of our algorithm to others, as in [28].
Figure 3 displays several daily profiles of $E_{bsrn}^{min}$ at Payerne, on which symbols representing the clear-sky instants detected by our algorithm (red dots) are superimposed. The sequence of five days displayed in Figure 3 includes various sky conditions varying from overcast to cloud-free conditions. As expected, instants with overcast (25 April 2015) and broken clouds (26 April) are rejected by the clear-sky detector. Our algorithm detected clear-sky periods during mixed conditions on 23 and 24 April 2015, while instants with a high variability were excluded. Note that on these two days there was a latency between a period with fluctuating SIS and a detected cloud-free condition. This is intentional to ensure that the selected instants are cloud-free conditions for sure. Finally, on 22 April 2015, the complete day was detected as being cloud-free when the SIS was greater than 200 W m⁻². The correct behaviour of our clear-sky detection algorithm was further verified by a visual inspection of other periods and stations.

<figure><img src="image_5.png" /><figcaption>**Figure 3.** Illustration of the detection of clear-sky instants. Daily profiles of $E_{bsrn}^{min}$ at Payerne on which the blue line represents E measurements and the red dots the clear-sky instants detected by our algorithm.</figcaption></figure>

## 2.6. Hourly Values

From the time series of $E_{clear\_bsrn}^{min}$, hourly means, $E_{ref}^{h}$, adjusted onto the time-stamping in UT used by the three reanalyses were constructed by a linear interpolation in $KT_{clear\_bsrn}^{min}$ for filling in the gaps, provided there are at least 20 known $KT_{clear\_bsrn}^{min}$ in the hour. Otherwise, the hour is said to be invalid. The hourly clearness indices, $KT_{ref}^{h}$, were computed as the ratio of $E_{ref}^{h}$ to $E_{0}^{h}$. These series $E_{ref}^{h}$ and $KT_{ref}^{h}$ are the reference against which the estimates of the reanalyses are compared. Their means and variances are termed reference means and variances hereafter.

Figure 4 exhibits the number of hours at each station and all stations merged for each period and all periods. Only periods with at least 50 valid hourly values were retained. The number of stations times the number of valid trimesters makes a total of 333 trimesters. Means of $E_{ref}^{h}$ and $KT_{ref}^{h}$ at each station and all stations merged are shown in Figures 5 and 6 for each period and all years. The minimum and maximum of $E_{ref}^{h\_mean}$ are 306 W m⁻² (DJF207 at Palaiseau (PAL)) and 846 W m⁻² (DJF2017 at Florianopolis (FLO)); those of $KT_{ref}^{h\_mean}$ are 0.654 (SON2018 at Toravere) and 0.900 (DJF2015 at Concordia Dome C (DOM)). The greatest $KT_{ref}^{h\_mean}$ are observed at the three Antarctic stations. The snow-covered ground reflects part of the solar irradiance and contributes to the measured SIS, thus increasing the clearness index. In addition, the Concordia Dome C is located at
3233 m above the mean sea level, which means a thinner atmospheric layer and, therefore, a greater clearness index.

Number of cloudless hourly values

<figure><img src="image_6.png" /><figcaption>**Figure 4.** Number of reference hourly means, $E_{ref}^h$, at each station for each trimester. Station codes are given in Table 1; bold italic: stations located beyond the polar circles (latitudes > 66.5° or <−66.5°); and bold: stations located in the tropical belt [−23.5°, 23.5°]. The quantities represented are dimensionless and, therefore, have no units.</figcaption></figure>
<figure><img src="image_7.png" /><figcaption>**Figure 5.** Means of $E_{ref}^h$, at each station and all stations merged (All) for each trimester and all trimesters (All periods). Station codes are given in Table 1; bold italic: stations located beyond the polar circles (latitudes > 66.5° or <−66.5°) and bold: stations located in the tropical belt [−23.5°, 23.5°].</figcaption></figure>
<figure><img src="image_8.png" /><figcaption>**Figure 6.** Means of $KT_{ref}^h$ at each station and all stations merged (All) for each trimester and all trimesters (All periods). Station codes are given in Table 1; bold italic: stations located beyond the polar circles (latitudes > 66.5° or <−66.5°) and bold: stations located in the tropical belt [−23.5°, 23.5°]. The quantities represented are dimensionless and, therefore, have no units.</figcaption></figure>

**3. Protocol for Comparison**

There are two kinds of comparison, one deals with the extraterrestrial irradiance, the other with the SIS. The protocol is fairly similar in both cases and applies to each reanalysis in turn. The comparison for the SIS was performed at each BSRN site and for all sites merged, for each trimester and for the ensemble of trimesters.

For each reanalysis, the comparison was made on the estimates provided by this reanalysis and either the extraterrestrial irradiance from the CAMS McClear service $E_0^h$, or reference values $E_{ref}^h$ and $KT_{ref}^h$ for coincident instants only. The Pearson correlation coefficients, slopes, and offsets of the least-squares fitting lines were computed as well as the differences between coincident estimates and reference values, which were summarized
by their mean, known as the bias or mean bias error, their standard deviation, and their root mean square error (RMSE). Relative biases, relative standard deviations, and relative RMSEs were expressed with respect to the reference means. Two-dimensional histograms of reference values and estimates were drawn. In addition, results were analysed as a function of the station, latitude, climate, elevation, year, $\theta_S^{eff-h}$, reference means, and variances to evidence possible trends. All produced graphs are available as Supplementary Materials.

The present protocol is not new and is similar to many used in comparable studies. It presents several drawbacks that may affect the comparison.

First of all, the reanalyses provide total irradiance, i.e., integrated over the whole spectrum, while BSRN measurements taken by pyranometers are based on a more limited range, as shown in Table 1. The often-used BSRN spectral range [285, 2800] nm accounts for about 99% of the total irradiance with slight dependence on the optical properties of the clear atmosphere constituents and the air mass. The 1% irradiance difference resulting from the limited spectral sensitivity of the pyranometer is partly compensated for during the assessment of the sensitivity coefficient of the instrument, which is the conversion factor between measured thermopile voltage and solar irradiance. The sensitivity coefficient is determined during the calibration of the instrument by comparing readings of the pyranometer and the World Radiometric Reference given by cryogenic solar absolute radiometers that measure total irradiance with a very high accuracy (see [39,40] for details). Consequently, the sensitivity coefficient implicitly accounts for the average difference between the irradiance in the interval [285, 2800] nm and the total irradiance. This suggests that a slight difference, with a zero mean and dependence on air mass and atmospheric content, may exist between the reanalyses and measurements due to the difference in spectral ranges (see [19] for further discussion on the difference in spectral ranges between measurements and models).

A second drawback is due to the mismatch in spatial representation of the field of *E* by the reanalyses and the BSRN pyranometers. The reanalyses provide *E* on cells several kilometers in size and, thus, do not contain the sub-cell variability of *E* that is partly perceived by the pyranometers. This adds to the discrepancy expected between the reanalyses' outputs and $E_{ref}^h$, whose magnitude cannot be predicted as it depends on the atmospheric conditions experienced by each station every 1 h. This is a well-known issue that encompasses the present case and applies to many domains in Earth sciences and others. Several examples are provided in [19] showing the large spatial variability of the SIS.

Also to be taken into account is the site elevation as well as the presence of seacoasts or mountains and, more generally, of orographic discontinuities. Urraca et al. [26] reported an increased positive bias in *E* due to the large size of the grid cells of ERA5 or MERRA-2 over the coasts where the spatial variability in *E* may be large. This effect may be minimized by a proper choice of BSRN stations located in homogeneous areas with respect to orography, vegetation, or the presence of large water bodies. The limited horizontal and vertical resolutions of the elevation model used in each reanalysis can also bring about a systematic error when the assumed elevation is less than the actual elevation of the station as it may imply, for example, a greater aerosol optical depth or water vapor content.

A final point is the uncertainty affecting the references $E_{ref}^h$ and $KT_{ref}^h$, which should not be neglected. This uncertainty is a combination of the uncertainty of the measurements themselves and of the construction of the hourly values with gaps.
## 4. Verification of the Extraterrestrial Irradiance Received on a Horizontal Surface

There is no measurement of the extraterrestrial irradiance at hand. This section aims at verifying the consistency over a year between the extraterrestrial irradiance received on a horizontal surface from the CAMS McClear service, $E_0^h$, and that from each reanalysis, $E_{0\_era5}^h$, $E_{0\_jra3q}^h$, and $E_{0\_merra2}^h$. Differences ($E_{0\_era5}^h-E_0^h$), ($E_{0\_jra3q}^h-E_0^h$), and ($E_{0\_merra2}^h-E_0^h$) were computed and analyzed. Null values were excluded from the analysis.

The nominal total solar irradiances of the CAMS McClear service and each reanalysis are very close: 1361 W m⁻², 1361.6 W m⁻², 1365 W m⁻², and 1365.2 W m⁻², respectively, for McClear service [18], ERA5 [10], JRA-3Q (according to ancillary metadata), and MERRA-2 [41]. This would be reflected in a small bias in the extraterrestrial irradiance. The standard deviation of the differences should be close to 0 and the correlation coefficients should be close to 1, meaning that $E_0^h$ and reanalyses vary similarly in time with similar amplitudes of the variations. As the extraterrestrial irradiance received on a horizontal surface depends on the latitude, the bias and standard deviations of differences depend on latitude too and exhibit symmetry in latitude, i.e., they do not depend on the hemisphere. They are greatest in the tropical belt and tend to decrease as the latitude increases northward or decreases southward. There is very little variation with the year. Table 2 provides results at latitudes 0° and 45° and year 2017 for illustration. Although the actual numbers depend on latitude, conclusions do not.

<table><caption>**Table 2.** Comparisons between the hourly extraterrestrial irradiance received on a horizontal surface from the CAMS McClear service, $E_0^h$, and that from each reanalysis, $E_{0\_era5}^h$, $E_{0\_jra3q}^h$, and $E_{0\_merra2}^h$ at two latitudes 0° and 45° over the year 2017. Relative values relate to the mean of $E_0^h$.</caption><thead><tr><th colspan="4">Latitude 0°</th><th colspan="4">Latitude 45°</th></tr><tr><th></th><th>Nb Values</th><th>Correl. Coeff</th><th>Mean or Bias W m<sup>-2</sup> (rel. %)</th><th>Stand. Dev. W m<sup>-2</sup> (rel. %)</th><th>Nb Values</th><th>Correl. Coeff</th><th>Mean or Bias W m<sup>-2</sup> (rel. %)</th><th>Stand. Dev. W m<sup>-2</sup> (rel. %)</th></tr></thead><tbody><tr><td>E<sup>h</sup><sub>0</sub></td><td>4729</td><td>-</td><td>770</td><td>440</td><td>4764</td><td>-</td><td>565</td><td>364</td></tr><tr><td>E<sup>h</sup><sub>0\_era5</sub></td><td>4733</td><td>-</td><td>769</td><td>440</td><td>4763</td><td>-</td><td>565</td><td>364</td></tr><tr><td>E<sup>h</sup><sub>0\_jra3q</sub></td><td>4627</td><td>-</td><td>789</td><td>431</td><td>4731</td><td>-</td><td>570</td><td>363</td></tr><tr><td>E<sup>h</sup><sub>0\_merra2</sub></td><td>5110</td><td>-</td><td>712</td><td>469</td><td>4775</td><td>-</td><td>563</td><td>364</td></tr><tr><td>(E<sup>h</sup><sub>0\_era5</sub> - E<sup>h</sup><sub>0</sub>)</td><td>4723</td><td>1.000</td><td>-1 (0)</td><td>0 (0)</td><td>4762</td><td>1.000</td><td>0 (0)</td><td>0 (0)</td></tr><tr><td>(E<sup>h</sup><sub>0\_jra3q</sub> - E<sup>h</sup><sub>0</sub>)</td><td>4627</td><td>1.000</td><td>19 (2)</td><td>-9 (-2)</td><td>4731</td><td>1.000</td><td>5 (1)</td><td>-1 (0)</td></tr><tr><td>(E<sup>h</sup><sub>0\_merra2</sub> - E<sup>h</sup><sub>0</sub>)</td><td>4729</td><td>0.997</td><td>-58 (-7)</td><td>29 (6)</td><td>4733</td><td>0.998</td><td>-2 (0)</td><td>0 (0)</td></tr></tbody></table>

The values for ERA5 and JRA-3Q are close to expectations. It is concluded that $E_0^h$, $E_{0\_era5}^h$, and $E_{0\_jra3q}^h$ are very similar, except for the nominal total solar irradiance. Note that this difference in the nominal total solar irradiance does not affect the computation of the clearness index (Equations (2) and (3)).

As for MERRA-2, the bias is close to 0 and the correlation coefficient is close to 1, although a little less. However, the standard deviations of the differences are much greater than expected: 34 W m⁻² and 21 W m⁻² (Table 2), instead of 0 W m⁻². To further investigate this discrepancy between $E_{0\_merra2}^h$ and $E_0^h$, Figure 7 displays the differences ($E_{0\_merra2}^h-E_0^h$) at 09:00, 12:00, and 15:00 UT at longitude 0° and two latitudes, 0° (top) and 45° (bottom), across the year. During a given day, the differences ($E_{0\_merra2}^h-E_0^h$) can vary from one hour to the next in a non-negligible way. They are close to 0 W m⁻² at 12:00 UT and increase in absolute value with $\theta_s$ for a given day, although they decrease as latitudes increase or decrease towards the poles. They also vary throughout the year for a given hour.
<figure><img src="image_9.png" /><figcaption>**Figure 7.** Yearly profiles (year 2017) of the differences ($E_{0\_merra2}^h - E_0^h$), labelled as (MERRA2–SG2), at 09:00, 12:00, and 15:00 UT at longitude 0° and two latitudes, 0° (top) and 45° (bottom), on which are superimposed the differences ($E_{0\_MST}^h - E_{0\_TST}^h$) (dotted lines), labelled as (E0_MST–E0_TST). Also drawn are the differences (MST–TST) in dashed line.</figcaption></figure>

Such variations are not expected and are reminiscent of those between the mean solar time (MST) and the true solar time (TST), which are drawn with the secondary axis in Figure 7, by means of equations taken from [9]. It is visible that ($E_{0\_merra2}^h - E_0^h$) vary over the year in correlation or anticorrelation with the difference (MST–TST). The MST is defined in such a way that it is equal to 12:00 on average over the year when the sun reaches its highest during a day. It is equivalent to the UT time at longitude 0°. Due to slight changes in the angular speed of the Earth along its elliptical orbit, the sun does not reach its highest at 12:00 MST every day. The peak time defines the TST, which may differ from 12:00 MST by up to ±15 min approximately. This difference is often termed “equation of time” though
it is not an equation in the usual sense. It depends on the eccentricity and obliquity of the orbit. Using the MST instead of the TST has an impact on the calculation of the solar zenith angle $\theta_s$ as the latter depends on the TST, and, therefore, it impacts the calculation of the extraterrestrial irradiance, and further the SIS. $E_0^h, E_{0\_era5}^h$, and $E_{0\_jra3q}^h$ are calculated with the TST and do not exhibit such behavior.

Despite their efforts, the present authors found no explicit mention of MST or TST in the calculation of radiation in MERRA-2 in publicly available documentation. Trolliet et al. [16] have suspected that the MST is used in MERRA-2 instead of the TST for computing $\theta_s$ and $E_{0\_merra2}^h$. This hypothesis, already partially confirmed by Figure 7, is studied in more detail in the following. The authors acknowledge that this development is somewhat lengthy but emphasize the need to understand the origin of the error in MERRA-2, which deeply influences the MERRA-2 SIS under clear-sky conditions and, moreover, under all sky conditions at intraday timescales.

The hourly extraterrestrial irradiance in MST $E_{0\_MST}^h$ and in TST $E_{0\_TST}^h$ were calculated over a year using equations in [9]. These equations are not very accurate but were used here for illustration purposes only. The differences ($E_{0\_MST}^h - E_{0\_TST}^h$) were superimposed onto those ($E_{0\_merra2}^h - E_0^h$) in Figure 7. It is clear that, at a given hour and latitude, the differences are similar. They have the same shape and the same amplitude; their averages differ by $-1 \text{ W m}^{-2}$ at both latitudes. The standard deviations of their differences are $2 \text{ W m}^{-2}$ at $0^\circ$ and $3 \text{ W m}^{-2}$ at $45^\circ$ and the correlation coefficient is equal to 0.997 and 0.988, respectively. This gives credence to the use of the MST instead of the TST in MERRA-2 and to a miscalculation of the solar zenith angle, which implies subsequent errors in $E_{0\_merra2}^h$ and SIS $E_{merra2}^h$.

A cross-correlation analysis between $E_0^h$ and $E_{0\_era5}^h, E_{0\_jra3q}^h$, and $E_{0\_merra2}^h$, respectively, was conducted to further verify that the lower correlation and higher bias observed for MERRA-2 is related to the omission of the equation of time. Indeed, if the time reference is correct, a maximum correlation will be obtained for a time lag equal to zero and only for this time lag. Otherwise, if the equation of time is omitted, the error in time is diagnosable by a maximum correlation for a non-zero lag evolving over the year. $E_0^h, E_{0\_era5}^h, E_{0\_jra3q}^h$ and $E_{0\_merra2}^h$ over the year 2021 at latitude $0^\circ$ and longitude $0^\circ$ were collected from the McClear service and the three reanalyses. They were resampled at 1 min resolution by linearly interpolating between two adjacent hourly values, yielding $E_0^{\min\_resampled}, E_{0\_era5}^{\min}, E_{0\_jra3q}^{\min}$, and $E_{0\_merra2}^{\min}$. More advanced interpolation methods would require accurate knowledge of time reference. The cross-correlation, i.e., the correlation between $E_0^{\min\_resampled}(t)$ and temporally lagged $E_{0\_merra2}^{\min}(t - \tau), E_{0\_era5}^{\min}(t - \tau)$, and $E_{0\_jra3q}^{\min}(t - \tau)$, has been calculated on a moving window of five days and the temporal lag $\tau$ ranges from -30 to 30 min. In contrast to the standard cross-correlation analysis, our approach, relying on a moving window, captures the temporal evolution of the results.

Figure 8 exhibits the calculated cross-correlation as a function of time (horizontal axis) and time lag (vertical axis), where yellow colours indicate correlations less than or equal to 0.995, while blue colours indicate correlations close to unity. Maximal cross-correlation is obtained for a zero lag throughout the year for ERA5 and JRA-3Q, indicating that the temporal reference is correct. In contrast, the maximal cross-correlation of MERRA-2 is observed for a time lag different from zero and varying over the year 2021. The shape of the domain corresponding to the maximal cross-correlation is very similar to the yearly evolution of the equation of time (Figure 7), supporting the hypothesis that the equation of time is missing in MERRA-2.
<figure><img src="image_10.png" /><figcaption>**Figure 8.** Cross-correlation between $E_0^{min}$ and $E_{0\_era5}^{min}$ (upper panel), $E_{0\_jra3q}^{min}$ (middle panel), and $E_{0\_merra2}^{min}$ (lower panel) at latitude 0° and longitude 0° as a function of time (horizontal axis) and time lag (vertical axis). The correlation is calculated on a five-day moving window. It is expected that the maximal correlation is obtained for a zero lag (black line).</figcaption></figure>

**5. Results for Clear-Sky Irradiance from ERA5 and Discussion**

Figure 9 exhibits the 2D histograms between $E_{ref}^h$ and $E_{clear\_era5}^h$ for each trimester and all stations merged. A good overall fit is observed between the estimates and the reference: the points are well aligned along the 1:1 line with limited dispersion, except during SON2017. The correlation coefficients are greater than 0.979 at any station, for any trimester, and the slope of the fitting line ranges between 0.93 and 1.08 (as detailed in the Supplementary Materials). One may note in Figure 9 that there is a tendency to underestimation of $E_{ref}^h$ by $E_{clear\_era5}^h$: the bias for all data merged is −11 W m⁻².

<figure><img src="image_11.png" /><figcaption>**Figure 9.** ERA5: 2D histograms of $E_{ref}^h$ (horizontal axis) and $E_{clear\_era5}^h$ (vertical axis) for each trimester, all stations merged. The number of elements in a class increases from yellow to deep blue.</figcaption></figure>
Figure 10 exhibits the bias and its relative value with respect to the reference means in $E^h$ (Figure 5) at each station and all stations merged for each trimester and all trimesters merged. The bias ranges between −42 W m⁻² and 27 W m⁻². It is most often negative: among the 333 cases, 70% underestimate with a bias less than −5 W m⁻², 17% have a bias in the interval [−5, 5] W m⁻², and 13% a bias greater than 5 W m⁻². Thirteen stations have a bias per trimester always less than −5 W m⁻²: TOR, LIN, CAB, PAL, PAY, CNR, LRC, MNM, ASP, FLO, DAA, GVN, and DOM. The largest negative biases for all trimesters merged are found at LIN (−28 W m⁻²), CAB, and ASP (both −25 W m⁻²). The bias has low variations from year to year or from season to season, except at FPE, SAP, and ASP. Though there are large variations in bias among stations, there is no visible relationship between the bias in $E^h$ and the latitude, or with the effective solar zenith angle $\theta_S^{\text{eff-h}}$. Notable changes in bias are observed with Köppen–Geiger climatic classes (Table 1), but no trend is visible.

<figure><img src="image_12.png" /><figcaption>**Figure 10.** ERA5: bias ($E_{\text{clear\_era5}}^h - E_{\text{ref}}^h$) and its relative value at each station and all stations merged for each trimester and all trimesters merged. The text color of the y-axis tick labels corresponds to station latitude, according to the convention described in the main text.</figcaption></figure>

Figure 11 exhibits the standard deviation of the differences ($E_{\text{clear\_era5}}^h - E_{\text{ref}}^h$) and its relative value with respect to the reference means in $E^h$ (Figure 5) at each station and all stations merged for each trimester and all trimesters merged. The standard deviation per trimester ranges between 2 W m⁻² and 44 W m⁻². In 57% of the 333 cases, the standard deviation is less than 10 W m⁻², and less than 20 W m⁻² in 94%. The standard deviations at the two southernmost sites GVN and DOM are always less than 10 W m⁻², while those at NYA, TOR, CAB, PAL, PAY, CAR, MNM, PTR, FLO, LAU, and SYO are always less than 15 W m⁻². Standard deviations per trimester greater than 25 W m⁻² are found at SAP, BOS, and GOB. In relative value, the standard deviation is less than 2% in 65% of the cases, and less than 3% in 90% of the cases. During a given season, the standard deviation has low variations from year to year, except at SAP and GOB. The contrast between seasons is weak in the southern hemisphere and more marked in the northern hemisphere, where values are among the greatest during MAM and JJA.

Figure 12 shows the 2D histograms between $KT_{\text{ref}}^h$ and $KT_{\text{clear\_era5}}^h$ for each trimester and all stations merged. As expected, the fit along the 1:1 line is not as good as for $E_{\text{clear\_era5}}^h$ but still it is quite good and the scatter is limited, except during SON2017, for which the bias is 0.01 and the standard deviation of differences is 0.03. The correlation coefficient
between $KT_{\text{clear\_era5}}^h$ and $KT_{\text{ref}}^h$ per trimester is greater than 0.85 in most of the cases, with several lower values at SAP and DOM (see Supplementary Materials). One may note in Figure 12 that there is a tendency to the underestimation of $KT_{\text{ref}}^h$ by $KT_{\text{clear\_era5}}^h$: the bias for all data merged is -0.016. Only 17% of the 333 cases exhibit a positive bias per trimester (see Supplementary Materials). The features of the biases and standard deviations per trimester per station for $KT^h$ are similar to those for $E^h$ (Figures 10 and 11).

<figure><img src="image_13.png" /><figcaption>**Figure 11.** ERA5: Standard deviation of the differences ($E_{\text{clear\_era5}}^h - E_{\text{ref}}^h$) and its relative value at each station and all stations merged for each trimester and all trimesters merged.</figcaption></figure>

<figure><img src="image_14.png" /><figcaption>**Figure 12.** ERA5: scatterplots of $KT_{\text{ref}}^h$ (horizontal axis) and $KT_{\text{clear\_era5}}^h$ (vertical axis) for each trimester, all stations merged. The number of elements in a class increases from yellow to deep blue.</figcaption></figure>
Such 2D histograms were drawn for each station and offer a great number of similarities (see Supplementary Materials). Figure 13 exhibits the 2D histograms between $KT_{ref}^h$ and $KT_{clear\_era5}^h$ at Gobabeb for each trimester as an example. Of particular interest here is the 2D histograms for JJA and SON. One may observe sorts of yellow horizontal lines that indicate that while $KT_{ref}^h$ varies, $KT_{clear\_era5}^h$ remains constant. This could reflect the influence of using a climatology that cannot reproduce the actual fluctuations in aerosol properties during this trimester, resulting in fluctuations in $KT_{ref}^h$, which are not captured by ERA5.

<figure><img src="image_15.png" /><figcaption>**Figure 13.** ERA5: 2D histograms of $KT_{ref}^h$ (horizontal axis) and $KT_{clear\_era5}^h$ (vertical axis) for each trimester at Gobabeb. The number of elements in a class increases from yellow to deep blue.</figcaption></figure>

Figure 14 exhibits boxplots of the differences ($KT_{clear\_era5}^h - KT_{ref}^h$) for all stations merged for several classes of $\theta_S^{eff\_h}$ for each trimester. The bias in $KT^h$ is always negative in this graph. The bias and the standard deviation of the differences vary with $\theta_S^{eff\_h}$ and the trimester. The bias is more and more negative and the standard deviation increases as $\theta_S^{eff\_h}$ increases. This effect is most pronounced during JJA. Changes from year to year are limited and less marked than those from season to season.

No trend with climate can be detected visually on the bias and standard deviation. Whatever the climate, the differences are most often negative; the standard deviation remains fairly constant, and the bias exhibits more fluctuations. Boxplots for each climate, each season, and each year for the $\theta_S^{eff\_h}$ classes exhibit complex patterns with no obvious trend. A thorough analysis is necessary and is beyond our investigative means.
Köppen-Geiger class: all climates (28 stations) / reanalysis: era5

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

**Figure 14.** ERA5: boxplots of the differences ($KT_{clear\_era5}^h - KT_{ref}^h$) for all stations merged for several classes of $\theta_S^{eff\_h}$ for each trimester (SON, MAM, JJA and SON). Percentages on the right indicate the fraction of cases in each class.

Table 3 summarizes the bias, standard deviation of differences, and root mean square of differences (RMSD) in $E^h$ and $KT^h$ at each station and all stations merged for all trimesters merged. The relative bias in $E^h$ and $KT^h$, respectively, is most often negative and ranges from -5% to 1% and -6% to 1%. The relative standard deviation in $E^h$ or $KT^h$ has low variations as it is most often in the range [2, 3%].
<table><caption>**Table 3.** ERA5: Statistical indicators in $E^h$ and $KT^h$ at each station and all stations merged for all trimesters merged. "Stand. dev." stands for standard deviation of differences and RMSD stands for root mean square of differences.</caption><thead><tr><th rowspan="2">Station Code</th><th colspan="3">E<sup>h</sup></th><th colspan="3">KT<sup>h</sup></th></tr><tr><th>Bias W m<sup>-2</sup><br>(rel. %)</th><th>Stand. Dev. W<br>m<sup>-2</sup> (rel. %)</th><th>RMSD W m<sup>-2</sup><br>(rel. %)</th><th>Bias (rel. %)</th><th>Stand. Dev.<br>(rel. %)</th><th>RMSD (rel. %)</th></tr></thead><tbody><tr><td>NYA</td><td>-4 (-1)</td><td>11 (3)</td><td>12 (3)</td><td>-0.01 (-1)</td><td>0.02 (3)</td><td>0.02 (3)</td></tr><tr><td>TOR</td><td>-21 (-4)</td><td>10 (2)</td><td>23 (5)</td><td>-0.03 (-5)</td><td>0.02 (3)</td><td>0.04 (6)</td></tr><tr><td>LIN</td><td>-27 (-5)</td><td>13 (2)</td><td>30 (6)</td><td>-0.04 (-6)</td><td>0.02 (3)</td><td>0.05 (7)</td></tr><tr><td>CAB</td><td>-25 (-5)</td><td>12 (2)</td><td>28 (5)</td><td>-0.04 (-5)</td><td>0.02 (3)</td><td>0.04 (6)</td></tr><tr><td>PAL</td><td>-19 (-4)</td><td>12 (2)</td><td>22 (4)</td><td>-0.03 (-4)</td><td>0.02 (3)</td><td>0.03 (5)</td></tr><tr><td>FPE</td><td>7 (1)</td><td>21 (3)</td><td>22 (4)</td><td>0.01 (1)</td><td>0.03 (4)</td><td>0.03 (4)</td></tr><tr><td>PAY</td><td>-20 (-3)</td><td>12 (2)</td><td>23 (4)</td><td>-0.03 (-4)</td><td>0.02 (3)</td><td>0.03 (5)</td></tr><tr><td>CAR</td><td>-10 (-2)</td><td>12 (2)</td><td>16 (3)</td><td>-0.02 (-2)</td><td>0.02 (3)</td><td>0.03 (4)</td></tr><tr><td>SAP</td><td>-9 (-2)</td><td>29 (5)</td><td>31 (5)</td><td>-0.01 (-2)</td><td>0.04 (6)</td><td>0.04 (6)</td></tr><tr><td>CNR</td><td>-14 (-2)</td><td>13 (2)</td><td>19 (3)</td><td>-0.02 (-3)</td><td>0.02 (3)</td><td>0.03 (4)</td></tr><tr><td>BOS</td><td>2 (0)</td><td>15 (3)</td><td>16 (3)</td><td>0.00 (1)</td><td>0.02 (3)</td><td>0.02 (3)</td></tr><tr><td>LRC</td><td>-16 (-3)</td><td>9 (2)</td><td>18 (3)</td><td>-0.02 (-3)</td><td>0.02 (2)</td><td>0.03 (4)</td></tr><tr><td>DRA</td><td>-5 (-1)</td><td>17 (3)</td><td>18 (3)</td><td>0.00 (0)</td><td>0.02 (3)</td><td>0.02 (3)</td></tr><tr><td>BIL</td><td>5 (1)</td><td>13 (2)</td><td>13 (2)</td><td>0.01 (1)</td><td>0.02 (3)</td><td>0.02 (3)</td></tr><tr><td>E13</td><td>-2 (0)</td><td>15 (3)</td><td>15 (3)</td><td>0.00 (0)</td><td>0.02 (2)</td><td>0.02 (2)</td></tr><tr><td>TAT</td><td>-12 (-2)</td><td>14 (3)</td><td>18 (4)</td><td>-0.02 (-3)</td><td>0.02 (3)</td><td>0.03 (4)</td></tr><tr><td>FUA</td><td>-15 (-3)</td><td>20 (3)</td><td>25 (4)</td><td>-0.02 (-3)</td><td>0.03 (4)</td><td>0.04 (5)</td></tr><tr><td>MNM</td><td>-16 (-2)</td><td>9 (1)</td><td>18 (3)</td><td>-0.02 (-3)</td><td>0.02 (3)</td><td>0.03 (4)</td></tr><tr><td>PTR</td><td>4 (1)</td><td>16 (2)</td><td>16 (3)</td><td>0.00 (0)</td><td>0.02 (3)</td><td>0.02 (3)</td></tr><tr><td>DWN</td><td>-4 (-1)</td><td>15 (3)</td><td>15 (3)</td><td>-0.01 (-1)</td><td>0.02 (3)</td><td>0.02 (3)</td></tr><tr><td>GOB</td><td>-7 (-1)</td><td>18 (3)</td><td>19 (3)</td><td>-0.01 (-1)</td><td>0.03 (3)</td><td>0.03 (4)</td></tr><tr><td>ASP</td><td>-25 (-4)</td><td>13 (2)</td><td>28 (4)</td><td>-0.03 (-4)</td><td>0.02 (3)</td><td>0.04 (5)</td></tr><tr><td>FLO</td><td>-20 (-3)</td><td>11 (2)</td><td>23 (4)</td><td>-0.03 (-4)</td><td>0.02 (3)</td><td>0.03 (5)</td></tr><tr><td>DAA</td><td>-19 (-3)</td><td>11 (2)</td><td>22 (3)</td><td>-0.03 (-3)</td><td>0.02 (2)</td><td>0.03 (4)</td></tr><tr><td>LAU</td><td>1 (0)</td><td>10 (2)</td><td>10 (2)</td><td>0.00 (0)</td><td>0.02 (2)</td><td>0.02 (2)</td></tr><tr><td>SYO</td><td>1 (0)</td><td>8 (2)</td><td>9 (2)</td><td>0.00 (0)</td><td>0.01 (2)</td><td>0.01 (2)</td></tr><tr><td>GVN</td><td>-21 (-4)</td><td>8 (2)</td><td>22 (4)</td><td>-0.03 (-4)</td><td>0.01 (1)</td><td>0.04 (4)</td></tr><tr><td>DOM</td><td>-14 (-3)</td><td>7 (1)</td><td>16 (3)</td><td>-0.02 (-3)</td><td>0.01 (2)</td><td>0.03 (3)</td></tr><tr><td>All stations</td><td>-11 (-2)</td><td>17 (3)</td><td>20 (3)</td><td>-0.02 (-2)</td><td>0.03 (3)</td><td>0.03 (4)</td></tr></tbody></table>

## 6. Results for Clear-Sky Irradiance from JRA-3Q and Discussion

Figure 15 exhibits the 2D histograms between $E^h_{ref}$ and $E^h_{clear\_jra3q}$ for each trimester and all stations merged. A good overall fit is observed between the estimates and the reference: the points are well aligned along the 1:1 line with limited dispersion. The correlation coefficients between $E^h_{clear\_jra3q}$ and $E^h_{ref}$ are greater than 0.980 at any station for any trimester, and the slope of the fitting line ranges between 0.95 and 1.09, with a mean equal to 1.01 (see Supplementary Materials).
<figure><img src="image_17.png" /><figcaption>**Figure 15.** JRA-3Q: 2D histograms of $E_{ref}^h$ (horizontal axis) and $E_{clear\_jra3q}^h$ (vertical axis) for each trimester, all stations merged. The number of elements in a class increases from yellow to deep blue.</figcaption></figure>

Figure 16 exhibits the bias and its relative value with respect to the reference means in $E^h$ (Figure 5) at each station and all stations merged for each trimester and all trimesters merged. The bias at Desert Rock (DRA) is always negative and very low, being between -13 W m⁻² and -63 W m⁻². Otherwise, the bias ranges between -36 W m⁻² and 33 W m⁻². The bias is in the range [-5, 5] W m⁻² in 114 cases (34%), in the range [5, 33] W m⁻² in 84 cases (25%), and less than -5 W m⁻² in 135 cases (41%). BIL, E13, and PTR exhibit the greatest biases. At a given station and during a given trimester, the bias has low variations from year to year except at the stations in Northern America (FPE, LRC, DRA, BIL, and E13) and TAT. The biases exhibit large variations from trimester to trimester, and from station to station. There is no visible relationship between the bias in $E^h$ and the latitude, or with the effective solar zenith angle $\theta_S^{eff\_h}$. Significant changes with Köppen-Geiger climatic classes (Table 1) are observed but no trend has been identified.

Figure 17 exhibits the standard deviation of the differences ($E_{clear\_jra3q}^h - E_{ref}^h$) and its relative value with respect to the reference means in $E^h$ (Figure 5) at each station and all stations merged for each trimester and all trimesters merged. The standard deviation per trimester ranges between 3 W m⁻² and 46 W m⁻². In 54% of the 333 cases, the standard deviation is less than 10 W m⁻², and less than 20 W m⁻² in 96% of cases. MNM and the
two southernmost sites GVN and DOM have standard deviations always less than 10 W m⁻². Those at NYA, TOR, LIN, CAB, PAL, PAY, CAR, PTR, FLO, and LAU, are always less than 15 W m⁻². Values greater than 25 W m⁻² are found at SAP and GOB. In relative value, the standard deviation is less than 2% in 75% of the cases and less than 3% in 96% of the cases. During a given season, the standard deviation has low variations from year to year, except at SAP and GOB. The contrast between seasons is weak in the southern hemisphere and more marked in the northern hemisphere, where values are among the greatest during MAM and JJA.

<figure><img src="image_18.png" /><figcaption>**Figure 16.** JRA-3Q: bias ($E_{clear\_jra3q}^h$ – $E_{ref}^h$) and its relative value at each station and all stations merged for each trimester and all trimesters merged.</figcaption></figure>

<figure><img src="image_19.png" /><figcaption>**Figure 17.** JRA-3Q: standard deviation of the differences ($E_{clear\_jra3q}^h$ – $E_{ref}^h$) and its relative value at each station and all stations merged for each trimester and all trimesters merged.</figcaption></figure>

Figure 18 shows the 2D histograms between $KT_{ref}^h$ and $KT_{clear\_jra3q}^h$ for each trimester and all stations merged. The fit along the 1:1 line is quite good though not as good as for $E_{clear\_jra3q}^h$. The scatter varies from trimester to trimester and is systematically high during
JJA. The correlation coefficient between $KT_{clear\_jra3q}^h$ and $KT_{ref}^h$ per trimester is greater than 0.85 in most of the cases. Very low values are found at CNR during DJF2017 (0.672) and DOM during MAM2017 (0.481) (see Supplementary Materials). The slope per trimester (see Supplementary Materials) is comprised between 0.95 and 1.05 in 24% of the cases. It is greater than 1.05 in 33% of the cases and less than 0.95 in 43% of the cases.

<figure><img src="image_20.png" /><figcaption>**Figure 18.** JRA-3Q: scatterplots of $KT_{ref}^h$ (horizontal axis) and $KT_{clear\_jra3q}^h$ (vertical axis) for each trimester, all stations merged. The number of elements in a class increases from yellow to deep blue.</figcaption></figure>

The bias in $KT^h$ per trimester (see Supplementary Materials) is close to 0 overall: 41% of the cases are in the range [-0.01, 0.01], 18% are greater than 0.01, 18% in the range [-0.02, -0.01], and 23% are less than -0.02. The bias is always less than -0.01 at CAB, DRA, and SYO and always greater than 0.01 at BIL and PTR. It varies from station to station, trimester to trimester, and year to year. The features of the biases and standard deviations per trimester per station for $KT^h$ are similar to those for $E^h$ (Figures 16 and 17).

Two-dimensional histograms were drawn at each station and offer a great number of similarities (see Supplementary Materials). Figure 19 exhibits the 2D histograms between $KT_{ref}^h$ and $KT_{clear\_jra3q}^h$ at Gobabeb for each trimester as an example. During JJA and SON, one may observe sorts of yellow quasi-horizontal lines (see e.g., JJA2015 for $KT_{clear\_jra3q}^h = 0.6$) that indicate that while $KT_{ref}^h$ varies, $KT_{clear\_jra3q}^h$ remains constant.
This could reflect the influence of using a climatology for aerosols, which cannot reproduce the actual fluctuations in aerosol optical properties during this trimester, resulting in fluctuations in $KT_{ref}^h$, which are not captured by JRA-3Q.

Figure 20 exhibits boxplots of the differences ($KT_{clear\_jra3q}^h - KT_{ref}^h$) for all stations merged for several classes of $\theta_S^{eff\_h}$ for each trimester. During a given year and a given season, changes in bias with $\theta_S^{eff\_h}$ are limited while the standard deviation increases as $\theta_S^{eff\_h}$ increases. During DJF and MAM, differences are most often negative whatever the value of $\theta_S^{eff\_h}$, meaning an underestimation of $KT_{ref}^h$ by $KT_{clear\_jra3q}^h$. Changes from year to year are limited during these seasons. During JJA and SON, differences are more centered on 0 for all years but 2016, which tends to be more negative. Other graphs show that the bias and standard deviation for each class of $\theta_S^{eff\_h}$ fluctuate between stations; no trend with latitude or climate can be detected visually.

<figure><img src="image_21.png" /><figcaption>**Figure 19.** JRA-3Q: 2D histograms of $KT_{ref}^h$ (horizontal axis) and $KT_{clear\_jra3q}^h$ (vertical axis) for each trimester at Gobabeb. The number of elements in a class increases from yellow to deep blue.</figcaption></figure>

Boxplots for each climate, each season, and each year for the $\theta_S^{eff\_h}$ classes exhibit complex patterns with no obvious trend. For each possible class, the differences may be low or high, negative or positive. Biases fluctuate significantly; standard deviations tend to increase as $\theta_S^{eff\_h}$ increases.
Table 4 summarizes the bias, standard deviation of differences, and root mean square of differences (RMSD) in $E^h$ and $KT^h$ at each station and all stations merged for all trimesters merged. The relative bias in $E^h$ and $KT^h$, respectively, ranges from -6% to 4% and -7% to 4%. The relative standard deviation in $E^h$ or $KT^h$ has low variations as it is often in the range [2, 3] %.

<table><caption>**Table 4.** JRA-3Q: Statistical indicators in $E^h$ and $KT^h$ at each station and all stations merged for all trimesters merged. "Stand. dev." stands for standard deviation of differences and RMSD stands for root mean square of differences.</caption><thead><tr><th rowspan="2">Station Code</th><th colspan="3">E<sup>h</sup></th><th colspan="3">KT<sup>h</sup></th></tr><tr><th>Bias W m<sup>-2</sup><br>(rel. %)</th><th>Stand. Dev. W<br>m<sup>-2</sup> (rel. %)</th><th>RMSD W m<sup>-2</sup><br>(rel. %)</th><th>Bias (rel. %)</th><th>Stand. Dev.<br>(rel. %)</th><th>RMSD (rel. %)</th></tr></thead><tbody><tr><td>NYA</td><td>-4 (-1)</td><td>12 (3)</td><td>12 (3)</td><td>-0.01 (-1)</td><td>0.02 (3)</td><td>0.03 (3)</td></tr><tr><td>TOR</td><td>-2 (0)</td><td>11 (2)</td><td>11 (2)</td><td>-0.01 (-1)</td><td>0.02 (3)</td><td>0.02 (3)</td></tr><tr><td>LIN</td><td>-9 (-2)</td><td>13 (3)</td><td>16 (3)</td><td>-0.02 (-2)</td><td>0.02 (3)</td><td>0.03 (4)</td></tr><tr><td>CAB</td><td>-14 (-3)</td><td>13 (2)</td><td>19 (4)</td><td>-0.02 (-3)</td><td>0.02 (3)</td><td>0.03 (4)</td></tr><tr><td>PAL</td><td>-11 (-2)</td><td>12 (2)</td><td>16 (3)</td><td>-0.02 (-3)</td><td>0.02 (3)</td><td>0.03 (4)</td></tr><tr><td>FPE</td><td>-7 (-1)</td><td>21 (4)</td><td>22 (4)</td><td>-0.01 (-1)</td><td>0.03 (4)</td><td>0.03 (4)</td></tr><tr><td>PAY</td><td>5 (1)</td><td>11 (2)</td><td>13 (2)</td><td>0.00 (1)</td><td>0.02 (2)</td><td>0.02 (2)</td></tr><tr><td>CAR</td><td>5 (1)</td><td>12 (2)</td><td>14 (2)</td><td>0.00 (1)</td><td>0.02 (3)</td><td>0.02 (3)</td></tr><tr><td>SAP</td><td>-4 (-1)</td><td>25 (4)</td><td>25 (4)</td><td>-0.01 (-2)</td><td>0.04 (5)</td><td>0.04 (5)</td></tr><tr><td>CNR</td><td>-4 (-1)</td><td>13 (2)</td><td>13 (2)</td><td>-0.01 (-1)</td><td>0.02 (3)</td><td>0.02 (3)</td></tr><tr><td>BOS</td><td>-1 (0)</td><td>17 (3)</td><td>17 (3)</td><td>0.00 (0)</td><td>0.02 (3)</td><td>0.02 (3)</td></tr><tr><td>LRC</td><td>-13 (-2)</td><td>15 (3)</td><td>20 (4)</td><td>-0.02 (-3)</td><td>0.02 (3)</td><td>0.03 (4)</td></tr><tr><td>DRA</td><td>-41 (-6)</td><td>24 (4)</td><td>47 (7)</td><td>-0.05 (-7)</td><td>0.02 (3)</td><td>0.06 (7)</td></tr><tr><td>BIL</td><td>20 (4)</td><td>14 (2)</td><td>24 (4)</td><td>0.03 (4)</td><td>0.02 (3)</td><td>0.03 (5)</td></tr><tr><td>E13</td><td>13 (2)</td><td>17 (3)</td><td>21 (4)</td><td>0.02 (3)</td><td>0.02 (3)</td><td>0.03 (4)</td></tr><tr><td>TAT</td><td>-15 (-3)</td><td>16 (3)</td><td>22 (4)</td><td>-0.03 (-4)</td><td>0.03 (4)</td><td>0.04 (5)</td></tr><tr><td>FUA</td><td>-10 (-2)</td><td>19 (3)</td><td>22 (4)</td><td>-0.02 (-2)</td><td>0.03 (4)</td><td>0.03 (5)</td></tr><tr><td>MNM</td><td>-9 (-1)</td><td>9 (1)</td><td>12 (2)</td><td>-0.01 (-2)</td><td>0.02 (2)</td><td>0.02 (3)</td></tr><tr><td>PTR</td><td>23 (4)</td><td>13 (2)</td><td>27 (4)</td><td>0.03 (4)</td><td>0.01 (2)</td><td>0.03 (4)</td></tr><tr><td>DWN</td><td>7 (1)</td><td>15 (3)</td><td>17 (3)</td><td>0.01 (1)</td><td>0.02 (4)</td><td>0.03 (4)</td></tr><tr><td>GOB</td><td>-17 (-2)</td><td>17 (2)</td><td>24 (3)</td><td>-0.02 (-3)</td><td>0.02 (3)</td><td>0.03 (4)</td></tr><tr><td>ASP</td><td>12 (2)</td><td>10 (2)</td><td>16 (2)</td><td>0.01 (2)</td><td>0.01 (2)</td><td>0.02 (3)</td></tr><tr><td>FLO</td><td>-2 (0)</td><td>12 (2)</td><td>12 (2)</td><td>0.00 (0)</td><td>0.02 (2)</td><td>0.02 (3)</td></tr><tr><td>DAA</td><td>-3 (0)</td><td>13 (2)</td><td>13 (2)</td><td>0.00 (-1)</td><td>0.02 (2)</td><td>0.02 (2)</td></tr><tr><td>LAU</td><td>6 (1)</td><td>10 (2)</td><td>11 (2)</td><td>0.01 (1)</td><td>0.01 (2)</td><td>0.02 (2)</td></tr><tr><td>SYO</td><td>3 (1)</td><td>9 (2)</td><td>10 (2)</td><td>0.00 (1)</td><td>0.02 (2)</td><td>0.02 (2)</td></tr><tr><td>GVN</td><td>-15 (-3)</td><td>9 (2)</td><td>17 (3)</td><td>-0.02 (-3)</td><td>0.01 (2)</td><td>0.03 (3)</td></tr><tr><td>DOM</td><td>-9 (-2)</td><td>7 (1)</td><td>11 (2)</td><td>-0.02 (-2)</td><td>0.01 (2)</td><td>0.02 (2)</td></tr><tr><td>All stations</td><td>-4 (-1)</td><td>20 (3)</td><td>20 (3)</td><td>-0.01 (-1)</td><td>0.03 (4)</td><td>0.03 (4)</td></tr></tbody></table>
Köppen-Geiger class: all climates (28 stations) / reanalysis: jra3q

<figure><img src="image_22.png" /><figcaption>**Figure 20.** JRA-3Q: boxplots of the differences ($KT_{clear\_jra3q}^h - KT_{ref}^h$) for all stations merged for several classes of $\theta_s^{eff\_h}$ for each trimester (SON, MAM, JJA and SON). Percentages on the right indicate the fraction of cases in each class.</figcaption></figure>

**7. Results for Clear-Sky Irradiance from MERRA-2 and Discussion**

The results for MERRA-2 regarding $E_s^h$ are affected by the mismatch in time system at intraday scales and, therefore, by a miscalculation of the solar zenith angle $\theta_s$. This miscalculation is propagated to the irradiance at the top of the atmosphere, $E_{0\_merra2}^h$, and also to the clear-sky SIS $E_{clear\_merra2}^h$ and clearness index $KT_{clear\_merra2}^h$ due to the dependence of absorption and scattering by atmospheric constituents on $\theta_s$. Figure 21 exhibits the 2D histograms between $E_{ref}^h$ and $E_{clear\_merra2}^h$ for each trimester and all stations merged. Overall, the points are well aligned along the 1:1 line. The slope of the fitting line per trimester per station ranges between 0.91 and 1.11 but ranges between 0.95 and 1.05 in most cases (see Supplementary Materials). The dispersion of the most populated
classes (green and blue) along the 1:1 line is noticeable during DJF and SON and is limited
during MAM and JJA. This is only a general feature when all stations are merged and the
actual situation varies from station to station. Similar 2D histograms at each station are
given in the Supplementary Materials. As an example, Figure 22 exhibits the 2D histograms
between $E_{ref}^h$ and $E_{clear\_merra2}^h$ for each trimester at Gobabeb. One may observe during DJF
that the dots are not well distributed along the 1:1 line and form horseshoe-shaped spots
along this line that are typical of the time system mismatch.

<figure><img src="image_23.png" /><figcaption>**Figure 21.** MERRA-2: 2D histograms of $E_{ref}^h$ (horizontal axis) and $E_{clear\_merra2}^h$ (vertical axis) for each trimester, all stations merged. The number of elements in a class increases from yellow to deep blue.</figcaption></figure>
<figure><img src="image_24.png" /><figcaption>**Figure 22.** MERRA-2: 2D histograms of $E_{ref}^h$ (horizontal axis) and $E_{clear\_merra2}^h$ (vertical axis) for each trimester at Gobabeb (GOB). The number of elements in a class increases from yellow to deep blue.</figcaption></figure>

Figure 23 displays the correlation coefficients per trimester per station. They are greater than 0.967 but many of them (14%) are less than 0.985. The coefficients exhibit a pattern that has seasonal and latitudinal dependency. The smallest are observed in the northern hemisphere during DJF and SON to a lesser extent and during MAM at the South Pole. This pattern in time and space originates from the mismatch in time system and can affect intraday calculations.
<figure><img src="image_25.png" /><figcaption>**Figure 23.** MERRA-2: correlation coefficient between $E_{clear\_merra2}^h$ and $E_{ref}^h$ at each station and all stations merged for each trimester and all trimesters merged.</figcaption></figure>

The bias and its relative value with respect to the reference means in $E^h$ (Figure 5) at each station and all stations merged for each trimester and all trimesters merged are displayed in Figure 24. The bias ranges between −39 W m⁻² and 34 W m⁻². It is in the interval [−5, 5] W m⁻² in 35% of the cases; 16% have a bias less than −5 W m⁻², 49% a bias greater than 5 W m⁻², and 17% a bias greater than 15 W m⁻². As a whole, the bias at a given station during a given trimester does not vary much with the year. It varies more from trimester to trimester and from station to station. The two southernmost stations have large negative biases, less than −15 W m⁻², during all trimesters. Except for these stations, there is no clear relationship between the bias in $E^h$ and the latitude, though one may note that the stations between 49° and 58° (TOR, LIN, CAB, and PAL) and between −23° and −45° (GOB, ASP, FLO, DAA, and LAU) have low relative biases, mostly comprised between −2% and 2%, with no clear link with Köppen-Geiger climatic classes (Table 1).
<figure><img src="image_26.png" /><figcaption>**Figure 24.** MERRA-2: bias ($E_{clear\_merra2}^h - E_{ref}^h$) and its relative value at each station and all stations merged for each trimester and all trimesters merged.</figcaption></figure>

Figure 25 exhibits the standard deviation of the differences ($E_{clear\_merra2}^h - E_{ref}^h$) and its relative value with respect to the reference means in $E_{ref}^h$ (Figure 5) at each station and all stations merged for each trimester and all trimesters merged. The standard deviation per trimester ranges between 5 W m⁻² and 43 W m⁻². It is large as a whole: it is less than 10 W m⁻² in only 5% of the 333 cases. It lies in the range [10, 20] W m⁻² in 57% of the cases, in the range [20, 30] W m⁻² in 24% of the cases, and is greater than 30 W m⁻² in 14% of the cases. In relative value, the standard deviation is less than 2% in only 9% of the cases. It is in the range [2, 3] % in 34% of the cases, in the range [3, 4] % in 20% of the cases, and greater than 4% in 23% of the cases. There is no clear picture between the standard deviation and latitude or station or climate. What is clear is the increase in the standard deviation during DJF and SON compared to MAM and JJA. This is to be related to the similar pattern discussed in Figure 23, which is attributed to the mismatch in time system.

<figure><img src="image_27.png" /><figcaption>**Figure 25.** MERRA-2: standard deviation of the differences ($E_{clear\_merra2}^h - E_{ref}^h$) and its relative value at each station and all stations merged for each trimester and all trimesters merged.</figcaption></figure>
Figure 26 shows the 2D histograms between $KT_{ref}^h$ and $KT_{clear\_merra2}^h$ for each trimester and all stations merged. As expected, the fit along the 1:1 line is not as good as for $E_{clear\_merra2}^h$ but it is still quite good and the scatter is limited. One may note during DJF, JJA, and SON, a group of dots of high $KT_{ref}^h$ that lies well below the 1:1: line. This group is due to the southernmost stations: SYO, GVN, and DOM, for which the bias is largely negative, being, respectively, -0.02, -0.05, and -0.06 when averaged over all trimesters. The northernmost station NYA, also in a polar climate, also has a negative bias of -0.01 on average. These four stations have the most noticeable negative biases averaged over all trimesters; the other stations exhibit biases close to 0 or with positive values of up to 0.02 (Table 5), as illustrated in Figure 26. It may be hypothesized that these negative biases are partly due to incomplete consideration of the high albedo surfaces surrounding the sites, though the authors do not have the details of the calculation of $KT_{clear\_merra2}^h$, and other reasons are possible. Among them, one can think of measurement issues due to frozen pyranometers [42] not detected by quality control, although measurement issues cannot be the major reason as they should also affect the results for ERA5 and JRA-3Q, which do not exhibit behaviours at polar stations similar to those for MERRA-2.

<figure><img src="image_28.png" /><figcaption>**Figure 26.** MERRA-2: scatterplots of $KT_{ref}^h$ (horizontal axis) and $KT_{clear\_merra2}^h$ (vertical axis) for each trimester, all stations merged. The number of elements in a class increases from yellow to deep blue.</figcaption></figure>
<table><caption>**Table 5.** MERRA-2: Statistical indicators in $E^h$ and $KT^h$ at each station and all stations merged for all trimesters merged. "Stand. dev." stands for standard deviation of differences and RMSD stands for root mean square of differences.</caption><thead><tr><th rowspan="2">Station Code</th><th colspan="3">E<sup>h</sup></th><th colspan="3">KT<sup>h</sup></th></tr><tr><th>Bias W m<sup>-2</sup><br>(rel. %)</th><th>Stand. Dev. W<br>m<sup>-2</sup> (rel. %)</th><th>RMSD W m<sup>-2</sup><br>(rel. %)</th><th>Bias (rel. %)</th><th>Stand. Dev.<br>(rel. %)</th><th>RMSD (rel. %)</th></tr></thead><tbody><tr><td>NYA</td><td>-6 (-2)</td><td>11 (3)</td><td>13 (3)</td><td>-0.01 (-2)</td><td>0.02 (3)</td><td>0.02 (3)</td></tr><tr><td>TOR</td><td>2 (1)</td><td>13 (3)</td><td>13 (3)</td><td>0.00 (1)</td><td>0.01 (2)</td><td>0.02 (2)</td></tr><tr><td>LIN</td><td>6 (1)</td><td>14 (3)</td><td>15 (3)</td><td>0.01 (1)</td><td>0.01 (2)</td><td>0.02 (2)</td></tr><tr><td>CAB</td><td>2 (0)</td><td>14 (3)</td><td>14 (3)</td><td>0.00 (1)</td><td>0.01 (2)</td><td>0.01 (2)</td></tr><tr><td>PAL</td><td>0 (0)</td><td>14 (3)</td><td>14 (3)</td><td>0.00 (0)</td><td>0.02 (2)</td><td>0.02 (2)</td></tr><tr><td>FPE</td><td>11 (2)</td><td>19 (3)</td><td>22 (4)</td><td>0.02 (2)</td><td>0.02 (3)</td><td>0.03 (4)</td></tr><tr><td>PAY</td><td>13 (2)</td><td>16 (3)</td><td>20 (4)</td><td>0.02 (2)</td><td>0.01 (2)</td><td>0.02 (3)</td></tr><tr><td>CAR</td><td>16 (3)</td><td>18 (3)</td><td>24 (4)</td><td>0.02 (3)</td><td>0.01 (2)</td><td>0.02 (3)</td></tr><tr><td>SAP</td><td>13 (2)</td><td>20 (3)</td><td>23 (4)</td><td>0.02 (2)</td><td>0.02 (3)</td><td>0.02 (4)</td></tr><tr><td>CNR</td><td>8 (1)</td><td>18 (3)</td><td>20 (3)</td><td>0.01 (1)</td><td>0.01 (2)</td><td>0.02 (2)</td></tr><tr><td>BOS</td><td>15 (3)</td><td>23 (4)</td><td>27 (5)</td><td>0.02 (3)</td><td>0.02 (3)</td><td>0.03 (4)</td></tr><tr><td>LRC</td><td>-1 (0)</td><td>22 (4)</td><td>22 (4)</td><td>0.00 (0)</td><td>0.01 (2)</td><td>0.01 (2)</td></tr><tr><td>DRA</td><td>0 (0)</td><td>25 (4)</td><td>25 (4)</td><td>0.01 (1)</td><td>0.02 (3)</td><td>0.02 (3)</td></tr><tr><td>BIL</td><td>19 (4)</td><td>24 (4)</td><td>30 (6)</td><td>0.03 (4)</td><td>0.02 (3)</td><td>0.03 (5)</td></tr><tr><td>E13</td><td>13 (2)</td><td>24 (4)</td><td>27 (5)</td><td>0.02 (3)</td><td>0.02 (3)</td><td>0.03 (4)</td></tr><tr><td>TAT</td><td>-2 (0)</td><td>24 (5)</td><td>24 (5)</td><td>0.00 (-1)</td><td>0.02 (3)</td><td>0.02 (3)</td></tr><tr><td>FUA</td><td>16 (3)</td><td>25 (4)</td><td>30 (5)</td><td>0.02 (3)</td><td>0.02 (3)</td><td>0.03 (4)</td></tr><tr><td>MNM</td><td>2 (0)</td><td>21 (3)</td><td>21 (3)</td><td>0.01 (1)</td><td>0.01 (2)</td><td>0.01 (2)</td></tr><tr><td>PTR</td><td>22 (4)</td><td>32 (5)</td><td>39 (6)</td><td>0.02 (3)</td><td>0.02 (2)</td><td>0.03 (4)</td></tr><tr><td>DWN</td><td>13 (2)</td><td>23 (4)</td><td>26 (5)</td><td>0.02 (2)</td><td>0.02 (3)</td><td>0.03 (4)</td></tr><tr><td>GOB</td><td>-4 (-1)</td><td>28 (4)</td><td>28 (4)</td><td>0.00 (-1)</td><td>0.02 (3)</td><td>0.02 (3)</td></tr><tr><td>ASP</td><td>4 (1)</td><td>28 (4)</td><td>28 (4)</td><td>0.01 (1)</td><td>0.02 (2)</td><td>0.02 (2)</td></tr><tr><td>FLO</td><td>7 (1)</td><td>22 (4)</td><td>23 (4)</td><td>0.01 (1)</td><td>0.02 (2)</td><td>0.02 (3)</td></tr><tr><td>DAA</td><td>-1 (0)</td><td>28 (4)</td><td>28 (4)</td><td>0.00 (0)</td><td>0.02 (2)</td><td>0.02 (2)</td></tr><tr><td>LAU</td><td>5 (1)</td><td>25 (4)</td><td>25 (4)</td><td>0.01 (1)</td><td>0.02 (2)</td><td>0.02 (2)</td></tr><tr><td>SYO</td><td>-14 (-3)</td><td>17 (3)</td><td>22 (4)</td><td>-0.02 (-2)</td><td>0.02 (2)</td><td>0.02 (3)</td></tr><tr><td>GVN</td><td>-32 (-6)</td><td>17 (3)</td><td>37 (7)</td><td>-0.05 (-6)</td><td>0.01 (2)</td><td>0.05 (6)</td></tr><tr><td>DOM</td><td>-35 (-7)</td><td>15 (3)</td><td>38 (7)</td><td>-0.06 (-7)</td><td>0.01 (2)</td><td>0.06 (7)</td></tr><tr><td>All stations</td><td>2 (0)</td><td>27 (4)</td><td>27 (4)</td><td>0.00 (0)</td><td>0.03 (4)</td><td>0.03 (4)</td></tr></tbody></table>

The correlation coefficient between $KT_{clear\_merra2}^h$ and $KT_{ref}^h$ per trimester is greater than 0.85 in most of the cases (see Supplementary Materials). This coefficient does not exhibit the same pattern as that between $E_{clear\_merra2}^h$ and $E_{ref}^h$ (Figure 23) as the effect of the diurnal cycle of $E_{0\_merra2}^h$ is removed in $KT_{clear\_merra2}^h$ and this mitigates the effects of the mismatch in time system.

Similar 2D histograms were drawn at each station and offer a great number of similarities (see Supplementary Materials). As an example, Figure 27 exhibits the 2D histograms between $KT_{ref}^h$ and $KT_{clear\_merra2}^h$ at Gobabeb for each trimester. Of particular interest here is the 2D histogram for SON2017, which can be compared to those obtained for ERA5 (Figure 13) and JRA-3Q (Figure 19) on which sorts of yellow quasi-horizontal lines indicate
that while $KT_{ref}^h$ varies, $KT_{clear\_era5}^h$ and $KT_{clear\_jira3q}^h$, respectively, remain constant, reflecting the influence of using a climatology. MERRA-2 does not exhibit such a pattern and captures the actual fluctuations in $KT_{ref}^h$, though the scatter cannot be neglected. This could be attributed to the assimilation of aerosols' properties in MERRA-2.

<figure><img src="image_29.png" /><figcaption>**Figure 27.** MERRA-2: 2D histograms of $KT_{ref}^h$ (horizontal axis) and $KT_{clear\_merra2}^h$ (vertical axis) for each trimester at Gobabeb. The number of elements in a class increases from yellow to deep blue.</figcaption></figure>

Figure 28 exhibits boxplots of the differences ($KT_{clear\_merra2}^h - KT_{ref}^h$) for all stations for several classes of $\theta_S^{eff\_h}$ for each trimester. Changes from year to year are limited for a given trimester. The bias and the standard deviation of the differences vary with $\theta_S^{eff\_h}$ and the trimester. During DJF, the bias is negative and decreases as $\theta_S^{eff\_h}$ increases. The differences are most often negative, meaning a very frequent underestimation of $KT_{ref}^h$ by $KT_{clear\_merra2}^h$. The pattern is fairly similar during SON, and, to a lesser extent, MAM, though the bias is positive, except at the greatest $\theta_S^{eff\_h}$, and differences are more balanced between negative and positive values. During JJA, the bias is always positive with no influence of $\theta_S^{eff\_h}$; $KT_{clear\_merra2}^h$ very frequently overestimates $KT_{ref}^h$. During each trimester, the standard deviation increases with $\theta_S^{eff\_h}$.
Köppen-Geiger class: all climates (28 stations) / reanalysis: merra2

<figure><img src="image_30.png" /><figcaption>**Figure 28.** MERRA-2: boxplots of the differences ($KT_{clear\_merra2}^h - KT_{ref}^h$) for all stations merged for several classes of $\theta_s^{eff\_h}$ for each trimester. Percentages on the right indicate the fraction of cases in each class.</figcaption></figure>

Boxplots for each climate, each season, and each year for the $\theta_s^{eff\_h}$ classes exhibit complex patterns; no trend can be visually detected, except at the polar stations, as already reported. Whatever the climate, the differences are frequently positive. They sometimes fluctuate with the year or with $\theta_s^{eff\_h}$. Standard deviations tend to increase with $\theta_s^{eff\_h}$. A deeper analysis of the links between differences, climate, and $\theta_s^{eff\_h}$ is beyond our investigative means.

Table 5 summarizes the bias, standard deviation of differences, and root mean square of differences (RMSD) in $E^h$ and $KT^h$ at each station and all stations merged for all trimesters merged. The relative biases in $E^h$ and $KT^h$ are often positive and range from -7% to 4%.
The relative standard deviation in $E^h$ or $KT^h$ has fairly low variations as it is often in the range [3, 4] %.

8. Conclusions

The present work has compared estimates of the hourly means in clear-sky conditions
provided by the ERA5, JRA-3Q, and MERRA-2 reanalyses against coincident ground-based
measurements from 28 BSRN stations covering four years (2015–2018) in an attempt to
identify strengths and potential drawbacks in the modelling of $E_{clear}^h$ and to inform the
agencies in charge of the reanalyses of further developments as uncertainties in $E_{clear}^h$ affect
the estimates of the SIS in all conditions. It is stressed here that the presented work deals
with hourly values $E_{clear\_era5}^h$, $E_{clear\_jra3q}^h$, and $E_{clear\_merra2}^h$, thus complementing the studies
already mentioned, dedicated to the analysis of daily and monthly means of SIS under
clear-sky conditions.

The 28 stations are located worldwide but still cover a limited number of climates and
a limited number of geographical and meteorological conditions. The BSRN measurements
are 1-min measurements and were processed to identify clear-sky conditions. The selected
measurements were then aggregated into hourly clear-sky irradiance $E_{ref}^h$. These processes
decrease the information content of the initial 1-min BSRN measurements to a limited
extent, but they are necessary to match the temporal resolution of the reanalysis.

It was found that ERA5 very frequently underestimates $E_{ref}^h$. The bias per trimester per station ranges from -7% to 5%, relative to the corresponding average in $E_{ref}^h$. The relative biases are most often negative: 70% of them are less than -1% and 52% are less than -2%. The bias varies little from one year to another and from one trimester to another at a given station, but, significantly, from one station to another. As a whole, the differences ($E_{clear\_era5}^h - E_{ref}^h$) tend to be more and more negative as the effective solar zenith angle $\theta_S^{eff\_h}$ increases. The fidelity or consistency to $E_{ref}^h$ is satisfactory as fluctuations in time are well captured by ERA5 with little dispersion: the correlation coefficients are very high and the standard deviations of differences per trimester relative to the corresponding average in $E_{ref}^h$ are fairly constant in time and across stations and are small: they are less than 2% in 60% of the cases, and less than 3% in 90% of the cases. From our results, we conclude that ERA5 can be considered as fairly reliable in time, despite its use of an aerosol climatology, i.e., it offers the same level of uncertainty because the statistical indicators do not fluctuate significantly with time at a given station. However, we observed that this is not true at certain locations, such as Gobabeb, where the use of a climatology leads to notable fluctuations in statistical indicators over time. In contrast, ERA5 is not reliable in space because of the large fluctuations in bias from station to station.

Our results suggest some areas for improving estimates of the hourly means in clear-sky conditions provided by ERA5. Further analyses are necessary to understand why $E_{clear\_era5}^h$ so frequently underestimates $E_{ref}^h$ and the reasons for the spatial variability of the bias. In particular, experiments could be conducted to assess whether the aerosol climatology alone explains the underestimation and its fluctuations in space. Another hypothesis to be tested is whether the underestimation is partly due to a systematic underestimation of the ground albedo.

As for JRA-3Q, the bias per trimester per station ranges from -9% to 6%, relative to
the corresponding average in $E_{ref}^h$. Values are fairly well balanced between negative and
positive values: 28% of the cases are in the range [-3, -1]%, 34% in [-1, 1]%, and 19% in
[1, 3]%. The bias most often exhibits low variations from year to year but large variations
across trimesters and stations are observed. Fluctuations in time are well captured by
JRA-3Q: the correlation coefficients are very high. The standard deviations of differences
per trimester relative to the corresponding average in $E_{ref}^h$ are small: they are less than
2% in 75% of the cases and less than 3% in 96% of the cases. Most often, the standard deviation has low variations across stations and years during a given season. In contrast, changes across seasons are weak in the southern hemisphere and significant in the northern hemisphere. Despite the low values of the standard deviations, these possible changes with the season do not guarantee a good fidelity to $E_{ref}^h$. No significant overall trend with $\theta_s^{eff\_h}$ was observed. From our results, we conclude that JRA-3Q is not reliable in time or space because of the large changes in bias and, sometimes, standard deviation from one station to another and from one trimester to another.

Our results suggest, as for ERA5, that additional analyses should be conducted to understand the origin of spatial and seasonal variations in bias and potentially identify ways to improve the model. One could also assess the impact of using aerosol climatology on the uncertainty of clear-sky estimates, especially to evaluate whether the observed differences in standard deviation and bias between regions and seasons can be explained by the use of a climatology.

The results for MERRA-2 regarding $E^h$ are strongly affected by the use of the mean solar time instead of the true solar time, which leads to a miscalculation of the solar zenith angle $\theta_s$. This miscalculation is propagated to the irradiance at the top of the atmosphere, $E_{0\_merra2}^h$, and also to the clear-sky SIS $E_{clear\_merra2}^h$ and clearness index $KT_{clear\_merra2}^h$ due to the dependence of absorption and scattering by atmospheric constituents on $\theta_s$. Among many examples of this influence on the statistical indicators, the correlation coefficients exhibit a pattern that has seasonal and latitudinal dependency that can affect intraday calculations. At a given station during a given trimester, the bias does not vary much with the year. It varies more across trimesters and stations. Only 43% of the cases have a relative standard deviation less than 3% and 23% of them have values greater than 4%. The standard deviation exhibits a marked seasonal pattern with an overall increase during DJF and SON compared to MAM and JJA. The standard deviation increases with $\theta_s^{eff\_h}$ while the link between the bias and $\theta_s^{eff\_h}$ is more complex and depends on the season. From our results, we conclude that MERRA-2 is not reliable in time or space because of the large changes in statistical indicators from one station to another and from one trimester to another.

Our recommendations for improving estimates of the hourly means in clear-sky conditions provided by MERRA-2 are as follows. Correct the time system. This should improve all statistical indicators. A study similar to this one may then reveal some features whose analysis can lead to further improvements. Attention should be paid to the extreme latitudes as the four polar stations exhibit the most noticeable negative biases compared to the other stations.

This work exposes the strengths and weaknesses of each reanalysis in clear-sky conditions and formulates suggestions for further improvements in the estimates of the hourly means of irradiance $E_{clear}^h$ in such conditions. Currently, none of the three studied reanalyses is satisfactory in cloud-free conditions and all are far from matching the quality of advanced models dedicated to the estimates of $E_{clear}^h$ such as CAMS-McClear and others, as reported in recent papers [19,28].

In addition to the assessment of the hourly means in clear-sky conditions provided by the ERA5, JRA-3Q, and MERRA-2 reanalyses, this paper has presented a new algorithm for detecting cloud-free instants that produces similar results to two other equivalent methods. This new algorithm requires further analysis to rigorously assess its strengths and limitations. This will be addressed in a separate article specifically dedicated to this investigation.
**Supplementary Materials:** The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/atmos16080949/s1, Figure S1: Bias, standard deviation of error, root mean square error, and slope of the three reanalyses for the variable E, analysed by trimester and station. Figure S2: Bias, standard deviation of error, root mean square error, and slope of the three reanalyses for the variable KT, by trimester and station. Figure S3: Scatterplots of reference versus analysis values of E for the three reanalyses, presented separately for individual stations as well as for the merged dataset of all stations. Figure S4: Scatterplots of reference versus analysis values of KT for the three reanalyses, presented separately for individual stations as well as for the merged dataset of all stations. Figure S5: Boxplots of the differences between reanalysis and reference values of E, shown for each individual station and for all stations combined, categorized by several classes of effective solar zenith angle $\theta_s^{\text{eff-h}}$ for each trimester. Figure S6: Boxplots of the differences between reanalysis and reference values of KT, shown for each individual station and for all stations combined, categorized by several classes of effective solar zenith angle $\theta_s^{\text{eff-h}}$ for each trimester.

**Author Contributions:** Conceptualization, methodology, validation, formal analysis, and visualization, Y.-M.S.-D. and L.W.; software, investigation, resources, data curation, supervision, project administration, and funding acquisition, Y.-M.S.-D.; writing—original draft preparation, L.W.; and writing—review and editing, Y.-M.S.-D. and L.W. All authors have read and agreed to the published version of the manuscript.

**Funding:** This research received no external funding.

**Data Availability Statement:** The BSRN measurements used in this study are available either via the BSRN FTP server (https://bsrn.awi.de/data/data-retrieval-via-ftp/, accessed on 11 July 2025) or through our institutional THREDDS data server (https://tds.webservice-energy.org/thredds/catalog/bsrn-stations/catalog.html, accessed on 11 July 2025). The reanalysis datasets are openly accessible: ERA5 was downloaded from the ECMWF Copernicus Climate Data Store (CDS), JRA-3Q from the THREDDS data server operated by UCAR, and MERRA-2 from the NASA Earthdata website.

**Acknowledgments:** The authors thank all operators of the ground stations of the BSRN network for their valuable measurements and the Alfred-Wegener Institute for hosting the BSRN website. We also acknowledge the institutions responsible for producing and maintaining the reanalysis datasets used in this study, namely ECMWF for ERA5, NASA GMAO for MERRA-2, and the Japan Meteorological Agency (JMA) for JRA-3Q.

**Conflicts of Interest:** The authors declare no conflicts of interest.

## References

1. Bojinski, S.; Verstraete, M.; Peterson, T.C.; Richter, C.; Simmons, A.; Zemp, M. The Concept of Essential Climate Variables in Support of Climate Research, Applications, and Policy. *Bull. Am. Meteorol. Soc.* **2014**, *95*, 1431–1443. [CrossRef]

2. Lean, J.; Rind, D. Climate Forcing by Changing Solar Radiation. *J. Clim.* **1998**, *11*, 3069–3094. [CrossRef]

3. Group on Earth Observations. *GEO Task US-09-01a: GEO User Requirements for Air Quality Report*; Group on Earth Observations: Geneva, Switzerland, 2013.

4. Dantas De Paula, M.; Gómez Giménez, M.; Niamir, A.; Thurner, M.; Hickler, T. Combining European Earth Observation Products with Dynamic Global Vegetation Models for Estimating Essential Biodiversity Variables. *Int. J. Digit. Earth* **2020**, *13*, 262–277. [CrossRef]

5. Group on Earth Observations. *The GEOSS Water Strategy: From Observations to Decisions*; Group on Earth Observations: Geneva, Switzerland, 2014.

6. Ranchin, T.; Trolliet, M.; Ménard, L.; Wald, L. Which Variables Are Essential for Renewable Energies? *Int. J. Digit. Earth* **2020**, *13*, 253–261. [CrossRef]

7. Juzeniene, A.; Brekke, P.; Dahlback, A.; Andersson-Engels, S.; Reichrath, J.; Moan, K.; Holick, M.F.; Grant, W.B.; Moan, J. Solar Radiation and Human Health. *Rep. Prog. Phys.* **2011**, *74*, 066701. [CrossRef]

8. Lefevre, M.; Blanc, P.; Espinar, B.; Gschwind, B.; Menard, L.; Ranchin, T.; Wald, L.; Saboret, L.; Thomas, C.; 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. Top. Appl. Earth Obs. Remote Sens.* **2014**, *7*, 1745–1753. [CrossRef]

9. Wald, L. *Fundamentals of Solar Radiation*, 1st ed.; CRC Press: Boca Raton, FL, USA, 2021; ISBN 978-1-003-15545-4.
10. Hersbach, H.; Bell, B.; Berrisford, P.; Hirahara, S.; Horányi, A.; Muñoz-Sabater, J.; Nicolas, J.; Peubey, C.; Radu, R.; Schepers, D.; et al. The ERA5 Global Reanalysis. *Q. J. R. Meteorol. Soc.* **2020**, *146*, 1999–2049. [CrossRef]

11. Kobayashi, S.; Ota, Y.; Harada, Y.; Ebita, A.; Moriya, M.; Onoda, H.; Onogi, K.; Kamahori, H.; Kobayashi, C.; Endo, H.; et al. The JRA-55 Reanalysis: General Specifications and Basic Characteristics. *J. Meteorol. Soc. Jpn.* **2015**, *93*, 5–48. [CrossRef]

12. Kosaka, Y.; Kobayashi, S.; Harada, Y.; Kobayashi, C.; Naoe, H.; Yoshimoto, K.; Harada, M.; Goto, N.; Chiba, J.; Miyaoka, K.; et al. The JRA-3Q Reanalysis. *J. Meteorol. Soc. Jpn.* **2024**, *102*, 49–109. [CrossRef]

13. Gelaro, R.; McCarty, W.; Suárez, M.J.; Todling, R.; Molod, A.; Takacs, L.; Randles, C.A.; Darmenov, A.; Bosilovich, M.G.; Reichle, R.; et al. The Modern-Era Retrospective Analysis for Research and Applications, Version 2 (MERRA-2). *J. Clim.* **2017**, *30*, 5419–5454. [CrossRef]

14. Oumbe, A.; Qu, Z.; Blanc, P.; Lefèvre, M.; Wald, L.; Cros, S. Decoupling the Effects of Clear Atmosphere and Clouds to Simplify Calculations of the Broadband Solar Irradiance at Ground Level. *Geosci. Model Dev.* **2014**, *7*, 1661–1669. [CrossRef]

15. Wang, Z.; Zhang, M.; Wang, L.; Feng, L.; Ma, Y.; Gong, W.; Qin, W. Long-Term Evolution of Clear Sky Surface Solar Radiation and Its Driving Factors over East Asia. *Atmos. Environ.* **2021**, *262*, 118661. [CrossRef]

16. Trolliet, M.; Walawender, J.P.; Bourlès, B.; Boilley, A.; Trentmann, J.; Blanc, P.; Lefèvre, M.; Wald, L. Downwelling Surface Solar Irradiance in the Tropical Atlantic Ocean: A Comparison of Re-Analyses and Satellite-Derived Data Sets to PIRATA Measurements. *Ocean Sci.* **2018**, *14*, 1021–1056. [CrossRef]

17. Lefèvre, M.; Oumbe, A.; Blanc, P.; Espinar, B.; Gschwind, B.; Qu, Z.; Wald, L.; Schroedter-Homscheidt, M.; Hoyer-Klick, C.; Arola, A.; et al. McClear: A New Model Estimating Downwelling Solar Radiation at Ground Level in Clear-Sky Conditions. *Atmos. Meas. Tech.* **2013**, *6*, 2403–2418. [CrossRef]

18. Gschwind, B.; Wald, L.; Blanc, P.; Lefèvre, M.; Schroedter-Homscheidt, M.; Arola, A. Improving the McClear Model Estimating the Downwelling Solar Radiation at Ground Level in Cloud-Free Conditions—McClear-v3. *Meteorol. Z.* **2019**, *28*, 147–163. [CrossRef]

19. Wandji Nyamsi, W.; Saint-Drenan, Y.-M.; Arola, A.; Wald, L. Further Validation of the Estimates of the Downwelling Solar Radiation at Ground Level in Cloud-Free Conditions Provided by the McClear Service: The Case of Sub-Saharan Africa and the Maldives Archipelago. *Atmos. Meas. Tech.* **2023**, *16*, 2001–2036. [CrossRef]

20. Babar, B.; Graversen, R.; Boström, T. Solar Radiation Estimation at High Latitudes: Assessment of the CMSAF Databases, ASR and ERA5. *Sol. Energy* **2019**, *182*, 397–411. [CrossRef]

21. Mabasa, B.; Lysko, M.D.; Moloi, S.J. Validating Hourly Satellite Based and Reanalysis Based Global Horizontal Irradiance Datasets over South Africa. *Geomatics* **2021**, *1*, 429–449. [CrossRef]

22. Sawadogo, W.; Bliefernicht, J.; Fersch, B.; Salack, S.; Guug, S.; Diallo, B.; Ogunjobi, K.O.; Nakoulma, G.; Tanu, M.; Meilinger, S.; et al. Hourly Global Horizontal Irradiance over West Africa: A Case Study of One-Year Satellite- and Reanalysis-Derived Estimates vs. in Situ Measurements. *Renew. Energy* **2023**, *216*, 119066. [CrossRef]

23. Tahir, Z.U.R.; Azhar, M.; Blanc, P.; Asim, M.; Imran, S.; Hayat, N.; Shahid, H.; Ali, H. The Evaluation of Reanalysis and Analysis Products of Solar Radiation for Sindh Province, Pakistan. *Renew. Energy* **2020**, *145*, 347–362. [CrossRef]

24. Tahir, Z.U.R.; Amjad, M.; Asim, M.; Azhar, M.; Farooq, M.; Ali, M.J.; Ahmad, S.U.; Amjad, G.M.; Hussain, A. Improving the Accuracy of Solar Radiation Estimation from Reanalysis Datasets Using Surface Measurements. *Sustain. Energy Technol. Assess.* **2021**, *47*, 101485. [CrossRef]

25. Sianturi, Y.; Marjuki, Sartika, K. Evaluation of ERA5 and MERRA2 Reanalyses to Estimate Solar Irradiance Using Ground Observations over Indonesia Region. In Proceedings of the International Energy Conference Astechnova 2019, Yogyakarta, Indonesia, 30–31 October 2019; p. 020002.

26. Urraca, R.; Huld, T.; Gracia-Amillo, A.; Martinez-de-Pison, F.J.; Kaspar, F.; Sanz-Garcia, A. Evaluation of Global Horizontal Irradiance Estimates from ERA5 and COSMO-REA6 Reanalyses Using Ground and Satellite-Based Data. *Sol. Energy* **2018**, *164*, 339–354. [CrossRef]

27. Reno, M.J.; Hansen, C.W. Identification of Periods of Clear Sky Irradiance in Time Series of GHI Measurements. *Renew. Energy* **2016**, *90*, 520–531. [CrossRef]

28. Bright, J.M.; Sun, X.; Gueymard, C.A.; Acord, B.; Wang, P.; Engerer, N.A. Bright-Sun: A Globally Applicable 1-Min Irradiance Clear-Sky Detection Model. *Renew. Sustain. Energy Rev.* **2020**, *121*, 109706. [CrossRef]

29. Ohmura, A.; Gilgen, H.; Hegner, H.; Müller, G.; Wild, M.; Dutton, E.G.; Forgan, B.; Fröhlich, C.; Philipona, R.; Heimo, A.; et al. Baseline Surface Radiation Network (BSRN/WCRP): New Precision Radiometry for Climate Research. *Bull. Amer. Meteor. Soc.* **1998**, *79*, 2115–2136. [CrossRef]

30. Tegen, I.; Hollrig, P.; Chin, M.; Fung, I.; Jacob, D.; Penner, J. Contribution of Different Aerosol Species to the Global Aerosol Extinction Optical Thickness: Estimates from Model Results. *J. Geophys. Res.* **1997**, *102*, 23895–23915. [CrossRef]

31. Tanré, D.; Geleyn, J.-F.; Slingo, J.-M. First Results of the Introduction of an Advanced Aerosol-Radiation Interaction in the ECMWF Low Resolution Global Model. In *Aerosols and Their Climatic Effects*; A. Deepak Publishing: Hampton, VA, USA, 1984; pp. 133–177.

32. Sato, M.; Hansen, J.E.; McCormick, M.P.; Pollack, J.B. Stratospheric Aerosol Optical Depths, 1850–1990. *J. Geophys. Res.* **1993**, *98*, 22987–22994. [CrossRef]
33. Randles, C.A.; Da Silva, A.M.; Buchard, V.; Colarco, P.R.; Darmenov, A.; Govindaraju, R.; Smirnov, A.; Holben, B.; Ferrare, R.; Hair, J.; et al. The MERRA-2 Aerosol Reanalysis, 1980 Onward. Part I: System Description and Data Assimilation Evaluation. *J. Clim.* **2017**, *30*, 6823–6850. [CrossRef]

34. Blanc, P.; Wald, L. The SG2 Algorithm for a Fast and Accurate Computation of the Position of the Sun for Multi-Decadal Time Period. *Sol. Energy* **2012**, *86*, 3072–3083. [CrossRef]

35. Mamajek, E.E.; Prsa, A.; Torres, G.; Harmanec, P.; Asplund, M.; Bennett, P.D.; Capitaine, N.; Christensen-Dalsgaard, J.; Depagne, E.; Folkner, W.M.; et al. International Astronomical Union General Assembly Resolution B3 on Recommended Nominal Conversion Constants for Selected Solar and Planetary Properties. *arXiv* **2015**, arXiv:1510.07674.

36. Blanc, P.; Wald, L. On the Effective Solar Zenith and Azimuth Angles to Use with Measurements of Hourly Irradiation. *Adv. Sci. Res.* **2016**, *13*, 1–6. [CrossRef]

37. Long, C.N.; Dutton, E.G. BSRN Global Network Recommended QC Tests, V2. 2002. Available online: https://epic.awi.de/id/eprint/30083/1/BSRN_recommended_QC tests_V2.pdf (accessed on 11 July 2025).

38. Mueller, R.; Dagestad, K.F.; Ineichen, P.; Schroedter, M.; Cros, S.; Dumortier, D.; Kuhlemann, R.; Olseth, J.A.; Piernavieja, G.; Reise, C.; et al. Rethinking Satellite-Based Solar Irradiance modellingThe SOLIS Clear-Sky Module. *Remote Sens. Environ.* **2004**, *91*, 160–174. [CrossRef]

39. ASTM G167-15; Test Method for Calibration of a Pyranometer Using a Pyrheliometer. G03 Committee ASTM International: West Conshohocken, PA, USA, 2015. [CrossRef]

40. ISO 9847:2023; Solar Energy—Calibration of Pyranometers by Comparison to a Reference Pyranometer. International Organization for Standardization (ISO): Geneva, Switzerland, 2023. [CrossRef]

41. Hinkelman, L.M. Differences Between Along-Wind and Cross-Wind Solar Irradiance Variability on Small Spatial Scales. *Sol. Energy* **2013**, *88*, 192–203. [CrossRef]

42. Cox, C.J.; Morris, S.M.; Uttal, T.; Burgener, R.; Hall, E.; Kutchenreiter, M.; McComiskey, A.; Long, C.N.; Thomas, B.D.; Wendell, J. The De-Icing Comparison Experiment (D-ICE): A Study of Broadband Radiometric Measurements under Icing Conditions in the Arctic. *Atmos. Meas. Tech.* **2021**, *14*, 1205–1224. [CrossRef]

**Disclaimer/Publisher's Note:** The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.