Aerosol Direct Radiative Effects under Cloud-Free Conditions over Highly-Polluted Areas in Europe and Mediterranean: A Ten-Years Analysis (2007–2016)

: This study investigates changes in aerosol radiative effects on two highly urbanized regions across the Euro-Mediterranean basin with respect to a natural desert region as Sahara over a decade through space-based lidar observations. The research is based on the monthly-averaged vertically-resolved aerosol optical depth (AOD) atmospheric proﬁles along a 1 ◦ × 1 ◦ horizontal grid, obtained from the Cloud-Aerosol Lidar with Orthogonal Polarization (CALIOP) instrument measurements aboard the Cloud-Aerosol lidar and Infrared Pathﬁnder Satellite Observation (CALIPSO). To assess the variability of the anthropogenic aerosols on climate, we compared the aerosol vertical proﬁle observations to a one-dimensional radiative transfer model in two metropolitan climate sensible hot-spots in Europe, namely the Po Valley and Benelux, to investigate the variability of the aerosol radiative effects and heating rate over ten years. The same analysis is carried out as reference on the Sahara desert region, considered subject just to natural local emission. Our ﬁndings show the efﬁcacy of emission reduction policies implemented at government level in strongly urbanized regions. The total atmospheric column aerosol load reduction (not observed in Sahara desert region) in Po Valley and Benelux can be associated with: (i) an increase of the energy ﬂux at the surface via direct effects conﬁrmed also by long term surface temperature observations, (ii) a general decrease of the atmospheric column heating rate, and likely (iii) an increase in surface temperatures during a ten-year period. Summarizing, the analysis, based on the decade 2007–2016, clearly show an increase of solar irradiation under cloud-free conditions at the surface of +3.6 % and +16.6% for the Po Valley and Benelux, respectively, and a reduction of − 9.0% for the Sahara Desert.


Introduction
Europe, the so-called old continent, includes also the Western part of the Eurasia. Four times smaller than America and Asia, it hosts 11% of the World population, including some high densely populated areas. Marked by cultural diversity and origins, some of European regions are historically related to social and economical progress, e.g., the industrial revolution in the XVII century that led to an unavoidable increase in anthropogenic emissions. In the last decades, air-quality and climate related research activities identified some European regions as climate hot-spots. In these areas, resulting climate changes happen faster then respect to other areas [1][2][3]. Among them, we can count the Mediterranean basin, the cradle of western civilization that developed its culture and historical background with the contribution of many diverse populations living and moving across the region, which are linked together by traditions and history. This cultural development was dramatically influenced over the last few centuries by the Mediterranean climate, which pushed local populations to develop different agricultural and farming techniques with respect to other European populations. The unique local biological heritage of this large region has a very high economic value. However, the region is also extremely sensitive to climate change [4] with some heavily industrialized areas, like the Po-Valley and Catalan region. Coastal cities also suffer relevant pollution episodes due to marine vessels [5]. On other hand, the Benelux region (Belgium, The Netherlands and Luxembourg) is considered to be another climate hotspot because it is a densely populated area, home of one of the major industrial district in Europe [6,7]. In Benelux, the recorded levels of particulate matter (PM) pollution, based on the EU's air quality standards, makes the air quality "very poor" which is the second-worst level on the air quality index, just before the "extremely poor" quality category.
The aerosol emissions in those two regions, anthropogenic and/or natural, besides deteriorating the air-quality [8], can locally warm or cool the atmospheric column, both directly through their effect in modulating incoming solar radiation and indirectly through their influence on lifetime of clouds and precipitations [9][10][11]. The aerosol species differently interact with incoming shortwave (SW) solar radiation and the outgoing longwave (LW) radiation emitted by the Earth surface. Shortwave reflecting aerosols, e.g., sulfates, marine aerosols, increase the planetary albedo cooling the surface [12], while black-carbon aerosol layers absorb the incoming SW radiation warming the surrounding air around them. Dust aerosol layers can either reflect the incoming SW radiation and trap the outgoing LW radiation (at a lesser extent). Both effects can alter regional atmospheric stability, vertical motions, and affect large-scale circulation and regional hydrological cycles, causing significant regional climate sensitivity [13]. The net aerosol radiative effect (RE) at ground level changed in the past few decades, especially due to anthropogenic contribution caused by the increased human industrial activities. The increase of pollutant concentrations in the atmosphere peaked in the 1970s/early 1980s followed by a substantial decrease in aerosols concentrations, resulting from mitigation policies adopted by European governments to reduce air pollution, is very likely the cause of the well documented change of surface solar radiation over the past half century known as "global dimming" (less solar radiation reaching the Earth surface), observed from the 1950s to the early 1980s. The global dimming is followed by the so-called "global brightening" (more solar radiation reaching the Earth surface) which is still happening [14,15]. Differently with respect to the global dimming, in the global brightening the clouds (which are indirectly affected by the aerosol presence) still played a dominant role on the variations of downward surface shortwave radiation (SSR) during the last few decades [16]. While the effects of greenhouse gases on climate are assessed with great accuracy, the radiative effects of the short-lived climate forcers, as those from aerosol emissions, show a much larger uncertainty. This is highlighted by the Fifth Assessment Report of the Intergovernmental Panel on Climate Change (IPCC) [17].
Surface air temperature, together with moisture, wind, and turbulent diffusion, is one of the most important meteorological parameters influencing particles dilution within the Mixing Layer Height (MLH [18]). The presence of different aerosol mixing may result in a net radiative heating/cooling of the lower atmosphere, thus promoting/inhibiting the vertical mixing and causing more/less air pollutants being transported to upper levels as a result of the increasing/decreasing the boundary layer height [19,20].
The IPCC estimates that (i) absorbing and non-absorbing aerosols have very different impacts on shortwave surface radiative cooling, and that (ii) their sizes also play a crucial role in radiative balance. A first species of fine mode aerosols are those originating from sulphates and nitrates. This kind of aerosols are non-absorbing smaller fine-mode particles that are strongly reflecting incoming solar radiation and then with a net effect of increasing the Earth-atmosphere albedo. Other fine-mode aerosol species are constituted by organic matter (e.g., organic carbon) and by strongly-absorbing soot (e.g., black carbon). Regarding coarse-mode aerosol species, we find non-absorbing sea salt and weakly-absorbing mineral dust [21]. As a net effect, the shortwave cooling exerted by both fine and coarse-mode aerosols, e.g., sulphates, nitrates, sea salt, organic matter, and dust over oceans is on average outweighting the net warming effects by soot and dust over continents [22]. As stated by Kinne [22], if the aerosol fine-mode is mainly considered as the net anthropogenic contributor to RE, investigations of anthropogenic impacts can be limited to the cooling effect exerted by non-absorbing or weakly-absorbing species and to the warming by strongly-absorbing ones.
In quantitative terms, as stated by the latest IPCC report, the total surface anthropogenic radiative effect for 2011 relative to 1750 is +2.29 [+1.13 to +3.33] W/m 2 , with a stronger acceleration starting from the 1970s compared to prior decades. The total anthropogenic RE best-estimate for 2011 is 43% higher with respect to year 2005 figures. This is caused by two combined effects: (i) continuous increase in greenhouse gas concentrations and (ii) an increase of aerosol concentrations at the surface causing a weaker net cooling effect (negative RE). Thus, atmospheric aerosol RE, which includes cloud adjustments, is overall −0.9 (−1.9 to −0.1) W/m 2 (medium confidence), which is the net result from a cooling effect from most anthropogenic aerosol species and a warming effect from black carbon solar radiation absorption.
In this study, we examine how the anthropogenic emissions influence the aerosol radiative effect at the surface in three different regions: two densely-populated European regions, i.e., Northern Italy and Benelux, and a desert region, i.e., the Bodele depression in Africa. The latter being chosen as reference because uninfluenced by human activity. Our objective is to corroborate links between atmospheric turbidity, aerosol RE at the surface and surface temperature trends. Other studies quantitatively highlighted the aerosol radiative effects over Europe [23] and their correlation with the surface temperature [24]. Nevertheless, none of these studies focused on anthropogenic aerosol emissions on regions that are critical for climate change comparing the results with a region without significant human footprints over a decade (2007)(2008)(2009)(2010)(2011)(2012)(2013)(2014)(2015)(2016). In particular, this study is based on computing the trends in atmospheric turbidity and direct radiative balance at the surface, either for anthropogenic (i.e., Po Valley and Benelux, anthropogenic aerosols as black-carbon) or natural (i.e., desert area, dust) scenarios. Vertically-resolved (CALIOP, Cloud-Aerosol Lidar with Orthogonal Polarization) and columnar (AERONET, AErosol RObotic NETwork, Level 2.0) measurements of aerosols optical depth and a radiative transfer model (Fu-Liou-Gu) [25][26][27] are used to quantify the variability over time of the vertically-resolved profile of atmospheric aerosol load and their radiative effect. In addition, surface temperature trends (i.e., before and after the 1980), calculated for the regions of interest, are compared with atmospheric turbidity observations over a ten-year period (i.e., 2007-2016). There is high confidence that short-lived climate effects exerted by aerosols and their interactions with clouds have offset a substantial portion of global mean radiative effect from well-mixed greenhouse gases. Still, the aerosols continue to represent the largest uncertainty to total RE estimates.

Regions of Interest
The European Environment Agency (EEA) [28] supports implementation policies to limit regional pollutants and respect air quality standards. EU legislation contributes to the definition of air pollution policies and long-term development strategies to improve air quality in Europe. The Po-Valley and Benelux are notoriously highly-populated and industrialized areas in Europe, prone to poor surface air quality [28]. On other hand, the selected area over northern Africa (basically centered over the Bodele Depression) is recognized as one of most important sources of natural aerosols [29]. Consequently, these three different regions were used to evaluate and confirm our hypothesis of an existing link between human activities and the surface temperature variability in the Euro-Mediterranean region. We investigate trends in atmospheric turbidity over a ten-year period and correlate aerosol radiative effect characteristics and inter-annual variability with regional surface temperatures.

Surface Air Temperature: Trend Calculation
Long-term analysis of Climate Research Unit (CRU) TS 4.02 surface temperature records derived at 0.5 latitude-longitude degree resolution by CRU [30] was produced for two periods: the first from 1951 to 1980, and the second from 1981 to 2016. The different boxes, representative of the three areas of interest (Po Valley, Benelux, and Sahara desert), are defined as shown in Figure 1. Anomalies of monthly mean temperature were calculated for the period from 2007 to 2016, the same period covered by CALIPSO measurements used for this study. Figure 2 illustrates the inter-annual anomalies relative to the ten-year mean (continuous lines), linear trends (dashed lines) computed before and after the 1980, and the period of CALIPSO observations (light-green area) considered for the comparisons reported in the Section 3. The ten-year linear trend was estimated applying the Theil-Sen method [31,32] estimating its significance through the Mann-Kendall non-parametric test [33]. We split the trend calculations into two different periods-before and after the 1980-to best isolate the relatively rapid temperature increase in the last four decades. As pointed out by Figure 2, the trends have similar slopes for the most recent period. On the other hand, between 1950 and 1980, linear trends are inconsistent: the temperature behavior in the Po Valley indicates a decrease, and an increase in Benelux and desert area. This can be partially explained by the implementation of the major industrial plants in Northern Italy, e.g., automotive, steelmakers. Instead, for Benelux, the analysis needs further investigation.

Aerosol Optical Depth: CALIPSO and AERONET Measurements
The standard Cloud-Aerosol Lidar with Orthogonal Polarization (CALIOP) data products are classified by processing level according to the definitions given by NASA's Earth Observing System (EOS) project. The CALIOP Version 4.10 Level 1 data consist of profiles of attenuated backscatter coefficients at 1064 nm and 532 nm, while the Level 2 data report such quantities as layer location, layer type, and a number of derived optical parameters, including optical depths.
In this paper, we aggregate the CALIOP lidar observations of the atmospheric verticallyresolved optical depth as monthly mean profiles over a uniform spatial grid of 1 • × 1 • up to an altitude of 10 km. The quality of the data is obtained through the process described by Campbell et al. [34] that replicates very closely the method described by Winker et al. [35]. The latter methodology is formally endorsed by the CALIPSO (Cloud-Aerosol Lidar and Infrared Pathfinder Satellite Observations) project team development. Nevertheless, the screening metrics in Campbell et al. [34] help to improve the retrieval method robustness for this investigation purposes. First, if the presence of a cloud is initially detected, the profile is not taken into consideration. This choice minimizes the relative attenuation of the lidar signal that can have an impact on the aerosol extinction coefficient retrieval. In the second instance, CALIOP profiles are not taken into consideration if the retrieval fails to resolve aerosol from surface up to 250 m. Lastly, based on Toth et al.'s [36] analysis, those profiles without aerosol loading (i.e., aerosol optical depth; AOD = 0) are disregarded. Applying those restrictions as a whole, it helps to limit the signal attenuation effects that could bias the averaged solutions of the aerosol extinction profiles derived for this study. An average of 1800 qualifying profiles were used each year to build the profiles used during our analysis. As final step, a Gaussian weighting function was used to solve the averages to increase the retrievals significance relative to the 1 • × 1 • study domain center. Moreover, the domain-average CALIOP 532 nm aerosol extinction profiles are built with a vertical resolution of 100 m that corresponds to each of the native Level 2 species resolved for use in a radiative transfer model experiment described in Section 2.4. In Reference [34], among others, the authors characterize CALIOP AOD measurements together with their accuracy. Eventually, lidar measurements provide critical information on vertical distribution of aerosol optical properties that are especially related to air quality, climatology and climatic effects. In addition, the column-integrated aerosol optical depth derived from Level 2 CALIOP 532nm observations is also widely used in combination with passive-based L2 aerosol retrievals for a comprehensive understanding of regional and global aerosol optical properties [37].
The NASA Aerosol Robotic Network (AERONET) [38], a federated network of homogeneous commercially-available sunphotomters, has collected long-term continuous and readily accessible public domain data of aerosol optical, microphysical and radiative prop-erties [38]. For this study, we consider the new Version 3 (V3) of Level 2 (Quality Assured) monthly AOD averages, a product directly downloadable at NASA's AERONET website.

The Fu-Liou-Gu Radiative Transfer Model
We assess the direct aerosol radiative effects and heating rate using the Fu-Liou-Gu (FLG) Radiative Transfer (RT) model [19,25,26,39,40]. Following the methodology reported in Reference [41], the FLG RT code was initialized using the vertically-resolved lidar CALIPSO extinction profiles (described above) of the different aerosol species, which were matched to the Optical Properties of Aerosol and Clouds (OPAC) catalog [42]-based physical and optical models written into the code [40]. The total number of the RT model vertical levels was adjusted to match the lidar extinction spatial resolution of 100 m from ground to 10 km (100 levels). The total aerosol radiative effect is computed adding all contributions (in terms of optical depth) from the different aerosol species identified by the CALIPSO classification algorithm in the considered region and matched with the corresponding FLG aerosol species, as shown in Reference [41]. The aerosols present in the atmosphere at highly urbanized sites and desert is very different. For this reason is very important to correctly chose the different aerosol types. From CALIOP data analysis on highly urbanized sites, among the eighteen aerosol species parameterized into the FLG RT model from the OPAC catalog, we considered only the main four species that match the CALIPSO classification: transported dust (CALIPSO: "dust"), urban (CALIPSO: "polluted continental/smoke"), 50% dust, and 50% urban (CALIPSO: "polluted dust") and black carbon (CALIPSO: "smoke"). Those aerosol species account for more than 98% of the total AOD value in each region. For the Sahara desert instead, dust is the unique aerosol species present in the atmosphere. A potential source of uncertainty in parameterizing Version 4 CALIOP aerosol typing into FLG RT is that the classification algorithm cannot distinguish between local urban pollution and an advected smoke aerosol layer from distant sources below 2.5 km. Even if it is reasonable to to assume that those episodes are not frequent, they can introduce a bias in the anlysis. To compute the radiative effects at the top of the atmosphere and at the surface, the FLG RT model, which account also for aerosol hygroscopicity, solves the radiative fluxes at each level for 18 spectral bands (12 short-wave, 6 long-wave) [25,26]. For each annually averaged lidar extinction profile used as input in the FLG code, the aerosol direct radiative effect (DRE) at the bottom of the atmosphere and the vertically-resolved heating rate (HR) is computed by subtracting the net radiative flux when the aerosols are present in the atmosphere from the net radiative flux obtained with pristine conditions, as shown in Equation (1): The ancillary data needed as input by the FLG model (i.e., the atmospheric thermodynamic variable profiles as the temperature, the atmospheric pressure, the mixing ratio and ozone concentration) are obtained from a mid-latitude standard atmosphere (USS976) for Po Valley and Benelux regions, while, for Sahara, we use the tropical standard atmosphere USS976 profiles, except for the relative humidity that was taken from the arctic atmosphere. The Solar Zenith Angle (SZA), used in FLG computations, is calculated for each box at noon local time of the 15th day of the month. While the anthropogenic aerosol contribution to LW radiative fluxes is negligible [19], dust aerosol layers (Sahara region) trap the outgoing longwave radiative fluxes (which strongly depend on water vapor). Nevertheless, the dust net radiative balance still remains negative at the surface. Following the approach reported in Reference [43], we applied a wavelength-independent constant albedo value of 0.12 for the two urban environments and 0.37 for the Sahara Desert (obtained integrating the hemispherical directional reflectance [44]), while the infrared surface emissivity is set to a constant value of 0.98.

Correlations between the Atmospheric Turbidity and Surface Air Temperature
Comparisons between long term time series of air temperature and aerosol optical depth (AOD) observations, as measured by AERONET and CALIOP, were performed to investigate trends throughout the considered period for this work (see Figure 1). CALIPSO measurements, covering the ten-year period (2007)(2008)(2009)(2010)(2011)(2012)(2013)(2014)(2015)(2016), were analyzed to understand temporal variability of aerosol vertical distribution. The CALIOP lidar AOD is retrieved inverting the lidar equation (Equation (2)): where P λ (r)r 2 is the range corrected received power, C λ is related to the system calibration coefficient, and T 2 λ (r) is the total two-way signal transmission containing the integration of α tot (r) the total extinction coefficient. β mol (r) and β aer (r) are the molecular and aerosol backscattering coefficients, respectively. The total two-way transmission can be expressed as (Equation (3))): while the AOD is defined as (Equation (4) where α aer is the aerosol extinction coefficient in km −1 , and r is the top of the atmosphere. Equation (2) has two unknowns (the extinction and backscattering aerosol coefficients) because the molecular optical properties can be retrieved using ancillary data, e.g., pressure and temperature profiles from the standard atmosphere. CALIPSO retrieval algorithm assumes a range independent ratio between the extinction and backscattering aerosol coefficient constant, the so-called lidar ratio (in sr): S shows a very high variability (20-120 sr), depending on the aerosol species [45]. The used values in CALIPSO retrieval can be found in Reference [46]. The authors decided to split, by arbitrary choice, the lower troposphere into three (3) levels to investigate the inter-annual variability of the atmospheric turbidity within the mixing layer height. For this study, we consider both the total AOD, and the AOD calculated at different vertical intervals from 40 m above the ground up to 2 km as follows: level 1 from 40 m to 500 m, level 2 from 500 m to 1 km, and level 3 from 1 km to 2 km. Together with AOD and air temperature, the aerosol radiative effects at the bottom-of-the-atmosphere (BOA) were computed by using the FLG Radiative Transfer Model (see Section 2.4).
In order to provide a quantitative assessment of the behavior of such parameters on the polluted regions subject to a strong human activity, as well as in natural desert zones, trend calculations were carried out based on the monthly averages of the lidar data aggregated from the original CALIPSO Level 2 products previously described. Table 1 reports the percentage variation of the computed radiative fluxes estimated per decade of AOD and BOA during 2007-2016 as inferred by CALIPSO observations through the Sen method. According to Toth et al. [47], this method is preferable when data are noisy as its output is a nonparametric estimate of the slope after computing all points dataset median of pairwise slopes. Nevertheless, our analysis shows some limitations because data are really noisy and often the correlation is not significant (p > 0.1). The same analysis on temperature records (Section 2.2), which, due to the higher temporal resolution, are much less noisy is much more statistical significant. On other hand, analysis on CALIOP data still gives indication on the observed trend during the decade, confirming the negative trend shown in Europe by Toth et al. [47]. During the analyzed decade, the AOD trend and, consequently, the total atmospheric turbidity for the Po Valley and Benelux is negative (i.e., −3.4% and −16.5%, respectively), which is associated with an increase of the incoming solar irradiation at the surface under cloud-free conditions (+3.6% and +16.6%, respectively). Because more solar radiation is reaching the ground, as a consequence, a net warming effect of the surface air masses is likely.
The authors speculate that high concentrations of fine and/or absorbing particles originated from combustion processes might be one of the most important causes of surface temperature increase. In fact, highly urbanized/industrialized areas, such as Po Valley and Benelux, are characterized by important sources of black carbon [48]. It is well-known that such type of particles absorb the incoming solar radiation influencing the vertical heating rate close to the major emission sources, while, far from emission areas, it can contribute to the surface temperature increase as a diabatic heat insulator. For instance, in 2004, Johnson et al. [49] demonstrated the warming effect of absorbing particles within the mixing layer. Recently, Sand et al. [50] found a non-linear relationship between emissions rate and temperature variation at ground level: the higher is the BC emission rate, the lower is the air temperature sensitivity because the convection in the lower atmospheric levels promote convection and turbulent motion moving BC from the surface to the mixing layer top. In other words, the net motion results in enhancing aerosol upward vertical transport. Such conditions favor the pollutant dilution within the atmospheric mixing layer resulting in lower AOD values at lower levels (e.g., level 1 and level 2 of Table 1). But, on the other hand, the positive aerosol load increase above 1km found in the Po-Valley might be also ascribed to the northern Africa mineral dust transport [51]. Dust outbreaks in Po-Valley occur with a frequency that spans between 10 and 15 events per year [51]. Over the Benelux, the impact of desert mineral dust is significantly lower as it is inversely related with latitude [52]. In fact, the AOD values decrease overall in the considered layers [53]. This translates into higher solar energy reaching the surface as suggested by the FLG computed BOA. This is consistent with the more rapid air temperature increase over the considered decade in Benelux compared to Po Valley (see Figure 2).
Over the desert regions, our results highlight an increase of the total AOD. A thicker aerosol dust layer over the years causes a reduction of the solar irradiance reaching the surface. This effect is partially compensated by the outgoing longwave radiation trapped by the dust layer, but still the net balance at the surface is negative. Such analysis points out that the net solar radiation reaching the ground decreased during the last decade, but the surface air temperature (see Figure 2) increased: this net warming effect is a consequence of the greenhouse gases concentration increase over the decade. Figure 3 illustrates the long-term temperature measurement record (1980-2020) and AOD both from space (CALIOP) and from ground measurements (AERONET; see Figure 3) for the Po Valley (upper panel), Benelux (middle), and Sahara (lower). Starting from the monthly averages, monthly air temperature deviations from the ten-year mean are also reported, coupled with the monthly averages of both CALIOP and AERONET observations. Looking at the surface temperature behavior, we find an unequivocal increase for the three considered areas. However, this increase has different causes depending on the considered region as previously described. The temperature trends put in evidence about 0.4 degrees/10 yr increment for urbanized areas (i.e., 0.43 degrees/10 yr in Benelux and 0.37 degrees/10 yr in Po Valley) and about 0.3 degrees/10 yr (i.e., 0.28 degrees/10 yr) for the desert region .
From lidar vertically-resolved observations, it is clear how changes in human activity differently shaped the surface and vertically-resolved AOD observations during the analyzed decade in the Po Valley and Benelux with respect to the desert Bodele depression. The AOD reduction is significant for the high-urbanized areas, on the other hand, as stated above, an increase of atmospheric turbidity is observed over the desert region (see Table 1). Such discrepancies are even more evident if we consider the trend calculated from CALIOP measurements. As reported in Table 1, AOD decrement is observed in both the Po Valley (up to 1 km) and Benelux (overall the column); on the contrary, the AOD enhances (overall the column) for Sahara.

Conclusions
Through our study, we finally give quantitative evidence to the changes in direct radiative effect and heating rate under cloud-free conditions associated with anthropogenic aerosol emissions over a decade (2007)(2008)(2009)(2010)(2011)(2012)(2013)(2014)(2015)(2016) in two highly-urbanized Euro-Mediterranean areas, i.e., the Po-Valley in Southern Europe and the Benelux in central Europe. We compared our findings with those of a wild and remote region not affected by human activity as the Sahara Desert, taken as reference. Through this study, using interpolated surface temperature measurements and remote sensing observations of the Aerosol Optical Depth (AOD) both vertically-resolved (from lidar) and columnar (from AERONET sunphotometers), for the first time, we quantitatively assess how much the implementation of emission mitigation policies reversed the so-called "aerosol dimming" over the polluted regions.
This research was carried out taking into account long term temperature anomalies, the satellite observations of atmospheric turbidity and the calculation of direct radiative effects. In higher urbanized regions, we considered four main aerosol species (through CALIPSO lidar observations) accounting for 98% the total Aerosol Optical Depth (AOD), dust (natural source), urban (anthropogenic source), black carbon (anthropogenic source), and polluted dust (anthropogenic source). As expected, in the desert region, the aerosol is 100% dust (natural source).
The aerosol direct radiative effects were estimated through the Fu-Liou-Gu Radiative Transfer Model using the CALIOP extinction profiles as input. The analysis carried on the Po-Valley and Benelux regions highlights a reduction of aerosol concentration in the atmosphere after year 2000 that yields an increase of solar radiation at the surface, which is compatible with the observed AOD yearly reduction rate of 0.005-0.01 yr −1 .
Results can be summarized as follows: (i) in the Po Valley and Benelux, two highly urbanized/industrialized regions (where anthropogenic aerosols are dominant) the columnar aerosol loading and the surface air temperature trends are inversely related (i.e., consistent with a reduction in solar dimming), while (ii) in the Sahara region (where the anthropogenic pressure is absent and crustal aerosol is dominant) the atmospheric turbidity does not appear to be temporally correlated with changes in surface air temperature.
These results reveal that the aerosol loading should be taken into consideration when assessing the temporal variability of the surface radiative flux budget. In particular, the concomitant increment of atmospheric transparency, and the enhancement of atmospheric greenhouse gases concentrations, can likely be considered as the two most significant effects that promote surface temperature increase. With respect to these considerations, additional quantitative investigations are required in future experiments and studies. A further appropriate assessment of spatial-temporal variations requires deeper investigations, as suggested by the differences between the highly populated areas and desert areas pointed out by this study.
More in-depth studies are also required to investigate the real effect of the solar irradiance (i.e., +3.6% in Po Valley and +16.6% in Benelux) in increasing the air temperature at the surface, considering the concurrent effects due to the presence of absorbing and reflecting aerosol layers at higher altitudes.