High-Resolution Solar Climate Atlas for Greece under Climate Change Using the Weather Research and Forecasting (WRF) Model

In the context of climate change and growing energy demand, solar technologies are considered promising solutions to mitigate Greenhouse Gas (GHG) emissions and support sustainable adaptation. In Greece, solar power is the second major renewable energy, constituting an increasingly important component of the future low-carbon energy portfolio. In this work, we propose the use of a high-resolution regional climate model (Weather Research and Forecasting model, WRF) to generate a solar climate atlas for the near-term climatological future under the Representative Concentration Pathway (RCPs) 4.5 and 8.5 scenarios. The model is set up with a 5 × 5 km2 spatial resolution, forced by the ERA-INTERIM for the historic (1980–2004) period and by the EC-EARTH General Circulation Models (GCM) for the future (2020–2044). Results reaffirm the high quality of solar energy potential in Greece and highlight the ability of the WRF model to produce a highly reliable future climate solar atlas. Projected changes between the annual historic and future RCPs scenarios indicate changes of the annual Global Horizontal Irradiance (GHI) in the range of ±5.0%. Seasonal analysis of the GHI values indicates percentage changes in the range of ±12% for both scenarios, with winter exhibiting the highest seasonal increases in the order of 10%, and autumn the largest decreases. Clear-sky fraction fclear projects increases in the range of ±4.0% in eastern and north continental Greece in the future, while most of the Greek marine areas might expect above 220 clear-sky days per year.


Introduction
The usage of solar energy applications in Greece has increased significantly during the last decade, providing both heat energy and electricity, meeting the demands for more efficient and more affordable production of solar technologies. Globally, the cumulative capacities for solar thermal and photovoltaics (PV) installations were around 480 and 506 GW, respectively, by the end of 2018 [1]. Between the top 10 countries who share their electricity production from solar energy in 2018, Greece holds the 9th place for the added capacity of solar water heating collectors and the 2nd place for solar PV, with a percentage of 8.2%, while the renewable power target for solar PV installed capacity has been set at 2.0 GW, by the end of 2030 [1].
With respect to the country's commitments for climate change mitigation and to implement European (EU) policy actions to decrease Greenhouse Gas (GHG) emissions, the upward trend of solar technologies installations is anticipated to further strengthen. In order to support coherent planning for Figure 1 presents the overall methodology of the development of the high-resolution solar climate atlas for Greece under climate change, using the WRF model that is based on extensive validation against the ECA dataset and bias correction of historic simulations and future projections. The solar Atmosphere 2020, 11, 761 3 of 17 resource assessment of the historic climatology period is performed using both measured and simulated solar irradiance data. As there are not sufficient weather stations to provide historical solar irradiance measurements in Greece, the use of high-resolution regional climate simulations to generate consistent estimations in future climate, at the seasonal to 25-year time periods, is of high priority [14].
Atmosphere 2020, 11, x FOR PEER REVIEW 3 of 17 and simulated solar irradiance data. As there are not sufficient weather stations to provide historical solar irradiance measurements in Greece, the use of high-resolution regional climate simulations to generate consistent estimations in future climate, at the seasonal to 25-year time periods, is of high priority [14].

Data Sources
The solar climate atlas is established using the Advanced Weather Research and Forecasting numerical model (v3.6.1) [22], forced by ERA-INTERIM [23] for the historic period  and EC-EARTH General Circulation Models (GCM) developed by the European Centre for Medium-Range Weather Forecasts (ECMWF) [24], for the same historic and future (2020-2044) periods, dynamically downscaled to the region of Greece at a scale of 5 × 5 km 2 [18,19]. EC-Earth is an Earth System Model (ESM) closely aligned with the ECMWF seasonal forecasting system, hence benefiting from its continuous updates. The model is widely used in a range of international climate modeling and service research projects. EC-Earth has been designed to cover all timescales (from subseasonal to long-term) and includes components that introduce more complexity into the system such as land-use changes and interactive vegetation in the simulations [24]. WRF was set up with two nested grids ( Figure 2): the outer domain, which is centered in the Mediterranean basin and has been configured using a mesh size of 265 × 200 grid cells with 20 km horizontal resolution (d01-EU), and the internal domain, which uses 185 × 185 grid cells of 5 km horizontal grid spacing covering the area of Greece (d02-GR). This setup has used 40 levels in the vertical dimension, arranged according to terrain-following hydrostatic pressure vertical coordinates, and one-way nesting has been applied, to avoid possible noise during feedback from the inner domain to the coarse domain. The last four days of the previous month were regarded as model spin-up for the following month, so the model was reinitialized every month [19]. According to [19], the ERA-INTERIM reanalysis dataset was used to provide initial and boundary conditions. The lateral boundary conditions and the sea surface temperature were both updated every 6 h, from ERA-INTERIM. Seven simulations were run using spectral nudging above the planetary boundary layer (PBL) and only over the coarse domain. The model physics configuration and extensive validation for temperatures and precipitation are presented in [18].
Climate simulations with the WRF model are performed for two periods (a) from 1980 to 2004, used for model validation purposes, and (b) 2020-2044, for the two RCPs (4.5 W/m 2 and 8.5 W/m 2 ), which is linked to the projected lifespan of modern solar systems [3].

Data Sources
The solar climate atlas is established using the Advanced Weather Research and Forecasting numerical model (v3.6.1) [22], forced by ERA-INTERIM [23] for the historic period  and EC-EARTH General Circulation Models (GCM) developed by the European Centre for Medium-Range Weather Forecasts (ECMWF) [24], for the same historic and future (2020-2044) periods, dynamically downscaled to the region of Greece at a scale of 5 × 5 km 2 [18,19]. EC-Earth is an Earth System Model (ESM) closely aligned with the ECMWF seasonal forecasting system, hence benefiting from its continuous updates. The model is widely used in a range of international climate modeling and service research projects. EC-Earth has been designed to cover all timescales (from subseasonal to long-term) and includes components that introduce more complexity into the system such as land-use changes and interactive vegetation in the simulations [24]. WRF was set up with two nested grids ( Figure 2): the outer domain, which is centered in the Mediterranean basin and has been configured using a mesh size of 265 × 200 grid cells with 20 km horizontal resolution (d01-EU), and the internal domain, which uses 185 × 185 grid cells of 5 km horizontal grid spacing covering the area of Greece (d02-GR). This setup has used 40 levels in the vertical dimension, arranged according to terrain-following hydrostatic pressure vertical coordinates, and one-way nesting has been applied, to avoid possible noise during feedback from the inner domain to the coarse domain. The last four days of the previous month were regarded as model spin-up for the following month, so the model was reinitialized every month [19]. According to [19], the ERA-INTERIM reanalysis dataset was used to provide initial and boundary conditions. The lateral boundary conditions and the sea surface temperature were both updated every 6 h, from ERA-INTERIM. Seven simulations were run using spectral nudging above the planetary boundary layer (PBL) and only over the coarse domain. The model physics configuration and extensive validation for temperatures and precipitation are presented in [18].
Climate simulations with the WRF model are performed for two periods (a) from 1980 to 2004, used for model validation purposes, and (b) 2020-2044, for the two RCPs (4.5 W/m 2 and 8.5 W/m 2 ), which is linked to the projected lifespan of modern solar systems [3].

Derived Solar Climatology Parameters
In our study, data of the GHI, from the ECA dataset [21], which has been provided by the Hellenic National Meteorological Service (HNMS), were used to validate the high-resolution climatology simulations generated by WRF. The following section introduces the main outputs from the WRF regional climate simulations. WRF model output data (GHI and fclear) were obtained at a frequency of 6 h intervals and then transformed into GHI climatology annual and seasonal values of winter (DJF), spring (MAM), summer (JJA), and autumn (SON).

Mean Daily Global Horizontal and Direct Normal Irradiation
The amount of radiation on a terrestrial surface at a given site, for a given time, depends on the orientation and slope of the surface. A flat surface receives the global radiation, which is the sum of DNI and Diffuse Horizontal Irradiation (DHI), while θ is the zolar zenith angle, between the sun's rays and the vertical (0 ≤ θ ≤90).
GHI is the amount of direct and diffuse irradiance absorbed on a horizontal plane and is reference irradiance for the assessment of climatic zones and a critical component for estimating the radiation on a tilted surface. Despite the fact that DNI is a key factor for energy yield calculation and performance assessment of CSP, concentrator solar photovoltaic (CPV) technologies, and tilted or sun-tracking photovoltaic modules, for the purposes of this study, the GHI has been preferred, as it is the most crucial factor for energy yield calculation and performance assessment of flat-plate

Derived Solar Climatology Parameters
In our study, data of the GHI, from the ECA dataset [21], which has been provided by the Hellenic National Meteorological Service (HNMS), were used to validate the high-resolution climatology simulations generated by WRF. The following section introduces the main outputs from the WRF regional climate simulations. WRF model output data (GHI and f clear ) were obtained at a frequency of 6 h intervals and then transformed into GHI climatology annual and seasonal values of winter (DJF), spring (MAM), summer (JJA), and autumn (SON).

Mean Daily Global Horizontal and Direct Normal Irradiation
The amount of radiation on a terrestrial surface at a given site, for a given time, depends on the orientation and slope of the surface. A flat surface receives the global radiation, which is the sum of DNI and Diffuse Horizontal Irradiation (DHI), while θ z is the zolar zenith angle, between the sun's rays and the vertical (0 ≤ θ z ≤ 90).
GHI is the amount of direct and diffuse irradiance absorbed on a horizontal plane and is reference irradiance for the assessment of climatic zones and a critical component for estimating the radiation on a tilted surface. Despite the fact that DNI is a key factor for energy yield calculation and performance assessment of CSP, concentrator solar photovoltaic (CPV) technologies, and tilted or sun-tracking photovoltaic modules, for the purposes of this study, the GHI has been preferred, as it is the most crucial factor for energy yield calculation and performance assessment of flat-plate photovoltaic (PV) technologies [25]. Average GHI (W/m 2 ) was obtained at 6 h intervals and then transformed into mean annual and seasonal values of the 25-year time period.

Sunshine Duration
Sunshine duration can be calculated using the following equation: where T is the sunshine duration in hours, ϕ is the latitude in degrees, and δ is the declination of the sun in degrees. The declination of the sun (δ) is given by Equation (3) As first, the assessment of the WRF model tended to heavily overestimate the mean GHI (W/m 2 ) values over the main parts Greece when compared to the ECA dataset, due to the fact that the WRF model provides GHI data (W/m 2 ) at 6 h intervals, thus it was preferred to transform the GHI from 6 h intervals (W/m 2 ) to continuous distribution per hours of daytime. The maxGHI value per day was then calculated through the sin −1 (t), where t corresponds to the 12:00 pm at each grid point.
where t2 corresponds to the time of the day and takes values from the range 0 < π/2 < π. Thus, the mean daily, annual, and seasonal solar GHI (kWh/m 2 ) was calculated, with the use of the mean daily GHI (W/m 2 ) and the sunshine duration (hours), for all the grid points of the Greek domain.

Clear-Sky Days and f clear
In general, clear-sky models estimate the terrestrial solar radiation under a cloudless sky, as a function of the solar elevation angle, site altitude, aerosol concentration, water vapor, and various atmospheric conditions. Apparently, the solar downward radiation is mostly affected by cloudy conditions. Clouds significantly decrease the incoming solar irradiance reaching the panels' surface, causing losses in the overall energy production. Indeed, clouds are able (i) to block a portion of the solar irradiance, decreasing mostly the direct radiation on the surface level [26]; (ii) to enhance the diffusive irradiation that arrives to the ground; and (iii) to regulate the global temperature and to contribute to global warming by capturing the long-wave radiation [17].
In this study, the clear-sky days were calculated at 6 h intervals from the WRF model output, through the difference of the total insolation for the actual cloud cover (GHI) and the clear-sky total insolation for when there were no clouds present (GHIcs). Also, the fraction of time for which no significant clouds blocked the sun (f clear ) was calculated using Suehrcke's Equation [27]. Following the methodological approach of [14], f clear is found by: where K is the mean clearness index and Kclear is the mean clear-sky clearness index.

WRF Solar Climatology Performance Validation
The implemented regional climate model WRF forced by the ERA-INTERIM performance was systematically assessed via a reliable evaluation procedure with the use of historical observations. The "nearest neighbor" method to the position of the monitoring stations was chosen. Section 3.1 describes the comparison against the available ECA dataset [21], from as many Greek meteorological stations as possible, of the same temporal period each time. Moreover, a qualitative comparison of GHI (kWh/m 2 ) with the annual climatological solar energy map for Greece for GHI values (kWh/m 2 ) of the period of 2002-2012, from the Atmosphere 2020, 11, 761 6 of 17 work of [16], was performed (Section 3.2). Also, discussion and comments were introduced, related to the qualitative comparison of WRF model results with the 35-year climatology simulations of the Solargis solar atlas published by the World Bank Group [28], funded by the Energy Sector Management Assistance Program (ESMAP), and with the 15-year climatology of the mean monthly GHI solar energy map of Greece, from the European Organization for the Exploitation of Meteorological Satellites (EUMETSAT) Satellite Application Facility on Climate Monitoring (CM SAF), at a spatial resolution of about 5 km.

ECA Dataset
Although Greece has a long history of routine meteorological observations mainly provided by the HNMS, the publications of high-quality data are not always available, as a matter of several issues (data policy, etc.). European Climate Assessment (ECA) concentrates its efforts on the implementation of holistic analysis methodologies to daily measurements datasets. Greece is one of the 34 countries who joined in the ECA project providing data from 10 stations. From them, those of Trikala, N. Filadelfeia, and Tymbakion have been excluded, as they had extensive gaps in their datasets [21]. Meanwhile, ECA has already applied a uniform comparison of the ECA dataset with existing gridded datasets, showing correlation coefficients greater than 0.8, which strengthens the ECA validity and confidence.
Under the purpose of this study, the climatology solar irradiance (W/m 2 ) of the high-resolution WRF climate simulation over Greece has been compared to the available climatology dataset of the ECA project [20] using the following statistical measures, presented in Table 1: (i) correlation coefficient (r), Equation (6); (ii) mean bias, Equation (7); (iii) mean absolute percentage error (MAPE), Equation (8); and (iv) the root mean square error (RMSE), Equation (9).
where s i and s represent the model simulations and mean monthly model simulated values, respectively, m i and m are the measurements (observations) and mean monthly measurements, respectively, and N the total number of data points. The evaluation of simulated mean surface GHI values (black line) against ECA observations (grey line) provided by the HNMS network to ECA yields that the stations of Arta, Aliartos, Argos, and Thessaloniki indicate the greatest correlation factors. Also in the locations of Athens, Aliartos, Andravida, and Sikiwna, model biases are within 40-45 W/m 2 (Table 1). MAPE is below 25% for the location of Athens, while the lowest RMSE values have been noticed for the locations of Athens and Aliartos. However, as Table 1 demonstrates, during the statistical comparison between WRF climatology and ECA dataset, the error analysis results yielded that the WRF model (black line) tended to overestimate constantly the GHI values of the ECA dataset (grey line), presenting a systematic bias error against observed values, with average values of the order of 57 W/m 2 for the whole dataset, resulting in the need for bias adjustment. Hence, the use of bias correction was applied in order to adjust the WRF model output according to the existing climate regime, which is necessary for the production of valuable outcomes in climate applications [29,30]. Thus, the bias in mean monthly model values is corrected by subtracting the value of 57 W/m 2 , which is the difference between modeled and observed mean monthly GHI values of the time period; see Table 1. The WRF direct output is bias-corrected according to the methodology proposed in [31]. The use of bias correction improved the performance of the WRF datasets to reasonably represent the GHI climatology data (Figure 3), except for the Athens observational dataset, which shows the greatest overall agreement among the two datasets, with the highest bias error, MAPE, and RMSE values in comparison to the other stations. Apparently, the improved simulated mean GHI values (black dashed line) against ECA observations show MAPE of the range of 9.0-16.5% and RMSE from 12 to 27 W/m 2 , while the greatest overall agreement between the observations and the corrected simulated datasets was found for the stations of Andravida (1980)(1981)(1982)(1983)(1984)(1985)(1986)(1987)(1988)(1989), Aliartos (1980)(1981)(1982)(1983)(1984)(1985)(1986)(1987)(1988)(1989)(1990), and Sikiwna (1980)(1981)(1982)(1983)(1984)(1985)(1986)(1987)(1988)(1989)(1990).

Qualitative Comparison
The skill of the WRF to simulate the GHI in Greece was further compared to the annual climatological solar energy map for Greece for GHI values (kWh/m 2 ) of the period of 2002-2012, of the work of [16], who studied the challenges in solar energy resource and forecasting in Greece. The

Qualitative Comparison
The skill of the WRF to simulate the GHI in Greece was further compared to the annual climatological solar energy map for Greece for GHI values (kWh/m 2 ) of the period of 2002-2012, of the work of [16], who studied the challenges in solar energy resource and forecasting in Greece. The analysis in Figure 4 shows that WRF simulations produce a mean climatology that is very close to the [16] solar GHI map, which has also been validated to the maps provided by SOLARGIS (www.solargis.com) for the period of 1994-2013, and found differences in the order of 25 kWh/m 2 . WRF values of the GHI Atmosphere 2020, 11, x FOR PEER REVIEW 9 of 17 analysis in Figure 4 shows that WRF simulations produce a mean climatology that is very close to the [16] solar GHI map, which has also been validated to the maps provided by SOLARGIS (www.solargis.com) for the period of 1994-2013, and found differences in the order of 25 kWh/m 2 . WRF values of the GHI concerning the period of 1994-2004 (kWh/m 2 ) show an excellent agreement compared to the long-term average GHI dataset of the annual climatological solar energy map of the period of 2002-2012, with the observed differences, while the spatial variability and the geographical distribution are consistent. Also, the WRF GHI solar map was further compared to the mean solar resource map of Greece, which was released previously with Global Solar Atlas and published by the World Bank Group [28], funded by ESMAP and prepared by SOLARGIS over a 35-year period (1994-2018), showing that model simulations produce a mean climatology that is very close to the mean SOLARGIS generated climate. Results also indicate the close relationship of the WRF model simulations (W/m 2 ) to the 15-year climatology of the mean monthly GHI solar energy map of Greece from the EUMETSAT's CM SAF at a spatial resolution of about 5 × 5 km 2 [32]. These qualitative comparisons provided evidence of the robustness of the WRF model for simulating the climate conditions during the recent period  in Greece and further confidence and reliance on the future scenario simulations.

Reconstructing a Very High-Resolution Solar Climate Atlas for Greece
Mean annual and daily GHI values, along with mean seasonal annual and daily GHI values, were calculated from WRF simulations of the present climate and are shown together with the differences between future and historical climate scenarios, based on 25-year mean values. Also, the clear-sky days are presented through the difference of the total insolation for the actual cloud cover (GΗΙ) and the clear-sky total insolation for when there were no clouds present (GHIcs).

Mean Annual GHI Analysis
The computed mean annual and daily GHI for the historical period 1980-2004, using the WRF regional climate simulations at surface level, with the adjusted bias (already discussed in Section 3.1), for the Greek domain and the projected differences in the mean daily and annual sum GHI between the future (RCPs 4.5/8.5) and historical scenarios are depicted in Figure 5. In general, the Also, the WRF GHI solar map was further compared to the mean solar resource map of Greece, which was released previously with Global Solar Atlas and published by the World Bank Group [28], funded by ESMAP and prepared by SOLARGIS over a 35-year period (1994-2018), showing that model simulations produce a mean climatology that is very close to the mean SOLARGIS generated climate. Results also indicate the close relationship of the WRF model simulations (W/m 2 ) to the 15-year climatology of the mean monthly GHI solar energy map of Greece from the EUMETSAT's CM SAF at a spatial resolution of about 5 × 5 km 2 [32]. These qualitative comparisons provided evidence of the robustness of the WRF model for simulating the climate conditions during the recent period  in Greece and further confidence and reliance on the future scenario simulations.

Reconstructing a Very High-Resolution Solar Climate Atlas for Greece
Mean annual and daily GHI values, along with mean seasonal annual and daily GHI values, were calculated from WRF simulations of the present climate and are shown together with the differences between future and historical climate scenarios, based on 25-year mean values. Also, the clear-sky days are presented through the difference of the total insolation for the actual cloud cover (GHI) and the clear-sky total insolation for when there were no clouds present (GHIcs).

Mean Annual GHI Analysis
The computed mean annual and daily GHI for the historical period 1980-2004, using the WRF regional climate simulations at surface level, with the adjusted bias (already discussed in Section 3.1), for the Greek domain and the projected differences in the mean daily and annual sum GHI between the future (RCPs 4.5/8.5) and historical scenarios are depicted in Figure 5. In general, the geographical pattern and the different climatic zones beyond the Greek domain justify the variances among continental Greece and marine regions or coastlines at the same latitude [33]. geographical pattern and the different climatic zones beyond the Greek domain justify the variances among continental Greece and marine regions or coastlines at the same latitude [33]. The WRF model for the historical period  presented spatial distribution of the GHI and demonstrated significant agreement with similar works [15][16][17]34], with values above 1650 kWh/m 2 in the broader area of central and northern Greece and GHI above 1800 kWh/m 2 in the broader area of the southern marine Greek parts of the Ionian and Aegean Sea, which correspond to 4.5 and to 5.1 kWh/m 2 daily GHI values, respectively. Over the southern and southeastern parts of the Aegean and Libyan Sea, the WRF model under the historical period shows values above 1850 kWh/m 2 . However, lack of adequate average solar potential was found at the locations of the highest altitudes, commonly observed at the mountain ridges, with values in the range of 1500 kWh/m 2 , which correspond to mean daily GHI values of 4.0 kWh/m 2 ( Figure 5).
On average, analysis of the mean annual GHI variations showed values of similar magnitude and pattern between the two future scenarios. The WRF model near-term climatology changes compared to mean historic climatology are found to be in the range of −100 to 100 kWh/m 2 for both examined RCPs scenarios, which is translated to a percentage range of ±5.0%. The mean GHI values, both in future scenarios, are projected to decrease in the most marine regions of the Aegean and Ionian Seas, as well as in specific locations of the western continental parts of Greece, by the same level of magnitude (from 1.0 to 3.5%). However, both future scenarios anticipated greater increases in higher altitude locations. The most noticeable percentage increases of the GHI from 2.0 to 4.5% have been projected in the Pindus mountain range, Mount Olympus, the mountains of Crete, and so forth, where they have already been found with poor to fair quality of GHI.

Mean Seasonal GHI Analysis
During the evaluation of the GHI seasonality in Greece, a general increase is anticipated from winter (DJF) to spring (MAM) and then to summer (JJA), as a result of the increased length of the day. Apparently, an increase of sunshine ratios from winter to summer happens as the time length of days increases and the percentage variability of sunshine also grows [33]. The smallest values of GHI hours take place in the winter months, where the cyclonic (depressions) and frontal activity increase their magnitude and frequency [33]. On average, spring (MAM) and summer are meaningfully the seasons with the greater mean seasonal GHI sum, compared to those of winter and autumn.
Regarding the seasonal analysis of the mean seasonal daily GHI, as presented in Figure 6, the most intense solar radiation is expected to occur during summer with daily GHI values to be around The WRF model for the historical period  presented spatial distribution of the GHI and demonstrated significant agreement with similar works [15][16][17]34], with values above 1650 kWh/m 2 in the broader area of central and northern Greece and GHI above 1800 kWh/m 2 in the broader area of the southern marine Greek parts of the Ionian and Aegean Sea, which correspond to 4.5 and to 5.1 kWh/m 2 daily GHI values, respectively. Over the southern and southeastern parts of the Aegean and Libyan Sea, the WRF model under the historical period shows values above 1850 kWh/m 2 . However, lack of adequate average solar potential was found at the locations of the highest altitudes, commonly observed at the mountain ridges, with values in the range of 1500 kWh/m 2 , which correspond to mean daily GHI values of 4.0 kWh/m 2 ( Figure 5).
On average, analysis of the mean annual GHI variations showed values of similar magnitude and pattern between the two future scenarios. The WRF model near-term climatology changes compared to mean historic climatology are found to be in the range of −100 to 100 kWh/m 2 for both examined RCPs scenarios, which is translated to a percentage range of ±5.0%. The mean GHI values, both in future scenarios, are projected to decrease in the most marine regions of the Aegean and Ionian Seas, as well as in specific locations of the western continental parts of Greece, by the same level of magnitude (from 1.0 to 3.5%). However, both future scenarios anticipated greater increases in higher altitude locations. The most noticeable percentage increases of the GHI from 2.0 to 4.5% have been projected in the Pindus mountain range, Mount Olympus, the mountains of Crete, and so forth, where they have already been found with poor to fair quality of GHI.

Mean Seasonal GHI Analysis
During the evaluation of the GHI seasonality in Greece, a general increase is anticipated from winter (DJF) to spring (MAM) and then to summer (JJA), as a result of the increased length of the day. Apparently, an increase of sunshine ratios from winter to summer happens as the time length of days increases and the percentage variability of sunshine also grows [33]. The smallest values of GHI hours take place in the winter months, where the cyclonic (depressions) and frontal activity increase their magnitude and frequency [33]. On average, spring (MAM) and summer are meaningfully the seasons with the greater mean seasonal GHI sum, compared to those of winter and autumn.
Regarding the seasonal analysis of the mean seasonal daily GHI, as presented in Figure 6, the most intense solar radiation is expected to occur during summer with daily GHI values to be around 6.0-7.5 kWh/m 2 and immediately after the spring season (4.5-6.0 KWh/m 2 ), which are translated to 550-700 kWh/m 2 and to 400-550 kWh/m 2 of seasonal GHI, respectively. On the contrary, the mean seasonal daily GHI is considerably less in autumn, with daily GHI values of the order of 2.7-4.0 kWh/m 2 , as well as in winter, with values in the range of 1.7-2.7 kWh/m 2 , which correspond to 200-350 kWh/m 2 and to 150-250 kWh/m 2 of seasonal GHI, respectively ( Figure 6).
Atmosphere 2020, 11, x FOR PEER REVIEW 11 of 17 6.0-7.5 kWh/m 2 and immediately after the spring season (4.5-6.0 KWh/m 2 ), which are translated to 550-700 kWh/m 2 and to 400-550 kWh/m 2 of seasonal GHI, respectively. On the contrary, the mean seasonal daily GHI is considerably less in autumn, with daily GHI values of the order of 2.7-4.0 kWh/m 2 , as well as in winter, with values in the range of 1.7-2.7 kWh/m 2 , which correspond to 200-350 kWh/m 2 and to 150-250 kWh/m 2 of seasonal GHI, respectively ( Figure 6). Comparing now the calculated differences in the mean seasonal GHI averaged over the 25 years, between the historical  and the future period (2020-2044), the projected changes were in the range of −13% to +12% (Figure 6).
During winter months (DJF), WRF projects for both scenarios (RCPs 4.5 and 8.5) increases in the mean seasonal GHI values for the eastern Greek mainland and in the northwestern Aegean Sea, as well as in Euboea and Crete, which might go up to 11%. However, partial increases in the mean seasonal GHI values are expected to occur further under RCP 8.5.
Throughout the spring months (MAM), a similar pattern of the seasonal mean GHI changes is anticipated in most parts of the Greek domain during the two future scenarios. In fact, slight increases are anticipated in most regions of the Greek domain, except for the marine areas in the Aegean and Ionian Sea, where minor decreases are projected. Comparing now the calculated differences in the mean seasonal GHI averaged over the 25 years, between the historical  and the future period (2020-2044), the projected changes were in the range of −13% to +12% (Figure 6).
During winter months (DJF), WRF projects for both scenarios (RCPs 4.5 and 8.5) increases in the mean seasonal GHI values for the eastern Greek mainland and in the northwestern Aegean Sea, as well as in Euboea and Crete, which might go up to 11%. However, partial increases in the mean seasonal GHI values are expected to occur further under RCP 8.5.
Throughout the spring months (MAM), a similar pattern of the seasonal mean GHI changes is anticipated in most parts of the Greek domain during the two future scenarios. In fact, slight increases are anticipated in most regions of the Greek domain, except for the marine areas in the Aegean and Ionian Sea, where minor decreases are projected.
During the summer months (JJA), both scenarios estimated increases in higher altitude locations (Pindus mountain range, Mount Olympus, etc.) and northeastern Greece, while those increases are anticipated to be higher during the RCP 4.5 scenario (7.0-10%). For the period of autumn (SON), both scenarios estimated slight decreases in the mean seasonal daily and annual GHI down to −5.0% in most parts of the Greek domain.

Clear-Sky Fraction-f clear -and Number of Clear-Sky Days
Aerosols and clouds are the key components in the atmosphere that reflect, transmit, or trap solar radiation and play a key role in the understanding of climate change. According to [16], in cloudy conditions, the high temporal and spatial fluctuation of cloudiness is the prime cause for the assessment of solar radiation. Also, changes in cloud coverage affect the diurnal temperature (DTR) variation, especially the maximum temperature during summer months and the minimum temperatures throughout winter, as [35] stated. However, cloud cover is a factor that cannot be smoothly simulated in climate models, as it provides a high level of uncertainty in the climate change predictions [36].
Under the prospect of this study, the mean annual f clear for the historical period 1980-2004, along with the anticipated differences between the future (RCPs 4.5/8.5) and historical scenarios, are presented in Figure 7. The WRF model for the historical period , as presented in Figure 7a, indicates lower mean f clear values in the high-altitude places by the order of 0.6-0.7, while marine and coastline locations show values above 0.85. As f clear is equivalent to the sunshine fraction, future changes between RCPs 4.5 and 8.5 are estimated to be in the range of ±4.5%, with the most increases up to 4.0% found in eastern and north continental Greece, as well in Euboea, Sporades, Rhodes, and Crete, while decreases down to 3.0% are projected in the western Greek mainland and in the Aegean Sea.
During the summer months (JJA), both scenarios estimated increases in higher altitude locations (Pindus mountain range, Mount Olympus, etc.) and northeastern Greece, while those increases are anticipated to be higher during the RCP 4.5 scenario (7.0-10%). For the period of autumn (SON), both scenarios estimated slight decreases in the mean seasonal daily and annual GHI down to −5.0% in most parts of the Greek domain.

Clear-Sky Fraction-fclear -and Number of Clear-Sky Days
Aerosols and clouds are the key components in the atmosphere that reflect, transmit, or trap solar radiation and play a key role in the understanding of climate change. According to [16], in cloudy conditions, the high temporal and spatial fluctuation of cloudiness is the prime cause for the assessment of solar radiation. Also, changes in cloud coverage affect the diurnal temperature (DTR) variation, especially the maximum temperature during summer months and the minimum temperatures throughout winter, as [35] stated. However, cloud cover is a factor that cannot be smoothly simulated in climate models, as it provides a high level of uncertainty in the climate change predictions [36].
Under the prospect of this study, the mean annual fclear for the historical period 1980-2004, along with the anticipated differences between the future (RCPs 4.5/8.5) and historical scenarios, are presented in Figure 7. The WRF model for the historical period , as presented in Figure  7a, indicates lower mean fclear values in the high-altitude places by the order of 0.6-0.7, while marine and coastline locations show values above 0.85. As fclear is equivalent to the sunshine fraction, future changes between RCPs 4.5 and 8.5 are estimated to be in the range of ±4.5%, with the most increases up to 4.0% found in eastern and north continental Greece, as well in Euboea, Sporades, Rhodes, and Crete, while decreases down to 3.0% are projected in the western Greek mainland and in the Aegean Sea.  [33], which shows that the sunshine duration in most coastline and marine areas of Greece ranges  [33], which shows that the sunshine duration in most coastline and marine areas of Greece ranges above 50 and 60%, which corresponds to clear-sky days above 170 and 210 days and 2800 and 3000 h of sunshine, respectively.
Atmosphere 2020, 11, x FOR PEER REVIEW 13 of 17 above 50 and 60%, which corresponds to clear-sky days above 170 and 210 days and 2800 and 3000 h of sunshine, respectively.  Atmosphere 2020, 11, x FOR PEER REVIEW 13 of 17 above 50 and 60%, which corresponds to clear-sky days above 170 and 210 days and 2800 and 3000 h of sunshine, respectively.  On average, the differences in annual mean clear-sky days between the historic climatology are found to be in the range of ±20 days, which corresponds to percentage changes of about ±5.0% for both examined RCPs scenarios (Figure 8), which are consistent with the percentage differences in the annual mean f clear fraction (Figure 7). Overall, both scenarios indicated slight increases in the clear-sky days per year for most parts of continental Greece, except for western Greece, where decreases were found under the RCP 4.5. Results are in agreement with those of [17], who found that over western Greece, lower cloud coverage is expected, despite the fact that RCM cannot efficiently project the cloud cover conditions [38].

Discussion and Conclusions
This study focused on the evaluation of the anticipated changes in mean daily, annual, and seasonal GHI and in the development of a high-resolution solar climate atlas for Greece under climate change, with the use of a very high-resolution-5 × 5 km 2 -regional climate model. Thus, a new solar atlas was produced, which can provide useful information to the country's climate change adaptation policy, stakeholders' long-term planning, and RET community activities under the growing demand for more efficient solar technologies. This work was a first step towards obtaining a high-resolution solar climate atlas in Greece, with some limitations due to the limited number of formally available observational data.
The analysis was conducted for the historical period 1980-2004 and the future period 2020-2044, under the RCPs 4.5 and 8.5 emissions scenarios that represent intermediate and worst-case scenarios, respectively. In general, the geographical pattern and the different climatic zones beyond the Greek domain justify the variances across continental Greece and marine areas or coastlines at the same latitude. In Greece, noteworthy changes have already been indicated across continental country, coastal, or marine regions, even located in similar latitudes, affecting the cloud temporal and spatial variability, as a matter of the different topographical characteristics and of the diverse climate pattern over the Greek domain [16]. The WRF model for the historical period  presented spatial distribution of the GHI and demonstrated significant agreement with similar works [15][16][17]34], with values above 1600 kWh/m 2 in the broader area of central and northern Greece, and above 1800 kWh/m 2 in the broader area of the southern marine Greek parts of the Ionian and Aegean Sea. However, lack of adequate average solar potential was calculated at the locations of the highest altitudes, commonly found at the mountain ridges.
On average, analysis of the mean annual GHI variations showed values of similar magnitude and pattern between the two future scenarios, while differences among future and historical scenarios were found to be in the range of ±100 kWh/m 2 , which is translated to a percentage range of ±5.5%. Moreover, seasonal analysis indicated increases in the mean seasonal GHI values for the eastern Greek mainland and in the northwestern Aegean Sea, as well as in Euboea and Crete, which were up to 11% for both scenarios (RCPs 4.5 and 8.5). Slight decreases in the mean seasonal daily and annual GHI down to −5.0%, under the future scenarios, were projected for most parts of the Greek domain.
As the impact of clouds decrease the GHI and thus the energy generation of a solar system, the f clear and the mean annual and the seasonal clear-sky days have been studied, indicating that the differences in annual mean clear-sky days between the historic and future climatology could be in the range of ±4.5 and ±5.0%, respectively. This finding helps to anticipate mainly increases in the fraction of time, with no significant clouds blocking the sun in the eastern and northeastern parts of the Greek mainland under both future scenarios. In addition, increases in clear-sky days under RCP 8.5 and decreases for the western parts of Greece under RCP 4.5 were found, which are in agreement with the results of [17].
Concluding, the study indicates that GHI might be affected from climate change in the future. The annual changes were found to be small and in the order of ±5.0%, while the seasonal changes were shown to be in the range of ±10%. As the PV and solar systems' cumulative capacities are increased further and further across Greece, future research should additionally focus on the climate effects of the projected changes in GHI and temperature in the efficiency of a broad number of PV and solar systems. To meet the growing demand in solar energy applications, the installation of more accurate monitoring systems could help to validate more efficiently the future projection changes obtained with models.