Next Article in Journal
Energy Poverty and Its Indicators: A Multidimensional Framework from Literature
Previous Article in Journal
Developing the Design of Single-Axis Sun Sensor Solar Tracking System
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Prediction of Global Solar Irradiance on Parallel Rows of Tilted Surfaces Including the Effect of Direct and Anisotropic Diffuse Shading

1
Institute of Earth Sciences, University of Évora, Rua Romão Ramalho 59, 7000-671 Évora, Portugal
2
Department of Mechatronics Engineering, University of Évora, Rua Romão Ramalho 59, 7000-671 Évora, Portugal
3
Physics Department, University of Évora, Rua Romão Ramalho 59, 7000-671 Évora, Portugal
*
Author to whom correspondence should be addressed.
Energies 2024, 17(14), 3444; https://doi.org/10.3390/en17143444
Submission received: 6 June 2024 / Revised: 5 July 2024 / Accepted: 10 July 2024 / Published: 12 July 2024
(This article belongs to the Topic Advances in Solar Technologies)

Abstract

:
Solar photovoltaic power plants typically consist of rows of solar panels, where the accurate estimation of solar irradiance on inclined surfaces significantly impacts energy generation. Existing practices often only account for the first row, neglecting shading from subsequent rows. In this work, ten transposition models were assessed against experimental data and a transposition model for inner rows was developed and validated. The developed model incorporates view factors and direct and circumsolar irradiances shading from adjacent rows, significantly improving global tilted irradiance (GTI) estimates. This model was validated against one-minute observations recorded between 14 April and 1 June 2022, at Évora, Portugal (38.5306, −8.0112) resulting in values of mean bias error (MBE) and root-mean-squared error (RMSE) of −12.9 W/m2 and 76.8 W/m2, respectively, which represent an improvement of 368.3 W/m2 in the MBE of GTI estimations compared to the best-performing transposition model for the first row. The proposed model was also evaluated in an operational forecast setting where corrected forecasts of direct and diffuse irradiance (0 to 72 h ahead) were used as inputs, resulting in an MBE and RMSE of −33.6 W/m2 and 169.7 W/m2, respectively. These findings underscore the potential of the developed model to enhance solar energy forecasting accuracy and operational algorithms’ efficiency and robustness.

1. Introduction

In recent decades, there has been a significant surge in the installed capacity of PV systems. In 2022 alone, solar photovoltaics comprised two-thirds of the new renewable energy capacity added to the grid, totaling 239 GW [1]. This growth is driven by global efforts to decarbonize the economy and achieve net-zero greenhouse gas emissions by 2050. The International Energy Agency (IEA) estimates that PV systems will continue to expand, with new capacity projected to exceed 700 GW by 2028 [2]. With the rapid increase in the installed capacity of solar energy systems worldwide and the characteristic variability of solar radiation, accurate estimation and forecasting of power output is becoming increasingly important. The factor that most affects this output is the irradiance incident on the solar energy collectors. However, weather stations, satellites, and numerical weather prediction (NWP) models usually provide global (GHI) and diffuse (DIF) irradiance data on the horizontal plane and direct irradiance on a plane normal to the sun’s rays (DNI). Therefore, transposition models that can compute the irradiance incident on a tilted surface from the available variables are essential.
Various transposition models have been developed and studied over the years. They can be divided into analytical, semi-empirical, and empirical models [3]. Analytical models are based on the laws of physics and only require the geometric characteristics and position data of the surface, while semi-empirical and empirical models also rely on observation data.
Physical models are usually based on the sum of the direct, diffuse, and reflected irradiances on the tilted surface. The computation of the direct and reflected components is often the same for all models, with the computation of the diffuse irradiance being the distinct factor [4]. Among these models, some assume a uniform sky dome radiance with the same intensity in all directions, termed isotropic. Anisotropic models, on the other hand, define indices representing the irradiance intensity from different regions of the sky dome, such as the circumsolar region (an annular region surrounding the sun disk), the horizon-brightening region (a band along the horizon), and the background or remaining region often called the isotropic region.
Transposition models also tend to show varying degrees of performance depending on sky conditions, whether clear, cloudy, or overcast [5]. The complexity and variability of the shape and position of clouds make the accurate instantaneous evaluation of the diffuse component extremely difficult without knowing the distribution of sky radiance. A model considering this would be more complex and require not only the simple measurement of the irradiance on the tilted plane but many more variables that are not commonly measured, such as the spectral radiance. Thus, physical models must rely on assumptions that can result in high transient errors [6]. Some of the most referenced analytical physical models in the literature are reviewed in this work, such as the Klucher [7], Hay–Davies [8], and Modified Bugler [9] models. On the other hand, abundant research on the comparison of the different models at different locations and positions of the surfaces has been published [3,4,10,11]. From these, no single model emerges as the best for all studied locations/positions, but it has been widely demonstrated that anisotropic models tend to outperform isotropic models when compared to measured data because they consider the angular dependence of sky radiance, thus reflecting real-world conditions more closely compared to isotropic models.
Semi-empirical and empirical models, including machine learning approaches which are data-driven, tend to be site-dependent, which means that they are biased towards specific locations and/or tilt and azimuth angles [12,13,14,15,16], and will not be considered here. While these models are gaining relevance and can provide substantial benefits in some situations, their applicability is often limited because extensive local experimental data are required for the model calibration. This increases the challenge and difficulty of generalizing them to different geographic locations or environmental conditions. Given the focus of this study on developing a transposition model that is applicable across various conditions without the need for site-specific parameter adjustments, these models were excluded to maintain a broader applicability and robustness of proposed approach.
Most transposition models have been developed for a tilted surface in an open field. However, a solar energy power plant, such as photovoltaic power plants, comprises several rows of modules, where rows other than the first have a different view of the sky dome and the ground compared to the first row. Therefore, when computing the global tilted irradiance (GTI) on a panel not located in the front row, the transposition model should be adjusted to account for the fraction of the sky dome obscured by the rows at the front, as well as the possible blocking of the direct beam and circumsolar irradiance from the sun disk and the reflected irradiance coming from the front row and the ground between them. Since in large solar power plants the number of panels in the first row can be much smaller than the number of panels in the remaining rows, having a more accurate transposition model for these panels can improve the estimation of power output. Few authors have focused on this aspect, with the most notable and recent works being those by Applebaum et al. [17], Varga and Mayer [18], and Tschopp et al. [19].
As mentioned above, one aspect of higher complexity in transposition models is the parametrization of diffuse irradiance. This component of the GTI is strongly related to the diffuse horizontal irradiance (DIF) and the global horizontal irradiance (GHI), which are typically measured in radiometric stations. However, if only GHI is measured, a separation method is needed to estimate DNI and DIF. If there are no measurements at the location of interest, estimations or forecast values need to be used instead. This increases the uncertainty associated with the estimation of GTI and, consequently, of power output. Of particular interest in this work is the usage of forecast values of solar radiation components, with the goal of predicting power output in a time horizon up to three days with an acceptable accuracy. In this regard, when utilizing forecast DNI and DIF as inputs, a transposition model allows for the forecasting of GTI and power output of photovoltaic systems if coupled with thermal and electric models of photovoltaic systems.
These forecasts can be obtained from a global operational numerical weather prediction model (NWP) such as the Integrated Forecasting System (IFS) of the European Centre for Medium-range Forecasts (ECMWF) [20]. By incorporating constitutive and state equations describing physical phenomena in the atmosphere, subject to boundary and initial conditions, NWP models provide the evolution of the atmospheric state, which include surface-level GHI. In the case of IFS/ECMWFS, DNI is also included, allowing for the study and assessment of solar resources and operational forecasts [21,22].
This work presents a comprehensive review of transposition models, categorizing models that do not include shading or irradiance masking as “first-row models”, and those that include these aspects as “inner-row models”. The development of a transposition model for inner rows is also presented, which can be applied to a desired first-row model. The proposed model is validated using experimental values, which are also used for the assessment and comparison of all models. The use of operational forecasts of solar irradiance is assessed to estimate the impact of the accuracy of the transposition model on GTI prediction at different time horizons. The primary aim is to develop a transposition model for inner rows of solar power plants, validate the proposed model against experimental measurements, and evaluate its performance using operational solar radiation forecasts. This study addresses gaps in the existing literature, which mainly focuses on the first row of panels, neglecting the impact of direct and anisotropic diffuse shading on inner rows. The hypotheses focus on the expected improvement in the accuracy of the developed model’s results over existing models for the first row, as well as its robustness when used with forecasted irradiance data as input.
This paper is organized as follows: An initial review of the most used analytical models and recent models developed for rows that are not the first is presented in Section 2. Section 3 details the development of a transposition model based on the works of Tschopp et al. and Varga and Mayer [18,19] for surfaces not located in the first row. Section 4 presents the experimental setup and procedure for testing the different models, and Section 5 discusses the evaluation and results. Section 6 establishes a connection with solar irradiance forecasting by applying the developed model to forecasted values of direct normal and diffuse horizontal irradiance, integrating it as an essential part of an operational algorithm for forecasting solar irradiance on a tilted surface. Finally, Section 7 presents the conclusions of this work.

2. Solar Irradiance Transposition Models for Tilted Surfaces

This section provides a review of analytical transposition models, ranging from the most widely used and commonly referenced transposition models that do not consider shading and obscuration to the few and most recently proposed models designed for inner rows in solar power plants.

2.1. Transposition Models for the First Row of Collectors in Solar Power Plants

For modeling solar irradiance on a tilted surface (GTI), nine of the most commonly used analytical transposition models are reviewed, where one is isotropic and the remaining are anisotropic models. In all the following equations, α is the solar height, Φ is the solar zenith angle, γ s is the solar azimuth, θ is the incidence angle on the surface, β is the tilt angle of the surface, and γ p is the azimuth of the surface, taking the horizontal plane and the local meridian as references.
Analytical transposition models include the computation of tilt factors for determining the beam ( R b ), diffuse ( R d ), and reflected ( R r ) irradiances on the tilted surface. These factors are used in Equation (1) for determining G T I , where D H I is the direct horizontal irradiance as defined in Equation (2):
G T I = R b D H I + R d D I F + R r G H I
D H I = D N I × c o s ( Φ )
For the first row of photovoltaic panels, assuming no obstruction of the sun’s rays, the direct irradiance on the tilted plane can be calculated at any given instant by multiplying D H I with the beam radiation tilt factor described in Equation (3):
R b = c o s ( θ ) c o s ( Φ )
The irradiance reflected onto the tilted surface by the surrounding surfaces is assumed isotropic and is modeled by multiplying G H I with the reflected irradiance tilt factor R r (Equation (4)), which consists of the product of the albedo of the surrounding ground ( ρ g ) and the view factor F g (Equation (5)).
R r = ρ g F g ,
F g = 1 c o s ( β ) 2
Numerous models have been proposed for determining the diffuse tilt factor R d since the diffuse radiation is highly variable due to its dependency on atmospheric conditions. Models that consider the diffuse irradiance on a tilted surface coming from the entire sky dome to have the same intensity are isotropic, while others that divide the sky dome into regions with different intensities of diffuse radiance are called anisotropic models. These regions are usually the circumsolar region, which constitutes an annular region surrounding the sun, i.e., the horizon-brightening region, which is the region of the sky dome near the horizon, and the remaining sky dome region is termed the isotropic region.
In the following, several models are presented, including the widely used isotropic model that was developed by Liu and Jordan [23]. It assumes that the diffuse solar radiation is uniformly spread across the sky, which simplifies the calculation of solar radiation on inclined surfaces by considering the sky as a homogenous light source. This model only considers the tilt angle of the surface as input, as presented in Equations (6) and (7), where the diffuse tilt factor is simply a view factor between the surface and the sky. Its advantages lie in its simplicity and widespread adoption for preliminary solar energy designs and studies. However, the model’s accuracy can be compromised due to its isotropic assumption, leading to potential inaccuracies in non-homogenous sky conditions and environments with significant obstructions. Despite these limitations, it remains a significant tool in solar energy research and engineering, often serving as a benchmark for more complex anisotropic models that account for the directional distribution of diffuse radiation.
Liu and Jordan model (LJ):
R d = F s
F s = 1 + c o s ( β ) 2
The remaining models analyzed are anisotropic models, which are often more representative of real sky conditions.
Bugler model (B):
The Bugler model [24] includes a circumsolar region in which its contribution to the diffuse irradiance is assumed to be 5% of DNI and can be computed using Equation (8), where F s is the sky view factor given by Equation (7). Bugler derived this model by analyzing solar radiation data and developing equations that incorporate these anisotropic factors, resulting in more accurate predictions of solar irradiance on tilted surfaces. However, this model overlooks the contribution of the circumsolar region to the isotropic irradiance and was later modified.
R d = F s + 0.05 D H I × R b D I F
Temps and Coulson model (TC):
The Temps and Coulson [25] model, defined by Equation (9), considers the circumsolar region through the factor 1 + c o s 2 θ s i n 3 Φ and the horizon brightening through 1 + s i n 3 β 2 . This model was obtained through measurements of direct and diffuse solar flux incident on slopes of various orientations for clear-sky conditions through a pyranometer.
R d = F s 1 + c o s 2 θ s i n 3 Φ 1 + s i n 3 β 2
Klucher model (K):
The Klucher model [7], given by Equation (10), is similar to the Temps and Coulson model, except for the inclusion of the clearness index f , defined by Equation (11). Under overcast conditions, this model reduces to the isotropic model (where f tends to 0). It was derived from a 6-month dataset of measured hourly diffuse and total solar radiation on a horizontal plane and total radiation on surfaces with 0° azimuth and tilts of 37° and 60°.
R d = F s 1 + f s i n 3 β 2 1 + f c o s 2 θ s i n 3 Φ
f = 1 D I F D H I + D I F 2
Hay–Davies model (HD):
The Hay–Davies model [8], which does not consider horizon brightening, is given by Equation (12). Here, A is named the anisotropy index, also termed direct or beam irradiance index, and is defined in Equation (13), where E N I is the extraterrestrial normal irradiance. This index represents the total transmittance of the atmosphere for beam radiation and is used to define the portion of circumsolar ( A × R b ) and isotropic ( F s 1 A ) diffuse irradiance on the tilted surface.
R d = A × R b + F s 1 A
A = D N I E N I
Ma Iqbal model (MI):
The Ma Iqbal model [26], described through Equation (14), is similar to the Hay–Davies model but uses the clearness index k T (Equation (15)) to define the circumsolar portion of diffuse irradiance instead of the direct irradiance index. The clearness index is computed as the ratio of global horizontal irradiance (GHI) to extraterrestrial horizontal irradiance (EHI). This formulation allows the Ma Iqbal model to account for varying sky conditions and solar geometry, enhancing its applicability in solar energy studies and irradiance forecasting.
R d = k T R b + F s ( 1 k T )
k T = G H I E H I
Modified Bugler model (MB):
The Modified Bugler model [9] improves the Bugler model by considering the contribution of the circumsolar region to the background isotropic diffuse radiance, computed through Equation (16). This model was validated using extensive datasets from various locations around the world. These datasets included measured solar radiation data collected under different sky conditions, such as clear sky, cloudy sky, and overcast sky conditions. Additionally, the model was validated against experimental data that accounted for various surface orientations, including horizontal, tilted, and vertical surfaces.
R d = 1 0.05 D H I D I F F s + 0.05 D H I × R b D I F
Modified Ma Iqbal model (MMI):
The modified Ma Iqbal model (Equation (17)) [3] is similar to the original model proposed by this author but uses a clearness index (Equation (18)) adjusted with an optical air mass as defined in Equation (19) [27].
R d = k T R b + ( 1 k T ) F s   a = 1 ,
k T = k T 1.031 e x p 1.4 0.9 + 9.4 M + 0.1
M = 1 cos Φ + 0.15 ( 93.885 Φ ) 1.253
Reindl model (R):
The Reindl model [28], based on the Hay–Davies model, uses the same anisotropic index, A (Equation (13)), but includes an additional term to account for horizon brightening (Equation (20)). This model was developed based on datasets from locations across the USA.
R d = A × R b + F s 1 A 1 + D H I D H I + D I F s i n 3 β 2
Perez model (P):
The Perez model, one of the most widely used transposition models in the literature, is included in this work even though it cannot be classified as an analytical model since it is built upon measured data and is subject to continuous revisions for this reason. The version of the Perez model [29] used in this work, often referred to as Perez3, is considered the most widely accepted version of this model [4]. This model categorizes the sky in three zones: the circumsolar zone, horizon band, and isotropic zone. The diffuse irradiance tilt factor is computed using Equations (21)–(27), with the assumption that the circumsolar irradiance originates from a point source and irradiance from the horizon from an infinitesimally thin region. The zenith angle is exceptionally used in radians for Equation (26). The coefficients F i j were determined using data from 10 American and 3 European locations and can be consulted in [4,29].
R d = F 1 a b + 1 F 1 F s + F 2 sin β
a = m a x ( 0 , cos θ )
b = m a x cos 85 ° , sin 90 ° θ z
F 1 = m a x 0 , F 11 ε + F 12 ε + F 13 ε θ z
F 2 = F 21 ε + F 22 ε + F 23 ε θ z
ε = G H I / D I F + 1.041 θ z 3 1 + 1.041 θ z 3
= D I F E N I cos θ z

2.2. Transposition Models for Surfaces Not Located in the Front Row of a Solar Power Plant

The models presented in Section 2.1 were developed for a single row. In the case of photovoltaic power plants, for instance, which are composed of various parallel rows, there is a partial obscuring of the sky radiance and, eventually, the blocking of direct sun rays by the front rows over the second and subsequent rows, thus affecting the global tilted irradiance on these surfaces. An obscuring of the ground between rows can also occur, which affects the reflected irradiance.
Some authors have considered this, namely Appelbaum et al. [17], who adjusted the Klucher transposition model [7] by computing the sky view factor for the second row and adjusting the indices of the circumsolar and horizon-brightening diffuse regions. For this, the authors considered the circumsolar irradiance coming from a point source (the sun) which would be completely obscured when the obscuring angle from the sun to the base of the second row caused by the front row is higher than the sun elevation angle. The correction for the horizon-brightening region was considered near sunrise and sunset. This adjusted model was validated with data obtained in Tel Aviv, showing that for rows other than the first, the effect of the anisotropic region defined as horizon brightening is negligible.
Varga and Mayer [18] modified the Hay–Davies model [8], which does not include a horizon-brightening region, to calculate the distribution of solar irradiance along a tilted surface for rows behind the first. Besides the adjustment of the view factors, the circumsolar irradiance was corrected through the modeling of the fraction of the region obscured by the front row by considering a circumsolar region with an apparent angular radius of 15°, and the ground was divided into two sections: a shaded section that reflects only isotropic diffuse irradiance and an unshaded section that reflects global irradiance. The model was validated against power generation data from Hungary showing a clear improvement when compared to a model that only considers shading of the direct irradiance.
Tschopp et al. [19] also adjusted the Hay–Davies model for surfaces that are not on the front row by dividing the length of the surface of the row being evaluated, the back surface of the row at its front, and the ground between them in segments and then computing diffuse and direct irradiance in each segment. For this, the view factors between all segments and between each segment and the sky are computed, while the obscuring of the direct and circumsolar irradiances are determined based on a simple formula for the shadow height on the surface. The model showed a significant improvement in the estimation of GTI for rows other than the first compared to the original transposition model.
Considering this, it is expected that combining the higher detail in the modeling of shadow and circumsolar irradiance obscuration presented in [18] and the inclusion of the back surface of the row at the front and the ground between rows presented in [19] will provide better estimations of GTI for rows that are not at the front of a solar power plant.

3. Development of Transposition Model for Surfaces Not Located in the Front Row of a Solar Power Plant

The model proposed in this work for determining solar irradiance on tilted rows adjacent to the first row builds upon the models presented in [18,19].Similar to these models, it assumes rows of panels with lengths much greater than their heights, resulting in a 2D representation as commonly used in the literature [17,18,19,30,31].
While the aforementioned models only consider cases when the sun is positioned at the front of the rows, here, the modeling of the direct and circumsolar irradiance shading for all involved surfaces is also included, considering any position of the sun, resulting in a more realistic model. Additionally, the developed model was made generic and can be applied to any first-row transposition model as long as it clearly considers direct, circumsolar, and isotropic irradiance, instead of relying on a predefined first-row model. Extra detail was also included in the modeling of ground shading by considering rows beyond the two main rows being modeled.
The three surfaces considered in this model, namely the front of the panel being evaluated, the back of the panel at its front, and the ground between rows, are divided into segments, as shown in Figure 1, where L , D , h 0 , and β are the length of the panels, the horizontal distance between rows, the vertical distance between the ground and the base of the rows, and the tilt angle, respectively. The value of G T I is computed for each segment c of the panel being evaluated and considers the reflected solar irradiance from each segment of the ground u and the back of the front panel v . The number of segments into which these surfaces are divided is defined by the user and can be adjusted, namely the number of segments of the panel being evaluated i = 1 : n , the number of segments of the ground j = 1 : p , and the number of segments of the back of the front panel k = 1 : q .
For each segment of both panels and ground, the G T I is computed as in [19] using Equation (28), where G T I , I r , and D are vectors with the values of the global tilted irradiance, the direct normal irradiance, and the diffuse horizontal irradiance, respectively, on each segment of the panel, ground, and front panel. I is the identity matrix while F is the view factor matrix between all segments and R is the reflectivity matrix for each segment. The view factors between all segments are computed using the Hottel crossed string rule [32], and the reflectivity of the front of the panel is assumed to be 0 as in [19].
G T I = I F R   1 I r + D ,
I r = D H I × R b
D = D I F × R d
R b = m a x 0 , c o s θ c o s Φ ( 1 S )
R d = F s X i + X c s ( 1 S c s )
For the computation of the transposed direct normal (Equation (29)) and diffuse (Equation (30)) irradiance on each segment, the vectors R b and R d are obtained from Equations (31) and (32) which consist of the tilt factors for the direct (beam) and diffuse irradiances. The shading of the direct irradiance in each segment is taken into consideration in the direct tilt factor through the vector S , whose values are either 0, when direct irradiance is not obscured, or 1, when the direct irradiance is obscured, depending on the geometrical characteristics of the installation and the apparent position of the sun in the sky. Regarding the diffuse tilt factor, an isotropic component, X i , and a circumsolar component, X c s , are modeled as in a first-row transposition model, which, as mentioned above, can be any first-row model as long as it includes direct, circumsolar, and isotropic irradiance components. In the present work, a set of analytical transposition models for first rows are firstly assessed against experimental values and then one is selected, which will then be used in this model evaluation, as reported in the following sections. The sky view factors are computed for each segment considering the summation and reciprocity rules. The fraction of circumsolar irradiance that is obscured is modeled by the vector S c s ranging from 0, when there is no obscuration, to 1, when all irradiance from the considered circumsolar region is obscured.
The vectors S and S c s (Equations (33) and (34), respectively) are determined based on the model proposed by Varga and Mayer [18] with various modifications for the inclusion of cases in which the sun is positioned at the back of the row and the computations for the different segments of the back surface of the row in front of the one being evaluated and the ground. Depending on the surface of each segment, namely the surface of the panel being evaluated, c , the back surface of the panel of the row at its front, v , or the surface of the ground between them, u , the way the shading vectors are modeled differs. Since this is a two-dimensional model, firstly, the projection of the solar elevation angle to the azimuth of the panels, α , is needed, which is obtained through Equation (35), where γ s is the solar azimuth and γ p is the azimuth of the surfaces.
S = S c S u S v
S c s = S c s , c S c s , u S c s , v
α = tan 1 tan α cos γ s γ p
The shading of direct beam and the fraction of obscured circumsolar irradiance are determined for each segment of the panel, ground, and back of the front panel according to the projected solar elevation angle relative to each of the angles shown in Figure 2. These angles are computed using Equations (36)–(43), where i , j , and k are the indices of the segments in the panel being evaluated, at the back of the front panel, and on the ground between the rows of panels, respectively. The angles with subscript c are obtained for each segment of the panel being evaluated, those with subscript v are obtained for each segment of the back of the front panel, and those with subscript u are obtained for each segment of ground between the two rows. The angles δ and ε result from the view of the top of an adjacent panel from a segment, the angles ζ and ξ result from the view of the bottom of an adjacent panel from a ground segment, and the angles λ and σ result from the view of the top of a subsequent panel from a ground segment.
δ c ( i ) = tan 1 L c i + 1 c ( i ) 2 s i n β D + c i + 1 c ( i ) 2 L c o s β Π a = 1 ,
δ v ( k ) = tan 1 L v k + 1 v ( k ) 2 s i n β D
δ u ( j ) = tan 1 L s i n β + h 0 u j + 1 u ( j ) 2 L c o s β
ε u ( j ) = tan 1 L s i n β + h 0 D u j + 1 u j 2 + L c o s β
ζ u ( j ) = tan 1 h 0 u j + 1 u j 2
λ u ( j ) = tan 1 L s i n β + h 0 D + u j + 1 u j 2 L c o s β
ξ u ( j ) = tan 1 h 0 D u j + 1 u j 2
σ u ( j ) = tan 1 L s i n β + h 0 2 D u j + 1 u j 2 + L c o s β
The shading of direct irradiance for each segment of the panel surface, i , is given by Equation (44). Shading can occur ( S c i = 1 ) when the sun is positioned either in front of or behind the rows. Specifically, shading occurs when the projected solar elevation angle is positive and lower than the angle from the middle of the segment to the top of the front row (indicating that the sun is behind the front row), or when the projected solar elevation angle is negative and lower than the tilt angle of the surfaces (indicating that the sun is behind the panel). For the remaining cases, the segments are not shaded.
S c i = 1 ,   δ c i > α γ s γ p < 90 °   β > α γ s γ p 90 ° 0 ,   δ c i α γ s γ p < 90 °   β α γ s γ p 90 °
The shading of direct irradiance for each segment of the back surface of the front row, denoted by index k , is computed through Equation (45). Each segment is always shaded if the sun is positioned in front of the rows. Additionally, if the sun is not in front of the rows, shading occurs when the projected solar elevation angle exceeds the tilt angle of the surfaces or when it is lower than the angle from the middle of the segment to the top of the panel being evaluated.
S v k = 1 , γ s γ p < 90 ° γ s γ p 90 ° δ v k > α β < α   0 , γ s γ p 90 ° δ v k α β α  
The computation of the shading of direct irradiance on the ground between the rows is performed through Equation (46), which involves comparing the projected solar elevation angle with six other angles as depicted in Figure 2 and calculated through Equations (38) to (43).
Each ground segment, denoted by index j , is shaded under the following conditions: when the sun is positioned in front of the rows and behind the row in front, or when its projected solar elevation angle is lower than the angle of the middle of the segment to the top of the subsequent row in front. Additionally, ground segments are shaded if the sun is behind the row being evaluated or when the projected solar elevation angle is lower than the angle from the middle of the segment to the top of the subsequent row behind. Moreover, for ground segments positioned below the front row, shading occurs when the projected solar elevation angle exceeds the tilt of the surfaces.
S u j = 1 , γ s γ p < 90 ° ζ u j < α < δ u j δ u j > 0   λ u j > α δ u j 0 α > ζ u j γ s γ p 90 ° ε u j > ξ u j ε u j > α > ξ u j ε u j ξ u j ε u j < α < ξ u j δ u j < 0   δ u j < α α σ u j 0 , γ s γ p < 90 ° δ u j 0 δ u j α   λ u j α ζ u j γ s γ p 90 ° δ u j 0 ε u j α ξ u j α   σ u j α ξ u j α ε u j δ u j < 0   ε u j α ξ u j α δ u j α
The modeling of the circumsolar irradiance obscuring follows a similar approach, albeit with the utilization of the vector S c s to represent the fraction of circumsolar irradiance that is obscured. Circumsolar irradiance is assumed as being uniformly distributed within the annular region surrounding the sun disk, with an apparent external angular radius, r , of 15°. The obscured area of the circumsolar region is given by Equation (47) as in [18]:
C x = s i n 1 r 2 α x 2 r 180 r 2 α x 2 2 π r
In the case of the segments of the panel being evaluated (Equation (48)), several conditions dictate the obscuring of circumsolar irradiance. When the sun is positioned in front of the rows, total obscuration ( S c s , c i = 1 ) occurs if the entire circumsolar region is below the angle defined by the middle of the segment to the top of the front row, denoted as δ c i . Substantial obscuration (more than 50%, S c s , c i = 1 C ( δ c i ) ) occurs if this angle is higher than the projected solar elevation angle and the difference between these two angles is smaller than the angular radius. It is less obscured (less or equal to 50%, S c s , c i = C δ c i ) if the projected solar elevation angle is higher than δ c i and the difference between these two angles is smaller than the angular radius. Finally, no obscuration ( S c s , c i = 0 ) occurs if the projected solar elevation angle exceeds δ c i and the difference between these two angles is higher than the angular radius. A similar principle applies when the sun is behind the rows, with the difference that instead of using the angle δ c i , the comparison is made with the tilt angle β .
S c s , c i = 1 , γ s γ p < 90 ° α δ c k r γ s γ p 90 °   α β r 1 C δ c i , γ s γ p < 90 ° r < α δ c i < 0 1 C β , γ s γ p 90 °   r < α β < 0 C δ c i γ s γ p < 90 ° 0 α δ c i < r C β , γ s γ p 90 °   0 α β < r 0 , γ s γ p < 90 ° α δ c i r γ s γ p 90 °   α β r
The fraction of circumsolar irradiance obscured for each segment of the back surface of the front panels is given by Equation (49). It is assumed that all circumsolar irradiance is obscured when the sun is positioned in front of the rows. When the sun is positioned behind the rows, the projected solar elevation angles plus or minus the circumsolar angular radius are compared in a similar manner to Equation (48), but with reference to the tilt angle of the surface and the angle of the middle of the back surface segment to the top of the panel being evaluated.
S c s , v k = 1 ,     γ s γ p < 90 ° γ s γ p 90 °   α δ v k r α β r 1 C δ v k , γ s γ p 90 °   r < α δ v k < 0   1 C β , γ s γ p 90 °     0 < α β < r C δ v k γ s γ p 90 °   0 α δ v k < r C β , γ s γ p 90 °   r < α β 0   0 , γ s γ p 90 °   α δ v k r   α β r
The computation of circumsolar irradiance obscuration for each segment of the ground between rows is more complex (Equation (50)). Complete obscuration is assumed when the entire circumsolar region is positioned either behind the front row panel or the panel being evaluated, considering the middle of the ground segment. Additionally, complete obscuration occurs when the angle of the top of the circumsolar region is lower than the angle of the top of the subsequent rows of panels, whether in front or behind. Similarly to the equations above (Equations (48) and (49)), the projected solar elevation angles plus or minus the circumsolar angular radius are compared with the six different angles shown in the scheme at the bottom of Figure 2 for the computation of the fraction of circumsolar irradiance that is obscured.
S c s , u j = 1 ,     γ s γ p < 90 ° δ u j 0 α ζ u j r α δ u j r α λ u j r δ u j < 0 α ζ u j r α r λ u j α + r ζ u j γ s γ p 90 °   α ξ u j r α ε u j r α σ u j r δ u j 0 α δ u j r α r σ u j α + r ε u j α r σ u j α + r ξ u j 1 C ( δ u j ) , γ s γ p < 90 ° δ u j 0 r < α δ u j < 0 γ s γ p 90 °     δ u j < 0 0 α δ u j < r   1 C ( ε u j ) , γ s γ p 90 °   r < α ε u j < 0 1 C ( ζ u j ) , γ s γ p < 90 °   α r > λ u j 0 < α ζ u j < r 1 C ( λ u j ) , γ s γ p < 90 °   α + r < ζ u j r < α λ u j < 0   1 C ( ξ u j ) γ s γ p 90 °   α r > σ u j 0 < α ξ u j < r 1 C ( σ u j ) , γ s γ p 90 °     α + r < ξ u j r < α σ u j < 0 C δ u j γ s γ p < 90 ° δ u j 0 0 α δ u j < r γ s γ p 90 °     δ u j < 0 r α δ u j < 0 C ε u j , γ s γ p 90 °   0 α ε u j < r C ζ u j ,   γ s γ p < 90 °   α r > λ u j r < α ζ u j < 0 C λ u j ,   γ s γ p < 90 °   α + r < ζ u j 0 < α λ u j < r C ξ u j , γ s γ p 90 °   α r > σ u j r < α ξ u j < 0 C σ u j ,   γ s γ p 90 °   α + r < ξ u j 0 < α σ u j < r 0 ,     γ s γ p < 90 ° δ u j 0 α δ u j r α ζ u j r α λ u j r γ s γ p 90 °   α ε u j r α ξ u j r α σ u j r α ε u j r

4. Experimental Setup and Procedure

In order to obtain observational data for both a first row and subsequent rows of panels with varying tilt angles and inter-row distances, a structure featuring a pyranometer for measuring GTI was constructed in an open field near Évora, Portugal (38.5306°, −8.0112°), as shown in Figure 3, Figure 4 and Figure 5.
The experimental setup consists of three frames: a base, a front frame, and a rear frame, with a pyranometer installed on the rear frame. The base was leveled, and two transversal bars on the sides ensured that both front and rear frames maintained the same tilt angle. The apparatus allows for adjustment with three degrees of freedom: tilt angle of the front and rear frames ( β , from 20° to 90°) through the solidary adjustment of the inclination of both frames; distance between frames ( D , from 0.80 m to 1.10 m) through three positions of where the front frame can be fixed to the base; and position of the pyranometer along the length of the rear frame ( c s , from 0.08 m to 1.08 m) by sliding the instrument along its supporting bar. The uncertainty on the measurements of each of these variables is 1 cm in the cases of distances and 1° in the case of the tilt angles.
To represent the adjacent row, three Alveopan bilaminate white polypropylene boards, with a total width W of 3.03 m and length L of 1.08 m, were installed in the front frame. The reflectivity of these boards was measured using a FieldSpec HandHeld 2 spectroradiometer (ASD, Inc., Boulder, CO, USA) [33], yielding an average reflectivity of 0.921. Although potential edge effects were acknowledged due to board sizes, these were not factored into the general model. For data collection purposes relevant to the assessment of transposition models applied to the first rows, these boards were removed.
Global tilted irradiance was measured using a Kipp & Zonen CMP11 pyranometer (Kipp & Zonen, Delft, The Netherlands), while global horizontal irradiance and reflected irradiances (for the computation of ground albedo) were measured using a Kipp & Zonen CM7B albedometer (Kipp & Zonen, Delft, The Netherlands). Both sensors were connected to a CR300 datalogger from Campbell Scientific (Shepshed, Loughborough, UK). Additionally, DNI, DIF, and GHI observation data were obtained from the Évora–PECS station of the DNI-ALENTEJO project network [34], located 5 m from the experimental setup.
The internal clock of the CR300 data logger used in the apparatus was synchronized with the data logger of the Évora–PECS station, both set to UTC time. Sensor outputs were sampled at 1 Hz and mean, maximum, minimum, and standard deviation values were recorded every minute. Observations were corrected following the best practices in the field, namely the WMO recommendations and the BSRN (Baseline Solar Radiation Network) guide, including corrections for sensor zero offset and filters according to the BSRN quality control procedure, considering the extremely rare limits [35] and removing measurements for zenith angles equal or above 85°.
Prior to the field measurements, a calibration procedure was conducted specifically for the CMP11 pyranometer and CM7B albedometer using a reference CMP21 pyranometer (Kipp & Zonen, Delft, The Netherlands), according to the ISO 9847:1992 standard [36].
For the accurate application of transposition models, the ground albedo value is needed. While a standard value of 0.2 is often used, this study estimated the ground albedo using Equation (51) [37], where B S A and W S A represent the black-sky albedo and white-sky albedo, respectively. B S A is defined as the albedo in the absence of the diffuse component and is a function of the solar zenith angle, while the W S A is the albedo considering the diffuse component as isotropic and in the absence of the direct component. The mean values of B S A and W S A were obtained by fitting Equation (44) to experimental data, where the albedo, ρ , was computed by applying the ratio of reflected (GRI) to global horizontal (GHI) irradiance observations from the albedometer. Following data treatment, including filtering and removal of the records for solar zenith angles lower than 70°, B S A and W S A were found to be 0.206 and 0.208, respectively, across all recorded data periods.
ρ = B S A + ( W S A B S A ) D I F G H I
Observations were conducted between 14 April and 1 June 2022. For first-row tests, 5 datasets or periods were generated, each corresponding to a specific tilt angle, as shown in Table 1, with the pyranometer positioned at c s = 0.50 m. Testing of the developed model for other rows resulted in 19 periods with various tilt angles, distance between rows, and pyranometer positions, as shown in Table 2. Periods 4, 9, and 19 include instances in which the pyranometer is shaded.
It should be noted that measurements were conducted during a period of relatively high solar elevation, minimizing shading effects. To better capture shading effects, the inter-row distance during experimental tests were shorter than typical photovoltaic power plant configurations. Nonetheless, the developed model is designed to encompass diverse real-world conditions in the field.

5. Results and Discussion

The different transposition models, including the developed model for rows that are not the first, were applied to the observations of DNI and DIF from the Évora–PECS station. Subsequently, the model outputs were compared with GTI observations from the experimental setup using multiple evaluation metrics developed in the software MATLAB R2018b. These metrics comprised the coefficient of determination (R2), mean bias error (MBE), and root-mean-squared error (RMSE) along with a global performance index (GPI) based on R2, MBE, and RMSE, where a higher value represents the better accuracy of the model [22].
Given the significant influence of atmospheric conditions, particularly cloud cover, on global irradiance, the mean and standard deviation of the clearness index were computed for each period based on 1 min observations. The clearness index, typically ranging from 0 to 1, represents the ratio of global horizontal irradiance measured at ground level to its counterpart estimated at the top of the atmosphere [34] (Equation (15)). It serves as an indicator of the total transmittance of the atmosphere, reflecting higher values under clear-sky conditions and lower values under overcast conditions. The subsequent subsections provide detailed tables presenting these clearness index values for each period.
Some clearness index values exceeded unity, with a maximum value of 1.079, attributed to cloud enhancement events. These phenomena occur when partly cloudy skies lead to a temporary increase in local GHI above the extraterrestrial irradiance, facilitated by multiple scatterings and reflections by clouds [35]. These values were kept in the analysis, as they capture the transient nature of atmospheric conditions during the observation periods.

5.1. Results for the First Row

The transposition models presented in Section 2.1 were computed for the periods shown in Table 1, with the resulting mean and standard deviation of the clearness index and GTI, alongside various evaluation metrics, summarized in Table 3. Across all models, GTI values were generally underestimated, with the exception of the Bugler and Modified Bugler models, where GTI was overestimated. This discrepancy could stem from the treatment of the direct normal component within the factor R b used for modeling the diffuse component. Nevertheless, the Modified Bugler model showed better results compared with the other models. The Modified Bugler model tends to perform best except for vertical surfaces (period 5), where the Klucher model seems to show better results. Given that the global performance index (GPI) for the overall data presented the Modified Bugler model as the best performing model, it was chosen as the primary model for first-row applications in the proposed model. Despite its tendency to overestimate GTI, with an overall MBE of 23.8 W/m2 and RMSE of 30.5 W/m2, it delivered optimal results for period 3, characterized by a tilt angle of 50° and small variation in sky conditions (standard deviation of clearness index of 0.087).

5.2. Results for Other Rows

For the evaluation of the developed transposition model for rows other than the first, some adjustments were implemented to accommodate the experimental setup (refer to Figure 5). Given the absence of panels in the second row, this row was not considered in the modeling. Instead, G T I was computed for each segment of the back of the front panel and ground, followed by the calculation of G T I at a designated point which represents the sensor. In this case, the angles ε u , λ u , ξ u , σ u , and δ v used for shadow computation and circumsolar irradiance obscuration were not applicable. Another modification involved the length of ground considered in the model. Since the setup comprised only one panel, reflections from the ground beyond the modeled rows could significantly impact the measured GTI and were thus incorporated into the model validation process (depicted in Figure 5).
GTI estimation was performed using both the developed model and the Modified Bugler model for reference and comparison, which is a common practice in the absence of a specific model for other rows. It is important to note that the configurations used in periods 4, 9, and 19 result in direct shading of the pyranometer for a certain time span of the day. As example, periods 1 and 12, when there was no shading, and period 19, when there was shading and obscuration of the pyranometer by the front row, are shown in Figure 6, Figure 7 and Figure 8 (the small data gaps during the day are a result of the filtering procedure mentioned in Section 4). Despite the fact that slightly lesser improvements are observed for period 12 (Figure 6), which is attributed to partially cloudy conditions, the effectiveness of the developed model over the Modified Bugler is evident across the evaluated periods.
The results for each period, all periods, and for the data when the pyranometer is shaded or unshaded are presented in Table 4. When compared with the original Modified Bugler model, which overestimates the GTI, the proposed model improves the MBE for most periods, albeit with a slight underestimation. Typically, the Modified Bugler performs better in periods characterized by higher pyranometer positioning and greater frame-to-frame distances, resembling first-row irradiance conditions. During periods of direct irradiance shading (periods 4, 9, and 19), the Modified Bugler model, which does not consider shading, shows significantly higher errors. Another aspect to highlight is the impact of clouds in the performance of the model. In periods 9 and 14, for example, when the mean clearness index is lower and its standard deviation is higher (indicating cloudier skies) the metrics show lower performance of both models.
Due to the variation in sky conditions along the different periods, the comparison between the different positioning of the setup proved challenging and thus, more importance is given to the overall results instead of each period. In this regard, the developed model for rows affected by the presence of rows in front showed an MBE of −12.9 W/m2 and a root-mean-squared error of 76.8 W/m2. As expected, this model outperforms the first-row model when the pyranometer is shaded. Even under unshaded conditions, the developed model is better than the model for the first row, showing the impact of the obscuring of the sky dome due to the other rows, considering the sky radiance anisotropy, namely the circumsolar region, and of the reflections from the front row and ground on the GTI.
To quantify the impact of the proposed model on reducing the error for each irradiance component on a tilted surface compared with the Modified Bugler model, a weight ( w i ) for each component i of GTI was computed through Equation (52):
w i = 1 n i D i M B G T I D G T I M B  
Here, D stands for developed model and M B for the Modified Bugler model. The results are presented in Table 5 for the direct, I b , diffuse circumsolar, D c s , diffuse isotropic, D i s o , and reflected, I r e f l , components and for both unshaded and shaded conditions. For the proposed model, the reflected component includes reflections in the ground and in the back side of the front panel. The mean bias error of each model is also included in the table for reference.
Overall, the masking of the isotropic diffuse irradiance has the highest weight in the difference between the Modified Bugler and the developed model followed by the modeling of reflected irradiance. As expected, when the direct irradiance is shaded, this becomes the most impactful component, while for unshaded conditions, it has no impact.

6. Operational Algorithm for Forecasting of Solar Irradiance on Tilted Surfaces

This section presents the application of the developed model to predict irradiance on tilted surfaces using operational direct normal and diffuse irradiance forecasts as inputs.

6.1. Forecast Input Data

The transposition models presented in this work, in particular the developed model for rows that are not the first, require several geometrical details on the solar panels of the power plant as input. However, only two meteorological variables, specifically DNI and DIF, are necessary for the determination of the GTI. The European Centre for Medium-range Weather Forecasts (ECMWF) developed the Integrated Forecasting System (IFS), which issues global hourly forecasts every day at 00 UTC up to 10 days ahead. These forecasts include the GHI and DNI variables, while DIF can be computed through the closure function G H I = B H I + D I F (where B H I is given by Equation (2)).
In the work by Pereira et al. [22], a method based on artificial neural networks (ANNs) was developed to generate improved spatial and temporal downscaled DNI forecasts using operational forecast outputs from the ECMWF/IFS and the Copernicus Atmospheric Monitoring Service (CAMS). This work employs the same models and procedure, while also including the analysis and usage of GHI as compared with previous work. Figure 9, Figure 10 and Figure 11 provide a direct comparison between ECMWF/IFS forecasts and observations made through calibrated sensors at Évora–Verney station in Portugal (38.5678, −7.9114). These figures show data for DNI, GHI, and DIF for each forecast day of the forecast horizon (day 0: 1 h to 24 h; day 1: 25 h to 48 h; day 2: 49 h to 72 h) over the period from 1 December 2016 to 31 May 2021, with a timestep of 10 min.
As expected, shorter forecast horizons correspond to smaller errors in the forecasts. Improved DNI forecasts were obtained using artificial neural network models, accounting for the nonlinear relationships between the DNI and various meteorological forecasted variables, such as GHI, cloud cover, and aerosol optical depth, as well as the variation in DNI over a given period before the forecasted instant. A brief explanation of this model and procedure can be found in Appendix A. The application of this procedure results in the metrics shown in Figure 12, where the diffuse component was computed using the GHI forecasted by the ECMWF/IFS and the improved DNI forecasts generated by the ANN models.
The application of the ANN model shows improvements in the forecasting of both DNI and DIF over a dataset that encompasses different sky conditions. The same procedure was applied by retrieving weather forecasts for the measurement periods used in this work to evaluate the different transposition models (Section 4). The resulting 10 min improved DNI and computed DIF forecasts were then evaluated against the 10 min mean observations obtained from the experimental data presented in Section 4.
Table 6 shows the results for first-row and inner-row tests. Detailed results for each testing period can be found in Appendix B. Similar to the larger dataset, the longer the forecast horizon, the higher the errors. However, the errors obtained tend to be greater than the ones aforementioned. This discrepancy is due to the smaller dataset used for this analysis, where even minor differences between forecasts and observations can lead to significant errors. Forecasts of DNI and DIF tend to have larger errors compared to global variables such as GHI or GTI, since accurate forecasting of clouds, specifically their location, is critical for this temporal resolution. The GTI results from the transposition models, based on the output of these forecast models, are presented in the following.

6.2. Operational Analysis of Transposition Models

To understand how the selected transposition models perform in an operational forecast context, the transposition models were applied using the resulting forecasts of DNI and DIF for the testing location with a 10 min timestep. The results were then compared with the observations.
Table 7 presents the results for the first row using the Modified Bugler model and for inner rows using the developed model for each day of the forecast horizon.
Regarding first-row results, when compared with the results using DNI and DIF measurements instead of forecasts (Table 3), it is evident that for the analyzed periods, the use of forecasts shows better results for days 0 and 1 regarding MBE. However, there is a deterioration for day 2 and across all forecast horizons regarding the MAE and RMSE.
As for the inner-row test results, when compared with GTI based on DNI and DIF measurements instead of forecasts (Table 4), there is a general deterioration in the results, except for the MBE in unshaded conditions for forecast days 0 and 1. Detailed results for each testing period can be found in Appendix B. The use of the developed transposition model resulted in overall MBE and RMSE values of 33.6 W/m2 and 169.7 W/m2, respectively, for the entire forecast horizon. These values show that the model is eneficial for linking irradiance forecast models with energy generation modeling in solar power plants, aiding in the production of power output forecasts. These forecasts are essential for better decision-making by operators of such energy systems due to the variability of resource and energy demand.
To understand how the mean bias error of DNI and DIF forecasts affects the bias of the computed GTI values, Figure 13 was created. This figure shows the difference between the MBE of GTI from the developed model using forecasts and experimental data as input as a function of the MBE of DNI and DIF predictions for the forecast at day 0. Positive values of this difference are represented with blue circles, while negative values are represented with white circles. For better readability, only periods with an MBE difference of less than 100 W/m2 are represented, which resulted in excluding only period 1.
Since the MBE of the proposed model is typically negative for the tests carried out (Table 7 and Table A5), and considering the relationship between the MBE values for the forecasted DNI and DIF components (when one tends to higher positive values, the other tends to more negative values (Table 6 and Table A3) as a consequence of the closure equation and knowing that the MBE of the GHI forecasts is lower), the bias of the GTI using prediction values may decrease as shown in Figure 13 due to a favorable combination of forecast and model errors. This is the reason why the MBE for unshaded conditions using forecast values (Table 7) is lower than the bias error for the same conditions using experimental values (Table 4). However, this trend is not observed for the shaded conditions, where the mean bias error is still lower when using ground-based measurements. A more detailed analysis of these aspects is needed in the future, which extends beyond the scope of the present work. Also, due to spatial resolution limitations, the forecasted variables from ECMWF (current operational horizontal resolution: ~9 km) are identical for the first and second rows in this experiment. In future work, this algorithm should be validated for a larger power plant that encompasses multiple grid points of a high-resolution forecast model at a hectometric scale.

7. Conclusions

This work presented a comprehensive analysis of nine analytical transposition models based on physics alongside the Perez transposition model to compute the global tilted irradiance on photovoltaic module surfaces. Additionally, it presented a model for the computation of this variable in rows of modules other than the first, which usually comprises most rows of solar power plants. The developed model can be applied to any first-row transposition model, provided it considers direct, circumsolar, and isotropic diffuse irradiance. This model computes the GTI for different longitudinal segments of the surfaces of the row of modules, the back of the row in front, and the ground between the rows. It takes into consideration the different view factors and the obscuring of direct and circumsolar irradiance for each of the segments for any apparent solar position and includes the shading effect of the succeeding rows on the ground segments.
The evaluation of these models utilized data collected in Évora, Portugal, for different tilt angles for first-row tests and also for different inter-row distances, including shading conditions, to assess the performance of the developed model. The clearness index helps address potential confounding variables by providing a baseline for sunshine conditions. However, we acknowledge that factors such as wind, relative humidity, precipitation, and aerosols also affect the experiment to some extent. In future research, these factors should be considered. Also, our dataset is limited to our experimental setup. While expanding to diverse locales would enhance generalizability, the focus of this study was on addressing challenges specific to inner rows of solar panels, rather than aiming for global applicability. Thus, while our findings offer valuable insights, they may not directly apply to all regions.
Results showed that the best analytical transposition model for the first row is the Modified Bugler model, showing an overall MBE of 23.8 W/m2 and RMSE of 30.5 W/m2. Conversely, for other rows, the developed model showed an MBE of −12.9 W/m2 and RMSE of 76.8 W/m2, resulting in an improvement of 368.3 W/m2 and 224.4 W/m2, respectively, compared to using the selected reference transposition model for first rows. This shows the importance of considering the direct shading and obscuring of the sky dome when computing GTI for surfaces in rows that are not the first.
Furthermore, the operational performance of transposition models was evaluated for GTI forecasting, using improved irradiance forecast values instead of measurements of DNI and DIF. These forecast values were obtained from artificial neural network models using numerical weather prediction and aerosol forecast data. Results of the first-row tests showed a MBE and RMSE for all data of 6.7 W/m2 and 48.7 W/m2 for forecast day 0, 15.5 W/m2 and 55.0 W/m2 for forecast day 1, and 71.5 W/m2 and 128.5 W/m2 for forecast day 2. This shows an increased error compared to results using observations which are a mean MBE and RMSE increase across the three days of forecast of 7.4 W/m2 and 46.9 W/m2. It also shows how the forecast performance tends to deteriorate with time. The same is visible for the tests performed for other rows, which show an overall MBE and RMSE of −31.5 W/m2 and 164.0 W/m2 for forecast day 0, −32.1 W/m2 and 165.6 W/m2 for forecast day 1, and −37.2 W/m2 and 179.4 W/m2 for forecast day 2.
This work demonstrated that transposition models that neglect shading and irradiance obscuration are not suitable for the accurate estimation of GTI in surfaces that are not in the front row of a solar power plant. The use of a dedicated model for these conditions, such as the one presented in this work, is of great importance, given that GTI is the main factor influencing the energy generation of solar photovoltaic systems.

Author Contributions

Conceptualization, S.P. and P.C.; methodology, S.P. and P.C.; validation, S.P. and P.C.; formal analysis, S.P. and P.C.; investigation, S.P., P.C. and R.S.; data curation, S.P. and P.C.; writing—original draft preparation, S.P.; writing—review and editing, S.P., P.C. and R.S.; visualization, S.P. and P.C.; supervision, P.C. and R.S. All authors have read and agreed to the published version of the manuscript.

Funding

This work was funded by National funds through FCT—Fundação para a Ciência e Tecnologia, I.P. (projects UIDB/04683/2020 and UIDP/04683/2020). S. Pereira acknowledges the support of FCT (the Portuguese Science and Technology Foundation) through the grant with reference SFRH/BD/145378/2019.

Data Availability Statement

The raw data supporting the conclusions of this article will be made available by the authors on request. Restrictions apply to the availability of these data. Data were obtained from ECMWF.

Acknowledgments

The authors are thankful to Afonso Cavaco, Josué Figueira, and Samuel Ramos Bárias for their technical support.

Conflicts of Interest

The authors declare no conflicts of interest.

Glossary

Nomenclature
A Anisotropy index
BHIBeam horizontal irradiance (W/m2)
C Fraction of circumsolar area obscured by the adjacent row
c Length along the panel (m)
c s Position of the pyranometer along the length of the panel (m)
D Distance between rows in the horizontal plane (m)
D Vector of diffuse horizontal irradiance vector (W/m2)
D c s Circumsolar diffuse component of GTI (W/m2)
DIFDiffuse horizontal irradiance (W/m2)
D i s o Isotropic diffuse component of GTI (W/m2)
DNIDirect normal irradiance (W/m2)
EHIExtraterrestrial horizontal irradiance (W/m2)
ENIExtraterrestrial normal irradiance (W/m2)
FView factor matrix
F g Ground view factor
F s Sky view factor
f Clearness index as defined in the Klucher model
GTIGlobal tilted irradiance (W/m2)
GRIGlobal reflected irradiance (W/m2)
h 0 Vertical distance between the ground and the panel base (m)
IIdentity matrix
I b Direct component of GTI (W/m2)
I r Vector of direct horizontal irradiance vector (W/m2)
I r e f l Reflection component of GTI (W/m2)
k T Clearness index
L Length of the panel (m)
MAEMean absolute error (W/m2)
MBEMean bias error (W/m2)
MAPEMean absolute percentage error (%)
RReflectivity matrix
R b Beam irradiance tilt factor
R d Diffuse irradiance tilt factor
RMSERoot-mean-squared error (W/m2)
RMSPERoot-mean-squared percentage error (%)
R2Coefficient of determination
R r Reflected irradiance tilt factor
r Apparent angular radius of circumsolar irradiance (°)
S Shading of direct component coefficient
S c s Shading of circumsolar component coefficient
u Length along the ground (m)
v Length along the back of the front panel (m)
W Row width (m)
w i Weight of irradiance component i to the bias reduction (%)
X i Isotropic component of first-row model
X c s Circumsolar component of first-row model
Acronyms
ANNArtificial neural networks
CAMSCopernicus Atmospheric Monitoring Service
ECMWFEuropean Centre for Medium-range Weather Forecasts
IFSIntegrated forecasting system
NWPNumerical weather prediction
Greek symbols
ΦSolar zenith angle (°)
αSolar elevation angle (°)
αProjection of α in the vertical plane of the local meridian (°)
βSurface tilt angle (°)
γpSurface azimuth (°)
γsSolar azimuth (°)
δAngle between the horizontal and the top of the front panel (°)
εAngle between the horizontal and the top of the panel being considered (°)
ζAngle between the horizontal and the bottom of the front panel (°)
θAngle of incidence (°)
λAngle between the horizontal and the bottom of the panel being considered (°)
ξAngle between the horizontal and the top of the panel in front of the front panel (°)
ρGround albedo
σAngle between the horizontal and the top of the panel behind the panel being considered (°)

Appendix A

The flowchart of the model used to generate forecast data of direct normal irradiance (DNI) [22] used in this work is shown in Figure A1.
Figure A1. Flowchart of the model used to generate DNI forecasts [22].
Figure A1. Flowchart of the model used to generate DNI forecasts [22].
Energies 17 03444 g0a1
The model takes as input data the variables shown in Table A1 from the operational Integrated Forecasting System (IFS) of the European Centre for Medium-range Weather Forecasts (ECMWF) and the Copernicus Atmospheric Monitoring Service (CAMS). These are run every day at 00 UTC, providing hourly forecast values up to 90 h ahead at discrete points of a global grid with a horizontal spatial resolution of 0.125° × 0.125°.
Table A1. Input variables obtained from numerical weather prediction models.
Table A1. Input variables obtained from numerical weather prediction models.
Variables Obtained from IFS/ECMWFVariables Obtained from CAMS
Direct normal irradianceTotal aerosol optical depth at 670 nm
Global horizontal irradianceTotal aerosol optical depth at 865 nm
Low cloud coverTotal aerosol optical depth at 1240 nm
Medium cloud coverSea salt aerosol optical depth at 550 nm
High cloud cover
Total cloud cover
Wind speed
Air temperature
Solar zenith angle
A temporal and spatial downscaling is performed on these variables, which results in forecasts for a specific location and higher temporal resolution. This downscaling is obtained through a bi-linear interpolation of the values in the four surrounding grid points of the desired location for the spatial downscaling and a piecewise cubic hermite interpolation of the hourly values into smaller timesteps (in this work, it was 10 min values) for the temporal downscaling.
The downscaled variables are then fed to an artificial neural network (ANN model A), which is a feed-forward network with one hidden layer comprising seven neurons and uses (i) a backpropagation learning function, namely Bayesian regularization backpropagation; (ii) a linear layer output with an initialization function that initializes the weights and biases of the layers according to the Nguyen–Widrow initialization algorithm; (iii) the hyperbolic tangent sigmoid transfer function; and (iv) the mean-squared error as a performance function. The input and output data are processed by removing rows with constant values and scaling the mean of each row to 0 and deviations to 1. This ANN takes into consideration the nonlinear relationships between the different atmospheric and aerosol variables and DNI resulting in improved DNI forecasts for the specified location and temporal resolution.
A second artificial neural network (ANN model B) is similar to ANN model A, but uses the Levenberg–Marquardt backpropagation algorithm and has eight neurons in the hidden layer, taking as input a time series of 12 timesteps of the improved DNI forecasts prior to the forecast moment (from the output of the ANN model A) along with the season and time of day. This model takes into consideration the temporal variation in DNI and further improves the DNI forecasts from ANN model A.

Appendix B

Table A2. Metrics of DNI and DIF forecasts for the measuring periods used for first-row tests (MBE, MAE, and RMSE in W/m2).
Table A2. Metrics of DNI and DIF forecasts for the measuring periods used for first-row tests (MBE, MAE, and RMSE in W/m2).
Period and
Variable
Day 0Day 1Day 2
R2MBEMAERMSER2MBEMAERMSER2MBEMAERMSE
1DNI0.5114−33.157.087.60.6120−14.938.367.90.4664−15.372.299.5
DIF0.122818.526.539.10.28746.011.013.70.0912−2.832.037.9
2DNI0.9082−68.768.972.70.7794−30.439.558.60.6228−36.161.685.1
DIF0.811940.841.146.50.00636.815.921.10.00197.031.736.5
3DNI0.3542−50.992.0115.00.0803−17.7175.5204.50.0014−13.2131.0164.0
DIF0.6937−8.759.770.40.52900.678.986.80.6843−9.470.477.7
4DNI0.247047.474.0127.70.2273−33.279.1128.70.1022−322.1322.1391.3
DIF0.2873−57.558.494.00.4425−24.540.269.20.488411.743.269.3
5DNI0.4456−86.6149.2220.40.1410−121.6206.9283.70.2980−132.1180.7256.1
DIF0.5380−3.232.539.90.051333.557.970.80.019536.965.176.6
Table A3. DNI and DIF forecast results for the measuring periods used for testing of the developed model (MBE, MAE, and RMSE in W/m2).
Table A3. DNI and DIF forecast results for the measuring periods used for testing of the developed model (MBE, MAE, and RMSE in W/m2).
Period and
Variable
Day 0Day 1Day 2
R2MBEMAERMSER2MBEMAERMSER2MBEMAERMSE
1DNI0.195726.3101.8141.40.112714.192.3146.10.0151−28.4128.6162.6
DIF0.1687−59.264.276.00.4663−36.436.952.90.00064.759.567.6
2DNI0.955255.455.666.80.476828.559.880.90.91539.421.733.9
DIF0.5153−48.548.552.30.0240−20.431.340.10.2740−23.830.739.8
3DNI0.675324.0121.2190.30.6811−16.6134.1187.70.69302.9134.2184.3
DIF0.6857−19.448.772.40.6798−15.551.371.10.5439−15.460.785.3
4DNI0.087946.2178.6273.20.193459.5162.1253.60.227751.8162.0251.3
DIF0.2560−55.363.197.30.5052−60.965.290.10.4056−55.664.691.0
5DNI0.1872−58.2106.8140.40.0197−45.0133.4184.80.00650.483.2149.6
DIF0.199921.240.346.00.0189−11.838.558.50.0143−4.135.656.4
6DNI0.3289−147.5169.6210.30.367346.1847.289.70.301338.046.689.3
DIF0.015288.096.9129.00.0550−33.533.546.30.0019−25.428.544.2
7DNI0.8260−24.436.744.50.898111.931.037.30.839437.649.065.9
DIF0.620514.622.628.40.4374−3.923.926.60.0125−25.938.149.2
8DNI0.867225.853.560.20.816930.574.985.30.865657.561.869.2
DIF0.5586−3.211.815.80.3760−13.726.830.60.0822−26.526.530.9
9DNI0.1415160.6213.8270.00.163599.9184.8223.00.0117190.2290.4341.6
DIF0.3964−102.6115.7143.90.0298−96.5126.6157.60.0196−102.5131.1161.4
10DNI0.339693.7196.6266.00.266133.3204.1268.50.4108−20.3201.3235.9
DIF0.2564−42.390.2122.90.2894−20.085.2113.80.3719−17.387.1105.8
11DNI0.2717−129.9170.9207.20.6058−42.883.0111.20.625722.069.2103.2
DIF0.083255.770.896.50.092922.937.445.30.3765−7.923.030.9
12DNI0.042045.0197.2228.60.106995.6258.8310.00.3188197.5295.7363.3
DIF0.5633−63.897.1113.20.6490−114.3128.1152.10.7377−108.5121.3150.0
13DNI0.2652200.3217.9265.30.2650231.0252.1300.20.2762262.2299.9341.5
DIF0.1240−77.5104.1127.20.1153−99.3119.9141.80.1656−119.1135.3154.5
14DNI0.3747101.8137.5163.50.4259156.6162.3210.60.2201107.0145.8186.7
DIF0.455526.289.6111.40.436914.885.0112.20.431455.295.9122.8
15DNI0.1945−102.7125.5161.50.2203−148.3173.5190.00.2292−116.1139.4170.9
DIF0.387266.868.793.00.440081.781.797.40.453373.076.7102.6
16DNI0.080520.7230.0265.30.1606−34.3209.9246.70.0197−197.8326.9405.4
DIF0.567518.953.065.60.479231.264.277.20.7111−15.646.052.8
17DNI0.3496−15.592.7146.00.4274−37.6119.2158.30.4238−96.0151.4205.5
DIF0.147914.329.438.50.274511.231.139.70.378231.241.150.6
18DNI0.02625.2145.8211.60.098319.4139.8198.70.015439.0142.2210.7
DIF0.09949.876.4102.30.19213.772.695.80.0518−4.081.5111.6
19DNI0.7739−27.365.699.70.89400.142.267.50.7421−33.489.0121.4
DIF0.43848.120.827.40.8150−2.112.918.20.065434.845.460.6
Table A4. GTI metrics of first-row tests for each day of forecast horizon and for each period using the Modified Bugler transposition model with DNI and DIF forecasts (MBE, MAE, and RMSE in W/m2).
Table A4. GTI metrics of first-row tests for each day of forecast horizon and for each period using the Modified Bugler transposition model with DNI and DIF forecasts (MBE, MAE, and RMSE in W/m2).
PeriodDay 0Day 1Day 2
R2MBEMAERMSER2MBEMAERMSER2MBEMAERMSE
10.979529.233.546.40.979829.535.646.40.971923.041.448.7
20.99808.913.816.60.994112.722.625.60.993311.023.827.7
30.9224−29.445.463.90.951131.645.765.90.934516.336.456.3
40.91779.844.759.40.871519.851.874.60.0893−203.9210.1277.7
50.963321.634.940.60.932921.841.748.40.974247.029.934.6
Table A5. GTI results of the developed model for rows that are not the front row for each day of forecast horizon and for each period using DNI and DIF forecasts (MAE and RMSE in W/m2, MAPE and RMSPE in %).
Table A5. GTI results of the developed model for rows that are not the front row for each day of forecast horizon and for each period using DNI and DIF forecasts (MAE and RMSE in W/m2, MAPE and RMSPE in %).
PeriodDay 0Day 1Day 2
R2MBEMAERMSER2MBEMAERMSER2MBEMAERMSE
10.0043−163.0207.1364.90.0059−163.5195.8364.00.0048−179.4215.5365.9
20.9913−30.734.938.80.9835−41.042.757.20.9768−50.353.270.5
30.6449−17.4117.5175.00.6302−36.0129.9181.50.5980−25.8128.2187.2
40.3640−3.781.8137.80.3764−3.084.0138.30.36301.185.3142.8
50.4550−87.7108.9214.50.4017−99.6131.8234.10.4038−72.3105.8217.9
60.9196−105.8113.3141.00.9651−29.245.257.80.9636−30.846.859.7
70.9938−43.343.848.80.9981−31.131.635.00.9971−26.527.030.9
80.9958−48.050.052.50.9963−50.251.553.80.9949−49.751.754.6
90.1442−82.6112.5166.30.1479−93.7114.7166.20.0815−84.5127.3179.2
100.1607−62.4142.8242.50.1907−78.1143.9239.90.2246−105.8147.0242.8
110.8884−92.897.4125.50.9728−62.167.478.00.9823−37.544.151.3
120.5099−43.6108.8150.60.4232−36.7136.7181.20.548844.0111.4167.7
130.911119.957.171.50.915428.654.771.80.902541.864.987.6
140.114295.4105.1135.20.1097107.8113.6155.50.0710111.3121.3155.7
150.8651−61.974.891.60.8613−68.587.496.90.8723−64.578.091.8
160.47645.2166.3204.40.4233−30.8172.6212.80.0174−227.0295.8380.3
170.9180−43.970.885.80.9211−59.276.791.80.9147−79.190.5107.8
180.8772−41.781.499.70.8803−39.978.197.90.8741−41.983.2100.8
190.6672−55.365.9114.30.6605−51.862.1114.60.6384−39.451.7109.7

References

  1. SolarPower Europe. Global Market Outlook for Solar Power 2023–2027; SolarPower Europe: Brussels, Belgium, 2023; Available online: https://www.solarpowereurope.org/insights/market-outlooks/global-market-outlook-for-solar-power-2023-2027-1 (accessed on 5 June 2024).
  2. International Energy Agency Renewables 2023; Analysis and Forecast to 2028; International Energy Agency: Paris, France, 2023.
  3. Nassar, Y.F.; Hafez, A.A.; Alsadi, S.Y. Multi-Factorial Comparison for 24 Distinct Transposition Models for Inclined Surface Solar Irradiance Computation in the State of Palestine: A Case Study. Front. Energy Res. 2020, 7, 1–19. [Google Scholar] [CrossRef]
  4. Yang, D. Solar Radiation on Inclined Surfaces: Corrections and Benchmarks. Sol. Energy 2016, 136, 288–302. [Google Scholar] [CrossRef]
  5. Gueymard, C.A. Direct and Indirect Uncertainties in the Prediction of Tilted Irradiance for Solar Engineering Applications. Sol. Energy 2009, 83, 432–444. [Google Scholar] [CrossRef]
  6. Gueymard, C. An Anisotropic Solar Irradiance Model for Tilted Surfaces and Its Comparison with Selected Engineering Algorithms. Sol. Energy 1987, 38, 367–386. [Google Scholar] [CrossRef]
  7. Klucher, T.M. Evaluation of Models to Predict Insolation on Tilted Surfaces. Sol. Energy 1979, 23, 111–114. [Google Scholar] [CrossRef]
  8. Hay, J.E. Calculating Solar Radiation for Inclined Surfaces: Practical Approaches. Renew. Energy 1993, 3, 373–380. [Google Scholar] [CrossRef]
  9. Hay, J.; McKay, D. Final Report IEA Task IX-Calculation of Solar Irradiances for Inclined Surfaces: Verification of Models Which Use Hourly and Daily Data; International Energy Agency: Paris, France, 1988. [Google Scholar]
  10. Nassar, Y.F.; Abuhamoud, N.M.; Miskeen, G.M.; El-Khozondar, H.J.; Alsadi, S.Y.; Ahwidi, O.M. Investigating the Applicability of Horizontal to Tilted Sky-Diffuse Solar Irradiation Transposition Models for Key Libyan Cities. In Proceedings of the 2022 IEEE 2nd International Maghreb Meeting of the Conference on Sciences and Techniques of Automatic Control and Computer Engineering (MI-STA 2022), Sabratha, Libya, 23–25 May 2022; pp. 9–14. [Google Scholar] [CrossRef]
  11. Mayer, M.J.; Gróf, G. Extensive Comparison of Physical Models for Photovoltaic Power Forecasting. Appl. Energy 2021, 283, 116239. [Google Scholar] [CrossRef]
  12. Khan, M.M.; Ahmad, M.J. Estimation of Global Solar Radiation Using Clear Sky Radiation in Yemen. J. Eng. Sci. Technol. Rev. 2012, 5, 12–19. [Google Scholar] [CrossRef]
  13. Padovan, A.; Del Col, D. Measurement and Modeling of Solar Irradiance Components on Horizontal and Tilted Planes. Sol. Energy 2010, 84, 2068–2084. [Google Scholar] [CrossRef]
  14. Muneer, T.; Gueymard, C.; Kambezidis, H. Hourly Slope Irradiation and Illuminance. In Solar Radiation and Daylight Models; Butterworth-Heinemann: Oxford, UK, 2004; p. 143. [Google Scholar] [CrossRef]
  15. Badescu, V. 3D Isotropic Approximation for Solar Diffuse Irradiance on Tilted Surfaces. Renew. Energy 2002, 26, 221–233. [Google Scholar] [CrossRef]
  16. Tian, Y.Q.; Davies-Colley, R.J.; Gong, P.; Thorrold, B.W. Estimating Solar Radiation on Slopes of Arbitrary Aspect. Agric. Meteorol. 2001, 109, 67–74. [Google Scholar] [CrossRef]
  17. Appelbaum, J.; Massalha, Y.; Aronescu, A. Corrections to Anisotropic Diffuse Radiation Model. Sol. Energy 2019, 193, 523–528. [Google Scholar] [CrossRef]
  18. Varga, N.; Mayer, M.J. Model-Based Analysis of Shading Losses in Ground-Mounted Photovoltaic Power Plants. Sol. Energy 2021, 216, 428–438. [Google Scholar] [CrossRef]
  19. Tschopp, D.; Jensen, A.R.; Dragsted, J.; Ohnewein, P.; Furbo, S. Measurement and Modeling of Diffuse Irradiance Masking on Tilted Planes for Solar Engineering Applications. Sol. Energy 2022, 231, 365–378. [Google Scholar] [CrossRef]
  20. ECMWF IFS Documentation. Available online: https://www.ecmwf.int/en/publications/ifs-documentation (accessed on 10 December 2021).
  21. Pereira, S.; Abreu, E.F.M.; Iakunin, M.; Cavaco, A.; Salgado, R.; Canhoto, P. Method for Solar Resource Assessment Using Numerical Weather Prediction and Artificial Neural Network Models Based on Typical Meteorological Data: Application to the South of Portugal. Sol. Energy 2022, 236, 225–238. [Google Scholar] [CrossRef]
  22. Pereira, S.; Canhoto, P.; Salgado, R. Development and Assessment of Artificial Neural Network Models for Direct Normal Solar Irradiance Forecasting Using Operational Numerical Weather Prediction Data. Energy AI 2023, 15, 100314. [Google Scholar] [CrossRef]
  23. Liu, B.Y.H.; Jordan, R.C. The Long-Term Average Performance of Flat-Plate Solar-Energy Collectors. Sol. Energy 1963, 7, 53–74. [Google Scholar] [CrossRef]
  24. Bugler, J.W. The Determination of Hourly Insolation on an Inclined Plane Using a Diffuse Irradiance Model Based on Hourly Measured Global Horizontal Insolation. Sol. Energy 1977, 19, 477–491. [Google Scholar] [CrossRef]
  25. Temps, R.C.; Coulson, K.L. Solar Radiation Incident upon Slopes of Different Orientations. Sol. Energy 1977, 19, 179–184. [Google Scholar] [CrossRef]
  26. Iqbal, M. An Introduction to Solar Radiation; Elsevier: Vancouver, BC, Canada, 1983; ISBN 9780123737502. [Google Scholar]
  27. Kasten, F. A New Table and Approximation Formula for the Relative Optial Air Mass. Arch. Für. Meteorol. Geophys. Und. Bioklimatol. Ser. B 1965, 14, 206–223. [Google Scholar] [CrossRef]
  28. Reindl, D.T.; Beckman, W.A.; Duffie, J.A. Evaluation of Hourly Tilted Surface Radiation Models. Sol. Energy 1990, 45, 9–17. [Google Scholar] [CrossRef]
  29. Perez, R.; Ineichen, P.; Seals, R.; Michalsky, J.; Stewart, R. Modeling Daylight Availability and Irradiance Components from Direct and Global Irradiance. Sol. Energy 1990, 44, 271–289. [Google Scholar] [CrossRef]
  30. Saint-Drenan, Y.-M.; Barbier, T. Data-Analysis and Modelling of the Effect of Inter-Row Shading on the Power Production of Photovoltaic Plants. Sol. Energy 2019, 184, 127–147. [Google Scholar] [CrossRef]
  31. Maor, T.; Appelbaum, J. View Factors of Photovoltaic Collector Systems. Sol. Energy 2012, 86, 1701–1708. [Google Scholar] [CrossRef]
  32. Eckert, E.R.G. Radiative Transfer, H. C. Hottel and A. F. Sarofim, McGraw-Hill Book Company, New York, 1967. 52 Pages. AIChE J. 1969, 15, 794–796. [Google Scholar] [CrossRef]
  33. ASD FieldSpec® HandHeld 2 User Manual 2010, 1–140. Available online: https://www.geo-informatie.nl/courses/grs60312/material2017/manuals/600860-dHH2Manual.pdf (accessed on 5 June 2024).
  34. Cavaco, A.; Canhoto, P.; Collares Pereira, M. Procedures for Solar Radiation Data Gathering and Processing and Their Application to DNI Assessment in Southern Portugal. Renew. Energy 2021, 163, 2208–2219. [Google Scholar] [CrossRef]
  35. Long, C.N.; Dutton, E.G. BSRN Global Network Recommended QC Tests, V2.0. Available online: https://epic.awi.de/id/eprint/30083/1/BSRN_recommended_QC_tests_V2.pdf (accessed on 28 May 2024).
  36. ISO 9847:1992; Solar Energy—Calibration of Field Pyranometers by Comparison to a Reference Pyranometer. International Standard Organization: Geneva, Switzerland, 1992.
  37. Betti, A.; Blanc, P.; David, M.; Saint-Drenan, Y.-M.; Driesse, A.; Freeman, J.; Fritz, R.; Gueymard, C.; Habte, A.; Höller, R.; et al. Best Practices Handbook for the Collection and Use of Solar Resource Data for Solar Energy Applications, Report IEA-PVPS 16-04:2021, 3rd ed.; International Energy Agency: Paris, France, 2021; ISBN 978-3-907281-19-2. [Google Scholar]
Figure 1. Schematic for modeling GTI in rows that are not the front row.
Figure 1. Schematic for modeling GTI in rows that are not the front row.
Energies 17 03444 g001
Figure 2. Schematic of the various angles for the computation of shadows and obscuring of circumsolar radiation for (a) a segment of the panel being evaluated, (b) the back of the front panel, and (c) the ground between the rows of panels.
Figure 2. Schematic of the various angles for the computation of shadows and obscuring of circumsolar radiation for (a) a segment of the panel being evaluated, (b) the back of the front panel, and (c) the ground between the rows of panels.
Energies 17 03444 g002
Figure 3. Experimental setup for measuring global tilt irradiance for different positions.
Figure 3. Experimental setup for measuring global tilt irradiance for different positions.
Energies 17 03444 g003
Figure 4. Overview of the experimental setup including (a) the Évora–PECS station and (b) the pyranometers used for albedo computations.
Figure 4. Overview of the experimental setup including (a) the Évora–PECS station and (b) the pyranometers used for albedo computations.
Energies 17 03444 g004
Figure 5. Schematic for modeling GTI on the sensor in the experimental apparatus.
Figure 5. Schematic for modeling GTI on the sensor in the experimental apparatus.
Energies 17 03444 g005
Figure 6. Global tilted irradiance observed and modeled by the Modified Bugler model (first-row model for reference and comparison) and the developed model for period 1 (1 min timestep).
Figure 6. Global tilted irradiance observed and modeled by the Modified Bugler model (first-row model for reference and comparison) and the developed model for period 1 (1 min timestep).
Energies 17 03444 g006
Figure 7. Global tilted irradiance observed and modeled by the Modified Bugler model (first-row model for reference and comparison) and the developed model for period 12 (1 min timestep).
Figure 7. Global tilted irradiance observed and modeled by the Modified Bugler model (first-row model for reference and comparison) and the developed model for period 12 (1 min timestep).
Energies 17 03444 g007
Figure 8. Global tilted irradiance observed and modeled by the Modified Bugler model (first-row model for reference and comparison) and the developed model for period 19 (1 min timestep).
Figure 8. Global tilted irradiance observed and modeled by the Modified Bugler model (first-row model for reference and comparison) and the developed model for period 19 (1 min timestep).
Energies 17 03444 g008
Figure 9. Comparison between downscaled 10 min forecasts of the ECMWF/IFS and observations made at Évora–Verney of DNI, DIF, and GHI for forecast day 0 (the colormap represents the number of data points in each bin; bin size: 20 × 20 W/m2).
Figure 9. Comparison between downscaled 10 min forecasts of the ECMWF/IFS and observations made at Évora–Verney of DNI, DIF, and GHI for forecast day 0 (the colormap represents the number of data points in each bin; bin size: 20 × 20 W/m2).
Energies 17 03444 g009
Figure 10. Comparison between 10 min downscaled forecasts of the ECMWF/IFS and observations made at Évora–Verney of DNI, DIF, and GHI for forecast day 1 (the colormap represents the number of data points in each bin; bin size: 20 × 20 W/m2).
Figure 10. Comparison between 10 min downscaled forecasts of the ECMWF/IFS and observations made at Évora–Verney of DNI, DIF, and GHI for forecast day 1 (the colormap represents the number of data points in each bin; bin size: 20 × 20 W/m2).
Energies 17 03444 g010
Figure 11. Comparison between 10 min downscaled forecasts of the ECMWF/IFS and observations made at Évora–Verney of DNI, DIF, and GHI for forecast day 2 (the colormap represents the number of data points in each bin; bin size: 20 × 20 W/m2).
Figure 11. Comparison between 10 min downscaled forecasts of the ECMWF/IFS and observations made at Évora–Verney of DNI, DIF, and GHI for forecast day 2 (the colormap represents the number of data points in each bin; bin size: 20 × 20 W/m2).
Energies 17 03444 g011
Figure 12. Comparison between improved 10 min forecasts and observations made at Évora–Verney of DNI and DIF for forecast days 0, 1, and 2 (the colormap represents the number of data points in each bin; bin size: 20 × 20 W/m2).
Figure 12. Comparison between improved 10 min forecasts and observations made at Évora–Verney of DNI and DIF for forecast days 0, 1, and 2 (the colormap represents the number of data points in each bin; bin size: 20 × 20 W/m2).
Energies 17 03444 g012
Figure 13. Difference between mean bias errors of GTI from the developed model results using forecasted or experimental data as input for day 0 (negative differences in white and positive differences in blue).
Figure 13. Difference between mean bias errors of GTI from the developed model results using forecasted or experimental data as input for day 0 (negative differences in white and positive differences in blue).
Energies 17 03444 g013
Table 1. Structure position and number of 1 min data points for each testing period of the first-row tests.
Table 1. Structure position and number of 1 min data points for each testing period of the first-row tests.
Periodβ (°)Number of Data Points
120551
238582
350620
470545
590508
All-2806
Table 2. Structure position and number of 1 min data points for each testing period of the proposed model for rows other than the first.
Table 2. Structure position and number of 1 min data points for each testing period of the proposed model for rows other than the first.
PeriodD (m)β (°) c s (m)Number of Data Points
11.10380.08585
21.10380.58604
31.10381.081528
40.81380.182345
50.81380.38605
60.81380.58615
70.81380.78608
80.81380.98585
90.80500.15581
100.80500.352412
110.80500.55617
120.80500.75632
130.80500.95651
140.84200.261807
150.84200.46365
160.84200.66648
170.84200.86612
180.84201.06633
190.81380.18718
All---17,151
Shaded---2163
Unshaded---14,988
Table 3. Mean and standard deviation of the clearness index and GTI and metrics of GTI from the transposition models analyzed for each and all measuring periods of the first-row tests. The best performing model is represented in bold for each indicator and period and the clearness index data are repeated for each metric for better readability.
Table 3. Mean and standard deviation of the clearness index and GTI and metrics of GTI from the transposition models analyzed for each and all measuring periods of the first-row tests. The best performing model is represented in bold for each indicator and period and the clearness index data are repeated for each metric for better readability.
MetricPeriod k t GTILiu
Jordan
KlucherHay
Davies
ReindlMa IqbalModified Ma IqbalBuglerModified BuglerTemps CoulsonPerez
MeanStdMean
(W/m2)
Std
(W/m2)
R210.6970.100700.0322.00.99950.99920.99950.99950.99950.99950.99940.99950.99870.9994
20.7200.081702.1330.00.99910.99880.99910.99920.99910.99910.99900.99910.99850.9991
30.6740.087585.5297.90.99550.99410.99650.99670.99590.99590.99460.99570.83410.9975
40.6860.141474.8288.40.99460.99360.99450.99380.99440.99440.99420.99440.93170.9934
50.6760.102327.7202.80.99820.99630.99780.99530.99790.99780.99910.99870.98370.9993
All0.6900.105558.1323.60.99590.99690.99740.99750.99700.99700.99380.99470.92160.9978
MBE (W/m2)10.6970.100700.0322.0−44.4−28.8−36.6−36.4−38.6−38.458.129.8−43.7−31.5
20.7200.081702.1330.0−53.7−40.7−44.6−43.9−46.9−46.751.223.6−34.6−36.5
30.6740.087585.5297.9−66.5−32.4−48.6−42.4−53.9−53.523.05.3−157.4−24.1
40.6860.141474.8288.4−51.4−28.4−45.1−38.4−46.5−46.537.219.7−75.0−24.0
50.6760.102327.7202.8−23.2−2.0−27.3−18.3−25.8−26.257.944.6−19.6−10.3
All0.6900.105558.1323.6−48.7−27.1−40.9−36.4−42.9−42.844.823.8−68.7−25.6
RMSE (W/m2)10.6970.100700.0322.049.736.141.341.243.543.359.130.445.837.5
20.7200.081702.1330.056.143.546.345.648.848.553.925.836.138.9
30.6740.087585.5297.971.938.051.745.157.657.227.415.2181.826.9
40.6860.141474.8288.456.534.848.742.750.450.342.926.794.830.8
50.6760.102327.7202.824.911.828.321.727.027.362.648.241.911.3
All0.6900.105558.1323.654.935.044.540.747.147.050.230.5100.230.9
GPIAll0.6900.105558.1323.61.0681.8481.4121.5691.3241.3291.1951.959−1.0001.953
Table 4. Mean and standard deviation of the clearness index and metrics of GTI from the Modified Bugler model and proposed model for other rows.
Table 4. Mean and standard deviation of the clearness index and metrics of GTI from the Modified Bugler model and proposed model for other rows.
Period k t GTIModified Bugler ModelDeveloped Model for Other Rows
MeanStdMean (W/m2)Std (W/m2)R2MBE (W/m2)RMSE (W/m2)R2MBE (W/m2)RMSE (W/m2)
10.6870.077566.1342.10.9905100.0106.10.9921−31.139.5
20.6750.108646.2324.40.998875.676.60.9987−38.544.8
30.5840.204456.2338.00.941357.996.20.9410−24.782.0
40.6410.191176.2126.30.0873521.5651.90.9844−15.722.7
50.7040.095627.8335.90.998476.478.50.9990−38.642.3
60.6900.103599.6356.20.999568.470.80.9991−39.442.3
70.6720.098582.9345.20.999669.371.00.9996−39.443.8
80.6840.104634.2335.20.999036.637.60.9988−69.577.3
90.5400.206248.4187.00.1385290.2376.00.3225−81.1161.5
100.5900.190420.5307.10.9769103.4112.00.9944−42.350.6
110.6740.145548.9307.60.998661.262.10.9962−51.657.8
120.6220.157432.0270.40.994952.455.20.9946−54.062.1
130.5610.134436.4272.70.996237.940.60.9961−44.758.7
140.3510.192124.579.550.4669250.6304.90.461087.0141.2
150.6360.154440.7300.00.998028.330.30.9987−62.769.6
160.5240.275481.7370.50.979235.161.70.9793−50.378.2
170.6950.172679.9360.30.990741.550.60.9909−65.978.0
180.6770.147643.9334.10.996637.442.40.9964−63.469.0
190.6440.177248.6232.20.0372373.2570.20.9063−12.475.5
All0.6380.188323.8304.60.4698381.2301.20.9543−12.976.8
Shaded0.6820.197203.7169.20.0028717.3775.00.2585−28.984.2
Unshaded0.6320.186421.6336.70.904788.9131.00.9552−26.575.7
Table 5. MBE and weight of each component of GTI on the reduction in bias error when comparing the developed model with the Modified Bugler model.
Table 5. MBE and weight of each component of GTI on the reduction in bias error when comparing the developed model with the Modified Bugler model.
PeriodMBE Modified Bugler (W/m2)MBE Developed Model (W/m2) w I b (%) w D c s (%) w D i s o (%) w I r e f l (%)
Unshaded88.9−26.50.00.556.642.9
Shaded717.3−28.966.52.628.72.2
All381.2−12.98.40.853.137.7
Table 6. DNI and DIF forecast results for all data of first-row tests and for all shaded and unshaded data of the developed model testing (MBE, MAE, and RMSE in W/m2).
Table 6. DNI and DIF forecast results for all data of first-row tests and for all shaded and unshaded data of the developed model testing (MBE, MAE, and RMSE in W/m2).
Data and
Variable
Day 0Day 1Day 2
R2MBEMAERMSER2MBEMAERMSER2MBEMAERMSE
First-row tests
AllDNI0.4859−38.587.2132.10.2568−41.5107.4169.20.2373−99.8151.2225.7
DIF0.5258−1.644.262.10.55404.041.360.30.53827.748.662.4
Inner-row tests
AllDNI0.538634.0146.3208.70.545839.4145.1208.60.493031.8156.0224.6
DIF0.425266.466.594.60.472064.064.091.50.431668.368.395.2
ShadedDNI0.6166126.4185.6267.40.6488134.8183.6263.80.5840121.6196.2272.8
DIF0.4670−49.792.8126.80.5025−53.197.1126.90.4255−44.4101.0131.0
UnshadedDNI0.535116.9139.0196.00.538921.8138.0196.70.487615.2148.5214.6
DIF0.3884−8.961.687.40.4437−17.657.883.30.4038−15.662.387.1
Table 7. GTI results of the developed model for rows that are not the front row for each day of forecast horizon using DNI and DIF forecasts (MAE and RMSE in W/m2; MAPE and RMSPE in %).
Table 7. GTI results of the developed model for rows that are not the front row for each day of forecast horizon using DNI and DIF forecasts (MAE and RMSE in W/m2; MAPE and RMSPE in %).
DataDay 0Day 1Day 2
R2MBEMAERMSER2MBEMAERMSER2MBEMAERMSE
First-row tests
All0.96466.734.548.70.960715.539.455.00.800171.567.3128.5
Inner-row tests
All0.7197−31.5100.6164.00.7162−32.1101.1165.60.6727−37.2106.8179.4
Shaded0.0074−115.3152.4285.10.0242−116.7152.4281.80.0078−113.9151.3284.5
Unshaded0.8031−16.091.0129.80.7907−16.591.6133.60.7318−23.098.6152.2
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

Pereira, S.; Canhoto, P.; Salgado, R. Prediction of Global Solar Irradiance on Parallel Rows of Tilted Surfaces Including the Effect of Direct and Anisotropic Diffuse Shading. Energies 2024, 17, 3444. https://doi.org/10.3390/en17143444

AMA Style

Pereira S, Canhoto P, Salgado R. Prediction of Global Solar Irradiance on Parallel Rows of Tilted Surfaces Including the Effect of Direct and Anisotropic Diffuse Shading. Energies. 2024; 17(14):3444. https://doi.org/10.3390/en17143444

Chicago/Turabian Style

Pereira, Sara, Paulo Canhoto, and Rui Salgado. 2024. "Prediction of Global Solar Irradiance on Parallel Rows of Tilted Surfaces Including the Effect of Direct and Anisotropic Diffuse Shading" Energies 17, no. 14: 3444. https://doi.org/10.3390/en17143444

APA Style

Pereira, S., Canhoto, P., & Salgado, R. (2024). Prediction of Global Solar Irradiance on Parallel Rows of Tilted Surfaces Including the Effect of Direct and Anisotropic Diffuse Shading. Energies, 17(14), 3444. https://doi.org/10.3390/en17143444

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop