Daily averages of AOD, EAE, SSA, ASYM, AAE, and the TPW and aerosol volume size distribution were calculated firstly based on the number of days with data, number of instantaneous data and then the monthly mean were obtained by calculating the average of the total number of days of a given month independently from the year. Finally, the monthly and seasonal means were analyzed to characterize the column-integrated aerosol properties Figure 1. The relationships of AAE vs. EAE were applied to infer the dominant aerosol types. AAE were calculated by using absorption AOD at 440, 675, and 870 nm versus wavelength at the logarithmic scale whereas EAE were also calculated from a linear regression of the same three mentioned wavelengths. The combination spectral difference between EAE (440 - 675) nm and EAE (675 - 870) nm were used as a graphic method for tracking mixtures of pollution containing dust based on AODs at a few wavelengths which was a good indicator to fine mode effective radius [5] . The seasonal aerosol size distributions binned by AOD was used to determine the variation of fine and coarse aerosols with increasing AOD. A straightforward graphical framework spaced on δα vs. EAE were employed to track mixtures of pollution aerosol with dust, to distinguish aerosol growth from cloud contamination and to examine aerosol humidification and finally a 120 hour back-trajectory calculation using the NOAA HYSPILT model was also used [6] [7] .

ICIPE Mbita station represents the most active aerosols monitoring station in Kenya with long-term data spanning from 24^th December 2006 to till date is located near the Southwestern border of Kenya and Uganda along the shores of Lake Victoria at the latitudes and longitudes of 0.42˚S, 34.20˚E respectively and 1125 m above sea level. Next to this station is also one of the third big city in Kenya, Kisumu City with a population of 409,928 [8] . Mbita is not only under the influence of lake-terrestrial air mass exchange [9] , but also are affected by dust aerosols caused by the transport highways and biomass burning from the nearby rural areas.

The data for the above mentioned station was obtained from the Air Recourses Laboratory HYSPLIT MODEL http://www.ready.noaa.gov/index.php and AERONET dataset (https://aeronet.gsfc.nasa.gov/new_web/units.html) and level 1.5 all points’ format data were used in this study. The level 1.5 data derived from the version 2 direct sun algorithm were used to analyze the AOD, EAE and TPW. The level 1.5 SSA, ASYM, and aerosol volume size distributions were derived from the version 2 inversion algorithm [10] .

Table 1 and Table 2 are the statistics of sun direct and inversion data for the AERONET ground station which shows the distribution of the data counts over different months and years. From Table 1, it shows some months with higher recoded days with data, for example in December 2016, the days with data were 27/31 whereas some days were with no recorded data. Because of daily changes of the local weather, the presence and the passage of clouds frequently interrupts clear-sky data leading to fewer days with recoded data. There is a direct relationship between the days with observed data and days with instantaneous

data as it is clearly shown in Table 1 and Table 2. As the number of observed days with data increases, there is also an increase days with instantaneous data though it’s not 1:1 ratio. During the hot dry season, most of the atmosphere is sky clear which leads to higher recorded of data as compared to the rainy season with overcast atmosphere. It was observed from Table 2 that total amount of the instantaneous data was 5287.

The daily variation of AOD440 nm was performed using a time series analysis. The annual mean of AOD for the ten years from 2007 to 2016 was about 0.25, 0.28, 0.18, 0.27, 0.27, 0.29, 0.29, 0.32, 0.26 and 0.56 respectively. It can be seen in Figure 3 that the mean value of AOD for the year 2014 was higher compared to other years because most of the observations were made during the dry season of July to August which is considered as the high aerosol loading period while 2009 was recorded the lowest AOD values of about 0.18 since observations were made for shorter period. The high variability in AOD440 nm indicates a variety of aerosol types over this study area due to the transported air mass of different origins and characteristics.

Figure 4 shows a plot of the monthly statistics of aerosol optical properties and total precipitable water content which are presented as box-and-whiskers diagrams. The variability of the median values in each month is labeled with a circle. It can be clearly shown that AOD showed a distinct seasonal variation for this rural area with high values that mainly occurred in June-July-August (dry season) while the lowest amount was occurred during March-April-May (wet season). The seasonal means of AOD 440 nm were approximately 0.25 ± 0.09 in DJF, 0.21 ± 0.06 in MAM, 0.38 ± 0.13 in JJA and 0.23 ± 0.08 in SON. The highest monthly mean values of AOD was in July (0.41 ± 0.07) while the lowest was in November 0.16 ± 0.03 and May 0.19 ± 0.07 Figure 4(a). The appearance of high values of AOD 440 nm in dry season was related with the contribution of mixed and marine aerosols as shown in Figure 7(b) though the previous literatures

work indicates that the dominant aerosols originates from dust and biomass burning particles emitted from surrounding rural areas [11] .

The monthly mean values of AOD was highest in July (0.41 ± 0.07) and lowest in November 0.16 ± 0.03 and May 0.19 ± 0.07. It can be seen that AOD, AEA and AAE graphs have similar patterns of variation such that when there is a rising of AOD parameter the other parameters also increase Figure 4(a), Figure 4(b) and Figure 4(e) whereas the relationship between AOD and TPW has negative correlation between them.

It was observed that low EAE values occurred during the wet season while high values were recorded during the dry season Figure 4(b). This may be related to high amount of TPW during wet season and low during dry season Figure 4(f). The annual mean values of SSA is 0.85 ± 0.10 while the relatively high monthly values of SSA occurred in the months of June to July (Dry season) Figure 4(c). This is due to the dominance of fine mode aerosols and low amount of water vapor in the atmosphere during this dry period. Ångström exponent (α) was a measure of the wavelength dependence of AOD and a good indicator of aerosol particle-size. It is clearly depicted from Figure 4(e) that the values of AAE are high in the dry months and low in the wet months. These results signify the presence of greater contribution of aerosols in the fine-mode to the extinction for the present study period [12] . Lastly the variations in ASYM Figure 4(d) are not apparent; the values only vary from 0.70 to 0.93.

Figure 5 shows the frequency percentage occurrence which is categorized by individual season for the AOD at 440 nm. Due to the different emission sources and atmospheric conditions, Figure 5 indicates that there is a large difference in the AOD distribution for each season although all the seasons indicate a unimodal distribution which can be seen from log normal blues curve. During the DJF period, the modal value was 0.18 of 60% frequency of occurrence while for the MAM season the modal values was 0.20 of 50% frequency of occurrence meanwhile during the dry season of JJA, the modal values was 0.35 with 40% frequency of occurrence and lastly for SON season, the modal value was at 0.25 with 40% also for frequency of occurrence.

It was also observed that there was inverse proportionality between the wavelength and the AOD, as the wavelength increased there was a decreasing in the values of AOD. During the dry season the values of AOD were higher than other seasons, the same trend was observed for the AE with higher values during JJA and lower in the wet season which was similar with what was observed by [13] .

Due to the breakup of large particles which have originated from multiple sources of particles, it is evident that these particles have grown growth mechanisms in the system which leads to bimodal graphs. From Figure 6, it is evident that the seasonal variation of aerosols volume size distribution at 440 nm in the rural area of Mbita are mainly of natural origin but with a moderate of influence of anthropogenic sources. The number of distribution is characterized by two modes at a radius of 0.25 to 0.3 µm for fine mode while the mass distribution for coarse mode is centered on 7 to 7.5 µm for all the four seasons. Figure 6(a) which is a transitional period between long rain season and short rain periods clearly shows dominance coarse particle in size range > 0.02 µm centered at a radius of 8 µm. During the long rain period Figure 6(b), the graph almost flattened at the radius of 0.2 µm remaining with only one peak centered at 8 µm. This may have resulted from growth mechanism of the fine particles from rain water droplets leading to coarse mode particles.

During the hot dry season Figure 6(c), both fine mode and coarse mode particles are clearly shown. The amount of water vapor in the atmosphere during this period is less. This may have contributed to the occurrence of both aerosols particles, though it can be seen that fine particles were more dominant with the highest peak of VSD at a value of 0.03 µm³/µm² than all other seasons. Figure 6(d), during the short rain season, the concentration of fine particles diminishes as it can be seen from the graph due to the interaction with rain water. The high values of VSD in the fine mode which is observed during the dry season is due to the contribution of more biomass burning during the land preparations and

cultivations for the crops and long range transport of dust and urban-industrial wastes [9] .

Different aerosols in the atmosphere can be classified into different types by using scatter plot techniques centered on their abundance such as AOD against α. AOD vs α which are both the functions of wavelength was preferred for this work because there was large instantaneous availability of data as compared to other cases like FMF, AAE and SSA which have averaged daily data as it was suggested by [14] . Five different classes of aerosols were identified by these techniques namely; Marine, Desert dust, Biomass burning/Urban Industrial and Clean Continental as seen in Figure 7. Different threshold values were used to classify the aerosol, for example marine influence aerosols ranges between 0 0.3, 0<α>0.85 which corresponds to 26.1% of total contribution Figure 7(b). Aerosols originates from desert dust and BB/UI ranges between 0.6 1.8, 0<α>0.6 and 0.6 1.8, 1<α>1.75 respectively with a small percentage of 1.6 while clean continental threshold values ranges between 0.6 0.3, 1<α>1.8 (26.1%). Lastly a large portion corresponding to 44.2%

of scattered points were grouped as mixed aerosols. It can be seen that aerosols originated from marine environment corresponding to (26.1%) of 3,554 because ICIPE Mbita station is situated next to Lake Victoria.

Figure 8 demonstrates the seasonal variations of aerosol classification at four different local seasonal periods over ICIPE Mbita namely; DJF, MAM, JJA and SON. From Figure 8(b), it can be seen that very few data points at low AOD 550 nm values during MAM season were outside the classification scheme which may be caused due to small uncertainties errors in dAE vs AE estimation while evaluating AOD 675 whereas other plots all over the data are within the classification scheme.

In JJA season Figure 8(c), AOD > 1.2 nm of (R_f ∼0.1 - 0.15 µm) is of η 90% however most aerosols clustering falls on the positive dAE with large clustering of fine mode AOD > 1 nm concentration found within the range of R_f (0.2 - 0.3) μm with η(30 - 50)%. The same scenario with large clustering of fine mode AOD > 1 nm concentration is found within the range of R_f (0.2 - 0.3) μm with η 30% occurring during MAM season Figure 8(b) however coarse mode occur more frequently in this season than in other seasons with (0.4 1.2) nm of η (30 - 50) % and R _f (∼0.05 - 0.10) µm were on the positive dAE. Both DJF and SON season show similar results although during SON there is large clustering of (0 0.6) nm and R _f (∼0.1 - 0.15) µm which are concentrated in dAE (0 - 0.2) positive side with η (50 - 70) % within the range of AE (0.8 - 1.6) nm, Figure 8(d).

The positive dAE values indicates that these aerosols are of bimodal distribution having a large coarse-mode fraction between η (50 - 70)%, whereas, the dAE values near to zero Figure 8(a) represents the lesser dominance of coarse-mode particles relative to fine-mode with η values of 70%.

This work emphasizes on obtaining the area of interest-scale characterization of the aerosol properties. It can be seen that small differences in the dAE vs α plots revealed establishing homogeneity in aerosol load although there was also present of small percentage of heterogeneity in aerosol load. The more positive dAE with increasing AOD is an indication of coarse-mode dominance under high AODs. As seen in Figure 8, even though several conclusions regarding these aerosol characteristics over ICIPE Mbita are apparent, the small scatter in the data leads to small visual uncertainties.

In order to understand and monitor the major sources of pollutants arriving at ICIPE Mbita, The Hybrid Single Particle Lagrangian Integrated Trajectory Model (HYSPLIT) [6] [7] of NOAA for backward trajectories analysis. The data was chosen based from Table 1 and Table 2 for the month of July in 2015. This month recorded the longest time duration data for 23 days with its corresponding

number of instantaneous data of 75. A 120 hours (5 days) back trajectories at different altitudes of 1,130, 2,000 and 3,000 meters above mean sea level was used since this station is at an altitude of 1,125 m. Generally, observations show Figure 9 shows most of the trajectories were originated from the southern hemisphere from south east direction. Figure 9(a), the trajectories at the altitude of 1,130 m amsl (red cycles) starts from the coastline of Mozambique in Indian Ocean and propagates along the coastline of Tanzania as it enters into Kenya, located at the border of Tanzania and the coast side of Kenya [15] [16] . This trajectory passes through industrial, agricultural and urban areas of Kenya toward western region of Kenya which might have carried large amounts of fine mode particles and industrial pollutants. For the trajectories, the altitudes of 2,000 and 3,000 originate from the central Ethiopian highland and Saudi Arabian respectively.

These trajectories propagate through the arid and semi-arid regions of Kenya, which brings in dust aerosols with them. Figures 9(a)-(d) shows air masses except level of 3,000 m from the 2^nd week until the 5^th week and for the

Figure 9(c) it indicates that the air masses had a strong maritime influence originating from the Indian Ocean and Madagascar Island, which are one of the major sources of causing sea salts aerosols and fine mode smoke plumes in East Africa [4] .

During summer in the northern hemisphere, the strong south easterly surface winds flow is predominantly as it can be seen from Figure 10. Through the application of a simple technique of wind analysis it can be seen that a persistent high-speed of winds of about (15 - 30) Knots in the form of low level jet streams, exists in the lower layers which originates from the tip of Madagascar which is then blown to Kenyan coastline as southerly. These southerly winds carry sea salts and dust from Indian Ocean and Madagascar Island respectively. Due to topographical intensification of the flow, the speed is generally reduced as winds

converges into western parts of Kenya leading to more concentration of sea salts and dust particles to the study area.