European Union Green Deal and the Opportunity Cost of Wastewater Treatment Projects

: The European Union Green Deal aims at curbing greenhouse gas emissions and introducing clean energy production. But to achieve energy efﬁciency, the opportunity cost of different energies must be assessed. In this article, two different energy self-sufﬁcient systems for wastewater treatment are compared. On the one hand, high-rate algal ponds system (HRAP) is considered; on the other hand, a conventional activated sludge system (AS) which uses photovoltaic power (PV) is studied. The paper offers a viability analysis of both systems based on the capacity to satisfy their energetic consumption. This viability analysis, along with the opportunity cost study, will be used in the article to compare these two projects devoted to the treatment of wastewater. In order to assess viability, the probability of not achieving the energy consumption threshold at least one day is studied. The results point that the AS+PV system self-sufﬁciency is achieved with much lesser land requirements than the HRAP system (for the former, less than 6500 m 2 , for the latter 40,000 m 2 ). However, the important AS capital cost makes still the HRAP system more economic, although storage provides a great advantage for using the AS+PV in locations where a lot of irradiance is available.


Introduction
In the last decade, international institutions have taken on a strong commitment to achieve a climate neutral economy with the Horizon 2030 [1]. For this commitment, the European Union (EU) passed the EU Green Deal in 2019 aiming at promoting the clean energy production [2,3]. In particular, growing energy demand and water consumption have increased concerns about efficient wastewater treatment [4][5][6]. But for this basic utility at least two possible systems are used. One is the activated sludge process (AS), which entails the biological reduction of organic matter in the wastewater. The oxidation and digestion of carbonaceous biological matter using aerobic micro-organisms in wastewater remove the non-active microbes from the system and provide an effective reduction of pollutant parameters [7,8]. This system has been installed for various municipalities and industries [9]. Another system of wastewater treatment is the high rate algal ponds system (HRAP) that have great potential for biofuel production [10][11][12] when climate is favourable, since the costs of algal cultivation and harvest for biofuel production are covered by the wastewater treatment function. It can be used to provide community-level energy supply. Microalgae are promising alternative renewable sustainable energy sources as they produce large amounts of biomass which can be used for production of "third-generation Besides, more environment-concerned storage design is needed and then, natural biomass derived carbons are excellent alternatives for substituting conventional carbon materials toward a wide range of applications [48,49]. Some biomass is even recycled from the agricultural daily wastes [50]. In particular, the HRAP system can achieve energetic self-sufficiency from its own biomass. But the AS system needs an additional facility that provides energy and, within the Green Deal commitment, this energy is better provided with a renewable technology (such as photovoltaic or wind energy). In this work, the option of complementing the AS system with photovoltaic technology is addressed, leaving for future works the possibility of studying other renewable sources, such as the small wind energy system. The availability of photo-synthetically active radiation (PAR) ground measurement is scarce worldwide. But the Project CGL2016-79284-P AEI/FEDER/UE has made it possible to have an ample network of such measurements in the peninsular Spanish territory, covering different types of climates, such as the Arid Mediterranean and the Oceanic climates.
According to the previous considerations, the main aim of this work is to compare the viability of wastewater treatment projects under outdoor conditions by HRAP and by AS conventional systems, reinforcing this last system with a photovoltaic facility (PV). Indeed, as previously mentioned, there are several studies devoted to the comparison between both systems, but in these studies the inclusion of a PV system is not considered. Thus, the novelty of the article is the comparison between both systems, which must have selfconsumption capacity. The different needs of space and energy consumption of these two wastewater treatment systems determines the assessment of their viability considering land requirement and consumption. Another contribution of the paper is the conceptualization of viability as the capacity of the project self-sufficiency (when the energy consumption associated to the system is covered by the energy supplied by the system). This viability is assessed from the probability of meeting certain consumption thresholds considering two cases: absence of storage and full capacity of storage. Obviously, self-sufficiency does not always imply cheaper energy prices, but energy needs for wastewater treatment projects must be maximised.
The work is organized as follows: first, the models performed to estimate the algal productivity and the PV power have been described, as well as the stations where the study takes place. Second, the methodology of the work is presented. Next, the results are shown and discussed. In this sense, a large analysis is developed, according to the PV installation area requirement, where systems with and without storage are compared. Finally, an opportunity cost analysis is included.

Materials and Methods
The two wastewater systems compared in this work (HRAP and AS+PV) are designed to serve a population equivalent to 10,000. The HRAP and AS systems used have been taken from the scenarios defined in [51]. According to these scenarios, HRAP only consumes 0.06 kWh m −3 , while the AS conventional system consumes 0.89 kWh m −3 . The viability of covering these consumptions in a self-sufficient way is analysed: the necessary energy for the AS system is provided by a PV system, while the energy for the HRAP system will come from algal productivity. In this sense, two aspects must be taken into account: on the one hand, the areas which must be covered by algal ponds and by PV panels in order to satisfy the respective demands; on the other hand, the economic aspect of both systems.

Site of Data Acquisition/Setting
GHI, PAR and temperature collected from two measurement stations are used for the work. The coordinates and type of climate corresponding to both stations are shown in Table 1.  The station located in Tabernas (Spain) has an arid Mediterranean climate (BShs in Köppen classification) with cool winters and very warm summers. Temperatures range from −5 • C to 45 • C. Rainfall is very low, accumulating an average of only 243 mm.
On the other hand, Lugo is a city in north-western Spain. It has a humid oceanic climate with dry summers, Cfb in the Köppen climate classification. Due to its remoteness from the Atlantic, its annual precipitation of 1084 mm can be considered low compared to the near areas.

Methodology
Below the steps followed to carry out this work are detailed: (a) PAR and temperature measurements are collected. These measurements are used for the GEOPAR project (Project CGL2016-79284-P AEI/FEDER/UE) developed by Centro de Investigaciones Energéticas, Medioambientales y Tecnológicas (CIEMAT). Lugo station includes an Apogee SQ-110 pyranometer for PAR measuring, as well as a HOBO S-THB Temperature/Humidity M002 psychrometer that measure humidity and temperature. PAR and temperature measurements in Tabernas are obtained by an Eko ML-020P device, and by a Vaisala HMP60 psychrometer to measure humidity and temperature. (b) Algal and PV productivities are estimated from PAR and temperature data in the study stations. These estimations are carried out for the two systems studied (HRAP and AS + PV), considering the results both with storage and without storage. In turn, the PV facility is performed for different land requirements. (c) Probabilities of not achieving the energy consumption threshold (0.89 kWh m −3 ) one day in the AS system, according to the areas occupied by PV panels, are determined for each station and for different PV installation areas. (d) Probabilities of not achieving the corresponding HRAP energy consumption threshold (0.06 kWh/m 3 ) are also determined for each station. (e) A PV facility viability assessment by exceedance probabilities is also included. For it, P50, P90 and P99 are estimated according to the PV installation areas for both cases, with and without storage. (f) Finally, an opportunity cost analysis is included. For this analysis, the costs associated to both studied systems (HRAP and AS + PV) are compared. For HRAP case, two costs are considered: the capital cost and the terrain cost. For AS + PV, besides the capital cost of AS, the capital and operation cost of PV, the terrain cost and, in the case of storage, the economic storage pack cost, are considered. Other possible opportunity costs are pointed out.

Calculations/Performance Models
The algal productivity depends on two meteorological variables: photo-synthetically active radiation (PAR) and temperature. In the case of the PV system, its productivity depends on the Global Horizontal Irradiance (GHI). Two scenarios described in [46] to perform two productivity models are used, one for the HRAP system and the other for the AS system.

Algal Productivity Model
A revision of models developed to estimate the algal productivity can be seen in [52]. In this work, the model performed in [53] is used, according to which the productivity Energies 2021, 14, 1994 5 of 18 P bio_est can be estimated from the difference between the specific growth rate, G, and the specific respiration rate, R. In turn, these rates can be assessed as follows [54].
The parameters of these expressions are shown in Table 2.
being ∅ T the temperature-dependent function: Table 3. Experimental parameters.
Besides, the HRAP system described in [51] is studied. The scenario 1 of this work considers a system with a total installation area of 40,000 m 2 and a flow rate of 1950 m 3 d −1 , whose total energy consumption is 0.06 kWh m −3 .

PV Productivity Model
PV power can be estimated from GHI and temperature using the following expression [56,57]: where PV est is the power generation estimation, PV ins is the PV installation capacity, PCS loss are the power losses due to the power conditioning system, System loss are the losses associated to PV system, G s is the GHI at standard test conditions and T loss is a reduction parameter related to the PV module temperature. This last parameter can be estimated as follows: where α pmax is the temperature dependency of PV power generation, T air is the atmospheric temperature and ∆T is the difference in PV module temperature.
The efficiency of the system as a function of the operating time can also decrease due to a decrease in the energy efficiency of photovoltaic panels and in the efficiency of energy storage (battery consumption). However, being the study conducted for one year time, the losses of efficiency of the PV system with the operation time can been considered negligible. If the study were for more years, these losses should be included-the average annual decrease in electricity yield could be near 0.60%·yr −1 [58].
The values taken from [57] for the different parameters are shown in Table 4. Table 4. Estimation of parameters for PV productivity model.

Parameter Units Estimation
The model above described has been performed again from a scenario drawn in [51] (scenario 3 of this work, corresponding to an AS system). The flow rate is the same than in the HRAP system (1950 m 3 d −1 ), but the total installation area is much lesser than in the algal pond system (only 900 m 2 vs. 40,000 m 2 ). The total energy consumption in this system is 0.89 kWh m −3 (much higher than that of the first system). It needs to be stressed that the daily energy demand of wastewater treatment is determined by an uneven wastewater inflow and pollutant load. The weather conditions and the wastewater temperature also affect the efficiency of the wastewater aeration system (solubility of oxygen in water). So, the amount of energy needed to remove a unit of load of organic pollutants from wastewater could have been calculated in kWh per kg BOD (biological oxygen demand) that is removed from wastewater. However, in the scenarios described in [51], energy consumption is considered constant. Actually, as shown in [59], in recent biological processes, technology has lowered energy requirements, making energy consumption much more stable.

Installation Area Analysis in Absence of Storage
The intra-annual variability of the HRAP system productivity, along with that of the meteorological variables (PAR and temperature) affecting this productivity are shown in Figure 1. This productivity has been estimated from the model and scenario previously described (40, The seasonal variability of PAR and temperature can also be observed in the power curves. Regarding PAR, Lugo station shows more fluctuations than Tabernas station, due to the oceanic climate of the former. In this climate, precisely, the seasonal variation is less pronounced. Thus, as mentioned, the power curves under outdoor conditions are clearly affected by these climatological aspects. On the other hand, the difference between the areas occupied by both systems (HRAP and AS system) is very large (40,000 m 2 vs. 900 m 2 ). The energy consumption of the HRAP and of the AS system is also very different (0.06 kWh m −3 vs. 0.89 kWh m −3 ), so the power need is much lower in the case of the algal system. In fact, the power derived from HRAP The seasonal variability of PAR and temperature can also be observed in the power curves. Regarding PAR, Lugo station shows more fluctuations than Tabernas station, due to the oceanic climate of the former. In this climate, precisely, the seasonal variation is less pronounced. Thus, as mentioned, the power curves under outdoor conditions are clearly affected by these climatological aspects.
On the other hand, the difference between the areas occupied by both systems (HRAP and AS system) is very large (40,000 m 2 vs. 900 m 2 ). The energy consumption of the HRAP and of the AS system is also very different (0.06 kWh m −3 vs. 0.89 kWh m −3 ), so the power need is much lower in the case of the algal system. In fact, the power derived from HRAP Energies 2021, 14, 1994 8 of 18 covers this consumption most of the year. However, covering the 900 m 2 of the AS system with PV panels of 2 m 2 with an installation capacity of 400 W (450 solar panels), the PV power never reaches this threshold.
Thus, in order to achieve self-sufficiency, the PV power must be increased and thus, the area occupied by panels. In Figure 2, the PV power obtained for different areas is shown.
Energies 2021, 14, x FOR PEER REVIEW 9 of 19 covers this consumption most of the year. However, covering the 900 m 2 of the AS system with PV panels of 2 m 2 with an installation capacity of 400 W (450 solar panels), the PV power never reaches this threshold. Thus, in order to achieve self-sufficiency, the PV power must be increased and thus, the area occupied by panels. In Figure 2, the PV power obtained for different areas is shown.  As areas occupied by the PV facility increases, the number of days in which the energy consumption threshold is overcome also increases. The same as the case of power derived from biomass, the seasonal evolution, associated to climate characteristic, is evident. This evolution is more pronounced in Tabernas and has more daily fluctuations in Lugo.
In order to compare the PV powers ( Figure 2) with the energy consumption threshold (0.89 kWh m −3 ), the probability of not achieving that threshold at least one day is estimated, according to the areas occupied by PV panels (Figure 3). This probability can serve to assess the viability of the project when a PV facility supplies the energy needed by the AS system. As observed, the decreasing slope is very pronounced for smaller areas, but later this slope strongly diminishes. This change in the slope is more abrupt in the case of Tabernas. In Figure 3, the probabilities of not achieving the corresponding HRAP energy consumption threshold (0.06 kWh m −3 ) have been included. These probabilities are 0.0546 for Tabernas case and 0.2077 for Lugo. In this situation, the behaviour of both systems could be considered similar, as the probabilities of not achieving the corresponding thresholds, according to the system considered (HRAP or AS system with PV power), are the same. These intersection points are obtained for 6030 m 2 in Tabernas and for 6552 m 2 in Lugo. Thus, these areas of PV panels are needed to obtain the same probability than in the 40,000 m 2 HRAP system. As areas occupied by the PV facility increases, the number of days in which the energy consumption threshold is overcome also increases. The same as the case of power derived from biomass, the seasonal evolution, associated to climate characteristic, is evident. This evolution is more pronounced in Tabernas and has more daily fluctuations in Lugo.
In order to compare the PV powers ( Figure 2) with the energy consumption threshold (0.89 kWh m −3 ), the probability of not achieving that threshold at least one day is estimated, according to the areas occupied by PV panels (Figure 3). This probability can serve to assess the viability of the project when a PV facility supplies the energy needed by the AS system. As observed, the decreasing slope is very pronounced for smaller areas, but later this slope strongly diminishes. This change in the slope is more abrupt in the case of Tabernas. In Figure 3, the probabilities of not achieving the corresponding HRAP energy consumption threshold (0.06 kWh m −3 ) have been included. These probabilities are 0.0546 for Tabernas case and 0.2077 for Lugo. In this situation, the behaviour of both systems could be considered similar, as the probabilities of not achieving the corresponding thresholds, according to the system considered (HRAP or AS system with PV power), are the same. These intersection points are obtained for 6,030 m 2 in Tabernas and for 6,552 m 2 in Lugo. Thus, these areas of PV panels are needed to obtain the same probability than in the 40,000 m 2 HRAP system.

Installation Area Analysis Including Storage
The former study does not consider the possibility of storing the energy surplus during those days in which the corresponding energy consumption threshold is exceeded. However, this possibility must be taken into account, as it is very difficult to reach the consumption threshold every day. Thus, in this section the self-sufficiency in both analysed systems is studied assuming that the possible energy surplus of one day is available for subsequent days.

Installation Area Analysis Including Storage
The former study does not consider the possibility of storing the energy surplus during those days in which the corresponding energy consumption threshold is exceeded. However, this possibility must be taken into account, as it is very difficult to reach the consumption threshold every day. Thus, in this section the self-sufficiency in both analysed systems is studied assuming that the possible energy surplus of one day is available for subsequent days.

Results for HRAP System
In the case of the HRAP system, the estimated power from biomass shown in Figure 1e,f can be modified by storing the energy that overcomes 0.06 kWh m −3 and by supplying this energy those days in which the power is lower than this threshold. By this storage system, the power intra-annual evolution no longer shows values below 0.06 kWh m −3 (Figure 4).

Results for HRAP System
In the case of the HRAP system, the estimated power from biomass shown in Figure  1 (e) and f)) can be modified by storing the energy that overcomes 0.06 kWh m −3 and by supplying this energy those days in which the power is lower than this threshold. By this storage system, the power intra-annual evolution no longer shows values below 0.06 kWh m −3 (Figure 4).

Results for AS System with PV Panels
Similar to the variability without storage represented in Figure 2, the PV power intraannual variability in case of storage for different PV installation areas is shown in Figure  5. The graphs show that, as area increases, the possibility that the threshold is achieved all days also increases. Indeed, for the area represented in Figure 5, this achievement is accomplished from 4000 m 2 both in Tabernas and in Lugo. The daily fluctuations are again much more pronounced in Lugo than in Tabernas.

Results for AS System with PV Panels
Similar to the variability without storage represented in Figure 2, the PV power intraannual variability in case of storage for different PV installation areas is shown in Figure 5. The graphs show that, as area increases, the possibility that the threshold is achieved all days also increases. Indeed, for the area represented in Figure 5, this achievement is accomplished from 4000 m 2 both in Tabernas and in Lugo. The daily fluctuations are again much more pronounced in Lugo than in Tabernas.
Similarly to the case without storage (Figure 3), storage may be included and the probability of not achieving the corresponding threshold at least one day according to the area occupied by PV panels may be obtained ( Figure 6). This probability can serve again to assess the viability of the AS system based on PV facility. In this case, the null probability is quickly reached (for 2500 m 2 in Tabernas and for 3500 m 2 in Lugo) as shown in Figure 6, unlike the case without storage where this null probability is never reached, as represented in Figure 3. That is, the storage allows us to ensure that the demand is always satisfied by using much less PV area, making the project more viable. On the other hand, the probabilities of not achieving the corresponding threshold for one day in the HRAP system, when storage is used, are 0 in both stations, as all values overcome this threshold (0.06 kWh m −3 ) (Figure 4). Thus, 2500 m 2 in Tabernas and 3500 m 2 in Lugo are, precisely, the areas needed to reach the same probability than in case of 40,000 m 2 HRAP system with storage. Therefore, thanks to storage, the land requirement has been greatly reduced in both locations.  Similarly to the case without storage (Figure 3), storage may be included and the probability of not achieving the corresponding threshold at least one day according to the area occupied by PV panels may be obtained ( Figure 6). This probability can serve again to assess the viability of the AS system based on PV facility. In this case, the null probability is quickly reached (for 2500 m 2 in Tabernas and for 3500 m 2 in Lugo) as shown in Figure  6, unlike the case without storage where this null probability is never reached, as represented in Figure 3. That is, the storage allows us to ensure that the demand is always satisfied by using much less PV area, making the project more viable. On the other hand, the probabilities of not achieving the corresponding threshold for one day in the HRAP system, when storage is used, are 0 in both stations, as all values overcome this threshold (0.06 kWh m −3 ) (Figure 4). Thus, 2500 m 2 in Tabernas and 3500 m 2 in Lugo are, precisely, the areas needed to reach the same probability than in case of 40,000 m 2 HRAP system with storage. Therefore, thanks to storage, the land requirement has been greatly reduced in both locations.

Comparison between the Storage and the Absence of Storage Cases
The project viability based on PV power has also been assessed by three recommended exceedance probabilities (P50, P90 and P99). In Figure 7, exceedance probabilities vs. the PV installation area have been represented for both cases, with and without storage. Since the PV power increases with the area used, the exceedance probabilities also increase as this area increases.

Comparison between the Storage and the Absence of Storage Cases
The project viability based on PV power has also been assessed by three recommended exceedance probabilities (P50, P90 and P99). In Figure 7, exceedance probabilities vs. the PV installation area have been represented for both cases, with and without storage. Since the PV power increases with the area used, the exceedance probabilities also increase as this area increases.
As observed, the slopes are higher in Tabernas than in Lugo, as the power achieved by PV is clearly higher in Tabernas, with a Mediterranean climate, than in Lugo, and this is more evident as the PV installation area increases. Storage brings probabilities higher or equal than those obtained in absence of storage, as power is increased by storage. Besides, storage affects small PV power, so the lower percentiles increase with storage. Thus, there is a greater difference between both cases (storage and absence of storage) for higher exceedance probabilities, associated to lower percentiles. That is, the differences between P99s are higher than the differences between P50s.
These probabilities provide information about the system, but they do not allow us to qualify a project as viable or not viable, unlike the study based on the project self-sufficiency. To achieve self-sufficiency, the estimate of the size of a power facility is of great interest. Indeed, the PV facility installation area required to satisfy the need of energy consumption in the AS system has been above estimated (Figure 3 in absence of storage and Figure 6 with storage). According to the results, storage leads to a large diminishing in this area: in absence of storage, the probability of not achieving the energy consumption threshold is always higher than zero in both stations, while this probability is zero for 2500 m 2 at Tabernas and for 3500 m 2 at Lugo, when storage is considered. As observed, the slopes are higher in Tabernas than in Lugo, as the power achieved by PV is clearly higher in Tabernas, with a Mediterranean climate, than in Lugo, and this is more evident as the PV installation area increases. Storage brings probabilities higher or equal than those obtained in absence of storage, as power is increased by storage. Besides, storage affects small PV power, so the lower percentiles increase with storage. Thus, there is a greater difference between both cases (storage and absence of storage) for higher exceedance probabilities, associated to lower percentiles. That is, the differences between P99s are higher than the differences between P50s.
These probabilities provide information about the system, but they do not allow us to qualify a project as viable or not viable, unlike the study based on the project self-sufficiency. To achieve self-sufficiency, the estimate of the size of a power facility is of great interest. Indeed, the PV facility installation area required to satisfy the need of energy consumption in the AS system has been above estimated (Figure 3 in absence of storage and Figure 6 with storage). According to the results, storage leads to a large diminishing in this area: in absence of storage, the probability of not achieving the energy consumption

Opportunity Cost Analysis
In the previous sections, space needs for self-sufficiency of two projects were assessed. However, the European Union Green Deal needs also to assess the opportunity cost of different energy sources, in special in terms of the availability and possible alternative uses of land or capital. Some green energy sources have great inconveniencies in terms of the needs of land and the loss of its alternative uses. Besides, the wastewater treatment systems proposed in the work need specific orographic conditions, such as non-steep slopes and they must not be installed in protected areas due to environmental reasons. Finally, other terrain uses must be considered, such as the profitability of a possible agriculture or livestock exploitation, which also implies benefits for farmers drawn from the diversification of income. In this sense, the fact that the optimal social and private solutions may differ must be taken into account.
The EU Green Deal aims at reaching the target of climate neutrality. Then, the energy system must be decarbonised, energy efficiency prioritised and a power sector based largely on renewable resources developed. But this objective will only be achieved by an opportunity cost analysis between alternative renewable energies. For instance, the extensive need for terrain in the HRAP projects limits their viability and convenience as compared to AS+PV-the need for land is 40,000 m 2 versus less than 6500 m 2 . However, the prices of both systems also need to be compared. A von Thünen framework [60] is suitable for analysing economic decisions when distance from different uses of energy matters. Transport is costly and reduces the profitability of bioenergy production when the distance to the power plant increases. In this case, distance is zero, but the result obtained with the alternative supplies of energy from other possible locations could be compared, as well as the possible uses of the energy produced in these locations in other facilities or industries. Some studies compare the process which uses the combustion of biomass with that of bio crude extraction [61]. Besides, as previously mentioned, different possible uses of the installation areas could be assessed [62,63]. These price comparisons are very specific of each location and moment and, so, they need to be done by each project.
Being aware of the circumstantial character of this price analysis, an economic analysis comparing the costs associated to both analysed systems, HRAP and AS system with PV facility, has been undertaken. For HRAP case, two type of costs are considered: the capital cost and the terrain cost. The operation and maintenance cost (energy and flocculants consumption) is not included as the energy consumption is supplied by the algae itself. According to [48], the capital cost is 192.55 €/p.e. Since the system was designed to serve a population equivalent to 10,000 p.e. the total capital cost will be 1,925,500 €. The terrain cost is also different according to the site considered (71.9 €/m 2 at Lugo and 150.8 €/m 2 at Tabernas according to the Spanish Official Statistics, http://www.fomento.es/be2/, accessed on 1 September 2020). As regards to the AS system with PV facility, according to [24], the AS construction costs are typically 70% higher than those of HRAPs, so the capital costs may be established in 6,418,333 €. On the other hand, the capital and operation cost of PV is 350 €/m 2 (an average estimation of the one offered in [64]). Additionally, the terrain cost must be added. Finally, the economic storage pack cost is included, that according to [65] can be estimated in 200 €/kWh of the battery energy storage systems. For the case of Lugo, the need of storage for 3500 m 2 is of 94.7 kWh, while in Tabernas for 2500 m 2 this need is of 81.1 kWh. The costs obtained for both systems are shown in Table 5: As shown, the total cost in the HRAP system is clearly lesser than in the AS + PV system in absence of storage. The difference is mainly due to the important capital cost of AS, although the difference in the terrain cost is, especially in Tabernas, very high, something that affects more the cost of the HRAP system. On the other hand, the presence of storage allows for a significant reduction of the costs of the PV facility. For example, in the case of Tabernas, although the HRAP system is 15.7% less costly than the AS + PV system, storage makes it more costly in a 3.5%. Indeed, storage provides a great advantage for using the AS + PV in locations where a lot of irradiance is available, as is the case with Tabernas. The existence of a threshold of irradiance may be needed so as to improve the costs of the AS + PV as compared to HRAP.
Then, the previous considerations about disadvantages linked to large land requirements of HRAP system seem to be balanced out by the little requirement of capital linked to HRAP. As shown, this opportunity cost analysis is then essential for a better assessment of different energy sources.

Conclusions
The present study concludes that the energy consumption is much smaller for HRAPs than for AS systems, where PV power is needed to supplement energy. However, HRAPs has much larger land requirements. Besides, storage reduces significantly the cost of the PV facility, providing a great advantage for the use of AS + PV in locations where a lot of irradiance is available. This study shows that climate characteristics are crucial for projects viability (the evolution of PV power according to seasonal evolution is more pronounced in Tabernas and has more daily fluctuations in Lugo). However, land prices should be compared, as well as the distance from different uses of energy. Here, a comparison of the total cost of the HRAP system with that of AS + PV has been made, showing that the latter is higher due to the important capital costs associated to AS, in spite of the fact that land requirements lead to a significant increase in the cost of the HRAP system. Another insight of the article is that PV power exceedance probabilities increase with the PV installation area, and that the probability differences between situations with storage and without storage increase for higher exceedance probabilities. These probabilities do not allow to qualify a project as viable or not, unlike the study based on the project self-sufficiency. In short, the opportunity cost analysis is essential for a better assessment and understanding of different energy sources. For future works, other renewable sources, such as the small wind energy, will be studied.