Joe Adabouk Amooliab,
Kwame Oppong Hackmanc,
Bernard Nanad and
Daniel M. Westervelt
*b
aUniversité Joseph Ki-Zerbo, Ouagadougou, Burkina Faso
bLamont-Doherty Earth Observatory, Columbia University, New York, NY, USA. E-mail: danielmw@ldeo.columbia.edu
cCompetence Centre, West African Science Service Centre on Climate Change and Adapted Land Use, Ouagadougou, Burkina Faso
dEcole Normale Superieure de l'Université de Koudougou, Koudougou, Burkina Faso
First published on 22nd July 2024
This study estimates PM2.5 concentrations in Ouagadougou using satellite-based aerosol optical depth (AOD) and meteorological parameters such as temperature, precipitation, relative humidity, wind speed, and wind direction. First, Simple Linear Regression (SLR), Multiple Linear Regression (MLR), Decision Tree (DT), Random Forest (RF), and eXtreme Gradient Boosting (XGBoost) models were developed using the available labeled data (AOD and meteorological parameters with corresponding PM2.5 values) in the city. The XGBoost model outperformed all other models that were used, with a coefficient of determination (R2) of 0.87 and a root-mean-square error (RMSE) of 15.8 μg m−3 after a five-fold cross-validation. The performance of the supervised XGBoost model was upgraded by incorporating a semi-supervised algorithm to use large amounts of unlabeled data in the city and allow for a more accurate and extensive estimation of PM2.5 for the period 2000–2022. This semi-supervised XGBoost model had an R2 of 0.97 and an RMSE of 8.3 μg m−3 after a five-fold cross-validation. The results indicate that the estimated 24 hour mean PM2.5 concentrations in the city are 2 to 4 times higher than the World Health Organization (WHO) 24 hour guidelines of 15 μg m−3 in the rainy season and 2 to 22 times higher than the WHO 24 hour guideline in the dry season. The results also reveal that the average annual estimated PM2.5 concentrations are 11 to 14 times higher than the WHO average annual guideline of 5 μg m−3. Finally, we find higher PM2.5 concentrations in the city's center and industrial areas than in the other areas. The results indicate a need for future air pollution policy and mitigation in Burkina Faso to achieve desired health benefits such as reduced respiratory and cardiovascular problems, which will, in turn, lead to decreased PM2.5 mortality rates.
Environmental significanceOuagadougou is a city with PM2.5 levels likely exceeding the WHO guidelines but with limited air monitoring stations to fully understand this. This paper addresses the limited data coverage by developing an accurate semi-supervised XGBoost model for estimating PM2.5 from satellite-based AOD and meteorological parameters using a combination of labeled and unlabeled data in Ouagadougou. The study also analyzed spatiotemporal variations of PM2.5 concentrations in Ouagadougou. The study reveals that the estimated PM2.5 concentrations in the city are higher than the WHO 24 hour guideline. It also highlights that the industrial and central areas are the major polluting areas of the city. This research underlines the potential of satellite-based AOD and meteorological parameters to increase data coverage of estimates of PM2.5 in cities with limited air monitoring stations while highlighting the need for a comprehensive modeling approach for better air quality decision-making. |
PM2.5 are particles with an aerodynamic diameter of 2.5 μm or less.12,13 The diameters of the larger particles in the PM2.5 size range are approximately 30 times smaller than that of human hair.14 The main sources of primary and secondary fine particles include electricity generation, industrialization, urbanization, and other processes that involve the burning of fuels like wood, heating oil, or coal, and natural sources like the dust storms from the Sahara Desert and forest fires.13,15–17
Ouagadougou is a city in Sub-Saharan Africa with PM2.5 levels likely exceeding the WHO guidelines due to dust from the Sahara Desert, re-suspension of road dust, emissions from traffic, widespread biomass burning, and frequent stable nocturnal atmospheric conditions but with limited air monitoring stations18–20 to fully understand this. This is due to the limited financial resources available to install these stations, which also contributes to the limited study of PM2.5 concentrations in the city. While it is true that many urban areas around the world struggle to meet WHO guidelines for PM2.5, Ouagadougou's specific environmental and socio-economic conditions exacerbate this issue, leading to consistently high pollution levels.21 Lindén et al.21 investigated the characteristics of the connections between Ouagadougou's metropolitan climate and air pollution. The author observed that Ouagadougou's air pollution situation was marked by significant spatial variations, high overall pollution levels, and extreme levels of coarse particles, frequently exceeding WHO air quality guidelines in all areas. Furthermore, Nana et al.20 quantified the concentrations of Nitrogen Dioxide (NO2), Sulfur Dioxide (SO2), BTEX (Benzene, Toluene, Ethylbenzene, and Xylene), and PM10 in Ouagadougou. Their findings showed that aside from downtown, where levels were frequently above the WHO Guideline, NO2 and SO2 concentrations in the city continued to be below the WHO-set guideline with ranges of 22 μg m−3 to 27 μg m−3 and 0.5 μg m−3 to 10.5 μg m−3 respectively. It was also found that the city had high quantities of BTEX (27.9 μg m−3) and PM10 (119 μg m−3 to 227 μg m−3) surpassing WHO guidelines.
Understanding the relationship between satellite data, particularly aerosol optical depth (AOD), meteorological parameters, and the concentration of surface PM2.5 will increase data coverage of estimates of PM2.5.22–28 Wang and Christopher22 investigated the relationship between hourly PM2.5 measured at the surface at seven locations in Jefferson county, Alabama for 2002 and column AOD derived from the Moderate Resolution Imaging Spectroradiometer (MODIS) on the Terra and Aqua satellites. Their findings showed that the satellite-derived AOD and PM2.5 had a strong correlation (r = 0.7), indicating that most of the aerosols were in the well-mixed lower boundary layer during the satellite overpass times. In addition, Van Donkelaar et al.23 estimated global ambient fine particulate matter concentrations from satellite-based AOD. They discovered that estimates of long-term average PM2.5 concentrations (1 January 2001 to 31 December 2006) at a spatial resolution of roughly 10 km × 10 km pointed to a population-weighted geometric mean PM2.5 concentration of 20 μg m−3 for the entire world. They found significant geographic agreement between the satellite-derived estimate and ground-based in situ measurements (r = 0.77; slope = 1.07) as well as between the satellite-derived estimate and noncoincidental observations elsewhere (r = 0.83; slope = 0.86). Tian and Chen29 developed a semi-empirical model using MODIS AOD, PM2.5, and meteorological parameters including specific humidity, air pressure, air temperature, and boundary layer height (BLH) to estimate, at a regional level, the hourly concentration of ground-level PM2.5 concurrent with satellite overpass. Their model was able to explain 65% of the variation in ground-level PM2.5 concentrations with a root-mean-square error (RMSE) of 6.1 μg m−3. Furthermore, Wang et al.30 and Xue et al.24 used a MODIS AOD, PM2.5, and meteorological data to develop a linear model for estimating surface PM2.5 in China. Their findings demonstrated good agreement between model-estimated PM2.5 and the observational data (RMSE of 23.0 μg m−3 and R2 of 0.72).
Given that surface PM2.5 data in Ouagadougou are sparse with the U.S. embassy site being the only measurement site, developing a model that makes use of the small amount of labeled data (AOD, temperature, precipitation, relative humidity, wind speed, and wind direction with corresponding PM2.5 values) and the large amount of unlabeled data (AOD, temperature, precipitation, relative humidity, wind speed, and wind direction without corresponding PM2.5 values) for estimating PM2.5 would be useful for extensive and accurate study of PM2.5 and help in air quality decision-making in the city.
This study aims to develop an accurate model for estimating fine particulate air pollution from satellite-based aerosol optical depth and meteorological parameters using a combination of labeled and unlabeled data in Ouagadougou. We also analyzed spatiotemporal variations in fine particulate air pollution concentrations in Ouagadougou. We hypothesized that an effective model could be developed for estimating fine particulate air pollution from satellite-based aerosol optical depth and meteorological parameters using a semi-supervised algorithm. We also hypothesized that Ouagadougou's fine particulate air pollution concentrations vary seasonally and are higher in central and industrial areas. According to the literature, Ouagadougou's air quality data challenges have not yet been analyzed using machine learning methods. To our knowledge, this is the first application of satellite AOD to surface PM2.5 in Ouagadougou, and one of the first on the African continent.
Satellite parameter | Spatial resolution | Temporal resolution | Source |
---|---|---|---|
AOD | 1 km | Daily | MODIS |
Precipitation | Resampled to 1 km | Daily | CHIRPS |
Temperature | Resampled to 1 km | Daily | Era5-Land |
Relative humidity | Resampled to 1 km | Daily | Era5-Land |
Wind speed | Resampled to 1 km | Daily | Era5-Land |
Wind direction | Resampled to 1 km | Daily | Era5-Land |
Daily AOD data were extracted from the Moderate Resolution Imaging Spectroradiometer (MODIS) Terra and Aqua Multi-Angle Implementation of Atmospheric Correction (MAIAC) Daily Level 2 Aerosol Product at spatial resolutions of 1 km × 1 km and band 550 nm (AOD550) for 16 areas in Ouagadougou for the period 2000–2022 using Google Earth Engine (GEE) and Python. The data availability of the MODIS satellite AOD in Table 1 is displayed in Table S2,† revealing significant gaps during the early 2000s (2000–2002) but more than 50% availability thereafter. With more than 50% of the data available since 2002, this does not pose any problems. The 16 areas in Ouagadougou were chosen based on the geographical distribution of the city, taking into account traffic, industrial, and commercial areas. The locations are on a regular grid due to the city's planned neighborhood layout. This method was chosen to cover the entire city. Because the ground-based weather observations from the ANAM station at Ouagadougou International Airport alone might not be a good representation of all the different areas in Ouagadougou, daily satellite weather observations including precipitation from the Climate Hazards Group InfraRed Precipitation with Station data (CHIRPS), temperature, relative humidity, wind speed, and wind direction from the ERA5-Land Daily Aggregated-ECMWF Climate Reanalysis were used as complement data for the same period as the observed weather parameters. These daily satellite weather observations were extracted for the 16 areas in the city where the AOD was extracted, which included areas that the ANAM station did not cover and the area covered by the station (Ouagadougou International Airport). The temperature was extracted directly from the temperature at 2 m band of Era5-Land. Relative humidity was extracted using the following equation.34
![]() | (1) |
The Wind speed was also extracted using the equation below;
![]() | (2) |
The wind direction was also computed using the following trigonometric expression;
![]() | (3) |
The ERA5 and CHIRPS data were 100% available for the collection period.
Satcorr,i = βi × Sati + β0,i | (4) |
The performance of the models was based on their coefficient of determination (R2) and the root mean square error (RMSE) as shown in Fig. S1.† Given the same climatic zone, these models were then applied to correct all satellite data for the other fifteen (15) locations in Ouagadougou, where ground observation data were not available.
The hourly PM2.5 data were averaged according to the local time for the period it was collected and used to obtain the hourly profile of PM2.5 concentrations at the U.S. embassy in Ouaga 2000. The 24 hour PM2.5 concentrations were averaged to obtain daily PM2.5 concentrations and used for the model development.
PM2.5 = β0 + βAOD × AOD | (5) |
Song et al.36 modified the simple linear eqn (5) by introducing meteorological parameters to obtain a multivariable equation for estimating PM2.5 as follows:
PM2.5 = (α + ε1) + (β1 + ε2) × AOD + (β2 + ε3) × TEMP + (β3 + ε4) × RH + (β4 + ε5) × WS | (6) |
PM2.5 = (α + ε1) + (β1 + ε2) × AOD + (β2 + ε3) × TEMP + (β3 + ε4) × RH + (β4 + ε5) × WS + (β5 + ε6) × Precip + (β6 + ε7) × WD | (7) |
Before each model development, Pearson correlation was performed to determine the Pearson correlation coefficient between PM2.5 and the parameters and to determine their strength and direction. In each model, 80% of the data was used for model development and 20% for model testing. The models were evaluated using the R2, RMSE, and significant F values. The F-value in regression determines if your model of linear regression fits the data set more closely than a model without any predictor variables. A regression table that includes the F-statistic and associated p-value is what you will get as an output when you fit a regression model to a dataset. If the regression model fits the data better, the p-value would be smaller than the selected significance threshold.
The model was trained using the inductive training principle. In the inductive training principle, the semi-supervised model is trained to keep the rules observed during the training process so it can generalize well to new unseen data.37 Eighty (80) % of the labeled data was combined with the unlabeled data, and 20% of the labeled data was used for model testing. The model was evaluated using the R2 and RMSE metrics obtained from a five-fold cross-validation. The semi-supervised XGBoost model was applied to estimate PM2.5 in the remaining fifteen (15) areas of the city where PM2.5 data were not available because these areas are in the same climatic zone with the location for which the model was trained. Estimation was performed for the period in which the model was trained (2000–2022) to study the growth and distribution of PM2.5 for this period.
Parameters | r | RMSE | MAE |
---|---|---|---|
Observed precipitation and CHIRPS precipitation (resampled to 1 km resolution) | 0.87 | 7.09 mm day−1 | 2.67 mm day−1 |
Observed relative humidity and Era5-Land relative humidity (resampled to 1 km resolution) | 0.96 | 5.65% | 4.53% |
Observed temperature and Era5-Land surface temperature | 0.92 | 1.24 °C | 0.92 °C |
Observed wind speed and Era5-Land wind speed (resampled to 1 km resolution) | 0.93 | 0.33 m s−1 | 0.25 m s−1 |
Observed wind direction and Era5-Land wind direction (resampled to 1 km resolution) | 0.89 | 22.36° | 14.04° |
The observed precipitation and CHIRPS satellite precipitation followed the same trends and were strongly correlated with a Pearson correlation coefficient of 0.87, an RMSE of 7.09 mm day−1, and an MAE of 2.67 mm day−1. These findings are similar to those of Plessis and Kibli,43 who reported a Pearson correlation coefficient of 0.77 between observed precipitation and CHIRPS precipitation over South Africa. The observed weather parameters and Era5-Land reanalysis parameters showed similar trends and were strongly correlated, with Pearson correlation coefficients ranging from 0.89 to 0.96, RMSE ranging from 0.33 to 22.36, and MAE ranging from 0.25 to 14.04. These Pearson correlation coefficients are similar to those found in Assamnew and Mengistu44 and Gleixner et al.45 who reported Pearson correlation coefficients ranging from 0.90 to 0.96 between Era5 weather parameters and their corresponding observed weather parameters over East Africa. The strong correlation between the observed weather parameters and the satellite weather parameters is due to the fact that Ouagadougou is mostly cloud-free except during the rainy season. The correlations in the rainy season are lower, ranging from 0.79 to 0.89, compared to the dry season where correlations range from 0.83 to 0.97, as shown in Table S1.† The reduced correlations in the rainy season are likely due to the clouds, which hinder satellite accuracy in capturing surface conditions.
Two peaks with different magnitudes were observed, one in the morning and one in the evening. The evening peak was higher with PM2.5 concentrations at about 105 μg m−3 while the morning peak was lower with PM2.5 concentrations at about 70 μg m−3. The morning peak was observed between 5:00 am and 9:00 am, and the evening peak was between 4:00 pm and 11:00 pm. The morning peak is likely due to morning vehicular traffic when people are going to work, and the evening peak is likely due to evening vehicular traffic when people return home after work. Cooking activities in the morning and evening also contribute to these peaks. Nana et al.20 and Ouarma et al.19 observed similar peaks in administrative sites in Ouagadougou. These peaks are also similar to those observed by McFarlane et al.46 in Kinshasa-Brazzaville and Bahino et al.47 in Accra and Abidjan. Govender et al.48 also found that PM2.5 concentrations increased during early mornings and late afternoons in the Vaal Triangle Area of South Africa. The high concentrations of PM2.5 in Ouagadougou in the evening is very likely due to the dynamics of the boundary layer (BL), which produces high dilution rates during the day and low dilution rates after sunset.49 Fig. S2† illustrates the temporal changes between weekends and weekdays in Ouagadougou. PM2.5 concentrations increase from Friday to Monday, followed by a decrease from Monday to Friday. Saturdays, Sundays, and Mondays consistently exhibit the highest pollution levels in the city, with PM2.5 concentrations exceeding 65 μg m−3. This trend could be attributed to increased weekend traffic associated with leisure activities, shopping, and social events, resulting in higher vehicle emissions, though further activity data are needed. Elevated PM2.5 levels on Mondays may stem from accumulated emissions over the weekend and heightened commuter traffic as the workweek begins. Notably, Saturdays in Ouagadougou are not considered full weekends, with work activities extending until 2 pm, which contributes to lower pollution levels compared to Sundays and Mondays, as depicted in Fig. S2.†
Fig. S3† presents an analysis of the potential source regions and contributions of PM2.5 using Trajstat models and the Hybrid Single-Particle Lagrangian Integrated Trajectory (HYSPLIT) model, developed by the National Oceanic and Atmospheric Administration's (NOAA) Air Resources Laboratory (ARL).50,51 We examined HYSPLIT back trajectories with meteorological data from the Global Data Assimilation System (GDAS) for extremely polluted days (hourly PM2.5 concentrations greater than 350 μg m−3 and daily mean PM2.5 concentrations greater than 200 μg m−3) in Ouagadougou. These days include January 28, February 7, March 13, 14, 19, and December 17. Back trajectories were calculated for 24 hour periods at three receptor heights (500 m, 1000 m, and 1500 m). The analysis reveals that on most of the extremely polluted days, the wind fields at these different heights predominantly originated from the east and northeast, extending to regions as far as Niger and northern Nigeria. An exception was noted on January 28, when back trajectories arriving at 500 m originated from the south. The West African region is influenced by dust from the northeast, where the Sahara Desert is located,19,47 suggesting that the dust likely originated from the Sahara Desert but was slightly redirected by atmospheric dynamics.52 This is also consistent with the conditional bivariate analysis in Fig. S4,† which shows that wind from the east and north has some association with high PM2.5 concentrations. These insights underscore the influence of meteorological conditions on PM2.5 levels in Ouagadougou.
Parameters | r |
---|---|
AOD-PM2.5 | 0.72 |
Relative humidity-PM2.5 | −0.55 |
Temperature-PM2.5 | 0.11 |
Precipitation-PM2.5 | −0.26 |
Wind speed-PM2.5 | −0.04 |
Wind direction-PM2.5 | −0.33 |
PM2.5 shows a weak negative correlation (r = −0.26) with CHIRPS precipitation. The negative correlation means that as precipitation increases, PM2.5 decreases and vice versa. Wet deposition (not in the dry season) serves as the primary sink for atmospheric particulate matter, hence increases in precipitation result in a decline in particle concentrations.53,54 The weak correlation means that most of the variations in PM2.5 are not directly explained by precipitation. All the Era5-Land variables showed a negative correlation with surface PM2.5 except temperature, which showed a positive correlation. Relative humidity had the strongest negative correlation (r = −0.55) followed by wind direction (r = −0.33) and wind speed (r = −0.04). The negative correlations imply that an increase in any of these variables will result in a decrease in PM2.5.55 High Relative humidity leads to hygroscopic growth of PM2.5, which can increase mass and increase deposition rates similar to wet deposition by precipitation hence decreasing the concentrations of PM2.5.19,56 However, it is worth noting that high relative humidity can also lead to more hydroxyl radical (OH) driven by sunlight in Ouagadougou and more aqueous formation of aerosols.54 Low wind speed leads to a stagnant atmosphere favoring PM2.5 accumulation.57
In contrast, high wind speeds promote PM2.5 dissipation. High wind speeds can generate and transport dust. Similarly, the direction of the wind can determine the concentration of PM2.5. If the wind direction is away from an area with PM2.5 sources, the PM2.5 particles are transported away from that area, resulting in lower PM2.5 concentrations in that area and hence a negative correlation.57 From the conditional bivariate analysis of the relationship between wind speed, wind direction, and PM2.5 concentrations displayed in Fig. S4,† slow-moving winds (less than 1.5 m s−1) from the north and medium-moving winds (1.5–3.0 m s−1) from the east and north are associated with high PM2.5 concentrations. Conversely, high-moving winds greater than 4 m s−1 are associated with lower PM2.5 concentrations. However, the weak correlations of the parameters with PM2.5 imply that most of the changes in PM2.5 are not directly explained by the parameters and their relationship might be more complex. The correlation between temperature and surface PM2.5 was 0.11. An increase in temperature leads to some increase in PM2.5 concentrations.53,54 However, the weak correlation observed in our findings indicates that most of the variations in PM2.5 are not directly explained by temperature.
![]() | ||
Fig. 3 Performance of the statistical regression models and the non-linear machine learning models at the U.S. Embassy site (Ouaga 2000). |
The SLR model developed is:
PM2.5 = 75.61 × AOD + 38.36 | (8) |
The model has an R2 of 0.52 after a five-fold cross-validation, indicating that MODIS AOD explains about half of the variations in surface PM2.5. The RMSE of the model after cross-validation was 38.3 μg m−3 and the significance F was 1.23 × 10−22 (less than α = 0.05). These findings are similar to those of by Van Donkelaar et al.58 They obtained R2 of 0.59 between MODIS AOD and PM2.5 over Eastern China. Wang and Christopher22 obtained an R2 of 0.49 at seven locations in Jefferson County, Alabama. Koelemeijer et al.59 also obtained an R2 of 0.36 between AOD and PM2.5 at some locations in Europe.
The MLR model developed is:
PM2.5 = 174.11 + 69.27 × AOD − 1.65 × T − 1.00 × RH − 14.04 × WS − 1.54 × Precip − 0.08 × WD | (9) |
The MLR model has an R2 of 0.67 and an RMSE of 33.7 μg m−3 after a five-fold cross-validation with a significance F of 8.85 × 10−30 (less than α = 0.05). One explanation for this is that the inclusion of meteorological parameters leads to better prediction of surface PM2.5, the MLR model explains 0.67 of the variations in surface PM2.5 with a smaller RMSE compared to the SLR model. These findings are consistent with the findings of Tian and Chen.29 Their regression model explained 0.65 of the variations in observed PM2.5.
The decision tree model explains 0.70 of the variations in surface PM2.5 with an RMSE of 34.0 μg m−3 after a five-fold cross-validation. The decision tree model performed the least in the non-linear models. The decision tree model uses a single tree and offers less hyperparameter tuning and hence was not able to better capture complex relationships. The random forest model explains 0.85 of the variations in observed PM2.5 with an RMSE of 16.6 μg m−3 after a five-fold cross-validation. McFarlane et al.60 used a random forest model for correcting low-cost sensors PM2.5 data in Kampala, Uganda, and had a similar performance (R2 of 0.86). The XGBoost model explains 0.87 of the variations of surface PM2.5 with a lower RMSE of 15.8 μg m−3 after a five-fold cross-validation outperforming all the models. XGBoost offers more hyperparameters tuning and this results in its outstanding performance. These findings clearly show that non-linear models perform better than linear models in estimating surface PM2.5.61,62
In all the models, MODIS AOD is the most important feature in PM2.5 estimation, followed by relative humidity and temperature as shown in Fig. S6.† Precipitation is a less important parameter in the models' estimation. This means that the impact of precipitation in explaining the variability of PM2.5 is less significant than that of AOD, relative humidity, temperature, wind speed, and wind direction.
Emission source and meteorological attribution data in Ouagadougou is scarce. A study by Lindén18 found that extremely stable nocturnal atmospheric conditions occurred on 80% of the days at the beginning of the dry season in Ouagadougou. The frequent occurrence of these stable conditions and the intra-urban thermal wind fields result in limited ventilation and dispersion of urban pollutants.18 Increased air pollution concentrations during stable conditions with low wind speeds have been identified as a factor in air quality issues in semi-arid African cities.63 We confirm this is also the case in our dataset in the conditional bivariate analysis of the relationship between wind speed, wind direction, and PM2.5 concentrations in Fig. S4† where low to moderate wind speeds tend to be associated with high PM2.5. Lindén18 also highlighted that key pollution sources include road dust re-suspension, transported dust, traffic emissions, and biomass burning. In developing countries, many households rely on biomass burning and using open fires, often in poorly ventilated areas, leading to severe air pollution.64 While biomass is more commonly used in rural areas, it is also prevalent in low-income urban areas.63 A study on household energy in Ouagadougou, Burkina Faso, found that over 80% of urban residents used biomass as their primary energy source.65 Additionally, in developing countries, traffic often consists of many old, poorly maintained vehicles and numerous two-stroke engines, which are highly polluting.20,21
The model has an R2 of 0.97, and an RMSE of 8.3 μg m−3 after a five-fold cross-validation, indicating that the model explains 0.97 of the variations in PM2.5 with the lowest RMSE. The model is validated using an independent dataset for the period of August 2023 to October 2023 from the tapered element oscillating microbalance (TEOM 1400a, a federal equivalent method gravimetric PM2.5 monitor)66 located at Université Joseph Ki-Zerbo in Ouagadougou. The model effectively explains 91% of the measurements recorded by the TEOM (R2 = 0.91) with minimal error (RMSE = 3.70 μg m−3) as shown in Fig. S8.† These results significantly enhance the reliability and confirm the generalizability of the model. These findings are similar to those of Bougoudis et al.67 because their semi-supervised ANN model explained approximately 0.9 of the variations in air pollutants in Athens, Greece. Similarly, the semi-supervised k-nearest neighbours (KNN) model proposed by Zhao et al.68 in China had an R2 of 0.97. Fig. 4 also shows the estimation of PM2.5 by the model on the whole data (labeled and unlabeled) after testing. It is observed that in the whole city, days in the dry season are associated with high PM2.5 concentrations (2 to 22 times higher than the WHO 24 hour guideline of 15 μg m−3), and days in the rainy season are associated with low PM2.5 concentrations (2 to 4 times higher than the WHO 24 hour guideline of 15 μg m−3), and the same trend is observed every year. This trend is similar in other Sahelian cities that share the same climatic conditions and experience similar Harmattan winds as Ouagadougou. Garrison et al.69 found that the mean 24 hour PM2.5 concentrations in Bamako are 43 μg m−3 (2.9 times higher than the WHO 24 hour guideline of 15 μg m−3) during the rainy season and up to 500 μg m−3 (33.3 times higher than the WHO guideline of 15 μg m−3) during the dry season. In Dakar, the mean 24 hour PM2.5 concentrations are 26 μg m−3 (1.7 times higher than the WHO guideline of 15 μg m−3) during the rainy season and up to 330 μg m−3 (22 times higher than the WHO guideline of 15 μg m−3) during the dry season.70 This implies that aside from internal sources of PM2.5 pollution within these cities, the Harmattan winds are a significant external contributor to the high levels of PM2.5 pollution. This natural phenomenon, combined with local pollution sources, results in the extremely high PM2.5 concentrations observed during the dry season.19
The yearly trend of PM2.5 in the city is not direct. The PM2.5 levels fluctuate irregularly, showing variations that are not easily predictable or following a simple upward or downward trend throughout the years in the study period. A direct trend means that PM2.5 levels show a predictable or linear increase, decrease, or stability over the years in the study period without significant fluctuations or irregularities. The yearly PM2.5 trend increases and decreases but, in general, there is a slight increasing trend which is not significant (p-value of 0.65 > 0.05). The average annual estimated PM2.5 concentrations range from 58.2 μg m−3 to 72.1 μg m−3, which is 11 to 14 times higher than the WHO guidelines of 5 μg m−3 annually and 6 to 8 times higher than the U.S. Environmental Protection Agency (EPA) guidelines of 9 μg m−3 annually. Significant increases in the trend were observed in 2004, 2010, and 2015 whilst significant decreases were observed in 2002, 2003, and 2020. Interannual variability and weather can affect year-to-year changes in PM2.5. Generally, the city is growing in population and economic activity,33 which can lead to enhanced PM2.5 concentrations. In 2020, however, when there was a lockdown due to COVID-19, it was observed that the PM2.5 trend decreased significantly. This is because of the reduction in industrial activities and heavy traffic due to the lockdown. However, the overall trend largely depends on the intensity of dust from the Sahara Desert and the variability of weather conditions21 as well as the intensity of emissions from unpaved roads, heavy traffic, and industrial activities in a given year.19,20,63 March has the highest PM2.5 concentrations with an average of 107 μg m−3 because the dust from the Sahara Desert transported by Harmattan winds reaches its peak in March. Dust from unpaved roads and biomass burning are also major contributors to higher concentrations in the dry season.20 August has the lowest PM2.5 concentrations with an average of 16 μg m−3, attributed to the significant precipitation during the month. These findings align with those of Awokola et al.,71 who observed elevated PM2.5 levels in March and April and lower PM2.5 levels in August in Balkuy, a town in the same region as Ouagadougou in Burkina Faso. Ouarma et al.19 also found mean PM2.5 levels of 22 μg m−3 during the rainy season months and 87 μg m−3 during the dry season months in 2019 in Ouagadougou. The increased rainfall in August leads to more particle deposition and reduced urban activities and road dust, contributing to cleaner air.
The mean estimated PM2.5 concentrations in all areas of the city in the rainy season of 2000–2009 were lower (between 30 μg m−3 and 40 μg m−3) except at Gounghin, Ouagadougou Airport, Tanghin, Kamboinsi, and Kossodo where the PM2.5 concentrations were slightly above 40 μg m−3. Wet deposition in the rainy season removes PM2.5 particles from the atmosphere.19,53,56 Also, rainfall dampens the ground, reducing the amount of dust and other particles that can be resuspended into the air through wind or human activities.72 Additionally, the rainy season often coincides with increased vegetation growth, which can help trap particulates and prevent them from becoming airborne.73
In the dry season of the same period, the PM2.5 concentrations were between 65 μg m−3 and 70 μg m−3 except in Ouagadougou Airport, Patte d'Oie, Gounghin, Tanghin, Kamboinsi, Kossodo, and Tengadogo where the PM2.5 concentrations were slightly higher, between 70 μg m−3 and 75 μg m−3. The slightly higher mean PM2.5 concentrations in these areas in both seasons could depict emissions from vehicles. The high PM2.5 concentrations in every area of the city in the dry season are mainly due to the dust transported from the Sahara Desert by the Harmattan winds from November to March18,20 and also dust from unpaved roads.20 Also, frequent atmospheric stability caused by temperature inversions during the early dry season favors the accumulation of air pollutants.18 Additionally, during the dry season, there is an increase in open burning of household waste, further exacerbating air pollution levels.74 The mean estimated PM2.5 concentrations in the rainy season of 2012–2022 were high between 40 μg m−3 and 45 μg m−3 in all areas except in Gounghin, Kossodo, Ouagadougou Airport, Patte d'Oie, and Tanghin where the mean PM2.5 concentrations were between 50 μg m−3 and 55 μg m−3. The mean estimated PM2.5 concentrations in Yagma, Koumdanyore, Zagtouli, Ouaga 2000, Sakoula, and Roumtenga in the dry season of the same period were between 60 μg m−3 and 70 μg m−3 whereas in Ouagadougou Airport, Patte d'Oie, Gounghin, Tanghin, and Kossodo, the estimated PM2.5 concentrations were between 80 μg m−3 and 90 μg m−3. Ouagadougou Airport, Patte d'Oie, Kamboinsi, and Tanghin are home to many commercial activities with high volumes of traffic. Kossodo and Gounghin are the most established industrial zones in the city with several gas-oil power plants, cement and brick factories, metal processing plants, and textile manufacturing units.32 Additionally, these areas are frequented by heavy trucks fueled by diesel for transporting raw materials and finished products; hence, higher levels of PM2.5 in these areas. For Patte d'Oie, Kossodo, Ouagadougou Airport, Pissy, Tanghin, Dassasgo, Gounghin, and Roumtenga, estimated PM2.5 concentrations in the rainy season increased to about 10–20 μg m−3 in 2012–2022 from their levels in 2000–2009. For Kossodo, Gounghin, and Pissy, estimated PM2.5 concentrations in the dry season increased to about 10–15 μg m−3. PM2.5 concentrations in Yagma, Ouaga 2000, Sakoula, and Tengadogo in the dry season decreased to about 3–5 μg m−3 compared to their levels in 2000–2009. These decreases may be attributed to improvements in the number of paved roads in these areas.
Footnote |
† Electronic supplementary information (ESI) available. See DOI: https://doi.org/10.1039/d4ea00057a |
This journal is © The Royal Society of Chemistry 2024 |