Updated GOES-13 Heliosat-2 Method for Global Horizontal Irradiation in the Americas

Increasing the utilization of renewable energy is at the center of most sustainability policies. Solar energy is the most abundant resource of this type on Earth, and optimizing its use requires the optimal estimation of surface solar irradiation. Heliosat-2 is one of the most popular methods of global horizontal irradiation (GHI) estimation. Originally developed for the Meteosat satellite, Heliosat-2 has been modified in previous work to deal with GOES-13 data and named here GOES_H2. This model has been validated through the computation of indicators and irradiation maps for the Guiana Shield. This article proposes an improved version of GOES_H2, which has been combined with a radiative transfer parameterization (RTP) and the McClear clear-sky model (MC). This new version, hereafter designated RTP_MC_GOES_H2, was tested on eight stations from the Baseline Surface Radiation Network, located in North and South America, and covered by GOES-13. RTP_MC_GOES_H2 improves the hourly GHI estimates independently of the type of sky. This improvement is independent of the climate, no matter the station, the RTP_MC_GOES_H2 gives better results of MBE and RMSE than the original GOES_H2 method. Indeed, the MBE and RMSE values, respectively, change from −11.93% to −2.42% and 23.24% to 18.24% for North America and from −4.35% to 1.79% and 19.97% to 17.37 for South America. Moreover, the flexibility of the method may allow to improve results in the presence of snow cover and rainy/variable weather. Furthermore, RTP_MC_GOES_H2 results outperform or equalize those of other operational models.


Introduction
In light of the energy transition challenge, electricity production using renewable energy has received growing interest [1,2]. For instance, the Energy Sector Management Assistance Program shows that the annual photovoltaic power potential is quite high for most countries [3]. However, maximizing solar radiation exploitation requires managing surface irradiation, namely, global horizontal irradiation (GHI) [4,5]. These data can be obtained in two main ways: in situ measurements [6] and estimation using satellite images [7]. Measurements are captured by instruments that require specific installation and continuous calibration and maintenance [8]. As a consequence, satellite-based estimation is receiving increasing interest [9].
Surface irradiation estimation via satellite falls into two main categories: direct and indirect models. Direct approaches use a radiative transfer model to evaluate the radiative loss between the top-of-atmosphere and the surface [10,11]. Indirect hybrid approaches combine physical models and empirical parameters to facilitate the estimation procedure [12,13]. Hence, while approaches in the first group are complex and time-consuming, parameter was introduced to deal with cloud cover in tropical climates, particularly in French Guiana. The introduced radiative transfer parameterization has shown promising results under cloudy skies [28].
The first objective of this paper is to propose an updated version of the GOES_H2 that would be efficient under all types of skies: cloudy and clear. This would be done by replacing the ESRA model by the McClear model in order to benefit from fully physical approach of the clear-sky modeling. In addition, by introducing a radiative transfer parameterization so that cloudy-sky irradiation would be efficiently taken into account and improve GHI estimates.
GHI estimation methods on the American continent are, either regional methods (BRASIL-SR [13], Heliosat family [7,20,22], PSM [15]), either methods with low spatial resolution (POWER [17]), which limits the availability of GHI data in this area. Hence, the second objective of this article is to propose the updated method named hereafter RTP_MC_GOES_H2 for providing accurate GHI estimates at fine resolution (1 × 1 km) and for GOES-13 spatial coverage (81.32°N to 81.32°S and 6.29°E to 156.29°W). Therefore, RTP_MC_GOES_H2 was validated on different areas of the American continent . GHI estimation was here performed on stations of the Baseline Surface Radiation Network (BSRN). The BSRN [29] is a global network of radiation sensors with high accuracy. Eight stations were chosen in North (four) and South (four) America to match the geographical coverage of GOES-13. Moreover, to assess the performance of the method, RTP_MC_GOES_H2 estimates were compared to other operational models (Heliosat-4, POWER and PSM) estimates.
Herein, the method RTP_MC_GOES_H2 and the data are first presented. Then, the results obtained with RTP_MC_GOES_H2 method are compared with those obtained with GOES_H2 and to results from free access web services generated by other methods. These results are then discussed before ending with conclusions.

GOES-13 Data
The geostationary operational environment satellite of the second generation, GOES-13, also known as GOES-N [30], was located at 75°W above the Equator from April 2010 to December 2017. The satellite had one band in the visible range centered on 0.63 µm for a spatial resolution of 1 × 1 km.
Moreover, GOES-13 covered the considered regions twice per hour at approximately minutes 15 and 45.
GOES-13 images from 1 January to 31 December 2015, were downloaded. The National Oceanic and Atmospheric Administration provided GOES-13 images via the Comprehensive Large-Array Stewardship System. The downloaded data were encoded in 16 bits.

GOES-13 Image Requests per Region
The Comprehensive Large-Array Stewardship System allows GOES-13 data to be downloaded for specific regions. Thus, instead of downloading full disk images, requests were made for the regions of interest.
In particular, the BSRN stations located in South America cover a large area. Hence, for computational efficiency, three areas were defined for downloading GOES-13 images: H2 is a hybrid model in which the key step is the computation of a cloud index, n, where n is a factor that quantifies the cloud cover above the specific area of interest. The operating principle is that cloudiness affects the solar irradiation that reaches the surface. Hence, the surface irradiation is effectively estimated by assessing the clearness of the sky. Thus, the model describes a method that, starting from digital counts, allows the estimation of n and then the global horizontal irradiation. The principal steps are described in the following.
GOES images were downloaded and converted into digital counts in 10 bits, followed by calibration for radiance. The best method is provided by the National Aeronautics and Space Administration (NASA) [31]: where L is radiance in W/(m 2 .sr) with sr standing for the steradian unit; d the number of days since launch; g 0 is the initial gain; g 1 , g 2 and g 3 are, respectively, the first-, second-, and third-order polynomial terms of the gain trend that account for the sensor degradation since the launch; CN 10 and CN 0 , are, respectively, the digital counts in 10 bits and an offset (https://www-pm.larc.nasa.gov/cgi-bin/site/showdoc?mnemonic=SAT_CALIB_ USER, 6 December 2021). For further details about the calibration, readers are redirected to [24]. The next step is to convert radiance into reflectance values through: where ρ sat is the satellite-derived reflectance, I 0 sat is total irradiance in the visible channel of the GOES-13 satellite (W/m 2 ), ε is the Earth-Sun distance correction factor and θ s is the solar zenith angle. Apparent albedo, ρ, is deduced from reflectance values using: where ρ atm is the apparent albedo of the atmosphere, θ v and ψ are the satellite viewing angle and the difference between the sun and satellite azimuth angles, respectively, and T(θ s ) and T(θ v ) are the downward and upward transmittances, respectively.
A key step in the H2 method is the computation of a cloud index, n, which allows the quantification of the cloudiness of the sky, the formula for which is: where ρ c and ρ g are the cloud and ground albedo, respectively. The ground albedo is chosen as the second minimum of the albedo in a time range, which, in this study, is a month. The second minimum is selected instead of the absolute minimum to avoid, among other issues, values generated by undetected defects in satellite images. Different methods have been proposed for the cloud albedo computation [24]; however, in this study, the original method introduced in [7] was used. The cloud albedo has been defined as the typical value for the brightest clouds. This value is quite difficult to estimate. The authors of [32] proposed an approach based on the computation of an effective cloud albedo (ρ eff ) using the sun zenith angle.
with ρ eff = 0.85 − 0.13 · 1 − exp −4 cos(θ s ) 5 (6) Empirical studies allow the computation of a clear-sky index, Kc, based on the cloud index [7]: Finally, an estimate of the global horizontal irradiation is given by: where Gc is the clear-sky irradiation-the irradiation that would reach the surface if the atmosphere were cloud-free-whose computation is discussed in the following part.

The MC_GOES_H2 Method
The European Solar Radiation Atlas (ESRA) [23] is the clear-sky model used in the original H2. It is a hybrid approach, combining both physical equations and empirical considerations, based on the Linke turbidity factor.
The Linke turbidity factor aims to quantify the attenuation of the radiation as it passes through the atmosphere [33]. An important problem with this factor is that it is considered fixed for each month. Hence, it does not consider changes in the atmosphere that occur in the same month based on geographical position, climate, etc.
Therefore, a new clear-sky model, namely, the McClear model, was proposed in a previous study [19] using recent results obtained through the Monitoring Atmosphere Composition and Climate (MACC) projects. The model is purported to more rapidly reproduce irradiation estimations obtained by a radiative transfer model by using lookup tables.
The McClear model is entirely physical and meteorological input data are continuously updated every three hours. Indeed, its inputs include atmospheric composition information provided by MACC models that use satellites and in situ data.
Clear-sky irradiation data are available at different time steps through the Copernicus Atmosphere Monitoring Service (CAMS) via a dedicated web interface [34]. In this study, the chosen time step is 1 min at download. Hence, to obtain the clear-sky irradiation every 30 min in the given interval, the integral (sum) of data was calculated on an interval of 30 min centered on the time of satellite acquisition at the corresponding point.
Let us name the MC_GOES_H2 the GOES_H2 method with the use of the McClear clear-sky model.

The RTP_MC_GOES_H2 Method
The H2 method is known to be less accurate for solar surface irradiation estimation under cloud cover. Furthermore, the thicker the clouds, the more difficult it is to estimate the irradiation. Indeed, in [7], the authors show that H2 gives less accurate estimations under cloudy skies in Europe and in [24], the authors have the same conclusions with a study in French Guiana using GOES_H2.
To overcome this issue, reference [28] proposed a modification, for tropical climate, of the equation of global horizontal irradiation (8) by combining clear-and cloudy-sky irradiation terms. Then, the equation for global horizontal irradiation becomes: The main objective was to improve GHI estimations in French Guiana, a region subject to significant cloud coverage, without negatively affecting clear-sky estimations or being time-consuming. Our work aims to generalize this approach to all types of climates and show that the GHI estimation is improved under all skies. Based on the work by [35], the cloudy-sky radiative transfer parameterization (RTP) term, namely, G b , is written as: where α c is cloud absorption (unitless). Radiative transfer models are the best models to estimate this term; however, to allow near-real-time estimation, a linear approximation function of the effective albedo (6) is proposed: The slope a is the one to be tuned. Authors in [36] propose to choose the value of coefficient a (11) between 0.04 and 0.17. For our study, different values were tested, and the best results, presented hereafter, were obtained with a = 0.085.
RTP_MC_GOES_H2 is an updated method that combines the H2 algorithm optimized for GOES-13 data [24] with the McClear clear-sky model [19] and a RTP term.
A flowchart of the methods is presented in Figure 1.

Choice of Ground Stations
The baseline surface radiation network is a global network of radiation measurement instruments [29]. The network allows the acquisition of surface solar irradiation measurements for stations that belong to it. However, maintaining the instruments is not one of its activities. This network is the best database of in situ data: the stations are chosen after rigorous screening.
Eight stations in the BSRN network were considered for this work. First, only stations located in America are relevant for our work since the objective is to validate the RTP_MC_GOES_H2 method in the geographical areas covered by the satellite. Second, we decided to analyze an entire year, so only stations that were active throughout this period were considered. Thus, eight stations were chosen: four in North America and four in South America for the year 2015. That year was chosen because, for all stations but one, measured data were available for all months. In North America, the stations are located in Oklahoma (Billings and Southern Great Plains), Illinois (Bondville), and Mississippi (Goodwin Creek). In South America, only four stations were available for the study period, all of which are in Brazil and located in Brasilia City, South Atlantic Ocean (Florianopolis), Petrolina, and São Martinho da Serra. The stations' information is summarized in Table 1.
The selected stations in North America are in UTC-6, and those in South America are in UTC-3.
Only data for which the cosine of the sun zenith angle is greater than zero are selected to discard the ones for which the sun is completely low in the sky. The study is pursued in all skies condition.
A quality control on extreme values was made for the extreme values according to [37] and the solar radiation Data (SoDa) website (solar radiation data. http://soda-pro.com/, 6 December 2021). Measured hourly surface irradiation was considered valid when: where G TOA is the top-of-atmosphere hourly irradiation, θ s is the sun zenith angle and I 0 the solar constant (1367 W m −2 ). Figure 1. Flowchart of the methods. GOES_H2 method is in the red block, in the green block is the MC_GOES_H2 method and blue block illustrates the radiative transfer parameterization. Here, we used the Köppen Geiger climate classification [38]. This system classifies climates worldwide by using names containing two or three letters. The three letters refer to the main climate, precipitation information, and temperature, respectively.
The stations are grouped into three types of climates: Aw, BSh, and Cfa. Aw corresponds to an equatorial Savanna climate with dry winters and applies to the Brasilia station. BSh is a steppe climate with hot and arid temperatures (that is, the annual temperature is greater than 18°C) and applies to the Petrolina station. Cfa is a warm, very humid climate with hot summers (the maximum temperature is greater than 22°C); the Billings, Bondville, Florianopolis, Goodwin Creek, São Martinho da Serra and Southern Great Plains stations are situated in this type of climate.

Configurations
Data are considered valid for the study if both BSRN measurements and estimations are available. That is, data are treated as missing if either the in situ measurement or estimation is unavailable. A missing estimation is encountered when no satellite image is available or the image is corrupted at the station position.
A general function is the computation of the hourly estimate based on the estimates available per hour. As two acquisitions are available per hour, two surface irradiation estimations should be available for each hour. Consequently, the hourly estimate is the sum of the two 30-min estimates for the previous hour; that is, where GHI h−1 15 is the mean irradiation computed for a 30-min time range, centered on the acquisition at minute 15 of h − 1; the same applies to GHI h−1 45 but at minute 45. The meteorological data (rainfall, snowfall, and temperature) [39] used in the below data analysis were provided by the Modern-Era Retrospective analysis for Research and Applications, Version 2 [40], and are available from the SoDa website (http://www.sodapro.com/web-services/meteo-data/merra, 6 December 2021).

Metrics
Three metrics were computed to validate the estimations compared to the in situ measurements: the mean bias error (MBE), the root-mean-square error (RMSE), and r-squared (R 2 ).

Mean Bias Error
The MBE (14), as its name implies, is the average difference between the estimation and the true value; it is equivalent to the systematic error of the model.
where N is the number of samples, y i is the hourly estimated GHI, and x i is the hourly measured GHI.

Root-Mean-Square Error
The RMSE (15) is a metric that penalizes large differences between estimates and measurements.

R-Squared
The coefficient of determination (16) computes how well the estimations fit the regression model generated with the ground data. The factor takes a value between 0 and 1, with 1 corresponding to an exact fit.
where x is the mean of the measured values. Table 2 shows a comparison of the metrics obtained with GOES_H2, MC_GOES_H2 and RTP_MC_GOES_H2 methods. The MBE of the overall stations is best with RTP_MC_GOES_H2 (−0.74%) than the one obtained with GOES_H2 (−8.91%).

Comparison between GOES_H2, MC_GOES_H2 and RTP_MC_GOES_H2 Estimates
Moreover, the RTP_MC_GOES_H2 outperforms the GOES_H2 for the overall stations regarding the RMSE (17.92% versus 21.94%). Figure 2 represents the scatter plots of the estimates versus the measured data for the overall stations as well as the mesh grids of the MBE depending on the sun zenith angle cosine and the clear-sky index K c for GOES_H2 and RTP_MC_GOES_H2. The scatter plots show that with the RTP_MC_GOES_H2 method, the line of regression is closer to the ideal case (estimates match measures) than with GOES_H2. The mesh grids show that with GOES_H2, the MBE are very low and the RMSE very high when K c is low and the cosine of the sun zenith angle is great. With, the RTP_MC_GOES_H2, these values are improved.
The MBE of the overall stations with MC_GOES_H2 (−11.93%) is slightly better than the one obtained with GOES_H2 (−10.07%) for stations located in North America. However, in South America, the MBE values obtained with MC_GOES_H2 are degraded versus those obtained with GOES_H2. Hence, the McClear clear-sky model has a low impact on the MBE, approximately 1% for North and South America.
RTP_MC_GOES_H2 significantly improves the estimates of GOES_H2 compared to MC_GOES_H2. Indeed, for the MBE, the RTP_MC_GOES_H2 provides an improvement up to ten points to the MBE while the MC_GOES_H2 improvement does not exceed three points. Regarding the RMSE, with the RTP_MC_GOES_H2, a gain up to 6 points is made while MC_GOES_H2 this value drops to 1.5.    On Figure 3, the slope of the line of regression is closer to one with RTP_MC_GOES_H2 than GOES_H2, indeed, the scatter is more compact. In addition, the color map shows that there is greater concentration of points around the ideal line.

Comparison between GOES_H2 and RTP_MC_GOES_H2 for North America
On the mesh grids on Figure 4, it can be seen that with GOES_H2 model, high (absolute) values of MBE are observed for K c < 0.6 regardless of the cosine of the sun zenith angle. Conversely, the RTP_MC_GOES_H2 gives better MBE values for all cases and even more when the cosine of the sun zenith angle is greater than 0.4.
The RMSE mesh grids confirm that the RTP_MC_GOES_H2 gives more accurate estimates than GOES_H2. Indeed, a decrease is visible for all cases.
Next, we focus on the Bondville station, as it produced the least accurate results. Figure 5 shows the monthly RMSE against the snowfall and temperature for this station. It can be seen that the RMSE follows the snowfall evolution: the greater the snowfall, the higher the RMSE. Moreover, for January and February, the average temperature is negative; that is, snow is not only present but also remains on the ground.
(a) (b) Figure 5. Representation of (a) the monthly RMSE and snowfall and (b) average temperature for the Bondville station.
In Figure 3, the scatter plot for the Bondville station contains several points at which in situ GHI values are greater than 300 and the satellite-based estimates are below 200; that is, the model highly underestimates the surface solar irradiation. These points form a sort of spike in the scatter plot. Identifying the data in this spike can be informative. To ensure that only outliers are selected, an ellipse is defined to exclude points that may just correspond to dispersion around the regression line, as shown in Figure 6. An analysis of this spike is summarized in Table 3. There are more than 100 data points involved, and more than half of them are observed during February, which is the month with the greatest snowfall. Almost 30 percent of the underestimated values are in March, followed by November. In fact, the months that have the greatest deviations correspond to those with a relatively large amount of snowfall. In November, only two days are of concern, and for these, the satellite images show persistent snow cover.  In Figure 7, the color map shows that there is a higher concentration of the points around the ideal regression line with the RTP_MC_GOES_H2 than with GOES_H2. Figure 8 shows an improvement for the MBE and the RMSE. Indeed, the metrics have higher values-in absolute values for the MBE-with GOES_H2 than with RTP_MC_GOES_H2. Even if the improvement is more noticeable when the sun is high in the sky, it can be seen that the add of the RTP term improves the estimation for all skies.
The performance is poorer for the Florianopolis station. However, the coefficient of determination is very good. Figure 9 shows that the monthly RMSE increases with rainfall, which can be quite high.
In addition, the station is on an island and in front of the ocean, where sunglint effects can occur. Sunglint is the reflection of the sun on the water surface at sunrise and sunset, which can distort the data acquired by the satellite. Table 4 shows the results obtained for different periods. On the one hand, better results are obtained when omitting the 6 p.m. hour, corresponding to sunset, and the 8 a.m. hour, during which sunrise occurs. On the other hand, the minimum relative RMSE is obtained for the period from 9 a.m. to 5 p.m., but with much less data.

Comparison to Other Estimations
In this part, GOES_H2 and RTP_MC_GOES_H2 metric values are compared to those obtained with GHI estimates obtained with other methods. All data for the latter are free and provided via web services.
The following databases were used: • The Copernicus Atmosphere Monitoring Service (CAMS) radiation service. The spatial coverage is the same as that of Meteosat. Thus, GHI estimates were only compared to BSRN stations in South America. Furthermore, Meteosat has a resolution of 3 km at nadir. Moreover, Heliosat-4 is the algorithm used for the estimation [18]. • The National Solar Radiation Database (NSRDB). NSRDB [41] provides GHI estimation via a physical solar model that uses information from multiple satellites. The resolution of the satellite data is about 4 km.
• The NASA Prediction Of Worldwide Energy Resources (POWER). POWER [17] uses a specific NASA algorithm to estimate GHI. The method combines many types of climatic and meteorological information from different satellites. The data resolution is approximately 50 km. Figure 9. Representation of the monthly rainfall and the monthly RMSE for the Florianopolis station. Figure 10 shows a comparison of the metrics obtained using the five datasets when considering all available data for each.
Regarding the MBE, RTP_MC_GOES_H2 provides more accurate results for all stations except Bondville and Florianopolis, while GOES_H2 estimates are the less accurate for five stations over height.
For the RMSE, the best results are obtained with RTP_MC_GOES_H2 estimates and GOES_H2 is the best second method in South America.
Following the R 2 metric, RTP_MC_GOES_H2 and CAMS give more accurate results followed closely by GOES_H2, except for Bondville.
Hence, RTP_MC_GOES_H2's performance is better of equal to those of CAMS, NSRDB, and POWER models for all the stations but Bondville.

Discussion
The GOES_H2 method was updated and coupled with the McClear clear-sky model and a radiative transfer parameterization. The method was validated on eight BSRN stations located in North and South America.
RTP_MC_GOES_H2 method outperforms the original GOES_H2 method. The introduction of the radiative transfer parameterization and the McClear model shows real improvement comparing to the original method.
RTP_MC_GOES_H2 gives better results than GOES_H2 under all skies, notably when the sky is clear and the sun is high in the sky. The method is also effective, regardless of the climate type, at least for the three types encountered in this study; that is, Aw, BSh, and Cfa in the Köppen Geiger classification.
However, less accurate results were obtained for the Bondville and Florianopolis stations. At Bondville, the RMSE results are highly correlated with the snowfall. Indeed, it is well known that snow cover negatively affects the estimation of irradiation. Snow cover reflects a high percentage of light in the visible range of the spectrum and can be considered as cloud cover in the estimation process. This situation negatively affects the estimation of the cloud index (4) and thus the surface irradiation estimation.
The discrimination between snow and cloud is generally based on the observation that, although snow and cloud are difficult to distinguish in the visible range of the spectrum, they behave differently in other ranges [42,43]. The approach would then be to combine the information of different bands of the satellite to isolate the hours where there is snow and exclude them from the computation of the cloud albedo ρ c . As theH2 method is only designed for the visible range [23], this process would be implemented in a third-party algorithm.
Two main aspects of the Florianopolis station are highlighted: the region is rainy, and it is an island. Thus, Florianopolis, as in all regions on the seafront, is affected by sunglint [44]. At sunrise and sunset, the reflection of the sun at the sea surface artificially induces a high clear-sky index value and, thus, the surface solar irradiation is overestimated.
Different methods are proposed in the literature to address sunglint effects [44]. As the best approach is to exclude data containing sunglint, many methods attempt to estimate whether sunglint will appear in a given hour. However, these estimations are frequently based on wind speed and, thus, it is necessary to be able to retrieve that information.
The metric values of the RTP_MC_GOES_H2 method estimates are better than those of CAMS, NSRDB, and POWER estimations for almost all stations: NSRDB and POWER data produce better results for the Bondville station. The meteorological and climatic data used as inputs in these models make it easier to deal with snow cover.
The RTP approach was introduced for tropical climates, but its effectiveness was demonstrated under all types of climates (encountered in the study). However, improvements can be made by first selecting the best proportionality parameter a from the radiative transfer parameter for each station individually. As the best parameter depends on the types of clouds that may be encountered in the region [28], the best a value should not be the same for all regions.
GOES-13 moved into on-orbit storage in 2018; a new satellite, GOES-16 [45], was launched in 2016 and has been in operation since 2018. Hence, to produce data estimates after 2018, it is necessary to adapt RTP_MC_GOES_H2 to GOES-16 data.

Conclusions
This article proposes an improvement to the GOES_H2 by updating the clear-sky model to the McClear one and introducing a radiative transfer parameterization.
The RTP_MC_GOES_H2 performs better than the GOES_H2 for surface solar irradiation estimation in North and South America.
Adding the RTP parameter improves the GHI estimates regardless of the clear-sky index value. In addition, RTP_MC_GOES_H2 gives lesser values of MBE and RMSE than the original GOES_H2 method no matter the station and the related climate (for the overall stations: −0.74% MBE and 17.92% for RTP_MC_GOES_H2 versus −8.91% MBE and 21.94% RMSE for GOES_H2). Hence, the method is effective regardless of the climate type, at least for the three types encountered in this study; that is, Aw, BSh, and Cfa in the Köppen Geiger classification. Moreover, as expected, the best estimations were obtained when the sky was clear and the sun was high in the sky.
Moreover, RTP_MC_GOES_H2 method outperforms operational CAMS, NSRDB, and POWER models for almost all stations.
Further work will focus on optimizing a method of estimating the proportionality coefficient of the radiative transfer parameter. Thus, the parameter value will adapt to different types of clouds encountered in different climates. Moreover, to enable the application of the method to more locations, establishing a means of discriminating snow from clouds is essential. In addition, addressing sunglint effects will increase the accuracy of estimations on islands, for example. These issues can be addressed by considering other bands of acquisition for the satellite or information from other satellites.
In future work, we also plan to optimize RTP_MC_GOES_H2 for the current operational satellite-GOES-16. Data Availability Statement: The satellite images are available at the Comprehensive Large-Array Stewardship System at https://www.avl.class.noaa.gov/index.html. CAMS data are available at http://www.soda-pro.com/web-services/radiation/cams-radiation-service. NSRDB data are available at https:maps.nrel.gov/nsrdb-viewer/. POWER data are available at https://power.larc. nasa.gov/data-access-viewer/, 6 December 2021.