Ambient air particulate total lung deposited surface area (LDSA) levels in urban Europe

This study is supported by the RI-URBANS project (Research Infrastructures Services Reinforcing Air Quality Monitoring Capacities in European Urban & amp; Industrial Areas, European Union's Horizon 2020 research and innovation program, Green Deal, European Commission, contract 101036245). This study is also supported by National Natural Science Foundation of China (42101470, 72242106) and in part by the Chunhui Project Foundation of the Education Department of China under Grant HZKY20220053. This study benefited from the Aerosol, Clouds and Trace Gases Research Infrastructure (ACTRIS), especially the so-called ACTRIS-2 H2020 research project (grant no 654109), and the authors would like to thank ACTRIS (The Aerosol, Clouds and Trace Gases Research Infrastructure), especially the ACTRIS in situ EBAS Data Centre (EBAS), for providing datasets to the study. This study is also partly funded by the National Institute for Health Research (NIHR) Health Protection Research Unit in Environmental Exposures and Health, a partnership between UK Health Security Agency (UKHSA) and Imperial College London, and the UK Natural Environment Re-search Council, and the views expressed are those of the author(s) and not necessarily those of the NIHR, UKHSA or the Department of Health and Social Care. The research was also supported by the Hungarian Research, Development and Innovation Office (grant no. K132254). We thank also the support from "Agencia Estatal de Investigacion" from the Spanish Ministry of Science and Innovation, and FEDER funds under the projects CAIAC (PID2019-108990RB-I00); and the Generalitat de Catalunya (AGAUR 2017 SGR41) and the Direccio General de Territori. IMT Nord Europe and LOA acknowledge financial support from the Labex CaPPA project, funded by the French National Research Agency (ANR-11-LABX-0005-01), and the CLIMIBIO and ECRIN projects, both financed by the Regional Council "Hauts-de-France" and the European Regional Development Fund (ERDF).This study aims to picture the phenomenology of urban ambient total lung deposited surface area (LDSA) (including head/throat (HA), tracheobronchial (TB), and alveolar (ALV) regions) based on multiple path particle dosimetry (MPPD) model during 2017-2019 period collected from urban background (UB, n = 15), traffic (TR, n = 6), suburban background (SUB, n = 4), and regional background (RB, n = 1) monitoring sites in Europe (25) and USA (1). Briefly, the spatial-temporal distribution characteristics of the deposition of LDSA, including diel, weekly, and seasonal pat-terns, were analyzed. Then, the relationship between LDSA and other air quality metrics at each monitoring site was investigated. The result showed that the peak concentrations of LDSA at UB and TR sites are commonly observed in the morning (06:00-8:00 UTC) and late evening (19:00-22:00 UTC), coinciding with traffic rush hours, biomass burning, and atmospheric stagnation periods. The only LDSA night-time peaks are observed on weekends. Due to the variability of emission sources and meteorology, the seasonal variability of the LDSA concentration revealed sig-nificant differences (p = 0.01) between the four seasons at all monitoring sites. Meanwhile, the correlations of LDSA with other pollutant metrics suggested that Aitken and accumulation mode particles play a significant role in the total LDSA concentration. The results also indicated that the main proportion of total LDSA is attributed to the ALV fraction (50 %), followed by the TB (34 %) and HA (16 %). Overall, this study provides valuable information of LDSA as a predictor in epidemiological studies and for the first time presenting total LDSA in a variety of European urban environments.RI-URBANS project (Research Infrastructures Services Reinforcing Air Quality Monitoring Capacities in European Urban amp; Industrial Areas, European Union's Horizon 2020 research and innovation program, Green Deal, European Commission, 101036245)National Natural Science Foundation of China (NSFC)Chunhui Project Foundation of the Education Department of ChinaAerosol, Clouds and Trace Gases Research Infrastructure (ACTRIS)National Institute for Health Research (NIHR) Health Protection Research Unit in Environmental Exposures and HealthUK Research & Innovation (UKRI) Natural Environment Research Council (NERC)National Research, Development & Innovation Office (NRDIO) - Hungary"Agencia Estatal de Investigacion" from the Spanish Ministry of Science and InnovationGeneralitat de Catalunya 42101470, 72242106Direccio General de Territori HZKY20220053Agence Nationale de la Recherche (ANR) 654109Regional Council "Hauts-de-France"European Union (EU) K132254, PID2019-108990RB-I00, AGAUR 2017 SGR41, ANR-11-LABX-0005-01ERD

Dynamics of the Atmospheric Boundary Layer over two middle-latitude rural sites with Doppler lidar

The Atmospheric Boundary Layer (ABL) over two middle-latitude rural sites was characterized in terms of mean horizontal wind and turbulence sources using a standard classification methodology based on Doppler lidar. The first location was an irrigated olive orchard in Úbeda (Southern Spain), representing one of the most important crops in the Mediterranean basin and a typical site with Mediterranean climate. The second location was PolWET peatland site in Rzecin (Northwestern Poland), representing one of the largest natural terrestrial carbon storages that have a strong interaction with the climate system. The results showed typical situations for non cloud-topped ABL cases, where ABL is fully developed during daytime due to convection, with high turbulent activity and strong positive skewness indicating frequent and powerful updrafts. The cloud-topped cases showed the strong influence that clouds can have on ABL development, preventing it to reach the same maximum height and introducing top-down movements as an important contribution to mixing. The statistical analysis of turbulent sources allowed for finding a common diurnal cycle for convective mixing at both sites, but nocturnal wind shear driven turbulence with marked differences in its vertical distribution. This analysis demonstrates the Doppler lidar measurements and the classification algorithm strong potential to characterize the dynamics of ABL in its full extent and with high temporal resolution. Moreover, some recommendations for future improvement of the classification algorithm were provided on the basis of the experience gained.Fundacion Ramon ArecesEuropean Space Agency 4000119961/16/NL/FF/mgPolish National Science Centre (NCN) 2021/40/C/ST10/00023Spanish Government CGL2015-73250-JIN CGL201681092-R CGL2017-83538-C3-1-R CGL2017-90884-REDT PID2020117825GB-C21 PID2020-120015RB-100Andalusian Regional Government P18-RT-3820FEDER-UGR program ARNM-430-UGR20University of GranadaACTRIS-2 Research Infrastructure Project of the European Union's Horizon 2020 research and innovation program 654109European Cooperation in Science and Technology (COST) CA18235Universidad de Granada/CBU

Impact of primary NO2 emissions at different urban sites exceeding the European NO2 standard limit

A large part of the European population is still exposed to ambient nitrogen dioxide (NO2) levels exceeding the European Union (EU) air quality standards, being a key challenge to reduce NO2 concentrations across many European urban areas, particularly close to roads. In this work, a trend analysis of pollutants involved in NO2 processes was done for the period 2003–2014 in traffic sites fromthree Spanish cities (Barcelona,Madrid and Granada) that still exceed the European NO2 air quality standard limits. We also estimated the contributions of primary NO2 emissions and photo-chemically formed NO2 to the observed ambient NO2 concentrations in order to explore their possible role in the observedNO2 concentration trends. TheNOx andNOconcentrations at these traffic sites showed significant decreasing trends during the period 2003–2014, especially at Barcelona (BARTR) andMadrid (MADTR) traffic stations. The NO2 concentrations showed statistically significant downward trends at BARTR and MADTR and remained unchanged at Granada traffic station (GRATR) during the study period. Despite the significant decrease in NO2 concentrations in BCNTR and MADTR during the analysed period, the NO2 concentrations observed over these sites still above the annual NO2 standard limit of 40 μg m−3 and, therefore, more efficient measures are still needed. Primary NO2 emissions significantly influence NO2 concentrations at the three analysed sites. However, as no drastic changes are expected in the after-exhaust treatment technology that can reduce primary NO2 emissions to zero in the near future, only a substantial reduction in NOx emissions will help to comply with the NO2 European air quality standards. Reduction of 78%, 56% and 16% on NOx emissions in Barcelona,Madrid and Granada were estimated to be necessary to comply with the NO2 annual limit of 40 μg m−3

Aerosol number fluxes and concentrations over a southern European urban area

This work was supported by the Spanish Ministry of Economy and Competitiveness through projects PID2020-120015RB-100, CGL201681092-R, and CGL2017-90884-REDT, by the Andalusia Regional Government through project P18-RT-3820 and P20-00136, by the European Union's Horizon 2020 research and innovation program through project ACTRIS-2 (grant agreement No 654109). This research was partially supported by Project RTI2018.101154.A.I00 funded by MCIN/AEI/10.13039/501100011033/FEDER "Una manera de hacer Europa". The authors thank the Parque de la Ciencias for making this research possible. Juan Andr ' es Casquero-Vera is supported by BES-2017-080015 funded by MCIN/AEI/10.13039/501100011033 and FSE "El FSE invierte en tu futuro". Funding for open access charge: Universidad de Granada/CBUA.Although cities are an important source of aerosol particles, aerosol number flux measurements over urban areas are scarce. These measurements are however important as they can allow us to identify the different sources/sinks of aerosol particles and quantify their emission contributions. Therefore, they can help us to understand the aerosol impacts on human health and climate, and to design effective mitigation strategies through the reduction of urban aerosol emissions. In this work we analyze the aerosol number concentrations and fluxes for particles with diameters larger than 2.5 nm measured by eddy covariance technique at an urban area (Granada city, Spain) from November 2016 to April 2018. This is the first study of particle number flux in an urban area in the Iberian Peninsula and is one of the few current studies that report long-term aerosol number flux measurements. The results suggest that, on average, Granada urban area acted as a net source for atmospheric aerosol particles with median particle number flux of 150 x 10(6) m(-2) s(-1). Downward negative fluxes were observed in only 12% of the analyzed data, and most of them were observed during high aerosol load conditions. Both aerosol number fluxes and concentrations were maximum in winter and 50% larger than those measured in summer due to the increased emissions from domestic heating, burning of residual agricultural waste in the agricultural area surrounding the site, as well as to the lower aerosol dilution effects during winter. The analysis of the seasonal diurnal variability of the aerosol number concentration revealed the significant impact of traffic emissions on aerosol population over Granada urban area in all seasons. It also shows the impact of domestic heating and agricultural waste burning emissions in winter as well as the influence of new particle formation processes in summer and spring seasons. Closer analysis by wind sector demonstrated that both aerosol concentrations and fluxes from urban sector (where high density of anthropogenic sources is located) were lower than those from rural sector (which includes agricultural area but also the main highway of the city). This evidences the strong impact of aerosol emissions from traffic circulating on the highway on aerosol population over our measurement site.Spanish Government PID2020-120015RB-100 CGL201681092-R CGL2017-90884-REDTAndalusia Regional Government P18-RT-3820 P20-00136European Commission 654109MCIN/AEI/FEDER "Una manera de hacer Europa" RTI2018.101154.A.I00FSE "El FSE invierte en tu futuro"Universidad de Granada/CBUA MCIN/AEI BES-2017-08001

Lidar and Radar Signal Simulation: Stability Assessment of the Aerosol–Cloud Interaction Index

This work was supported by the Spanish Ministry of Economy and Competitiveness through projects CGL2016-81092-R, PID2020-120015RB-I00 and RTI2018-101154-A-I00, the Regional Government of Andalusia through project AEROPRE (P18-RT-3820), and by the Spanish Ministry of Education, Culture and Sports and Spanish Ministry of universities through grant FPU19/05340. The financial support for EARLINET in the ACTRIS.IMP 871115 (H2020-INFRADEV-2018-2020) is gratefully acknowledged. This work is related to activities within the COST Action CA18235 PROBE (PROfiling the atmospheric Boundary layer at European scale). The authors thank the University of Granada, Programa Operativo FEDER Andalucia 2014-2020 through project DEM3TRIOS (A-RNM-430-UGR20). Juan Antonio Bravo-Aranda received funding from the Marie Skodowska-Curie Action Cofund 2016 EU project-Athenea3i under grant agreement no. 754446. Maria J. Granados-Munoz project has received funding from the European Union's Horizon 2020 research and innovation programme under the Marie Skodowska-Curie grant agreement No 796539. The financial support for UGR and FEDER funds through project B-RNM-496-UGR18 is gratefully acknowledged.Aerosol-cloud interactions (ACI) are in the spotlight of atmospheric science since the limited knowledge about these processes produces large uncertainties in climate predictions. These interactions can be quantified by the aerosol-cloud interaction index (ACI index), which establishes a relationship between aerosol and cloud microphysics. The experimental determination of the ACI index through a synergistic combination of lidar and cloud radar is still quite challenging due to the difficulties in disentangling the aerosol influence on cloud formation from other processes and in retrieving aerosol-particle and cloud microphysics from remote sensing measurements. For a better understanding of the ACI and to evaluate the optimal experimental conditions for the measurement of these processes, a Lidar and Radar Signal Simulator (LARSS) is presented. LARSS simulate vertically-resolved lidar and cloud-radar signals during the formation process of a convective cloud, from the aerosol hygroscopic enhancement to the condensation droplet growth. Through LARSS simulations, it is observed a dependence of the ACI index with height, associated with the increase in number (ACINd) and effective radius (ACIreff) of the droplets with altitude. Furthermore, ACINd and ACIreff for several aerosol types (such as ammonium sulfate, biomass burning, and dust) are estimated using LARSS, presenting different values as a function of the aerosol model. Minimum ACINd values are obtained when the activation of new droplets stops, while ACIreff reaches its maximum values several meters above. These simulations are carried out considering standard atmospheric conditions, with a relative humidity of 30% at the surface, reaching the supersaturation of the air mass at 3500 m. To assess the stability of the ACI index, a sensitivity study using LARSS is performed. It is obtained that the dry modal aerosol radius presents a strong influence on the ACI index fluctuations of 18% cause an ACI variability of 30% while the updraft velocity within the cloud and the wet modal aerosol radius have a weaker impact. LARSS ACI index uncertainty is obtained through the Monte Carlo technique, obtaining ACIreff uncertainty below 16% for the uncertainty of all LARSS input parameters of 10%. Finally, a new ACI index is introduced in this study, called the remote-sensing ACI index (ACIRs), to simplify the quantification of the ACI processes with remote sensors. This new index presents a linear relationship with the ACIreff, which depends on the Angstrom exponent. The use of ACIRs to derive ACIreff presents the advantage that it is possible to quantify the aerosol-cloud interaction without the need to perform microphysical inversion retrievals, thus reducing the uncertainty sources.Spanish Government CGL2016-81092-R PID2020-120015RB-I00 RTI2018-101154-A-I00Junta de Andalucia P18-RT-3820Spanish Government FPU19/05340EARLINET in the ACTRIS.IMP 871115University of Granada, Programa Operativo FEDER Andalucia through project DEM3TRIOS A-RNM-430-UGR20European Commission 754446 796539UGREuropean Commission B-RNM-496-UGR1

Impact of urban aerosols on the cloud condensation activity using a clustering model

The indirect effect of aerosols on climate through aerosol-cloud-interactions is still highly uncertain and limits our ability to assess anthropogenic climate change. The foundation of this uncertainty is in the number of cloud condensation nuclei (CCN), which itself mainly stems from uncertainty in aerosol sources and how particles evolve to become effective CCN. We analyze particle number size distribution (PNSD) and CCN measurements from an urban site in a two-step method: (1) we use an unsupervised clustering model to classify the main aerosol categories and processes occurring in the urban atmosphere and (2) we explore the influence of the identified aerosol populations on the CCN properties. According to the physical properties of each cluster, its diurnal timing, and additional air quality parameters, the clusters are grouped into five main aerosol categories: nucleation, growth, traffic, aged traffic, and urban background. The results show that, despite aged traffic and urban background categories are those with lower total particle number concentrations (Ntot) these categories are the most efficient sources in terms of contribution to the overall CCN budget with activation fractions (AF) around 0.5 at 0.75%supersaturation (SS). By contrast, road traffic is an important aerosol sourcewith thehighest frequency of occurrence (32 %) and relatively high Ntot, however, its impact in the CCN activity is very limited likely due to lower particlemean diameter and hydrophobic chemical composition. Similarly, nucleation and growth categories, associated to new particle formation (NPF) events, present large Ntot with large frequency of occurrence (22%and 28%, respectively) but the CCN concentration for these categories is about half of the CCN concentration observed for the aged traffic category, which is associated with their small size. Overall, our results show that direct influence of traffic emissions on the CCN budget is limited, however, when these particles undergo ageing processes, they have a significant influence on the CCN concentrations and may be an important CCN source. Thus, aged traffic particles could be transported to other environments where clouds form, triggering a plausible indirect effect of traffic emissions on aerosol-cloud interactions and consequently contributing to climate change.BioCloud project - MCIN/AEI RTI2018.101154.A.I00FEDER "Unamanera de hacer Europa" European Commission 871115 ATMO_ACCESS 101008004Ministry of Science and Innovation, Spain (MICINN) Spanish Government PID2020-12001-5RB-I00 GL2016-81092-R CGL2017-90884REDTJunta de AndaluciaUGREuropean Commission B-RNM- 474-UGR18NIMBUS B- RNM-496UGR18Junta de Andalucia P2000136 AEROPRE P-18-RT-3820University of Granada Plan Propio PPVS2018-04 LS2022-1Spanish Government FPU19/05340Ministry of Science and Innovation, Spain (MICINN)Spanish Government PRE2019-09082

Spatial and temporal variability of carbonaceous aerosols: assessing the impact of biomass burning in the urban environment

Biomass burning (BB) is a significant source of atmospheric particles in many parts of the world. Whereas many studies have demonstrated the importance of BB emissions in central and northern Europe, especially in rural areas, its impact in urban air quality of southern European countries has been sparsely investigated. In this study, highly time resolved multi-wavelength absorption coefficients together with levoglucosan (BB tracer) mass concentrations were combined to apportion carbonaceous aerosol sources. The Aethalometer model takes advantage of the different spectral behaviour of BB and fossil fuel (FF) combustion aerosols. The model was found to be more sensitive to the assumed value of the aerosol Ångström exponent (AAE) for FF (AAEff) than to the AAE for BB (AAEbb). As result of various sensitivity tests the model was optimized with AAEff = 1.1 and AAEbb = 2. The Aethalometer model and levoglucosan tracer estimates were in good agreement. The Aethalometer model was further applied to data from three sites in Granada urban area to evaluate the spatial variation of CMff and CMbb (carbonaceous matter from FF or BB origin, respectively) concentrations within the city. The results showed that CMbb was lower in the city centre while it has an unexpected profound impact on the CM levels measured in the suburbs (about 40%). Analysis of BB tracers with respect to wind speed suggested that BB was dominated by sources outside the city, to the west in a rural area. Distinguishing whether it corresponds to agricultural waste burning or with biomass burning for domestic heating was not possible. This study also shows that although traffic restrictions measures contribute to reduce carbonaceous concentrations, the extent of the reduction is very local. Other sources such as BB, which can contribute to CM as much as traffic emissions, should be targeted to reduce air pollution.This research was partially supported by the Andalusia Regional Government through projects P10-RNM-6299 and P12-RNM-2409, by the Spanish Ministry of Economy and Competitiveness and FEDER through project CGL2013_45410-R; by EUREKA and the Slovenian Ministry of Economic Development and Technology grants (Eurostars grant E!4825 FC Aeth, JR-KROP grant 3211-11-000519); and by European Union's Horizon 2020 Research and Innovation Programme under grant agreement No. 654109, ACTRIS-2. The authors would like to thank Air Quality Service from Junta de Andalucía (Consejería de Medio Ambiente y Ordenación del Territorio) and Vicerrectorado de Política Científica e Investigación from the University of Granada for their support in the installation of the Aethalometer at PC and GV, respectively.G. Titos was partially funded by Programa del Plan Propio de Investigación “Contrato Puente” of the University of Granada and by the Spanish Ministry of Economy and Competitiveness under postdoctoral program Juan de la Cierva – Formación (FJCI-2014-20819)

Extinction-related Angström exponent characterization of submicrometric volume fraction in atmospheric aerosol particles

The AEAOD– ΔAEAOD grid proposed by Gobbi et al. (2007) is a graphical method used to visually represent the spectral characterization of aerosol optical depth (AOD), i.e. Angström exponent (AE) and its curvature, in order to infer the fine mode contribution (η) to the total AOD and the size of the fine mode aerosol particles. Perrone et al. (2014) applied this method for the wavelengths widely used in lidar measurements. However, in neither case does the method allow for a direct relationship between η and the fine mode fraction contribution to the total aerosol population. Some discussions are made regarding the effect of shape and composition to the classical AE-ΔAE plot. The potential use of particle backscatter measurements, widely used in aerosol characterization methods together with extinction measurements, is also discussed in the AE-ΔAE grid context. A modification is proposed that yields the submicron contribution to the total volume concentration by using particle extinction data, and a comparison to experimental measurements is made. Our results indicate that the use of a modified AE-ΔAE grid plot to directly obtain submicrometric and micrometric mode fraction to the total aerosol population is feasible if a volume-based bimodal particle size distribution is used instead of a number-based one.Andalusia Regional Government through project P12-RNM-2409Spanish Ministry of Sciences, Innovation and Universities (CGL2016-81092 and CGL2017 -90884 - REDT

Cloud condensation nuclei activation properties of Mediterranean pollen types considering organic chemical composition and surface tension effects

Supplementary data to this article can be found online at https://doi.org/10.1016/j.atmosenv.2023.119961This work was supported by BioCloud project (RTI2018.101154.A. I00) funded by MCIN/AEI/10.13039/501100011033, FEDER “Una manera de hacer Europa” and NUCLEUS project (PID2021-128757OB- I00) funded by MCIN/AEI/10.13039/501100011033 and NextGener- ationEU/PRTR. This work received support from the European Union’s Horizon 2020 research and innovation program through projects ACT- RIS.IMP (grant agreement No 871115) and ATMO_ACCESS (grant agreement No 101008004), by the Spanish Ministry of Science and Innovation through projects ELPIS (PID2020-120015RB-I00) and ACT- RIS-Espa˜na (CGL2017-90884REDT)). By the Junta de Andalucía Excel- lence, project ADPANE (P20-00136), AEROPRE (P-18-RT-3820) and by University of Granada Plan Propio through Visiting Scholars (PPVS2018-04), Singular Laboratory (LS2022-1) programs and Pre- Competitive Research Projects Pre-Greenmitigation3 (PP2022.PP34). Funding for open access charge, University of Granada/CBUA. Andrea Casans is funded by Spanish ministry of research and innovation under the predoctoral program FPI (PRE2019-090827) funded by MCIN/AEI/ 10.13039/501100011033, FSE “El FSE invierte en tu futuro”. Fernando Rejano is funded by Spanish ministry of universities through predoctoral grant FPU19/05340. Juan Andr´es Casquero-Vera is funded by FJC2021- 047873-I, MCIN/AEI/10.13039/501100011033 and NextGener- ationEU/PRTR. Elisabeth Andrews is funded in part by NOAA cooper- ative agreements NA17OAR4320101. Thanks to the NOAA Global Monitoring Laboratory for the use of the CCN counterWind-dispersed pollen grains emitted from vegetation are directly injected into the atmosphere being an important source of natural aerosols globally. These coarse particles of pollen can rupture into smaller particles, known as subpollen particles (SPPs), that may act as cloud condensation nuclei (CCN) and affect the climate. In this study, we characterize and investigate the ability of SPPs of 10 Mediterranean-climate pollen types to activate as CCN. A continuous flow CCN counter (CCNC) was used to measure the activation of size-selected (80, 100 and 200 nm dry mobility diameter) particles at different supersaturations (SS). Hygroscopicity parameter (κ) for each SPP type and size has been calculated using κ-K¨ohler theory. Organic chemical speciation and protein content has been determined to further characterize pollen solutions. Furthermore, the surface activity of SPPs has also been investigated by using pendant drop tensiometry. All studied SPP samples show critical supersat- uration (SSCrit) values that are atmospherically relevant SS conditions. Hygroscopicity κ values are in the range characteristic of organic compounds (0.1–0.3). We found that organic speciation and protein content vary substantially among pollen types, with saccharides and fatty acids being the only organic compounds found in all pollen types. A clear relationship between SPP activation and its organic composition was not observed. This study also reveals that all SPPs investigated reduce the surface tension of water at high concentrations but at diluted concentrations (such as those of activation in the CCNC), the water surface tension value is a good approximation in K¨ohler theory. Overall, this analysis points out that pollen particles might be an important source of CCN in the atmosphere and should be considered in aerosol-cloud interactions processes.BioCloud project (RTI2018.101154.A. I00) funded by MCIN/AEI/10.13039/501100011033, FEDER “Una manera de hacer Europa” and NUCLEUS project (PID2021-128757OB- I00) funded by MCIN/AEI/10.13039/501100011033NextGenerationEU/PRTREuropean Union’s Horizon 2020 research and innovation program through projects ACT- RIS.IMP (grant agreement No 871115)European Union’s Horizon 2020 research and innovation program through project ATMO_ACCESS (grant agreement No 101008004Spanish Ministry of Science and Innovation through projects ELPIS (PID2020-120015RB-I00) and ACT-RIS-España (CGL2017-90884REDT)Junta de Andalucía Excel- lence, project ADPANE (P20-00136), AEROPRE (P-18-RT-3820)University of Granada Plan Propio through Visiting Scholars (PPVS2018-04), Singular Laboratory (LS2022-1) programs and Pre- Competitive Research Projects Pre-Greenmitigation3 (PP2022.PP34)Funding for open access charge, University of Granada/CBUASpanish ministry of research and innovation under the predoctoral program FPI (PRE2019-090827) funded by MCIN/AEI/ 10.13039/501100011033Spanish ministry of universities through predoctoral grant FPU19/05340FJC2021- 047873-I, MCIN/AEI/10.13039/501100011033 and NextGenerationEU/PRTRNOAA cooperative agreements NA17OAR432010