Exploring Wind and Solar PV Generation Complementarity to Meet Electricity Demand

: Understanding the spatiotemporal complementarity of wind and solar power generation and their combined capability to meet the demand of electricity is a crucial step towards increasing their share in power systems without neglecting neither the security of supply nor the overall cost e ﬃ ciency of the power system operation. This work proposes a methodology to exploit the complementarity of the wind and solar primary resources and electricity demand in planning the expansion of electric power systems. Scenarios that exploit the strategic combined deployment of wind and solar power against scenarios based only on the development of an individual renewable power source are compared and analysed. For each scenario of the power system development, the characterization of the additional power capacity, typical daily proﬁle, extreme values, and energy deﬁcit are assessed. The method is applied to a Portuguese case study and results show that coupled scenarios based on the strategic combined development of wind and solar generation provide a more sustainable way to increase the share of variable renewables into the power system (up to 68% for an annual energy exceedance of 10% for the renewable generation) when compared to scenarios based on an individual renewable power source. Combined development also enables to reduce the overall variability and extreme values of a power system net load.


Introduction
Renewable energy sources (RES) continue to grow and gain increased relevance in modern electric power. The main driver of this growth was based on subsidies, typically, and feed-in tariffs that aim to reduce the air pollution through the replacement of fossil energy sources by clean and safe RES [1][2][3]. Within the different types of RES, wind and solar photovoltaic (PV) are the most promising technologies to jointly produce a large share of renewable electricity needed to reach the EU's 2030 ambitious targets [4,5]. Indeed, most of the power system scenarios to achieve 100% RES in European countries strongly rely on these two power sources that differ from conventional energy sources, due to the stochastic nature of their primary resource and the non (or poorly) dispatchability of the generated power, being normally referred to as variable renewable energies (VRE).
The stochastic behavior of VRE introduces new challenges to a transmission system operator (TSO) that can exacerbate the power system's flexibility requirements to deal with their variability, especially in periods with severe power ramp events [6,7]. The added flexibility required to maintain a stable operation will increase the overall operational costs of the power system and may reduce the environmental benefits of VREs. Therefore, to decarbonize the electric power systems in a sustainable manner, new capacity deployment planning strategies addressing the variability and uncertainty associated with the large-scale integration of VRE are needed [8][9][10]. One of these strategies can methodology permits to optimize the wind and solar PV power capacity ratio within a TSO control region (in this case a whole country) by considering the electricity demand needs and different optimization perspectives that enable to strategically identify the most adequate regions and install each power source. Using high resolution data and Portugal as a case study, the following deployment scenarios are considered and compared, strategically adding to the existing VRE energy mix: (1) only wind power, (2) only solar PV power, (3) adding wind and solar PV power. For each scenario, a characterization of the additional power capacity, typical daily profiles, extreme values, and the VRE energy surplus/deficit are assessed. The introduction of the VRE complementarity [28] and demand criteria in new VRE deployment constitutes a crucial step towards the planning of a secure and reliable nearly 100% renewable power system by reducing the flexibility needs to deal with the variability and uncertainty related to these technologies, and additional potential benefits are expected for the power system and the final consumers. This work does not replace the necessity of additional technologies capable of providing the power system's flexibility needs to deal with VRE variability, namely, the use of energy storage technology, which will always be necessary, and whose optimized dimensioning is outside the scope of this paper. Nevertheless, this work contributes significantly to the minimization of storage capacity and, as such, to the reduction of system costs, as this type of solution still has a high investment cost [1,29].
The methodology presented in this work represents a new planning paradigm, since until now, the only criteria to deploy the VRE capacity has been to find locations without any environmental/technical constraints that maximize the generation, which understandably allows higher revenues for the power producers although it can raise some concerns regarding the power systems' flexibility.

Holistic Renewable and Complementarity Assessment Methodology
A schematic flowchart with the main steps of the methodology developed is depicted in Figure 1. It should be noted that the proposed methodology only addresses the temporal and spatial integration aspects of wind and solar PV within a region of control (of the power system) and the benefits of the complementarity assessment. It does not cover other relevant aspects such as economic costs, land availability, social acceptance, storage technologies, and so forth that should be also analyzed in future works [30]. Below, a detailed explanation of the methodology is provided that is able to be transferred and applied in other regions/power systems.

Identification of Power Generation Patterns
The identification of power generation patterns based on high-resolution data allows synthesizing the information of each point of the grid by extracting the main pattern features in the region under analysis. With this approach, it is possible to filter the local effects of the individual grid points and group grid points with a similar power production behavior providing a better understanding of generation features [31,32]. Moreover, instead of detecting specific grid points that could be unfeasible due to constraints (e.g., protected zones), using this approach, the most suitable It should be noted that the proposed methodology only addresses the temporal and spatial integration aspects of wind and solar PV within a region of control (of the power system) and the benefits of the complementarity assessment. It does not cover other relevant aspects such as economic costs, land availability, social acceptance, storage technologies, and so forth that should be also analyzed in future works [30]. Below, a detailed explanation of the methodology is provided that is able to be transferred and applied in other regions/power systems.

Identification of Power Generation Patterns
The identification of power generation patterns based on high-resolution data allows synthesizing the information of each point of the grid by extracting the main pattern features in the region under analysis. With this approach, it is possible to filter the local effects of the individual grid points and group grid points with a similar power production behavior providing a better understanding of generation features [31,32]. Moreover, instead of detecting specific grid points that could be unfeasible due to constraints (e.g., protected zones), using this approach, the most suitable regions for the installation of future VRE capacity and plants are identified. Thus, the identification of power production patterns that distinguish regional circulations is a key step to explore the wind and solar PV power production complementarity. This step can be obtained through a clustering analysis.
Cluster analysis is a classic technique employed to divide data among a set of clusters according to the similarity among observations in such a way that cluster elements are always more similar (in some sense) to each other than to those of other groups, with higher dissimilitude between the clusters [33,34]. With this step, it is possible to detect statistical patterns of data that can often be associated with different physical processes. Although cluster analysis does not always reveal the results of a specific physical process, if there are prevailing behaviors in a group of elements, they will be reflected in one of the clusters. Therefore, in this work, the power production patterns were identified using the K-medoids clustering technique [34], Equation (1), applied to the power production input matrix (X n ), Equation (2).
In the previous equations, D is the distance between the observations and the medoid of each cluster, C k is the cluster medoid of the K-th cluster, n represents the records available, while i represents the spatial points available. The clustering algorithm is applied separately for wind and solar PV production spatial fields. Despite the subjectivity on the exact definition of the cluster centers, this methodology enables a consistent classification for most of the data. The K-medoids technique is a non-hierarchical clustering algorithm to group the data in K clusters, where K is beforehand defined. The optimal value of K is determined using the Calinski-Harabasz (CH) criterion [35]. This criterion allows detecting the appropriate number of clusters by calculating the Euclidean distance between the clusters and comparing the internal sum of squared errors of each cluster. Hereafter, the production patterns obtained using this technique are designated as wind production patterns (WPP) and solar PV production patterns (SPP).

Optimization Greedy Algorithm
Currently, due to the increasing environmental, energy security, and policy goals, the deployment of renewable capacity is constrained by different objectives functions. In multi-objective problems, there is no solution that simultaneously optimizes all the being necessary to decide a priori the relative importance (the weight) of each objective to create a suitable (prescribed) optimal objective function [36]. Thus, despite the importance, for example, of the storage systems and interconnections to accommodate energy surpluses/deficits, this work focus only in the synergy between wind and solar generation to attain the electricity demand. In this sense, a greedy algorithm to identify the future VRE deployment pathways was implemented [37]. This algorithm is especially useful for optimization problems with linear objective functions, due to its simplicity when compared with other heuristic search algorithms that present a higher computational effort and require the definition of path cost functions.
The greedy algorithm chooses the "most attractive" alternative in every iteration [37]. Broadly speaking, using the current net load power capacity (computed on the basis of the demand and the current VRE generation) as a baseline, at each iteration, 10 MW are added to the previous net load time series by selecting the WPP or SPP that minimizes a predetermined criterion. The algorithm stops when Energies 2020, 13, 4132 5 of 21 a determined annual VRE energy generation exceeds the annual electricity consumption-Surplus Max . Different VRE energy exceedance (i.e., surplus) values were tested-from 1% to 10%. The maximum value criterion-10%-is related to the maximum VRE energy curtailment observed during an entire year [38]. An additional restriction for the installation of new parks in each region was imposed by considering only a maximum of 10% (Max area ) of the area available for additional wind and solar power capacity deployment (WPP area and SPP area ), i.e., after the identification of the spatial points associated to each cluster, only 10% of the area available can be used. For the wind turbine used in this work, a wind park with 10 MW and an interturbine distance of three rotor wind turbine diameter requires an area of 1.0 km 2 (WindPark area ), while the area needed (based on data from Portugal) was a typical fixed solar PV park with 10 MW that corresponds to 0.3 km 2 (SolarPark area ).
Two different objective functions are dedicated to identifying the pathways regarding the integration of these technologies into power systems. These functions are the minimization of: (i) the net load (NL) annual variability by computing the annual standard deviation (Equation (3)) and (ii) one-hour net load step-change standard deviation (Equation (4)). The selection of these objective functions is related to the operational impact of VRE variability into the power system, namely, the required operational capacity that TSO has to adjust, using the production of dispatchable power plants (running or start-up/shut-down), storage, and interconnections, for maintaining the power system balance [8,39].
The proposed model has two instances according to the definition of Z(t) that can either be (5) or (6): and are subject to restrictions related to the annual VRE energy generation surplus regarding the annual electricity consumption (7) and the area available for additional wind (8) and solar (9) power capacity deployment within each region: w it N WPP (it, w) · WindPark area ≤ w Max area · WPP area (w) (8) s it N SPP (it, s) · SolarPark area ≤ s Max area · SPP area (s) (9) with: Z, NL, VRE, VRE surplus , N W , N S ∈ N (13) In Equations (7) to (13), N WPP and N SPP represent the number of wind and solar power parks identified through the optimization algorithm for each WPP and SPP, VRE surplus denotes the power generation surplus. Demand represents the hourly electricity consumption while WP and SP denote the existing wind and solar PV hourly generation.

Characteristics of the Portuguese Power System and Weather Conditions
In Figure 2a, the geographical wind and solar installed power capacity are depicted for 2015 (reference year in this work, as in the Portuguese national energy and climate plan for 2030). At the end of 2015, the wind power capacity was 5034 MW. At the end of 2019, the wind power installed capacity only increased 391 MW to 5429 MW. Contrasting with the solar PV power capacity, most of the wind power sites are located in the mountain regions of the center/north of Portugal, with some exceptions on the west and southwest coast. This concentration is explained by the wind resource available, and as described in [40], it can introduce additional challenges to TSOs. Indeed, the authors relate the occurrence of strong wind power ramp events, unleashed by certain weather conditions that usually extend over several hundreds of kilometers (e.g., cold fronts), with these locations. Moreover, as depicted in Figure 3, due to local effects such as the mountain/valley breezes, on average, the highest wind power generation levels occur during nighttime, showing no correlation with the electricity demand. Thus, and taking into account the wind energy sector maturity in Portugal, a strategical wind power deployment is needed to (i) complement the current wind power production promoting the statistical power smoothing effect and (ii) meet the electricity demand.  In Figure 3 is also depicted the demand and the net load profile (the consumption minus the wind and solar production). The net load corresponds to the load that needs to be suppressed by the remaining power generation fleet. On average, the highest demand consumption values are observed between 7-9 p.m., due to domestic users. To some extent, the current net load profile allows already to observe the complementarity between wind and solar PV generation, which contributes to a daily average value of nearly 1800 MW. Nevertheless, this complementarity is incapable of following the electricity demand. The annual VRE energy surplus percentage value regarding the electricity demand (i.e., the instantaneous sum of the wind and solar generation above the consumption) during the year 2015 was only 9.2e−04.  remaining power generation fleet. On average, the highest demand consumption values are observed between 7-9 p.m., due to domestic users. To some extent, the current net load profile allows already to observe the complementarity between wind and solar PV generation, which contributes to a daily average value of nearly 1800 MW. Nevertheless, this complementarity is incapable of following the electricity demand. The annual VRE energy surplus percentage value regarding the electricity demand (i.e., the instantaneous sum of the wind and solar generation above the consumption) during the year 2015 was only 9.2e−04. As identified for other phenomena (e.g., precipitation [42]), due to the geographical location, along with its topographical features and weather conditions, different regional wind and solar power production patterns can be observed in Portugal [43]. The atmospheric circulation in Portugal is strongly influenced by the seasonal migration of the mid-latitude high pressure system [44]. In winter, when the subtropical high pressure is usually centered at lower latitudes, warm and cold fronts and other baroclinic synoptic perturbations moving eastward from the Atlantic Ocean are observed. In late spring to early autumn, the western perturbations are blocked by this subtropical system that expands to higher latitudes (typically centered in the Azores). Under these conditions, The solar PV power capacity is still low in Portugal (731 MW at the end of 2019). At the end of 2015 (reference year), the total capacity was 454 MW. Most of this capacity is located in the south region of the country, Figure 2a. Nevertheless, with the rapid growth of the solar PV industry and the substantial drop in the price of solar panels, it is foreseen that solar PV generation will play a crucial role in the Portuguese power system in the near future. This expectative, strongly supported by recent governmental policies, is associated with: (i) a very high resource potential, due to the geographical position and climate conditions of Portugal, (ii) the opportunity to decentralize electricity production enabling to reduce transmission losses, and (iii) the daily profile of the solar resource patterns that are correlated with the electricity demand ( Figure 3).
In the last years, the annual share of electricity generation from renewable energy sources has ranged from nearly 40% to 60%, Figure 2b. The wind power generation contributed with 24% in 2015 and 28 % in 2019. These values are even more significant, as the electricity demand also increased from 49.0 TWh in 2015 to 50.3 TWh [41]. The contribution share from solar PV power is still residual −2%. From Figure 2b it is possible to verify that the dependence on fossil fuels to balance the demand is essentially supplied by coal and natural gas. The reduced hydro contribution is associated with extremely dry years. Indeed, the role of this technology in the Portuguese power system has been characterized by a strong variability in consecutive years (severe dry year followed by strong raining year). As in many countries, hydropower reservoirs (with and without pumped technology) in Portugal are also responsible for counterbalancing the electricity demand and the VREs variability. However, due to the increased environmental restrictions, in the coming years, a significant deployment capacity of this type of technology in Portugal is not foreseen. Thus, these limitations associated, coupled with the foreseen decommissions of the old coal-based power plants, highlight the need for a careful assessment of the deployment of the new VRE power capacity.
In Figure 3 is also depicted the demand and the net load profile (the consumption minus the wind and solar production). The net load corresponds to the load that needs to be suppressed by the remaining power generation fleet. On average, the highest demand consumption values are observed between 7-9 p.m., due to domestic users. To some extent, the current net load profile allows already to observe the complementarity between wind and solar PV generation, which contributes to a daily average value of nearly 1800 MW. Nevertheless, this complementarity is incapable of following the electricity demand. The annual VRE energy surplus percentage value regarding the electricity demand (i.e., the instantaneous sum of the wind and solar generation above the consumption) during the year 2015 was only 9.2e−04.
As identified for other phenomena (e.g., precipitation [42]), due to the geographical location, along with its topographical features and weather conditions, different regional wind and solar power production patterns can be observed in Portugal [43]. The atmospheric circulation in Portugal is strongly influenced by the seasonal migration of the mid-latitude high pressure system [44]. In winter, Energies 2020, 13, 4132 8 of 21 when the subtropical high pressure is usually centered at lower latitudes, warm and cold fronts and other baroclinic synoptic perturbations moving eastward from the Atlantic Ocean are observed. In late spring to early autumn, the western perturbations are blocked by this subtropical system that expands to higher latitudes (typically centered in the Azores). Under these conditions, this synoptic circulation is weak and diffused over Portugal, which promotes small scale (e.g., valley/mountain and sea/land breezes) and regional locally-forced atmospheric features as a thermal low pressure system inland in the middle of the Iberian Peninsula [45]. During this period, when a sharp contrast between high temperatures over land and inferior temperatures over the sea is observed, strong parallel winds and changes in the cloudiness (and radiation) [43] are triggered. These seasonal winds, also known as Nortada, are observed mainly in western coastal regions [46]. The inner center/south and southern regions show a complex topography, with several mountains. Thus, the displacements and changes in the intensity of the subtropical high pressure, along with their interaction with the complex orography, generate a range of differentiated regional wind and solar climates [43].

Data
For the purpose of this study, one-hour resolution data (electricity demand, wind, and solar PV power production time series) for the year 2015 were used. As supported by [9,47], an hourly resolution permits modelling the energy system flexibility requirements with a sufficient level of detail.
The year 2015 was used as a reference period for the Portuguese national energy and climate plans and, according to the Portuguese system operator, was considered a typical year, regarding wind power production.
The national aggregated demand and production data were gathered from the national TSO website [41]. In this study, the electricity demand was considered inelastic and its potential increase in the next years with electricity consumption was not considered.
The wind speed data were derived from a numerical weather prediction (NWP) model, with a high spatial resolution (5 km × 5 km). NWP allows performing regional and/or national wind speed characterization without resorting to an extensive and costly network of anemometric stations. These models have the ability to describe, in an accurate way, the behavior and evolution of air masses and to treat explicitly the inherent phenomena of turbulence and atmospheric stratification [48]. In this work, the fifth-generation mesoscale model (also known as MM5) model [48] was used, and it was fed with the NCEP Climate Forecast System Version 2 (NCEP-CFSv2) reanalysis [49], with a horizontal grid spacing of 0.5 • × 0.5 • . On the basis of sensitivity tests of mesoscale model initial and boundary data and parameterizations and using meteorological stations the model was calibrated for the region under analysis. The model was configured for recording data every hour during 2015 in three domains. The calibration steps follow a similar approach to the one described in [50]. The last domain has a spatial resolution of 5 km, encompassing the onshore and offshore regions of Portugal. Regarding the offshore data, only the spatial points within the bathymetries 0 to 300 m were considered. For each point, wind speed data at 80 m above ground level were obtained and converted to wind power using a power curve of the most typical wind turbine in Portugal-Enercon E82, with 2.0 MW nominal capacity. According to the power curve used, the operational wind speeds range from 2 m/s (cut-in) to 25 m/s (cut-out), reaching the nominal power at 13 m/s. To avoid the installation of wind turbines in locations without economic viability, the spatial points that show a yearly wind speed below 5.5 m/s were excluded.
The hourly solar power data used in this work were gathered from the Photovoltaic Geographical Information System (PVGIS) web-based tool, which was developed at the Renewables & Energy Efficiency Unit. The solar PV power production is based on solar radiation data from the Satellite Application Facility on Climate Monitoring (CMSAF) [51,52]. This database is being extensively validated over Europe, showing a high level of accuracy and being a suitable tool for solar radiation and photovoltaic uses at any place over the globe [53,54]. In the present work, the data were extracted for the same geographical points considered in the wind resource step that considers a crystalline Energies 2020, 13, 4132 9 of 21 silicon panel and the following configuration: 10 MW of installed peak PV power, system overall losses equal to 10%, and the optimal azimuth and inclination angles for each point when considering the influence of shadows from the terrain.

Wind Power
In Figure 4a, the average power capacity profiles considering a wind park with 10 MW for each spatial point is shown. Based on the CH criterion, a suitable number of WPP is ten. Figure 4b depicts the spatial distribution of each WPP, while Figure 4c exhibits the current nominal power capacity located within each WPP. In Figure 5, the average daily and monthly wind power patterns are depicted to provide further insights regarding the features of each WPP. These profiles are computed using all points associated with each cluster.

Solar PV Power
In Figure 6a, the average power capacity considering a solar PV park with 10 MW for each spatial point is shown. Based on the CH criterion, a suitable number of SPP is six. Figure 6b depicts the spatial distribution of each SPP, while Figure 6c exhibits the current nominal power capacity located within each SPP. In Figure 7, the average daily and monthly solar PV power patterns are depicted to provide further insights regarding the features of each SPP.

Solar PV Power
In Figure 6a, the average power capacity considering a solar PV park with 10 MW for each spatial point is shown. Based on the CH criterion, a suitable number of SPP is six. Figure 6b depicts the spatial distribution of each SPP, while Figure 6c exhibits the current nominal power capacity located within each SPP. In Figure 7, the average daily and monthly solar PV power patterns are depicted to provide further insights regarding the features of each SPP.   contrast between land/sea regions combined with their exposure to severe meteorological events, explain the WPPs #1, #2, #5, #7, and #10 daily and annual profiles. In these regions, the temperature contrast between the ocean and the land increases in the summer allowing to observe higher wind power production values during these months when compared to winter months, Figure 5b. In the winter, due to their location on the western coast, where most of the windstorms reach Portugal [55], significant wind power productions can also be observed. A decrease in the level of production is expected during the spring and autumn. Some small differences can be found in the amplitude of the daily profile and peak hour of these WPPs. Thus, WPPs #2, #7, and #10 have their peak production in the middle of the afternoon, and the reaming patterns tend to reach the maximum in the late afternoon. These dissimilarities can be associated with local thermal effects as discussed previously. Figure 4c evidences that more than 3 GW of the current wind power capacity is located in only four onshore regions-WPPs #3, #4, #6, and #8, which are spatially concentrated in the center/north of Portugal. All these WPPs show a similar daily pattern associated with the atmospheric thermal stratification [56,57]. On average, the highest wind power values are observed during the night-time and a reduction in their intensity is expected during the day, leading to a "U" shape daily profile, Figure 5a. In the case of WPP #8, the minimum value is reached at nearly noon, while in the remaining WPPs, the minimum is obtained in the first hours of the morning. WPP #9 also exhibits a similar daily profile, although the hourly gradients tend to show a smooth behavior contrary to previous WPPs. The daily amplitude of the power production is very reduced for WPP #4 compared to the remaining WPPs. This result can be partially explained by the widely spatial distribution of this WPP that can mitigate severe variations. The analysis of the monthly profile highlights the similarities of these WPPs, showing a clear annual cycle, Figure 5b. The high wind season is during the winter months.

Solar PV Power
In Figure 6a, the average power capacity considering a solar PV park with 10 MW for each spatial point is shown. Based on the CH criterion, a suitable number of SPP is six. Figure 6b depicts the spatial distribution of each SPP, while Figure 6c exhibits the current nominal power capacity located within each SPP. In Figure 7, the average daily and monthly solar PV power patterns are depicted to provide further insights regarding the features of each SPP.
Energies 2020, 13, x FOR PEER REVIEW 11 of 21 and as shown in Figure 6a, the south/center regions of Portugal present more favorable conditions regarding the exploitation of this technology compared with the north region. By comparing the WPP ( Figure 5) and the SPP (Figure 7) profiles, a high temporal complementarity degree can be observed, namely, between WPPs #3, #4, #6, and #8 and all SPP profiles. In these cases, the wind production profiles show the lowest production period during the day, when solar production is higher, and vice versa.

Scenarios and Optimization Perspectives
The VRE deployment scenarios and optimization perspectives analyzed in this work are presented in Table 1. The cluster analysis enabled the identification of six statistical different spatial patterns in solar PV production, Figure 6. These patterns are strongly related to interactions amid the atmospheric conditions and the orography. As in the wind power production, solar power generation shows a seasonal behavior, with high positive gradients from sunrise on, reaching a plateau with negligible gradients around noon (11-14 h) and high negative gradients until sunset, Figure 7a. Figure 7b highlights the seasonal contrast of weather conditions between the southern and northern regions. During the March to October months, no substantial differences are found in power production. In contrast, in the winter months, depending on the region, the dissimilarity in power production can reach 0.5 MW. This behavior leads to significant differences in the average solar PV power daily profile and in the current capacity installed. From an annual energy power production point of view, and as shown in Figure 6a, the south/center regions of Portugal present more favorable conditions regarding the exploitation of this technology compared with the north region. By comparing the WPP ( Figure 5) and the SPP (Figure 7) profiles, a high temporal complementarity degree can be observed, namely, between WPPs #3, #4, #6, and #8 and all SPP profiles. In these cases, the wind production profiles show the lowest production period during the day, when solar production is higher, and vice versa.

Scenarios and Optimization Perspectives
The VRE deployment scenarios and optimization perspectives analyzed in this work are presented in Table 1. The Wind STD and Wind HSC are based on the wind power deployment, taking into consideration the wind resource complementarity with the presently installed wind power plants. In the PV STD and PV HSC scenario, a similar approach is taken by considering only the PV technology and the SPP identified. In the remain scenarios, a joint analysis of both technologies' generation to attain Energies 2020, 13, 4132 12 of 21 the demand is performed by taking advantage of the complementarity between the wind and solar generation. For each scenario, the two objective functions are applied.

Additional VRE Deployment Capacity and Daily Profiles
Through the greedy algorithm and the power production time series associated with each WPP or SPP, it was possible to identify the additional power capacity needed to comply with an annual VRE energy exceedance of 5%, Figure 8. In Figure 9, the average net load profiles associated with each scenario are depicted.

HSC WindHSC
Minimize one-hour net load step change standard deviation by adding: The WindSTD and WindHSC are based on the wind power deployment, taking into consideration the wind resource complementarity with the presently installed wind power plants. In the PVSTD and PVHSC scenario, a similar approach is taken by considering only the PV technology and the SPP identified. In the remain scenarios, a joint analysis of both technologies' generation to attain the demand is performed by taking advantage of the complementarity between the wind and solar generation. For each scenario, the two objective functions are applied.

Additional VRE Deployment Capacity and Daily Profiles
Through the greedy algorithm and the power production time series associated with each WPP or SPP, it was possible to identify the additional power capacity needed to comply with an annual VRE energy exceedance of 5%, Figure 8. In Figure 9, the average net load profiles associated with each scenario are depicted. Scenarios that use the NL standard deviation as the function objective (STD scenarios) permit to install a higher power capacity when compared with the HSC scenarios. Although, strong one-hour step changes are expected, especially in the PV and PV + Wind STD scenarios, Figure 9a. According to Figure 9a, a strong reinforcement of PV power capacity can induce, on average, a net load below zero at 1 and 2 p.m. Moreover, the highest amplitudes between the minimum and maximum of the net load are also expected for the PV scenarios. Thus, with a strong reinforcement of the PV capacity, the so-called "duck curve" is also observed [58]. This effect can raise several concerns regarding the power systems flexibility, namely the capability of the conventional power plants to accommodate the ramp rate range needed to fully exploit solar energy without (i) increasing the overall costs of the system and (ii) reducing the environmental benefits of VRE [58,59]. The regions with a higher solar PV potential are the ones with a high additional capacity-SPP #2 and #5, Figure 8. While in the STD scenario, the capacity is installed mainly in SPP #5, in the HSC scenario, the required power capacity In the case of the Wind scenarios, and taking into account the current wind power capacity, Figure 4c, a strong deployment in regions 9 and 10 is needed to accomplish the criteria established (minimization of wind power variability and the one-hour step change). In this case, a strategic decentralization of the wind power deployment capacity can enable meeting the demand without a strong net load variability. This solution involves wind power deployment in coastal/offshore zones, which can bring further benefits, once those power plants will always be near the larger consumption centers, contrary to the current situation. In WPP10, one of the regions with the highest wind power capacity in all scenarios, the first floating Portuguese offshore floating park with a 27 MW nominal capacity (already in the planning phase) is foreseen. The continuing decline in costs of the offshore technologies can be conducive to further investments in this technology [60]. WPP 5, which is one of the regions with a high wind power potential, presents a low additional capacity value. Thus, a strategic wind power deployment strategy can pass to exploit locations with wind power profiles that allow following the demand patterns, even if this represents less annual energy production. Additional wind power capacity is also required in regions with high installed capacity, e.g., WPP6. In these cases, the additional capacity needed can be also achieved through repowering procedures. For wind scenarios, the minimum NL values are expected during the first hours of the day and the maximum values are expected between 10 a.m. and 1 p.m., where it was not possible to identify an adequate region to meet the electricity demand. This drawback can be easily overcome by exploring the power production complementarity between this technology and solar PV.
In the PV + Wind scenarios, a strong complementarity between the VRE generation is observed. The solar PV plants can overcome the limitations observed in the wind power generation during the period between 10 a.m. and 1 p.m. The results of PV + WindSTD scenarios are similar to PV scenario results concerning the total power capacity, notwithstanding in the first one (i) a shift in the intensity of the net load is introduced by the wind power capacity during the night time and (ii) in average, the net load values are always higher than 0 MW, since less solar PV power capacity is installed. In the case of PV + WindHSC scenario, the complementarity enables to obtain a flat net load profile during the day. The daily amplitude is nearly 1500 MW with reduced hourly step-changes enabling to decrease the power system flexibility requirements needs. In this case, the highest net load values are expected by the end of the day (7 to 10 p.m.). Comparing both PV + Wind scenarios, the use of HSC objective functions requires less installed capacity.
Further details of the PV + WindHSC scenario are depicted in Figure 9b. According to this figure, it is possible to increase the VRE power capacity in other hand, the peak between 7 to 10 p.m. is still observed. Notably, this peak hour, even for large scale integration of VRE, is only reduced in the Wind scenarios.  The reduced NL values achieved means that is possible to replace conventional generation by renewable generation, bringing environmental benefits, due to a reduction in the CO2 emissions [61]. Impacts are also expected in the electricity markets on the basis of the marginal cost of each technology (the cost necessary to produce a megawatt) in the formation of its hourly price. Since VRE technologies have very low marginal costs (of the order of 0 €/MWh), compared to conventional technologies, there is a tendency for a decrease in the values verified in the wholesale markets-the

Hour
Current 1% 2% 3% 4% 5% 6% 7% 8% 9% 10% Scenarios that use the NL standard deviation as the function objective (STD scenarios) permit to install a higher power capacity when compared with the HSC scenarios. Although, strong one-hour step changes are expected, especially in the PV and PV + Wind STD scenarios, Figure 9a. According to Figure 9a, a strong reinforcement of PV power capacity can induce, on average, a net load below zero at 1 and 2 p.m. Moreover, the highest amplitudes between the minimum and maximum of the net load are also expected for the PV scenarios. Thus, with a strong reinforcement of the PV capacity, the so-called "duck curve" is also observed [58]. This effect can raise several concerns regarding the power systems flexibility, namely the capability of the conventional power plants to accommodate the ramp rate range needed to fully exploit solar energy without (i) increasing the overall costs of the system and (ii) reducing the environmental benefits of VRE [58,59]. The regions with a higher solar PV potential are the ones with a high additional capacity-SPP #2 and #5, Figure 8. While in the STD scenario, the capacity is installed mainly in SPP #5, in the HSC scenario, the required power capacity is split by SPPs #2 and #5. This result can be partially explained by the need to reduce severe hourly power step-changes imposed by the objective function. Thus, the algorithm selects these two regions to promote the statistical power smoothing effect. Moreover, the northern region of Portugal is more susceptible to severe meteorological events, which impact the solar irradiation, unleashing severe hourly power step-changes.
In the case of the Wind scenarios, and taking into account the current wind power capacity, Figure 4c, a strong deployment in regions 9 and 10 is needed to accomplish the criteria established (minimization of wind power variability and the one-hour step change). In this case, a strategic decentralization of the wind power deployment capacity can enable meeting the demand without a strong net load variability. This solution involves wind power deployment in coastal/offshore zones, which can bring further benefits, once those power plants will always be near the larger consumption centers, contrary to the current situation. In WPP10, one of the regions with the highest wind power capacity in all scenarios, the first floating Portuguese offshore floating park with a 27 MW nominal capacity (already in the planning phase) is foreseen. The continuing decline in costs of the offshore technologies can be conducive to further investments in this technology [60]. WPP 5, which is one of the regions with a high wind power potential, presents a low additional capacity value. Thus, a strategic wind power deployment strategy can pass to exploit locations with wind power profiles that allow following the demand patterns, even if this represents less annual energy production. Additional wind power capacity is also required in regions with high installed capacity, e.g., WPP6. In these cases, the additional capacity needed can be also achieved through repowering procedures. For wind scenarios, the minimum NL values are expected during the first hours of the day and the maximum values are expected between 10 a.m. and 1 p.m., where it was not possible to identify an adequate region to meet the electricity demand. This drawback can be easily overcome by exploring the power production complementarity between this technology and solar PV.
In the PV + Wind scenarios, a strong complementarity between the VRE generation is observed. The solar PV plants can overcome the limitations observed in the wind power generation during the period between 10 a.m. and 1 p.m. The results of PV + Wind STD scenarios are similar to PV scenario results concerning the total power capacity, notwithstanding in the first one (i) a shift in the intensity of the net load is introduced by the wind power capacity during the night time and (ii) in average, the net load values are always higher than 0 MW, since less solar PV power capacity is installed. In the case of PV + Wind HSC scenario, the complementarity enables to obtain a flat net load profile during the day. The daily amplitude is nearly 1500 MW with reduced hourly step-changes enabling to decrease the power system flexibility requirements needs. In this case, the highest net load values are expected by the end of the day (7 to 10 p.m.). Comparing both PV + Wind scenarios, the use of HSC objective functions requires less installed capacity.
Further details of the PV + Wind HSC scenario are depicted in Figure 9b. According to this figure, it is possible to increase the VRE power capacity in a sustainable way without introducing strong net load variations. Compared with the present net load profile, the highest differences are observed during the day time, due to the solar PV generation. On the other hand, the peak between 7 to 10 p.m. is still observed. Notably, this peak hour, even for large scale integration of VRE, is only reduced in the Wind scenarios.
The reduced NL values achieved means that is possible to replace conventional generation by renewable generation, bringing environmental benefits, due to a reduction in the CO 2 emissions [61]. Impacts are also expected in the electricity markets on the basis of the marginal cost of each technology (the cost necessary to produce a megawatt) in the formation of its hourly price. Since VRE technologies have very low marginal costs (of the order of 0 €/MWh), compared to conventional technologies, there is a tendency for a decrease in the values verified in the wholesale markets-the order of merit effect [61,62]. For certain levels of VRE penetration, this effect is positive for consumers, due to a reduction in electricity market prices. On the other side, this effect leads to a reduction in producers profitability ("self-cannibalization effect") [62]. This situation may lessen the incentive to invest in new capacity deployment, preventing to achieve high shares of VRE generation [3]. Thus, adequate policy schemes are needed to promote the nearly 100 % renewable power system without increasing the bill costs paid by consumers [1,3,9]. Figure 10 shows the load duration curve for each scenario, considering 5% and 10% annual energy exceedance values, respectively.

Net Load Duration Curve and Extreme Values
The reduced NL values achieved means that is possible to replace conventional generation by renewable generation, bringing environmental benefits, due to a reduction in the CO2 emissions [61]. Impacts are also expected in the electricity markets on the basis of the marginal cost of each technology (the cost necessary to produce a megawatt) in the formation of its hourly price. Since VRE technologies have very low marginal costs (of the order of 0 €/MWh), compared to conventional technologies, there is a tendency for a decrease in the values verified in the wholesale markets-the order of merit effect [61,62]. For certain levels of VRE penetration, this effect is positive for consumers, due to a reduction in electricity market prices. On the other side, this effect leads to a reduction in producers profitability ("self-cannibalization effect") [62]. This situation may lessen the incentive to invest in new capacity deployment, preventing to achieve high shares of VRE generation [3]. Thus, adequate policy schemes are needed to promote the nearly 100 % renewable power system without increasing the bill costs paid by consumers [1,3,9]. Figure 10 shows the load duration curve for each scenario, considering 5% and 10% annual energy exceedance values, respectively.

Net Load Duration Curve and Extreme Values
From Figure 10, it is possible to observe that the PV + Wind scenarios are the most appropriate to reduce and smooth the net load values. By taking advantage of the generation complementarity, these scenarios show the lowest net load: (i) amplitude between the maximum and minimum values observed and (ii) percentage of time with negatives values, i.e., less surplus generation during the period under analysis is observed. On the other hand, the PV scenario shows (i) the lowest capability to reduce the net load and (ii) the highest net load amplitude values. These results, associated with the intrinsic features of PV generation, i.e., no production during the night time and peak production during the middle of the day, are more pronounced as the solar PV share into the power system increases. From Figure 10, it is possible to observe that the PV + Wind scenarios are the most appropriate to reduce and smooth the net load values. By taking advantage of the generation complementarity, these scenarios show the lowest net load: (i) amplitude between the maximum and minimum values observed and (ii) percentage of time with negatives values, i.e., less surplus generation during the period under analysis is observed. On the other hand, the PV scenario shows (i) the lowest capability to reduce the net load and (ii) the highest net load amplitude values. These results, associated with the intrinsic features of PV generation, i.e., no production during the night time and peak production during the middle of the day, are more pronounced as the solar PV share into the power system increases.
The Wind scenarios present a slightly better performance to reduce NL values when compared with the PV scenarios. Nevertheless, similar extreme negative values are also observed. This result can be partly explained by the weather conditions, where the power smoothing effect associated with the aggregation of wind parks, even when widely dispersed, is not enough to mitigate some extreme meteorological events, such cold fronts [40]. This situation usually occurs when the entire country is immersed under weather conditions with a coherent structure, which usually extends over several hundreds of kilometers (e.g., cold fronts). Therefore, both scenarios based on a single VRE show the lowest performance, meaning that the natural complementarity of the wind and solar PV is essential to mitigate extreme values. Therefore, the complementarity enables to reduce the power system flexibility requirements associated with a large-scale integration of VRE.
The percentage of time with positive net load values is very similar among the different scenarios, 88% for the 5% annual energy exceedance, while for the 10% scenario, this time is reduced to 83%. In the solar scenarios, this value is reduced to 96% of the time. Remarkably, for 98% of the time of the year, the existing pumped hydro storage (PHS) capacity in Portugal could already absorb the surplus generated in the Wind and PV + Wind scenarios with a 10% energy surplus.
In Table 2, further insights regarding the seasonal NL extreme values are presented. According to Table 2, for the PV scenarios, the maximum values are equal to the present day NL and are observed at 8 p.m. This scenario emphasizes that even with a large penetration of solar PV technology, due to its intrinsic behavior, the required capacity from other technologies to satisfy the demand needed is above 8 GW in the winter and in the summer is nearly 6 GW. These high values, that need to be compensated using other technologies, raise special concerns during the summer, due to the expected reduction of the available resource for hydropower plants and the foreseen decommissioning of the Portuguese conventional coal power plants. The Wind and PV + Wind scenarios show that is possible to reduce the current seasonal 98th NL percentile value. Considering a Surplus Max. equal to 5%, the minimum values of this statistical parameter are typically observed in the PV + WIND HSC scenario. However, for a Surplus Max. equal to 10%, the Wind STD scenario is the one that shows the best performance. Thus, a strategic synergy between wind and solar PV generation enables to reduce the power system dependence from other technologies. Notwithstanding, meaningful NL values are observed. Consequently, for a large-scale integration of renewable power system, the results suggest that a long-term holistic approach considering (i) batteries and hydro pumping, (ii) demand-side response, (iii) interconnection capacity and (iv) sector coupling is crucial to dully allocate the energy surplus/deficit from VRE [9,12].
Alongside the power production, the hourly gradients of power production are of particular concern during the electric power operation [8]. They represent a key challenge for balancing the energy system by adjusting the production of dispatchable power plants (running or start-up/shut-down) and storage. In Figure 11, the duration curve of the net load hourly step change is shown to assess the variability obtained in each scenario.
On the basis of Figure 11, the highest hourly step-changes are associated with PV scenarios, with the most severe ramp-down reaching −2650 MW and −2834 MW for 5% and 10% surplus values, respectively. Regarding hourly ramp-up, the maximum values are 3290 MW and 3545 MW. Due to the daily variability of solar PV generation, the scenarios based on this technology show the highest extreme values of this parameter. In the remain scenarios, a high percentage of time with hourly net load step change near 0 MW is observed. Notably, the PV + Wind HSC scenario shows a behavior very similar to the current duration curve. Thus, the results suggest that is possible to increase the VRE capacity without exacerbating the current flexibility to deal with the variability associated with these renewable sources. On the basis of Figure 11, the highest hourly step-changes are associated with PV scenarios, with the most severe ramp-down reaching −2650 MW and −2834 MW for 5% and 10% surplus values, respectively. Regarding hourly ramp-up, the maximum values are 3290 MW and 3545 MW. Due to the daily variability of solar PV generation, the scenarios based on this technology show the highest extreme values of this parameter. In the remain scenarios, a high percentage of time with hourly net load step change near 0 MW is observed. Notably, the PV + WindHSC scenario shows a behavior very similar to the current duration curve. Thus, the results suggest that is possible to increase the VRE capacity without exacerbating the current flexibility to deal with the variability associated with these renewable sources.

VRE Share Values in the Final Consumption
In Figure 12, the VRE share in the final consumption is depicted considering different Surplus . values. As expected, due to the wind and PV generation synergy, the PV + Wind scenarios enable to obtain the highest VRE penetration values for the same annual energy exceedance value. The highest value (scenario PV + WindSTD) enables to obtain a VRE contribution of 68%, while in the PV + WindHSC scenario, the contribution is 67%. Taking into account this slight difference found and the previous results achieved (e.g., typical power profile and hourly step-change), results suggest that the PV + WindHSC scenario can provide a more sustainable way to increase the VRE penetration.
In the PV scenarios, only a maximum of 56% share of VRE is estimated for a Surplus . equal to 10%-HSC scenario. Notably, the Wind scenarios show that a strategic deployment can enable to achieve a VRE share of 65%. After the 3% annual energy surplus threshold, the results highlight that the wind power generation complementarity across the different spatial locations allows to achieve

VRE Share Values in the Final Consumption
In Figure 12, the VRE share in the final consumption is depicted considering different a Surplus Max. values.
(a) (b) Figure 11. Duration curve of the net load hourly step change using a Surplus . equal to (a) 5% and (b) 10%.
On the basis of Figure 11, the highest hourly step-changes are associated with PV scenarios, with the most severe ramp-down reaching −2650 MW and −2834 MW for 5% and 10% surplus values, respectively. Regarding hourly ramp-up, the maximum values are 3290 MW and 3545 MW. Due to the daily variability of solar PV generation, the scenarios based on this technology show the highest extreme values of this parameter. In the remain scenarios, a high percentage of time with hourly net load step change near 0 MW is observed. Notably, the PV + WindHSC scenario shows a behavior very similar to the current duration curve. Thus, the results suggest that is possible to increase the VRE capacity without exacerbating the current flexibility to deal with the variability associated with these renewable sources.

VRE Share Values in the Final Consumption
In Figure 12, the VRE share in the final consumption is depicted considering different Surplus . values. As expected, due to the wind and PV generation synergy, the PV + Wind scenarios enable to obtain the highest VRE penetration values for the same annual energy exceedance value. The highest value (scenario PV + WindSTD) enables to obtain a VRE contribution of 68%, while in the PV + WindHSC scenario, the contribution is 67%. Taking into account this slight difference found and the previous results achieved (e.g., typical power profile and hourly step-change), results suggest that the PV + WindHSC scenario can provide a more sustainable way to increase the VRE penetration.
In the PV scenarios, only a maximum of 56% share of VRE is estimated for a Surplus . equal to 10%-HSC scenario. Notably, the Wind scenarios show that a strategic deployment can enable to achieve a VRE share of 65%. After the 3% annual energy surplus threshold, the results highlight that the wind power generation complementarity across the different spatial locations allows to achieve  As expected, due to the wind and PV generation synergy, the PV + Wind scenarios enable to obtain the highest VRE penetration values for the same annual energy exceedance value. The highest value (scenario PV + Wind STD ) enables to obtain a VRE contribution of 68%, while in the PV + Wind HSC scenario, the contribution is 67%. Taking into account this slight difference found and the previous results achieved (e.g., typical power profile and hourly step-change), results suggest that the PV + Wind HSC scenario can provide a more sustainable way to increase the VRE penetration.
In the PV scenarios, only a maximum of 56% share of VRE is estimated for a Surplus Max. equal to 10%-HSC scenario. Notably, the Wind scenarios show that a strategic deployment can enable to achieve a VRE share of 65%. After the 3% annual energy surplus threshold, the results highlight that the wind power generation complementarity across the different spatial locations allows to achieve a higher penetration share compared to the PV scenarios. Therefore, it is possible to observe a 1% difference between the Wind and PV scenarios for a 3% energy exceedance, while for the 10% energy exceedance, this difference increases to 4%, considering the best result of each scenario.

Final Remarks
In this study, a methodology to assess the complementarity of wind and solar power production and their contribution to meet the national demand was developed and applied using Portugal as a case study.
Through the clustering analysis of high-resolution data, the spatiotemporal variability of the wind and solar PV power production in Portugal was identified. Ten regions with different wind power profile features were detected, while for solar PV, six regions were identified. On the basis of these profiles and a greedy algorithm, different deployment scenarios are considered by strategically adding to the current energy mix: (1) only wind power, (2) only solar PV power, (3) adding wind and solar PV power. PV scenarios show the highest values regarding the net load hourly step change and the lowest capability to reduce the net load values throughout the year. Scenarios based on wind power technology show a higher VRE penetration share in the annual consumption when compared with strong PV solar reinforcement scenarios. However, the highest penetration share was attained by exploring the wind and solar PV generation complementarity. In the PV + Wind HSC scenario, a sustainable way to increase the VRE penetration until 68%, considering and annual energy surplus regarding the consumption equal to 10%, was identifying and considering a wind power capacity reinforcement of 4240 MW (through new wind parks but also considering repowering procedures) and 2200 MW of solar PV. These numbers, together with the current share of dispatchable hydropower production (nearly 20%), show that it is possible to obtain a sustainable pathway toward decarbonization of the Portuguese electric power system, without exacerbating the power system's flexibility requirements to deal with VRE variability. The methodology developed can be generalized and directly transferred to other regions/power systems.
The proposed methodology only addresses the temporal and spatial integration aspects of VRE variability, and the benefits of the complementarity assessment form a power system point-of-view. Thus, some limitations of the analysis should be highlighted. First, it does not cover other relevant aspects such as economic costs, land availability, social acceptance, environmental benefits, use of storage technologies, and the electricity market. Second, the electricity demand was considered inelastic. These aspects will be analyzed in future works. Notwithstanding these limitations, the result of this work clearly highlights the added value of using wind and solar PV complementarity and electricity criteria as a planning strategy for new VRE capacity deployment aiming to reduce the power flexibility needs, namely, the use of expensive energy storage systems.
Although a long period is needed to obtain results with high statistical significance, the results suggest that high net load values are still expected, even under conditions of large-scale integration of VRE. Thus, in this energy transition phase, a long-term holistic approach with a cost analysis (e.g., strategic use of storage, hydro pumping, sector coupling, and hydrogen generation) is required to sustainably (i) deal with the expected energy surplus/deficit from VRE and (ii) create public policies to fairly promote the deployment of wind and solar PV towards the mitigation of climate change.
Author Contributions: A.C. performed the literature review, collected, and processed all the data, developed the mathematical model, implemented the model, and wrote the manuscript. A.E. discussed the work to be developed, supported the model developed and the identification of the most relevant results to be shown, and also reviewed the manuscript. All authors have read and agreed to the published version of the manuscript.