Long-Term Patterns and Trends of Shortwave Global Irradiance over the Euro-Mediterranean Region

: The spatiotemporal patterns and trends of shortwave global irradiance (SWGI) are a crucial factor affecting not only the climate but also sectors of the economy. In this work, the ERA5-Land reanalysis dataset is employed and evaluated against in situ measurements from a dense network of surface stations operated by the National Observatory of Athens over Greece, revealing a good agreement between the two datasets. Then, the spatiotemporal variability of SWGI is investigated over the Euro-Mediterranean region (10 ◦ W–42 ◦ E and 30 ◦ N–52 ◦ N) for a 40-year period (1981–2020). SWGI exhibits a smooth latitudinal variability from north to south of − 5.4 W/m 2 /degree on an annual scale, while it varies signiﬁcantly on a seasonal basis and is almost four times lower in the winter than in the summer. The SWGI trend during the analyzed period was found to be positive and statistically signiﬁcant at the 95% conﬁdence level. Spring and summer are the periods where positive and the strongest rates of SWGI trends are evident, while in the winter and autumn, negative or neutral trends were found. The increasing SWGI trend shows a slowdown during the beginning of the 2000s in all seasons, except autumn. The SWGI trend decreases by about − 0.06 W/m 2 /decade every 100 m of elevation increase.


Introduction
Shortwave global irradiance (SWGI), the irradiance from direct sunlight and diffuse sky radiation that reaches the Earth's surface (with a wavelength range 300-3000 nm), plays a primary role in the Earth's climate system, as it is a key component of the surface energy balance. It has significant implications for hydrological components (e.g., evaporation), plant photosynthesis, carbon cycle, cloud formation, climate, and climate change [1][2][3][4]. In addition, SWGI constitutes a renewable source of energy; thus, the accurate analysis of its trends is critical for the siting and operation of renewable energy plants [5]. Hence, in recent decades, high priority has been placed on understanding the changes in SWGI, as they have profound environmental, societal, and economic implications.
In order to better understand SWGI, it is important to have as many observations available as possible. Such observations may come from either in situ measurements or satellite platforms, which have better coverage than ground-based observations. These two methods may complete each other so as to provide improved datasets. On the other hand, atmospheric modeling is a widely used alternative since SWGI is a main component of numerical weather prediction models. The main drawback of these model-based datasets is that they cannot always reproduce local conditions in different locations.
On a seasonal scale, Manara et al. [10] pointed out that the brightening period is particularly pronounced in autumn, while in winter, the seasonal trend is negative (not statistically significant) in northwest Italy during the period of 1990-2016. Furthermore, Chiacchio and Wild [7] and Sanchez-Lorentzo et al. [22] showed a pronounced increasing tendency in spring over Europe. Additionally, Kambezidis et al. [22] found an increasing SWGI trend over the Mediterranean, which is the largest in spring. On the other hand, Perdigao et al. [14] and Bais et al. [6] reported that the maximum positive SWGI trend is observed in summer. The seasonal SWGI trends found in previous studies are presented in Table 1. These interannual and seasonal variations are a complex issue that has been attributed to either direct and indirect aerosol effects or to cloud fluctuations [13]. Despite the fact that the SWGI climatology is a well-documented field of research, only a few studies have investigated the temporal evolution of SWGI beyond 2000 [16][17][18]23,24] and especially during the 2010s over Europe. These studies have reported that there is an indication of stabilization or a slowdown of the upward trend in SWGI during the beginning of the 2000s [6,16,18,21,25,26]. Thus, there is lack of detailed knowledge regarding SWGI trends over the last decade, which motivated the present study. The present study has a twofold objective. First, we investigate the ability of the ERA5-Land database to represent the actual SWGI measured by ground-based stations and to reproduce local conditions in different areas (e.g., continental, seaside, islands, etc.) over the Greek territory.
Second, we provide a long-term analysis of the patterns and trends of SWGI for the Euro-Mediterranean region, which compared to the existing literature, (a) has a more recent time span, as it covers the last decade, (b) provides a long term overview of SWGI (1981-2020), and (c) is of enhanced spatial resolution since it is based on ERA5-Land reanalysis dataset with a 9 km spatial resolution, compared to the 31 km of ERA5 or 80 km spatial resolutions of the ERA-Interim data sets.
The paper is structured as follows: Section 2 provides a detailed presentation of the methodology followed and the datasets used, Section 3 discusses the results, and finally, Section 4 presents the conclusions and the future prospects of this study.

Materials and Methods
In this paper, data form ground measurements (for the period 2010-2019) and data from ERA5-Land re-Analysis (for the period 1981-2020) were used. First, these data were used to validate the ability of ERA5-Land to represent the local conditions of SWGI. Secondly, analysis regarding (a) annual and seasonal spatiotemporal distribution, (b) annual and seasonal trends, and (c) the validation of SWGI with elevation were performed based on ERA5-Land. The workflow chart of the analysis conducted in this paper is presented in Figure 1, while details about the datasets and methodology used are given in next Sections 2.1-2.3.

Ground-Based Measurements
Observations provided by surface networks of meteorological stations are necessary for model validation and verification (e.g., Kotroni and Lagouvardos [28]; Akylas et al. [29]). Figure 2 shows the locations of the surface weather stations that were used to evaluate the ERA5-Land dataset. Indeed, from the~100 surface weather stations that are equipped with solar sensors, we selected those that have the longest time-series available (as the network has grown during the last 15 years), had no data gaps, presented no drift, and were distributed across Greece so as to cover different geographical regions from north to south, with various geomorphological characteristics (coastal, continental and insular areas). Hourly data are provided by the dense network of automated surface meteorological stations operated by the METEO unit at the National Observatory of Athens [30]. The validation period spans 2010-2019. The evaluation is restricted to shortwave global irradiance as measured by the selected ground stations. The weather station type that was used was the Davis Vantage Pro2, whichmeasures ambient air temperature, relative humidity, wind speed and direction, rainfall, atmospheric pressure, and solar and UV radiation. It also measures indoor air temperature and humidity at the location of the station's console/datalogger and calculates a suite of bioclimatic indices and derived meteorological parameters. Table 3 provides the technical characteristics of the various sensors incorporated within each station [31,32]. Table 3. Technical characteristics of the weather station used in NOAAN [30].

ERA5-Land Re-Analysis Dataset
The ERA5-Land is a global atmospheric reanalysis dataset produced by the European Centre for Medium-Range Weather Forecast's (ECMWF) and has data data available from 1981 [33], providing a consistent view of the evolution of land variables. It has an enhanced resolution of 0.1 • × 0.1 • compared to ERA5 climate reanalysis (31 km resolution), while the temporal frequency of the model output is hourly. ERA5-Land has been produced by implementing a series of improvements with respect to ERA5 (representing the land component). The shortwave radiative variables were computed by the Fouquart and Bonnel [33] scheme using a six-band shortwave parameterization model. The variable used in this study was the "surface solar radiation downwards" from the ERA5-Land dataset, which is available to download from https://cds.climate.copernicus.eu/ (accessed on 25 October 2021) [33]. The surface solar radiation downwards represents the shortwave irradiation (in units of J/m 2 ) that is accumulated over a time period of 24 h. The daily mean SWGI values (in units of W/m 2 ) were calculated by dividing surface solar radiation downwards by 86,400 s (i.e., the length of a day in seconds).

Methods
The comparison between the ERA5-Land data and the measurements from the meteorological stations was performed by using common statistical indices, such as bias, absolute bias, and root mean square error (RMSE). The grid (model) points were the nearest neighbours to the stations. The descriptions of those indices are provided below: • Bias or mean error (ME) of forecasted-model values: • Absolute bias or mean absolute error (MAE): the absolute mean value of the above differences. • Root mean square error (RMSE): • Mean percentage error (MPE): the computed average of mean errors by which forecasts of a model differ from actual values of the quantity being forecast.
• Mean absolute percentage error (MAPE): the computed average of mean absolute errors by which forecasts of a model differ from actual values of the quantity being forecast.
Here obs(i) denotes the recorded (observed) value by each meteorological station at time i for(i) the respective ERA5-Land value of the one nearest to the meteorological station grid point and k the size of the sample (the number of observations that are available for each station).
For the climatology and trend analysis of the SWGI, a much larger area than the area of validation was used that covered the wider Mediterranean and Southern Europe, extending from 10 W to 42 E and from 30 N to 52 N. The annual and seasonal mean values were computed using the ERA5-Land mean daily values of SWGI for the period from 1981 to 2020 in a grid resolution of 0.1 • × 0.1 • . Annual and seasonal anomalies were calculated as the difference between the corresponding value and the climatic average values by considering the period 1981-2020 as the reference period. The seasons were defined according to the World Meteorological Organization (WMO) i.e., winter (December-January-February), spring (March-April-May), summer (June-July-August,) and autumn (September-October-November).
The slope of the trends was estimated using the Theil-Sen method [33,34], and the statistical significance was assessed using the Mann-Kendall non-parametric test [34,35] at the 95% confidence level. The slope (b) of a trend in sample data could be estimated as follows: where x t and x t * are the sequential data series at the interval t and t * (t > t * ). The Mann-Kendall test confirms the existence of a positive or negative trend for a given confidence level and has been widely used in previous studies (e.g., [2,36]), as it has a lot of advantages. The Mann-Kendall test does not require that the data are normally distributed, and it is not affected by missing data, by outliers, or by the length of the time series. Specifically, for a time series containing a set of observations (x i , x j ), the Mann-Kendall test is given by where The value of Mann-Kendall test Z MK is computed by Where q is the number of tied groups, and t p is the number of observations in the pth group. The value of Z MK (Equation (7)) determines if there is a monotonic trend at the time-series. A monotonic trend, which is significant at the 95% confident level, is accepted if |Z MK | ≥ 1.96. A positive value in Equation (6) indicates that there is a positive trend in the observations and vice versa.

Validation of ERA5-Land SWGI
In order to evaluate the ability of ERA5-Land to represent SWGI, radiation measurements from meteorological stations located in different geographical areas of Greece were employed, as described in Section 2.1.
In the frame of the present work, the comparisons were performed on a seasonal basis and are presented in Figure 3, Figure 4, Figure 5, Figure 6. During the winter, ME ranged from 5-38 W/m 2 , MAE ranged from 27 to 47 W/m 2 , and RMSE ranged from 37-79 W/m 2 ; MPE ranged from 11-40%, and MAPE ranged from 30-42% In the spring, the deviations increased for ME (23-57 W/m 2 ), for MAE (54-78 W/m 2 ), and for RMSE (74-123 W/m 2 ) and decreased for MPE (4-22%) and MAPE (21%-35%). During summer, the statistical scores were similar to the ones for spring, and more specifically, ME ranged from 7 to 56 W/m 2 , MAE ranged from 39 to 74 W/m 2 , and RMSE ranged from 51 to 124 W/m 2 , but MPE reduced to 0-19%, and MAPE reduced to 14-26%. Finally, during autumn, ME values were between 1 and 37 W/m 2 , MAE values ranged from 34 to 55 W/m 2 , RMSE ranged from 45 to 96 W/m 2 , MPE ranged from −5% to 18%, and MAPE ranged from 20% to 37%. In Sianturi et al. [37], the ERA5 (31 km horizontal resolution) and MERRA2 (approximately 50 km horizontal resolution) reanalysis data were evaluated against in situ measurements. The results showed that for ERA5, the model overestimates SWGI, which is in agreement with our results. The aforementioned study focused over an area of islands in Indonesia (which is closer to the Earth's equator than Greece, so larger absolute deviations are to be expected), and the statistical indices for MAE revealed higher deviations (50-96 W/m 2 ) than our study did for the Greek islands (49-79 W/m 2 ). The values of the statistical measures of the present study are close to the ones presented by Perdigao et al. [15], who compared the downscaled output of the Weather Research and Forecasting model (WRF) using NCEP/NCAR reanalysis data at a 5 km horizontal resolution against ground station measurements for the areas of Portugal and Spain.    During the spring, the errors are higher and are distributed more evenly around Greece ( Figure 4). This result is in agreement with Perdigao et al. [15], where the statistics for the area of the Iberian Peninsula reveal higher deviations during spring and summer. In the present analysis, during the summer months, higher differences are found in the mainland stations, while at the coastal and island areas, they are lower ( Figure 5). This fact may be attributed to the atmospheric thermal instability during summer, which mainly occurs over the Greek peninsula at around noon, causing increased cloud cover compared to other areas. This localized summertime instability cannot be easily reproduced by the model at the available resolution (9 km × 9 km). In autumn, MAEs and RMSEs are lower than they are in spring, but they are also evenly distributed around the area of Greece ( Figure 6). It should be noted that as expected, during spring and summer, the values of SWGI are higher because of the longer duration of the daylight as well as the lower cloud cover duration. Finally, we have to note that the stations that are close to large cities such as Athens are much more affected by smog, scattering and absorbing more aerosols than other stations on the islands.
In addition to the above statistics, we also present the scatter plots (with the relevant trendlines) for three representative stations at different areas of Greece, one in the northern part (Notio Pedio Kozani), the second in the central Greece (Spata), and the third inthe southeast (island of Rhodes) (Figure 7). There is a good correlation between observational and ERA5-Land points, which confirms the good quality of the model data.
Here, we should note that the results from the presented validation suffer from the fact that in situ observations have been compared to gridded model data that are representative of an area that is 10 km × 10 km. Observational and model time-series are both valuable depending on their availability, quality, homogeneity, and representativeness. Keeping in mind the limitations of each source of data, in the following, a climatological analysis of SWGI is presented based on the ERA5-Land reanalysis data, which is a very recent data set with higher resolution than previously available reanalysis provided in the literature. Figure 8 shows the spatial distribution of multiyear mean SWGI (in W/m 2 ) in the extended Mediterranean area from 1981 to 2020 acquired from ERA5-Land reanalysis data. The multiyear mean of SWGI was computed from mean daily values in a grid resolution of 0.1 • × 0.1 • . The multiyear mean SWGI shows a latitudinal gradient from north to south, ranging from 100 W/m 2 in the northern regions to 260 W/m 2 in the southern regions, corresponding to a spatial average value of 177 W/m 2 . The annual latitudinal gradient from south to north is −5.4 W/m 2 per degree ( Figure 9).   The mean values of SWGI are in line with previous studies over Europe. Specifically, Alexandri et al. [17] studied the eastern Mediterranean region using satellite data from 1983 to 2013 and found mean SWGI values greater than 250 W/m 2 in southern Africa, while in Greece and Turkey, the mean SWGI ranged from 150 to 220 W/m 2 . Perdigao et al. [14,15], using data obtained from the WRF model for the period of 1950-2010 and reanalysis data from ERA-40, found SWGI ranging from 150 to 250 W/m 2 over the Iberian Peninsula. The calculated mean values in the present study are consistent with those obtained from the analysis of observational data throughout Europe. Specifically, Sanchez-Lorenzo et al. [19] analysed data from 47 surface stations over the entire European region from 1983 to 2005 and found that the mean annual SWGI ranged from 87 to 174 W/m 2 , while Ohmura [13], who used observational data for the surface station in Davos, Switzerland, reported that the annual SWGI in the period from 1990-2005 ranged from 140 to 180 W/m 2 , which compares well with the respective values in the same region in our analysis, which range from 150 to 175 during the same period (not shown).

Annual and Seasonal Average Distribution
The seasonal spatial distribution of SWGI (in W/m 2 ) are depicted in Figure 10. As it was expected, the seasonal mean of SWGI reaches the minimum values in winter (mean spatial average: 89 W/m 2 ) and the maximum values in summer (mean spatial average: 265 W/m 2 ). The month with the lowest SWGI values is December (mean spatial average: 70 W/m 2 ), while the months with the highest SWGI values are June and July (mean spatial average: 275 and 273 W/m 2 , respectively) due to the high solar elevation and decrease in cloudiness over the Mediterranean [18]. The interannual variability of the mean SWGI in spring and summer are similar to the annual series, as these seasons contribute 68% to the annual mean SWGI, which is in line with other studies (e.g., Alexandri et al. [18]; Perdigao et al. [15]).
In the winter, the mean SWGI ranges from 28 W/m 2 in North and central Europe to 180 W/m 2 in African regions. Spring, as a transition period to the maximum mean SWGI values of the summer, exhibits increased mean SWGI values that range from 144 W/m 2 in the United Kingdom to 303 W/m 2 in Algeria. In summer, the mean SWGI values reach their maxima in all regions, showing an increase of 200% compared to winter. In autumn, the mean SWGI is reduced by 50% compared to summer. The mean SWGI seasonal values are in line with other studies provided in the literature focusing on the European region. Specifically, Alexandri et al. [17] found average values of SWGI over the eastern Mediterranean of 100.3 W/m 2 in winter, 230.5 W/m 2 in spring, 309.6 W/m 2 in summer, and 168.9 W/m 2 in autumn within a period cover 1983-2013. Additionally, they attributed the highest/lowest values of the mean SWGI over the southern/northern parts of the eastern Mediterranean to the lower/higher latitudes and cloudiness. Furthermore, Perdigao et al. [15] studied the SWGI seasonality over the Iberian Peninsula and found the average values of 75-125 W/m 2 in winter,~230 W/m 2 in spring,~320 W/m 2 in summer, and~180 W/m 2 in autumn in the period from 1950 to 2010. The interannual variability of SWGI, in all seasons, shows a latitudinal gradient from north to south as the spatial distribution of the multiyear mean of SWGI ( Figure 10). The latitudinal gradient is −5.6, −5.2, −4.9, and −6.1 W/m 2 per degree towards the North in winter, spring, summer, and autumn, respectively (Figure 9).

Annual and Seasonal SWGI Trend
The interannual variability of SWGI has been linked to changes in aerosols (e.g., Silveira et al. [38]; Kazadzis et al. [9]; Nabat et al. [11]; Folini and Wild [8]), in the cloud cover (e.g., Sanchez-Lorenzo et al. [22]; Chiacchio and Wild [7]) as well as to changes in the cloud types and optical cloud properties (Hatzianastassiou et al. [39]). Aerosols influence solar radiation in two ways: directly, by reflecting incoming sunlight (e.g., by sulfate and sulfuric acid particles), and indirectly, by changing the number of cloud condensation nuclei particles. Indeed, as aerosols act as cloud condensation nuclei particles, the increasing aerosol load increase the average number and size of the cloud droplets and hence, affect the scattering properties and cloud albedo [12]. Thus, the reason for the interannual variation of SWGI is a complex problem, which is still subject to ongoing research (e.g., Sfica et al. [40]; Chiacchio and Wild [7]). Figure 11 shows the annual trend of SWGI in the extended Mediterranean area along with the statistical significance at the 95% level, which was estimated over the period from 1981-2020 from the ERA5-Land data. Overall, the annual SWGI trend ranges from −0.9 to 4.5 W/m 2 per decade (mean regional value 1.7 W/m 2 per decade, as shown in Table 4) and is characterized by positive values almost throughout the entire study area. It must be highlighted that regions near the Mediterranean Sea, such as Greece, Croatia, Albania, Turkey, Italy, Spain, and northern African countries, for which higher values of mean SWGI were depicted (Figure 8), are associated with lower values of positive annual trend values ( Figure 11). In contrast, in Central Europe, the positive annual SWGI trend is higher and ranges from 1.8 to 4.5 W/m 2 per decade, except in Great Britain, which exhibits lower trend values that are comparable to those in the southeastern Mediterranean. The highest trend values are observed in the northeastern part of the study area (region between 46 • N to 53 • N and 26 • E to 40 • E). Almost neutral to negative trend values (but not statistically significant) are observed in some limited areas over Greece, Turkey, and in Algeria. The spatial distribution of the annual SWGI trend is consistent with previous studies performed in various regions of Europe, e.g., in Spain [14], in Switzerland and Germany [26], and in Central Europe [22], although the magnitude of these trends is dependent on both the study area and the analysed period. Table 2 summarizes some of the results of previous studies dealing with the SWGI trends over various regions in Europe.  The increasing trend of SWGI is in line with the reduction of sulfate emissions and other anthropogenic aerosols since 1980 [41][42][43] and the observed decrease of aerosol optical depth (AOD) reported in many studies (e.g., Founda et al. [44]; Nabat et al. [45]; Zerefos et al. [46]; Ruckstul et al. [27]; Mishchenko et al. [47]). The spatial distribution of SWGI trends reveals that the brightening period presents geographical variations. The Euro-Mediterranean "hot spot" occurs in the northeastern part of the study area, which is in accordance with other studies (e.g., Sanchez-Lorenzo et al. [22]; Kambezidis et al. [18]; Chiacchio and Wild [7]). Indeed, in Eastern Europe, anthropogenic emissions were strongly reduced in the late 1980s [41], while reductions in Western Europe started in the 1980s. Additionally, Chiacchio et al. [48] found that the annual change in sulfate AOD for Europe was −69%, while Western Europe had a decline of 63%, and Eastern Europe had a maximum of 75% for the period from 1985-2007. Moreover, Sfica et al. [40] showed that the cloud cover decreased over Central Europe during the period of 1981-2014 and that the maximum decline occurred in the east of Europe.
It is worth mentioning that the southern part of Europe has lower SWGI trends compared to Central Europe, which could be attributed to the higher AOD values found in southern Europe [6,49,50]. Kazadzis et al. [9] found that the brightening trend of Athens (Greece) is lower than the reported SWGI trends around Europe for the period from 1980-2012. Furthermore, Kambezidis et al. [18] reported higher rates of SWGI increase over the western Mediterranean in comparison to the eastern part for the period of 1979-2012. Figure 12 shows the spatial distribution of the seasonal SWGI trends during the period of 1981-2020. In winter, the SWGI trends range from −4.0 to 1.8 W/m 2 per decade, while negative trend values prevail in the area of study, except in some sub regions south of 42 • N and over Great Britain. However, the SWGI trend is not statistically significant in most of the region (87% of the whole area). During this season, the highest negative trend is calculated, which is in agreement with Sanchez-Lorenzo et al. [16]. In spring, the sign of the SWGI trend becomes positive, ranging from 0 to 7.4 W/m 2 per decade, with the exception of small regions in Algeria and Spain, where the trend is slightly negative. The SWGI trend is higher in Central Europe, with the maxima occurring in the area north of the Black Sea and on the border between France and Germany. Increased positive trends are found in the summer in the regions of Central Europe where the maximum annual values range from 2.4 to 10.1 W/m 2 per decade (in the continental European region north of 49 N). In areas such as Great Britain, Cyprus, Israel, and part of the Algeria-Morocco region, slightly negative trend values are observed, but they are not statistically significant. In autumn, a more neutral to positive SWGI trend prevails in most of the study area (70% of the area). Slightly negative values are found mainly in the Balkan Peninsula and in the Algeria-Morocco region. It should be noted that this trend is statistically significant in very few areas (10% of the whole area) in autumn. Overall, the seasonal variability is in line with previous studies reported by Sanchez-Lorenzo et al. [16,22] and Chiacchio and Wild [7], who found that spring and summer are the seasons with positive and the strongest rates of SWGI trends, while in winter, there is a tendency for negative or neutral trends, and in autumn, the trends are weaker, with a tendency for negative rates over the European region. Table 1 summarizes the results of previous studies dealing with the seasonal SWGI trends over various regions in Europe. In previous studies, the strongest increase in the SWGI trend values were reported in spring [7,16,22], in autumn [10], and in summer [6,14], while in the present study, the maximum positive trend is found in summer. The magnitude of the trends found in previous studies is dependent on various parameters such as the used dataset, the selected studied region, the length of the analysed period, and the selected reference period used to compute the anomaly. Nevertheless, the maximum SWGI trend found in the present study is observed in regions where the maximum decline of aerosols in summer has been reported in previous studies (e.g., Filonchyk et al. [51]; Parding et al. [52]).
Chiacchio et al. [48] found that the temporal decrease in AOD correlates well with the increase in spring and summer mean SWGI, but it does not explain the trends found in SWGI during winter and autumn. Chiacchio and Wild [7] proposed that during winter and autumn, the North Atlantic Oscillation (NAO) and the associated cloud cover could be contributing to the decadal changes in SWGI. Indeed, partial warming in Europe in the last few decades has been attributed to a positive NAO trend [53]. Furthermore, Sfica et al. [40] showed that in the warm season (from May to September), the occurrence and intensity of the Azores High had slightly decreasing ridges during the period of 1981-2014, resulting in a reduction in the cloud cover over Eastern Europe. During the cold season (from November to March), the increasing trend of an anticyclonic circulation over Eastern Europe induced the blocking of the Mediterranean cyclones along their track over the Balkan Peninsula, resulting in less cloudiness over the eastern part of Europe and enhanced cloudiness over the southern part of the continent. Indeed, in autumn, the contrasting trend between the east and west Mediterranean is pronounced (Figure 12d). Figures 13 and 14 show the mean annual and seasonal series of the SWGI, which have been regionally averaged for the entire study area, together with a smoothed line obtained with a 5 year Gaussian low-pass filter and with the trend lines, which have been calculated for the periods 1981-2000 and 2001-2020. As it was expected, the annual SWGI time series show positive trends (brightening) during 1981-2020, with a statistically significant increase of +1.7 W/m 2 per decade. Interestingly, the smoothed time series of SWGI shows a local maximum near 2000, followed by a small decline between 2000 and 2010. Indeed, the SWGI trend over 1981-2000 is almost twice as large as the SWGI trend from 2001-2020 (Table 4). This is in line with previous studies, which reported indications of less distinct solar brightening after 2000 compared to the 1990s [6,21,26] or stabilizing SWGI [16,18,25,54]. The transition of solar dimming to brightening in Europe during the 1980s was mainly attributed to the sharp decrease of aerosols (e.g., Streets et al. [43]). The largest decrease in the concentrations of anthropogenic aerosols occurred before 2000 (e.g., Stjern et al. [55]; Ruckstuhl et al. [27]); thus, the AOD and its impact on solar radiation could not have changed dramatically after 2000. Indeed, Wild et al. [23]), who compared sunshine duration data with the average surface solar radiation from 2000 to 2005, argued that decreases in cloudiness rather than in aerosol concentration have primarily affecting the SWGI trend in Europe since 2000.
In winter, the SWGI trend is close to zero (0.5 W/m 2 per decade), with an indication for a slight overall decrease; however, it is non-significant at the 95% confidence level. The spring and summer SWGI time series shows a strong and significant upward SWGI trend throughout the whole period (+2.5 and +3.6 W/m 2 per decade, respectively) and shows a similar evolution to the annual time series. Interestingly, SWGI trends over the period from1981-2000 for spring and summer are almost 4 and 2 times greater than the SWGI trends over 2001-2020 (Table 4), respectively. In both seasons, the mean SWGI temporal evolution (filtered time series) shows a peak near 2000 followed by a period with little change. In autumn, the SWGI series show a slight overall increase of +0.6 W/m 2 per decade, which is significant at the 95% confidence level.

Variation of SWGI Trend with Elevation
In developed regions such as Europe, trends in aerosol concentrations are influenced by anthropogenic emissions, which mostly originate from areas close to sea level. As a consequence, AOD trends and their associated SWGI trends are lower in mountainous areas than they are in low-elevated areas (e.g., Stjern et al. [55]; Dong et al. [56]). Conversely, the higher frequency of thick orographic clouds over mountains shields a significant fraction of the solar radiation, contributing to local minima at mountainous areas [17]. For the present analysis, the SWGI trend values were paired with the gridded elevation fields (Figure 15). The resulting number of pairs was divided into 25 bins (with equal number of pairs), and for each bin, the mean SWGI trend was estimated. The box and whisker plots presented in Figure 15 provide information for the distribution of the median, lower (10% and 25%), and upper (75% and 90%) percentiles of the SWGI trend for various bins of the elevation. The annual SWGI trend decreases by about 0.06 W/m 2 per decade for every 100 m of increase in elevation. It is worth mentioning that the evolution of the SWGI trend shows a slightly increasing trend for those grid-points located at elevations lower than 200 m and a decreasing trend for grid-points located at elevations higher than 200 m. A similar behavior is observed on seasonal scales, especially in spring and summer, where the variation of the SWGI trend with elevation shows that the trend decreases by about −0.06 (spring) and −0.21 W/m 2 per decade (summer) for every 100 m of increase in elevation ( Figure 16). In winter, the trend in SWGI increases slightly, (+0.03 W/m 2 per decade every 100 m), while in autumn, the trend decreases slightly (−0.01 W/m 2 per decade every 100 m). All of the aforementioned trends are statistically significant at the 95% level. These findings are in line with Manara et al. [10], Zeng et al. [57], and Philipona [58], who also found an overall decreasing SWGI trend with elevation.

Discussion and Concluding Remarks
In this paper, a climatological and trend analysis of SWGI in the Euro-Mediterranean area (10 W-42 E and 30 N-52 N) was performed. For this purpose, we used the ERA5-Land reanalysis dataset from the period of 1981-2020 with an enhanced resolution of 0.1 • × 0.1 • . The comparison of this dataset with observational data provided by a network of surface weather stations in Greece showed a relatively good ability to represent the SWGI in the ERA5-Land dataset, with MAE ranging from 30-50 W/m 2 during the autumn and winter periods and rising to 50-70 W/m 2 in spring and summer. These values appear to be more accurate than those presented in other similar studies (which have used ERA5 data) for other areas (e.g., Indonesia in Sianturi et al. [29], Iberian Peninsula in Perdigao et al. [30]).
The suitability of using the ERA-5 land data set for this climatological analysis is further supported by the high correlations between this dataset and the observational dataset. Further, in order to support the suitability of the ERA-5 land dataset for the trend analysis provided in this work, we note that the comparison of trends between the observational and model data sets showed the same slopes and similar trend values for the validation period (not shown). Thus, there is confidence that the trend slopes discussed in this analysis are correct in terms of sign, while the values might deviate from those of an observational network. Although the validation area is restricted over Greece for the purposes of this study, we assume that the statistics are similar over the much larger area used for the climatology and trend analysis. This is supported by the fact that the statistics over Greece are comparable to those provided in the literature over other areas of the studied domain. Nevertheless, the aforementioned limitations should be taken into account when exploring the results of the present work.
Analysis of the annual patterns of SWGI over the area revealed a latitudinal variability from north to south of −5.4 W/m 2 per degree. The mean spatial average value of SWGI is 177 W/m 2 . The seasonal mean of SWGI reaches the minimum values in winter (89 W/m 2 ), while spring and summer contribute 68% to the annual mean SWGI. The maximum mean SWGI values in summer found in the present study are consistent with previous SWGI studies SWGI (e.g., Alexandri et al. [18]; Perdigao et al. [15]) and AOD values reported in previous studies (e.g., Filonchyk et al. [51]; Parding et al. [52]). Indeed, Parding et al. [52] linked the highest AODs values in summer to the high relative humidity and temperature, which contribute to the hygroscopic growth of thin aerosols [59] and strengthen photochemical reactions, leading to the formation of secondary aerosols [60]. In general, high AOD values in spring and summer are associated with the intensification of anthropogenic aerosol emissions from crop biomass burning [61], forest fires [52], and dust transported from remote sources (e.g., Saharan dust) [62,63].
The SWGI trend analysis revealed an overall positive trend (brightening) for the period of 1981-2020, with a mean regional value of 1.7 W/m 2 per decade. The spatial analysis revealed that the brightening period presents geographical variations. The highest SWGI trend appears in Central Europe, especially in the northeastern part of the study area, with values reaching up to 4.5 W/m 2 per decade. Almost neutral to negative trend values are observed in some limited regions of Greece, Turkey, and in Algeria, although these trends are not statistically significant. The spatial pattern of the SWGI trend is strongly related to the reductions of anthropogenic aerosol emissions in the 1980s in Europe (e.g., Nabat et al. [17]; Chiacchio et al. [17]) and the decreasing cloud cover over the inner European continent (e.g., Sfica et al. [17]; Mateos et al. [3]; Hatzianastassiou et al. [3]).
On a seasonal level, spring and summer are the seasons with the highest positive SWGI trend, while in winter, there is a tendency for negative or neutral trends, and in autumn, the trends are slightly positive or neutral over the Euro-Mediterranean region. Previous studies have indicated that the brightening in spring and summer is possibly associated to the temporal evolution of AOD (e.g., Chiacchio et al. [3]), while the trends in winter and autumn could be attributed to cloud cover variation (Chiacchio and Wild [7]). In this study, the largest positive SWGI trends were found in summer and not in spring or autumn, as reported in previous studies [7,10,16,22]. Nevertheless, it must be highlighted that the magnitude of the trend is strongly dependent on various parameters such as the length of the analysed period and the selected reference period used to compute the anomaly. Moreover, the maximum SWGI trend found in the present study is in line with the maximum decline of aerosols in summer in previous studies [51,52]. Interestingly, the increasing SWGI trend shows a slowdown during the beginning of the 2000s in all seasons, except autumn, confirming the reported indications of previous studies (e.g., Kambezidis et al. [7]; Wild et al. [23]). This slowdown in the increasing SWGI trend has been attributed to the cloud cover beyond 2000 [23] and not to aerosol variability [6]. The analyses of the dependence of trends in SWGI on elevation revealed that the SWGI trend decreases about −0.06 W/m 2 per decade (p value ≤ 0.05) for every 100 m of increase in elevation. This finding is in line with previous studies that have found a decreasing AOD trend with elevation (e.g., Dong et al. [56]). This elevational gradient is negative in all seasons, except in autumn. The largest negative gradient is found in summer.
The present study provides an update on SWGI trends for Europe by extending previous analyses beyond the year 2010 and using reanalysis datasets with a better resolution than those used in the past. Interestingly, it has been shown that the solar dimming masked the greenhouse effect up until 1980s, while the subsequent brightening contributed to global warming (Wild et al. [64]). Future investigation should seek to disentangle the contribution of cloudiness and aerosol emissions to the variability of SWGI. Additionally, the contribution of each type of cloud (high, mid, and low) to the variations of SWGI should be investigated. This would lead to a better understanding of the complex interactions between aerosols, cloud cover, and irradiance and their influence on climate change.
Funding: This research received no external funding.
Institutional Review Board Statement: Not applicable.

Data Availability Statement:
The raw ERA5-Land data required to reproduce these findings are available to download from https://cds.climate.copernicus.eu/ (accessed on 25 October 2021). The raw data of the observations and all of the datasets generated during the current study are available from the corresponding author upon reasonable request.