key: cord-0982453-j0c9krav authors: Orak, Nur H.; Ozdemir, Ozancan title: The impacts of COVID-19 lockdown on PM(10) and SO(2) concentrations and association with human mobility across Turkey date: 2021-03-19 journal: Environ Res DOI: 10.1016/j.envres.2021.111018 sha: c17fa36c0502abd3e68ec898d982b9f6009352e7 doc_id: 982453 cord_uid: j0c9krav The SARS-CoV-2 virus pandemic (COVID-19) has caused 2.25 million deaths worldwide by February 3, 2021 (JHU, 2021) and still causing severe health and economic disruptions with increasing rates. This study investigates the impact of lockdown measures on ambient air pollution and its association with human mobility in 81 cities of Turkey. We conducted a countrywide analysis using PM(10) and SO(2) measurement data by the Turkish Ministry of Environment and Urbanization and mobility data derived from cellular device movement by Google. We observed the most significant change in April 2020. PM(10) and SO(2) concentrations were lower in 67% and 59% of the cities, respectively in April 2020 compared to the previous five years (2015–2019). The correlation results show that Restaurant/Café, Transit, and Workplaces mobility is significantly correlated with PM(10) and SO(2) concentration levels in Turkey. This study is the first step of a long-term investigation to understand the air quality impacts on population susceptibility to COVID-19. COVID-19 is a global infectious disease (Zhu et al., 2020) caused by severe acute respiratory syndrome coronavirus 2 (WHO, 2020a), which was first reported in Wuhan, China in late 2019 (Lu et al., 2020) . On March 11, 2020, the World Health Organization named COVID-19 as a pandemic with more than 118,000 confirmed cases in 114 countries, and 4291 deaths (WHO 2020). COVID-19 has caused 2.25 million deaths worldwide by February 3, 2021 (JHU, 2021 . Different countries implemented different lockdown measures to protect public health by minimizing the movement of citizens. As a result, the pandemic caused an estimated global consumption loss of $3.8 trillion, job losses equivalent to 147 million full-time positions, and a loss of $2.1 trillion in wages and salaries (Lenzen et al., 2020) . The Ministry of Health of Turkey confirmed the first case of COVID-19 on March 11, 2020 in Istanbul, Turkey (TRT 2020). The Turkish government announced the first restrictions on March 12, 2020, to limit the mobility of citizens, schools were closed starting from March 16, 2020, sports events were canceled and state officers were limited to travel internationally until further notice (BBC 2020) . According to the Stringency Index of Oxford COVID-19 Government Response Tracker (Hale et al., 2020) , Turkey started with a stringency level of 23.2 on March 11th (Fig. S1 ). This index is driven by policy indicators (Table S1 ). The value of the index on a specific day is the average of nine sub-indices. Fig. S1 shows the gradual increase of the stringency index over time until March 28 with 75.9 and decreased to 63.9 on June 2, 2020, in Turkey. The ease on lockdown continued during summer and the index dropped to 42.6 by July 25, 2020. However, restrictions came back by September due to a high number of new cases, the index was as high as 73.6 with the new restriction on September 12, 2020. Several studies investigated the impact of air quality and meteorological variables on transmission dynamics of COVID-19 (Bashir et al., 2020; Coccia, 2020a Coccia, , 2021a Coccia, , 2021b Srivastava, 2021) . Although the environmental impacts of COVID-19 lockdowns are still under investigation by several research groups all around the world, several studies have investigated the change in atmospheric pollution during the pandemic (Archer et al., 2020; Baldasano, 2020; Berman and Ebisu, 2020; Chen et al., 2020; Kerimray et al., 2020; Krecl et al., 2020; Muhammad et al., 2020; Pei et al., 2020) . A major decline in Nitrogen dioxide (NO₂) emissions over early 2020 were observed by satellite images from the Centre for Research on Energy and Clean Air (CREA), NASA, and the European Space Agency (ESA) (NASA, 2020) . Another study showed a decline in PM 2.5 (11.3%) and NO 2 (25.5%) during the pandemic in the US (Berman and Ebisu, 2020) . Similarly, 30% reduction in NO 2 during the month of March in the urban northeastern US was observed by satellite data (Blumberg, 2020). Bao and Zhang (2020) showed that the air quality index decreased by 7.80% and SO 2 , PM 2.5 , PM 10 , NO 2 , and CO decreased by 6.76%, 5.93%, 13.66%, 24.67%, and 4.58%, respectively as an average of 44 cities in China due to lockdown measures. Another study estimated a decrease in the carbon footprint burden by 20% compared with 2015-2019 in Italy (Rugani and Caro, 2020) . In Brazil, researchers observed significant reductions (30.3-48.5) in CO levels and limited reductions for PM 10 levels during the first week of the lockdown (Dantas et al., 2020) . Another group analyzed the air quality impacts in São Paulo, where drastic reductions observed for NO, NO 2 , and CO by 77.3%, 54.3%, and 54.3%, respectively compared to the five-year monthly mean. On the other hand, the researchers observed an approximately 30% increase in ozone concentrations in urban areas (Nakada and Urban, 2020) . In Kazakhstan, PM 2.5 , CO, and NO 2 concentrations were reduced by 21%, 49%, and 35%, respectively during the first lockdown period between mid-March and mid-April compared to the average on the same days in 2018-2019 (Kerimray et al., 2020) . In India, where a nationwide lockdown was imposed initially for three weeks from 24th March to 14 th April 2020, PM 10 and PM 2.5 were reduced more than 50% compared to the pre-lockdown phase (Mahato et al., 2020) . Similarly, NO 2 and CO concentration levels were reduced by 52.68% and 30.35%, respectively. In Salé City, Morocco, 75%, 49%, and 96% reductions for PM 10 , SO 2 , and NO 2 were observed, respectively (Otmani et al., 2020) . In the United Arab Emirates, Nitrogen Dioxide (NO 2 ), Aerosol Optical Depth, and Surface Urban Heat Island Intensity variables were examined during lockdown periods and the levels declined by 23.7%, 3.7%, and 19.2%, respectively, compared to the same period in 2019 (Alqasemi et al., 2021) . The majority of these studies focus on the first lockdown period between March and April of 2020. In literature, mobility data is mainly used to investigate the spread of COVID-19 in societies (Chang, 2020) . In this study, we explore the continuous change in air quality between January 2020 and November 2020, especially March and April, in 81 cities of Turkey by an extensive database of ground monitoring stations. Besides, we investigate the associations between air quality parameters and mobility trends to understand the impact of lockdown measures on air quality. We obtained 24-h average PM 10 , PM 2.5, SO 2 , CO, O 3 , NOx, and NO 2 concentration data for 2015 through 2020 from the public air quality monitoring network designed by the Turkish Ministry of Environment and Urbanization (UHKIA, 2020). Each city has at least one station for air quality measurement (Fig. 1) . However, only PM 10 and SO 2 parameters were selected for the analysis due to being the only available parameters for 81 cities between January 2015 and November 2020. There are several missing data for other air quality parameters in smaller cities. We aim to investigate the change in all cities in TR compared to five years before lockdown; therefore, we only selected two parameters for this study. Mobility measure is the most important statistical parameter to define the distance and movement of a person in a society. The aim of mobility data is to capture general patterns in human mobility based on mobile device utilization. Even though the government set social distancing guidance and limitations, which is the same for all 81 cities of Turkey, we observed different behavior in each region. Therefore, we aim to capture actual patterns of human mobility by using datasets from Google (GoogleLLC, 2020). The importance of the dataset is showing a relative trend of how the movement of people in Turkey changes day after day. Google's mobility dataset is based on users' location information, which is classified into Workplaces, Residential, Retail and Recreation, Grocery and Pharmacy, Transit, Stations, Parks, and Restaurant/Café. We have used six categories, Workplaces, Residential, Grocery and Pharmacy, Transit, Stations, Parks, and Restaurant/Café, which can have direct effects on air quality. We have also obtained car-purchasing data from the Turkish Statistical Institute reports (TUIK, 2020) to interpret the results. We performed unpaired Welch's two-sample t-test analysis by Rsoftware (RStudio Team, 2020) in order to determine the significant reductions in PM 10 and SO 2 concentrations in 2020. In this analysis, we N.H. Orak and O. Ozdemir Environmental Research 197 (2021) 111018 assumed the samples follow normal distribution due to the sufficient number of observations in each sample. Then, we performed F-test to check the variance assumption of the t-test. F-statistic is defined as: where n1 and n2 are the sample sizes of the samples. The formula for t-statistics varies based on the variance equality. If the variance of both samples are statistically equal is the pooled variance. In this equation, X 1 and X 2 are the means of the samples, n1 and n2 are the sample sizes, s 2 1 and s 2 2 are the variances. If the variance equality is not satisfied, Here X 1 and X 2 are the means of the samples, n1 and n2 are the sample sizes, s 2 1 and s 2 2 are the variances. We divided the data set into two groups, mean of 2015-2019 and 2020. H 0 represents there is no significant effect of COVID-19 on air quality, where H 1 represent the significant difference between the means of the two groups. We performed Pearson's correlation analysis to measure the relationship between air pollutants and mobility data sets. Pearson's correlation coefficient is a number ranging between − 1 and +1, which shows a negative to a positive correlation between variables. The formula for Pearson's correlation is: is the covariance between two variables in the scenario, σ is the standard deviation of x 1 and x 2 . In this study, we aim to test the difference between PM 10 and SO 2 concentrations in 2020 and 2015-2019. Therefore, we performed unpaired or independent Welch's two-sample t statistic, which is based on comparing two groups based on their average. Another reason for selecting this parametric test rather than the nonparametric test is the sample size of the dataset allows us to approach normality by Central N.H. Orak and O. Ozdemir Environmental Research 197 (2021) 111018 Limit Theorem. The second aim is to check the relationship between PM 10 and SO 2 concentrations. For this purpose, we performed Pearson correlation coefficients between related variables to find out the strength and direction of the relationship if it exists. We compared daily average concentrations of PM 10 and SO 2 of 2020 to the average of the daily concentrations of the years 2015 through 2019 (Fig. S2 ). We did not take monthly averages to be able to provide the most accurate variability that represents the lockdown restrictions, which were updated weekly not monthly. In addition, we used the average of five-year data (2015-2019) to provide meaningful baseline data to be able to compare with the associated dates in 2020. Second, we investigated the association between air quality parameters and mobility data in each city. There is a general lower concentration trend for PM 10 throughout 2020 compared to the 2015-2019 average at the ground monitoring sites in Turkey (Fig. 2a) . However, we observe a sharp decrease in March as a result of lockdown measures. The average concentration for March 2020 is 43.75 μg/m 3 while the average concentration is 53.90 μg/m 3 for the baseline. We see the real impact of lockdown in the month of April, approximately 67% of cities had lower average PM 10 concentrations in April 2020 compared to the associated dates in all of the previous five years. The gap between the two trends narrows down during the summer months, which can be a result of ease in lockdown restrictions starting from June 2020. Fig. 2b shows the difference between the daily average concentrations of PM 10 (μg/m 3 ) between 2020 and 2015-2019. The median concentration for 2020 is significantly lower than 2015-2019. There are outliers for 2020 due to occasional very high concentrations related to sudden changes in restriction policy. Table S4 ). The variability for average 2020 decreased compared to the previous month. In May, we see a similar concentration level for 2020 but higher variability due to starting fewer restrictions by the end of the month (Fig. 3) . The stringency level decreased to 63.9 by June 2, 2020. As a result, we see an increasing trend in PM 10 concentrations in June with lower variability. Fig. 4 shows the spatial variability of average PM 10 concentrations. We see a general betterment in PM 10 concentrations across the country. Kastamonu city is excluded from 2015 to 2019; Kastamonu, Eskişehir, Uşak, Mersin, and Osmaniye cities were excluded from the 2020 data set due to data availability. The (Table 1) . We observe a more drastic change in April. In Hakkari city, the average PM 10 concentration decreased by 521% compared to the previous years in April (Fig. S2) . However, we see a significant increase in a few cities, such as Batman city (+53%). One of the reasons can be agricultural burning practices, which is common for example in Batman city. In metropolitan cities, the impact of stay-at-home order on PM 10 concentration is slightly lower but causes a significant decrease such as in Istanbul city (− 55%), which can be a direct result of reduced road traffic. On the other hand, we do not see this impact in the capital city, Ankara, where people's actual mobility is different from the state ordinances. This can be related to characteristics of the city where more state offices locate than private businesses. Fig. 5 shows the average daily SO 2 concentrations throughout 2020 compared to 2015-2019 mean across Turkey. We see a gradual decrease from March to May for both trends, which is parallel to the domestic heating season. The average SO 2 concentration for 2020 is slightly higher compared to the previous five years based on a two-tailed t-test (p < .02) where both datasets are right-skewed with outliers. The occurrence of more outliers is closely related to the sudden changes in mobility restrictions, which corresponds to high variation in data. In March and April 2020, the average concentrations of SO 2 are slightly higher than 2015-2019; however, the differences are not significant (Table S2) . On the other hand, we observe significantly higher concentrations in May and June of 2020 as shown in Fig. 6 (p < .001 ). The month of May serves as a transition period, as a result, we see higher Table 1 Average concentration of PM 10 (μg/m 3 ) from largest to smallest change by city for the years 2015-2019 and 2020. Change indicates the percent change in 2020 with respect to the 2015-2019 mean. N.H. Orak and O. Ozdemir Environmental Research 197 (2021) 111018 variability in data. The significant increase in May and June can be a result of behavioral change. People tend to use more private vehicles during the COVID-19 pandemic. Vehicle exhaust emissions are a major contributor to SO 2 emissions. According to Turkey Statistical Institute (TUIK), the number of the newly registered automobiles increased 55% between January-August in 2020 compared to 2019 (Table S3) . 48.5% is petrol and 42.9% is diesel-fueled vehicles out of the newly registered cars in 2020. As a consequence, we see higher average daily SO 2 concentration levels throughout 2020 (Fig. 5a ). There might be two reasons for high variability in May and June; first, weekly changes in lockdown measures, second is the movement of air masses before arriving at the monitoring sites. In terms of spatial variability, SO 2 reductions were recorded in 59% of the cities in Turkey (Fig. 7) , the highest decreases were observed in Afyon (− 291%), Ardahan (− 222%), Kars (− 203%), Denizli (− 161%), Yozgat (− 137%), Muş (− 127%), Karaman (− 126%), and Sinop (− 114%) cities (Table 2) . Noticeable increases occurred in Hakkari (+74%), Bursa (+72%), Tunceli (+71%), and Ordu (+69%) cities. Fig. 7 shows the average SO 2 concentrations between January-November in 2020 and 2015-2019 mean. We see a stronger change in April compared to the average concentration levels in 2020 (Table S5, Fig. S4 ). It is important to investigate the actual mobility of people in Turkey to capture the effect of COVID-19 on local air quality because there was no complete compliance with the government restrictions. We used six classes to analyze the relationship between mobility and PM 10 concentration levels; Restaurant/Café, Grocery and Pharmacy, Park, Transit, Workplaces, and Residential. We calculated the Pearson's productmoment correlation coefficient for the concentration levels and mobility data in each category. Table S5 shows the average mobility values of the six categories for Turkey and each city. The negative values represent the decrease in the percentage of mobility while the positive ones show the increase in the percentage of mobility with respect to baseline denoted by 0. NA values in the table mean that those cities have no information for the related category. Fig. S5 shows the average change in mobility for the whole country. Results show that the mobility in Restaurant/Café, Grocery and Pharmacy, Transit, Workplaces decreased in Turkey during the given time period while the mobility in both Parks and Residential increased as expected. On the other hand, we can see that the highest decrease in the mobility at Restaurant/Café is in Hakkari city while the lowest one is in Mugla city. The main differences between these two cities are population and characteristics. Mugla city is a holiday town with a population of more than three times of Hakkari city. If we look at the mobility in Grocery and Pharmacy, Tunceli has the most reduction in mobility while Agrı has the most increase. The most decrease in the mobility at Parks is in Gümüşhane while the highest increase is in Mugla. The results for the Transit category shows that the highest decrease occurred in Karabük while the lowest was in Table 2 Average concentration of SO 2 (μg/m 3 ) from largest to smallest change by city for the years 2015-2019 and 2020. Change indicates the percent change in 2020 with respect to the 2015-2019 mean. the eastern city of Van. The mobility in workplaces for provinces shows that Istanbul has the highest decrease while the Hakkari has the lowest one. Lastly, Istanbul has the highest positive change in Residential mobility while Artvin has the highest decrease. Fig. S6 shows that the Restaurant/Café mobility and the average concentration of PM 10 are positively correlated (r = 0.39, p < .001). We see a similar trend for the average concentration and Grocer-y&Pharmacy mobility (r = 0.32, p < .001), Parks (r = 0.18, p = .002), Transit (r = 0.43, p < .001), and Workplaces (r = 0.40, p < .001). On the other hand, we detect a negative correlation between Residential and mobility (r = − 0.32, p < .001), which is expected. The results show that the strongest positive impact on average SO 2 concentration is Restaurant/Café mobility (r = 0.32, p < .001), Transit mobility (r = 0.26, p < .001), and Workplaces mobility (r = 0.39, p < .001) (Fig. S7 ). We observe a negative correlation between stay-at-home and average SO 2 concentration (r = − 0.16, p = .007), which shows that the main source of SO 2 is related to traffic emissions. However, we need data on industrial activities for further investigation, because SO 2 is a common air pollutant to describe industry emissions. This study investigated the effects of COVID-19 restrictions on PM 10 and SO 2 concentrations across Turkey. All cities of Turkey (81) implemented the same lockdown measures. The state set the first mobility restrictions in March 2020; however, the rules have been updated weekly throughout the year. First, we analyzed the daily average concentrations of both pollutants between January and November of 2020, compared to the previous five years, 2015-2019 based on public air quality monitoring network data. We analyzed the correlation between mobility and air quality by obtained mobility data from Google. The presented results in this study are critically important for future investigations of air quality and mortality rates. Several studies explored the impact of air quality on mortality rates all around the world (Buonanno et al., 2020; Coker et al., 2020; Giani et al., 2020; Holshue et al., 2020; Liang et al., 2020; Liu et al., 2020; Ogen 2020) . Mortality rates by the city are not available yet in Turkey; however, we will be able to continue the next step when this information is available in 2021. As we know, this is the first study that investigates the effect of COVID-19 on ambient air quality across Turkey and the relationship with mobility. Also, 2020 provides a rare opportunity to experiment and analyze how human behavior can affect PM 10 and SO 2 concentrations. The results can help policymakers to understand the impact of lockdown measures on air quality in different cities with different characteristics, which provides a comprehensive air quality picture of the region that can be used in future studies on air quality and infectious diseases. Nur H. Orak: Conceptualization, Methodology, Formal analysis, Writingoriginal draft, Writingreview & editing, Ozancan Ozdemir: Methodology, Formal analysis, Visualization, Writingreview & editing. Impact of COVID-19 lockdown upon the air quality and surface urban heat island intensity over the United Arab Emirates Changes in air quality and human mobility in the USA during the COVID-19 pandemic Covid-19 lockdown effects on air quality by NO 2 in the cities of Barcelona and Madrid (Spain) Correlation between environmental pollution indicators and COVID-19 pandemic: a brief study in Californian context Does lockdown reduce air pollution? Evidence from 44 cities in Northern China Koronavirüs -okullar 16 mart'tan itibaren 1 hafta tatil, spor karşılaşmaları Nisan sonuna dek seyircisiz Changes in U.S. Air pollution during the COVID-19 pandemic Estimation of airborne viral emission: quanta emission rate of SARS-CoV-2 for infection risk assessment Mobility network models of COVID-19 explain inequities and inform reopening Air pollution reduction and mortality benefit during the COVID-19 outbreak in China The effects of atmospheric stability with low wind speed and of air pollution on the accelerated transmission dynamics of COVID-19 Effects of the spread of COVID-19 on public health of polluted cities: results of the first wave for explaining the dejà vu in the second wave of COVID-19 pandemic and epidemics of future vital agents An Index to Quantify Environmental Risk of Exposure to Future Epidemics of the COVID-19 and similar Viral Agents: Theory and Practice The effects of air pollution on COVID-19 related mortality in northern Italy The impact of COVID-19 partial lockdown on the air quality of the city of Rio de Janeiro Short-term and long-term health impacts of air pollution reductions from COVID-19 lockdowns in China and Europe: a modelling study COVID-19 Community Mobility Reports Variation in Government Responses to COVID-19 First case of 2019 novel coronavirus in the United States Johns Hopkins University Coronavirus Resource Center Assessing air quality changes in large cities during COVID-19 lockdowns: the impacts of traffic-free urban conditions in Almaty Drop in urban air pollution from COVID-19 pandemic: policy implications for the megacity of Sao Paulo Global socio-economic losses and environmental gains from the coronavirus pandemic Urban air pollution may enhance COVID-19 case-fatality and mortality rates in the United States Effects of COVID-19 lockdown on global air quality and health Genomic characterisation and epidemiology of 2019 novel coronavirus: implications for virus origins and receptor binding Effect of lockdown amid COVID-19 pandemic on air quality of the megacity Delhi COVID-19 pandemic and environmental pollution: a blessing in disguise? Covid-19 pandemic: impacts on the air quality during the partial lockdown in Sao Paulo state Airborne nitrogen dioxide plummets over China. Earth Observatory. North American space agency Assessing nitrogen dioxide (NO 2 ) levels as a contributing factor to coronavirus (COVID-19) fatality Impact of COVID-19 lockdown on PM 10 , SO 2 and NO 2 concentrations in sale city (Morocco) Response of major air pollutants to COVID-19 lockdowns in China RStudio Team, 2020. RStudio. Integrated Development for Impact of COVID-19 outbreak measures of lockdown on the Italian carbon footprint COVID-19 and air pollution and meteorology-an intricate relationship: a review Türkiye'de ilk koronavirüs vakası tespit edildi TUIK, 2020. Turkey Statistical Institute Turkish Ministry of Environment and Urbanization Director-general's opening remarks at the media briefing on covid-19-11 WHO, 2020b. Naming the Coronavirus Disease (Covid-19) and the Virus that Causes it 2020 A novel coronavirus from patients with pneumonia in China The authors declare that they have no actual or potential competing financial interests. Supplementary data to this article can be found online at https://doi. org/10.1016/j.envres.2021.111018.