Improving the Irradiance Data Measured by Silicon-Based Sensors

: Silicon-based sensors are widely used for monitoring solar irradiance, in particular, in the ﬁeld of Photovoltaic (PV) applications. We present a method to correct the global horizontal irradiance measured by silicon-based sensors that reduces the difference to the standard thermopile sensor measurements. A major motivation to use silicon-based sensors for the measurements of irradiance is their lower cost. In addition, their response time is much lower, and their spectral response is much closer to that of the PV systems. The analysis of the differences is based on evaluating four parameters that inﬂuence the sensor measurements, namely the temperature, cosine error, spectral mismatch, and calibration factor. Based on the analysis, a correction model is applied to the silicon sensors measurements. The model separates measurements under a clear sky and cloudy sky by combining the clearness index and the solar zenith angle. By applying the correction model on the measurements of the silicon-based sensor, the differences between sensor readings have been reduced signiﬁcantly. The relative root mean squared difference (rRMSD) between the daily solar irradiation measured by both sensors decreased from 10.6% to 5.4% after applying the correction model, while relative mean absolute difference (rMAD) decreased from 7.4% to 2.5%. The difference in total annual irradiation decreased from 70 KWh/m 2 (6.5%) to 15 kWh/m 2 (1.5%) by the correction. The presented correction method shows promising results for a further improvement in the accuracy of silicon-based sensors. Silicon sensors in this monocrystalline solar cells from Ingenieurbüro Mencke and Tegtmeyer (IMT) [25] (hereafter referred to simply as SiS). The solar cells are in an aluminum casing with a glass cover and a Pt100 temperature sensor to measure the cell temperature. According to manufacturer speciﬁcations, the sensor’s reading should agree within 5% with the thermopile pyranometer’s reading within an ambient temperature between − 20 to 70 ◦ C and normal incident irradiance.


Introduction
Solar irradiance provides the energy that powers the earth's climate and biosphere. This energy is the primary source for many processes on earth, including processes that sustain living systems and circulation of the atmosphere and oceans [1]. This energy is also readily available for providing electricity and heat for industrial and domestic applications [2]. Accurate measurements of solar irradiance are needed for understanding the primary source of energy input to the earth-atmosphere-ocean system and for evaluation of photovoltaic (PV) power profiles [3]. Moreover, solar irradiance is the most important input parameter for PV software to predict the expected solar energy.
The global horizontal irradiance (GHI) is the most commonly measured component of solar radiation [4] and includes both direct beam and diffuse radiation. GHI is of particular interest to photovoltaic installations and is defined as the total solar radiation per unit area that is intercepted by a flat, horizontal surface.
Due to the increasing interest in solar energy and climate changes, solar irradiance measurements are gaining higher importance compared to recent years. According to IPCC [5], a small variation in solar irradiance can produce natural forcing of earth's climate with global and regional-scale responses. Therefore, accurate irradiance measurements are essential for the detection and attribution of climate change [6].
Instruments measuring solar irradiance may be classified as thermal sensors (thermopile pyranometers) and silicon-based (PV) sensors. Radiometers equipped with thermal sensors are widely used to measure broadband solar irradiance due to their nearly constant spectral response over the whole solar spectral range [7]. However, thermal sensors are costly in terms of hardware and calibration.
On the other hand, silicon-based radiation sensors provide the simplest and cheapest alternative. However, these sensors only respond to wavelengths between 300 and 1200 nm, and their spectral response within this interval is not uniform. The response to the red and near-infrared light is noticeably higher than for blue and ultraviolet light. This limited and non-uniform spectral response causes a spectral mismatch of the broadband irradiance measurement [8].
Many studies have highlighted the differences between both sensors and evaluated the uncertainties for solar radiation measurements. Understanding these differences is important because PV system performance analysis often depends on accurate solar irradiance data, and sensors of different technologies may be used. Dunn et al. [9] calculated typical measurement uncertainties for PV sensor and thermopile pyranometer measurements. The calculations were performed for a fixed-tilt system under clear sky conditions. They found the uncertainties in irradiance measurements to be in the order of ±2.4% for PV reference devices and ±5% for thermopile pyranometers. Meydbray et al. [10] considered the spectral effect and calculated the deviation of the two sensors in four different locations in the USA. They found that the daily solar irradiation deviation can amount up to 3% and that this deviation is highly variable over days, weeks, and months. Azouzoute et al. [11] compared the Plane of Array irradiance data, measured at a tilt plan of 32 • by a first-class thermopile sensor and a reference PV cell. They found that the monthly deviation between the solar irradiation measurements from both devices ranges between −4.8% in June to −0.7% in December. Several researchers have developed correction methods that reduce the systematic errors of silicon sensors. The temperature correction is almost similar in all versions, while the methods used to correct the spectral effects vary between the publications. The corrections depend on the sensor temperature, the solar zenith angle (SZA), air mass (AM), diffuse horizontal irradiance (DHI), and global horizontal irradiance (GHI).
In the following, we present different approaches for spectral corrections. Alados-Arboledas et al. [12] used tabular factors for different sky parameters and a functional correction depending on the solar incidence angle. King and Myers [13] proposed functional corrections based on air mass and the angle of incidence derived for global irradiation. Vignola [14] further developed this approach and included diffuse and subsequently direct beam irradiance. Moreover, the German Aerospace Center (DLR) has developed a method using functional corrections, including a particular spectral parameter. The method was developed in 2003 and based on global, diffuse, and direct irradiance [15]. Forstinger et al. [16] suggested a new correction and calibration method based on a physical approach. The method aimed to remove the systematic errors and is based on information of the sensor properties, which includes its directional response and the site's atmospheric conditions. This study aims to examine the differences in the readings of thermopile pyranometer and silicon sensor at a timescale of minutes, days, months, and years in Hannover. The goal is to evaluate the irradiance data from each sensor, identify contributing factors behind the variations between the readings, and create simple empirical correction methods that require a minimum number of inputs to improve the silicon sensor measurements and reduce the differences between both sensors.

Instruments
The measurement system installed on the rooftop of IMUK consists of different meteorological and radiation instruments. The systems have been described shortly in Mubarak et al. [17,18]. We will focus only on the irradiance devices, as their data will be used in this work, namely the thermopile pyranometers and silicon-based sensors.

Thermopile Pyranometer
Pyranometers are broadband instruments that measure global solar irradiance received from the whole skydome (2 Ω solid angle). Besides global irradiance, pyranometers can also measure diffuse irradiance. For this, a small shading disk can be mounted on an automated solar tracker. Alternatively, a shadow ring may prevent the direct component from reaching the sensor the whole day. The International Standard [19] and the World Meteorological Organization (WMO) distinguishes three classes of pyranometers based on performance characteristics and specifications: the best class is called secondary standard, the second-best is called first class, and the third is called second-class. The thermopile sensors used in this study are CMP11 pyranometers from Kipp & Zonen [20] (hereafter referred to simply as Pyr). The sensors were used to measure GHI and DHI.
A typical pyranometer (Figure 1) consists of a black-painted disk (detector) sealed by two glass domes, which protect the sensor from thermal convection and weather threat (e.g., rain and wind). The double glass domes also limit the instrument's spectral sensitivity in the wavelength range between 280 nm and 2800 nm and usually have a bubble for silicon-based. Moreover, the specially designed double glass domes also produce a more accurate angular response in the instrument and reduce thermal losses [21]. barak et al. [17,18]. We will focus only on the irradiance devices, as their data will in this work, namely the thermopile pyranometers and silicon-based sensors.

Thermopile Pyranometer
Pyranometers are broadband instruments that measure global solar irradia ceived from the whole skydome (2Ω solid angle). Besides global irradiance, pyrano can also measure diffuse irradiance. For this, a small shading disk can be mounte automated solar tracker. Alternatively, a shadow ring may prevent the direct com from reaching the sensor the whole day. The International Standard [19] and the Meteorological Organization (WMO) distinguishes three classes of pyranometer on performance characteristics and specifications: the best class is called secondary ard, the second-best is called first class, and the third is called second-class. The t pile sensors used in this study are CMP11 pyranometers from Kipp & Zonen [20] ( ter referred to simply as Pyr). The sensors were used to measure GHI and DHI.
A typical pyranometer (Figure 1) consists of a black-painted disk (detector) se two glass domes, which protect the sensor from thermal convection and weathe (e.g., rain and wind). The double glass domes also limit the instrument's spectral s ity in the wavelength range between 280 nm and 2800 nm and usually have a bub silicon-based. Moreover, the specially designed double glass domes also produce accurate angular response in the instrument and reduce thermal losses [21].

Silicon Sensors
During the last decade, silicon-based sensors have been used as low-cost ra instruments to monitor PV plants [22]. Silicon-based sensors have a similar spec sponse (300-1200 nm) as a solar panel. As a result, silicon sensors provide a more a representation of the energy available for conversion to electricity by a solar pa tenth of the cost of a thermopile pyranometer [23]. In general, there exist two mai of silicon sensors to monitor PV solar systems, specifically silicon-photodiode pyr ters and reference solar cells. A photodiode-based pyranometer essentially moni short circuit current of a solar cell under a diffusing lens. The pyranometer body a fusing lens are designed to minimize deviations from a true angular response [22 ever, the irradiance values indicated by these pyranometers, without correction, m fer from the "true" broadband solar irradiance by over 10% [24].

Silicon Sensors
During the last decade, silicon-based sensors have been used as low-cost radiation instruments to monitor PV plants [22]. Silicon-based sensors have a similar spectral response (300-1200 nm) as a solar panel. As a result, silicon sensors provide a more accurate representation of the energy available for conversion to electricity by a solar panel at a tenth of the cost of a thermopile pyranometer [23]. In general, there exist two main types of silicon sensors to monitor PV solar systems, specifically silicon-photodiode pyranometers and reference solar cells. A photodiode-based pyranometer essentially monitors the short circuit current of a solar cell under a diffusing lens. The pyranometer body and diffusing lens are designed to minimize deviations from a true angular response [22]. However, the irradiance values indicated by these pyranometers, without correction, may differ from the "true" broadband solar irradiance by over 10% [24].
Unlike photodiode-based pyranometers, reference cells do not contain a diffuser; they have glazing, which allows as much of the incident solar radiation as possible to pass inside, such as PV modules (see Figure 2). Reference solar cells are expected to have a similar spectral response as photodiode-based pyranometers since photodiodes and reference cells are both solar cells for which the output is monitored in a short circuit configuration [21]. Generally, reference cells have an internal temperature measurement that can be used to adjust the temperature dependence of the measured irradiance.
Silicon sensors used in this research are monocrystalline solar cells from Ingenieurbüro Mencke and Tegtmeyer (IMT) [25] (hereafter referred to simply as SiS). The solar cells are in an aluminum casing with a glass cover and a Pt100 temperature sensor to measure the cell temperature. According to manufacturer specifications, the sensor's reading should agree within 5% with the thermopile pyranometer's reading within an ambient temperature between −20 to 70 • C and normal incident irradiance.
Unlike photodiode-based pyranometers, reference cells do not contain a they have glazing, which allows as much of the incident solar radiation as possibl inside, such as PV modules (see Figure 2). Reference solar cells are expected t similar spectral response as photodiode-based pyranometers since photodiodes erence cells are both solar cells for which the output is monitored in a short circui uration [21]. Generally, reference cells have an internal temperature measurement be used to adjust the temperature dependence of the measured irradiance.
Silicon sensors used in this research are monocrystalline solar cells fro nieurbüro Mencke and Tegtmeyer (IMT) [25] (hereafter referred to simply as S solar cells are in an aluminum casing with a glass cover and a Pt100 temperatur to measure the cell temperature. According to manufacturer specifications, the reading should agree within 5% with the thermopile pyranometer's reading w ambient temperature between −20 to 70 °C and normal incident irradiance. Measurements of both sensors were carried out from January 2016 to Decem The data were stored in data loggers. The loggers recorded all data, such as irr and temperature, from the sensors every one minute. In the following, the mai ences between sensors will be discussed.

Differences between Sensors
As explained above, both sensors demonstrate differences in the used tech measuring principle, and reactions to incident solar irradiance, as well as makin ences in measured irradiance ( Table 1). The main differences between the sensor further demonstrated and discussed. The influence of temperature on the irradiance sensor signal is lower in the pyranometers [26,27] than in silicon devices [13,28]. The different reaction of se temperature is due to the different measurement principle. In the case of silico Measurements of both sensors were carried out from January 2016 to December 2019. The data were stored in data loggers. The loggers recorded all data, such as irradiation and temperature, from the sensors every one minute. In the following, the main differences between sensors will be discussed.

Differences between Sensors
As explained above, both sensors demonstrate differences in the used technology, measuring principle, and reactions to incident solar irradiance, as well as making differences in measured irradiance ( Table 1). The main differences between the sensors will be further demonstrated and discussed. The influence of temperature on the irradiance sensor signal is lower in thermopile pyranometers [26,27] than in silicon devices [13,28]. The different reaction of sensors to temperature is due to the different measurement principle. In the case of silicon-based sensors, the temperature affects the short circuit current of monocrystalline silicon cells; it records higher values at higher temperatures [29]. The used silicon sensors do not have a temperature compensation that reduces temperature dependence. Accordingly, the temperature change during the measurements should be taken into account and corrected. The available SiS have internal temperature sensors that can be used to adjust the temperature dependence of the measured irradiance.

Spectral Response
The spectral response (SR) describes the sensor's sensitivity to radiation of different wavelengths. As defined, thermopile pyranometers measure solar irradiance within a wide wavelength range (290-2800 nm). Unlike thermopile pyranometers, silicon sensors do not respond to all incident wavelengths equally. Thus, they have a non-flat spectral response ( Figure 3). The SR of silicon describes how well a material can utilize the light of a specific wavelength to generate an electric current. The SR of a silicon-based sensor has an upper limit of approximately 1200 nm, which is mainly determined by the spectral response of the used PV material. The spectral response of the used SiS has been shown in previous work [17]. sensors, the temperature affects the short circuit current of monocrystalline silicon cells; it records higher values at higher temperatures [29]. The used silicon sensors do not have a temperature compensation that reduces temperature dependence. Accordingly, the temperature change during the measurements should be taken into account and corrected. The available SiS have internal temperature sensors that can be used to adjust the temperature dependence of the measured irradiance.

Spectral Response
The spectral response (SR) describes the sensor's sensitivity to radiation of different wavelengths. As defined, thermopile pyranometers measure solar irradiance within a wide wavelength range (290-2800 nm). Unlike thermopile pyranometers, silicon sensors do not respond to all incident wavelengths equally. Thus, they have a non-flat spectral response ( Figure 3). The SR of silicon describes how well a material can utilize the light of a specific wavelength to generate an electric current. The SR of a silicon-based sensor has an upper limit of approximately 1200 nm, which is mainly determined by the spectral response of the used PV material. The spectral response of the used SiS has been shown in previous work [17].

Angular Response
The angular "cosine" response reflects the fact that the response of the irradiance sensor varies with the angle of incidence [30]. Global irradiance is measured with instruments assumed to have a true angular response. However, we know from different studies [31][32][33] that no instrument is perfect in this regard. The deviation of global irradiance measurements from the ideal cosine law is known as cosine error. The angular response of thermopile pyranometers is determined by the glass domes and the spectral and spatial uniformity of the thermopile detector. The used thermopile pyranometer has a good angular response, where the maximum deviation from the ideal angular response is less than 10 W/m 2 (up to an incidence angle of 80°) with respect to 1000 W/m 2 irradiance at normal incidence (0°) [34].
The angular response of a silicon sensor can be described as the reduction in sensor output when solar radiation impinges at angles deviating from the normal to the surface. This deviation affects the calibration of sensors and introduces energy losses in photovoltaic conversion. King et al. [13] attributed this decrease to two sources: The first one is a geometrical factor (cosine law). The second source results from the optical properties of the used sensor, primarily from the reflectance of input optics. The influence of incidence angle on the SiS response is shown in Figure 4. With an increasing incidence angle, the relative deviation rises significantly.

Angular Response
The angular "cosine" response reflects the fact that the response of the irradiance sensor varies with the angle of incidence [30]. Global irradiance is measured with instruments assumed to have a true angular response. However, we know from different studies [31][32][33] that no instrument is perfect in this regard. The deviation of global irradiance measurements from the ideal cosine law is known as cosine error. The angular response of thermopile pyranometers is determined by the glass domes and the spectral and spatial uniformity of the thermopile detector. The used thermopile pyranometer has a good angular response, where the maximum deviation from the ideal angular response is less than 10 W/m 2 (up to an incidence angle of 80 • ) with respect to 1000 W/m 2 irradiance at normal incidence (0 • ) [34].
The angular response of a silicon sensor can be described as the reduction in sensor output when solar radiation impinges at angles deviating from the normal to the surface. This deviation affects the calibration of sensors and introduces energy losses in photovoltaic conversion. King et al. [13] attributed this decrease to two sources: The first one is a geometrical factor (cosine law). The second source results from the optical properties of the used sensor, primarily from the reflectance of input optics. The influence of incidence angle on the SiS response is shown in Figure 4. With an increasing incidence angle, the relative deviation rises significantly.  Figure 5 shows the difference between GHI measured by the Pyr and that measured by the SiS under a clear sky and overcast conditions. Under clear sky conditions, the Pyr measured higher irradiance, and the difference between both sensors increases in the morning and evening hours when the solar zenith angle is high. This may mainly be due to the higher cosine error of the used SiS. Once the sun is totally covered by clouds and the present irradiance is only diffuse, the sensitivity of the SiS increases compared to calibration conditions. The SiS measures a higher solar irradiance than Pyr, which has nearly constant spectral sensitivity over the complete solar spectrum. This result confirms that the spectral distribution of incident solar irradiance significantly impacts the irradiance measured by silicon-based sensors. It is also worth investigating the influence of these differences between the sensors on the monthly and annual irradiation. Figure 6 shows the average total monthly and annual uncorrected global horizontal irradiation (2016-2019) measured by both sensors.  Figure 5 shows the difference between GHI measured by the Pyr and that measured by the SiS under a clear sky and overcast conditions. Under clear sky conditions, the Pyr measured higher irradiance, and the difference between both sensors increases in the morning and evening hours when the solar zenith angle is high. This may mainly be due to the higher cosine error of the used SiS. Figure 5 shows the difference between GHI measured by the Pyr and that measured by the SiS under a clear sky and overcast conditions. Under clear sky conditions, the Pyr measured higher irradiance, and the difference between both sensors increases in the morning and evening hours when the solar zenith angle is high. This may mainly be due to the higher cosine error of the used SiS. Once the sun is totally covered by clouds and the present irradiance is only diffuse, the sensitivity of the SiS increases compared to calibration conditions. The SiS measures a higher solar irradiance than Pyr, which has nearly constant spectral sensitivity over the complete solar spectrum. This result confirms that the spectral distribution of incident solar irradiance significantly impacts the irradiance measured by silicon-based sensors. It is also worth investigating the influence of these differences between the sensors on the monthly and annual irradiation. Figure 6 shows the average total monthly and annual uncorrected global horizontal irradiation (2016-2019) measured by both sensors. Once the sun is totally covered by clouds and the present irradiance is only diffuse, the sensitivity of the SiS increases compared to calibration conditions. The SiS measures a higher solar irradiance than Pyr, which has nearly constant spectral sensitivity over the complete solar spectrum. This result confirms that the spectral distribution of incident solar irradiance significantly impacts the irradiance measured by silicon-based sensors.
It is also worth investigating the influence of these differences between the sensors on the monthly and annual irradiation. Figure 6 shows the average total monthly and annual uncorrected global horizontal irradiation (2016-2019) measured by both sensors. GHI measured by the thermopile pyranometer was overall higher than the SiS irradiance. The monthly differences ranged between 3.7% in July and 17.1% in January. The annual Pyr irradiation was 6.5% higher. GHI measured by the thermopile pyranometer was overall higher than the SiS irradiance. The monthly differences ranged between 3.7% in July and 17.1% in January. The annual Pyr irradiation was 6.5% higher. The results demonstrate that both sensors measure close to each other. The difference between their readings is minimal during summer months (when solar irradiance level around midday is close to STC value). Once the solar elevation angle decreased in winter months and accordingly, the solar irradiance was lower compared to the STC value, the SiS measured lower than Pyr, and the difference between the sensors was higher. This result agrees with the results achieved in [17].

Measurement Site
The measurements, which this work is based on, were performed over the course of four years (January 2016-December 2019) on the roof of the Institute for Meteorology and Climatology (IMuK) of the Leibniz Universität Hannover (Hannover, Germany; 52.23° N, 09.42° E and 50 m above sea level). To better evaluate the results, it is necessary to give an overview of the prevailing solar radiation conditions in the measurement site. According to the German Weather Service (DWD), the average annual global horizontal irradiation in Germany (1981-2010) ranges between 951 kWh/m 2 and 1257 kWh/m 2 [35]. Solar radiation in Hannover as a north German city is close to the minimum value, representing a low radiation climate in this region. Ground measurements for GHI at IMUK in the last ten years (2010-2019) showed that the global irradiation ranged between 989 kWh/m 2 and 1162 kWh/m 2 (Figure 7, upper graphic). The results were within the range of DWD. The results demonstrate that both sensors measure close to each other. The difference between their readings is minimal during summer months (when solar irradiance level around midday is close to STC value). Once the solar elevation angle decreased in winter months and accordingly, the solar irradiance was lower compared to the STC value, the SiS measured lower than Pyr, and the difference between the sensors was higher. This result agrees with the results achieved in [17].

Measurement Site
The measurements, which this work is based on, were performed over the course of four years (January 2016-December 2019) on the roof of the Institute for Meteorology and Climatology (IMuK) of the Leibniz Universität Hannover (Hannover, Germany; 52.23 • N, 09.42 • E and 50 m above sea level). To better evaluate the results, it is necessary to give an overview of the prevailing solar radiation conditions in the measurement site. According to the German Weather Service (DWD), the average annual global horizontal irradiation in Germany (1981-2010) ranges between 951 kWh/m 2 and 1257 kWh/m 2 [35]. Solar radiation in Hannover as a north German city is close to the minimum value, representing a low radiation climate in this region. Ground measurements for GHI at IMUK in the last ten years (2010-2019) showed that the global irradiation ranged between 989 kWh/m 2 and 1162 kWh/m 2 (Figure 7, upper graphic). The results were within the range of DWD. The lower graphic of Figure 7 shows the monthly average of global-and d izontal irradiation in Hannover within the last ten years. GHI had its maxim due to the height of the sun and sunshine duration in this month compare months. It can be seen that the lowest irradiation values were registered in months when the sunshine duration at its minimum. The monthly average in t months was eight to nine times higher than in the winter months. The sun's h the main factor that determines the amount of diffuse irradiance, but the atmos ditions in the measurement location play an essential role. The monthly diffuse makes up about half of the global irradiation during summer months and mor thirds of global irradiation in winter months. These values can give an idea abo radiation climate of Hannover. The big difference between summer and winter and the distribution of diffuse radiation over the year influences the tilt angle tation at which the maximum annual solar radiation is collected.

Correction Model
One of the aims of this work is to improve the SiS solar irradiance measur The lower graphic of Figure 7 shows the monthly average of global-and diffuse horizontal irradiation in Hannover within the last ten years. GHI had its maximum in June due to the height of the sun and sunshine duration in this month compared to other months. It can be seen that the lowest irradiation values were registered in the winter months when the sunshine duration at its minimum. The monthly average in the summer months was eight to nine times higher than in the winter months. The sun's height is not the main factor that determines the amount of diffuse irradiance, but the atmospheric conditions in the measurement location play an essential role. The monthly diffuse irradiation makes up about half of the global irradiation during summer months and more than two-thirds of global irradiation in winter months. These values can give an idea about the solar radiation climate of Hannover. The big difference between summer and winter irradiation and the distribution of diffuse radiation over the year influences the tilt angle and orientation at which the maximum annual solar radiation is collected.

Correction Model
One of the aims of this work is to improve the SiS solar irradiance measurements and reduce the difference to the thermopile sensor. Therefore, we proposed a set of correlations and equations as a correction model to improve the measured SiS-irradiance and to approach the reading of Pyr. The model addressed the main measurement uncertainties that the SiS suffers from in relation to the thermopile pyranometer. The temperature correction was first applied to correct the temperature error caused by changing temperatures during the measurements and the deviation from the STC value. The model started with this correction because the temperature correction depends only on sensor temperature, regardless of sky conditions (clear or cloudy). In the second step, the clear sky data was separated from cloudy-sky conditions using the clearness index (K t ) and SZA. Then cosine, spectral, and calibration corrections were applied according to sky conditions. Figure 8 shows a block diagram of the calculation process of the model.
Energies 2021, 14, x FOR PEER REVIEW was separated from cloudy-sky conditions using the clearness index ( ) and cosine, spectral, and calibration corrections were applied according to sky con ure 8 shows a block diagram of the calculation process of the model.

Temperature Correction ( )
The daytime temperature of SiS as a reference PV cell differs from the am perature since cells are dark in color and accordingly absorb a greater amo energy [36]. During the day, the silicon sensor becomes hotter than the ambie ture by a factor that depends on incident solar irradiance. It is known that the current ( ) increases with increasing cell temperature [36]. The used silicon erate next to a short circuit [25], and accordingly, the sensor's sensitivity incr temperature increase. The influence of temperature on a sensor's response is rected using a temperature coefficient α, which represents the change in w ature changes related to the temperature value under STC. A typical value of ature coefficient for a crystalline silicon (c-Si) PV device is 0.05%/°C [37]. The t

Temperature Correction (C T )
The daytime temperature of SiS as a reference PV cell differs from the ambient temperature since cells are dark in color and accordingly absorb a greater amount of solar energy [36]. During the day, the silicon sensor becomes hotter than the ambient temperature by a factor that depends on incident solar irradiance. It is known that the short circuit current (I sc ) increases with increasing cell temperature [36]. The used silicon sensors operate next to a short circuit [25], and accordingly, the sensor's sensitivity increased as the temperature increase. The influence of temperature on a sensor's response is usually corrected using a temperature coefficient α, which represents the change in I sc with temperature changes related to the temperature value under STC. A typical value of the temperature coefficient for a crystalline silicon (c-Si) PV device is 0.05%/ • C [37]. The temperature coefficient is considered by the manufacturer (IMT), and the correction is done according to Equation (1): where E(T d ,SMM,CosEr,Cf) and E(SMM,CosEr,Cf) are the SiS Irradiance before and after applying the C T , T sen is the sensor temperature measured by the integrated Pt100 temperature sensor that was mounted to the back of the sensor, and α represents the temperature coefficient (0.05%/ • C). The C T depends only on the sensor's temperature, and therefore, the same correction can be applied to the SiS measurements under all weather conditions. In the next step, the model separated the clear sky from cloudy measurements. In this work, a combination of the clearness index (K t ) and the SZA was used to separate the clear sky measurements from those measured under cloudy conditions. Figure 9 shows the change in (K t ) with SZA for clear sky and cloudy global irradiances.
where E 0 is the extraterrestrial irradiance, and E is the measured global irradiance.  = where is the extraterrestrial irradiance, and E is the measured global irradi The extraterrestrial irradiance was calculated as described by Duffie and Beckm [1]. The values of were different in both cases for the same SZA. We used th ences in the clearness index for each SZA range to separate the clear sky irrad those measured under cloudy conditions. Therefore, we defined a critical clear value ( ) for every SZA-range according to Table 2.  The extraterrestrial irradiance was calculated as described by Duffie and Beckman (1980) [1]. The values of K t were different in both cases for the same SZA. We used these differences in the clearness index for each SZA range to separate the clear sky irradiance from those measured under cloudy conditions. Therefore, we defined a critical clearness index value (K tc ) for every SZA-range according to Table 2.

SZA Range
Clearness Index

Cosine Correction (C c )
Cosine error affects mainly the beam irradiance, while diffuse irradiance has less dependence on SZA. Accordingly, the correction model treated the measurement data differently according to sky conditions. The SiS angular response that was measured in the laboratory of manufacture (Figure 3) was used to correct the clear sky global irradiance. In the case of cloudy days, GHI was corrected for the CosEr using a correction factor, which was obtained as a fit function through empirical measurements. The ratio (E SiS /E Pyr ) vs. the SZA is plotted (Figure 10), and a polynomial fitting to the cloudy irradiance tendency was obtained and used to correct the GHI measured in cloudy days according to Equation

Spectral Mismatch Correction ( )
The spectral distribution of the incident solar irradiance varies during the d is also affected by the atmospheric conditions. As a result, it is generally differ

Spectral Mismatch Correction (C s )
The spectral distribution of the incident solar irradiance varies during the day, and it is also affected by the atmospheric conditions. As a result, it is generally different from the STC spectrum (AM1.5). These spectral mismatches lead to a measurement uncertainty sometimes noted as spectral mismatch error that is quantified by the spectral mismatch factor (SMM) [38]. The SMM depends on the semiconductor materials used in PV sensors and the atmospheric conditions in terms of cloud cover and aerosol content [13,39]. The spectral distribution of the beam irradiance differs from that of the diffuse irradiance. Accordingly, spectral mismatch affects both irradiance components differently. The model treated the spectral mismatch of clear sky measurements in a way that differed from cloudy conditions. One way to determine the magnitude of SMM of the sensor under clear sky conditions is to measure the spectral distribution over the day and compare it to the AM1.5 spectral distribution that is used for the calibration. The air mass is calculated using the formula of Kasten and Young [40]. Due to the unavailability of simultaneous spectral measurements, we calculated the spectral irradiance in this work using the UVSPEC model in the Libradtran package [41,42]. As an input parameter, horizontal visibility of 20 km, perceptible water of 15 kg/m 2 , an albedo of 0.02, and a total ozone column of 300 DU were used (Figure 11, left). First, the pyranometer irradiance (E Pyr ) was calculated, and then the SiS irradiance (E SiS ) was determined according to Equations (4) and (5). Then the spectral irradiances at AM1.5 were calculated and used as reference values for both sensors (E Pyr.cal , E SiS.cal ). The spectral mismatch factor was calculated according to Equation (8), 1200 360 E AM1.5 (λ)SR(λ)dλ (8) where, E AM1.5 is the spectral irradiances at AM1.5. The corrected irradiance was then calculated by dividing the uncorrected value by the SMM: where E(SMM,Cf), E(Cf) are the solar irradiance before and after applying the spectral correction. For cloudy days, the SMM was obtained as a fit function through empirical measurements. The ratio (E SiS /E Pyr ) vs. the air mass is plotted in Figure 11 (right), and a polynomial fitting to the cloudy irradiance tendency was obtained and used to correct the GHI measured in cloudy days according to Equation (10). For cloudy days, the SMM was obtained as a fit function through empirical measurements. The ratio ( / ) vs. the air mass is plotted in Figure 11 (right), and a polynomial fitting to the cloudy irradiance tendency was obtained and used to correct the GHI measured in cloudy days according to Equation (10). It is essential to refer here that the SMM may be strongly influenced by the atmospheric conditions and the PV semiconductor materials [39].

Calibration Correction (C cal )
The calibration of silicon sensors was carried out under STC: 1000 W/m 2 of irradiance, AM1.5 spectrum, and 25 • C of cell temperature. In fact, these reference conditions were hardly obtainable in the outdoor measurements as they combined the irradiance of a clear summer day with the sensor temperature of a clear winter day and the spectrum of a clear spring day. Operating the sensors under real atmospheric conditions that deviate from the STC introduces a calibration error. This deviation depends on the real atmospheric conditions that were changing continuously, making the exact evaluation of the error more complicated. We calculated the calibration deviation factor for both clear and cloudy conditions differently in this work. The most influential parameter that affected the calibration factor's magnitude was the solar irradiance prevailing in the field. Calibration factor correction for clear sky conditions (Cf-cl) was calculated by comparing the irradiance measured by both sensors around noontime (about ± 60 min). The calculations were performed after applying the temperature, cosine, and spectral mismatch corrections.
The calibration correction factor (C cal ) was calculated by plotting the ratio (E SiS /E Pyr ) vs. the SZA min (as the irradiance is maximum) and obtaining the correlation equation as a linear fit. The obtained equation was then multiplied by SiS irradiance according to Equation (11): where E(Cf-cl), E corr are the solar irradiance before and after applying the calibration correction for clear sky measurements and SZA min is the minimum daily SZA value.
In the case of cloudy conditions, there was, unfortunately, no reliable cloudy diffuse reference to which the measurements could be compared. Thus, the calibration factor was obtained as a fit function through empirical measurements. The ratio (E SiS /E Pyr ) vs.
the SZA was plotted (after applying the other corrections), and a correlation equation was applied as a correction factor to correct the GHI measured under cloudy conditions Equation (12).
where E(Cf-cd), (E corr ) are the overcast irradiance before and after applying calibration correction.

Results and Discussion
This section presents the results of applying the correction model on SiS measurements. We will start by investigating the effect of individual corrections, namely temperature, cosine, spectral, and calibration corrections, as well as evaluating the magnitude of each correction on the measured irradiance. Finally, we apply all corrections on the SiS measurements and compare them to the Pyr irradiance before and after applying the corrections. Figure 12 describes the influence of individual correction factors on the measured SiS irradiance for two clear sky days. The two case-study days on 29 June (upper graphic) and 16 April (lower graphic) were chosen to investigate the influence of the model on clear sky conditions in different temperatures and solar heights. The value of temperature correction affected the measured irradiance according to sensor temperature with relation to calibration value (green curve). When sensor temperature exceeded the calibration value (25 • C), the SiS measured higher values than it should be, and then the purpose of C T was to reduce the SiS irradiance.  June 2019) before (gray) and after applying each step of the correction model. The corrections were applied to the raw data one after another. First, the temperature correction was applied (green), then cosine correction (gold), then spectral correction (blue), and finally, the calibration correction (brown).
The reduction had its maximum at about noon when the irradiance and thus the temperatures were at maximum. The influence of decreased with increasing SZA, where the temperatures decreased. Once the temperature was below 25 °C, had an increasing influence on SiS irradiance (Figure 10, right).
Due to a large cosine error of the used silicon sensors, the cosine correction played before (gray) and after applying each step of the correction model. The corrections were applied to the raw data one after another. First, the temperature correction was applied (green), then cosine correction (gold), then spectral correction (blue), and finally, the calibration correction (brown).
The reduction had its maximum at about noon when the irradiance and thus the temperatures were at maximum. The influence of C T decreased with increasing SZA, where the temperatures decreased. Once the temperature was below 25 • C, C T had an increasing influence on SiS irradiance (Figure 10, right).
Due to a large cosine error of the used silicon sensors, the cosine correction played the greatest role in correcting the SiS irradiance (golden curve). The losses of SiS irradiance due to cosine error were the highest under clear sky conditions, as shown in Figure 12. C c increased the SiS irradiance depending on the SZA, and its influence was strongest in early mornings and late afternoons.
The blue curve shows the influence of spectral mismatch correction on the SiS irradiance after applying the temperature and cosine corrections. C s increased the SiS irradiance as long as the AM was below the calibration value (AM1.5). Once the air mass exceeded 1.5, the path of solar radiation through the atmosphere was longer, and the incident radiation contained more red light. Accordingly, the response of SiS irradiance increased compared to STC, and then the purpose of C s was to decrease the SiS irradiance. Therefore, the blue curve on both days under study lay above the golden curve around noontime (Am < 1.5) and lay below in the morning and evening hours (Am > 1.5). In general, the spectral mismatch showed seasonal variation depending on the semiconductor material of the sensor [2]. The spectral effect of the SiS (C-si) was smaller, and there was a small seasonal variation [2,3]. Applying the correction model on the SiS measurements under clear sky conditions increased the SiS irradiance (independent of the time of the day) and thus decreased the difference to the Pyr irradiance.
The influence of the correction model on the monthly and annual values is shown in Figure 13. The deviation of the SiS irradiation from the Pyr values decreased significantly after applying the correction model on the SiS measurements. The differences in total annual irradiation decreased from 70 KWh/m 2 to 15 kWh/m 2 with the correction.  after applying the correction model on SiS measurements. The differences between the sensor readings decreased significantly after the correction. The difference in total annual irradiation between the sensor readings decreased from 6.5% to 1.5%.
The ability of the model to decrease the difference between the sensors under cloudy conditions was limited due to the anisotropic distribution of sky radiation and a combination of multiple influencing factors. Figure 14 shows the results of applying the correction model on an overcast day, 16 September 2019. It can be seen that the temperature and cosine corrections had a small increasing influence on the SiS irradiance due to the low sensor temperature (below 25 °C) and the cosine error of the sensor. On the other hand, the spectral correction decreased the irradiance due to the high AM values. However, applying the model to the measurements decreased the SiS irradiance and reduced the difference to the Pyr irradiance. after applying the correction model on SiS measurements. The differences between the sensor readings decreased significantly after the correction. The difference in total annual irradiation between the sensor readings decreased from 6.5% to 1.5%.
The ability of the model to decrease the difference between the sensors under cloudy conditions was limited due to the anisotropic distribution of sky radiation and a combination of multiple influencing factors. Figure 14 shows the results of applying the correction model on an overcast day, 16 September 2019. It can be seen that the temperature and cosine corrections had a small increasing influence on the SiS irradiance due to the low sensor temperature (below 25 • C) and the cosine error of the sensor. On the other hand, the spectral correction decreased the irradiance due to the high AM values. However, applying the model to the measurements decreased the SiS irradiance and reduced the difference to the Pyr irradiance.  Figure 14. The daily curve of the ratio ( / ) for an overcast day (16 September 2019) before (gray) and after applying each step the correction model. First, the temperature correction was applied (green), then cosine (gold) and spectral correction (blue), and finally the calibration correction (brown).
The ability of the correction model to correct the GHI measured by the SiS was also analyzed by means of the relative Root Mean Square Difference (rRMSD) and relative Mean Absolute Difference (rMAD). Table 3 shows the average deviation of the sensors' readings after each correction step. The rRMSD of daily irradiation between SiS and Pyr decreased from 10.6% to 5.4% after applying the correction model, while the rMAD values decreased from 7.4% to 2.5%.

Conclusions
In this work, we focused on the differences between two of the commonly used solar sensors: thermopile pyranometer and silicon-based sensor. These sensors were different in the used technology, measuring principle, and responses to incident solar irradiance. Therefore, the values of irradiance measured with these sensors were different. We aimed to analyze these differences in order to propose a method to improve the SiS measurements by minimizing the difference between SiS and Pyr readings using the minimum number of input parameters. Reducing the errors of silicon-based sensors can open a new avenue for using these sensors as low-cost irradiance sensors. We discussed four parameters that influence the sensor measurements, namely the temperature, cosine error, spectral mismatch, and calibration factor. In order to eliminate the uncertainty caused by these parameters, a correction model based on the analysis was applied to SiS measurements.
Applying the temperature correction on the SiS measurements increased the measured value of irradiance in the winter months when the cell temperature was below the calibration value (25 °C). On the other hand, the correction caused a decrease in the measured value of SiS irradiance on hot days when cell temperatures were above 25 °C.
In comparison to temperature and spectral corrections, the effect of cosine error cor- Figure 14. The daily curve of the ratio (E SiS /E Pyr ) for an overcast day (16 September 2019) before (gray) and after applying each step the correction model. First, the temperature correction was applied (green), then cosine (gold) and spectral correction (blue), and finally the calibration correction (brown).
The ability of the correction model to correct the GHI measured by the SiS was also analyzed by means of the relative Root Mean Square Difference (rRMSD) and relative Mean Absolute Difference (rMAD). Table 3 shows the average deviation of the sensors' readings after each correction step. The rRMSD of daily irradiation between SiS and Pyr decreased from 10.6% to 5.4% after applying the correction model, while the rMAD values decreased from 7.4% to 2.5%.

Conclusions
In this work, we focused on the differences between two of the commonly used solar sensors: thermopile pyranometer and silicon-based sensor. These sensors were different in the used technology, measuring principle, and responses to incident solar irradiance. Therefore, the values of irradiance measured with these sensors were different. We aimed to analyze these differences in order to propose a method to improve the SiS measurements by minimizing the difference between SiS and Pyr readings using the minimum number of input parameters. Reducing the errors of silicon-based sensors can open a new avenue for using these sensors as low-cost irradiance sensors. We discussed four parameters that influence the sensor measurements, namely the temperature, cosine error, spectral mismatch, and calibration factor. In order to eliminate the uncertainty caused by these parameters, a correction model based on the analysis was applied to SiS measurements.
Applying the temperature correction on the SiS measurements increased the measured value of irradiance in the winter months when the cell temperature was below the calibration value (25 • C). On the other hand, the correction caused a decrease in the measured value of SiS irradiance on hot days when cell temperatures were above 25 • C.
In comparison to temperature and spectral corrections, the effect of cosine error correction on the SiS irradiance was greater due to a large cosine error of the silicon sensors under investigation for SZAs larger than 60 • , which is shown in Figure 12. The influence of this correction appeared to be more marked at sunsets and sunrises of clear sky days. Cosine error affects mainly the beam irradiance, while diffuse irradiance had less dependence on SZA.
Spectral mismatch correction increases or decreases the SiS irradiance, depending on the air mass with respect to the AM1.5 spectrum used in STC. For AM values below AM1.5, the correction increased the SiS irradiance. For air mass higher than AM1·5, the spectral distribution of the global solar irradiance changed in favor of the spectral response of SiS. Accordingly, the correction decreased the measured SiS irradiance (Figure 12, blue curve). The spectral mismatch and the spectral correction were strongly influenced by the used PV material and atmospheric conditions.
The model applied as a last correction the calibration factor correction that stems from calibration of SiS under STC, which deviated from the real atmospheric conditions. The magnitude of this correction in clear sky days depends on the solar irradiance; its value is greater in winter months than in summer months. This correction increased the SiS irradiance in clear sky conditions by 2-5% and decreased by up to 3% for overcast conditions.
As a result of applying the correction model on the GHI of SiS, the differences between the values measured with both sensors were reduced significantly. The daily rRMSD between SiS and Pyr irradiances decreased from 10.6% to 5.4% after applying the correction model, while the rMAD decreased from 7.4% to 2.5%.
The model was used to correct the GHI measurements. The accuracy of the model for correcting tilted irradiance was not tested due to the unavailability of thermopile measurements on inclined surfaces. However, we assumed that this model can be used to correct the SiS irradiance measured on tilted surfaces. The assumption is based on the fact that the temperature factor depends only on the sensor temperature and the cosine error depends on the angle of incident (AOI). In addition, the spectral mismatch factor has a negligible dependence on the tilt angle of the sensor [2], and the calibration factor depends on the level of incident irradiance.
It can be concluded that applying the correction model on the SiS irradiance significantly reduces the differences between the readings of the sensors and makes the SiS more useful as a low-cost instrument to measure GHI.
Author Contributions: R.M. conceived and designed the study and wrote the draft paper; H.S. contributed to the discussion, collected the data, and supervised the measurements; G.S. contributed to the conception and discussion of the data. All authors significantly contributed to the final version of the manuscript. All authors have read and agreed to the published version of the manuscript.