Next Article in Journal
Prospects and Obstacles Associated with Community Solar and Wind Farms in Jordan’s Suburban Areas
Previous Article in Journal
A Multi-Stage Approach to Assessing the Echo-Tech Feasibility of a Hybrid SAM-CREST Model for Solar PV Power Plants in Maryland, USA
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

New Decomposition Models for Hourly Direct Normal Irradiance Estimations for Southern Africa

by
Francisca Muriel Daniel-Durandt
and
Arnold Johan Rix
*,†
Department of E&E Engineering, Stellenbosch University, Stellenbosch 7602, South Africa
*
Author to whom correspondence should be addressed.
These authors contributed equally to this work.
Solar 2024, 4(2), 269-306; https://doi.org/10.3390/solar4020013
Submission received: 4 March 2024 / Revised: 29 April 2024 / Accepted: 9 May 2024 / Published: 14 May 2024

Abstract

:
This research develops and validates new decomposition models for hourly direct Normal Irradiance (DNI) estimations for Southern African data. Localised models were developed using data collected from the Southern African Universities Radiometric Network (SAURAN). Clustered areas within Southern Africa were identified, and the developed cluster decomposition models highlighted the potential advantages of grouping data based on shared geographical and climatic attributes. This clustering approach could enhance decomposition model performance, particularly when local data are limited or when data are available from multiple nearby stations. Further, a regional Southern African decomposition model, which encompasses a wide spectrum of climatic regions and geographic locations, exhibited notable improvements over the baseline models despite occasional overestimation or underestimation. The results demonstrated improved DNI estimation accuracy compared to the baseline models across all testing and validation datasets. These outcomes suggest that utilising a localised model can significantly enhance DNI estimations for Southern Africa and potentially for developing similar models in diverse geographic regions worldwide. The overall metrics affirm the substantial advancement achieved with the regional model as an accurate decomposition model representing Southern Africa. Two stations were used as a validation study, as an application example where no localised model was available, and the cluster and regional models both outperformed the comparative decomposition models. This study focused on validating the model for hourly DNI in Southern Africa within a range of K t -intervals from 0.175 to 0.875, and the range could be expanded and validated for future studies. Implementing accurate decomposition models in developing countries can accelerate the adoption of renewable energy sources, diminishing reliance on coal and fossil fuels.

1. Introduction

Photovoltaic (PV) systems require accurate modelling and monitoring to ensure their profitability. The amount of irradiance at the site, the global plane-of-array irradiance (GPI), is the foundation of designing, modelling and monitoring PV systems. The GPI comprises the plane-of-array’s (POA) direct beam, ground and diffuse irradiance components. GPI is used to model and monitor PV systems, as this shows the amount of generated solar power, and, therefore, it is one of the most important contributing factors to designing a PV system. The global horizontal irradiance (GHI), direct normal irradiance (DNI), and diffuse horizontal irradiance (DHI) components are required to calculate these irradiance components.
Irradiance components with a transposition model calculate GPI ( G P O A ) as
G P O A = G B C + G D C + G R C .
G B C is the direct beam irradiance, G R C is the ground-reflected irradiance, and G D C is the diffuse irradiance component in the POA on the collector. The GHI, DNI, and DHI components are required to calculate G B C , G D C and G R C . The sum of the DNI projected onto the horizontal surface using the cosine of the solar zenith angle θ Z and DHI gives the GHI, as shown in Figure 1 [1]:
G H I = D N I · cos θ Z + D H I .
The units of GHI, DHI, and DNI are W / m 2 .
Most ground-based stations at least have measurements of GHI. Other measurements include radiometric data such as DNI, DHI, and ultra-violet and meteorologic data such as the temperature, pressure, rainfall, relative humidity, wind direction and wind speed. Pyranometers measure DHI and GHI, and the pyrheliometer measures DNI.
GHI is measured with a hemispherical view and is mounted horizontally. Similar in setup to other pyranometers, the DHI pyranometer includes the additional feature of being shaded from direct sunlight. The pyrheliometer has a narrow view that only measures the beam directly from the Sun and is usually a Sun tracker for increased accuracy [2]. The irradiance measurements are converted to W / m 2 and logged accordingly.
Calibrating the equipment to the ISO 9060:2018 standard [3] is necessary, and it is advisable to undergo recalibration every two years to ensure the reliability of measurements. The maintenance required is to clean the domes and regularly check and replace the desiccant, which keeps the instruments dry internally.
GHI, DNI, and DHI are interdependent; therefore, having only two irradiance measurements is sufficient to estimate the third using the decomposition models (also sometimes called separation models) [4]. If only the GHI is available, the DNI and DHI also are estimated using the decomposition models. The transposition models calculate GPI using the irradiance components. Therefore, GHI, DHI and DNI correlations are usually empirically expressed as a decomposition model [5].
Indices are relationships between different irradiance components. Decomposition and transposition models utilise these relationships.
The definition of the direct beam transmittance K n and diffuse transmittance K d is
K n = D N I G 0 n ,
K d = D H I G H I .
Ref. [6] defines the K t as
K t = G H I G 0 n cos θ Z .
All K-values ( K t , K n , and K d ) are unitless.
The extraterrestrial irradiance on a normal surface G 0 n depends on the day of the year
G 0 n = ( Solar Constant ) 1 + 0.033 · cos 360 · n 365 .
The solar constant is usually 1367 W / m 2 .
Determining the horizontal extraterrestrial irradiance G 0 h involves multiplying it by the cosine of θ Z as expressed in Equation (7):
G 0 h = G 0 n · cos θ Z .
Multipredictor decomposition models can improve accuracy compared to single predictor models [7]. However, the disadvantage is that multiple measurements must be available, which is not always the case for developing countries or brand-new sites of PV installations.
Refs. [8,9] developed a logistical model to estimate solar diffuse radiation [8,9]. Refs. [10,11,12] have proposed machine-learning-based models to predict solar diffusion and direct components [10,11,12]. Ref. [13] proposed a satellite-based decomposition model as an alternative to ground-based measurements [13,14] proposed statistical models for estimating diffuse radiation [14].
Decomposition models have been developed by assessing previous models and improving the accuracy of these estimations. As more data and measurements become available, researchers have the opportunity to develop models for different climates and temporal resolutions. Most models predominantly use K t . Some of the variables used in the decomposition models are the solar altitude angle β and dew point temperature T d . Using K t as the main predictor in decomposition models is popular because of its simplicity and applicability [7].
Ref. [15] developed a relationship between the K t and K d [15,16] extended the K t - K d relationship to latitudes from 31 to 42° North [16]. Ref. [17] established a GHI and DNI relationship for a Mediterranean site to estimate K n using K t [17].
The Direct Insolation Simulation Code (DISC) was developed by [18]. Refs. [18,19] developed the Dirint model with the hopes of increasing the performance of the DISC model [19]. The Dirint model of [19] has shown superior performance when estimating the DNI [20].
In Korea, ref. [21] developed a model using six Korean locations. The authors of [4,21] developed a new model using [18]’s DISC model by refitting the coefficients [4]. Ref. [22] developed a DNI estimation model using the solar elevation angle for Norway based on hourly GHI and DHI records [22].
Ref. [23] derived K d for Hong Kong [23]. Ref. [24] determined K d using two models with K t and β [24].
The main limitations of decomposition models are that some have a limited climate scope, and the dataset’s temporal resolution affects the irradiance estimation accuracy. A decomposition model in a tropical climate may be unsuitable for a desert climate and vice versa. Intra-hourly-based models perform differently from daily or monthly models, which is why many available decomposition models exist.
The accuracy of decomposition models has been evaluated in several regions, such as Belgium [5], China [25], the USA [20], and North Africa [26].
Ref. [27] provided an extensive study of 140 available decomposition models. The authors state that the predicted DNI’s accuracy highly depends on the decomposition model. Validation studies exist but are limited to a few models and test stations, i.e., biased to a specific location or climate [27]. Research indicates that no decomposition model has been developed and validated for South Africa.
Ref. [28] state that, in general, decomposition models tend to overestimate DHI and underestimate DNI, and typically, models tend to underestimate DHI in overcast periods and overestimate during clear-sky periods [20].
Higher resolution data include higher K t values, resulting in extreme overestimations of DNI. These hourly DNI estimates have higher accuracy than 1-min DNI estimates. Subhourly estimations would be highly beneficial for the real-time monitoring and forecasting of solar power [27].
Figure 2 visualises the testing and validation countries of common decomposition models in green for models such as [15,16,17,24], DISC ([18]), and Dirint ([4,19,21,22,23]) [4,15,16,17,18,19,21,22,23,24].
The decomposition model in South America was developed in Brazil in [29] and in Argentina and Brazil in [30]. Northern African models include Nigeria [31], Algeria [32], and Morocco [33].
Ref. [34] developed a model for Australia and observed that the model only slightly outperformed the Dirint model [34]. The BRL model developed by [9] included a method to construct multiple variable logistic models for the diffuse solar fraction, which includes Mozambique [9]. Figure 2 represents these discussed models [9,29,30,31,32,33,34] in red.
South African research on decomposition models includes the following: ref. [35] published the only Southern African-based study on the relationship between radiation and K t [35]. However, this relationship is with photosynthetically active radiation related to agricultural practices, not PV systems. Clear-sky model assessments and validation studies have been performed by [36,37] for Southern African countries. Clear-sky models simplify atmospheric attenuation to estimate solar irradiance under clear-sky conditions, do not represent decomposition models, and do not include these studies as comparison models, as they are irrelevant to the research.
Ref. [38]’s thesis assessed decomposition and transposition models in South Africa and showed that the models tend to overestimate the DHI but underestimate the DNI [38]. Furthermore, the DISC and Dirint decomposition models showed the most accurate estimations of the DNI and DHI for South African climatic conditions [39].
As discussed, decomposition models are empirical relationships between GHI, DHI and DNI. All three irradiance components are required to estimate GPI. Decomposition models are useful as they reduce the measurement equipment by decomposing one irradiance component into two others; for example, they use GHI to estimate DHI and DNI.
Most decomposition models are not universally applicable and localised to a specific climate, and the temporal resolution is not always transferable. There has not been extensive literature published representing the Southern African region in decomposition models, which this research article will attempt to address.

2. Model Development

The methodology to develop a novel decomposition model is based on selected data from the automated quality control (QC) procedure proposed in [40] and addresses three geographical models:
  • A localised decomposition model, which is site-specific;
  • A clustered decomposition model, which encapsulates several sites to group an area based on their geographical location;
  • A regional (Southern African) model, which encapsulates the data from the SAURAN network for developing a model specific to Southern Africa.
Figure 3 details an overview of the decomposition model development, which is discussed in this section. The first part is ensuring the data are quality-controlled, followed by the development of the decomposition model. The accuracy of the models are rigorously scrutinised to determine if the proposed models outperform the baseline models.

2.1. SAURAN Database

The Southern African Universities Radiometric Network (SAURAN) is a network that includes multiple stations across Southern Africa, collecting meteorological data such as irradiance and wind, among others [41]. Table 1 summarises the SAURAN stations’ corresponding geographical information, such as latitude, longitude, and elevation above sea level.
Table 2 summarises the SAURAN database and dataset sizes, followed by Table 3 which shows the subsequent data points available for model development. Further, the data points assessed are K t between 0.175 and 0.875.
The data points are hourly measurements of the GHI, DNI, and DHI. The split of the train–validation–test datasets is 50:25:25, with the exceptions of two datasets, ILA and MIN. The ILA and MIN have a 0:0:100 data split and are two unknown datasets as part of the test study.
Table 3 also shows each station’s mean GHI, DNI, and DHI, determined after applying the QC procedure.

2.2. Comparison Metrics

The comparison metrics are the root mean square error (RMSE), mean absolute error (MAE), and mean bias error (MBE).
RMSE = i = 1 N ( x i x i ^ ) 2 n , MAE = 1 n i = 1 n | x i x i ^ | , MBE = 1 n i = 1 n ( x i x i ^ ) ,
where x i is the measured value and x i ^ is the predicted value. A low RMSE and MAE indicate a good model, whereas an MBE should be closer to zero. RMSE indicates the concentration of data around the line of best fit. Therefore, a smaller RMSE is indicative of a more accurate model.
The Pearson correlation coefficient r indicates the correlation between data:
r = x i x ¯ y i y ¯ x i x ¯ 2 y i y ¯ 2 .
In Equation (9), x i and y i represent the individual points with index i, and x ¯ and y ¯ represent the mean of the x and y sample set. An r closer to −1 has a negative correlation, meaning if one variable increases, the other decreases. In contrast, if r is closer to 1, it has a positive correlation, meaning if one variable increases, the other would also [42].
Statistical indicators used for the comparison metrics are the MBE, RMSE, and MAE, which are all expressed as a percentage of the mean measured DNI [27] and R 2 . Further comparison metrics are two MAE K t -intervals: K t < 0.60 and K t 0.60 .
The MBE indicates whether a model over- or underestimates the DNI, and the RMSE indicates the deviation of the errors. A significant difference between MAE and RMSE indicates a larger variance in the data. Lower RMSE and MAE are ideal, whereas an MBE closer to zero is optimal. The MAE is an unbiased estimator and also evaluates the two K t intervals. Lower and higher K t indicate overcast and clear-sky conditions, respectively. Therefore, the two K t intervals assess the models under varying weather conditions.

2.3. Regression and Fitting

The relationship between two variables is quantified using statistical methods like regression. Regression techniques can be linear, multi-linear and non-linear.
The definition of a linear relationship is
y = b 0 + b 1 x ,
where y is the response, x is the regressor, b 0 is the intercept, and b 0 is the slope. A regression analysis quantifies the strength of a relationship between y and x [42].
The least squares method estimates b 0 and b 1 so that the sum of the squares of the residuals is at a minimum. The residual sum of squares is denoted as SSE and is the sum of squares of the errors about the regression line. Thus, the minimisation of
SSE = i = 1 n y i y i ^ 2 ,
where y ^ denotes the predicted or fitted value.
The coefficient of determination, R 2 , indicates how good the fit of a model is and is a number between zero and one.
R 2 = 1 y i y i ^ 2 y i y i ¯ 2 .
A higher R 2 value indicates that the model explains the variation in the response variable around its mean, and the regression model fits the observation better [42].
Polynomial regression is the modelling of a dependent, y, as an n t h -degree polynomial of x
y = b o + b 1 x + b 2 x 2 + + b r x 2 .
Exponential regression is where the best fit of an equation is an exponential function, like
y = a + b c x ,
or
y = a + b e c x .
Multi-linear regression has multiple variables, which is the outcome of a response variable
y = β 0 + b 1 x 1 i + b 2 x 2 i + + b k x k i .

2.4. Software Development Tools

The model development utilises a combination of data science applications and modelling. The primary tool is the open-source language Python with the anaconda interface [43], and various available libraries [44,45,46].

2.5. Baseline Models

Three comparative models are used as a baseline to compare the new models. Based on the literature, the DISC and Dirint models performed well for Southern African climates [39,47].
The Dirint [19] and Lee [4] models are also used for comparison because their foundation is similar to the DISC model [18].
Ref. [18]’s DISC quasi-physical approach has three assumptions [4]:
  • The relative air mass ( A M ) is the dominant parameter affecting the relationship between K n and K t ;
  • The physical model used to calculate K n will provide a physically based reference from which the changes in K n can be calculated (see Equation (20) below);
  • Seasonal, annual and climate variations in the relationship between K n and K t are fully accounted for by parametric functions in K t that relate Δ K n to A M , cloud cover, and PW vapour.
A M is defined as [48]:
A M = cos θ Z + 0.5057 · 96.080 θ Z 1.634 1 .
The absolute AM ( A M a ) is the pressurised normalisation of AM, expressed as
A M a = AM P P o ,
where P refers to the atmospheric pressure at the test site and P o is the atmospheric pressure at sea level.
The modelled DNI is determined using Equation (3):
D N I = G 0 n · K n ,
where
K n = K n c Δ K n ,
and
Δ K n = a DISC + b DISC e c DISC · A M .
The clear-sky limit K n c is a polynomial in A M :
K n c = 0.866 0.122 A M + 0.0121 A M 2 0.00653 A M 3 + 0.000014 A M 4 .
Two K t intervals determine the coefficients a DISC , b DISC and c DISC : K t 0.60 and K t > 0.60 .
For K t 0.60
a DISC = 0.512 1.56 K t + 2.286 K t 2 2.222 K t 3 , b DISC = 0.370 + 0.962 K t 3 , c DISC = 0.280 + 0.932 K t 2.048 K t 2 .
For K t > 0.60
a DISC = 5.743 + 21.77 K t 27.49 K t 2 + 11.56 K t 3 , b DISC = 41.4 118.5 K t + 66.05 K t 2 + 31.90 K t 3 , c DISC = 47.01 + 184.2 K t 222.0 K t 2 + 73.81 K t 3 .
Ref. [18]’s model possesses a different functional form because the quasi-physical approach is applied; therefore, it partially reflects the physics involved in the atmospheric transmission of solar radiation [4]. The a DISC , b DISC and c DISC parameters were fitted based on solar radiation data from Atlanta, Georgia, USA, 1981 [18]. Ref. [18] adopted the Bird clear-sky model for K n c (see Equation (22)). The parameters a DISC , b DISC and c DISC , as described in Equations (23) and (24), were then fitted based on the dataset.
The DISC model, termed ‘quasi-physical’, combines a clear-sky model with experimental fits for other sky conditions. The model is a clear-sky irradiance attenuated by a function of K t . The authors of [18] derived the empirical regressions from 12 years of recorded radiation data at 70 stations [5,18].
The Dirint model is based on the DISC model and was developed by [19]. The goal was to improve the accuracy of the DISC model by [18].
The Dirint model uses a clearness index variation parameter K t :
K t = K t 1.031 e 1.4 / ( 0.9 + 9.4 / A M ) + 0.1 .
Furthermore, a stability index parameter Δ K t ,
Δ K t = 0.5 | K t ( i ) K t ( i + 1 ) | + | K t ( i ) K t ( i 1 ) | ,
considers the previous ( i 1 ), current (i) and next hourly ( i + 1 ) record. When the preceding or hourly record is missing, Δ K t is
Δ K t = | K t ( i ) K t ( i ± 1 ) | .
A low Δ K t is a stable condition, whereas a high Δ K t characterises unstable conditions, which allows the distinction between hazy and partly cloudy conditions. The T d is an adequate atmospheric PW estimator [19]. The Dirint model’s atmospheric PW (W) is estimated using
W = exp ( 0.07 · T d 0.075 ) .
The Dirint is a four-dimension conditional model, with the θ Z , K t , Δ K t and W. Based on the four-dimensional model, the calculation of hourly DNI is
D N I = I 0 n · k b e 1.4 / ( 0.9 + 9.4 / A M ) 0.87291 .
where
k b = 0 for K t < 0.2 , k b = a Dirint K t + b Dirint .
Coefficients a Dirint and b Dirint are from a complex lookup table.
Ref. [4] created a new model for Korea with the same format as [18]’s DISC model.
a Lee = 0.342 0.3782 K t , b Lee = 0.5329 + 0.2676 K t 0.0216 k t 2 + 0.1584 K t 3 .
For K t 0.5
c Lee = 0.2117 0.0513 K t + 1.2976 K t 2 3.3222 K t 3
or K t > 0.5
c Lee = 0.7221 10.2801 K t + 30.3285 K t 2 27.9766 K t 3 .
The evaluation consists of comparing the localised, clustered and regional models against the three baseline models: DISC, Dirint and Lee. The DISC and Dirint models were selected based on their performance in estimating DNI for Southern African climates. The Lee and Dirint models have foundational similarities to the DISC model. These models consider whether the newly developed decomposition model improves the accuracy of hourly DNI estimations for Southern Africa. The accuracy evaluation uses the comparison metrics discussed in the next section.

2.6. Decomposition Model Development Methodology

The methodology builds on the DISC model. The DISC model stands out as one of the better-performing models for estimating DNI for South Africa [39]. Its simplicity is evident in its lack of a need for a complex four-dimensional lookup table, unlike the Dirint model.
The original DISC model uses Equation (21), an exponential function. However, the regression model for an exponential function, as discussed in Section 2.3, showed difficulty in finding optimal a, b and c coefficients in all cases. Instead, a second-order polynomial function of A M
Δ K n = a + b · A M + c · A M 2
is a suitable substitute with similar regression results.
The training set then fits a, b and c for intervals K t 0.60 and K t > 0.60 :
a = a 0 + a 1 K t + a 2 K t 2 + a 3 K t 3 , b = b 0 + b 1 K t + b 2 K t 2 + b 3 K t 3 , c = c 0 + c 1 K t + c 2 K t 2 + c 3 K t 3 ,
and the validation and testing sets evaluate the model’s accuracy.
Figure 3 summarises the development for each decomposition model in this study, where each model undergoes the following initial processing steps:
  • Empirical formulae estimate θ Z , A M , pressure, I 0 n , K t , and K n . From this, the assessment of available models aids in developing a new model.
  • Data are split into intervals of 0.05 K t , starting from 0.175 to 0.875.
  • Δ K n is then modelled as Equation (34);
  • The interval or intervals are then fitted against the function to determine Equation (34) to determine the a, b and c coefficients using a least squares regression analysis.
  • From the K t -interval function, the a 0 - a 3 , b 0 - b 3 and c 0 - c 3 coefficients are fitted to a polynomial of Equation (35) with regards to K t .
  • These equations can be used to determine Δ K n and K n , which, in turn, calculate the DNI (see Equations (19) and (20)).
For each SAURAN station, a localised decomposition model is developed. A clustered decomposition model describes an area with similar irradiance patterns using the clustered areas discussed in [40]. Ref. [49] first presented a two-cluster correlation map using the SAURAN database [49], and, by using this approach, this study formulated four clusters instead of two in Southern Africa, as shown in Figure 4a.
Figure 4a shows the clusters’ geographical location, and Figure 4b shows the penetration levels of GHI. Table 4 shows the different clusters’ training sets’ mean GHI, DNI, and DHI.
Cluster 1 receives the most GHI and DNI, and Cluster 3 receives the least, as evident from Figure 4b. The different climates are also evident in these clusters: Cluster 3 is more humid and receives, on average, more DHI than Cluster 1.
Figure 5 shows how the cluster data are combined. Each cluster and the regional (Southern African) model are combined with even distributions of datasets to avoid introducing a bias, as some stations are over-represented in the original dataset. Some stations, such as the SUN, UPR and RVD stations, have considerably more data available as they are either older stations or have not been closed down.
The different stations have varying climates, and therefore, a larger representation of one station will result in a biased model towards that station. The advantage of the even distribution is that every station is sufficiently represented and will not cause a model bias, but this reduces the amount of available data.
Cluster 2’s stations have higher elevation and summer humidity due to its warm, rainy summers and dry, cold winters. The expected annual irradiance levels are lower, as seen in Figure 4b. The stations have higher humidity because of their location and higher DHI levels.
The two stations in Cluster 2, UPR and CSIR, are expected to have more diffuse particles due to the higher air pollution levels and, therefore, higher DHI levels. Cluster 2 has a large bias of the data from Pretoria, South Africa, from the CSIR and UPR datasets.
Cluster 4 has lower annual irradiance levels, as seen in Figure 4b, and FRH and NMU are closer to the coastline, whereas GRT is inland.

3. Development of New Decomposition Models

This section discusses the newly developed a, b and c coefficients of Equation (34).
The section consists of three subsections:
  • The localised decomposition models, developed using the training dataset of the SAURAN station;
  • The clustered decomposition models, which are modelled on the training data of all the stations within the cluster, as discussed in Figure 5;
  • The regional model is modelled on all the stations’ training data (Table 3).

3.1. Localised Decomposition Models

The localised decomposition model equations for the a, b and c coefficients are presented in Appendix A.

3.2. Cluster Decomposition Models

Figure 6, Figure 7, Figure 8 and Figure 9 show the different corresponding clusters’ model coefficients.

3.2.1. Cluster 1

Cluster 1 comprises the HLO, NUST, RVD, SUN and VAN datasets, as shown in Figure 4a and Figure 5.
Figure 6 shows the Cluster 1 and five stations’ a, b and c coefficients. The discussion of the different stations is in Appendix A under Appendix A.5 (HLO), Appendix A.11 (NUST), Appendix A.12 (RVD), Appendix A.13 (SUN) and Appendix A.19 (VAN).
The RVD model is the only model showing difficulty fitting the coefficients with K t . Table 3 indicates that the RVD station has the highest mean DNI and GHI, with the lowest DHI measurements, compared to the rest of Cluster 1’s stations.
The coefficients for Cluster 1 are
a = 2.4134 15.428 K t + 50.7433 K t 2 50.3864 K t 3 for K t < 0.60 , 18.4363 61.7241 K t + 67.5365 K t 2 23.9963 K t 3 for K t 0.60 b = 1.4538 + 13.4628 K t 43.4278 K t 2 + 40.7081 K t 3 for K t < 0.60 , 7.7071 + 23.0993 K t 20.5561 K t 2 + 4.7883 K t 3 for K t 0.60 c = 0.2232 2.1593 K t + 6.6964 K t 2 6.0805 K t 3 for K t < 0.60 , 0.7064 1.9679 K t + 1.5214 K t 2 0.2153 K t 3 for K t 0.60

3.2.2. Cluster 2

Cluster 2 consists of the CSIR, CUT, UBG, UFS, UPR, and UNV datasets. Figure 7 shows the Cluster 2 and six stations’ a, b and c coefficients.
The discussion of the different stations are in Appendix A under Appendix A.1 (CSIR), Appendix A.2 (CUT), Appendix A.14 (UBG), Appendix A.15 (UFS), Appendix A.18 (UPR) and Appendix A.16 (UNV). The UFS have the greatest deviation from the Cluster 2 fit.
The coefficients for Cluster 2 are
a = 1.4234 6.6647 K t + 25.5673 K t 2 27.8939 K t 3 for K t < 0.60 , 10.5939 28.4856 K t + 21.7611 K t 2 3.348 K t 3 for K t 0.60 . b = 0.8033 + 7.723 K t 26.9074 K t 2 + 25.9995 K t 3 for K t < 0.60 , 6.1286 + 15.6674 K t 9.4832 K t 2 0.4949 K t 3 for K t 0.60 . c = 0.1604 1.5968 K t + 5.0219 K t 2 4.537 K t 3 for K t < 0.60 , 0.8966 2.6121 K t + 2.2392 K t 2 0.4802 K t 3 for K t 0.60 .

3.2.3. Cluster 3

Cluster 3 consists of the KZH, KZW and UNZ datasets. Figure 8 shows the Cluster 3 and three stations’ a, b and c coefficients.
The discussion of the different stations is in Appendix A under Appendix A.7 (KZH), Appendix A.8 (KZW) and Appendix A.17 (UNZ). The three models fit quite well and are similar to Cluster 3.
The coefficients for Cluster 3 are
a = 1.7678 9.8995 K t + 34.9076 K t 2 36.2495 K t 3 for K t < 0.60 , 41.1735 157.5062 K t + 201.1335 K t 2 85.5391 K t 3 for K t 0.60 . b = 1.0914 + 10.2302 K t 33.9924 K t 2 + 32.4023 K t 3 for K t < 0.60 , 31.151 + 121.964 K t 158.1413 K t 2 + 67.9737 K t 3 for K t 0.60 . c = 0.2256 2.1961 K t + 6.8092 K t 2 6.1997 K t 3 for K t < 0.60 , 3.4215 13.4513 K t + 17.5277 K t 2 7.5726 K t 3 for K t 0.60 .

3.2.4. Cluster 4

Cluster 4 consists of the NMU, FRH and GRT datasets. Figure 9 shows the Cluster 4 and three stations’ a, b and c coefficients.
The discussion of the different stations is in Appendix A under Appendix A.10 (NMU), Appendix A.3 (FRH) and Appendix A.4 (GRT). The GRT station’s c-coefficient does show difficulty in a fit determination.
The coefficients for Cluster 4 are
a = 0.7671 0.9387 K t + 9.7073 K t 2 14.2827 K t 3 for K t < 0.60 , 20.0495 67.0086 K t + 73.6919 K t 2 26.4669 K t 3 for K t 0.60 . b = 0.2382 + 2.3908 K t 11.3676 K t 2 + 12.3585 K t 3 for K t < 0.60 , 12.5095 + 42.6536 K t 47.1486 K t 2 + 16.8022 K t 3 for K t 0.60 . c = 0.0451 0.5012 K t + 1.8485 K t 2 1.7706 K t 3 for K t < 0.60 , 1.4976 5.1597 K t + 5.8075 K t 2 2.1264 K t 3 for K t 0.60 .

3.3. Regional Decomposition Model

The regional (Southern African) decomposition model data are an even distribution of the SAURAN stations regarding the number of data points used per station. Multiple climates, different elevations, and various pollution levels are represented within the dataset, leading to a better decomposition model for Southern Africa and a regional application.
Figure 10 shows the coefficients a, b and c of the regional model and the four clusters.
The coefficients for the regional model are
a = 1.2893 5.2531 K t + 21.5081 K t 2 24.5156 K t 3 for K t < 0.60 , 19.0295 63.9357 K t + 70.7485 K t 2 25.6524 K t 3 for K t 0.60 . b = 0.6327 + 5.891 K t 21.4431 K t 2 + 21.2593 K t 3 for K t < 0.60 , 11.7813 + 39.923 K t 43.6596 K t 2 + 15.3252 K t 3 for K t 0.60 . c = 0.1118 1.1105 K t + 3.6089 K t 2 3.3222 K t 3 for K t < 0.60 , 1.49 5.2009 K t + 5.9381 K t 2 2.2137 K t 3 for K t 0.60 .

4. Results

Each station is discussed individually by assessing the dataset’s comparison metrics: the R 2 -value, MBE, RMSE, MAE and MAE of two K t -intervals. The results compare the localised, clustered and regional (Southern African) models to the three baseline models, DISC, Dirint and Lee. The tables visualise the results for each station using red and green, with green denoting lower error and red denoting higher error.
Table 3 discusses the validation data. In the previous section, the localised, clustered, and regional models were empirically determined. Appendix A expands on the equations for the localised models.
Section 3.2 and Section 3.3 discussed the clustered and regional models. The test data also introduce two unknown datasets, the ILA and MIN datasets. These datasets assess the models with new data for the developed models. ILA and MIN have no localised model, but geographically, they fall within a cluster: ILA falls under Cluster 1 and MIN under Cluster 2.

4.1. Testing and Validation Results

4.1.1. CSIR

Appendix A.1 shows the decomposition model equations for the CSIR station. Table 5 shows the results of the CSIR station. The results show that the localised, Cluster 2 and regional models outperform the baseline models in all metrics. The localised model significantly improves for lower K t , reducing the MAE from around 60% to 36%.
The test results of the CSIR dataset are presented in Figure A1. As seen in the figure, the localised, cluster, and regional models outperform the baseline models, which is consistent with the validation results of the previous section in Table 5.

4.1.2. CUT

Appendix A.2 shows the decomposition model equations for the CUT station. Table 6 shows the CUT station results. The localised Cluster 2 and regional model significantly improve the comparison metrics over the three baseline models. The Lee model has a similar MBE to the regional model (±0.7) and has a higher K t -metric similar to Cluster 2. However, the Lee RMSE and MAE still do not outperform the new models.
Figure A2 presents the test results of the CUT dataset, where the localised, cluster and regional models outperform the baseline models. The test results are consistent with the validation results presented in Table 6.

4.1.3. FRH

Appendix A.3 shows the decomposition model equations for the FRH station. Table 7 shows the results of the FRH station. The localised model outperforms the baseline models by improving R 2 and MBE and reducing MAE and RSME. The Lee model shows the lowest MAE for higher K t -values; however, it does show an overestimation for DNI with a higher MBE. For most metrics, the localised, Cluster 4 and regional model outperforms the baseline models.
Figure A3 presents the test results of the FRH dataset. The localised Cluster 2 and regional models outperform the baselines, but no significant difference exists between the three new models. The test results presented in Figure A3 correspond with the validation results in Table 7.

4.1.4. GRT

Appendix A.4 shows the decomposition model equations for the GRT station. Table 8 shows the GRT station results. The localised model does show improvement over the DISC and Dirint model but does not significantly outperform the Lee model. The Lee model has a higher R 2 and lower MBE and RMSE, whereas the localised model has a lower MAE for the entire dataset and the two K t intervals. The Cluster 4 and regional models perform better than the DISC and Dirint models but do not significantly outperform all the baseline models.
Figure A4 shows the test results of the GRT dataset. The results correspond with the validation results in Table 8. The localised, Cluster 4 and regional model does outperform the DISC and Dirint model but only marginally outperforms the Lee model.

4.1.5. HLO

Appendix A.5 shows the decomposition model equations for the HLO station. Table 9 shows the HLO station results. The localised model performs better than the baseline models and improves all comparison metrics.
Figure A5 shows the test results of the HLO dataset. The validation results in Table 9 and the test results correspond well, indicating that the localised Cluster 1 and regional models outperform the baseline models.

4.1.6. ILA

Figure A6 presents the test results of the ILA dataset. The ILA dataset has no localised decomposition model; therefore, the testing only assesses the Cluster 1 and regional models. The results show that the Cluster 1 and regional models outperform the baseline models. The results highlight the substitution of using a Cluster model when no localised model is available, subject to the geographical location within the Cluster area.

4.1.7. KZH

Appendix A.7 shows the decomposition model equations for the KZH station. Table 10 shows the results of the KZH station. The localised, Cluster 3 and regional models all show significant improvements in reducing the error over the baseline models. The DISC has a lower MBE than the regional model.
Figure A7 shows the test results of the KZH dataset. The localised, Cluster 3 and regional models all outperform the baseline models. The regional model does not significantly outperform Cluster 3 or the localised model.

4.1.8. KZW

Appendix A.8 shows the decomposition model equations for the KZW station. Table 11 shows the results of the KZW station. The localised, clustered, and regional models show improvement over the baseline models with metrics that assess the entire dataset.
Figure A8 shows the test results of the KZW dataset. The validation and testing results from Table 11 correspond.

4.1.9. MIN

The test results of the MIN dataset are presented in Figure A9. MIN has no localised decomposition model and falls geographically under Cluster 2. The cluster model and localised model show improvement over the baseline models. Much like the ILA dataset, the MIN dataset demonstrates how the clustered and regional models can serve as alternatives to enhance DNI estimations in Southern Africa.
Appendix A.10 shows the decomposition model equations for the NMU station. Table 12 shows the NMU station results. The localised, Cluster 4 and regional models show significant improvement in reducing the errors from the baseline models. Based on the higher MBE, the Cluster 4 and regional models overestimate the DNI more than the DISC and Dirint models.
The test results of the NMU dataset are presented in Figure A10. Localised and cluster models outperform baseline models, which is consistent with the results in Table 12.

4.1.10. NUST

Appendix A.11 shows the decomposition model equations for the NUST station. Table 13 shows the results of the NUST station. The localised model shows superior performance over the baseline models, as well as the clustered and regional models. The metrics of the clustered model compared to the baselines indicate that the regional model slightly overestimates the DNI compared to the lowest baseline model (Lee), which slightly underestimates the DNI.
The test results of the NUST dataset are presented in Figure A11. Localised, clustered and regional models outperform the baseline models, which is consistent with the validation results presented in Table 13. The regional model shows marginal underperformance compared to the localised and Cluster 1 model, but not significant enough to make it unusable.

4.1.11. RVD

Appendix A.12 shows the decomposition model equations for the RVD station. Table 14 shows the RVD station results. The localised, clustered and regional models outperform the baseline models. The Lee model performs better than the regional model but does not outperform the localised and cluster models. The RVD station receives more irradiance on average than other stations in the SAURAN database.
Figure A12 shows the test results of the RVD dataset. The results indicate that the localised, cluster and regional models outperform the baseline models, which is consistent with the validation results of the previous section in Table 14. The localised model’s RMSE is higher than the Lee model; however, the localised model performs best in reducing the error for the other metrics. Though the regional model outperforms the baseline models, it does show the worst performance of the three newly developed models for RVD.

4.1.12. SUN

Appendix A.13 shows the decomposition model equations for the SUN station. Table 15 shows the SUN station results. The localised model outperforms the baseline models by improving R 2 and reducing the MBE, RMSE and MAE. The Cluster 1 and regional models show a slightly worse MBE than the Lee baseline model but otherwise outperform the baseline models. The Lee model also predicts higher K t points with a lower MAE than the regional model; however, the other metrics indicate that the regional model shows better results overall.
Figure A13 shows the test results of the SUN dataset. The results indicate that the localised, cluster and regional models outperform the baseline models, which is consistent with the testing results of the previous section in Table 15. As with the validation results, the regional model is the worst-performing new model, but it still outperforms the baseline models.

4.1.13. UBG

Appendix A.14 shows the decomposition model equations for the UBG station. Table 16 shows the UBG station results. The localised, clustered and regional models all outperform the baseline models. The Lee model has a lower MBE than the Cluster 1 and regional models. The Lee model also has a lower MAE for K t 0.60 ; however, the other metrics indicate that the model does not improve the R 2 , RMSE, overall MAE or K t < 0.60 MAE.
Figure A14 shows the test results of the UBG dataset. The test results show that the localised, Cluster 1 and regional models outperform the baseline models. The Lee model has a lower MBE than the regional model, which is consistent with the validation results in Table 16.

4.1.14. UFS

Appendix A.15 shows the decomposition model equations for the UFS station. Table 17 shows the UFS station results. The localised, Cluster 2 and regional models outperform the baseline models. The Lee model underestimates the DNI slightly better than the Cluster 2 model.
Figure A15 shows the test results of the UFS dataset. All three new decomposition models significantly improve the errors compared to the baseline models, which is consistent with the validation results in Table 17.

4.1.15. UNV

Appendix A.16 shows the decomposition model equations for the UNV station. Table 18 shows the UNV station results. The localised, Cluster 2 and regional models significantly improved over the baseline models. The Cluster 2 and regional model overestimates the DNI more than the DISC model, based on the MBE.
Figure A16 shows the test results of the UNV dataset. The test results correspond with the validation results in Table 18, where the localised, Cluster 2 and regional models outperform the baseline models. The only exception is the MBE, where the Cluster 2 and regional models perform worse than the DISC model. Considering all the metrics, the new models outperform the baselines in reducing the overall error of DNI estimations.

4.1.16. UNZ

Appendix A.17 shows the decomposition model equations for the UNZ station. Table 19 shows the results of the UNZ station. The localised, clustered and regional models all show improvement over the baselines. The Dirint model has a lower MBE than the regional model.
Figure A17 shows the test results of the UNZ dataset, which correspond with the validation results in Table 19.

4.1.17. UPR

Appendix A.18 shows the decomposition model equations for the UPR station. Table 20 shows the UPR station results. The localised, cluster and regional models outperform the baseline models.
Figure A18 shows the test results of the UPR dataset. The comparison metrics of the entire dataset indicate that the localised, cluster and regional models outperform the baseline models, which is consistent with the results of the validation dataset in Table 20.

4.1.18. VAN

Appendix A.19 shows the decomposition model equations for the VAN station. Table 21 shows the results of the VAN station. The localised model outperforms the baseline models by improving R 2 and reducing the MBE, RMSE and MAE. The new models significantly reduce the MAE in the higher K t compared to the baseline models. Even when outperforming the baseline models, the regional model performs worse than the new model. The VAN station receives a very high average DNI and DHI and lower DHI than the rest of the database’s stations. These results are similar to the RVD station, which has significantly higher irradiance levels than the other stations.
Figure A19 shows the test results of the VAN dataset. The new models all outperform the baseline models. The regional model shows the worst performance of the new models, even when outperforming the baseline models, similar to the RVD station that receives more irradiance on average compared to the other stations. These results are consistent with the validation results in Table 21.

4.2. Discussion

Table 22 summarises the performance of the localised, clustered and regional models for both the test and validation sets.
As expected, the localised models outperformed the baseline models for all station datasets because of the site-specific climatic training data. As discussed in the previous section, the cluster model combines multiple stations in a similar geographical area.
The clustered model will have significantly more data from which to train a model. A clustered model is ideal if a site has no data for localised model development using the discussed methodology. The regional (Southern African) model also shows improvement over the baseline models, indicating that this model may be appropriate for adoption as a new model for Southern Africa. The two models with no localised models (ILA and MIN) showed improvement using the clustered and regional models in the validation study.

5. Conclusions

This article presents the development of a new decomposition model of hourly DNI estimations for Southern Africa using the SAURAN database. The new models improved the DISC model in developing new decomposition models localised for Southern African climates. The new decomposition models significantly improved the DNI estimation errors over the baselines for all the SAURAN stations’ validation and test data sets. The proposed methodology can be helpful for the development of local decomposition models for other areas worldwide.
The results indicate that a localised model will improve the estimations of DNI. Clustered models also indicate that grouping data based on similar geographical and climatic properties can also improve the performance of decomposition models. This phenomenon could be helpful when using a clustered decomposition model if no local model or limited data are available but are available from two or more geographically close stations.
The overall model, the regional decomposition model, is encapsulated by different climatic regions and geographical locations. There are also some exceptions where the model over- or underestimates the DNI; however, the overall metrics indicate that the Southern African model significantly improves over the baseline models.
The validation study conducted at two stations further substantiates the superiority of the developed hourly decomposition models over comparative models, emphasising their robustness and applicability in realistic scenarios where the localised models are usually not unavailable. The study validates hourly irradiance data for K t intervals between 0.175 and 0.875. Recommendations for future work include developing models for higher and lower K t values and models for higher temporal resolutions with increased accuracy, which is ideal for the real-time monitoring and short-term forecasting of PV power.
Developing countries with an accurate decomposition model can open the path to expanding the use of renewable energy sources and reducing their dependence on coal and fossil fuels. Good-quality data are needed to ensure the progress of research and development related ot solar energy, and this paper addresses the shortcomings of the Southern African region’s lack of accurate solar resource representation in existing models. The next step is assessing the decomposition models with well-known transposition models to determine improved accuracy for PR estimations.

Author Contributions

Conceptualisation, F.M.D.-D.; methodology, F.M.D.-D.; software, F.M.D.-D.; validation, F.M.D.-D.; formal analysis, F.M.D.-D.; investigation, F.M.D.-D.; resources, F.M.D.-D. and A.J.R.; data curation, F.M.D.-D.; writing—original draft preparation, F.M.D.-D.; writing—review and editing, F.M.D.-D. and A.J.R.; visualisation, F.M.D.-D.; supervision, A.J.R.; project administration, F.M.D.-D. and A.J.R.; funding acquisition, F.M.D.-D. and A.J.R. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The original contributions presented in the study are included in the article, further inquiries can be directed to the corresponding author.

Conflicts of Interest

The authors declare no conflicts of interest.

Appendix A. Localised Decomposition Models

Appendix A.1. CSIR

The a, b and c coefficients for Equation (34) for the CSIR station are
a = 0.6472 1.0883 K t + 12.3429 K t 2 17.7027 K t 3 for K t < 0.60 , 6.1444 9.484 K t 5.138 K t 2 + 9.2371 K t 3 for K t 0.60 . b = 0.5032 + 5.67 K t 22.0707 K t 2 + 22.2654 K t 3 for K t < 0.60 , 2.8669 22.2584 K t + 43.4833 K t 2 24.9682 K t 3 for K t 0.60 . c = 0.1418 1.4904 K t + 4.7805 K t 2 4.3489 K t 3 for K t < 0.60 , 1.1287 + 5.8459 K t 9.4568 K t 2 + 4.8729 K t 3 for K t 0.60 .
The test results of the CSIR dataset are in Figure A1.
Figure A1. Hourly test results of decomposition model development for CSIR.
Figure A1. Hourly test results of decomposition model development for CSIR.
Solar 04 00013 g0a1

Appendix A.2. CUT

The a, b and c coefficients for Equation (34) for the CUT station are
a = 2.4122 15.382 K t + 48.3773 K t 2 45.8766 K t 3 for K t < 0.60 , 15.7689 + 86.1878 K t 141.6902 K t 2 + 73.2142 K t 3 for K t 0.60 . b = 1.6214 + 14.8897 K t 45.5364 K t 2 + 40.6731 K t 3 for K t < 0.60 , 17.2881 85.5492 K t + 134.0146 K t 2 67.4052 K t 3 for K t 0.60 . c = 0.2864 2.6906 K t + 7.8329 K t 2 6.7318 K t 3 for K t < 0.60 , 1.9147 + 9.6635 K t 15.3053 K t 2 + 7.7563 K t 3 for K t 0.60 .
The validation results of the CUT dataset are in Figure A2.
Figure A2. Hourly test results of decomposition model development for CUT.
Figure A2. Hourly test results of decomposition model development for CUT.
Solar 04 00013 g0a2

Appendix A.3. FRH

The a, b and c coefficients for Equation (34) for the FRH station are
a = 1.3199 6.0578 K t + 24.6953 K t 2 27.3825 K t 3 for K t < 0.60 , 33.3782 124.5221 K t + 155.5513 K t 2 64.7203 K t 3 for K t 0.60 . b = 0.5658 + 5.6247 K t 21.3334 K t 2 + 21.2952 K t 3 for K t < 0.60 , 26.2959 + 102.0318 K t 131.6271 K t 2 + 56.3592 K t 3 for K t 0.60 . c = 0.0963 1.024 K t + 3.5042 K t 2 3.2887 K t 3 for K t < 0.60 , 3.9117 15.5175 K t + 20.4822 K t 2 8.9754 K t 3 for K t 0.60 .
The test results of the FRH dataset are in Figure A3.
Figure A3. Hourly test results of decomposition model development for FRH.
Figure A3. Hourly test results of decomposition model development for FRH.
Solar 04 00013 g0a3

Appendix A.4. GRT

The a, b and c coefficients for Equation (34) for the GRT station are
a = 2.9602 19.6392 K t + 58.9738 K t 2 55.4505 K t 3 for K t < 0.60 , 29.0697 + 138.0498 K t 209.2486 K t 2 + 102.455 K t 3 for K t 0.60 . b = 2.2843 + 20.1509 K t 58.384 K t 2 + 51.362 K t 3 for K t < 0.60 , 34.8693 155.9554 K t + 227.6821 K t 2 108.7518 K t 3 for K t 0.60 . c = 0.4656 4.1668 K t + 11.5129 K t 2 9.7101 K t 3 for K t < 0.60 , 8.4664 + 36.6029 K t 51.9382 K t 2 + 24.2327 K t 3 for K t 0.60 .
Figure A4 shows the test results of the GRT dataset.
Figure A4. Hourly test results of decomposition model development for GRT.
Figure A4. Hourly test results of decomposition model development for GRT.
Solar 04 00013 g0a4

Appendix A.5. HLO

The a, b and c coefficients for Equation (34) for the HLO station are
a = 3.1156 21.0151 K t + 62.5312 K t 2 57.2369 K t 3 for K t < 0.60 , 50.8058 194.5525 K t + 248.2624 K t 2 105.3977 K t 3 for K t 0.60 . b = 2.0533 + 18.2104 K t 53.4859 K t 2 + 46.7754 K t 3 for K t < 0.60 , 33.1822 + 128.4612 K t 164.8983 K t 2 + 70.2185 K t 3 for K t 0.60 . c = 0.3508 3.2108 K t + 9.1319 K t 2 7.7849 K t 3 for K t < 0.60 , 3.7247 14.4408 K t + 18.5853 K t 2 7.9382 K t 3 for K t 0.60 .
Figure A5 shows the test results of the HLO dataset.
Figure A5. Hourly test results of decomposition model development for HLO.
Figure A5. Hourly test results of decomposition model development for HLO.
Solar 04 00013 g0a5

Appendix A.6. ILA

There is no decomposition model developed for the ILA station.
The test results of the ILA dataset are in Figure A6.
Figure A6. Hourly test results of decomposition model development for ILA.
Figure A6. Hourly test results of decomposition model development for ILA.
Solar 04 00013 g0a6

Appendix A.7. KZH

The a, b and c coefficients for Equation (34) for the KZH station are
a = 1.3444 6.1333 K t + 23.7139 K t 2 25.5604 K t 3 for K t < 0.60 , 56.4073 216.3047 K t + 276.4264 K t 2 117.5281 K t 3 for K t 0.60 . b = 0.7156 + 6.8854 K t 24.0569 K t 2 + 22.9604 K t 3 for K t < 0.60 , 45.6353 + 178.5713 K t 231.593 K t 2 + 99.6179 K t 3 for K t 0.60 . c = 0.1646 1.6569 K t + 5.2451 K t 2 4.7368 K t 3 for K t < 0.60 , 6.2156 24.5162 K t + 32.0741 K t 2 13.9205 K t 3 for K t 0.60 .
The test results of the KZH dataset are in Figure A7.
Figure A7. Hourly test results of decomposition model development for KZH.
Figure A7. Hourly test results of decomposition model development for KZH.
Solar 04 00013 g0a7

Appendix A.8. KZW

The a, b and c coefficients for Equation (34) for the KZW station are
a = 2.2627 14.4806 K t + 47.1629 K t 2 46.2386 K t 3 for K t < 0.60 , 35.0603 132.558 K t + 167.2728 K t 2 70.2359 K t 3 for K t 0.60 . b = 1.4332 + 13.4317 K t 42.4822 K t 2 + 39.1838 K t 3 for K t < 0.60 , 28.1708 + 109.9004 K t 141.9142 K t 2 + 60.6943 K t 3 for K t 0.60 . c = 0.2678 2.6116 K t + 7.9363 K t 2 7.1072 K t 3 for K t < 0.60 , 3.1759 12.4602 K t + 16.1949 K t 2 6.9726 K t 3 for K t 0.60 .
Figure A8 shows the test results of the KZW dataset.
Figure A8. Hourly test results of decomposition model development for KZW.
Figure A8. Hourly test results of decomposition model development for KZW.
Solar 04 00013 g0a8

Appendix A.9. MIN

There is no decomposition model developed for the MIN station.
The test results of the MIN dataset are in Figure A9.
Figure A9. Hourly test results of decomposition model development for MIN.
Figure A9. Hourly test results of decomposition model development for MIN.
Solar 04 00013 g0a9

Appendix A.10. NMU

The a, b and c coefficients for Equation (34) for the NMU station are
a = 0.0688 + 5.2168 K t 7.122 K t 2 + 0.2163 K t 3 for K t < 0.60 , 15.7683 48.8606 K t + 48.5956 K t 2 15.0438 K t 3 for K t 0.60 . b = 0.3483 2.9436 K t + 3.2778 K t 2 0.0145 K t 3 for K t < 0.60 , 11.0783 + 37.046 K t 39.8763 K t 2 + 13.6687 K t 3 for K t 0.60 . c = 0.0758 + 0.598 K t 1.1265 K t 2 + 0.7025 K t 3 for K t < 0.60 , 1.9782 7.2272 K t + 8.7285 K t 2 3.4834 K t 3 for K t 0.60 .
The test results of the NMU dataset are in Figure A10.
Figure A10. Hourly test results of decomposition model development for NMU.
Figure A10. Hourly test results of decomposition model development for NMU.
Solar 04 00013 g0a10

Appendix A.11. NUST

The a, b and c coefficients for Equation (34) for the NUST station are
a = 2.0021 11.1478 K t + 36.3216 K t 2 35.8544 K t 3 for K t < 0.60 , 17.9917 + 91.7087 K t 144.5222 K t 2 + 72.3519 K t 3 for K t 0.60 . b = 1.2658 + 11.3567 K t 35.6964 K t 2 + 32.4545 K t 3 for K t < 0.60 , 24.1237 111.1616 K t + 165.4905 K t 2 79.997 K t 3 for K t 0.60 , c = 0.2056 1.9277 K t + 5.7745 K t 2 5.0518 K t 3 for K t < 0.60 , 2.59 + 11.8847 K t 17.6199 K t 2 + 8.4894 K t 3 for K t 0.60 .
The test results of the NUST dataset are in Figure A11.
Figure A11. Hourly test results of decomposition model development for NUST.
Figure A11. Hourly test results of decomposition model development for NUST.
Solar 04 00013 g0a11

Appendix A.12. RVD

The a, b and c coefficients for Equation (34) for the RVD station are
a = 4.0191 28.309 K t + 83.0024 K t 2 74.0684 K t 3 for K t < 0.60 , 84.2546 327.2114 K t + 422.3587 K t 2 181.2345 K t 3 for K t 0.60 . b = 2.8895 + 25.4823 K t 74.043 K t 2 + 63.2378 K t 3 for K t < 0.60 , 65.5943 + 256.2577 K t 331.9301 K t 2 + 142.7166 K t 3 for K t 0.60 . c = 0.5456 4.9199 K t + 13.8139 K t 2 11.415 K t 3 for K t < 0.60 , 10.4557 41.161 K t + 53.7847 K t 2 23.3378 K t 3 for K t 0.60 .
The test results of the RVD dataset are in Figure A12.
Figure A12. Hourly test results of decomposition model development for RVD.
Figure A12. Hourly test results of decomposition model development for RVD.
Solar 04 00013 g0a12

Appendix A.13. SUN

The a, b and c coefficients for Equation (34) for the SUN station are
a = 1.4996 7.646 K t + 30.0235 K t 2 33.6216 K t 3 for K t < 0.60 , 32.6225 122.3879 K t + 152.4023 K t 2 62.9667 K t 3 for K t 0.60 . b = 0.835 + 7.9823 K t 28.5353 K t 2 + 28.4965 K t 3 for K t < 0.60 , 19.4403 + 73.2046 K t 90.6277 K t 2 + 36.9662 K t 3 for K t 0.60 . c = 0.1285 1.3014 K t + 4.3214 K t 2 4.0834 K t 3 for K t < 0.60 , 2.0866 7.7979 K t + 9.595 K t 2 3.8902 K t 3 for K t 0.60 .
The test results of the SUN dataset are in Figure A13.
Figure A13. Hourly test results of decomposition model development for SUN.
Figure A13. Hourly test results of decomposition model development for SUN.
Solar 04 00013 g0a13

Appendix A.14. UBG

The a, b and c coefficients for Equation (34) for the UBG station are
a = 0.5601 + 1.4348 K t + 1.9992 K t 2 6.4178 K t 3 for K t < 0.60 , 42.1623 152.8269 K t + 183.6564 K t 2 73.01 K t 3 for K t 0.60 . b = 0.0378 + 0.6077 K t 6.2716 K t 2 + 7.0084 K t 3 for K t < 0.60 , 44.4913 + 167.9403 K t 209.5344 K t 2 + 86.5114 K t 3 for K t 0.60 . c = 0.0159 0.2605 K t + 1.1245 K t 2 0.8877 K t 3 for K t < 0.60 , 9.2841 35.956 K t + 46.1639 K t 2 19.6586 K t 3 for K t 0.60 .
The test results of the SUN dataset are in Figure A14.
Figure A14. Hourly test results of decomposition model development for UBG.
Figure A14. Hourly test results of decomposition model development for UBG.
Solar 04 00013 g0a14

Appendix A.15. UFS

The a, b and c coefficients for Equation (34) for the UFS station are
a = 1.1152 5.4355 K t + 27.3687 K t 2 32.6276 K t 3 for K t < 0.60 , 18.8962 58.0528 K t + 56.4791 K t 2 16.8735 K t 3 for K t 0.60 . b = 0.5439 + 6.5875 K t 27.7113 K t 2 + 28.8812 K t 3 for K t < 0.60 , 19.1711 + 65.1123 K t 71.8324 K t 2 + 25.6916 K t 3 for K t 0.60 . c = 0.1395 1.6147 K t + 5.7356 K t 2 5.4904 K t 3 for K t < 0.60 , 5.238 19.7369 K t + 24.6715 K t 2 10.2415 K t 3 for K t 0.60 .
The test results of the SUN dataset are in Figure A15.
Figure A15. Hourly test results of decomposition model development for UFS.
Figure A15. Hourly test results of decomposition model development for UFS.
Solar 04 00013 g0a15

Appendix A.16. UNV

The a, b and c coefficients for Equation (34) for the UNV station are
a = 1.6679 8.8496 K t + 30.8258 K t 2 31.7801 K t 3 for K t < 0.60 , 17.0947 58.7362 K t + 67.3532 K t 2 25.6036 K t 3 for K t 0.60 . b = 0.9329 + 8.8795 K t 29.6546 K t 2 + 28.1554 K t 3 for K t < 0.60 , 11.0559 + 39.0859 K t 45.1859 K t 2 + 17.0921 K t 3 for K t 0.60 . c = 0.1744 1.744 K t + 5.4271 K t 2 4.9023 K t 3 for K t < 0.60 , 0.7738 2.5106 K t + 2.5888 K t 2 0.8303 K t 3 for K t 0.60 .
The test results of the UNV dataset are in Figure A16.
Figure A16. Hourly test results of decomposition model development for UNV.
Figure A16. Hourly test results of decomposition model development for UNV.
Solar 04 00013 g0a16

Appendix A.17. UNZ

The a, b and c coefficients for Equation (34) for the UNZ station are
a = 1.1129 5.2859 K t + 25.5284 K t 2 30.5998 K t 3 for K t < 0.60 , 33.4179 127.1391 K t + 161.8637 K t 2 68.7679 K t 3 for K t 0.60 . b = 0.6222 + 7.0597 K t 28.1112 K t 2 + 29.4187 K t 3 for K t < 0.60 , 24.1838 + 94.1721 K t 121.4204 K t 2 + 51.9098 K t 3 for K t 0.60 . c = 0.1398 1.5931 K t + 5.5902 K t 2 5.4749 K t 3 for K t < 0.60 , 2.3838 9.2877 K t + 11.9991 K t 2 5.1432 K t 3 for K t 0.60 .
The test results of the UNZ dataset are in Figure A17.
Figure A17. Hourly test results of decomposition model development for UNZ.
Figure A17. Hourly test results of decomposition model development for UNZ.
Solar 04 00013 g0a17

Appendix A.18. UPR

The a, b and c coefficients for Equation (34) for the UPR station are
a = 1.3766 6.4439 K t + 25.7243 K t 2 28.9372 K t 3 for K t < 0.60 , 2.5908 + 24.5466 K t 49.3217 K t 2 + 28.3294 K t 3 for K t 0.60 . b = 0.8308 + 8.2061 K t 29.0534 K t 2 + 28.667 K t 3 for K t < 0.60 , 9.0641 45.9866 K t + 73.7132 K t 2 37.7903 K t 3 for K t 0.60 . c = 0.1626 1.6496 K t + 5.2698 K t 2 4.8501 K t 3 for K t < 0.60 , 1.7989 + 8.4045 K t 12.7021 K t 2 + 6.2414 K t 3 for K t 0.60 .
The test results of the UPR dataset are in Figure A18.
Figure A18. Hourly test results of decomposition model development for UPR.
Figure A18. Hourly test results of decomposition model development for UPR.
Solar 04 00013 g0a18

Appendix A.19. VAN

The a, b and c coefficients for Equation (34) for the VAN station are
a = 1.8649 12.2743 K t + 47.2823 K t 2 50.5904 K t 3 for K t < 0.60 , 7.6031 11.9134 K t 6.9687 K t 2 + 12.3734 K t 3 for K t 0.60 . b = 1.1265 + 12.1632 K t 44.5245 K t 2 + 44.4602 K t 3 for K t < 0.60 , 3.6296 27.7921 K t + 54.161 K t 2 31.162 K t 3 for K t 0.60 . c = 0.2277 2.5008 K t + 8.4732 K t 2 8.0774 K t 3 for K t < 0.60 , 1.3247 + 7.0528 K t 11.5818 K t 2 + 6.025 K t 3 for K t 0.60 .
Figure A19 shows the test results of the VAN dataset.
Figure A19. Hourly test results of decomposition model development for VAN.
Figure A19. Hourly test results of decomposition model development for VAN.
Solar 04 00013 g0a19

References

  1. Masters, G. Renewable and Efficient Electric Power Systems, 2nd ed.; John Wiley and Sons, Inc.: Hoboken, NJ, USA, 2013; pp. 186–245, 316–398. [Google Scholar]
  2. Sengupta, M.; Habte, A.; Wilbert, S.; Gueymard, C.; Remund, J. Best Practices Handbook for the Collection and Use of Solar Resource Data for Solar Energy Applications, 3rd ed.; NREL: Golden, CO, USA, 2021. [Google Scholar] [CrossRef]
  3. ISO 9060:2018; Solar Energy—Specification and Classification of Instruments for Measuring Hemispherical Solar and Direct Solar Radiation. ISO: Geneva, Switzerland, 2018.
  4. Lee, H.; Kim, S.Y.; Yun, C. Comparison of Solar Radiation Models to Estimate Direct Normal Irradiance for Korea. Energies 2017, 10, 594. [Google Scholar] [CrossRef]
  5. Bertrand, C.; Vanderveken, G.; Journée, M. Evaluation of decomposition models of various complexity to estimate the direct solar irradiance over Belgium. Renew. Energy 2015, 74, 618–626. [Google Scholar] [CrossRef]
  6. Liu, B.Y.; Jordan, R.C. The interrelationship and characteristic distribution of direct, diffuse and total solar radiation. Sol. Energy 1960, 4, 1–19. [Google Scholar] [CrossRef]
  7. Yang, D.; Gueymard, C.A. Ensemble model output statistics for the separation of direct and diffuse components from 1-min global irradiance. Sol. Energy 2020, 208, 591–603. [Google Scholar] [CrossRef]
  8. Boland, J.; Huang, J.; Ridley, B. Decomposing global solar radiation into its direct and diffuse components. Renew. Sustain. Energy Rev. 2013, 28, 749–756. [Google Scholar] [CrossRef]
  9. Ridley, B.; Boland, J.; Lauret, P. Modelling of diffuse solar fraction with multiple predictors. Renew. Energy 2010, 35, 478–483. [Google Scholar] [CrossRef]
  10. Soares, J.; Oliveira, A.P.; Božnar, M.Z.; Mlakar, P.; Escobedo, J.F.; Machado, A.J. Modeling hourly diffuse solar-radiation in the city of São Paulo using a neural-network technique. Appl. Energy 2004, 79, 201–214. [Google Scholar] [CrossRef]
  11. Talvitie, T.; Eldali, F.; Pinney, D. Predicting Solar Diffuse and Direct Components Using Deep Neural Networks. In Proceedings of the 2021 IEEE Power & Energy Society Innovative Smart Grid Technologies Conference (ISGT), Washington, DC, USA, 16–18 February 2021; pp. 1–5. [Google Scholar]
  12. Kalyanam, R.; Hoffmann, S. A Novel Approach to Enhance the Generalization Capability of the Hourly Solar Diffuse Horizontal Irradiance Models on Diverse Climates. Energies 2020, 13, 4868. [Google Scholar] [CrossRef]
  13. Bessafi, M.; Oree, V.; Khoodaruth, A.; Chabriat, J.P. Impact of decomposition and kriging models on the solar irradiance downscaling accuracy in regions with complex topography. Renew. Energy 2020, 162, 1992–2003. [Google Scholar] [CrossRef]
  14. Janjai, S.; Phaprom, P.; Wattan, R.; Masiri, I. Statistical models for estimating hourly diffuse solar radiation in different regions of Thailand. In Proceedings of the International Conference on Energy and Sustainable Development: Issues and Strategies (ESD 2010), Chiang Mai, Thailand, 2–4 June 2010; pp. 1–6. [Google Scholar] [CrossRef]
  15. Orgill, J.; Hollands, K. Correlation equation for hourly diffuse radiation on a horizontal surface. Sol. Energy 1977, 19, 357–359. [Google Scholar] [CrossRef]
  16. Erbs, D.; Klein, S.; Duffie, J. Estimation of the diffuse radiation fraction for hourly, daily and monthly-average global radiation. Sol. Energy 1982, 28, 293–302. [Google Scholar] [CrossRef]
  17. Louche, A.; Notton, G.; Poggi, P.; Simonnot, G. Correlations for direct normal and global horizontal irradiation on a French Mediterranean site. Sol. Energy 1991, 46, 261–266. [Google Scholar] [CrossRef]
  18. Maxwell, E. A Quasi-Physical Model for Converting Hourly Global Horizontal to Direct Normal Insolation; Technical Report; Solar Energy Research Institute: Golden, CO, USA, 1987. [Google Scholar]
  19. Perez, R.; Ineichen, P.; Maxwell, E.; Seals, R.; Zelenka, A. Dynamic global-to-direct irradiance conversion models. ASHRAE Trans. 1992, 98, 354–369. [Google Scholar]
  20. Lave, M.; Hayes, W.; Pohl, A.; Hansen, C.W. Evaluation of Global Horizontal Irradiance to Plane-of-Array Irradiance Models at Locations Across the United States. IEEE J. Photovoltaics 2015, 5, 597–606. [Google Scholar] [CrossRef]
  21. Lee, K.; Yoo, H.; Levermore, G.J. Quality control and estimation hourly solar irradiation on inclined surfaces in South Korea. Renew. Energy 2013, 57, 190–199. [Google Scholar] [CrossRef]
  22. Skartveit, A.; Olseth, J.A. A model for the diffuse fraction of hourly global radiation. Sol. Energy 1987, 38, 271–274. [Google Scholar] [CrossRef]
  23. Lam, J.C.; Li, D.H. Correlation between global solar radiation and its direct and diffuse components. Build. Environ. 1996, 31, 527–535. [Google Scholar] [CrossRef]
  24. Reindl, D.; Beckman, W.; Duffie, J. Diffuse Fraction Correlations. Sol. Energy 1990, 45, 1–7. [Google Scholar] [CrossRef]
  25. Yao, W.; Li, Z.; Xiu, T.; Lu, Y.; Li, X. New decomposition models to estimate hourly global solar radiation from the daily value. Sol. Energy 2015, 120, 87–99. [Google Scholar] [CrossRef]
  26. Khalil, S.A.; Shaffie, A. A comparative study of total, direct and diffuse solar irradiance by using different models on horizontal and inclined surfaces for Cairo, Egypt. Renew. Sustain. Energy Rev. 2013, 27, 853–863. [Google Scholar] [CrossRef]
  27. Gueymard, C.A.; Ruiz-Arias, J.A. Extensive worldwide validation and climate sensitivity analysis of direct irradiance predictions from 1-min global irradiance. Sol. Energy 2016, 128, 1–30. [Google Scholar] [CrossRef]
  28. Laiti, L.; Giovannini, L.; Zardi, D.; Belluardo, G.; Moser, D. Estimating Hourly Beam and Diffuse Solar Radiation in an Alpine Valley: A Critical Assessment of Decomposition Models. Atmosphere 2018, 9, 117. [Google Scholar] [CrossRef]
  29. Oliveira, A.P.; Escobedo, J.; Machado, A.; Soares, J. Correlation models of diffuse solar-radiation applied to the city of São Paulo, Brazil. Appl. Energy 2002, 71, 59–73. [Google Scholar] [CrossRef]
  30. Salazar, G.; Pedrosa Filho, M. Analysis of the Diffuse Fraction from Solar Radiation Values Measured in Argentina and Brazil Sites. In Proceedings of the Solar World Congress 2019, Santiago, Chile, 3–7 November 2020. [Google Scholar] [CrossRef]
  31. Chendo, M.; Maduekwe, A. Hourly global and diffuse radiation of Lagos, Nigeria—Correlation with some atmospheric parameters. Sol. Energy 1994, 52, 247–251. [Google Scholar] [CrossRef]
  32. Chikh, M.; Mahrane, A.; Haddadi, M. Modeling the Diffuse Part of the Global Solar Radiation in Algeria. Energy Procedia 2012, 18, 1068–1075. [Google Scholar] [CrossRef]
  33. Benchrifa, M.; Tadili, R.; Idrissi, A.; Essalhi, H.; Mechaqrane, A. Development of New Models for the Estimation of Hourly Components of Solar Radiation: Tests, Comparisons, and Application for the Generation of a Solar Database in Morocco. Int. J. Photoenergy 2021, 2021, 8897818. [Google Scholar] [CrossRef]
  34. Engerer, N. Minute resolution estimates of the diffuse fraction of global irradiance for southeastern Australia. Sol. Energy 2015, 116, 215–237. [Google Scholar] [CrossRef]
  35. Tsubo, M.; Walker, S. Relationships between photosynthetically active radiation and clearness index at Bloemfontein, South Africa. Theor. Appl. Climatol. 2005, 80, 17–25. [Google Scholar] [CrossRef]
  36. Nijegorodov, N. Improved ashrae model to predict hourly and daily solar radiation components in Botswana, Namibia, and Zimbabwe. Renew. Energy 1996, 9, 1270–1273. [Google Scholar] [CrossRef]
  37. Mabasa, B.; Lysko, M.D.; Tazvinga, H.; Zwane, N.; Moloi, S.J. The Performance Assessment of Six Global Horizontal Irradiance Clear Sky Models in Six Climatological Regions in South Africa. Energies 2021, 14, 2583. [Google Scholar] [CrossRef]
  38. Mahachi, T. Energy Yield Analysis and Evaluation of Solar Irradiance Models for a Utility Scale Solar PV Plant in South Africa. Master’s Thesis, University of Stellenbosch, Stellenbosch, South Africa, 2016. [Google Scholar]
  39. Daniel, F.M.; Rix, A.J. The Evaluation of Decomposition Models Under Varying Weather Conditions for South Africa. In Proceedings of the Southern African Sustainable Energy Conference, Stellenbosch, South Africa, 17–19 November 2021; pp. 207–213. [Google Scholar]
  40. Daniel-Durandt, F.M.; Rix, A.J. Automating Quality Control of Irradiance Data with a Comprehensive Analysis for Southern Africa. Solar 2023, 3, 596–617. [Google Scholar] [CrossRef]
  41. SAURAN. 2022. Available online: https://sauran.ac.za/ (accessed on 23 June 2022).
  42. Walpole, R.; Myers, R.; Myers, S.; Ye, K. Probability & Statistics for Engineers & Scientists; Pearson Education Inc.: Boston, MA, USA, 2012. [Google Scholar]
  43. Anaconda Software Distribution. 2020Anaconda Inc. Available online: https://www.anaconda.com/ (accessed on 14 February 2020).
  44. Holmgren, W.; Hansen, C.; Mikofski, M. pvlib python: A python package for modeling solar energy systems. J. Open Source Softw. 2018, 3, 884. [Google Scholar] [CrossRef]
  45. McKinney, W. Data structures for statistical computing in python. In Proceedings of the 9th Python in Science Conference, Austin, TX, USA, 28–30 June 2010; Volume 445, pp. 51–56. [Google Scholar]
  46. Hunter, J.D. Matplotlib: A 2D graphics environment. Comput. Sci. Eng. 2007, 9, 90–95. [Google Scholar] [CrossRef]
  47. Mahachi, T.; Rix, A. Evaluation of irradiance decomposition and transposition models for a region in South Africa Investigating the sensitivity of various diffuse radiation models. In Proceedings of the IECON 2016—42nd Annual Conference of the IEEE Industrial Electronics Society, Florence, Italy, 24–27 October 2016; pp. 3064–3069. [Google Scholar] [CrossRef]
  48. Kasten, F.; Young, A.T. Revised optical air mass tables and approximation formula. Appl. Opt. 1989, 28, 4735–4738. [Google Scholar] [CrossRef] [PubMed]
  49. Farmer, W.; Rix, A.J. Mapping the spatial perturbations seen by the power system network due to intermittent renewable energy sources. In Proceedings of the Southern African Sustainable Energy Conference, Stellenbosch, South Africa, 17–19 November 2021; pp. 222–228. [Google Scholar]
  50. GHI Solar Map ©2017 Solargis. 2021. Available online: https://solargis.com/maps-and-gis-data/download/south-africa (accessed on 27 July 2022).
Figure 1. The irradiance relationships between GHI, DNI, DHI, and θ Z .
Figure 1. The irradiance relationships between GHI, DNI, DHI, and θ Z .
Solar 04 00013 g001
Figure 2. Validation sites of discussed decomposition models.
Figure 2. Validation sites of discussed decomposition models.
Solar 04 00013 g002
Figure 3. Decomposition model development.
Figure 3. Decomposition model development.
Solar 04 00013 g003
Figure 4. Clusters within the Southern African context. (a) SAURAN [41]; (b) GHI across South Africa [50].
Figure 4. Clusters within the Southern African context. (a) SAURAN [41]; (b) GHI across South Africa [50].
Solar 04 00013 g004
Figure 5. Distribution of data within clusters.
Figure 5. Distribution of data within clusters.
Solar 04 00013 g005
Figure 6. Cluster 1 coefficients in Δ K n = a + b · A M + c · A M 2 . (a) a coefficient; (b) b coefficient; (c) c coefficient.
Figure 6. Cluster 1 coefficients in Δ K n = a + b · A M + c · A M 2 . (a) a coefficient; (b) b coefficient; (c) c coefficient.
Solar 04 00013 g006
Figure 7. Cluster 2 coefficients in Δ K n = a + b · A M + c · A M 2 . (a) a coefficient; (b) b coefficient; (c) c coefficient.
Figure 7. Cluster 2 coefficients in Δ K n = a + b · A M + c · A M 2 . (a) a coefficient; (b) b coefficient; (c) c coefficient.
Solar 04 00013 g007aSolar 04 00013 g007b
Figure 8. Cluster 3 coefficients in Δ K n = a + b · A M + c · A M 2 . (a) a coefficient; (b) b coefficient; (c) c coefficient.
Figure 8. Cluster 3 coefficients in Δ K n = a + b · A M + c · A M 2 . (a) a coefficient; (b) b coefficient; (c) c coefficient.
Solar 04 00013 g008aSolar 04 00013 g008b
Figure 9. Cluster 4 coefficients in Δ K n = a + b · A M + c · A M 2 . (a) a coefficient; (b) b coefficient; (c) c coefficient.
Figure 9. Cluster 4 coefficients in Δ K n = a + b · A M + c · A M 2 . (a) a coefficient; (b) b coefficient; (c) c coefficient.
Solar 04 00013 g009
Figure 10. Regional model coefficients in Δ K n = a + b · A M + c · A M 2 . (a) a coefficient; (b) b coefficient; (c) c coefficient.
Figure 10. Regional model coefficients in Δ K n = a + b · A M + c · A M 2 . (a) a coefficient; (b) b coefficient; (c) c coefficient.
Solar 04 00013 g010aSolar 04 00013 g010b
Table 1. SAURAN station summary [40,41].
Table 1. SAURAN station summary [40,41].
Name (Location)Coordinates Lat (°S), Long (°E))Elevation (m)
CSIRCSIR Energy Centre (Pretoria, South Africa)25.747, 28.2791400
CUTCentral University of Technology (Bloemfontein, South Africa)29.121, 26.2161397
FRHUniversity of Fort Hare (Alice, South Africa)32.785, 26.845540
GRTGraaff-Reinet (Graaff-Reinet, South Africa)32.485, 24.586660
HLOMariendal (Mariendal, South Africa)33.854, 18.824178
ILAIlanga CSP Plant (Upington, South Africa)28.490, 21.520884
KZHUniversity of KwaZulu-Natal Howard College (Durban, South Africa)29.871, 30.977150
KZWUniversity of KwaZulu-Natal Westville (Durban, South Africa)29.817, 30.945200
MINCRSES Mintek (Johannesburg, South Africa)26.089, 27.9781521
NMUNelson Mandela University (Gqeberha, South Africa)34.009, 25.66535
NUSTNamibian University of Science and Technology (Windhoek, Namibia)22.565, 17.0751683
RVDRichtersveld (Alexander Bay, South Africa)28.561, 16.761141
SUNStellenbosch University (Stellenbosch, South Africa)33.935, 18.867119
UBGGaborone (Gaborone, Botswana)24.661, 25.9341014
UFSUniversity of Free State (Bloemfontein, South Africa)29.111, 26.1851491
UNVVenda (Vuwani, South Africa)23.131, 30.424628
UNZUniversity of Zululand (KwaDlangezwa, South Africa)28.853, 31.85290
UPRUniversity of Pretoria (Pretoria, South Africa)25.753, 28.2291410
VANVanrhynsdorp (Vanrhynsdorp, South Africa)31.617, 18.738130
Table 2. SAURAN database and dataset sizes from [41].
Table 2. SAURAN database and dataset sizes from [41].
StationDataset SizeStart DateEnd Date
Before QCAfter QC
CSIR46,43426,53911 March 201731 October 2022
CUT28,07714,61924 October 201731 October 2022
FRH40,89522,2337 February 201724 February 2022
GRT18,541  977427 November 201324 January 2016
HLO21,53211,7288 October 201527 October 2020
ILA  8832  467613 October 202131 October 2022
KZH52,32338,8987 December 201507 August 2022
KZW20,29110,7567 December 201512 December 2018
MIN  8185  442328 October 202131 October 2022
MRB  4201  246217 March 201722 October 2019
NMU39,96923,13010 December 201530 September 2022
NUST52,00427,40126 July 201631 October 2022
PMB  9773  541513 July 202131 October 2022
RVD63,71634,45727 March 201428 July 2021
SALT14,151  990821 July 201722 December 2020
STA40,25621,7517 December 201519 April 2021
SUN87,72047,73324 May 201031 October 2022
SUT17159028 February 201720 April 2017
UBG38,91720,64626 November 20146 November 2020
UFS31,66517,15216 January 201430 August 2017
UNV59,10033,14423 April 201531 October 2022
UNZ56,39930,37311 July 201431 October 2022
UPR78,79242,12819 September 201331 October 2022
VAN24,70113,23426 August 201610 July 2019
Table 3. Model development stations indicating the mean GHI, DNI, and GHI and sizes of training, validation, and testing sets.
Table 3. Model development stations indicating the mean GHI, DNI, and GHI and sizes of training, validation, and testing sets.
StationMean 1Dataset 2Cluster Allocation
GHI
[W/m2]
DNI
[W/m2]
DHI
[W/m2]
TotalTrainValidationTest
CSIR57559916714,9917495374837482
CUT60963915991614580229022912
FRH54458315112,2246112305630564
GRT57362415157882894144714474
HLO55060813870613530176517661
ILA58968013127090027091
KZH53351717987824391219521963
KZW53151118459452972148614873
NMU55654516510,5625281264026414
NUST61467014915,9017950397539761
MIN56457316127610027612
RVD63072912519,6249812490649061
SUN55664513328,50814,254712771271
UBG59160215812,1376068303430352
UFS56765413710,2575128256425652
UNV57952419715,8747937396839692
UNZ53052817610,0555027251425143
UPR56860916328,08914,044702270232
VAN59768312678603930196519651
1 Daylight average, 2 Dataset size after quality control as in [40] and 0.175 ≤ Kt ≤ 0.875.
Table 4. Mean cluster irradiances.
Table 4. Mean cluster irradiances.
Mean 4
GHIDNIDHI
[ W / m 2 ][ W / m 2 ][ W / m 2 ]
Cluster 1592669135
Cluster 2583604165
Cluster 3534523178
Cluster 4557579158
4 Mean values of training set.
Table 5. Hourly validation results of a decomposition model development for CSIR.
Table 5. Hourly validation results of a decomposition model development for CSIR.
Solar 04 00013 i001
Table 6. Hourly validation results of decomposition model development for CUT.
Table 6. Hourly validation results of decomposition model development for CUT.
Solar 04 00013 i002
Table 7. Hourly validation results of decomposition model development for FRH.
Table 7. Hourly validation results of decomposition model development for FRH.
Solar 04 00013 i003
Table 8. Hourly validation results of decomposition model development for GRT.
Table 8. Hourly validation results of decomposition model development for GRT.
Solar 04 00013 i004
Table 9. Hourly validation results of decomposition model development for HLO.
Table 9. Hourly validation results of decomposition model development for HLO.
Solar 04 00013 i005
Table 10. Hourly validation results of the decomposition model development for KZH.
Table 10. Hourly validation results of the decomposition model development for KZH.
Solar 04 00013 i006
Table 11. Hourly validation results of decomposition model development for KZW.
Table 11. Hourly validation results of decomposition model development for KZW.
Solar 04 00013 i007
Table 12. Hourly validation results of decomposition model development for NMU.
Table 12. Hourly validation results of decomposition model development for NMU.
Solar 04 00013 i008
Table 13. Hourly validation results of decomposition model development for NUST.
Table 13. Hourly validation results of decomposition model development for NUST.
Solar 04 00013 i009
Table 14. Hourly validation results of decomposition model development for RVD.
Table 14. Hourly validation results of decomposition model development for RVD.
Solar 04 00013 i010
Table 15. Hourly validation results of decomposition model development for SUN.
Table 15. Hourly validation results of decomposition model development for SUN.
Solar 04 00013 i011
Table 16. Hourly validation results of decomposition model development for UBG.
Table 16. Hourly validation results of decomposition model development for UBG.
Solar 04 00013 i012
Table 17. Hourly validation results of decomposition model development for UFS.
Table 17. Hourly validation results of decomposition model development for UFS.
Solar 04 00013 i013
Table 18. Hourly validation results of decomposition model development for UNV.
Table 18. Hourly validation results of decomposition model development for UNV.
Solar 04 00013 i014
Table 19. Hourly validation results of the decomposition model development for UNZ.
Table 19. Hourly validation results of the decomposition model development for UNZ.
Solar 04 00013 i015
Table 20. Hourly validation results of decomposition model development for UPR.
Table 20. Hourly validation results of decomposition model development for UPR.
Solar 04 00013 i016
Table 21. Hourly validation results of decomposition model development for VAN.
Table 21. Hourly validation results of decomposition model development for VAN.
Solar 04 00013 i017
Table 22. Summary of test and validation sets of stations outperforming baseline models.
Table 22. Summary of test and validation sets of stations outperforming baseline models.
DatasetLocalised Model Outperforms Baseline ModelsCluster Model Outperforms Baseline ModelsBaseline ModelsRegional Model Outperforms Baseline Models
TestValidationTestValidationTestValidation
CSIR
CUT
FRH
GRT
HLO
ILA---
KZH
KZW
MIN---
NMU
NUST
RVD
SUN
UBG
UFS
UNV
UNZ
UPR
VAN
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.

Share and Cite

MDPI and ACS Style

Daniel-Durandt, F.M.; Rix, A.J. New Decomposition Models for Hourly Direct Normal Irradiance Estimations for Southern Africa. Solar 2024, 4, 269-306. https://doi.org/10.3390/solar4020013

AMA Style

Daniel-Durandt FM, Rix AJ. New Decomposition Models for Hourly Direct Normal Irradiance Estimations for Southern Africa. Solar. 2024; 4(2):269-306. https://doi.org/10.3390/solar4020013

Chicago/Turabian Style

Daniel-Durandt, Francisca Muriel, and Arnold Johan Rix. 2024. "New Decomposition Models for Hourly Direct Normal Irradiance Estimations for Southern Africa" Solar 4, no. 2: 269-306. https://doi.org/10.3390/solar4020013

APA Style

Daniel-Durandt, F. M., & Rix, A. J. (2024). New Decomposition Models for Hourly Direct Normal Irradiance Estimations for Southern Africa. Solar, 4(2), 269-306. https://doi.org/10.3390/solar4020013

Article Metrics

Back to TopTop