Overview and Seasonality of PM10 and PM2.5 in Guayaquil, Ecuador

The focus of this study is the assessment of total suspended particles (TSP) and particulate matter (PM) with various aerodynamic diameters in ambient air in Guayaquil, a city in Ecuador that features a tropical climate. The urban annual mean concentrations of TSP (Total Suspended Particles), and particle matter (PM) with various aerodynamic diameters such as: PM10, PM2.5 and PM1 are 31 ± 14 µg m−3, 21 ± 9 µg m−3, 7 ± 2 µg m−3 and 1 ± 1 µg m−3, respectively. Air mass studies reveal that the city receives a clean Southern Ocean breeze. Backward trajectory analysis show differences between wet and dry seasons. During the dry season, most winds come from the south and southwest, while air masses from the peri urban may contribute as pollutant sources during the wet season. Although mean values of PM10 and PM2.5 were below dangerous levels, our year-round continuous monitoring study reveals that maximum values often surpassed those permissible limits allowed by the Ecuadorian norms. A cluster analysis shows four main paths in which west and southwest clusters account for more than 93% of the pollution. Total vertical column of NO2 shows the pollution footprint is strongest during the dry season, as opposed to the wet season. A microscopic morphological characterization of ambient particles within the city during the wet and the dry season reveals coarse mode particles with irregular and rounded shapes. Particle analysis reveals that samples are composed of urban dust, anthropogenic and organic debris during the dry season while mainly urban dust during the wet season.


ESPOL
Escuela Superior Politécnica del Litoral GIOVANNI Geospatial interactive online visualization and analyze infrastructure HR SEM High resolution scanning electronic microscopy image HYSPLIT Particulate matter with aerodynamic diameter < 1 µm PM 10 Particulate matter with aerodynamic diameter < 10 µm PM 2. 5 Particulate matter with aerodynamic diameter < 2.5 µm TOPAS Turnkey optical particle analysis system

Introduction
Aerosols originate both naturally or anthropogenically, and their composition, their effect on atmospheric processes and their physical and chemical properties are key components to understanding interactions with biogeochemical cycles that could affect rain and drought (Andreae 1995). Aerosol can be classified by sizing as: (a) fine mode particles corresponding to those with sizes < 1 µm, and (b) coarse mode particles with aerodynamic diameters > 2.5 µm and < 10 µm (PM 2.5-10 ) that vary from natural origin e.g., fungus, bacteria, pollen, or sea salt. Under natural conditions, fine particles, corresponding to the cloud condensation nuclei and ice nuclei, are suspended in the clouds and they may precipitate to the Earth's surface through raindrops in a process known as wet deposition (Seinfeld and Pandis 1998;Pöschl 2005). However, if no precipitation is present the concentration of airborne particle increases Fröhlich-Nowoisky et al. 2012). These airborne particles ultimately reach the Earth's surface through convective transport, diffusion, etc. (Pöschl 2005). Trade winds assist air masses to transport airborne particles longer distances, known as long range transport, which can occur through natural causes such as volcanic eruptions (Kahn and Limbacher 2012) and dust transport (Abdelkader et al. 2017;Di Biagio et al. 2017). The main drivers affecting air quality and aerosol composition in cities stems from incomplete combustion during biomass burning, a common occurrence in many industries. For instance, biomass burning is usually considered one of the main culprits of deteriorated air quality (Ault et al. 2012;Zhang et al. 2014). While biomass burning is typically carried out in rural locations, the air quality of peri-urban areas can be altered by anthropogenic emissions as well. This may be due to the proximity between the countryside and big cities. This is more evident over those places with no buffering zones and that passing through an expansion of agriculture areas that promote biomass burning. Further potential factor affecting air quality is the high levels of smog (García-Franco 2020) and recurrent dust events from nearby deserts (Andreae et al. 2015;Soto-Coloballes 2020).
To understand South American air pollution in cities, Table 1 summarizes some relevant works conducted in selected cities in the region. These studies analyzed the concentration of both PM 10 and PM 2.5 fractions. The concentrations found were related to conditions of human population, altitude, land use, and geographical position. For instance, early studies report that Cuenca, the third largest city in Ecuador, measures high daily particle concentrations of PM 2.5 (particulate matter < 2.5 µm) levels up to 30 µg m −3 (Astudillo-Alemán et al. 2015;Parra 2018). While in nearby cities such as Lima, Peru measure PM 2.5 concentrations of 26 µg m −3 which is accounted by the considerable number of industries in those regions (Silva et al. 2017). Hence, Table 1 Summary of air quality among different cities in South America NA not available a Vargas et al. (2012) b Astudillo-Alemán et al. (2015), Parra (2018) c Silva et al. (2017) d Zalakeviciute et al. (2018bZalakeviciute et al. ( , 2020  such levels pose a threat to those with sensitive respiratory systems (children and the elderly) which is corroborated by those findings carried out in a global context (Lelieveld et al. 2015). Aerosol concentrations of PM 10 are at high levels in cities like Rio de Janeiro and Sao Paulo, Brazil, which have a greater human population and major exchange markets (Martins et al. 2017), this is similar to Bogota, Colombia (Vargas et al. 2012). Recent studies on the effect of pollution in Andean populated cities, including Quito, revealed that environmental issues could be serious threats to human health (Bravo Alvarez et al. 2012;Cazorla 2016;Samek et al. 2018;Zalakeviciute et al. 2018aZalakeviciute et al. , 2020. One of the main forms of transportation in these Andean cities is the use of internal combustion (IC) engines. The atmospheric pressure in these cities is much lower (~ 15%) than the pressure at sea level (Zalakeviciute et al. 2018b). The IC engines operating in Andean cities are less efficient (as they were designed to operate at one atm), hence emitting greater amounts of hydrocarbons/PM than those IC engines operating at lower elevations (Bandowe et al. 2018).
Seasonality is a sub-analyze that is not common in Latin America urban pollution studies, as referred to in Table 1. There are many examples in northern hemisphere studies (Juneng et al. 2009;Hand et al. 2012); particularly due to the distinctive aerosol particle concentration during winter, spring, summer and autumn. Hence, this study will be focused on the intercomparison between wet and dry seasons which are common over the Tropical region.
The aim of the present study is (i) to provide a detailed characterization of ambient aerosol concentrations from total suspended particles (TSP), PM 10 , PM 2.5 and PM 1 based on 14 months of continuous monitoring in the city of Guayaquil, Ecuador; (ii) to relate the spatial and temporal patterns to meteorological conditions of the city; and (iii) to conduct fine/coarse particle morphology and determine the chemical composition of PM collected at various sampling locations. The present work is mainly focused on particulate matter and not on other pollution traces.

Study Area and Instrumentation
Guayaquil is a coastal city located in the southwest sector of Ecuador at 4 m asl with approximately 3 million inhabitants. It is one of the most populated cities in Ecuador and its major industry is international trade due to its harbor located at the end of the gulf of Guayaquil with an estimated area of 344.5 km 2 (Delgado 2013). The exponential growth of the city's population and trade market have made it a congested urban city with significant environmental threats (Rissler et al. 2014;Zalakeviciute et al. 2020). However, there are surprisingly few objective studies on the environmental effects and pollution stemming from population growth and rise in industry. During the early 1990's, an external consulting company, Espey & Houston Inc-COPADE, prepared the first report regarding the air quality of the city. The report defined hotspots from north to south and reported an estimated of 50 µg m −3 with heavy concentrations over the industrial site "Isla Trinitaria" due to electric power plants at the time (Borbor and Barriga 2007). In addition, several local studies from 2000 and 2016 reported values of PM 10 at approximately 50 µg m −3 (Eficacitas 2007) and of PM 2.5 at ~ 25 µg m −3 (FLACSO et al. 2008).
The geographical location of the studied area is located at: 2° 10′ 2.1′′ S, 79° 53′ 29.93′ W at 10 m asl; in which ambient aerosols were measured with a dust particle counter that was located outside a storage cellar. The data was accessed online from the AirQweb 2018 website (https:// www. airqw eb. com/ Main, last visited on June 20, 2018). Initially, the data were gathered using an external particle counter TOPAS from Turnkey Instruments (Cheshire UK from LATAN) installed from September 1, 2015 to November 17, 2016 UTC-5, with 10-min time resolution. The instrument allows the measurement of TSP and mass concentrations up to 6000 µg m −3 , PM 10 , PM 2.5 and PM 1 ; the relevant online data is summarized and presented in Table 2. Due to the fact that the TOPAS instrument went to maintenance and it was put offline on November 17 at 17:50 (UTC-5), the TOPAS was replaced with a new particle counter, Dylos DC1100 with PC interface (Gramsch et al. 2020). The performance of the low-cost aerosol sensor (Dylos DC1100) compared to the TOPAS instrument was adequately validated prior to the data collection. The intercomparison between TOPAS and Dylos DC1100 is presented in a graphic (included as supplementary S1 material) and it shows an acceptable correlation. Other researchers have reported the evaluation of the lowcost monitoring (Dylos) with more expensive instruments and correlation with similar results (Steinle et al. 2015;Vercellino et al. 2018).
Guayaquil weather correspond to the tropical savanna climate (Aw) corresponding to the Köppen climate classification; in which, mean temperature fluctuated from 23° to 26 °C and annual precipitation rate from 500 to 1000 mm and it varies from year to year. The seasonality split on two seasons: the wet season runs from the end of December to May and the dry season runs from July until November while December and June are considered as transitional periods. Annual seasonality has been characterized by many authors (Cañadas Cruz 1983;García-Garizábal 2017) to be due to trade winds related to the cold Humboldt current that converge near the equatorial. These periods may differ each year, but the trend patterns seem to remain consistent (Thiel et al. 2007).

Air Masses Trajectory Analysis
Air mass characterization was conducted using Hybrid Single Particle Lagrangian Integrated Trajectory (HYSPLIT) (Stein et al. 2015;Draxler et al. 2018) in which the air counter location was taken as the starting point and with historical data from January 1, 2015 to December 31, 2017. Using a combination of air mass studies encompassing a 96 h backward trajectory at a 200 m height above ground level (https:// ready. arl. noaa. gov/ hypub-bin/ trajt ype. pl? runty pe= archi ve) and a georeferenced list of the areas being studied (e.g. mining, industry, energy and water infrastructure), we were able to identify pollution hot-spots and emission sources, which can be compared to determine how the region is affected by different industries. Using this data, we have been able to identify how industrial pollution is translocated to Guayaquil.
Anthropogenic activities such as: burning fuel, vehicle emissions, power plants and others are among the different sources of NO 2 (Seinfeld and Pandis 1998). It is believed that higher NO 2 concentrations play a crucial role of anthropogenic sources. A remote atmospheric sounding instrument (TROPOspheric Monitoring Instrument, TRO-POMI) that provide spatial continuous information of the earth surface at a daily scale, was used to survey air pollution of the city. We use the OMI Total Vertical Column NO 2 product with level-2 NRTI:2018-10-24 (http:// www. tropo mi. eu/ data-produ cts/ nitro gen-dioxi de) last visited on 22 July 2020. The selected periods (Aug-Nov 2018; Aug-Nov 2019) were averaged daily map tiles for the dry season and selected periods (Jan-May 2019) were averaged daily map tiles for the wet season. The methodology used in our work closely follows the published works by other researchers (Krotkov et al. 2017;Zhao et al. 2020).
Additionally, some products from NASA Giovanni (Geospatial Interactive Online Visualization and Analyze Infrastructure) (Acker et al. 2014;Luna et al. 2018) web server were used to visualize aerosol presence, in which there were: (i) Aerosol Optical Depth at 550 nm (AOD 550nm ) at a wavelength of 550 nm from the moderate resolution imaging spectroradiometers (MODIS) on the satellites Terra and Aqua combined dark target deep blue AOD products (MOD08_D3_V6_ and (MYD08_D3_v6), and (ii) Tropical Rainfall Measurement Mission TRMM (TRM_3B42_Daily_v7 product), (https:// giova nni. gsfc. nasa. gov/ giova nni/), last visited on May 15 2018. Wet season (Jan-May, in blue color), dry season (Aug-Nov, in orange color) and transition season in white, as presented in Fig. 1

Sampling and Microscopy Analysis
Carbon membranes (disk of 12-mm diameter, brand Agar Scientific) were placed at three selected locations: (a) outside the city, (b) residential area and (c) city center. The sampling locations were uniquely selected as to represent the various types of airborne particle matters. Thus, the sampling encompassed a heavy-duty truck highway or the "highway to via costa" outside of the city, a residential area that is highly populated, and the final selected sampling location is city center or the "Archivo Historico". The wet season in Guayaquil is composed of heavy rains (precipitation or ~ 200 mm per month), flooding, high humidity (RH > 80%) and other related issues. Therefore, due to climate conditions the collection of samples during the wet season is not as practical as during the dry season. Additionally, the interpretation of samples collected during the wet season can be difficult and a good assessment may require an analysis (physical and chemically) of rainwater samples for comparison (Holecek et al. 2007). Due to the complexity of the wet season, only one location was considered. This location is named "Colon". The location was selected because of being in the vicinity of the TOPAS. City center is heavily trafficked by people, taxis, small passenger vehicles, etc.; and it is close to the airport cellar measurement station. During the sampling, carbon membranes were left for approximately 4-7 days in open environments of the selected locations, refer to Fig. 4c and Table 3. The carbon membrane disks were placed at 2 m height above the ground, facing SW winds and with a cover protection over the border. After the samples were collected, they were kept isolated in a cold environment at 7 °C in a refrigeration unit. Collected aerosol samples were exposed to gold sputtering for 120 s to render them stable under the electron microscope beam. The FEI-SEM with an integrated energy spectrometry (EDS) was used to study the chemical composition of selected samples using Genesis spectrum software (Baron and Willeke 2001;EDAX 2006). The methodology for chemical characterization of ambient aerosols by Moffet (2011), Burkhardt (2010) mentioned that SEM techniques are good for certain elements but not constituent (e.g. carbon, nitrogen), fine and coarse morphology. The SEM images were analyzed at a working distance of 10 mm, the detectors used were: Everhart-Thornley detector ETD and Back Scattering Electron Detector BSED with magnifications of: 3000x, 12,000 × until 100,000x; while performing EDS analysis, all analysis were taken by point to point to the sample on the surface of a selected structure. Elements like Au and Pd are considered traces of the gold sputtering. The sampling methodology on the collection of particles was conducted following closely the work in the literature (Morillas et al. 2016).
The data analysis was performed using the software Dylos Logger version 3.1.0, and R program version R 3.5.1 coupled with rstudio and packages openair, maps, maptools (Carslaw and Ropkins 2012;Carslaw and Beevers 2013;Grange et al. 2016) and Igor pro version 6.3.7 map script templates (Pöhlker et al. 2019) and ArcView for maps and data analysis and plotting (Zhao et al. 2020).

Meteorological Conditions and Particle Concentrations
Meteorological data showed that the wet season of 2016 received a cumulative precipitation of 1000 mm of rain (Fig. 1a). The majority of rainfall occurring in the first quarter of the year January to April, while the dry season ranged from August to November. Besides, certain months such as May, June, July and December are transitional months when sporadic rain episodes may occur. Average temperatures during the wet season and mean relative humidity (RH) of 27 °C and 78% occurred, while during the dry season, temperatures dropped to 24 °C with a mean RH of 73%, refer to Fig. 1b, c. As shown in Fig. 1a, the presented data agree with previous findings from Cañadas Cruz (1983); in which the tropical weather in the coast of Ecuador is described. Average wind speed during the daytime were 1 m s −1 during the wet season and increased to 2 m s −1 during the dry season, while nighttime average wind speeds were maintained at 1 m s −1 during the wet season and rose to 3 m s −1 (approximately two folds higher) during the dry season  (Fig. 1b, c). In regards of the wind direction, it showed a continuous path from south and south-west wind throughout the year. There are some occasional periods during the daytime, when winds from the north when wind speed remained quiet low (< 1 m s −1 ). During the rainy season, most of the large particles, particularly those of TSP and PM 10 were affected due to wet deposition, when rain helps TSP and PM precipitate to the ground rather than remain floating in the air (Pöschl 2005;Zalakeviciute et al. 2018a;Zhou et al. 2020). This effect is mainly noticeable with averaged values about 30 µg m −3 for TSP and 20 µg m −3 for PM 10 , respectively. Nonetheless, the concentrations of TSP and PM 10 sporadically spikes during the absence of rain Fig. 1a, f, and g and Table 2; in which values peak to 160 µg m −3 for TSP and 98 µg m −3 for PM 10 , respectively. Regarding PM 2.5 and PM 1 were at about 6 µg m −3 for PM 2.5 and 1 µg m −3 for PM 1 . Nonetheless, sporadic spikes in PM 2.5 and PM 1 can be found the whole year, no matter the seasonality, see Fig. 1a, h and i and Table 2; in which values peaked to 23 µg m −3 for PM 2.5 and 8 µg m −3 for PM 1 .
In contrast, during the dry season for TSP and PM 10 , Fig. 1a, f and g and Table 2, shows averaged values about 29 µg m −3 for TSP and 20 µg m −3 for PM 10 , respectively. Meanwhile, sporadic spikes in TSP and PM 10 can be found during the absence of rain during the dry season with peaked values of 140 µg m −3 for TSP and 93 µg m −3 for PM 10 . The average concentrations of PM 2.5 and PM 1 , Fig. 1a, h and i and Table 2, are 7 µg m −3 for PM 2.5 and 1 µg m −3 for PM 1 . Meanwhile, there are sporadic spikes in PM 2.5 and PM 1 and they take place mostly during the dry season with peaked values. Moreover, there are higher wind speeds and variations of wind direction, refer to Fig. 1b, c, during the dry season rather than the wet season, which may contribute to the transport of PM during the dry season.
According to the Ministry of Environment in Ecuador MAE air quality norm Book IV Annexes 3 and 3a, annual average PM 10 should not exceed 50 µg m −3 , while the reported data are 21 µg m −3 for the whole year (Table 2). PM 2.5 should not exceed 15 µg m −3 during a 24 h continuous measurement, while the reported data show 7 µg m −3 for the whole year. Furthermore, the maximum values reported show higher numbers particularly for the months of June and October of 2016 for PM 10 with values of 195 µg m −3 and 173 µg m −3 , respectively. Within the months of June and September of 2016 for PM 2.5 with values of 41 µg m −3 and 37 µg m −3 , respectively. Diurnal average data reveal that higher peaks tend to become more frequent in the afternoon and night times, refer to Fig. 2d, f. The windy season in Guayaquil starts on the second half of the year, while the windiest months in Guayaquil occurs from August to November (Johansson et al. 2018). According with our data, the average wind speed in February is 1 m s −1 while in September is 3 m s −1 ; throughout the year, the increase of the wind occurs in the afternoons around 14:00-18:00 and gradually increased until midnight. Seasonality plays a key role in here because high wind speed transport pollution during the dry season. Therefore, our study suggest that this may correspond to the double high peaks particularly during the dry months, refer to Figure Supplementary  S2. In regards of the wind directions, some studies suggested that higher PM 2.5 concentrations at night occur due the decrease of the planetary boundary layer; and consequently, enhance the presence of pollution in the air of megacities (Gago et al. 2013;Miao et al. 2019). According with Seinfeld and Pandis (1998) those events occur under lower temperatures at night and under steady wind speed conditions. In contrast, our study shows high wind speed events during the night (refer to Figure Supplement S2) no matter the seasonality. Thus, we imply the high PM concentration peaks are due to long-range transport events from the south gulf path. Moreover, it is important to remark that there is no mandatory reference in regards of TSP or PM 1 ; nonetheless, those values are also reported in Table 2 as reference.
Diurnal cycles of PM show that PM 10 and PM 2.5 average values are below those limits set by MAE air quality norm Book IV Annexes 3 and 3a. However, while considering the minimum and maximum data fluctuated substantially in all cases (TSP, PM 10 , PM 2.5 and PM 1 ) in terms of annual and seasonal occurrence, the collected data of aerosol particulates during the wet and dry season show a bimodal shape with a very prominent peak at 7:00-9:00, mainly for TSP and PM 10 during the wet season, refer to Fig. 2a, c. Thus, for PM 2.5 and PM 1 average concentrations, there are no distinguished peaks, rather small humps from mornings and afternoons, refer to Fig. 2e, g. In contrast, considering the minimum and maximum values during the dry season, the morning peak is shifted towards 10:00-11:00 h for TSP and PM 10 . A second peak occurs at 14:00-17:00 h for all particle size TSP, PM 10 , PM 2.5 and PM 1 during the dry season, refer to Fig. 2a-d. Hence, regardless of the season (wet/dry), the average values of PM 10 and PM 2.5 between 6:00 and 10:00 rise closer to the MAE limit; most likely due to the morning urban traffic no matter the seasonality, similar situations have been found by (Martins et al. 2017) studies at certain urban areas in Brazil and by in Quito (Zalakeviciute et al. 2018b); mainly during the wet season.

Backward Trajectories Analysis, Seasonality and NO 2
Backward trajectories of air masses analysis showed that during the wet season, there is an influx of air masses. Thus, air masses follow the trajectory from the gulf of Guayaquil, as well as the Guayas River, which steers air masses from south towards the north until reaching continental inland (Fig. 3a). During the dry season, air masses become south and south-west dominant (Fig. 3b). Besides, the precipitation rates during the wet season (see Fig. 3c) are prominent particularly in the coastland and shows a remarkable contrast to the dry season, as presented and correlated with the data presented in Fig. 1a.
There is an absence of strong winds during the wet season; in contrast, south and south-west winds are strong during the dry season. This can be observed in Fig. 1b, c, respectively.
Furthermore, the aerosol Optical Depth (AOD) shows the presence of aerosol in the atmosphere during the wet season and it shows a footprint of aerosol particle suspended in the air which fades out during the dry seasons. Despite the fact of high incidence of aerosol particles along the trajectory during the wet season, most PM are scavenged due to wet deposition before reaching the city (Fig. 3c, d) and Fig. 1a. Thus, higher aerosol peaks are seen during the absence of rain and increase in temperature, particularly during heavy traffic hours (07:00 in the mornings and 16:00 in the afternoons, UTC-5). Consequently, if there is more PM during the dry season it may be produced by two effects: (a) long range transport from south/west locations and (b) from the background urban pollution. The PM shown during the period 2015-2016 may correspond to Guayaquil's unique geographical conditions. During the wet season, the city experiences more regional pollution, while during the dry season, it may experience long range transport in addition to its local pollution. Similar findings agree with pollution and wind relationship as mentioned by Wang and Ogawa (2015) in Japan. Thus, we observe that any sign of external contamination may be occurring when strong winds increase (e.g. > 3 m s −1 ); possibly transporting external PM 2.5 and PM 1 from south and south west sources. The higher wind speeds and directions correlate well with the geographical location of the city, which is located on the west bank of a large river flowing into the gulf and finally into the Pacific Ocean. Ecuador's main port is located at the southern part of Guayaquil.
A cluster classification of all trajectories occurred from 2015 to 2016 showed four main paths along with their respective contribution: C1 (47%), C2 (28%), C3 (19%) and C4 (6%). The three main sources correspond to south (red dashed line, C1), southwest (blue dashed line, C2), west regions (green dashed line, C3) and east regions (purple dashed line, C4), refer to Fig. 4a, b. Thus, a temporal analysis shows the C1 and C2 are mainly occurring the whole year, while C3 and C4 are most common during the wet season. Main anthropogenic activities such as: mining, energy and water infrastructure, are located at the south and north of the city, while mining and energy, are both located at the west of the north of peri-urban area of the city, as presented in Figs. 4c/5a-c.
It seems that using the total vertical column of NO 2 as pollutant tracer, it sums up the local pollution during the dry season, refer to Fig. 5a, c, while during the wet season, pollution is mainly concentrated on the city center, refer to Fig. 5b. These findings correlate well with previously data presented in Fig. 2 (PM seasonality) and figure supplement S2. Total NO 2 column during the wet season, refer to Fig. 5b, shows an increase of the pollution spectrum from the north and northeast. The reason of this may be due to that most industries, energy and mining activities are concentrated nearby. They coincide with the cluster analysis C1, C2 and C3, which contribute with 94% of the air masses throughout the year. Consequently, part of those anthropogenic activities could be reflected as part of the pollution transport into the city center particularly at nighttime, according to Fig. 5c.
These trends are presented as guidance to understand critical pollution hotspot in the city. As previously presented by Zhao et al. (2020), total NO 2 column data are more precise over rural areas rather than urban and suburban areas in Canada in the city of Toronto. Moreover, ground measurements carried out by Cazorla (2016) in Ecuador in the city of Quito, showed low ozone and NO x levels during three months measurements from July to September, emphasizing that pollution ratio may imply further complex interactions at least for high altitude cities like in Quito at 2800 m asl.
Although the data presented provide good information of airborne particle concentrations within the city, this does not answer the question of the morphology and characteristics of the ambient particles found. To provide insight about the morphology and main characteristics of aerosol particles, some samples were taken during the dry season of 2018.

Morphology of Aerosol Particles During the Dry Season
The sampling period during the dry season of 2018 shows PM 2.5 concentrations of 9, 11 and 12 µg m −3 from outside the city, residential area and city center, respectively, refer to Table 3. Hence, these areas were selected for placing the carbon membranes for collecting the physical samples. Regarding the morphological characteristics the samples shown in Fig. 6a1-a3 are from the residential area; Fig. 6b1-b3 represents characteristics of samples collected from outside the city; and Fig. 6c, d shows samples collected from the city center. The low-resolution SEM images of lighter contrast (Fig. 6c1) represent processed raw SEM images to assist in visualization and comparison of the density of the "as-trapped" airborne PM present at the various studied locations. Higher resolution (HR) SEM and SEM-EDS analysis show that the structural morphology and chemical constituents of these urban dusts are quite complex and vary with location. SEM imaging shows that the dust samples are composed of very irregular shaped or anomalous structures, some have sharp-edges (dotted arrows in Fig. 6a2), as well as structures with much more well-defined shapes such as spherical or rod-like morphologies. Some of the dust particles appear to be composed of a collection or aggregates of smaller shapeless particles (Fig. 6b2, b3), while others appear to be composed of relatively larger solid angular and shapeless particles with a coarse surface texture having a high density of smaller particles attached to their surface (as highlighted with the solid arrows in Fig. 6a2, a3. The particle in Fig. 6a3) has a coarse-grained porous texture. The urban dust present in the city downtown area appears to be composed of particles with a variety of shapes and of a polydisperse size distribution (Fig. 6c1). The urban dust in Fig. 6c1 is composed of relatively large and sharp-edged (magnified view 1), spheroidal-shape particles (inmagnified view 2), straight and axial elongated rod-like structures (magnified view 3) and fiber-like with a radius of gyration (magnified view 4). The spherule in Fig. 6c2 is a "close-up" view of one of the smaller structures in Fig. 6c1. The spheroidal has a small fracture on its surface (solid arrow) revealing that the particle was formed of a collection of smaller shapeless dust particles (dotted arrow). However, the surface of the spheroidal is Fig. 6 a1-a3 Represents lower resolution SEM images of the "as-trapped" airborne PM after being strike on carbon tapes placed in open conditions from the residential area; b1-b3 from outside city; c1-c3 city center and d1-d6 city center during the wet season. The urban dust of Guayaquil city center appears to be composed an aerosol of a polydisperse size distribution and varying shapes (magnified view 1-4) rather smooth and does not contain deposition or absorbed materials attached to its surface as seen present in the samples collected in the residential area Fig. 6a2, a3. The HR-SEM shows submicron scaled particles form the larger structures Fig. 6c3. Figure 6d1-d5 represents SEM images collected on the surface of sample collectors positioned in the city center during the wet season. Figure 6d1 shows that large spherical shape structures are surrounded by other large irregular-in-shape structures. Higher resolution SEM imaging of the nearly spherical and irregular structures are shown in Fig. 6d2-d4. A higher resolution SEM analysis on a typical circular structure shows that its surface contains smaller irregular-in-shape deposits. Additionally, a closeup view on surface areas of the sphere where there is not accumulation of deposits show that the spherical particles have a very peculiar structure. It appears that the spherical structure is formed of smaller independent structures as shown by the presence of grainlike boundaries along the surface of the sphere (inset of Fig. 6d2). It is interesting to note that spherical structure that the spherical structure of the dry sample (Fig. 6c2) appears to have a similar morphology. SEM imaging analysis also show that some of the wet deposits are composed of large irregular and spherical like structures (Fig. 6d5). The inset of Fig. 6d5 shows the spherical structures of various diameters at the bottom of the cluster.
The localized pollution in Guayaquil is considered to be at a safe level by governmental mandated standards (MAE 2011). The reason for such safe levels of pollution is likely related to the unique geographical location of the city, which is surrounded by the gulf of Guayaquil and estuaries (Pinto and Slowey 2021), as well as some hills < 500 m in height with smooth slopes, refer to Fig. 4b, c. In contrast, wind coming from the Pacific Ocean, particularly during the dry season, may contribute to the long-range transport of pollutants to the most eastern parts of the city.
Single particle chemical analysis of the capture particles differ from their location; in which: (a) residential samples show high trace of N, Na, S, Ca, Si and Al; (b) city center samples show trace of N, Na, Mo, Ca during the dry season; (c) outside the city shows trace of Na, Mg, Al, Si, Ca and Fe, and d) N, Na, Fe, Mg, Al, during the wet season, refer to Figure Supplementary S3. Chemical elements like: C, Au and Pd correspond to the sampling methodology of the carbon pin and the gold sputtering. Na, Ca and S were the dominant elements found in most of the samples regardless the location and due to the size and their morphology we infer they can be some traces of anthropogenic emissions (e.g. urban dust, soil and industries). In that regards, similar findings agree which compare urban and rural areas in the Mediterranean area (Morillas et al. 2016;Genga et al. 2018;Perrone et al. 2018); in Poland by Samek et al. (2018); and in Quito, Ecuador by Zalakeviciute et al. (2020). Besides, sea salt aging plus the urban pollution of the city was found in Macao by Li et al. (2011) with traces of N, O, Na and S were found in certain samples at the residential areas and downtown below PM 2.5 ; our findings present similar pattern at certain locations but at particles corresponding to coarse mode PM 10 mostly.
Dust particles can be characterized by their unique surface morphology. It has been suggested that dust particles can have two exclusive and often coexisting types of surface roughness. One type is that the particle surface itself can be rough, while in the other type, the roughness is created by natural deposition or agglomeration of small particles on the surface of the large dust particle (Nousiainen 2009). In the context of particle morphology, dust particles are not spherical or nearly spherical, rather they have a very irregular shape and structure (Nousiainen 2009). The morphology (shape and size) of some of our structures resemble dust particles. Dust particles can have a variety of shapes, compactness, and be formed of mixtures of different mineral species (most notably birefringent) (Nousiainen 2009). Particles released from biomass burning or incomplete combustion have a nearly spherical shape and are composed mostly of carbon and oxygen (Chou et al. 2008); among other elements as: Ca-rich, Fe-rich, soot, and Si-rich (Fruhstorfer and Niessner 1994). In our collected SEM images, we can observe that irregular dust-like structures correspond to (PM < 2.5 µm) and (PM < 10 µm) structures that are nearly spherical. EDS on those spherical particles show them to be composed mostly of carbon and oxygen with traces of Ca (see Fig.  S3). It has been reported that structures with similar morphology, shape (nearly spherical) and composition (carbon, oxygen and traces of Ca) were generated from the biomass burning in the work of (Chou et al. 2008;Fruhstorfer and Niessner 1994). The deposition of compounds of many different mineral species on the surface of large dust particles is evident in our EDS analysis, refer to Fig. S3; in which, some particles related to mineral origins and some are mixed with anthropogenic origins (like biomass, burning or incomplete combustion) as mentioned by Pachauri (2013). It can be observed that the overall elemental composition of the dust particles is homogenous with only additions of chemicals on the surface. The EDS profiles enclosed in the orange (dry season) and blue (wet season) color box of Fig. S3.
In general, our microscopic analysis images showed a dominance of mineral and round shape alike fly ash components regardless of the seasonality. Two exceptions are Fig. 6b, d in which its rounded or oval shape and its chemical composition of C, O, N, Na and Ca are probably part of biological particles, similar shapes are described by Priyamvada et al. (2017). During the wet season, we noticed a predominance of larger particles (corresponding to TSP, except fig. S3-g) some shapeless and some with rounded shape and rough surfaces. From our early findings, we can 1 3 conclude that macro particles and C-rich, N-rich, Al-rich, Si-rich, Mg-rich, Ca-rich, Fe-rich and fly ash alike counted as dominant pollutants under the wet regime. This information is encompassed with the data presented in Fig. 5a-c (dry season) and 5b (wet season); in which, the pollution footprint is reduced during rainy periods.
We hypothesized that less variety and the absence of PM 10 and PM 2.5 particles may be due the current sampling technique that may be not suitable to work under adverse meteorological conditions, with a high precipitation rate (> 70 mm day −1 )and high relative humidity (> 90%); in such conditions, some ambient particles may be affected by hygroscopicity (Posfai and Buseck 2010;Pöhlker et al. 2014). Our preliminary findings on this matter are not conclusive but a close approach of the airborne particles found in the environment. Further studies must be conducted to fully understand the underlying mechanisms that the role of the city's landscape, meteorology, seasonality may have on urban aerosol composition.

Conclusions
After 14 months of continuous monitoring of particle concentrations in the city of Guayaquil, no environmental threat was found during the evaluated period. PM 10 and PM 2.5 mean values are within thresholds allowed by regulations from local government during 2015 and 2016.
Seasonality play an important role in the measurement of PM concentrations and their sources. We note that rains during the wet season, help to reduce the suspended particles within the city particularly during the months of February and March; while in January and May are more heterogeneous and PM increased at times during the absence of rain. Contrary, during the dry season, PM concentrations are more consistent from August to November regardless the particle size.
The south and south-west path helps to move fresh marine air masses into the city as long-range transport, which is corroborate with the meteorological data. Guayaquil city has a smooth topography, and it is protected by the gulf of Guayaquil and experiences a direct exposure of south and south-west trade winds.
Satellite data help to understand urban pollution with air masses backward trajectories, rain patterns, AOD and NO 2 . The total vertical column NO 2 helps: i) to visualize the local pollution footprint of the whole city, ii) to identify possible pollution hotspots and its relationship with anthropogenic activities, and iii) to show the pollution plume of the NO 2 differs through the year, broader during the dry season and narrower during the wet season. Our findings suggest that the complexity of urban pollution, may be part of the contributor enhanced by long-range transport pollutants and drought particularly during the second half of the year.
The microscopic analysis of the ambient aerosols shows that two main groups of samples: i) large and coarse mode particles (PM 10 and PM 2.5 ) with irregular shape and aggregate alike shapes, most likely corresponding to urban dust and mineral origins, and ii) other shapes like round shape and irregular shape particles with C-rich, N-rich, Al-rich, Si-rich, Mg-rich, Ca-rich and Fe-rich belonging to the coarse mode particles (PM 10 ). We note that although most of the sampled particles are uniform, some of them contain deposit or coating with materials on the surface that have additional chemical composition, which may be an indication of mixed or secondary organic aerosol particles. A comparison of SEM images of samples collected during the dry and wet seasons shows that there is a similarity in the morphological characteristics of the structures.
In summary, we considered the results presented in this study to be the first insight for further studies particularly over the interaction of a middle size south American city and its growth in the coming years. Hence, the microscopic analysis is the first step towards understanding the current state and nature of the suspended particles in the air.
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http:// creat iveco mmons. org/ licen ses/ by/4. 0/.