Comparison of Metropolitan Cities for Mortality Rates Attributed to Ambient Air Pollution Within the Context of SDGs Using the AirQ Model

Abstract

In the present study, the air pollution dynamics of the metropolitan cities of Balıkesir, Bursa, İstanbul, Kocaeli, Sakarya and Tekirdağ in the Marmara Region, which is the geographical region with the highest urban and industrial activity in Turkey, were examined for the time period between 2016 and 2019. Annual changes in the cities in terms of air pollution, which was examined with a focus on the PM_2.5 parameter as indicated by United Nations (UN) Sustainable Development Goals (SDGs), differences in the cities by years, and the seasonal changes in air pollution in the cities were investigated. Additionally, mortality rates attributed to air pollution were calculated with the AirQ+ software based on integrated exposure-response function recommended by the World Health Organization (WHO) and the UN using city-scale statistics of fatal disease cases that can be attributed to air pollution. It was determined that all cities in the Marmara Region study area exceeded the limit PM_2,5 values specified by the European Union (EU) in the years 2016, 2017 and 2018 while only Kocaeli and Tekirdağ were below the limit values in 2019. The limit values specified by the WHO were exceeded in all cities in each year. A total of 46.920 premature deaths attributed to the exceedance of WHO limit values were calculated for the years 2016, 2017, 2018 and 2019 with 11.895, 13.853, 11.748 and 9.429, respectively.

Determining national limit values for the PM_2.5 parameter, which is among the most important factors of air pollution, and monitoring it in a sustainable manner using a sufficient number of well-equipped stations is of great importance. This way, national, regional and urban action plans regarding the impact of air pollution on human health, as indicated by UN SDGs, can be prepared.

1. Introduction

Air pollution is of vital importance for many areas including energy, climate, transportation, trade, agriculture and biological diversity. In addition to posing risks for the environment, economy and food safety, its impact on human health is also critical (Pope et al., 1995; Rafaj et al., 2018). It has been revealed through the epidemiological studies conducted since the 1980s that the relationship between particulate pollution and mortality is consistent. In these studies, it was reported that air pollution at high concentrations has caused the death of hundreds of thousands of individuals in Europe and the United States (Ware et al., 1981; Katsouyanni et al., 1997; Dominici et al., 2003; Samoli et al., 2008; Fattore et al., 2011; Ghanbari Ghozikali et al., 2016). The negative effect of particulate matters, especially those with a grain diameter of less than 2.5 mm (PM_2.5), on human health was determined not only at high concentrations but also at low concentrations, and such matters are linked to respiratory morbidity and mortality (Pope et al., 1995; Sun et al., 2004; Bell, Peng and Dominici, 2006; Downs et al., 2007). According to previous epidemiological studies, PM_2.5 has a stronger impact in terms of premature deaths and potential health effects from prolonged exposure compared to PM₁₀ (Hamanaka and Mutlu, 2018). It was reported that particulate matter of this size bypasses the defence mechanisms of the upper respiratory tract and affects all body systems by first reaching the pulmonary alveoli and then the systemic circulation (Lee, Kim and Lee, 2014).

One of the important risk factors for particularly lung cancer (Vena, 1982) stroke (Hong et al., 2002), type II diabetes (Rao et al., 2015) and cardiovascular diseases (Yamamoto, Phalkey and Malik, 2014) is air pollution. The risk factor of air pollution varies at different levels of socio-economic development. For example, urban air pollution poses a greater risk factor in middle-income countries compared to high-income countries because significant improvements have been achieved in high-income countries in terms of controlling this risk with public health policies (WHO, 2009). Current calculations attribute 7 million premature deaths around the world to air pollution, and approximately 4.2 million of these result from outdoor air pollution (WHO, 2021). The economic cost of premature deaths attributed to outdoor particulate matter pollution and indoor air pollution to the European Nation in 2010 was approximately 1.5 trillion dollars (Prüss-Üstün et al., 2016). In a study conducted by the WHO on deaths attributed to outdoor air pollution for the year 2012, it was reported that approximately 36% of the deaths resulted from stroke while this rate was 36% for ischemic heart diseases, 14% for lung cancer, 8% for chronic obstructive lung diseases and 6% for acute lower respiratory tract infections (WHO, 2016a).

Findings obtained from epidemiological studies show that exposure to air pollution is essentially related to five diseases depending on the age group. These can be listed as: i) Acute low respiratory infections (ALRI) in young children, ii) Cerebrovascular diseases (stroke) in adults, iii) Ischaemic heart diseases (IHD) in adults, iv) Chronic obstructive pulmonary disease (COPD) in adults and v) Lung cancer in adults (WHO, 2016c). Various epidemiological studies determined a relationship between the daily increase in suspended particulate matter concentration and the increase in mortality rate, while also revealing the relationship between mortality and the number of hospitalizations due to respiratory and cardiovascular diseases (Miri et al., 2016).

For Turkey, which is a party to the SDGs, the Marmara Region is the most significant geographical region in terms of both urban and industrial activities. While the six metropolitan municipalities surrounding the Sea of Marmara host some of the most important industrial production points in Turkey, they also represent settlements with the highest urbanization rate and population density. Despite the fact that air pollution resulting from the use of coal for heating purposes in the Marmara Region has decreased over the years, the increase in population-related urban activities and industrial activities keep air pollution among significant health threats.

In the present study, the change in mortality rate attributed to particulate matter air pollution in the Marmara Region over the years (2016–2019) has been examined and the cities were compared. In this comparison, it was decided to use the AirQ + model proposed by the WHO. The results obtained were compared with those gathered using the same model for different cities around the world, and interpreted.

2. Materials And Methods

In the present study on the average level of air pollution and the impact of this level on human health, the years 2016, 2017, 2018 and 2019 are evaluated. 2020 data have been excluded in order to avoid the misleading interferences of the impact of the COVID-19 pandemic on air quality. The flow diagram in Figure 1 shows the steps followed in the study and the path followed in the assessment of the findings.

Study area and air quality monitoring stations

The Marmara Region consists of 11 cities located around the Sea of Marmara in the northwest of Turkey where approximately 31% of the country's population lives. There are 6 metropolitan municipalities in the Marmara Region where urban and industrial activity is relatively denser. By population, these cities can be listed as İstanbul (15.462.452), Bursa (3.101.833), Kocaeli (1.997.258), Balıkesir (1.240.285), Tekirdağ (1.081.065) and Sakarya (1.042.649) (TURKSTAT, 2021). The surface area of the 6 metropolitan municipalities selected as the study area is approximately 45.000 km² and these cities represent the region as they account for approximately 68% of its terrestrial area. Air quality measurements in Turkey are conducted regionally by the Clean Air Centres affiliated with the Ministry of Environment and Urbanization. As of 2019, there are 5 air quality measurement stations affiliated with the Marmara Clean Air Center in Balıkesir, 6 in Bursa, 12 in İstanbul, 12 in Kocaeli, 4 in Sakarya and 5 in Tekirdağ. Figure 2 shows the location of the study area and the air quality measurement stations located in the area.

The data of the stations that made measurements for more than 75% of the days in a calendar year were regarded as suitable within the scope of the present study. The daily average PM_2.5 concentrations of the stations that met these conditions were averaged, and the annual average PM_2.5 concentration was calculated by considering all of the resulting daily average PM_2.5 concentrations of each city. Within the scope of the present study, the PM_2.5 concentrations for the air quality measurement stations where no PM_2.5 measurement was performed were calculated by converting the PM₁₀ concentrations measured using the conversion coefficient (the conversion coefficient specified for Turkey by the WHO is 0,666128). Table 1 shows the number of stations that performed direct PM_2.5 measurements in the cities included in the present study, and the number of stations where the PM_2.5 concentration was calculated by converting PM_10, by years.

Table 1: Number of stations for measured and converted PM_2,5 concentrations

2.2 Age-specific mortality rate

Crude death rate (CDR) is defined as the number of deaths per 100.000 of the population in year (UN, 1991, 2001, 2021; WHO, 2016c). In the present study, the multi-year natural death rates for the cities were calculated to cover the period between 2009 and 2018. Deaths that could not technically be associated with air pollution, such as traffic accidents, murders and suicides, were not included in the total number of deaths (Revised Mortality Underlying Cause Instructions, 2014). Age-specific mortality rates and statistics related to the causes of death were obtained from the Turkish Statistical Institute by city.

Formula 1 is used to calculate the age-specific mortality rate (ASMR) for adult individuals aged over 30 in each city.

For this calculation, the age-specific mortality rate statistics published by the Turkish Statistical Institute and the Address-Based Population Registration System were used. N represents the number of natural deaths over the age of 30 while P represents the population in the same age group. It was acknowledged that all individuals aged over 30 were exposed to the average level of the measured or converted PM_2.5 concentrations.

2.3 Air Q model and calculation of mortality rate attributed air pollution

Version 2.0 o the AirQ+ software, which was developed by the World Health Organization (WHO) Regional Office for Europe and the Bonn Office of the European Centre for Environment & Human Health to calculate the mortality and morbidity of air pollution at a given settlement with reference to the General Principles of Health Risk Assessment of Air Pollution (WHO, 2016b), was used in the present study. Using AirQ+, the number of deaths attributed to long-term PM_2.5 exposure exceeding the ambient air level of 10 mg/m³ as specified in air quality guidelines was calculated.

The Integrated Exposure-Response (IER) function is used in the prediction of the burden of disease attributed to long-term exposure to PM_2.5 in ambient air. (Burnett et al., 2014). Despite the presence of certain limitations, the IER approach provides a beneficial framework to evaluate the consistency of findings and estimate the burden of disease attributed to air pollution (Pope, Cohen and Burnett, 2018). The AirQ+ tool was selected as it uses this function, which enables the effects of air pollution originating from different sources to be evaluated together.

The measurement of the health effects of exposure to air-pollutant materials is based on the concept of Attributable Proportion (AP) and the calculation of Years of Life Loss (YLL). AP, i.e. the fraction of the health outcome, exposure of a specific population within a certain time period can be calculated using Formula 2 (Ezzati et al., 2002):

Within this equation, RR(c) represents the relative risk for the health outcome of exposure in category c while P(c) represents the population rate in category c of exposure.

Relative Risk (RR) is a function used to predict the size of the relationship between exposure and disease. Relative risks originating from air pollution are generally modelled with log-linear function (Formula 3) (Ezzati et al., 2002; Lim et al., 2012).

In the equation, X represents the pollutant concentration in mg/m³. X₀ shows the cut-off or counterfactual concentration. In the present study, the annual average cut-off value of 10mg/m³ (as specified in the WHO air quality guidelines) was used for X₀ in the AirQ+ model. Finally, b represents the change in Relative Risk (RR) for a change of one unit at a concentration of X.

For all-cause mortality, the AirQ+ model takes the meta-analysis of 13 cohort studies conducted by Hoek et al. (2013) as reference and regards the all-cause mortality as 1.062 (95% CI: 1,040–1,083) at a confidence interval (CI) of 95% (Héroux et al., 2015; Hoek et al., 2013).

3. Results And Discussion

3.1 Evaluation of air quality measurements 

In the cities included in the present study, data was obtained from a total of 32 stations in 2016 while this number was 42 in 2019. Table 2 shows the annual PM_2.5 concentrations recorded in the cities and the descriptive statistics parameters Mean±Standard Deviation, median and Interquartile Range – IQR values.

Table 2: Descriptive statistics parameters for PM_2.5 concentrations of cities by year

Whether the PM_2.5 concentrations of the cities differed within themselves in the times of measurement was examined using repeated measures ANOVA. According to this, a significant difference (p<0.05) was detected between four temporal measurements of PM_2.5 values in all cities (2016, 2017, 2018, 2019). It can be said that air quality was relatively improved in 2019 compared to the other years.

It was recorded that 68% of a total of 1.327 stations in 32 countries reporting to the European Environment Agency exceeded the annual PM_2.5 concentration limit of 10 mg/m³ (as specified in WHO guidelines) in 2016. It is observed that with the exception of Iceland, Estonia, Finland, Ireland, Sweden, Switzerland, Norway, United Kingdom, Austria, Portugal and Germany, a significant portion of the countries with average PM_2.5 concentrations below the limit value of 25 mg/m³ specified in the EU legislation exceeded the limit determined by the WHO (EEA, 2018). The annual average PM_2.5 limit value is included in national legislation in certain countries. Based on this, this value was specified as 25 mg/m³ in the EU, 8 mg/m³ in Australia, 10 mg/m³ in Canada and Japan, 12 mg/m³ in the United States and Mexico, 15 mg/m³ in Ecuador, El Salvador and the Dominican Republic, 20 mg/m³ in South Africa, 15 mg/m³ in rural parts of China, 35 mg/m³ in urban parts of China and 40 mg/m³ in India. (THHP, 2019). It is observed that all cities included in the present study exceeded the limit value of 10mg/m³ specified by the WHO in each year (2016-2019). Additionally, it is observed that Balıkesir, İstanbul, Kocaeli and Tekirdağ had annual average PM_2.5 concentrations at EU limits in 2019, but the annual average PM_2.5 concentrations of Bursa and Sakarya were above the EU limits. It is thought that the air pollution occurring in Bursa, particularly in winter months, essentially originates from industrialization, and the significant pollution factor is coal, which is used for domestic heating and as an industrial fuel. It is known that the ring road passing through the north of the city creates traffic-related pollution based on the prevailing wind direction. The compaction of the particulate matter carried from the north through this type of wind also affects the air quality in Bursa negatively (MoEU, 2020a). Agricultural areas are located in the north of the city and the dust originating from agricultural activities is regarded as among the important factors of inner-city air pollution. The use of coal with low calorific value is the leading air pollution factor in Sakarya. In relation to this, the low air circulation in the city centre, particularly during winter months, causes pollutants to get concentrated in the low parts of the city. Additionally, the number of facilities that were granted environmental permits in terms of emission control was 94 in the year 2016 while this number has increased to 220 in the year 2020, which is an important fact in terms of industrial air pollution (MoEU, 2020b). 

The average PM_2.5 concentrations of the cities, recorded each day of the year between 2016 and 2019, were examined and Figure 3 shows the daily concentration accumulations. In all cities, accumulations were observed in the second quarter in the daily PM_2.5 concentrations recorded. 

Since the cities included in the present study were subjected to multiple comparisons within themselves (4 years), whether there was a statistically significant difference in pairwise comparisons between the years (two years) was examined using the Bonferroni Correction. According to this, significant differences were determined in all years except the dual years of 2016-2018 in Balıkesir, 2016-2017 and 2016-2018 in Bursa, 2016-2018 in İstanbul, 2016-2018 in Kocaeli, 2016-2017, 2016-2018 and 2017-2018 in Sakarya, and 2018-2019 in Tekirdağ. Table 3 shows the Bonferroni Correction results regarding the pairwise comparisons of the cities' annual PM_2.5 concentrations.

Table 3: Bonferroni Correction results of dual comparison of years

In the light of these determinations, it is acknowledged that the decrease in the average concentration in each city was recorded as an improvement in air quality since the difference in the dual comparisons of 2016-2019 and 2017-2019 was significant. Additionally, the decrease in the number of days in which the average PM_2.5 concentration limit value of 25mg/m³ specified by the WHO was exceeded within a calendar year is significant with the p<0.05 value in the 2016-2019 dual comparison according to the Bonferroni Correction. Figure 4 shows the change in the number of days in which the WHO limit values were exceeded in the cities by year. 

The average PM_2.5 concentration values of the cities included in the present study during periods where sunny and cloudy days were common in the years 2016, 2017, 2018 and 2019 were examined using the Unpaired T-Test. September, October, November, December, January and February constituted the cloudy period while March, April, May, June, July and August were regarded as constituting the sunny period. According to this, it was observed that the annual seasonal differences between the sunny/cloudy periods in all cities were statistically significant (p<0.05) with the exception of the sunny and cloudy periods in Tekirdağ in 2018 (t: 1.846, p:0.066), and that the air quality during the sunny periods was better compared to the cloudy period. Based on this, Figure 5 shows the seasonal differences in the sunny/cloudy periods of the cities by year.

When the change in the average PM_2.5 concentrations of the cities over 48 months is examined, the result obtained is in parallel with the Unpaired T-Test performed for seasonal change. Figure 6 shows the change in the monthly average PM_2.5 values of the cities recorded over 48 months and the increase in the PM_2.5 concentration during cloudy months.

When the monthly average PM_2.5 concentrations of the cities included in the present study are considered, the highest monthly average in 2016 was measured during the cloudy months while the lowest average was measured during the sunny months with the exception of Balıkesir and Tekirdağ. In 2017, the highest monthly average was measured during the cloudy months with the exception of Tekirdağ while the lowest monthly average was measured during the sunny months in all cities. In 2018, the highest monthly average was measured during the cloudy months with the exception of İstanbul and Kocaeli while the lowest monthly average was measured, once again, during the cloudy months with the exception of Bursa and Sakarya. In 2019, the highest monthly average was measured during the cloudy months while the lowest monthly average was measured during the sunny months. Table 4 shows the months in which the highest and lowest PM_2.5 values were recorded in the cities between 2016 and 2019, along with mean±sd values.

Table 4: Highest and lowest monthly average PM_2.5 concentrations

Meteorological factors are important in terms of understanding the change in air pollutant concentrations and the change in their distribution. In addition to other air pollutants, the tendencies of PM are also affected by topography and meteorological variables (Rodríguez, Dupont-Courtade and Oueslati, 2016; Mohtar et al., 2018; Liang and Gong, 2020). In the study conducted by Bodor et al. (2020) in the Transylvania region of Romania, the highest monthly average PM concentration was recorded in February at four times the annual average. Additionally, a negative correlation was observed between PM and temperature (Bodor et al., 2020). 

It was revealed by Cichowichz et al. (2017) that air pollutant concentrations in Poland are high during the winter months and heating seasons and low during the summer months. Mohtar et al. (2018) revealed that PM concentrations in Malaysia were high during these seasons in relation to burning biomass and coal with low calorific value. In the present study, statistically significant differences were determined in the average PM_2.5 concentrations during the seasons categories by the number of sunny and cloudy days, and it was revealed that PM-related air pollution was highest in all cities during November, December, January and February.

3.2 Natural mortality rates and air pollution attributable deaths

In the cities included in the present study, the natural mortality rates in adults aged over 30 were calculated in units of 100 thousand, taking into account all mortality cases except those in relation to external injuries and poisoning. Figure 7 shows the mortality rates in the cities.

The natural mortality rate in adults (30+) during the years 2013-2014 was recorded as 876 in Paris, 871 in Santiago, the capital of Chile, during 2013-2016, and 934 in the Lombardia region of Italy during 2013-2016 (WHO, 2018). In the present study, the city with the highest mortality rate was recorded as Balıkesir while Istanbul had the lowest mortality rate.

As a fraction of the health outcome, the exposure of the population aged over 30 in the cities was calculated in the form of attributable proportion and shown in Figure 8.

In the attributable proportion, upper and lower rates were calculated by considering a confidence interval of 95% (95% CI). Based on this, among the four years considered, a maximum attributable proportion of 26.32% is calculated in Bursa for the year 2017. In the study conducted by Toscano et al. (2021) covering the period between 2003 and 2019, it is reported that the sharp decrease in the NO₂ concentration in Italy’s Campania region and the relative decrease in the PM₁₀ and PM_2.5 concentrations also caused a decrease in the attributable proportion parameter.

The AirQ model was operated with the Integrated Exposure-Response Function reference using the age-based (30+) natural mortality rates and annual average PM_2.5 concentrations in the cities. Based on this, Figure 9 shows the premature deaths caused by long-term exposure to PM_2.5.

In relation to the statistically significant decrease (p<0.05) in PM_2.5 concentrations between 2016 and 2019, it is thought that the decreasing trend in the number of premature deaths in all cities is statistically significant. It is thought that premature deaths may be prevented if air quality is improved to a point where the annual average PM_2.5 concentration in cities is reduced to the annual average PM_2.5 concentration threshold specified in the WHO guidelines.

Pope & Dockery (2006) revealed that an increase of 10 μg/m³ in the 24-hour average PM_2.5 concentration increases the relative risk of daily cardiovascular mortality by 0.6-1.1%. Additionally, Garcia et al. (2016) revealed a significant relationship between PM_2.5 exposure and deaths caused by cardiovascular diseases, ischemic heart diseases and cardiopulmonary diseases in California. Approximately 39% of all deaths in Turkey in the years 2016, 2017 and 2018 resulted from vascular diseases. Approximately 65% of all vascular diseases in Turkey result from ischemic heart disease and other heart diseases (TURKSTAT, 2019). In the examination performed within the scope of the present study, it was calculated that 46% of the deaths in Balıkesir resulted from vascular diseases while this rate was 44% in Bursa and 45% in Sakarya. This rate was calculated as 34% for İstanbul, 34% for Kocaeli and 40% for Tekirdağ. It is thought that the average PM_2.5 concentrations measured, the mortality rates in the cities and national cause of death statistics are in a significant relationship.

Figure 10 shows the mortality rate (per 100 thousand) of preventable premature deaths in similar to attributable proportion. This indicator is of importance as it allows for a depopulated assessment on the health effects of air pollution. Although there was a decrease in the mortality rates in Bursa, Balıkesir and Sakarya between 2016 and 2019, it is observed that mortality rate in these cities is clearly higher compared to other cities. 

In addition to the significant decrease (p<0.001) in the mortality rate in Tekirdağ from 2016 to 2019, it is also acknowledged that the change between 2017 and 2018 is statistically significant (p<0.001) and very apparent. Additionally, İstanbul and Kocaeli had a very low mortality rate compared to the other four cities in 2016, with a steady decrease (p<0.05) up to 2019. The dual comparison results regarding the PM_2.5 value in the cities were found to be statistically significant (p<0.05). Based on this significant finding, the PM_2.5 values of the cities differ over the years. Only İstanbul and Kocaeli had similar PM_2.5 values over the years (p=0.493>0.05).

In a study conducted in 11 of the largest cities in Poland, which is the country with the highest annual average PM_2.5 concentration among European Union countries, it was reported that the annual average PM_2.5 concentration ranged between 14.3 and 52.5 μg/m³. The number of deaths attributed to long-term exposure to PM_2.5 concentration was estimated to be 9.6 to 22.8 of every 100.000 lung cancer deaths and 48.6 to 136.6 in cases resulting from heart and lung diseases (Badyda, Grellier and Dabrowiecki, 2017). It is estimated that in 25 cities in Iran where PM_2.5 pollution is common, the annual average PM_2.5 concentrations are 1.5-6.1 times higher than the WHO threshold value and, consequently, 13.321 of all annual deaths recorded in these cities are classified as premature death attributed to air pollution (Hadei et al., 2020). In the present study, a decrease was observed in the overall natural mortality rates between the years 2016 and 2019. In 2019, the cases in Balıkesir, Bursa and Sakarya were calculated as above 100 (per 100 thousand), higher than İstanbul, Kocaeli and Tekirdağ.

The impact of sharp increases and decreases in air pollutant concentrations on morbidity and mortality has been subject to many studies. In the study conducted by Turner et al. (2011), it is stated that an increase of 10 μg/m³ in the PM_2.5 level has caused an increase of 15-27% in deaths from cancer. Mahesvaran et al. (2010) examined the impact of air pollution on post-stroke survival rates with 3.320 patients and 1.856 death cases in London. Based on this, it was revealed that an increase of 10 μg/m³ in the concentration of particulate matter with a diameter of less than 10 μm increased the probability of poststroke mortality by 52%. In the study conducted by Brook et al. (2010) on the relationship between cardiovascular diseases and air pollution, it was reported that a decrease of 10 μg/m³ in PM_2.5 concentrations increased life expectancy by 0.61 ± 0.20 years. In the study conducted by Pala et al. (2021) in all of Turkey, it was stated that 44.617 cases of premature death in 2018 could have been prevented had the air quality been at the annual average PM_2.5 concentration level specified by the WHO. In this context, in the present study, 11.895 (95% CI, 7.951-15.409), 13.853 (95% CI, 9.281-17.906), 11.748 (95% CI, 7.838-15.243) and 9.429 (95% CI, 6.250-12.303) premature deaths were attributed to excessive air pollution values in 2016, 2017, 2018 and 2019, respectively. It is thought that such premature deaths can be prevented if a level of air quality that is in accordance with the WHO limits can be established.

4. Conclusion

With the modern industrialization and motorization around the world, the emission of anthropogenic aerosols is sharply increasing, which brings along not only regional but also global effects by means of long-range transportation. Aerosols are globally regarded for their potential effects on climate change and human health (Sun et al., 2004). For this reason, the goal of reducing the death and diseases caused by air pollution by 2030 is pointed out under the 3rd Goal titled Healthy and Quality Life in the United Nations Sustainable Development Goals published in 2015. Additionally, under the 11th Goal titled Sustainable Cities, the annual average of thin particulate matter concentration in cities was specified as a target indicator (WHO, 2016a). For these purposes, in the present study covering the period between 2016 and 2019 for 6 cities in the Marmara Region, premature deaths attributed to air pollution were calculated based on the PM_2.5 concentrations measured, and the relationship between air quality, human health and life expectancy was focused on. The most striking aspect of the present study is that the impact of COVID-19 on the premature deaths attributed to air pollution in the Marmara Region and findings and conclusions related to urban comparisons can be the subject of a separate study.

PM_2.5 is regarded as one of the most essential parameters in understanding air quality and its impact on human health (EEA, 2018; THHP, 2019; WHO, 2016a, 2021). Hamanaka & Mutlu (2018) revealed that PM_2.5 had a stronger effect on human health compared to PM₁₀. Additionally, Lee et al. (2014) emphasized that PM_2.5 affected the whole body system. For this reason, directly measuring the PM_2.5 parameter rather than calculating with any conversion coefficient is of importance in terms of monitoring this parameter in a sustainable manner. Additionally, the relatively small number of stations that directly measure the PM_2.5 parameter is thought to be related to the fact that no national-scale limit values were specified regarding the PM_2.5 parameter and that it was not included in air quality regulations. In this respect, specifying limit values for PM_2.5 in the Regulation on Air Quality Assessment and Management within the framework of WHO documents and EU regulations is thought to be of importance for the effective monitoring of its effects on human health.

It is known that almost all parameters related to air pollution originate from incineration processes, industrial activities and traffic. It is also known that atmospheric convection and meteorological conditions play a role in the distribution and effects of air pollutants. For this reason, it is necessary to limit the use of coal with low calorific value in order to reduce air pollution caused by heating. It may be effective for gas distribution firms and local governments to cooperate in providing incentives and support mechanisms for the transition from coal to natural gas. It would be beneficial to conduct studies towards adopting sustainable production and consumption principles with the aim of reducing industrial air pollution. These studies can focus on the main subjects of technological developments, sustainability of administrative sanctions and the effective use of economic tools. In order to reduce traffic-based emissions in urban areas, it would be beneficial to promote walking and cycling short distances through infrastructure improvements. It is also thought that it would be effective to introduce administrative regulations promoting the use of electric or hybrid vehicles rather than those using fossil fuels. Additionally, providing digital communication tools that can recommend alternative routes to drivers and notify pedestrians on what time to go out in order to reduce exposure to traffic-based emissions as a public service is thought to be beneficial from many aspects. Finally, it is of importance to plan residential areas in consideration of air pollution dynamics, meteorological conditions and topographical factors.   

Declarations

Acknowledgement

We would like to thank Mr. Emre Yaşar who supported us in checking the necessary statistical analyses within the study, the Marmara Clean Air Centre who provided the location information of the Air Quality Measurement Stations in the Marmara Region. Also, the authors special thank to Prof. Dr. Burcu Onat and Prof. Dr. Ülkü Alver Şahin for their valuable assessments and comments.

Authors’ contributions

ACK: literature search, study design, data acquisition, data analysis, statistical analysis, AirQ model results, manuscript preparation, corresponding author

NS: concept, design, data interpretation, definition of intellectual content, manuscript editing, 'guarantor'

Funding

There is not any funding body in the design of the study and collection, analysis, and interpretation of data.

Data availability

The datasets generated and/or analysed during the current study are available in the;

Ministry of Environment and Urbanization, National Air Quality Monitoring Network repository, [https://www.havaizleme.gov.tr/]

Turkish Statistical Institute, Death Statistics repository, [https://biruni.tuik.gov.tr/medas/?kn=114&locale=en]

Compliance with ethical standard

Ethics approval and consent to participate Not applicable

Consent for publication Not applicable

Competing interests The authors declare that they have no competing interests

References

Badyda, A. J., Grellier, J., & Dabrowiecki, P. (2017). Ambient PM2.5 Exposure and Mortality Due to Lung Cancer and Cardiopulmonary Diseases in Polish Cities. Advances in Experimental Medicine and Biology, 944, 9–17. https://doi.org/10.1007/5584_2016_55

Bell, M. L., Peng, R. D., & Dominici, F. (2006). The exposure-response curve for ozone and risk of mortality and the adequacy of current ozone regulations. Environmental Health Perspectives, 114(4), 532–536. https://doi.org/10.1289/ehp.8816

Bodor, Z., Bodor, K., Keresztesi, Á., & Szép, R. (2020). Major air pollutants seasonal variation analysis and long-range transport of PM10 in an urban environment with specific climate condition in Transylvania (Romania). Environmental Science and Pollution Research, 27(30), 38181–38199. https://doi.org/10.1007/s11356-020-09838-2

Brook, R. D., Rajagopalan, S., C. Arden Pope, I., Brook, J. R., Bhatnagar, A., Diez-Roux, A. V., Holguin, F., Hong, Y., Luepker, R. V., Mittleman, M. A., Peters, A., Siscovick, D., Sidney C. Smith, J., Whitsel, L., & Kaufman, J. D. (2010). Particulate Matter Air Pollution and Cardiovascular Disease. Circulation, 121(21), 2331–2378. https://doi.org/10.1161/CIR.0B013E3181DBECE1

Burnett, R. T., Arden Pope, C., Ezzati, M., Olives, C., Lim, S. S., Mehta, S., Shin, H. H., Singh, G., Hubbell, B., Brauer, M., Ross Anderson, H., Smith, K. R., Balmes, J. R., Bruce, N. G., Kan, H., Laden, F., Prüss-Ustün, A., Turner, M. C., Gapstur, S. M., … Cohen, A. (2014). An integrated risk function for estimating the global burden of disease attributable to ambient fine particulate matter exposure. Environmental Health Perspectives, 122(4), 397–403. https://doi.org/10.1289/ehp.1307049

Cichowicz, R., Wielgosiński, G., & Fetter, W. (2017). Dispersion of atmospheric air pollution in summer and winter season. Environmental Monitoring and Assessment, 189(12). https://doi.org/10.1007/s10661-017-6319-2

Dominici, F., McDermott, A., Zeger, S. L., & Samet, J. M. (2003). Airborne particulate matter and mortality: Timescale effects in four US cities. American Journal of Epidemiology, 157(12), 1055–1065. https://doi.org/10.1093/aje/kwg087

Downs, S. H., Schindler, C., Liu, L.-J. S., Keidel, D., Bayer-Oglesby, L., Brutsche, M. H., Gerbase, M. W., Keller, R., Künzli, N., Leuenberger, P., Probst-Hensch, N. M., Tschopp, J.-M., Zellweger, J.-P., Rochat, T., Schwartz, J., & Ackermann-Liebrich, U. (2007).  Reduced Exposure to PM 10 and Attenuated Age-Related Decline in Lung Function . New England Journal of Medicine, 357(23), 2338–2347. https://doi.org/10.1056/nejmoa073625

EEA. (2018). Air quality in Europe 2018 (Issue 12/2018). papers2://publication/uuid/1D25F41B-C673-4FDA-AB71-CC5A2AD97FDD

Ezzati, M., Lopez, A. D., Rodgers, A., Vander Hoorn, S., & Murray, C. J. L. (2002). Selected major risk factors and global and regional burden of disease. Lancet, 360(9343), 1347–1360. https://doi.org/10.1016/S0140-6736(02)11403-6

Fattore, E., Paiano, V., Borgini, A., Tittarelli, A., Bertoldi, M., Crosignani, P., & Fanelli, R. (2011). Human health risk in relation to air quality in two municipalities in an industrialized area of Northern Italy. Environmental Research, 111(8), 1321–1327. https://doi.org/10.1016/j.envres.2011.06.012

Garcia, C. A., Yap, P.-S., Park, H.-Y., & Weller, B. L. (2016). Association of long-term PM2.5 exposure with mortality using different air pollution exposure models: impacts in rural and urban California. International Journal of Environmental Health Research, 26(2), 145. https://doi.org/10.1080/09603123.2015.1061113

Ghanbari Ghozikali, M., Heibati, B., Naddafi, K., Kloog, I., Oliveri Conti, G., Polosa, R., & Ferrante, M. (2016). Evaluation of Chronic Obstructive Pulmonary Disease (COPD) attributed to atmospheric O3, NO2, and SO2 using Air Q Model (2011-2012 year). Environmental Research, 144(2016), 99–105. https://doi.org/10.1016/j.envres.2015.10.030

Hadei, M., Shahsavani, A., Krzyzanowski, M., Querol, X., Stafoggia, M., Nazari, S. S. H., Jafari, A. J., Yarahmadi, M., Kermani, M., & Khosravi, A. (2020). Burden of mortality attributed to PM2.5 exposure in cities of Iran; contribution of short-term pollution peaks. Atmospheric Environment, 224(February), 117365. https://doi.org/10.1016/j.atmosenv.2020.117365

Hamanaka, R. B., & Mutlu, G. M. (2018). Particulate Matter Air Pollution: Effects on the Cardiovascular System. Frontiers in Endocrinology, 9. https://doi.org/10.3389/FENDO.2018.00680

Héroux, M. E., Anderson, H. R., Atkinson, R., Brunekreef, B., Cohen, A., Forastiere, F., Hurley, F., Katsouyanni, K., Krewski, D., Krzyzanowski, M., Künzli, N., Mills, I., Querol, X., Ostro, B., & Walton, H. (2015). Quantifying the health impacts of ambient air pollutants: recommendations of a WHO/Europe project. International Journal of Public Health, 60(5), 619–627. https://doi.org/10.1007/s00038-015-0690-y

Hoek, G., Krishnan, R. M., Beelen, R., Peters, A., Ostro, B., Brunekreef, B., & Kaufman, J. D. (2013). Long-term air pollution exposure and cardio-respiratory mortality: A review. Environmental Health: A Global Access Science Source, 12(1). https://doi.org/10.1186/1476-069X-12-43

Hong, Y. C., Lee, J. T., Kim, H., & Kwon, H. J. (2002). Air pollution: A new risk factor in ischemic stroke mortality. Stroke, 33(9), 2165–2169. https://doi.org/10.1161/01.STR.0000026865.52610.5B

Revised mortality underlying cause instructions, WHO International Statistical Classification of Diseases (2014). https://icd.who.int/icd10updateplatform/PropD.aspx?prop=2109

Katsouyanni, K., Touloumi, G., Spix, C., Schwarte, J., Balducci, F., Medina, S., Rossi, G., Wojtyniak, B., Sunyer, J., Bacharova, L., Schouten, J. P., Ponka, A., & Anderson, H. R. (1997). Short term effects of ambient sulphur dioxide and particulate matter on mortality in 12 European cities: Results from time series data from the APHEA project. British Medical Journal, 314(7095), 1658–1663. https://doi.org/10.1136/BMJ.314.7095.1658

Lee, B., Kim, B., & Lee, K. (2014). Air pollution exposure and cardiovascular disease. Toxicological Research, 30(2), 71–75. https://doi.org/10.5487/TR.2014.30.2.071

Liang, L., & Gong, P. (2020). Urban and air pollution: a multi-city study of long-term effects of urban landscape patterns on air quality trends. Scientific Reports, 10(1), 1–13. https://doi.org/10.1038/s41598-020-74524-9

Lim, S. S., Vos, T., Flaxman, A. D., Danaei, G., Shibuya, K., Adair-Rohani, H., Amann, M., Anderson, H. R., Andrews, K. G., Aryee, M., Atkinson, C., Bacchus, L. J., Bahalim, A. N., Balakrishnan, K., Balmes, J., Barker-Collo, S., Baxter, A., Bell, M. L., Blore, J. D., … Ezzati, M. (2012). A comparative risk assessment of burden of disease and injury attributable to 67 risk factors and risk factor clusters in 21 regions, 1990-2010: A systematic analysis for the Global Burden of Disease Study 2010. The Lancet, 380(9859), 2224–2260. https://doi.org/10.1016/S0140-6736(12)61766-8

Maheswaran, R., Pearson, T., Smeeton, N. C., Beevers, S. D., Campbell, M. J., & Wolfe, C. D. (2010). Impact of outdoor air pollution on survival after stroke: population-based cohort study. Stroke, 41(5), 869–877. https://doi.org/10.1161/STROKEAHA.109.567743

Miri, M., Derakhshan, Z., Allahabadi, A., Ahmadi, E., Oliveri Conti, G., Ferrante, M., & Aval, H. E. (2016). Mortality and morbidity due to exposure to outdoor air pollution in Mashhad metropolis, Iran. The AirQ model approach. Environmental Research, 151, 451–457. https://doi.org/10.1016/j.envres.2016.07.039

MoEU. (2020a). Bursa Clean Air Action Plan 2020-2024.

MoEU. (2020b). Sakarya Clean Air Action Plan 2020-2024.

Mohtar, A. A. A., Latif, M. T., Baharudin, N. H., Ahamad, F., Chung, J. X., Othman, M., & Juneng, L. (2018). Variation of major air pollutants in different seasonal conditions in an urban environment in Malaysia. Geoscience Letters, 5(1). https://doi.org/10.1186/s40562-018-0122-y

Pala, K., Aykac, N., & Yasin, Y. (2021). Premature deaths attributable to long-term exposure to PM2.5 in Turkey. Environmental Science and Pollution Research, 2017. https://doi.org/10.1007/s11356-021-13923-5

Pope 3rd, C. A., & Dockery, D. W. (2006). Health effects of fine particulate air pollution: lines that connect. Journal of the Air & Waste Management Association (1995), 56(6), 709–742. https://doi.org/10.1080/10473289.2006.10464485

Pope, C. A., Cohen, A. J., & Burnett, R. T. (2018). Cardiovascular disease and fine particulate matter lessons and limitations of an integrated exposure-response approach. Circulation Research, 122(12), 1645–1647. https://doi.org/10.1161/CIRCRESAHA.118.312956

Pope, C. A., Thun, M. J., Namboodiri, M. M., Dockery, D. W., Evans, J. S., Speizer, F. E., & Heath, C. W. (1995). Particulate Air Pollution as a Predictor of Mortality in a Prospective Study of U.S. Adults. American Journal of Respiratory and Critical Care Medicine, 151(3_pt_1), 669–674. https://doi.org/10.1164/ajrccm/151.3_Pt_1.669

Prüss-Ustün, A., Wolf, J., Corvalan, R., Bos, R., & Neira, M. (2016). Preventing disease through healthy environments: A global assessment of the environmental burden of disease. In WHO. www.who.int/about/licensing/copyright_form/en/index.html

Rafaj, P., Kiesewetter, G., Gül, T., Schöpp, W., Cofala, J., Klimont, Z., Purohit, P., Heyes, C., Amann, M., Borken-Kleefeld, J., & Cozzi, L. (2018). Outlook for clean air in the context of sustainable development goals. Global Environmental Change, 53(April), 1–11. https://doi.org/10.1016/j.gloenvcha.2018.08.008

Rao, X., Patel, P., Puett, R., & Rajagopalan, S. (2015). Air Pollution as a Risk Factor for Type 2 Diabetes. Toxicological Sciences, 143(2), 231–241. https://doi.org/10.1093/TOXSCI/KFU250

Rodríguez, M. C., Dupont-Courtade, L., & Oueslati, W. (2016). Air pollution and urban structure linkages: Evidence from European cities. Renewable and Sustainable Energy Reviews, 53, 1–9. https://doi.org/10.1016/j.rser.2015.07.190

Samoli, E., Peng, R., Ramsay, T., Pipikou, M., Touloumi, G., Dominici, F., Burnett, R., Cohen, A., Krewski, D., Samet, J., & Katsouyanni, K. (2008). Acute effects of ambient particulate matter on mortality in Europe and North America: Results from the APHENA study. Environmental Health Perspectives, 116(11), 1480–1486. https://doi.org/10.1289/EHP.11345

Sun, Y., Zhuang, G., Wang, Y., Han, L., Guo, J., Dan, M., Zhang, W., Wang, Z., & Hao, Z. (2004). The air-borne particulate pollution in Beijing - Concentration, composition, distribution and sources. Atmospheric Environment, 38(35), 5991–6004. https://doi.org/10.1016/j.atmosenv.2004.07.009

THHP. (2019). Policy Note on PM2.5. https://www.temizhavahakki.com/pm25-politika-notu/

Toscano, D., Murena, F., Costabile, F., Hussein, T., & Massimi, L. (2021). The Historical Trend of Air Pollution and Its Impact on Human Health in Campania Region (Italy). https://doi.org/10.3390/atmos12050553

TURKSTAT. (2019). Distribution of causes of death by gender, 2009-2018. Turkish Statistical Institute. https://data.tuik.gov.tr/Bulten/Index?p=Olum-ve-Olum-Nedeni-Istatistikleri-2019-33710

TURKSTAT. (2021). Address Based Population Registration System Results 2020. Turkish Statistical Institute. https://biruni.tuik.gov.tr/medas/?kn=95&locale=en

Turner, M. C., Krewski, D., Pope 3rd, C. A., Chen, Y., Gapstur, S. M., & Thun, M. J. (2011). Long-term ambient fine particulate matter air pollution and lung cancer in a large cohort of never-smokers. American Journal of Respiratory and Critical Care Medicine, 184(12), 1374–1381. https://doi.org/10.1164/RCCM.201106-1011OC

UN. (1991). Handbook of Vital Satistics Systems and Methods (F/35). United Nations Publication. https://unstats.un.org/unsd/publication/SeriesF/SeriesF_35v1E.pdf

UN. (2001). Principles and Recommendations for a Vital Statistics System (ST/ESA/STA). UN Department of Economic and Social Affairs. https://unstats.un.org/unsd/demographic-social/Standards-and-Methods/files/Principles_and_Recommendations/CRVS/SeriesM_19rev2-E.pdf

UN. (2021). UN Data Glossary. United Nations Statistics Division. http://data.un.org/Glossary.aspx?d=DemogG

Vena, J. E. (1982). Air Pollution as a Risk Factor in Lung Cancer. American Journal of Epidemiology, 116(1), 42–56. https://doi.org/10.1093/OXFORDJOURNALS.AJE.A113401

Ware, J. H., Thibodeau, L. A., Speizer, F. E., Colome, S., & Ferris Jr., B. G. (1981). Assessment of the health effects of atmospheric sulfur oxides and particulate matter: Evidence from observational studies. Environmental Health Perspectives, Vol. 41(October), 255–276. https://doi.org/10.1289/ehp.8141255

WHO. (2009). Global Health Risks: Mortality and Burden of Disease Attributable to Selected Major Risks. http://www.who.int/healthinfo/global_burden_disease/GlobalHealthRisks_report_full.pdf

WHO. (2016a). Ambient air pollution: A global assessment of exposure and burden of disease. http://marefateadyan.nashriyat.ir/node/150

WHO. (2016b). Health Risk Assessment of air pollution. In WHO Regional Office for Europe. World Health Organization. https://www.euro.who.int/en/health-topics/environment-and-health/air-quality/activities/airq-software-tool-for-health-risk-assessment-of-air-pollution%0Ahttp://www.euro.who.int/en/health-topics/environment-and-health/air-quality/activities/airq-software-t

WHO. (2016c). Institutional information of SDG 3.9.1. Mortality Rate Attributed to Household and Ambient Air Pollution (Issue 1, p. 1). World Health Organization. https://www.who.int/data/gho/data/indicators/indicator-details/GHO/ambient-and-household-air-pollution-attributable-death-rate-(per-100-000-population-age-standardized)

WHO. (2018). AirQ+ software tool for health risk assessment of air pollution - Information documents on Mortality (2.0; pp. 1–6). World Health Organization. https://www.euro.who.int/en/health-topics/environment-and-health/air-quality/activities/airq-software-tool-for-health-risk-assessment-of-air-pollution

WHO. (2021). Air Pollution. World Health Organization. https://www.who.int/health-topics/air-pollution#tab=tab_1

Yamamoto, S. S., Phalkey, R., & Malik, A. A. (2014). A systematic review of air pollution as a risk factor for cardiovascular disease in South Asia: Limited evidence from India and Pakistan. International Journal of Hygiene and Environmental Health, 217(2–3), 133–144. https://doi.org/10.1016/J.IJHEH.2013.08.003