116
Views
0
CrossRef citations to date
0
Altmetric
Research Article

The effect of coal-fired power plants on ambient air quality in Mpumalanga province, South Africa, 2014–2018

ORCID Icon, ORCID Icon, ORCID Icon & ORCID Icon
Received 29 Nov 2023, Accepted 29 Apr 2024, Published online: 06 May 2024

ABSTRACT

Several coal-fired power plants (CFPPs) were built in South Africa, mainly in the central Mpumalanga Province, due to an increase in the demand for Eskom, the national power utility, to keep up with socio-economic growth. The CFPPs, of which 90% are owned by Eskom, generate a significant share of the country’s electricity but contribute to the air pollution experienced in the country. The paper discusses sulphur dioxide (SO2), nitrogen dioxide (NO2) and particulate matter of size less than 10 micrometre (μm) in diameter (PM10), using data from 2014 to 2018. The statistics revealed higher PM10 concentrations during winter than in summer and spring at the Kriel and Komati sites; associated with the higher contribution of domestic burning. The study’s results could influence legislation and policies and help to understand the source of poor ambient air quality by assessing the three pollutants within the area of the selected power plants.

Introduction

Air quality deterioration has become an increasing global concern due to its impacts on public health and the environment from the Industrial Revolution in the 18th century to modern times (Olufemi et al. Citation2018; Xiong and Xu Citation2021; Zhang et al. Citation2022). Today, industries require considerably more energy compared to the 1800s. This compounded with industry’s heavy reliance on fossil fuels, which significantly contributed to poor ambient air quality (Naidoo et al. Citation2014; Moletsane et al. Citation2021; Walton et al. Citation2021). Several coal-fired power plants (CFPPs) were built in South Africa (SA), mainly in the central Mpumalanga Province, due to an increase in the demand for the national power utility, Eskom, to keep up with socio-economic growth in the country (Rorich and Galpin Citation1998). These CFPPs are mostly (90%) owned and operated by the power utility, and 80% are located within Mpumalanga.

The National Environmental Management: Air Quality Act (NEMAQA) (No. 39 of 2004) air quality regulations provide the minimum emission standards (MES) for CFPPs, regulating the pollution that may be emitted from point sources, such as stacks. Compliance with existing and new emission limits was specified in 2015 and 2020, respectively, depending on the date of construction (DoE Department of Energy Citation2019). The emission reduction initiatives in SA are also in line with the 2015 adoption of the United Nations’ 2030 Agenda for Sustainable Development Goals (SDGs), which the country has signed as outlined in the study by Udeagha and Muchapondwa (Citation2023). To reduce the impact of air pollution due to CFPPs, the power utility has developed an Emission Reduction Plan, focusing on reducing SO2, among other pollutants (Eskom Citation2021). Particulate emission reduction initiatives have been implemented at some of the CFPPs to reduce emissions at the existing coal-fired fleet (as seen in –Matimolane et al. Citation2021), such as the retrofit of electrostatic precipitators (ESPs), flue gas desulphurisation (FGD) and fabric filter plants (FFPs). Currently, the FGD is used only at the Kusile Power Plant to reduce SO2 emissions, although it is currently not fully operational due to infrastructure breakdowns. The retrofitting with FGD at Medupi Power Plant will be implemented in 2025.

Figure 1. Timescale in relative emissions of particulates (ash) from Eskom’s power plants from the early 1980s to present (Matimolane et al. Citation2021).

Figure 1. Timescale in relative emissions of particulates (ash) from Eskom’s power plants from the early 1980s to present (Matimolane et al. Citation2021).

The power utility, in compliance with its atmospheric emission license (AEL), has strategically installed ambient air quality stations to attain continuous air pollution parameter measurements to represent the ambient air of the areas potentially impacted by the power plant emissions (Eskom Citation2013). The measuring stations are distributed and located in residential areas to measure and understand the impact of air pollution on health (Stolz et al. Citation2020) and the environment in general. The burning of solid fuels is the primary source of domestic energy for people living in low-income settlements in developing countries, which causes gaseous and particulate emissions (Xie et al. Citation2015).

Appreciable reductions in air pollutant concentrations and emission factors from CFPPs in China, India, Japan, Europe and the USA were observed over the past two decades (Wang et al. Citation2020). The combination of control measures such as the application of removal technologies, stringent emission standards, equipment improvement and the new perspective of policymaking as well as effectiveness of government policies to reduce emissions (Filonchyk and Peterson Citation2023a; Wan et al. Citation2023). However, from the listed control measures, technological upgrades are the most essential for addressing further emission control issues to improve air quality (Weng et al. Citation2023). The current study is of significance since there is continued pressure to move away from fossil fuel power generation. Furthermore, the assessment and understanding of the gaseous and particulate emission impacts from CFPPs relates to how the changes in air quality are managed.

The aim of the current study was to analyse the air quality data from 2014 to 2018 obtained from the power utility to determine ambient air quality trends in SO2, NO2 and PM10 from each community neighbouring the selected CFPPs. The present study, therefore, selected one of the previously mothballed CFPPs, namely the Komati Power Plant, which was returned to service in 2013 and one large 3000 MW Kriel Power Plant. The data from ambient monitoring stations near these two CFPPs were used to assess compliance of ambient pollution concentrations with the National Ambient Air Quality Standards (NAAQS), to determine ambient air quality trends of these pollutants and, lastly, to determine the correlation between the three selected criteria pollutants.

Materials and methods

Study area

Two ambient air quality monitoring sites were selected within a monitoring network of 600 km2 as observed in and were considered to be representative of the study area. These monitoring sites are located at the Kriel and the Komati villages corresponding to the two CFPPs selected for the current study. The Kriel village site commenced in June 2006 and the monitoring site, as observed in , is located ±7.8 km east of this power plant. The Komati site () is situated 2.2 km west-south-west of the Komati Power Plant and monitoring commenced on 1 June 2006 (Eskom Citation2013). The sites represented a range of highly industrialised areas such as Secunda, where the Kriel village monitoring site is located, to the more rural area of Komati, where the Komati site is located. The air pollution sources in the area include a nearby coal mining and handling company, along with haul truck exhaust emissions; domestic fuel use; agriculture and biomass burning as well as small industries. Both sites are affected by CFPPs and, further afield, a large petrochemical refinery.

Figure 2. Map of the South African Highveld including major towns, power stations and monitoring locations (Lourens et al. Citation2011)

Figure 2. Map of the South African Highveld including major towns, power stations and monitoring locations (Lourens et al. Citation2011)

Data set description

The Kriel village and the Komati sites are equipped to continuously monitor ambient concentrations of SO2, NO2, NOx, and coarse particulate matter of particulate size less than 10 μm in diameter (PM10). The standard specifications, as well as the equipment/techniques used for the measurement of SO2 and NOx, conform to the US-EPA equivalent method numbers EQSA-0486-060 and RFNA-1289-074, respectively (Komane Citation2019). The meteorological parameters of wind velocity, wind direction, ambient pressure, ambient temperature, humidity, rainfall and solar radiation are also recorded but not analysed in this discussion.

Data analysis

Various descriptive statistics were calculated from hourly (where possible), daily, weekly, monthly and annual measures of ambient SO2, NO2 and PM10. The reason for focusing on these three pollutants is because of their cosmopolitan nature and the large quantities released from the coal combustion process for the generation of electricity, as well as for industrial and domestic use. The power utility’s in-house personnel from Research, Testing and Development (RT&D) provided the SO2, NO2 and PM10 data used in the current study. Statistical analyses were performed using R (version 4.1.0), and the R Openair package version 2.11 (https://bookdown.org/david_carslaw/openair/). The R statistical package allows for easy importing of the data and descriptive analysis of the data (R Core Team Citation2021). Note that the data integrity was accepted since it was assumed that it was already completed by the RT&D team. Outliers were not removed, but those were incomplete data for the pollutants at various timestamps that were removed in order to complete the statistical analyses. Incomplete data were due to plant unavailability reported as the percentage of time that electrical power was unavailable to the monitoring station. In addition, the analysers installed during mid-month sometimes became faulty and were then removed from the site for maintenance. Considering hourly data, the percentages of missing values were about 49% for the Kriel village site (from a total of 37,944 observations), and 43% for the Komati site (from a total of 29,928 observations). The measured data from 2014 to 2018 were used to assess compliance of ambient pollution concentrations with the NAAQS to determine ambient air quality trends and to assess correlations between the three ambient pollutants. In this paper, the Theil-Sen analysis was used to estimate the trends of the pollutants for the given period as it tends to yield accurate confidence intervals and is robust against outliers as it is based on the median of the slopes (Carslaw and Ropkins Citation2012; Munir et al. Citation2013). Note that this study used the 95th percentile as an index that captures the range in the distribution of daily observed criteria pollutants selected during the study period.

Results and discussion

The results are presented for the three criteria pollutants of greatest concern in SA, namely SO2, NO2 and PM10, concerning the environment in general. The raw data for annual variations of ambient air pollutants were visualised with time-series plots for each station for the full five-year period for the Kriel village and the Komati sites, as shown in supplemental material, i.e. Figures A1 and A2 for PM10, NO2 and SO2. Note that in supplemental Figures A1 and A2, the pollutant is specified on each y-axis, and the y-axis range is in the units of the indicated pollutant (μg.m−3), with its minimum and maximum values specified.

Comparison with ambient air quality standards

National Ambient Air Quality Standards (NAAQS) have been established for criteria pollutants and are a combination of a limit for an average period and an allowable frequency of exceedance (supplemental material: Table A1). The study used verified daily concentrations over annual periods for each pollutant and at each site over the investigated period, and the NAAQS limit exceedances are highlighted in bold, as observed in supplemental Table A2. The PM10 levels exceeded the NAAQS 24-hour standard of 75 μg.m−3 at both the Kriel village and the Komati sites in all the investigated years when considering the daily 95th percentile pollutant concentrations over annual periods. The significance of mentioning the PM10 exceedance at both selected sites is mainly due to the domestic burning that increases in winter and therefore, the PM10 exceedance is not only due to the power station emissions (Lourens et al. Citation2011, Vente and Lourens Citation2021). In addition to the domestic burning during colder months are the veld fires that happen due to the dry conditions within the environment (Venter et al. Citation2012, Vente and Lourens Citation2021). Note the NAAQS for PM10’s one-year limit of 40 μg.m−3 for the 2014 data compared to 50 μg.m−3 for the 2015 to 2018 data, observed in supplemental material in Table A1.

Sulphur dioxide

The daily mean SO2 concentrations over annual periods at the two selected monitoring sites are presented in supplemental Table A2. The ambient SO2 concentrations for 2014 to 2018 were low at both sites with daily mean values in the range of 11.0 to 15.4 μg.m−3 at the Kriel village site and 11.4 to 21.0 μg.m−3 at the Komati site, during the investigated period. The daily mean SO2 concentrations over annual periods at all the sites did not exceed the SO2 24-hour NAAQS limit of 125 μg.m−3 as observed in supplemental Table A2. For SO2, the maximum concentration range was 79.1 μg.m−3 at the Kriel village site during the investigated period. There were zero (0) exceedances of the SO2 daily limit at both the Kriel village and the Komati sites. Levels of SO2 in urban areas have been linked to coal combustion in industries and the burning of coal is the main contributor of SO2 into the atmosphere (Xu et al. Citation2017; Morakinyo et al. Citation2020).

Nitrogen dioxide

The daily mean NO2 concentrations over annual periods at the monitoring sites are presented in supplemental Table 2. The ambient NO2 concentrations over the period were low with daily mean values in the range of 7.2 to 12.9 μg.m−3 at the Kriel village site and 8.0 to 12.6 μg.m−3 at the Komati site. The maximum daily mean NO2 concentration was 68.7 μg.m−3 at the Kriel village site during the investigated period. There are no 24-hour NAAQS limits for NO2 according to the South African NAAQS established for criteria pollutants (Government Gazette Republic of South Africa Citation2009) as shown in supplemental Table A1.

Particulate matter less than 10 micrometres in diameter

The daily mean PM10 concentrations over annual periods at the monitoring sites are presented in supplemental Table A2. The ambient PM10 concentrations over the period were moderate to high with mean values in the range of 36.8 to 60.9 μg.m−3 at the Kriel village site and 53.1 to 70.8 μg.m−3 at the Komati site. The maximum mean PM10 concentrations over the annual period were 248.1 μg.m−3 at the Komati site during the investigated period. At the Kriel village site, there were 236 exceedances of the PM10 daily limit, while at the Komati site, there were 423 exceedances during the entire monitoring period. The significance of mentioning these PM10 exceedances besides domestic burning and veld fires is that PM10 is especially problematic during dry season since rain not only scrub the atmosphere but also assists in suppression; hence, an increase in PM10 is typically observed between May and September annually (Lourens et al. Citation2011; Venter Citation2021).

Monthly, seasonal, weekly and daily variations of ambient air pollutants

The daily 95th percentile PM10 concentration monthly average observed in 2017 at the Kriel village site of 103.6 mg.cm−3 was among the five out of 10 lowest concentrations at both sites over the years, as observed in the box-and-whisker plot in . In contrast, the Komati site had the highest daily 95th percentile PM10 concentrations per month of 144.0 mg.cm−3 recorded in August 2017 (winter), as observed in the box-and-whisker plot in .

The Theil-Sen graphs of weekly mean PM10 concentrations have demonstrated decreasing trends at the Kriel village site of −0.67 mg.cm−3.y−1, with a 95% confidence interval of −4.72,3.73 mg.cm−3.y−1, as observed in . Note that if the 95% confidence interval does not contain zero then the estimated slope is statistically significant at the 0.05 level, i.e. significantly different from zero with 95% confidence. The daily mean PM10 concentrations at the Kriel village site have an estimated slope of −0.36 (−5,4.12) mg.cm−3.y−1, as observed in . However, the Theil-Sen analysis indicated increasing trends in weekly mean PM10 concentrations at the Komati site of 3.1 (−1.92,6.47) mg.cm−3.y−1 and daily mean PM10 concentrations of 3.17 (−2.34,5.79) mg.cm−3.y−1, as observed in . Despite the increase in PM10 emissions from the CFPPs, there were no temporal increasing trends in PM10 concentrations at the Kriel village site. Similarly, the daily SO2 mean and weekly concentrations at the Kriel village site shown by the Theil-Sen plots in compared to the SO2 at the Komati () site do not reflect, among other things, the sulphur contents of coal that have increased due to the poor coal quality procured by the power utility through its Primary Energy department (Eskom Citation2021).

Figure 3. The box-and-whisker plot of PM10 at the Kriel village (a) and Komati (b) sites showing the year and month with lowest and highest pollutant values.

Figure 3. The box-and-whisker plot of PM10 at the Kriel village (a) and Komati (b) sites showing the year and month with lowest and highest pollutant values.

Figure 4. The Theil-Sen plots of PM10 at the Kriel village (a) and Komati (b) sites showing trends in the daily and weekly mean pollutant concentrations.

Figure 4. The Theil-Sen plots of PM10 at the Kriel village (a) and Komati (b) sites showing trends in the daily and weekly mean pollutant concentrations.

The Theil-Sen plot showed that the weekly and daily mean NO2 concentrations have been decreasing at both the Kriel village and the Komati sites over the five-year period of assessment, as observed in . This could be attributed to very small changes in gaseous and particulate emissions or as a consequence of weather conditions that were not considered in the current study. Although the meteorological data was recorded but not discussed in the current article, TROPOspheric Monitoring Instrument (TROPOMI) has the capability to map a multitude of trace gases such as NO2, ozone, formaldehyde, SO2, methane, carbon monoxide and aerosols (Filonchyk and Peterson Citation2023a, Citation2023b). The latter would have affected PM10, and therefore, it is important to be aware of this aspect for all three pollutants considered in the current study. Moreover, in Kumar et al. (Citation2021) NOx concentration was high from vehicular sources and PM10 concentration was high from industrial sources at ambient concentration when using meteorological data. In the current study, the higher concentrations of PM10 primarily occurred in winter (especially August), with domestic burning a major contributor.

Figure 5. The Theil-Sen plots of SO2 at the Kriel village (a) and Komati (b) sites showing trends in the daily and weekly mean pollutant concentrations.

Figure 5. The Theil-Sen plots of SO2 at the Kriel village (a) and Komati (b) sites showing trends in the daily and weekly mean pollutant concentrations.

Figure 6. The Theil-Sen plots of NO2 at the Kriel village (a) and Komati (b) sites showing trends in the daily and weekly mean pollutant concentrations.

Figure 6. The Theil-Sen plots of NO2 at the Kriel village (a) and Komati (b) sites showing trends in the daily and weekly mean pollutant concentrations.

Interestingly, in winter and spring seasons the surface CO2 footprints are spatially expanded as a result of contributions of emissions from biomass and domestic fossil fuel combustion (Ncipha and Sivakumar Citation2022). Levels of NO2 and SO2 in urban areas have been linked to coal combustion in industries and biomass burning in residential areas (Morakinyo et al. Citation2020). In the current study, the NO2 concentrations peaked in spring and winter (supplemental material A1) with no significant trends otherwise. The increase in the number of vehicles in the country has contributed to changes in the emissions profiles, including an increase in vehicle emissions (Pretorius et al. Citation2015). The NO2 concentrations during the investigated period were variable from month to month, with lower concentrations during January, April and October as shown in supplemental material A1a (Kriel) and b (Komati). The average concentrations of NO2, SO2 and PM10 were in the order of winter > autumn > spring > summer. These seasonal variations reflect the effects of meteorological conditions occasioned by dawdling winds and a greater percentage of wind calm hours; lower mixing of boundary heights, limiting the diffusion capacity of the atmosphere and thus confining the emitted pollutants close to the ground, less rainfall, and lower temperatures (Xu et al. Citation2017; Morakinyo et al. Citation2020).

Correlation analysis between ambient air pollutants

The R Openair package version 2.11 (R Core Team Citation2021) was used to perform a correlation analysis between measured pollutants to understand the linear relationship among air pollutants as observed in supplemental Figures B1 to B10. Moderate linear correlations in data indicate that there are potential linear relationships between the three pollutants, but further statistical analyses are required to quantify the precise relationships between the pollutants. A correlation percentage of 20% to 40% in absolute value is considered to be a weak correlation. Correlations between 40% and 60% in absolute value are considered moderate, while those below 20% in absolute value are considered very weak. Correlation percentages whose magnitudes are between 60% and 80% in absolute value indicate variables that can be considered strongly correlated. The correlation percentage of 80% and above in absolute value is considered very strong.

The linear correlations between PM10, NO2, and SO2 were mostly moderate to high at all sites for all the years considered, as observed in supplemental Figures B1 to B10. Possible reasons for moderate-to-high linear correlations observed in this study’s pollutant data might be from coal combustion as a source of SO2, NO2 and PM10. Agriculture and wind-blown dust, on the other hand, consist of only particulate matter. Additionally, the lack of flue gas desulphurisation (FGD) in all the selected stations near the monitoring sites might be a result of the atmospheric chemistry driving the conversion of SO2 to secondary particulates (Eskom Citation2021).

Results showed that the concentrations of SO2 and NO2 changed with seasons, and they manifested the same trend, i.e. correlations in winter are greater than spring, which are greater than in autumn, which are greater than summer. In contrast, at the Komati site, in 2015 and 2016, there were very strong positive correlations between NO2 and SO2 in summer of 81% and 90%, respectively, as observed in supplemental Figures B1 to B10. Noticeably, a very strong linear correlation of 82% between NO2 and PM10 in September was also observed at the Komati site in 2017. However, there was a very weak negative correlation of −2% between SO2 and PM10 at the Komati site in the spring of 2016. In summary, from the current study, the main seasons of concern to control SO2 and NO2 were identified as winter and spring.

Conclusion

In the current study, the three criteria pollutants, namely SO2, NO2 and PM10 were analysed from the hourly (where possible), daily, weekly, monthly and annually monitored data obtained for the period 2014 to 2018. The results showed that only PM10 measured regularly exceeded the annual NAAQS at the Kriel village and the Komati sites when both daily mean and 95th percentile PM10 concentrations over annual periods were considered. The Theil-Sen and box-and-whisker plots revealed higher PM10 concentrations and higher exceedances during winter (especially August), and lower concentrations during summer in all the selected years, with the highest concentrations reached in 2017 at the Komati site. However, this trend could not be attributed to the operational activities at both power plants alone, as there are other CFPPs, various mines, larger industries, and other agricultural activities near these monitoring sites. Moreover, the linear correlations between PM10, NO2, and SO2 were moderate to high in all the investigated years and selected sites. Understanding such interactions is essential for urban climate studies and our study provides a basis for legislation and policies to estimate possible changes in ambient air quality for these two CFPPs to fully comply with the NEMAQA Act.

The decrease in the emissions of the three pollutants in 2018 compared to previous years, could be explained by the retrofitting of plants with pollution abatement technologies in meeting the 2015 MES target deadline (Pretorius et al. Citation2015). The ambient monitoring stations are helpful in showing the impacts of several other activities directly and indirectly associated with the operation of CFPPs. These should ideally be supplemented with either dispersion plume (Morakinyo et al. Citation2020) or Bayesian model approaches geared at assessing the impact of the CFPP operations using historical and current data (Ngamlana et al. Citation2024). The daily PM10 over annual period concentrations exceedances have shown that the study’s results can influence legislation and policies to estimate possible changes in particulate emissions.

Acknowledgements

This is contribution number 883 of the North-West University (NWU) Water Research Group.

Disclosure statement

No potential conflict of interest was reported by the author(s).

Data availability statement

Data will be made available on request.

Additional information

Funding

The author(s) reported there is no funding associated with the work featured in this article.

References

  • Carslaw D, Ropkins K. 2012. Openair—an R package for air quality data analysis. Environ Modell Software. 27–28:52–61. doi: 10.1016/j.envsoft.2011.09.008.
  • (DoE) Department of Energy. 2019, The South African energy sector report. Directorate: energy data collection, management and analysis. Pretoria, South Africa. Online: 2019-South-African-Energy-Sector-Report.pdf.[accessed 2023 August 15] DoE.
  • Eskom. 2013. Eskom brings mothballed stations back to life. [accessed 2020 Apr 8]. https://m.engineeringnews.co.za/article/eskom-brings-mothballed-stations-back-to-life-2013-12-13/rep_id:4433 Eskom.
  • Eskom. 2021. Eskom integrated report, 31 March. https://www.eskom.co.za/investors/integrated-results/[accessed 2022 February 03].
  • Filonchyk M, Peterson MP. 2023a. An integrated analysis of air pollution from US coal-fired power plants. Geosci Front. 14(2):101498. doi: 10.1016/j.gsf.2022.101498.
  • Filonchyk M, Peterson MP. 2023b. NO2 emissions from oil refineries in the Mississippi Delta. Sci Total Environ. 898:165569. doi: 10.1016/j.scitotenv.2023.165569.
  • Government Gazette Republic of South Africa. 2009. Dec 24. (32816), last access: 16 Aug 2023. http://faolex.fao.org/docs/pdf/saf122986.pdf.
  • Komane A. 2019. Eskom kriel village air quality report. Eskom research testing & development. Johannesburg: Eskom Holdings SOC.
  • Kumar A, Dikshit AK, Patil RS. 2021. Use of simulated and observed meteorology for air quality modeling and source ranking for an industrial region. Sustainability. 13(8):4276. doi: 10.3390/su13084276.
  • Lourens AS, Fourie GD, Burger JW, Pienaar JJ, Read CE, Jordaan JH, Van Zyl PG, Beukes JP. 2011. Spatial and temporal assessment of gaseous pollutants in the highveld of South Africa. S Afr J Sci. 107(1):1–8. https://hdl.handle.net/10520/EJC97109.
  • Matimolane M, McCourt B, Herbst D. 2021. Air quality offsets implementation plan for Gert Sibande district municipality: tutuka, majuba, Camden and grootvlei power stations: march 2021 update. ( No. ENV21-R052). Sunninghill: Eskom Holdings SOC Ltd.
  • Moletsane SD, Lindeque F, Language B, Nkosi NC, Adesina JA, Burger RP, Mkhatshwa G, Piketh SJ. 2021. Intra-urban variability of PM25 in a dense, low-income settlement on the South African highveld. CAJ. 31(1):1–9. doi: 10.17159/caj/2021/31/1.9413.
  • Morakinyo OM, Mukhola MS, Mokgobu MI. 2020. Concentration levels and carcinogenic and mutagenic risks of PM 2.5-bound polycyclic aromatic hydrocarbons in an urban–industrial area in South Africa. Environ Geochem Health. 42(7):2163–2178. https://link.springer.com/article/10.1007/s10653-019-00493-2.
  • Munir S, Habeebullah TM, Seroji AR, Gabr SS, Mohammed AM, Morsy EA. 2013. Quantifying temporal trends of atmospheric pollutants in Makkah (1997–2012). Atmos Environ. 77:647–655. doi: 10.1016/j.atmosenv.2013.05.075.
  • Naidoo S, Piketh S, Curtis C. 2014. Quantification of emissions generated from domestic burning activities from townships in Johannesburg. CAJ. 24(1):34–40. doi: 10.17159/caj/2014/24/1.7047.
  • Ncipha XG, Sivakumar V. 2022. The first national scale spatial and temporal analysis of surface CO2 over South Africa utilising satellite retrievals. S Afr Geogr J. 104(2):137–154. doi: 10.1080/03736245.2021.1934093.
  • Ngamlana Z, Malherbe W, Gericke G. 2024. Association of coal fired power plants with river water quality in South Africa. Limnologica. 104:8–9. doi: 10.1016/j.limno.2023.126140.
  • Olufemi AC, Bello PO, Mji A. 2018. Conflict implications of coal mining and environmental pollution in South Africa: lessons from niger delta, Nigeria. Afr J Conflict Resolut. 18(1):7–35. https://hdl.handle.net/10520/EJC-ff5423180.
  • Pretorius I, Piketh S, Burger R, Neomagus H. 2015. A perspective on South African coal fired power station emissions. J Energy South Afr. 26(3):27–40. doi: 10.17159/2413-3051/2015/v26i3a2127.
  • R Core Team. 2021. R: a language and environment for statistical computing. Vienna, Austria: R Foundation for Statistical Computing. https://www.R-project.org/.
  • Rorich RP, Galpin JS. 1998. Air quality in the mpumalanga highveld region, South Africa. S Afr J Sci. 94:109–114. https://journals.co.za/doi/abs/10.10520/AJA00382353_133.
  • Stolz T, Huertas ME, Mendoza A. 2020. Assessment of air quality monitoring networks using an ensemble clustering method in the three major metropolitan areas of Mexico. Atmos Pollut Res. 11(8):1271–1280. doi: 10.1016/j.apr.2020.05.005.
  • Udeagha MC, Muchapondwa E. 2023. Striving for the united nations (UN) sustainable development goals (SDGs) in BRICS economies: the role of green finance, fintech, and natural resource rent. Sustainable Dev. 31(5):3657–3672. doi: 10.1002/sd.2618.
  • Venter AD, Beukes JP, Van Zyl PG, Tiitta P, Josipovic M, Pienaar JJ, Kulmala M. 2012. An air quality assessment in the industrialised western Bushveld Igneous Complex, South Africa. S Afr J Sci. 108(9):1–10. doi: 10.4102/sajs.v108i9/10.1059.
  • Venter AD, Lourens AS. 2021. Ambient air quality data reported at Sasol Secunda monitoring stations during COVID-19 lockdown-Mpumalanga, South Africa. Clean Air Journal. 31(1):1–7. doi: 10.17159/caj/2021/31/1.8954.
  • Walton NM, Piketh SJ, van Zyl P, Maenhaut W, Burger R, Language B, Formenti P. 2021. Source apportionment of ambient fine and coarse aerosols in embalenhle and kinross, South Africa. CAJ. 31(2):1–13. doi: 10.17159/caj/2021/31/2.11980.
  • Wang G, Deng J, Zhang Y, Zhang Q, Duan L, Hao J, Jiang J. 2020. Air pollutant emissions from coal-fired power plants in China over the past two decades. Sci Total Environ. 741:1–9. doi: 10.1016/j.scitotenv.2020.140326.
  • Wan R, Tang L, Guo J, Zhai W, Li L, Xie Y, Bo X, Wu J. 2023. Cost-benefits analysis of ultra-low emissions standard on air quality and health impact in thermal power plants in China. J Environ Manage. 345:1–10. doi: 10.1016/j.jenvman.2023.118731.
  • Weng Z, Song Y, Cheng C, Tong D, Xu M, Wang M, Xie Y. 2023. Possible underestimation of the coal-fired power plants to air pollution in China. Resour Conserv Recycl. 198:1–10. doi: 10.1016/j.resconrec.2023.107208.
  • Xie Y, Zhao B, Zhang L, Luo R. 2015. Spatiotemporal variations of PM2.5 and PM10 concentrations between 31 Chinese cities and their relationships with SO2, NO2, CO and O3. Particuology. 20:141–149. doi: 10.1016/j.partic.2015.01.003.
  • Xiong J, Xu D. 2021. Relationship between energy consumption, economic growth and environmental pollution in China. Environ Res. 194:1–8. doi: 10.1016/j.envres.2021.110718.
  • Xu Y, Ren C, Ma P, Ho J, Wang W, Lau KKL, Lin H, Ng E. 2017. Urban morphology detection and computation for urban climate research. Landscape Urban Plann. 167:212–224. doi: 10.1016/j.landurbplan.2017.06.018.
  • Zhang X, Han L, Wei H, Tan X, Zhou W, Li W, Qian Y. 2022. Linking urbanization and air quality together: a review and a perspective on the future sustainable urban development. J Cleaner Prod. 346:1–14. doi: 10.1016/j.jclepro.2022.130988.