**Spatiotemporal Assessment of Air Quality and Heat Island Effect Due to Industrial Activities and Urbanization in Southern Riyadh, Saudi Arabia**

**Abeer Salman 1,\*, Manahil Al-Tayib 2, Sulafa Hag-Elsafi 1, Faisal K. Zaidi <sup>3</sup> and Nada Al-Duwarij <sup>1</sup>**


<sup>3</sup> Department of Geology and Geophysics, King Saud University, Riyadh 11564, Saudi Arabia;

**\*** Correspondence: abalsalman@ksu.edu.sa

**Abstract:** The aim of this paper is to evaluate the air and thermal pollution in the southern suburbs of Riyadh, where people are suffering from poor air quality due to the rapid development of the industrial facilities in the area. The study involved the distribution of questionnaires to 405 residents living in that area in order to obtain their opinions about air pollution. In addition, land surface temperature and 12 components of air were measured at 18 points. In addition, the air quality data from 2016 to 2020 for Al Khaldya and Southern Ring Road air stations were assessed. Al Misfat (Oil Refinery) and the Second Industry City are significant contributors to air pollution and also result in the urban heat island effect from high temperature due to factories and industrial activities. However, all the measured components of air quality are lower than the standard limits except the element particulate matter (PM)10, which exceeds the standard limits in several parts of the study area and during several months of the year. This can surely have a negative impact on the health of residents. At the same time, this study found that the preventive measures taken to stop the spread of COVID-19 infections have led to a positive impact in the area in terms of improvement in air quality.

**Keywords:** air pollution; Riyadh; industrial city; land surface temperature; geostatistic

#### **1. Introduction**

Urbanization and industrialization are the main sources of air pollution in the world. Consequently, air pollution has become the main topic of many environmental studies [1–5]. The World Health Organization states that air pollution is a threat to human health and is responsible for millions of deaths annually [6–9]. Annually, about 21,000 Canadians die due to exposure to polluted air [10].

Several studies have shown a positive correlation between cancer and exposure to polluted air in different parts of the world [11,12] and in the Arabian Peninsula [13]. The industrial areas in the world are the main contributors to air pollution and are responsible for the deteriorating health of the population.

In Riyadh, the capital of Saudi Arabia, industrialization has increased dramatically during the past couple of decades in response to the increase in population and urban development. Establishment of two industrial cities has resulted in better employment opportunities and, consequently, a high population growth of about 4% in Riyadh [14]. The First Industrial City was established in 1973 and has an area of 0.5 km2. The Second Industrial City was established in 1976 and has an area of 19 km2 [15].

Though the Second Industrial City and the Saudi Aramco Oil Refinery (Al-Misfat) were established in the south of Riyadh, far from the city, urban expansion and associated economical activities forced the people to live around the industrial cities. Accordingly,

**Citation:** Salman, A.; Al-Tayib, M.; Hag-Elsafi, S.; Zaidi, F.K.; Al-Duwarij, N. Spatiotemporal Assessment of Air Quality and Heat Island Effect Due to Industrial Activities and Urbanization in Southern Riyadh, Saudi Arabia. *Appl. Sci.* **2021**, *11*, 2107. https://doi.org/10.3390/ app11052107

Academic Editor: Tung-Ching Su

Received: 8 January 2021 Accepted: 19 February 2021 Published: 27 February 2021

**Publisher's Note:** MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

**Copyright:** © 2021 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (https:// creativecommons.org/licenses/by/ 4.0/).

fzaidi@ksu.edu.sa

overcrowded buildings, especially along roads; an increase in traffic and vehicular emissions; and frequent desert dust storms in the area were the principal sources that increased the level of environmental pollutants [16], leading to air pollution [17]. The cumulative effect of these factors resulted in deterioration of air quality indoors and outdoors at many sites in Riyadh [18].

During several visits to the study area in southern Riyadh, it was evident that the area is suffering from multiple aspects of pollution, reduction in visibility during the mornings, and dirty and eroded buildings. Accordingly, residents and workers complained of unhealthy air quality in the area.

Several efforts have been made to study the concentration of chemicals and particulate matter in the air in Saudi Arabi [19–21]; however, most studies were limited and not comprehensive due to insufficient data [18]. As a result, the Royal Commission for Riyadh City and the General Authority of Meteorology and Environment Protection are operating 32 air quality control stations in Riyadh, over the past few years, to monitor the air quality. Some of these monitoring stations are fixed, whereas others are mobile. In addition, the General Authority of Meteorology and Environment Protection launched a dashboard for Riyadh air quality [22] to visualize the air quality index (AQI).

Mapping techniques are useful to assist decision makers in quantifying and locating pollution sites, even with limited measurements. (Geographic information system) GIS and remote sensing are useful tools for spatial analysis and detection of thermal pollution. Geostatistics is a strong technique to assess the spatial analysis of the studied variables [23]. The kriging technique is widely used for investigating sources of pollution as it is a reliable technique with fewer errors, even with small samples [24–26]. In addition, remote sensing imageries, especially in the thermal bands, provide simpler and more accurate information about surface temperature and emissivity [27] to detect urban heat islands (UHIs) [28–30].

Accordingly, this study aims to shed light on the opinions of residents living in southern Riyadh, near the industrial city, regarding emissions of oxides and particulate matter from industrial areas and traffic congestion. Most studies when surveying about the quality of life in cities do not take into account people's opinions about air pollution aspects. The spatial and temporal changes in thermal and air pollution were evaluated from 2016 to 2020. Additionally, the impact of the lockdown on air quality during the coronavirus pandemic was also assessed.

#### **2. Materials and Methods**

#### *2.1. Study Area*

The study area is located in the southern suburbs of Riyadh. It lies between the coordinates 24.482–24.594 North and 46.759–46.946 East (Figure 1). It covers about 128.7 km2 and includes the following seven districts: Ad-Difa, Al-Iskan, the second Industrial City, Al-Misfat, Al-Mansuriah, Ad-Dar Al-Baida, and Taybah.

The total number of inhabitants in the study area is 116,510 [31]. The Second Industrial City, including Al Iskan and Ad Dar Al Baida, has a higher population due to better infrastructure and amenities (Figure 2).

**Figure 1.** Study area location.

**Figure 2.** Population of the study area.

In the study area, two major industrial factories were established, namely the Second Industrial City and Al-Misfat (the Saudi Aramco Oil Refinery) (Figure 3). The Second Industrial City was established in 1979 and has more than 823 factories with different activities [15]. The Al-Misfat refinery was established in 1974 and produces various kinds of petroleum derivatives. Moreover, there is heavy vehicular traffic throughout the day to support the industrial area.

**Figure 3.** Land-use map of the study area and location of air sampling.

The western part of the study area is bounded by Wadi Hanifa, where the treated wastewater from Riyadh City is disposed. This water flows along the valley, with some agricultural activities in the Mansuriah and Taybah districts, which are less densely populated.

#### *2.2. The Questionnaire*

The study involved a survey about people's opinions regarding air pollution. The survey was conducted through a questionnaire, which was prepared in accordance with relevant guidelines and regulations and was approved by the Geography Department at King Saud University. Random samples of the population were chosen from the seven suburbs, and the questionnaire was distributed to 405 persons. The sub-samples, determined by the study, were according to population, as shown in Table S1.

The questionnaire consists of three parts: inhabitants' conditions, recognizing pollution, and pollution impact on the environment and health. The study applied chi-square and Cramer's V tests to understand the relationship between districts and pollution aspects.

#### *2.3. Spatial and Temporal Analysis*

Thermal analysis was carried out using the Landsat 8 satellite sensor image, and the land surface temperature (LST) was calculated from band 10 and band 11 by using land surface emissivity (ε). In addition, the normalized difference vegetation index (NDVI) was calculated using band 4 and band 5 to determine the proportion of vegetation (*PV* and ε) to estimate the LST [32–35], as shown in Figure 4. Three images from May 2014, 2017, and 2020, with a ground resolution of 30 m, were selected.

**Figure 4.** Flowchart to determine land surface temperature.

At the same time, the study involved the measurement of particulate matter (PM10, PM4, PM2.5, and PM1) and gases (CO, CO2, NO, NO2, SO2, O3, H2S, and Volatile organic compounds VOCs) using portable instruments at 18 points from 18 to 25 May 2014. The measuring was according to Presidency of Meteorology and Environment Protection standards (PME) in Riyadh and performed by Geotechnical and Environmental Company (GECO), as shown in Table S2. The distance between sample sites was 2–3 km, as shown in Figure 3.

The study used descriptive analysis for air sampling, to picture and summarize the data [36], and used correlation coefficients to assess the strength of linear relationship between variables. Principal component analysis (PCA) was carried out to transform the original set of variables into a smaller set of linear combinations [37].

Data analysis was carried out geostatistically, using ArcGIS 10.1, while ordinary kriging was done by creating a semivariogram. The study estimated the central tendency, dispersion, and shape to determine the type of distribution of numerical variables. A semivariogram was created by building the spatial variability structure of each attribute, as in shown in Equation (1).

$$\gamma(h) = \frac{1}{2N(h)} \sum\_{i=1}^{N(h)} \left[ (Z(\mathbf{x}) - Z(\mathbf{x} + h))^2 \right] \tag{1}$$

where *γ*(*h*) is the semi-variance, *Z*(*x*) is the value of initial potential at site *x*, *Z*(*x+h*) is the value of potential at site (*h*) distance apart from (*x*), and *N* is the number of sample pairs [38].

Ordinary kriging maps were established to predict values at a non-sampled point that assumes a constant unknown meaning, Equations (2a) and (2b):

$$Z(X\_0) = \mu + \varepsilon(X\_0) \tag{2a}$$

$$Z(X\_0) = \sum \lambda\_i \gamma(x\_i) \,, \quad \sum \lambda\_i = 1 \,\tag{2b}$$

where *μ* is an unknown constant, *ε*(*X*0) is the the error associated with an unknown location *X*0, *Z*(*X*0) is the estimated value of *Z* at *X*0, and *λ<sup>i</sup>* is the the weight that gives the best-possible estimation from the surrounding points.

Finally, the trend of the monthly average of NO2 and PM10 concentrations was determined for the Southern Ring Road air station from January 2016 to September 2018 and the Al-Khaldiyah air station from September 2019 to May 2020. These two air stations were selected as they are the closest to the study area and have air quality conditions similar to those in the study area. Moreover, a continuous data set for the time period mentioned above was available from these two air-quality-monitoring stations.

#### **3. Results**

#### *3.1. The Questionnaire*

Most respondents (nearly 64.85%) were in the age group of 20–40 years. Of these, 85.11% are employees (private sector or government), and their housing is near their workplace. Half of them are professionals and are living in the study area for more than 10 years.

Most of the respondents reported that the buildings they live in are dirty, the plants are dusty, and fine dusty clouds extending about 1–3 km or more are present around the residential buildings, as shown in Figure 5.

**Figure 5.** Noticing the pollution aspects in the study area.

They reported that most of the aspects related to air quality are similar in many residential districts of Riyadh; however, the problem of dirty buildings is specific to the study area, as shown in Table S3. Al-Iskan, Ad-Dar Al-Baida, and the Second Industrial

City are the most heavily polluted districts in the study area; and as a consequence, their residents suffer from many respiratory and eye diseases, as shown in Figure 6.

**Figure 6.** Diseases the population suffers from due to air pollution.

The survey results confirm the existence of problems related to air pollution in the study area. Field measurements and spatial analysis of gases and particulate matters, suspended in the air, are required to determine the extent and concentration of pollutants.

#### *3.2. Thermal Analysis*

The study widely used the NDVI as an indicator of the vegetation index and to calculate the LST. Generally, the study area shows a high variation in temperature. The highest temperature is observed in Al-Misfat, the factories, and the Second Industry City. The NDVI imageries from 2014 to 2020 (Figure 7) showed low values in most of the study area, except in the agricultural area in Taybah and Al-Mansuriah. The NDVI showed a gradual decrease in the study area, showing a steady increase in temperature. However, in May 2020, a lower temperature was recorded as much of the activities had stopped due to the lockdown to prevent the spread of COVID-19. This observation proves that human activities and reduced vegetation cover are major reasons for the development of heat islands in the region.

**Figure 7.** Land surface temperature and the normalized difference vegetation index (NDVI) (2014–2020).

#### *3.3. Statistical and Geostatistical Analysis*

The results of the descriptive statistics are given in (Table 1). The mean values of the particulate matter observed in the study area during most days of the year exceed the National Ambient Air Quality Standards (NAAQS).


**Table 1.** Descriptive statistics of air parameters.

PM, particulate matter; NAAQS, National Ambient Air Quality Standards.

The correlation coefficient shown in Table S4 shows strong-to-moderate correlation between particulates, CO2, CO, and NO2 and indicates the same source of pollution.

Factor analysis (PC) helped in identifying four factors having eigen values of >1. The four factors explain 82.046% of the total variance, as shown in Table S5. PC1 accounts for 41.2% of the data variability and shows significant factor loadings of PM1, PM2.5, PM4, PM10, CO, and NO2, as shown in Table S6. These variables have strong-to-moderate correlation, as shown in Table S4. PC2 represents 14.5% of the total variability and has significant factor loadings of SO2 and CO2. PC3 and PC4 represent 13.6% for O3 and VOCs and 12.756% for H2S and NO, respectively.

Ordinary kringing maps were built, as shown in Figure 8, for the particulate matter and gases, except O3 because there was no significant spatial variation, as seen in the semivargrom model shown in Table S7. The distribution of all variables was lognormal for different types of the chosen model, depending on the least mean square error.

High concentrations of pollutants and particulate matter are in Al-Misfat, the Second Industrial City, and along the highway.

The concentration of PM varies in the study area, depending on land use and working conditions. However, a high concentration of coarse particulate matter (PM10) covers large areas, and its mean (214 ug/m3) exceeds the standard limit (150 ug/m3). Fine particulates (smaller than 2.5) are more critical and may cause more respiratory problems. They are concentrated in specific areas around Al-Misfat and the Second Industrial City. These fine particulates, derived mostly from combustion, may remain suspended for weeks and drift for many kilometers [39]. Moreover, spatial distributions of particulates are similar to those of CO and NO2 and consequently considered to be within PC1.

A high variation in the concentration of oxides, CO2, NO2, and SO2, is observed in the study area. The highest concentrations are around the Second Industrial City and Al-Misfat, where industries and automobile emissions are considered as the major sources of NO2 [20]; however, emissions remained below standard at all places. VOC emission is concentrated in Al-Misfat and in the southeastern part of the Second Industrial City. The oil and natural gases sector is the main source of VOCs. H2S is concentrated in the urban area in the northern part of the study area and near Al-Misfat. The main sources of H2S are petroleum refineries, Al-Misfat, and the sewage system. Ad-Dar Al-Baida suffers from the leakage of an old sewage system, and the population suffers greatly from stagnant sewage pools and bad smells.

PM10 and NO2 concentrations showed a continuous increase from 2016 to 2019. However, they clearly decreased during the spring of 2020 by 62% and 78%, respectively, due to the lockdown imposed during the coronavirus pandemic, as shown in Figure 9, wherever human activity decreased dramatically.

The southern part of the study area suffers from pollution, mainly particulate matter from factories and dust storms. In addition, vehicles are one of the main causes of air pollution and emissions of carbon from PM2.5 [21]. Ad-Dar Al-Baida, Taybah, and Al-Mansuriah are the least affected due to their distance from the industrial area. Despite that, Ad-Dar Al-Baida residents complain of pollution, and remeidial measures must be adopted to reduce the seriousness of the situation. More studies should be conducted in the study area to determine the general direction of the concentrations of oxides and to identify the chemical composition of the particulate matter.

Moreover, citizens may change their lifestyle to mitigate the sources of air pollution, such as transportation and energy consumption [40].

**Figure 8.** *Cont*.

**Figure 8.** Ordinary kriging map of particulate matter and gases.

**Figure 9.** Trend in PM10 and NO2 during 2016–2020 on the Southern Ring Road and the Al Khaldiyah air station.

#### **4. Conclusions**

This study based on the questionnaire responses of residents living in southern Riyadh and air quality monitoring from 2016 to 2020 from a couple of air-monitoring stations close to the area of interest reveals that the adverse effects of poor air quality due to urbanization and industrialization are mostly felt by the residents of Al-Iskan, Ad-Dar Al-Baida, and the Second Industrial City. Respiratory problems and irritation in the eyes were reported by 49% and 50% of the residents, respectively, during the survey. A high concentration of PM10 has been observed in many places in the study area, and it exceeds the permissible limits. Thermal analysis of the area showed the maximum temperature in the Second Industrial City and Al-Misfat, indicating the presence of the thermal island effect in these regions due to industrial activities. The anthropogenic influence on air quality was evident from the air quality data for 2020, which showed a marked decrease in parameters such as NO2 due to the impacts of lockdown on various activities during the coronavirus pandemic. The land surface temperature image for May 2020 showed a temperature range between 54 ◦C and 32.4 ◦C, which is less than the temperature recorded during the same period in 2014 and 2017 and clearly shows the positive impact of minimal commercial/industrial activity during the lockdown period in 2020.

**Supplementary Materials:** The following are available online at https://www.mdpi.com/2076-341 7/11/5/2107/s1, Table S1: Distribution of questionnaire in the study area; Table S2: Air sampling methods and instruments; Table S3: Relationship between districts and impact of pollution; Table S4: Correlation coefficients at 95% confidence level; Table S5: Eigen values of PC of the data sets; Table S6: Principal Component loadings of the data sets; Table S7: Semivariogram parameters.

**Author Contributions:** A.S. wrote the main manuscript and interpreted the data; M.A.-T. prepared questioners and analyzed the statistical data; S.H.-E. analyzed the Landsat images; F.K.Z. drafted and editing the work; N.A.-D. designed maps and prepared Figures 1–3. All authors have read and agreed to the published version of the manuscript.

**Funding:** This research was funded by the Deanship of Scientific Research at King Saud University through research group number RG-1438-023.

**Institutional Review Board Statement:** Not applicable.

**Informed Consent Statement:** Not applicable.

**Data Availability Statement:** Data sharing not applicable.

**Acknowledgments:** Sincere thanks to the Royal Commission for Riyadh City and the General Authority of Meteorology and Environment Protection for providing air quality data.

**Conflicts of Interest:** The authors declare no conflict of interest.

#### **References**


#### MDPI

St. Alban-Anlage 66 4052 Basel Switzerland Tel. +41 61 683 77 34 Fax +41 61 302 89 18 www.mdpi.com

*Applied Sciences* Editorial Office E-mail: applsci@mdpi.com www.mdpi.com/journal/applsci

Academic Open Access Publishing

www.mdpi.com ISBN 978-3-0365-8279-5