Spatial-temporal Variation and Local Source Identification of Air Pollutants in a Semi-urban Settlement in Nigeria Using Low-cost Sensors

Low-cost sensors were deployed at five locations in a growing, semi-urban settlement in southwest Nigeria between June 8 and July 31, 2018 to measure particulate matter (PM2.5 and PM10), gaseous pollutants (CO, NO, NO2, O3 and CO2), and meteorological variables (air temperature, relative humidity, wind speed and wind-direction). The spatial and temporal variations of measured pollutants were determined, and the probable sources of pollutants were inferred using conditional bivariate probability function (CBPF). Hourly PM2.5 and PM10 concentrations ranged from 20.7 ± 0.7 to 36.3 ± 1.6 μg m–3 and 47.5 ± 1.5 to 102.9 ± 5.6 μg m–3, respectively. Hourly gaseous pollutant concentrations ranged from 348 ± 132 to 542 ± 200 ppb CO, 21.5 ± 7.2 ppb NO2 and 57.5 ± 11.3 to 64.4 ± 14.0 ppb O3. Kruskal-Wallis ANOVA on ranks determined statistically significant spatial differences in the hourly-average pollutant concentrations. Diel variation analyses indicated that CO2, PM2.5, and PM10 peaked in the early hours of most days, O3 at noon while NO, NO2, and CO peaked in the evening. Most pollutants were of anthropogenic origins and exhibited the highest contributions from the southwest at most sampling locations. There were strong similarities between pollutants source contribution at two of the monitoring sites that were in residential areas with a frequently used paved road. Mitigation strategies need to be established to avoid further deterioration of ambient air quality that negatively affect public health.


INTRODUCTION
Environmental degradation is a major consequence of increasing population, economic development, rapid industrialization, and technological advancement (Morais et al., 2012;Mohsin et al., 2019). Air pollution is a well-known global problem in developed and developing nations due to its effects on public health and the environment at local, regional, and continental scales. Air pollution heightens the risk of respiratory and cardiovascular diseases, lowers the survival rates in the newborns, reduces the life expectancy of people with existing health issues, and produces acute effects on the exposed population (Dockery et al., 1993;Dockery and Pope, 1994;Pope, 2000). Clean air is a basic requirement for human comfort, health, and well-being and is being affected by air pollution, leading to human morbidity and mortality. Environmental issues such as acid rain, visibility impairment, alteration of the atmospheric radiation budget, and modification of cloud properties are also adverse effects of air pollution (Owoade et al., 2012;Tian et al., 2014). Air pollutants have been studied from both human and environmental health perspectives. Regulatory and advisory bodies such as the United States Environmental Protection Agency (U.S. EPA) and the World Health Organization (WHO) have established standards/guidelines for daily and annual concentrations of the different air pollutants. However, these values are exceeded in many locations.
Air pollutants may arise from natural or anthropogenic sources, and they can be transported over long distances, crossing both local and international boundaries (Owoade et al., 2012;Fawole et al., 2016). The major atmospheric pollutants include sulfur dioxide (SO2), oxides of nitrogen (NO/NO2/NOx), ozone (O3), oxides of carbon (CO/CO2) and particulate matter (PM). PM is particularly important because of its variable chemical composition, size, morphology, frequent occurrences of high concentrations, and its significant effects on human health and the environment (Sumesh et al., 2017). Ozone is formed in the troposphere through photochemical reaction of NOx and volatile organic compounds (VOCs) (Seinfeld and Pandis, 2016). Significant sources of NO in urban settings are motor vehicles and power generation systems (Ul-Haq et al., 2015). NOx is primarily emitted as NO with a smaller proportion of NO2. However, the reaction with O3 rapidly converts NO to NO2. The CO concentrations found in the lower atmosphere are typically associated with incomplete combustion activities as well as the oxidation of atmospheric methane.
In addition to the variable strength of the emission sources, variations in air pollutant concentrations are associated with the meteorological conditions, local environmental properties (land area, size, vegetation cover and population), and the availability of other pollutant precursors (Khiem et al., 2010;Wang et al., 2018). Notable spatial and temporal variations have been found for CO within the boundary layer up to scales of high inter hemispheric gradients (Yashiro et al., 2009). PM concentrations at any location vary with the local source strengths and other processes such as atmospheric chemical reactions, deposition, precipitation, and new particle formation (Sumesh et al., 2017). Hu et al. (2018) and Zhang and Jiang (2018) reported that local emission sources and potential effects of meteorology at meso -and synoptic scales might influence PM concentrations. The mixing layer depth and local weather dynamics drive the dispersion of source emissions at local scales while long-range transport affects the movement of PM over regional, transboundary, and continental scales. Multiple authors (Owoade et al., 2012;Guttikunda and Gurjar, 2012;Jayamurugan et al., 2013;Cheng et al., 2013;Gogikar and Tyagi, 2016;Gogikar et al., 2018) reported that meteorological variables such as wind speed and direction, air temperature, relative humidity, and rainfall govern the advection and dispersion of pollutants in any region.
A major challenge in monitoring air pollutant concentrations across an area and the influence of meteorology is the availability of measurement instruments. Most federal reference method (FRM or research grade) instruments are costly and difficult to move because of their size, power requirements, and ancillary systems such as sampling manifolds. Thus, they are typically stationed in fixed locations (Karagulian et al., 2019). Data obtained from the FRMs yield limited information on spatial variations of pollutants. Recent improvements in micro-technology have led to the development of portable monitoring systems for particulate matter and gaseous pollutants. These monitors are low-cost, very mobile, and can be deployed indoors and outdoors (Kumar et al., 2015;Sousan et al., 2016;Manikonda et al., 2016;Zikova et al., 2018;Crilley et al., 2018;Masiol et al., 2018;Li et al., 2020).
Existing studies of particulate air pollution in this area of Ile-Ife, Nigeria focused on the elemental characterization and receptor modeling to identify the likely particulate pollutant sources Owoade et al., 2015). The only major point source in the area is a scrap steel and iron smelting facility (Owoade et al., 2015;Ogundele et al., 2017). Owoade et al. (2016) identified four sources in both the PM2.5 and PM2.5-10 size fractions that were labeled: soil (44%), savannah burning (26%), scrap processing (18%), and vehicular emissions (12%) for PM2.5 and soil plus biomass burning (71%), sea salt (22%), scrap processing (5%) and vehicle emissions (including tire wear) (2%) for PM2.5-10. The biomass burning in the PM2.5 fraction would likely also include cooking emissions that include both the biomass burning for heat and the emitted cooking particles. All of these measurements were made either adjacent to the iron/steel facility or on the Obafemi Awolowo University campus. Thus, limited information was available on the PM concentration variability across the area and its influencing factors.
If the temporal and spatial variability and directional characteristics of the pollutants can be determined, they can be employed to identify hot spots, pollutant sources and provide the information needed to implement appropriate mitigation measures. Therefore, this study monitored particulate matter and gaseous pollutants concentrations in Ile Ife, Nigeria, analyzed their temporal variations and spatial variability, and provided directional information of their measured concentration using low-cost sensors. Omokungbe et al. (2020) made an initial report of the PM2.5 and PM10 results in which they examined the relationships between the measured PM values with temperature and relative humidity values. Higher PM concentrations were observed at lower temperatures and there was no variation of PM with RH when the RH < 80%. There were effects of RH for values > 80% likely due to the hygroscopic growth of the particles above the likely deliquesce relative humidity value. However, detailed spatial and temporal analyses were not made and those analyses and the resulting information about gaseous and particulate pollutions is the focus of this work.

Study Area
The study was conducted in a fast-growing city, Ile-Ife in southwestern Nigeria (7.33°N and 4.31°E, 294 m A.S.L.). It has an average population of about 502,000 (NPC, 2006) that was projected to have grown to 886,300. A major road passes through the city with several minor roads (paved and unpaved) linking different parts of the city for the easy movement of vehicles and passengers. The city hosts Obafemi Awolowo University (OAU), Obafemi Awolowo University Teaching Hospital Complex (OAUTHC), and a private polytechnic college. In addition to these institutions, a steel smelting factory is situated on the outskirts of the city along the Ife-Ibadan Expressway. Climatically, it experiences two major seasons (wet and dry). The dry season occurs from November to March the following year while the rainy season runs from April to October. During the dry season, there is little precipitation with high solar radiation, clear sky conditions, and high air temperatures. Cold, dust-laden northeasterly trade winds from the Bodele Depression in the Chad Basin in the Sahara Desert can blow across the region delivering significant quantities of dust for several days. In the wet season, the near surface flow is controlled by southwesterly winds emanating from the Atlantic Ocean and resulting into moderate temperatures and an abundance of rainfall (Falaiye et al., 2013).

Air Monitoring
Gaseous pollutants (CO, CO2, NO, NO2, and O3) and size segregated particulate matter fractions (PM2.5 and PM2.5-10) were monitored using low-cost monitoring sensors in a monitoring package at five locations in Ile Ife, southwestern Nigeria. The sampling locations were distributed across the city and the choices of the location were considered appropriate based on open space and the absence of any obstruction that would hinder the free circulation of air. Three sites (P1, P2 and P3) were located within the Obafemi Awolowo University campus, and two locations (P4 and P5) were within a fast-growing residential area of the city with one location generally downwind of an industrial area situated very close to the busy Ife-Ibadan Expressway. Fig. 1 shows the sampling site locations in the city. The five low-cost sensor packages were simultaneously deployed at all the chosen locations. Each low-cost monitor was placed 3 m above the ground level. Brief descriptions of each site are presented in Table 1. Air monitoring took place from June till July 2018. The duration was short as this is a preliminary study on air quality research in southwestern Nigeria.
The sensor packages used in this study were AlphaSense low-cost portable air monitoring devices. Optical particle counters (OPC-N2) measured particulate matter while gaseous pollutants were monitored with electrochemical sensors. The OPC sensor records particles per milliseconds in 16 size bins and estimates PM concentrations in the range of 0.38 to 17 microns with a maximum particle count of 10000 s -1 (Crilley et al., 2018). The OPC-N2 converts the particle concentration to mass concentration corresponding to the PM metrics: PM2.5 (for particulate matter ≤2.5 µm) and PM10 (for particulate matter ≤10 µm). The electrochemical sensors operate on amperometric  Eleweran Residential area further away from iron smelter and major road, closer to a minor road.
sensor principles that sense selected toxic gases at the parts per-million-level in an environment. Each unit also has temperature and relative humidity (RH) sensors as well as a 2-D sonic anemometer to measure wind speed and directions. Further details of the fundamental principles and operational procedures of the sensor and associated atmospheric variables sensors were given by Popoola (2012)  The major advantages of this sensor include low cost, low power use, high sensitivity, fast response, selectivity, multiple species monitoring, near real time measurement, and simple operating procedures (Popoola, 2012;Popoola et al., 2018). At each site, a box containing the sensor nodes were mounted at about 3 m above the ground on a 6 m pole and powered with a rechargeable battery. The low-cost sensors were operated simultaneously at the five sites for two months (June 8-July 30, 2018 in the wet season). The sensors recorded the gaseous and PM concentrations, and the meteorological variables every second and stores 20 s averages of these parameters. For the continuous measurements, the data were retrieved from a universal serial board (USB) stick every two weeks and then underwent further processing with suitable computer packages and software developed by the Department of Chemistry, University of Cambridge, UK. The PM mass concentrations and meteorological data were averaged to hourly values and further statistical analyses were performed using the Openair R software package (Carslaw, 2015;Carslaw and Ropkins, 2012) and SigmaPlot V14.

Statistical Analysis
The summary statistics were calculated. Normality was tested under the hypothesis: Ho: there is no normal distribution pattern with data obtained from the study area to determine to use either parametric or non-parametric methods. The confidence intervals for this study were set at 95% such that p-values ≤ 0.05 resulted in the acceptance of the null hypothesis. The normality test was the Kolmogorov-Smirnov test. Spatial variations in the measured variables were assessed using the Kruskal-Wallis analysis of variance (ANOVA) on ranks (Kruskal and Wallis, 1952) with confidence level set at 95% under the hypothesis that there was spatial variation (median values were not equal) in pollutant concentrations and meteorological parameters observed across the monitoring sites. Dunn's method (Dunn, 1964) was used to compare the site to site (pairwise comparison) variation at a 95% confidence level.

Source Directionality Using Conditional Bivariate Probability Function
The directional dependence of local sources of PM affecting each sampling sites was explored with conditional probability function (CPF) embedded in a bivariate polar plot (BPP) using the conditional bivariate probability function (CBPF). Bivariate polar plots provide information on the potential source influences by wind speed and direction on the receptor site. The CBPF estimates the probability that high pollutant concentrations (above a predetermined threshold value) are associated with a given wind direction and speed (Kim et al., 2003;Uria-Tellaetxe and Carslaw, 2014;Sumesh et al., 2017). It is based on the ratio of the number of samples in a wind sector and wind speed interval with concentrations greater than the given threshold to the total number of samples in that sector and wind speed interval. The CBPF values were calculated using Eq. (1): where m∆θ,∆u is the number of samples with wind direction from a given ∆θ and wind speed interval in the interval of ∆u having a concentration c greater than a threshold x and n∆θ,∆u is the total number of samples of that combination of wind speed and direction intervals (∆u, ∆θ). In this study, x was set at the 75 th percentile.

Coefficient of Divergence
Coefficient of divergence (CoD) is a parameter used for the evaluation of degree of divergence of two datasets. This analysis was performed to understand the variability or degree of uniformity of pollutants species measured simultaneously in two sites during a specified period. The CoD was evaluated as follows: where j and k represents the paired sampling sites; xij and xik stands for the average concentration for the same pollutants at the two sites and p is the number of pollutants under study. CoD is known to be self-normalizing and can be obtained from short and long-term averages (Wongphatarakul et al., 1998). CoD values range from zero to unity where values close to zero reveals strong similarity between the two sites while values close to unity depicts between the sites (Wongphatarakul et al., 1998;Liu et al., 2017). Wilson et al. (2005) reported that CoD values of 0.2 were consistent with spatial homogeneity of air pollutants.

Coefficient of Variation
The Coefficient of Variation (CoV or r-squared) was also used to assess the spatial variation of pollutant species between sites. The r-squared values shows the variance of a pollutant in a site that can be explained by the pollutants at another site. Like CoD, r-squared values range from 0-1. The r-squared values greater or equal to 0.7 suggests strong similarity in the pollutant temporal variations between the compared sites while values closer to zero suggests significant difference in the pollutant variation between the sites.

Average Hourly Concentration of the Measured Pollutants
The results of the average (± standard deviation) of the concentration of the measured pollutants at the five sites are presented in Table 2. The maximum 1-hour average concentration for PM2.5 and PM10 are 36 ± 1.5 and 111 ± 14 µg m -3 (P4), respectively. The high PM concentrations can be attributed to the location of this site in the fast-developing residential area with unpaved roads and a smelting industry. PM2.5/PM10 ratios have been used to aid source identification (Gogikar et al., 2018). In this study, moderate values of PM2.5/PM10 ratio were found at P5, P2, and P4, (0.44, 0.43 and 0.43 respectively) showing high anthropogenic contributions to the PM concentrations at these sites likely including suspended dust from the roads. A previous study on the Obafemi Awolowo University campus using gravimetric methods reported daily mean PM2.5 and PM10 mass concentrations of 25.4 ± 2.4 and 37.2 ± 32.4 µg m -3 (Owoade et al., 2012) were lower than the mass concentration obtained in this study, but with a much higher PM2.5/PM10 ratio. The daily average PM10 mass concentrations in this study exceeded the 50 µg m -3 limit value recommended by the WHO. Also, the daily average values of PM10 were 2.0, 1.8 and 1.6 times the threshold limit values of 50 µg m -3 stipulated by EU (2014) at P1, P4, and P3, respectively.
The maximum one-hour average CO concentration was 542 ppb at P3 that was attributed to its proximity to a busy road. The vehicular emissions and other anthropogenic activities like cooking with solid fuels (Obaseki et al., 2017) contributed to the values observed at this site. NO maximum concentration of hourly average was 7.34 ppb at site P3, similar to the value of CO. Ul-Haq et al.
(2015) associated concentration of NO in cities to vehicular emission. The minimum average concentration of CO2 was found at P3 (428 ppm) and the maximum was 471 ppm at P4. The maximum hourly concentration of NO2 occurred at P4 (21.5 ppb), an area downwind of an industrial plume while the lowest concentration was found at the farm settlement (13.7 ppb). The highest ozone concentration was found at P4 (64.4 ppb). These results showed that maximum hourly-mean concentrations of PM2.5, PM10, CO2, O3 and NO2 were observed at P4, a fast-growing residential area and downwind of a smelting industry. The range of concentrations of O3, NO and NO2 observed in this study were similar to those reported by Hagenbjörk et al. (2017) at three different sites (regional, traffic and background) in Sweden.

Observed Meteorology during the Study Period
The descriptive statistics for the meteorological variables (temperature, relative humidity and wind speed) at the five sites are presented in Table S1. Mean (±SD) temperature across the sites ranged from 23.4 ± 3.7°C (P2) to 25.6 ± 3.2°C (P5). The Kruskal-Wallis test indicated a statistically significant differences (p < 0.05) in temperature across the sites except for P1/P3, P2/P3, and P4/P5. The mean (± SD) of hourly relative humidity ranged from 74.9 ± 16.0% (P2) to 79.9 ± 15.7% (P3). There was a statistically significant difference in the hourly average relative humidity observed across the monitoring sites (p < 0.05) except for P1/P5. The hourly averaged wind speed range across study sites was 0.05 to 2.48 m s -1 . The Kruskal-Wallis analysis of the hourly mean wind speeds indicated statistically significant differences (p < 0.05) except P1/P2 and P3/P5. Given that these sites are in built areas, it is likely that local surface obstructions are diverting the wind and potentially reducing the observed speeds. Low wind speeds during calm period limit the dispersion of emissions within a locality leading to increased pollutant concentration (Dai et al., 2020).
The prevailing wind directions shown in Fig. S1 reflect winds primarily from the southwest with some northwesterly and limited southeasterly flows. Calm winds were observed about 16.7% of the time. There were prevailing low wind speeds that aided pollutant accumulation. It can be seen that there are some notable differences among the sites with P4 and P5 having northwesterly winds than at the other sites. P1/P2 and P1/P5 are the only pairs that have similar median directions.

Spatial Variation of Pollutants
The distributions of hourly-averaged concentration of the pollutants at the five monitoring sites are presented in Fig. 2 as box and whisker plots. The result showed the median (black line within the box), 25 th and 75 th percentiles of the pollutants. The lower and upper whiskers showed the minimum and maximum pollutant's concentration during the study. The results of the detailed statistical analyses are presented in Tables S1 and S2 of the Supplemental Information file. Those results are summarized in Table S3.
CO had the largest interquartile range and thus, its highest variation at site P3. Generally, there were significant differences in CO across all sites. This pattern could be attributed to a variety of activities including open burning and vehicular emissions. The highest concentrations of CO being at P3 could be attributed to its location near a bus park. For CO2, there were significant differences among most of the sites. However, P1/ P5 and P2/ P3 did not have significant differences. P1 is at the farm at one end of the domain while P5 is at the other end of the domain and is closest to the Iron and steel mill and in the residential district. P2 and P3 are quite close to one another with no significant sources that could affect these two sites differently.
There were significant differences in the variation of NO across most of the study area. However, the Dunn's result showed that there were again no statistically significant differences between the two sites within the university (P2 and P3). Similar periods of cooking and associated burning activities were likely to occur within the residential areas. Significant differences in the O3 and NO2 concentrations were also found across most of the sites. Sites P1 and P2 did not differ for either the NO2 or O3 medians. O3 medians did differ between P1/P5, P2/P4, and P2/P5. P4 and P5 are likely to be more highly affected by the major road south of the residential area, local combustion, and possibly emissions from the iron and steel mill.
There were also significant differences in PM2.5 and PM10 values across most of the study sites. Greater variability in the particulate matter fractions were observed at P4, a residential area with many unpaved roads and extensive two-wheel motor vehicle transport. For both fractions, it was observed that there were no significant differences observed between P1 and P3. P2 and P4 also had similar PM10 median values although P4 had a much greater range of values.

Diel variations
The diel variations of CO concentration (Fig. 3) had two major peaks that occurred in the morning (06:00-08:00 hr) and at night (18:00-20:00 hr). The variation shown in the plots provides the 5 th and 95 th confidence intervals represented by the shaded areas. These times are Volume 21 | Issue 10 | 200598 typically periods of higher traffic volumes (P3 and P5) while domestic combustion and electricity generators are mostly use at night at P2 and P4. Thus, these sources likely drive this observed variation. Higher morning peaks were observed at P3 and P5 while the higher peaks at the other locations (P1, P2, and P4) occurred in the evening. The diel variations of CO2 were similar across the study area. Boundary layer dynamics and source emission patterns are the major factors that influence the variability of CO2 concentration (Wang et al., 2010;Fang et al., 2016). The hourly peaks observed for CO and CO2 at all locations peaked in the evening when local biomass burning for cooking were likely driving the observed patters. There was a secondary CO peak in the morning likely from morning cooking. The diel variations of NO, NO2, and O3 (Fig. 3) at the monitoring sites were very similar to one another. There was an increase in NO from early morning 03:00 hr until 14:00-15:00 hr. NO peaked prior to the peak of NO2 and O3 which showed that photochemical transformation of NO through titration by ozone until the NO was sufficiently low that ozone concentrations could rise (Wang et al., 2018).
The similitude in the temporal variations can probably be attributed to general combustion (transport and residential combustion emissions) in the study area. The concentration of NO2 decreased during the convective hours (12:00-15:00 hr) of the day when mixing and secondary pollutant formation was maximal but dilution and vertical transport of primary pollutants predominated. The particulate matter fractions (PM2.5 and PM10) showed similar diel variations (Fig. 3) at all stations with a major peak at 05:00-06:00 hr except at P5 that peaked earlier than the other stations. P4 and P5 are close to the iron and steel smelting factory where its major operations take place at night (6:00 pm-5:00 am). Unlike these sites, P4 and P5 had two peaks between midnight and 06:00 hr that could represent the time for batch processes of the factory's operation. Additionally, the PM2.5 and PM10 temporal variation in the study area can be assumed to be influenced by atmospheric dispersion playing significant role in transporting the pollutants from the source to the sites.

Day of the Week Variation of Pollutants
The day of week variations ( Fig. S2) showed that CO exhibited similar daily variations at all locations. There were higher concentrations during the weekdays than on the weekends (Saturday and Sunday) for most sites. P3 showed the largest drop in CO concentration on weekends. This reduction was expected because it was close to a major university road that was used less on weekends. CO2 followed a similar variation pattern at all the monitoring stations. Unlike CO, CO2 concentration increased on weekend days (Saturday and Sunday) compared to weekdays. This result may reflect people not working and thus present at their residences as well as additional cooking being done.
The daily concentration of NO2 ranged from 10 to 20 ppb during the study and showed a uniform temporal variation at the sites. There were no differences in the day of week variations of NO and O3 at the monitoring sites. The day of week variations of PM2.5 and PM10 at the monitoring sites were similar. This result could be due to similar sources, largely transportation and the transport of particles from iron and steel smelting factory at few kilometres from the sites. For most of the sites, the highest PM concentrations were observed midweek (Tuesday and Wednesday).

Source Identification and Characterization
3.6.1 Conditional bivariate probability plots Fig. 4 shows the CBPF results for CO. The residential area sites (P4 and P5) show opposite directions at relatively low wind speeds mostly under 2 m s -1 suggesting that local activities in this area are the primary sources. Given that the winds are generally from the southwest (Fig. S1), there is clearly substantial local emissions. These sources would include the minor road bordering the area and local combustion for cooking. P1 shows the highest probabilities in the southwesterly direction toward the campus area for wind speeds under 1.5 m s -1 . P2 and P3 on the OAU campus again show preferred southerly directions. P2 appears to be seeing emissions from the road and the area around the Fire Station whereas P3 is affected by the campus area to the south southwest to the west southwest.   5 presents the CBPF plots for CO2. They have very different patterns than those observed for CO with all of the higher probability areas appearing at higher wind speeds. There are much smaller relative variations in the CO2 concentrations since local emissions are superimposed on the ~400 ppm global background. P1 showed a strong influence to the southwest under higher wind speeds (~5 m s -1 ), a lower probability area to the southeast under lower wind speeds, some influence from the north at speeds between 2.0 and 2.5 m s -1 , and at very little low wind speed probabilities. Thus, there are none or very weak sources in the immediate vicinity of P1. There are large residential areas to the southwest to the southeast (city of Ife). Sources affecting P2 are on the OAU campus area to the west southwest. P3 shows a weaker influence from the campus area and a strong focused area to the east southeast. A large source is the Obafemi Awolowo University Teaching Hospitals Complex (OAUTHC) that is 4.8 km from P3 at 97.5°. They would be cooking using LPG as well as drawing considerable traffic in the area. Another smaller potential source is the Moremi High School is 2.5 km at 96.9° from P3. However, there is no cooking done there and food is brought in as needed.
Figs. 6 and 7 show the CBPF plots for PM2.5 and PM10. All these plots show that the highest probabilities for high concentrations occur when wind speeds are low (typically < 1 m s -1 ). The P3 plot suggests influence from the road that curves around the Fire Station building. However, directionality is not well defined at low wind speeds. Thus, the differences in patterns observed between P2 and P3 may not be as substantial as it appears. The P4 plot suggests emissions in the residential area even with only weak flows from the westerly and northerly directions. The P5 plot also points to the residential area, but again all under low wind speed conditions. Neither the P4 or P5 plots suggest an effect of the iron and steel mill. The only obvious difference between the two size fraction plots is at P3 where the PM10 plot shows areas of influence to the northwest and east southeast at higher wind speeds.
The CBPF plots for NO, NO2, and O3 are provided in Figs. S3-S5. There are notable differences between the NO and NO2 plots suggesting that NO is more local in origins while NO2 is   transported. At P4, there is a strong node to the SW at higher wind speeds that may reflect the high temperature combustion occurring at the iron and steel mill. The P2 plot shows a likely influence of combustion sources on the OAU campus such as cooking with LPG that did not affect P3. P3 again shows an NO area to the ESE (Fig. S4) similar to those seen in the O3 and PM10 plots at higher wind speeds suggesting a source that is likely to be the OAUTHC. Most of the high probability areas in the NO2 plots are at higher wind speeds around the outside of the plots similar to the O3 plots (Fig. S5). These patterns suggests that NO2 and O3 are likely in photolytic steady-state and then interact with the local NO emissions to reduce the ozone concentrations. However, the NO titration does not produce sufficient NO2 relative to the transported concentrations to produce high probabilities at low wind speeds (Wang et al., 2018).

Coefficients of divergence
The CoD values are presented in the Table 3. CoD values < 0.2 are considered similar in concentration values (Wilson et al., 2005). The most similar sites were P2 and P5. These sites are the residential areas near a heavily travelled paved road. Overall, there were similarities among

Coefficients of variance
The r-squared values are presented in Table 4. The variance of pollutants at most sites cannot be explained with by the same pollutant at the other sites except for CO2 at P3 and P1 which are locations within the university and presence of sources of combustion within close proximity. This variation could also be influenced by prevailing meteorology of the study sites.

Limitations
Due to limited availability of the monitoring equipment, this study was only conducted over a 2-month period during the later spring and early summer. Thus, it does not reflect seasonal variations in emission or meteorological conditions. Given differences in wind directions, precipitation, and human activities over an entire year, this work only provides a snapshot of this period.

CONCLUSIONS
This study used low-cost sensor-based monitors to measure gaseous pollutants, PM, and meteorological variables within a fast-growing semi-urban settlement over a period of 8 weeks. The study found temporal and spatial variations of the pollutants across the study area. Wind in the study area had speeds from 0.5 to 2.5 m s -1 mostly from the southwest suggesting sources in the identified wind direction. The observed spatial variations for each pollutant were statistically significant at a 95% confidence level. The diel variations of pollutants revealed that CO, CO2, PM2.5 and PM10 had two prominent peaks at morning and night times while NO, NO2 and O3 were observed to peak between 12:00-15:00 hr. Average daily concentration showed that the all the pollutants peaked majorly during the weekdays (Tuesday or Wednesday) except CO2 that generally peaked on Sunday. Local sources such as heating/cooking and traffic were likely major sources as well as the scrap processing factory to the SW of the study area. The area is a fast-growing semi-urban settlement with increasing urban development. Thus, there is a need to implement mitigation strategies to reduce local emissions of pollutants. Also, limiting emission sources to the southwest of the study area such as from the scrap processing plant could lower the pollution load across the study area.

SUPPLEMENTARY MATERIAL
Supplementary material for this article can be found in the online version at https://doi. org/10. 4209/aaqr.200598