Investigation of the Effect of Albedo in Photovoltaic Systems for Urban Applications: Case Study for Spain

: Rooftop photovoltaic generation can help cities become key players in the transition to clean energy. The optimal solar photovoltaic production on rooftops depends on two angles: tilt angle and azimuth angle. It is accepted in all studies that the ideal orientation of photovoltaic modules is toward the south (north) in the northern hemisphere (south). In contrast, the determination of the optimum tilt angle is more complex, and there are different equations for its calculation. Most of these equations do not take albedo into account. In this work, 47 Spanish province capitals representing the most populated cities have been studied with different equations for the calculation of the optimum annual tilt angle (Technical report by the Spanish Institute for the Diversiﬁcation and Saving of Energy (IDAE), Lorenzo’s and Jacobson’s equation) and different types of albedo. Accounting for the geographical and the meteorological conditions of the cities, we analyzed the optimum tilt angle through a Mathematica© optimization code. The inﬂuence that different variables have on optimum tilt angle has been quantiﬁed by means of the term relative energy harvested. The use of the equations as a function of latitude increases the annual relative energy harvested by increasing the albedo. When the albedo is 0.2, the annual relative energy harvested is very similar in all equations. Comparing to the method that maximizes the total irradiation incident on a tilted surface, the minimum and maximum value of the percentage of relative energy harvested per year were 0.01 and 2.50% for the IDAE guideline, 0.00 and 2.38% for Lorenzo’s equation, 0.00 and 2.46% for Jacobson’s equation. A simpliﬁed polynomial regression model to estimate optimum tilt angle as a function of latitude, altitude and albedo has been proposed as well.


Introduction
Electricity generation is the engine of the economic and social development of society. The need for large amounts of electrical energy has led to unprecedented environmental issues due to the use of fossil fuels and therefore the increasing of greenhouse gas emissions. There is a global consensus to reduce the use of fossil fuels and replace it with renewable energy sources of energy [1]. It is the called clean energy transition. A long-term sustainable energy model is therefore necessary to meet the commitments made by governments at the global level.
Solar energy is a renewable energy source that can be used for electricity generation. Electricity can be generated directly by photovoltaic systems [2] or by converting the Sun's energy into high-temperature heat (concentrating solar power plants) [3]. From the point of view of their use in cities, the first technology is the right one.
Cities are major players in the environmental issues considered as examples of unsustainability [4]. Buildings in the European Union (EU) are responsible for 40% of final energy consumption and 36% of greenhouse gas emissions [5]. For this reason, cities are at the center of the implementation of EU policies [6] such as the promotion of low and zero-energy buildings. In this scenario, cities need profound changes toward an ecological future [7].
The decentralized electricity production (DEP) systems can help make cities key players in the clean energy transition. This new paradigm of electricity production is revolutionizing the power generation sector [8]. In this sense, Gómez et al. [9] introduce the concept of "prosumer", or a customer who consumes and can sell generated electricity. Electricity production by means of building-integrated photovoltaics systems and buildingapplied photovoltaics systems are two solutions used in DEP systems.
The installation of PV systems on rooftops is rapidly growing [10], as they produce benefits from various points of view: health [11], environmental [12], economic [13], and social [14]. These systems have been studied by several authors [15,16].
The optimal solar PV production on rooftops depends on two angles [17]: the angle between the plane of the PV module and the horizontal plane (called tilt angle β) and the angle between the projection on a horizontal plane of the normal to the PV module and south direction (called azimuth angle γ). These angles not only have a direct impact on the amount of solar irradiance incident on a PV module [18] but also influence the surface area required for the installation of the PV system [19]. The optimum azimuth angle for tilted surfaces is facing due south (north) in the northern (southern) hemisphere [17]. Many researchers have studied procedures to find the optimum tilt angle, grouping them into two categories.
The first category uses equations based on the latitude, in which the optimum tilt angle is taken as the latitude plus or minus a specific value obtained via analytic methods (e.g., regression analysis). These equations are very easy to use but prone to error. In addition, these equations have been determined for albedos close to 0.2. Yet in the building sector, materials with albedos greater than 0.2 can be used. Therefore, the use of these equations leads to errors in the estimation of the energy absorbed by the PV modules. These equations are used as an installation guide of PV modules in practical engineering applications.
A large number of studies belong to this first category. Ullah et al. [20] showed a detailed summary of these studies. These studies were carried out with the aim of obtaining formulas to determine the optimum tilt angle of PV modules. For locations with λ < 45 • (where λ is the latitude of the site), Lorenzo [21] obtained the optimal tilt of solar collector as a function of latitude: β opt = 3.7 + 0.69 · |λ|. Jacobson et al. [22] carried out an analysis that includes a 3rd degree polynomial: β opt = 1.3793 + λ(1.2011 + λ(−0.014404 + 0.000080509λ)). These two studies do not take into account the variation of the albedo. Nicolás-Martín et al. [23] obtained the annual optimum tilt angle as a function of latitude, diffuse fraction and albedo. To do this, 14, 468 sites spread across the globe from the One Building database were used. Their study did take into account the albedo variation.
The second category cover procedures that maximize the total solar irradiance incident on a tilted surface over a period of time. These procedures are more accurate, yet they are more complex. Soulayman et al. [24] presented a modified general algorithm to determine the annual optimum tilt angle over all mid-latitude zone. The latitude of these sites varied between 23.45 ( • ) N and 43.45 ( • ) N and between 23.45 ( • ) S and 43.45 ( • ) S. This work suggested following equation β opt = 1.171 + 0.916 · λ for the mid-latitudes zone of the northern hemisphere. Al Garni et al. [25] investigated the annual optimum tilt angle of PV modules for 18 sites of Saudi Arabia. The latitude of the cities varied between 18.22 ( • ) and 31.02 ( • ). For these sites, the optimal tilt was only slightly higher than the latitude. Ullah et al. [20] investigated the annual optimum tilt angle for PV modules for several cities in Pakistan. This study was experimentally validated for one of the cities studied. Barbón et al. [2], applying Cavaleri's principle, investigates the annual optimum tilt angles for PV modules for 39 locations in the northern hemisphere. These four studies do not take into account the variation of the albedo. Table 1 shows a summary of the advantages and disadvantages of procedures to find the optimum tilt angle. The objective of this paper is to study the influence that several variables have on the optimum tilt angle for 47 Spanish province capitals and to propose equations that estimate the optimum tilt angle that take into account the variations of the albedo. The notable difference with the methods in the first category is that our procedure allows the albedo value to be varied. The analysis is performed using an algorithm that maximizes the solar irradiation reaching the tilted surface for a given period of time. The annual optimum tilt angle for each of the cities is computed using Cavaleri's principle. The influence that the different variables have on optimum tilt angle is quantified by means of the term relative energy harvesting (REH). REH can be defined as the ratio of the energy collected through the direct application of each model to the amount obtained by application of Cavaleri's principle.
The specific contributions of this study can be summarized as follows: (i) We analyze the equations proposed by the Spanish regulations for the determination of the annual optimum tilt angle. (ii) We determine the error caused through the use of the equations based on the latitude that do not take into account the variation of the albedo. (iii) We propose a method to compute the annual optimum tilt angle that accounts for the albedo variations. (iv) We conduct a detailed analysis of the relative energy harvesting as a function of different variables. (v) Finally, a simplified polynomial regression model to estimate optimum tilt angle as a function of latitude, altitude and albedo is proposed.
This paper is organized as follows: the geographic characteristics of the Spanish province capitals are presented in Section 2. Section 3 presents the solar irradiance estimation model and the procedure followed to calculate the amount of total irradiation on a tilted plane. The methods used to calculate the optimum tilt angle are described in Section 4. Section 5 presents the results of the study. Finally, Section 6 summarizes the main contributions and conclusions of the paper.

A Case Study in Spain
A total number of 47 province capitals comprising 29.8% of the Spanish population have been analyzed, representing the most populous cities of Spain. These locations are given in Table 2 together with their main geographical characteristics. The geographical distribution of province capital cities analyzed in this paper is illustrated in the map presented in Figure 1.  These cities are located between latitudes 35 • N and 45 • N; therefore, there are significant seasonal variations. This is an important fact to keep in mind. Figure 2 shows a solar map provided by the SolarGIS organization of the global horizontal solar irradiation for Spain. SolarGIS offers a reliable solar resource database with local geographical data of high accuracy and a spatial resolution to the site specifications [26].

Solar Irradiance Estimation Model
The methods used to estimate the optimum tilt angle based on maximizing the solar irradiation on a tilted surface are highly sensitive to the solar irradiance estimation model. This is because the solar irradiance is specific for each geographical location, and it is affected by the distribution of the local cloud cover. Therefore, a set of precise data about the solar irradiance arriving at Earth for each particular site is needed. These data, provided by the ground-level meteorological stations, are the global and diffuse solar irradiance on horizontal surfaces. However, the number of these meteorological stations in the world is not very high. Therefore, models are used to estimate solar irradiance in the absence of meteorological stations. The results of these models are only approximate, but they are accepted by the scientific community. Nonetheless, it is necessary to be careful when using these models and take into account that the accuracy of a model varies depending on the specific site where it is used [27]. In recent decades, there are many different models published in the technical literature, such as the clear-sky models [28], the satellite-based models [29], the temperature-based methods [30], etc.
In the present study, it is necessary to estimate the hourly beam and diffuse solar irradiance on a horizontal surface of a geographical location. For this reason, the procedure presented in [31] has been used. This procedure uses the term 'adjusted solar irradiance' to refer to the solar irradiance that takes into account the weather conditions of the site. It adapts the clear sky models, Hottel's model [32] and Liu's and Jordan's model [33], to the weather conditions of the site to estimate the adjusted hourly beam and diffuse solar irradiance on a horizontal surface. Fourier series approximations are used to correct the clear-sky models. Several works have already applied this method [2,34,35], because this procedure demonstrates its accuracy and its application to different climates compared with the actual data obtained from ground-level meteorological stations (WRDC database [36]).

Estimation of the Amount of Total Irradiation on a Tilted Plane
Most authors use the same procedure to calculate the total solar irradiance (I t ) on a tilted surface. The procedure in obtaining its components separately: the beam (I bt ), the diffuse (I dt ), and the ground reflected (I rt ) irradiances. These components are defined in [17], and its calculation requires estimating beforehand the beam and the diffuse irradiance on a horizontal plane.
The amount of beam irradiance on a tilted surface can be estimated from the geometric relation between the horizontal plane and the tilted surface [17], which is expressed as: is the azimuth angle, n (day) is the day of the year and T (h) is the solar time. The incident angle θ i can be determined following [17].
The diffuse irradiance on a tilted surface has three components: isotropic, circumsolar, and horizon brightening irradiances [17]. Depending on the components considered for calculation, two models may be used: isotropic and anisotropic model. Liu and Jordan's model [37] is an isotropic model widely used in the specialized literature [38,39]. In the comparison between the isotropic and anisotropic models presented by Mehleri et al. [40], it is concluded that Liu and Jordan's model obtains accurate results. According to this model, the diffuse irradiance on a tilted surface can be defined as: where I dh (W/m 2 ) is the diffuse irradiance on a horizontal plane, β ( • ) is the tilt angle, n (day) is the day of the year and T (h) is the solar time.
Although it is impossible to exactly compute the ground-reflected irradiance due to the many factors contributing to it [41], most works assume that the reflection on the ground of the beam and diffuse solar irradiances is isotropic [17,37,42,43], and they use (3) to compute the ground-reflected irradiance: where I bh (W/m 2 ) is the beam irradiance on a horizontal plane, I dh (W/m 2 ) is the diffuse irradiance on a horizontal plane, β ( • ) is the tilt angle, n (day) is the day of the year, T (h) is the solar time, and ρ g (dimensionless) is the ground reflectance or albedo. Typical ρ g values for different ground surfaces have been computed by [44][45][46]. Some of these values are: for weathered concrete ρ g = 0.22, for dark surfaces of buildings (red brick, dark paints, etc.) ρ g = 0.27, and for light surfaces of buildings (light brick, light paints, etc.) ρ g = 0.60. A value of 0.2 is commonly adopted if no information is available about ground surface [47]. So, the total solar irradiance on a tilted surface can be calculated as: Equation (4) is fitted to each location and climate data and, therefore, fitted under real weather conditions.
In order to compute the adjusted total solar irradiation on a tilted surface for each day of the year, the ordinary method of adjusted solar irradiance integration between sunrise and sunset can be expressed as [2]: where H t (Wh/m 2 ) is the adjusted total solar irradiation on a tilted surface, I t (W/m 2 ) is the adjusted hourly distribution of total solar irradiance on a tilted surface, β ( • ) is the tilt angle, γ.( • ) is the azimuth angle, n (day) is the day of the year, T (h) is the solar time, T R (h) is the sunrise solar time, and T S (h) is the sunset solar time.
The adjusted annual solar irradiation on a tilted surface, H a t (β, γ), for different values of tilt and azimuth angles can be expressed as [2]: The ideal position of a PV module is defined by the optimum tilt and azimuth angles. As it is well known, the optimum azimuth angle is 0 • in the northern hemisphere [48,49]. Therefore, Equation (5) can be simplified, obtaining the two-variable function, H t (n, β, 0). The plot of this function is shown in Figure 3 for Madrid (Spanish capital). As the ideal position can be referred to a period of time (year, month, day), H t (n, β, 0) can be used to obtain the ideal position for a certain period of time.

Procedure
The following considerations have been made in this study: (i) The effect of the weather conditions for each city under study is taken into account with the method proposed by [31]. (ii) The Mathematica © code uses the satellite-derived PVGIS data [50] for each city under study as inputs of monthly-averaged beam and diffuse solar irradiation. (iii) Any albedo value can be used. (iv) The study was carried out in Spain. However, the procedure can be applied anywhere in the world.
The optimum tilt angle and the total incident energy are outputs of the algorithm.

Computation of the Optimum Tilt Angle
In what follows, we shall look at several proposals to compute the optimum tilt angle for non-tracker-based panels. Assuming the PV modules are south oriented, the technical literature provides several equations for calculating the annual optimum tilt angle based on the latitude of the site.

Technical Report by the Spanish Institute for the Diversification and Saving of Energy (IDAE)
A Technical Report [51] by IDAE requires that the annual optimum tilt angle be determined by Equation (7): where λ is the latitude ( • ).

Other Equations
For locations with λ < 45 • , Lorenzo [21] proposes the following equation for calculating the annual optimum tilt angle: β opt = 3.7 + 0.69 · |λ| (8) where λ is the latitude ( • ). This model is has been used in [52,53]. Jacobson's model [22] has also been used extensively [23,25,54,55]. This model proposes the following equation for calculating the annual optimum tilt angle: where λ is the latitude ( • ). Equations (7)-(9) do not take into account other variables that influence the annual solar irradiation on a tilted surface, such as albedo.

Cavaleri's Principle
The previous methods are based on statistical procedures and do not include the effect of albedo variation. In order to take into account the effect of albedo in the determination of the optimum tilt angle, we present an application of Cavaleri's principle.
Another procedure to obtain the optimum tilt angle is based on maximizing the solar irradiation falling on a tilted surface over a period of time. In this case, the period of time is a year. For this purpose, Cavaleri's principle can be applied [2,34].
The volume underneath the graph of the equation H t (n, β, 0) is given by the double integral [2]: where D is the rectangle D : [1, 365] × [0, 90]. The application of integral calculus to Cavaleri's principle yields [2]: The interval [0, 90] can be discretized as the integral for each of the values calculated. This optimization procedure seeks β y opt such that [2]: A flowchart outlining the proposed methodology is shown in Figure 4. Figure 5 shows the representation of optimum tilt angle in Madrid for each value of ρ g . This graph shows that the annual optimum tilt angle is higher as the albedo increases.
This procedure can be used in other parts of the world, both in the northern and southern hemispheres. The proposed procedure is not limited by the latitude of the site. Take for example Cairo, Egypt (Latitude: 30 • 29 24 N, Longitude: 31 • 14 38 W, Altitude: 41 (m)). Figure 6 shows the representation of optimum tilt angle in Cairo for each value of ρ g .

Energy Harvesting
Relative energy harvesting (REH) can be used in the evaluation of the equation for the determination of the optimum tilt angle [24]. In this context, it is reasonable to introduce the "maximum" value of the study, assuming the PV modules are south oriented. For this purpose, we provide a baseline "maximum" value using single-axis trackers aligned to the east-west axis. Duffie's [17] formula for this tracker is used, which yields: β = arctan(tan θ z | cos γ s |) (13) where, in this case, the azimuth angle of the surface depends on γ S as follows [17]:      Table 3 shows the maximum annual irradiation for single-axis trackers aligned to the east-west axis for each value of ρ g . These are maximum values which obtain a tilted surface assuming the PV modules are south oriented. Table 3. Maximum annual irradiation (MWh/m 2 ) on a tilted plane.

Id
City ρ g = 0.2 ρ g = 0.4 ρ g = 0.6 ρ g = 0.8  REH is then calculated as the difference between the energy harvested according to each specific model under study and the single-axis trackers aligned with the east-west axis, as a percentage of energy: where the subindex * stands, as above, for the corresponding equation (IDAE recommendation, Lorenzo's, or Jacobson's equations).
The REH can be also calculated as the difference between the energy harvested after application of the model under study and the method based maximizing the total irradiation incident on a tilted surface applying Cavaleri's principle, as a percentage of energy: where the subindex * stands, as above, for the corresponding equation (IDAE recommendation, Lorenzo's equation, Jacobson's equation). Based on (15) and (16), one can obtain reliable information regarding the model from a practical point of view. In addition, these equations will allow us to evaluate the validity of the models when the albedo is different from the usual.

Results and Discussion
Equations (8) and (9) are widely used to estimate the amount of total irradiation on a tilted plane. Equation (7) is recommended by a technical report by IDAE. The aim of this section is to assess these equations depending on the albedo variations. Based on the equations presented in Section 2, each type of solar irradiance is calculated with the effect of the weather conditions. A Mathematica © optimization code was used for computing the direct, diffuse, and reflected components of the solar irradiance. The satellite-derived PVGIS data [50] were used to obtain the monthly-averaged beam and diffuse solar irradiation of each city under study. The method proposed by [31] was used to account for the effect of weather conditions.

Equations as Function of Latitude
We compare the yearly solar irradiation using the models in (7)-(9) to determine the optimum tilt angle, taking as a baseline the single-axis tracker. Figure 7 shows a graphical representation of the yearly REH 1 between the use of the annual optimum tilt angle and a single-axis tracker for each value of ρ g for which Equation (15) has been used. The charts show the loss of harvested solar energy with respect to the baseline for four different albedo values. The following conclusions can be inferred: (i) The annual solar irradiation for each city is consistently higher as the albedo increases. (ii) For ρ g = 0.2, the three models closely match. Lorenzo's and Jacobson's equations obtain very similar total annual irradiation values. The maximum REH of total yearly irradiation is −6.76% in the case of IDAE's guidelines, −6.66% for Lorenzo's equation, and −6.65% for Jacobson's equation. The deviations are not greater than 0.11%. The reflectivity of most types of ground surfaces is rather low; therefore, the contribution of this type of solar irradiance falling on a PV module is low. When no specific albedo values are available, the value of 0.2 is typically assumed [47]. Therefore, in this case, these equations are valid.

Method Based on to Maximize the Total Irradiation Incident on a Tilted Surface
Equations (7)-(9) depend on latitude but not on albedo. In contrast, the procedure that maximizes the total irradiation incident on a tilted surface takes into account the albedo. Figure 8 shows the influence of albedo on the annual optimum tilt angle in Madrid. The annual optimum tilt angle is higher as albedo increases. The equations to determine the optimum tilt angle for Spain are compared with the method based on maximizing the total irradiation incident on a tilted surface applying Cavaleri's principle. Figure 9 shows a graphical representation of the yearly REH 2 comparing the use of different equations to determine the optimum tilt angle and the application of the Cavaleri's principle for each value of ρ g , for which (16) was used. Figure 9 suggests the following conclusions: (i) The method that maximizes the total irradiation incident on a tilted surface obtains the best results. (ii) For ρ g = 0.2, the three equations closely match with the method that maximizes the total irradiation incident on a tilted surface. (vi) The greater the albedo, the greater the error introduced by the three models, since they only account for the latitude of the place. The results suggest that the proposed optimization method is able to take into account the variation of the albedo. (vii) With increased albedo, the IDAE guideline obtains the worst results (viii) It is recommended to use Lorenzo's equation and Jacobson's equation for the cities studied for ρ g = 0.2. Figure 9. Yearly REH between the equations as function of latitude and the method based on maximizing the total irradiation incident on a tilted surface.

Simplified Model for the Annual Optimum Tilt Angle
Previous studies have clearly shown the need to take into account the influence of albedo when calculating the optimum tilt angle. The method based on Cavaleri's principle presents a high accuracy in the calculations, as it takes into account not only the albedo but also the direct and diffuse irradiance regarding the meteorological conditions of each location. However, we are aware that their use requires the use of software that may not be available to all users. Therefore, in this section, we present a simplified model (similar to Lorenzo's or Jacobson's) that produces good approximation to the value of the optimum annual tilt obtained with our method.
We proceeded with a Monte Carlo-like experiment, using the latitude and altitude data of the 47 cities-as well as f our cases of varying albedo-as samples. After applying Cavalieri's principle to these samples, we obtained the optimum tilt, which we employed as our response variable in a regression problem, in which the regressors or independent variables were precisely the latitude, altitude, and albedo. We then proceeded with a symbolic regression through genetic programming based on Python's project gplearn [56]. Symbolic regression is a machine learning technique employed to identify underlying mathematical expressions based on a catalog of possible functions. It begins by constructing a population of naive random formulae to represent a relationship between regressors and a response variable in order to predict new data. Each successive generation of programs then evolves from the previous one by selecting the fittest individuals from the population to undergo genetic operations. Genetic programming represents solutions to the regression problem as trees (by contrast, genetic algorithms represent them as strings). Using this tree representation is quite convenient in our application, since mathematical expressions can be represented by trees.
We performed a symbolic regression analysis by testing arithmetic functions along with trigonometrical and power functions as possible fittings. After ten thousand iterations, the best parsimonious fit turned out to be a simple linear regression model. We comment on two alternatives next.
(i) A first linear model, a function of latitude λ and albedo ρ, is: The goodness of fit obtained, as measured by the coefficient of determination, gives a value of R 2 = 0.9537. In Figure 10. we see the cloud of 4 × 47 = 188 sample points and the fitting surface (the plane in this case) β(λ, ρ). (ii) A second model, incorporating altitude A, is: The goodness of fit is very similar to the previous one, obtaining now R 2 = 0.9553. A noteworthy feature is that the albedo term is the same in both models. The altitude of the location only plays a marginal role in accuracy.
After extensive testing, we have found that it does not pay to use higher degree models (quadratic or cubic), as they unnecessarily complicate the model with little improvement in goodness of fit. For example, the quadratic model and the cubic model we obtain are: β(λ, ρ) = −110.245 + 6.48852λ − 0.0753977λ 2 + 11.0489ρ + 14.5745ρ 2 (19) β(λ, ρ) = 204.222 − 17.1264λ + 0.513803λ 2 − 0.00489214λ 3 + 14.9574ρ + 5.79787ρ 2 + 5.85106ρ 3 (20) obtaining R 2 = 0.9658 and R 2 = 0.9659, respectively. We have also investigated the possibility of finding non-linear models for the fit. By way of example, we present two of the models of this type that have given us the best results, combining precision and simplicity. On the one hand, the simplest non-linear model: β(λ, ρ) = 25.6067 + 0.634111λρ (21) which with the same variables (λ and ρ) gives an R 2 = 0.9541. On the other hand, we also tested a more sophisticated model that incorporates the ratio between the adjusted annual beam irradiance on a horizontal plane, I a bh (kWh/m 2 ) and the adjusted annual diffuse irradiance on a horizontal plane I a dh (kWh/m 2 ). In this case, we obtained as the best fit:  (22) and the goodness of fit is very similar, now obtaining R 2 = 0.9596.
To deepen this analysis, the key factor is to study the seasonality of the irradiances of each location. The study is very complex and again does not produce a substantial improvement of the fit.

Conclusions
In this paper, the optimum tilt angles of photovoltaic systems for urban applications have been analysed for 47 cities in Spain through a Mathematica © optimization code. These cities are the Spanish province capitals and represent the most populous cities. The analysis was carried out for fixed tilt angle and different types of albedo, taking into account both the geographical and the meteorological conditions of the sites. The yearly energy harvesting has been evaluated using: (i) models based on latitude (IDAE's guideline, and Lorenzo's and Jacobson's models); (ii) equations as function of latitude (IDAE guideline, Lorenzo's equation, Jacobson's equation) versus a method that maximizes the total irradiation incident on a tilted surface. The summary of the results is as follows: (i) The use of the equations as a function of latitude increases the annual REH by increasing the albedo. When the albedo was 0.2, the annual REH was very similar for all three equations. Lorenzo's and Jacobson's models did not appreciably change the annual REH in the range of albedos studied. (ii) With regard to single-axis tracker, the minimum and maximum value of the percentage of REH per year obtained were 6.76 (ρ g = 0.2) and 8.19% (ρ g = 0.8) for IDAE's guideline, 6.66 (ρ g = 0.2) and 7.66% (ρ g = 0.8) for Lorenzo's model, 6.65% (ρ g = 0.2) and 7.46% (ρ g = 0.8) for Jacobson's model. (iii) The annual REH was increased when the equations as a function of latitude were used, compared with the method that maximizes the total irradiation incident on a tilted surface. This increase is stronger for albedo values equal to 0.8. (iv) The yearly optimum tilt angle is higher as albedo increases. Therefore, equations that only depend on latitude are less accurate as the albedo increases. With regard to the method that maximizes the total irradiation incident on a tilted surface, the minimum and maximum value of the percentage of REH per year obtained were 0.01 and 2.62% for IDAE guideline, 0.00 and 2.38% for Lorenzo's equation, and 0.00% and 2.46% for Jacobson's equation. (v) We have obtained two linear models that account for the latitude, albedo and the altitude of the location to predict the optimum tilt angle. We believe that these formulae show great simplicity of use with acceptable accuracy.
We believe that our methodology can be used to make optimal decisions in the choice of tilt angle in photovoltaic modules as a function of albedo, obtaining important benefits from the point of view of total energy absorption. Future work will consist of an experimental study of the effect of albedo on the equations used in the determination of the optimum tilt angle.