Next Article in Journal
Assessing the Level of Knowledge, Beliefs and Acceptance of HPV Vaccine: A Cross-Sectional Study in Romania
Previous Article in Journal
Simple and Autonomous Sleep Signal Processing System for the Detection of Obstructive Sleep Apneas
Previous Article in Special Issue
Relationship between Meteorological and Air Quality Parameters and COVID-19 in Casablanca Region, Morocco
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
IJERPH
Volume 19
Issue 11
10.3390/ijerph19116937
ijerph-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Machine Learning and Meteorological Normalization for Assessment of Particulate Matter Changes during the COVID-19 Lockdown in Zagreb, Croatia †

by
Mario Lovrić
1,2,*,
Mario Antunović
3,
Iva Šunić
2,
Matej Vuković
4,
Simonas Kecorius
5,
Mark Kröll
1,
Ivan Bešlić
6,
Ranka Godec
6,
Gordana Pehnec
6,
Bernhard C. Geiger
1,
Stuart K. Grange
7,8 and
Iva Šimić
6,*
1
Know-Center, Inffeldgasse 13, 8010 Graz, Austria
2
Institute for Anthropological Research, Gajeva 32, 10000 Zagreb, Croatia
3
Ascalia d.o.o., Ulica Trate 16, 40000 Čakovec, Croatia
4
Pro2Future GmbH, Inffeldgasse 25F, 8010 Graz, Austria
5
Institute of Epidemiology, Helmholtz Zentrum München, Ingolstädter Landstr. 1, 85764 Neuherberg, Germany
6
Environmental Hygiene Unit, Institute for Medical Research and Occupational Health, Ksaverska cesta 2, 10000 Zagreb, Croatia
7
Empa, Swiss Federal Laboratories for Materials Science and Technology, 8600 Dübendorf, Switzerland
8
Wolfson Atmospheric Chemistry Laboratories, Department of Chemistry, University of York, York YO10 5DD, UK
*
Authors to whom correspondence should be addressed.
†
This paper is an extended version of the work presented at the 7th International Conference on Environmental Pollution, Treatment and Protection (ICEPTP’22).
Int. J. Environ. Res. Public Health 2022, 19(11), 6937; https://doi.org/10.3390/ijerph19116937
Submission received: 30 March 2022 / Revised: 30 May 2022 / Accepted: 2 June 2022 / Published: 6 June 2022
(This article belongs to the Special Issue The Role of Air Pollution during the COVID-19 Pandemic)
Browse Figures
Versions Notes

Abstract

:
In this paper, the authors investigated changes in mass concentrations of particulate matter (PM) during the Coronavirus Disease of 2019 (COVID-19) lockdown. Daily samples of PM1, PM2.5 and PM10 fractions were measured at an urban background sampling site in Zagreb, Croatia from 2009 to late 2020. For the purpose of meteorological normalization, the mass concentrations were fed alongside meteorological and temporal data to Random Forest (RF) and LightGBM (LGB) models tuned by Bayesian optimization. The models’ predictions were subsequently de-weathered by meteorological normalization using repeated random resampling of all predictive variables except the trend variable. Three pollution periods in 2020 were examined in detail: January and February, as pre-lockdown, the month of April as the lockdown period, as well as June and July as the “new normal”. An evaluation using normalized mass concentrations of particulate matter and Analysis of variance (ANOVA) was conducted. The results showed that no significant differences were observed for PM1, PM2.5 and PM10 in April 2020—compared to the same period in 2018 and 2019. No significant changes were observed for the “new normal” as well. The results thus indicate that a reduction in mobility during COVID-19 lockdown in Zagreb, Croatia, did not significantly affect particulate matter concentration in the long-term..

1. Introduction

Particulate matter (PM) is recognized as one of the major air pollutants affecting human health. Particle size plays an important role in determining pollutant respiratory deposition and thus potential health risks. Airborne particles PM10 (with aerodynamic diameter less than 10 μm) and especially its smaller fractions (e.g., PM2.5—with aerodynamic diameter less than 2.5 μm and PM1—with aerodynamic diameter less than 1 μm) are known to effectively enter the human body, e.g., trachea (upper throat) or bronchi, and even reach all the way down to the alveoli in the lungs, where it can penetrate from the lung alveoli into the blood [1,2]. In general, the smaller the particle size, the greater the adverse health effect [3,4,5]. Therefore, further reduction of PM pollution both in developed and developing countries has the potential to improve both life quality and expectancy. To better understand sources, as well as environmental and health impacts of air pollution, long-term measurement data sets are used in source appointment, epidemiological, and air quality studies [6]. On the other hand, short-term traffic bans can be used to pinpoint pollution contributors and raise awareness of air quality problems [7]. Ironically, besides causing worldwide health and economical disturbance, the current COVID-19 pandemic has also provided means to investigate air pollution [8]. Published evidence on the impact of the COVID-19 lockdown on the concentration of ambient air pollutants highlights the importance of transport and industrial activities [9,10]. For example, there is clear evidence for reduced gaseous (e.g., nitrogen dioxide (NO2)) and particulate pollutant concentrations in some urban areas, which can be linked to reduced transportation due to COVID-19 [11,12,13,14]. In contrast, the European Environment Agency reported that a consistent reduction of PM2.5 cannot be seen in European cities during the lockdown period [15]. The main reasons could be that local pollution sources are more various, including not only industrial activities and road traffic, but also the combustion of different fuels for the heating, as well as the formation of secondary aerosols [16]. Furthermore, it is not entirely clear how lockdown period pollutant concentrations depend on other effects/confounders that should be accounted for, e.g., weather effects [13]. The methods used in lockdown-related air pollution studies differ significantly. Many studies have tried to assess the influence of lockdown measures on air pollutant concentrations by a simple comparison of basic statistical parameters during the lockdown with the same period a year before or up to five years before [17,18,19,20,21,22,23]. However, in such studies, the variability of meteorological factors between years as well as long-term trends are not considered. There are not many studies that include more sophisticated statistical and modelling tools [13,24,25,26,27]. A summary of the review will be presented in the discussion section. In machine learning prediction of air pollutant concentration, one often assumes that the concentration (dependent variable) is a function of temporal and meteorological determinants (independent variables) [28,29,30]. With that in mind, one can employ complex algorithms which are considering possible non-linear relationships within the data and present true influencing determinates based on inference. Two commonly seen (non-linear) machine learning algorithms are Random Forest (RF) regression and neural networks (NN). In previous work [13] RF regression was used to predict pollutant concentrations during the lockdown in Graz and presented the advantages of utilizing such methods over the historical comparison of pollution. Similarly, RF was used to assess changes in pollutant levels during different stages of lockdown in Los Angeles by comparing predicted concentrations under different traffic emission scenarios [31]. A similar approach was used by Brancher [32] who refers to baseline models (non-lockdown periods) as “Business as usual” scenarios. The model describes hourly-averaged concentrations per pollutant and monitoring station to investigate air quality changes before and after lockdown and to verify the models’ predictive skill to reproduce the pollutant measurements. A NN approach was used to investigate whether changes in air quality in Nigeria were caused primarily by the lockdown. In this case, monthly average values of ground-level fine aerosol optical depth (AODf) across Nigeria from 2001 to 2020 were used [33]. Another method for the assessment of air pollution during the lockdown period is the difference-in-differences (DID) model. Xu et al. [34] used this method to evaluate air pollutants and air quality before and during the lockdown. The DID model calculates the effect of treatment (independent variable) on outcome (dependent variable) by comparing the average changes in each of the groups. In this case, the outcome is the level of air pollution. Control variables such as temperature, humidity, wind speed, etc. are also included. The model considers whether the lockdown was enforced or not for each date and based on this calculates relative changes in air pollution levels. The study by Gope et al. [35] used the Air Quality Index (AQI), which is calculated from the concentration of the pollutants, to analyse the impact of lockdown on the environment. Comparison of the AQI for these periods showed that most cities reduced their pollution. The pandemic of COVID-19 caused many changes in human activities, not only during the lockdown but also in the months following it. For example, working from home and virtual meetings caused lower mobility and, in some areas, even became a regular practice. There are mixed results published regarding the lockdown and “new normal” effects on particulate matter. While many methods are being used, there is a lack of a standardized approach for understanding these phenomena. In this work, the authors present an assessment of particulate matter in three mass fractions (PM1, PM2.5 and PM10) based on daily measurements over a long period of 12 years at one urban background location. Previous PM measurements at the same location have shown significant air pollution during winter months [36,37,38]. Data from air quality monitoring stations [39] shows that in a few years prior to the COVID-19 pandemic, levels of PM10 and PM2.5 were below regulatory limits set by Croatian and EU legislation for protection of human health. Annual averages of PM10 and PM2.5 were below 40 µg/m3 and 25 µg/m3, respectively. Daily limit value for PM10 is 50 µg/m3 and should not be exceeded more than 35 times during the calendar year. This criterion was met at that location since 2017. However, considering new WHO guidelines [40] which, in the light of recent scientific evidence, suggests much lower limit values for both, PM10 and PM2.5 fractions, it is necessary to apply new measures for reducing air pollution. With the aim to protect people’s wellbeing, it is important to know the main pollution sources and the efficiency of implemented measures. The intention of this paper was to examine whether a reduction in mobility during the COVID-19 lockdown caused changes in PM levels. The hypotheses are that the lockdown and the “new normal” both show reductions in particulate matter concentration. A reduction during the “new normal” is hypothesized due to a restriction on travel which affects Croatia’s tourism and more working from home. To test these hypotheses, these periods were evaluated by using ANOVA on the normalized (de-weathered) data.

2. Materials and Methods

2.1. Particulate Matter and Meteorological Measurements

Aerosol PM concentrations were measured in Zagreb, Croatia, at a sampling site located in the northern, residential part of the city (45°50′7″ N, 15°58′42″ E, 116 m a.s.l.,).
The area is characterized by modest traffic and population density. The household heating (gas and/or wood) season usually starts in October and lasts until April. The PM samplers (Sven Leckel, engineering office, Berlin, Germany) were positioned at about 20 m from the nearest street. Twenty-four-hour samples of PM1, PM2.5 and PM10 fractions have been collected continuously every day on quartz filters (47 mm in diameter) during 12-years period (2009–2020). PM mass concentrations were determined gravimetrically (Mettler TOLEDO MX5 balance, Greifensee, Switzerland) according to the EN 12341:1998 and EN 14907:2005 standards from 2009–2014 and EN 12341:2014 standard from 2015–2020. Before and after the sampling, filters were conditioned at a constant temperature (20 ± 1 °C) and relative air humidity (45–50% RH) for 48 h. Meteorological parameters (temperature, RH, wind speed and direction, pressure, and precipitation) were obtained from the Croatian Meteorological and Hydrological Service. For the input of data used in this study an explorative plot of the particle mass concentrations over time used in this study is shown in Figure 1. The collected data can be found at ref. [41].

2.2. Data Processing and Model Training

The dataset used in this analysis is air pollution data collected over the period of 12 years (2009–2020) in a daily frequency (daily average) in Zagreb, Croatia. The dataset contains the particulate matter mass concentrations of PM1, PM2.5 and PM10 and temporal information such as: day of week, Julian date (days counted from 1 January 1970), month, year, holiday tag, etc. In addition, temporally aligned meteorological influences were added into the dataset. These consist of maximum daily temperature (T), minimal daily T, difference of max and min T, average T, maximum and minimum daily pressure (p), difference of max and min p, average p, maximum daily relative humidity (RH), minimum daily RH, difference of max and min daily RH, average RH, wind speed and precipitation. Temporal and meteorological variables are given as independent or predictive variables. To retain a high amount of data for machine learning (ML), missing datapoints were imputed with backfill strategy (missing values are filled with the ones from the following day). Python programming language (www.python.org accessed on 1 February 2022, v3.7.10) was used for analysis, while data processing and model training follows the method and the process described in [13,30]. It is assumed that the concentrations of particulate matter (PM1, PM2.5, PM10; dependent variables) can be modelled through temporal and meteorological variables as independent ones, previously listed. To model the air pollution, RF [42] and LightGBM [43] methods were used. RF is an ensemble ML algorithm which consists of many individual decision trees and applies bootstrap aggregation (bagging) and feature randomness techniques in building each decision tree. Like RF, LightGBM is an ensemble method that relies on tree-based learning but utilizes gradient boosting techniques as well as different tree-building techniques. RF and LightGBM are non-parametric and as such require no formal distributional assumptions enabling these methods to deal with skewed and multi-modal data. Consequently, these methods are well suited for modelling challenging phenomena such as air pollution, but also other settings as outlined in several studies [13,28,30,44,45]. Air pollution ML models were trained for PM1, PM2.5 and PM10 respectively, with their daily concentrations representing target (dependent, predicted) variables. Following the method outlined in previous studies [46,47], hyperparameters of these regression models were optimized through 10-fold cross-validation with Bayesian optimization. The training dataset (TDS) consists of data between 1 January 2009 and 31 December 2019 while data from 2020 was split into several smaller datasets indicating different validation and interest periods. Finally, the models were tested on their generalization performance on MVS which is reported in the Results section.
The first coronavirus patient was confirmed in Croatia on 25 February 2020. The Coronavirus Disease of 2019 (COVID-19) disease in Croatia was announced by the government on the 11 March 2020. In the second half of March, all public events and gatherings were canceled and all non-essential activities (shopping centers, bars and nightclubs, restaurants, cinemas, libraries, gyms, sports centers and sport events, dance schools, children’s workshops, religious and other public gatherings) were closed. On 19 March the decision was passed to restrict staying on the streets, squares, and other public places. Two days later it was followed by suspensions of public transport, suspension of intercity lines for trains and buses, and shortened working hours of shops and post offices. Finally, on 23 March a ban on leaving one’s place of residence or permanent residence was set and only persons and activities important for the movement and movement of goods were exempt. Relaxation of strict measures started on 23 April and was carried out in three phases till 11 May 2020: were relaxed between 27 April 27 and 11 May 2020 after which the bans were lifted. Sources regarding the given information can be found following references [48,49,50,51]. These interest periods and overall timeframes for 2020 are depicted in Figure 2 for the sake of simplifying the timeframes.
A priori, the 2020 data on “out-of-ordinary” events which might have affected the pollutant concentrations were analyzed. “Out-of-ordinary” events include: a 5.3-magnitude earthquake in Zagreb (22 March 2020), long-range transport of desert dust events on 26–30 March [13,52] and construction works near the measurement site in August 2020.
As shown in Figure 3, these events clearly disrupt normal concentrations and introduce bias in the models. Therefore, the respective timeframes were excluded from the presented analyses. As a result, even though the lockdown timeframe lasted longer than the one given as LDS, several dust events needed to be excluded. A subset from MVS for comparison to LDS (comparison set or CS) was split, which is set between 3 January and 29 February 2020. The subset is shorter than MVS due to several construction activities at the site in March. LDS and NNS present the timeframe in focus of this lockdown pollution investigation. MVS was used to better understand the model generalization.

2.3. Meteorological Normalization (De-Weathering)

In this work, the authors followed the methodology from Grange et al. [28,29] for meteorological normalization of the daily particulate matter time series. Meteorological normalization was achieved by firstly creating an ML model per pollutant (particulate matter concentration) that generalizes well on unseen data. In the next step, all predictive variables (except Julian day) are repeatedly randomly sampled without replacement and used to predict pollutant concentration using the individual trained RF models.
The procedure of meteorological normalization removes the short-term variation in the time series. Reasoning for this procedure is that pollutant data must be corrected for meteorological and temporal effects which are changing over time and can therefore affect pollutant concentration. Herein, an example is shown by means of changing temperature and precipitation through the given years (Figure 4). One can observe that there was a trend in Zagreb, Croatia towards higher temperature and precipitation at the given monitoring station.
To normalize, the model predictions for each pollutant for 100 random samples were then averaged into the normalized time series (normalized PM1, PM2.5, PM10). The complete procedure of data processing, machine learning model training and meteorological normalization is presented in Figure 5.

3. Results

The model scores by means of root-mean-square-error, R2 scores and chosen models by means of MVS performance are shown in Table 1. Based on R2 scores, all three models show good predictive values (good generalization). The prediction quality in RMSE is similar for PM10 and PM2.5, while it increases for PM1. When comparing the R2 scores to the author’s previous work [13], the observed values in this study suggest a reasonably good generalization with R2 scores above 0.77.
Once models were trained, the data was normalized (de-weathered) as described in Section 2.3 and depicted in Figure 5. To evaluate change in airborne pollution concentrations due to the lockdown, yearly trends were assessed by means of median of the normalized time series. Three timeframes were compared (Figure 6), namely the months of January and February together (CS) and June together with July (NNS) which are considered to be the new normal (post-lockdown changes). Normalized time series during the months of April (LDS) every year is given in Figure 6. In Figure 6a, for CS (January and February, pre-lockdown reference) there is a continuous reduction from 2009 to 2017 for all size fractions of PM. However, starting with 2017 an unexpected increase for the normalized values can be observed for PM10 and a slight increase from 2019 to 2020 for PM1, while PM2.5 has a steady level from 2017 onwards. During NNS (the months of June and July, Figure 6b) a reduction in pollution compared to previous years was expected due to lower tourism rate and the travel ban as well as working from home and many isolations. Instead, the observed pattern for NNS is like CS, with PM2.5 showing similar levels from 2018–2020, while PM1 and PM10 show slight increases. In the case of April, when the lockdown took place (LDS; boxplots in Figure 7), PM2.5 shows a steady decrease from 2009 onwards with a visually insignificant decrease from 2019 to 2020. PM10 mass concentration shows a similar trend with visually insignificant changes from 2017 onwards. The results are slightly different from PM1 which shows mixed periods of decrease and increases from 2009 onwards. Even though there seems to be a huge drop from 2018 to 2019, there is an increase from 2019 to 2020. ANOVA was utilized to assess whether the changes observed in normalized PM-time series are significant.
For each of PM1, PM2.5 and PM10 a test was created with the year as categorical variable (2018, 2019, 2020) as the independent variable and the normalized concentrations as the dependent variable. ANOVA showed that yearly changes for April do not show statistical significance for any of the PM fractions—PM1 (p = 0.10), PM2.5 (p = 0.47) and PM10 (p = 0.76). Similar results were obtained for the NNS (new normal) with PM1 (p = 0.26), PM2.5 (p = 0.81) and PM10 (p = 0.72).

4. Discussion

Table 2 presents the literature findings on the influence of lockdown on air pollutant’s levels and methods applied. It is evident that most of the studies did not include long-term PM measurements and were mostly focused on the comparison of the lockdown period with the same period in 2019. Furthermore, they do not account for the effect of year-to year variability and seasonal variability caused by meteorology and weather conditions, which was the reason for using normalization [28,29]. However, even normalization can suffer from a lack of data or model quality. For this reason, we utilized, besides RF, also LightGBM and Bayesian optimization for tuning all models to improve model accuracy. Comparing to data sources listed in Table 2, this work has a fair amount of data with a collection starting from 2009, which the authors deem to be needed for such an undertaking.
Regarding the changes in air pollution, the results obtained in this study are in line with the other published work [13,14,34,73]. In a previous study carried out in Zagreb, average concentrations of NO2, PM1 and PAHs in PM1 during the lockdown period were compared with the average concentrations for the same period in 2019. It was found that during lockdown at the traffic measuring site concentrations decreased by 35% for NO2 and PM1 compared with the same period in 2019. However, at the urban background measuring site NO2 decreased by 27% while PM1 levels remained like the year before [18].
There is mixed evidence for whether PM mass concentrations were truly affected by the lockdown, i.e., change in mobility, at least for the sites not heavily affected by traffic. In the authors previous work [13], it was found that a drop in PM10 mass concentration during the lockdown in Graz, Austria is rather inconsistent when comparing it to gaseous pollutants (e.g., NO2). A drop in NO2 gas concentrations was found to be around 40%, while PM10 mass concentration decrease was in the range of 6–14% compared to previous years. However, these analyses were conducted without meteorological normalization. Given the comparison to the observed reduction in traffic, it can be assumed that the reduced traffic was not a dominant contributing factor to the changes in PM10 atmospheric load. Xu et al. [34] who were using the difference-in-differences method to compare air pollution before and during the lockdown in China, found no change in PM10 and PM2.5 concentrations. This finding implies that traffic cannot be considered the main PM source. Although particulate matter concentrations show short-term declines within three days after the lockdown policies were implemented, when compared to the measurements from different monitoring sites that acted as a control group, the lockdown effects were not as pronounced. The given results are in accordance also with a study by Etchie et al. [33] who observes no effect on PM by the lockdown. A study that also used RF and meteorological normalization shows only a moderate decrease for PM10 [73]. The results imply that in Zagreb, Croatia, traffic is not the main contributor in such a site to air pollution by means of particulate matter. This is essential to the discussion of which factors contribute to public health, such as traffic, residential heating, and urban planning. If traffic is not the main contributor to particulate matter, regardless of its contribution by means of nitrogenous oxides [13] and polycyclic aromatic hydrocarbon (PAH) [18] then local policymakers need to invest more efforts in understanding other contributors and prevention policies which will improve air quality.

Limitations and Future Work

First and foremost, there are many contributing factors to atmospheric PM levels, such as secondary particle formation and long-range transport of particulate matter which are here not considered. Beyond those, there are also solar irradiation, traffic density, contribution from resuspension of road dust, etc.
The height of the boundary layer changes significantly during the day. In summer, this is more pronounced because the sun warms the substrate, vertical currents appear, and the height of the boundary layer increases. As night approaches, the height of the boundary layer decreases. Since in this work the authors used the daily values of the concentrations of suspended particles, this daily course is not distinguishable in the given data on the concentration of PM. Furthermore, in the city of Zagreb, a high stability of the boundary layer of the atmosphere was observed during the cold part of the years when increased concentrations of PM were also observed. It should be emphasized that Zagreb is surrounded by Medvednica mountain on the north and the river Sava on the south, and there are no major industrial cities along these routes. East and west of Zagreb there are only suburban settlements (Zaprešić, Samobor, Sesvete, Dugo Selo) which are residential centers without pronounced dominant sources, and this contribution is not included in this paper.
Another limiting factor is also the data frequency and measurement method. This work is based on daily averages of particulate matter and meteorology which are rough estimates of pollution trends and their dependence on meteorology. Even though the models have good generalization with R2 scores above 0.77, using these predictions will lead to error propagation since the models do not cover all the variance.
In the future, the authors aim to shift the focus towards PM composition and source apportionment regarding changes during the lockdown which helps understand the contributions. Furthermore, the authors intend to apply recent in-house research on the intersection of physics-based and machine learning based models (so call physics-inspired machine learning) such as seen in recent research [45,74] where models’ accuracy profits from the combination of the two worlds.

5. Conclusions

In this work the effects of the lockdown on PM mass concentration in Zagreb, Croatia were evaluated. The authors hypothesized that the COVID-19 lockdown (April 2020) and the “new normal” (June, July 2020) would both exhibit a decrease in PM1, PM2.5, and PM10 mass concentrations due to changed human behavior and mobility. To investigate the anticipated decrease in PM mass concentrations, machine learning, by means of Random Forests (RF) and LightGBM (LGB) were utilized and combined with meteorological normalization. The RF and LGB models trained in this study exhibited a reasonably good generalization on the test set (R2 scores > 0.77). By using normalization, the trend component of the PM mass concentration was extracted and compared it pre-, during, and post-lockdown timeframes. The results by means of normalized concentrations show that over the course of 2009–2017/2018, for the city of Zagreb at an urban location, PM mass concentrations dropped, however, no significant changes were observed in PM mass concentrations due to the lockdown or post-lockdown events. Besides that, given that NO2 reductions were earlier observed at the same site, one can speculate that a reduction in mobility did not affect particulate matter to a significant extent at this specific site.

Author Contributions

Conceptualization, M.L. and I.Š. (Iva Šimić); Data curation, M.A., I.Š. (Iva Šunić), I.B., R.G. and I.Š. (Iva Šimić); Formal analysis, M.L., M.A., M.V. and R.G.; Funding acquisition, M.K. and G.P.; Investigation, I.Š. (Iva Šunić), I.B., R.G., S.K.G. and I.Š. (Iva Šimić); Methodology, S.K., G.P., B.C.G. and S.K.G.; Project administration, M.K. and G.P.; Supervision, B.C.G.; Validation, I.B.; Visualization, M.L., M.A. and M.V.; Writing—original draft, M.L., I.Š. (Iva Šunić), M.V., S.K., G.P., S.K.G. and I.Š. (Iva Šimić); Writing—review & editing, S.K., M.K., R.G., B.C.G. and S.K.G. All authors have read and agreed to the published version of the manuscript.

Funding

The modelling parts of the presented work were developed within two projects 1) HPC and BigData Technologies for Global Systems (HiDALGO), under grant agreement No. 824115 and 2) HEAt -Holistische digitale Infrastruktur durch optimierte Erfassungs- und Analysemethoden des Verkehrsgeschehens. M.V. has been supported by the FFG-COMET-K1 Center “Pro²Future” (Products and Production Systems of the Future), Contract No. 881844. S.K.G. is supported by the Swiss Federal Office for the Environment (FOEN) and the Natural Environment Research Council (NERC) while holding associate status at the University of York.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Data is available at [41].

Acknowledgments

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study, analyses, or interpretation of data; in the writing of the manuscript, or in the decision to publish the results.

References

  1. Kim, K.H.; Kabir, E.; Kabir, S. A review on the human health impact of airborne particulate matter. Environ. Int. 2015, 74, 136–143. [Google Scholar] [CrossRef] [PubMed]
  2. Anderson, J.O.; Thundiyil, J.G.; Stolbach, A. Clearing the Air: A Review of the Effects of Particulate Matter Air Pollution on Human Health. J. Med. Toxicol. 2012, 8, 166–175. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  3. Li, N.; Chen, G.; Liu, F.; Mao, S.; Liu, Y.; Hou, Y.; Lu, Y.; Liu, S.; Wang, C.; Xiang, H.; et al. Associations of long-term exposure to ambient PM1 with hypertension and blood pressure in rural Chinese population: The Henan rural cohort study. Environ. Int. 2019, 128, 95–102. [Google Scholar] [CrossRef] [PubMed]
  4. Yang, M.; Guo, Y.M.; Bloom, M.S.; Dharmagee, S.C.; Morawska, L.; Heinrich, J.; Jalaludin, B.; Markevychd, I.; Knibbsf, L.D.; Lin, S.; et al. Is PM1 similar to PM2.5? A new insight into the association of PM1 and PM2.5 with children’s lung function. Environ. Int. 2020, 145, 106092. [Google Scholar] [CrossRef] [PubMed]
  5. Jakovljević, I.; Sever Štrukil, Z.; Godec, R.; Bešlić, I.; Davila, S.; Lovrić, M.; Pehnec, G. Pollution Sources and Carcinogenic Risk of PAHs in PM1 Particle Fraction in an Urban Area. Int. J. Environ. Res. Public Health 2020, 17, 9587. [Google Scholar] [CrossRef] [PubMed]
  6. Stojić, A.; Stojić, S.S.; Reljin, I.; Čabarkapa, M.; Šoštarić, A.; Perišić, M.; Mijić, Z. Comprehensive analysis of PM10 in Belgrade urban area on the basis of long-term measurements. Environ. Sci. Pollut. Res. 2016, 23, 10722–10732. [Google Scholar] [CrossRef]
  7. Wiedensohler, A.; Andrade, M.; Weinhold, K.; Müller, T.; Birmili, W.; Velarde, F.; Moreno, I.; Forno, R.; Sanchez, M.F.; Laj, P.; et al. Black carbon emission and transport mechanisms to the free troposphere at the La Paz/El Alto (Bolivia) metropolitan area based on the Day of Census (2012). Atmos. Environ. 2018, 194, 158–169. [Google Scholar] [CrossRef]
  8. Lovrić, M.; Antunović, M.; Šunić, I.; Vuković, M.; Kecorius, S.; Kröll, M.; Bešlić, I.; Šimić, I.; Pehnec, G. Insignificant Changes in Particulate Matter during the COVID-19 Lockdown: A Machine Learning Study in Zagreb, Croatia. In Proceedings of the Proceedings of the 7th World Congress on Civil, Structural, and Environmental Engineering (CSEE’22); Naggar, H., El Barros, J., Eds.; INTERNATIONAL ASET INC: Lisbon, Portugal, 2022. [Google Scholar] [CrossRef]
  9. Liu, F.; Wang, M.; Zheng, M. Effects of COVID-19 lockdown on global air quality and health. Sci. Total Environ. 2021, 755, 142533. [Google Scholar] [CrossRef]
  10. Singh, R.P.; Chauhan, A. Impact of lockdown on air quality in India during COVID-19 pandemic. Air Qual. Atmos. Health 2020, 13, 921–928. [Google Scholar] [CrossRef]
  11. Donzelli, G.; Cioni, L.; Cancellieri, M.; Morales, A.L.; Suárez-Varela, M.M.M. The effect of the COVID-19 lockdown on air quality in three italian medium-sized cities. Atmosphere 2020, 11, 1118. [Google Scholar] [CrossRef]
  12. Menut, L.; Bessagnet, B.; Siour, G.; Mailler, S.; Pennel, R.; Cholakian, A. Impact of lockdown measures to combat COVID-19 on air quality over western Europe. Sci. Total Environ. 2020, 741, 140426. [Google Scholar] [CrossRef] [PubMed]
  13. Lovrić, M.; Pavlović, K.; Vuković, M.; Grange, S.K.; Haberl, M.; Kern, R. Understanding the true effects of the COVID-19 lockdown on air pollution by means of machine learning. Environ. Pollut. 2020, 274, 115900. [Google Scholar] [CrossRef] [PubMed]
  14. Grange, S.K.; Lee, J.D.; Drysdale, W.S.; Lewis, A.C.; Hueglin, C.; Emmenegger, L.; Carslaw, D.C. COVID-19 lockdowns highlight a risk of increasing ozone pollution in European urban areas. Atmos. Chem. Phys. 2021, 21, 4169–4185. [Google Scholar] [CrossRef]
  15. Agency, E.E. Air Quality and COVID-19—European Environment Agency. Available online: https://www.eea.europa.eu/themes/air/air-quality-and-covid19 (accessed on 1 May 2022).
  16. Pinto, J.N.; Garcia, M.; Borrega, J.; Coelho, L.M.R. Airborne PM Impact on Health, Overview of Variables, and Key Factors to Decision Making in Air Quality. In Air Pollution-Monitoring, Quantification and Removal of Gases and Particles; IntechOpen: London, UK, 2018. [Google Scholar] [CrossRef] [Green Version]
  17. Ljubenkov, I.; Haddout, S.; Priya, K.L.; Hoguane, A.M. SARS-CoV-2 epidemic: Changes in air quality during the lockdown in Zagreb (Republic of Croatia). Toxicol. Environ. Chem. 2020, 102, 302–303. [Google Scholar] [CrossRef]
  18. Jakovljević, I.; Štrukil, Z.S.; Godec, R.; Davila, S.; Pehnec, G. Influence of lockdown caused by the COVID-19 pandemic on air pollution and carcinogenic content of particulate matter observed in Croatia. Air Qual. Atmos. Health 2021, 14, 467–472. [Google Scholar] [CrossRef]
  19. Dimovska, M.; Gjorgjev, D. The effects of COVID-19 lockdown on air quality in Macedonia. Open Access Maced. J. Med. Sci. 2020, 8, 353–362. [Google Scholar] [CrossRef]
  20. Gualtieri, G.; Brilli, L.; Carotenuto, F.; Vagnoli, C.; Zaldei, A.; Gioli, B. Quantifying road traffic impact on air quality in urban areas: A COVID19-induced lockdown analysis in Italy. Environ. Pollut. 2020, 267, 115682. [Google Scholar] [CrossRef]
  21. Baldasano, J.M. COVID-19 lockdown effects on air quality by NO2 in the cities of Barcelona and Madrid (Spain). Sci. Total Environ. 2020, 741, 140353. [Google Scholar] [CrossRef]
  22. Wyche, K.P.; Nichols, M.; Parfitt, H.; Beckett, P.; Gregg, D.J.; Smallbone, K.L.; Monks, P.S. Changes in ambient air quality and atmospheric composition and reactivity in the South East of the UK as a result of the COVID-19 lockdown. Sci. Total Environ. 2021, 755, 142526. [Google Scholar] [CrossRef]
  23. Ropkins, K.; Tate, J.E. Early observations on the impact of the COVID-19 lockdown on air quality trends across the UK. Sci. Total Environ. 2021, 754, 142374. [Google Scholar] [CrossRef]
  24. Bera, B.; Bhattacharjee, S.; Shit, P.K.; Sengupta, N.; Saha, S. Significant impacts of COVID-19 lockdown on urban air pollution in Kolkata (India) and amelioration of environmental health. Environ. Dev. Sustain. 2020, 23, 6913–6940. [Google Scholar] [CrossRef] [PubMed]
  25. Granella, F.; Reis, L.A.; Bosetti, V.; Tavoni, M. COVID-19 lockdown only partially alleviates health impacts of air pollution in Northern Italy. Environ. Res. Lett. 2021, 16, 035012. [Google Scholar] [CrossRef]
  26. Tadano, Y.S.; Potgieter-Vermaak, S.; Kachba, Y.R.; Chiroli, D.M.G.; Casacio, L.; Santos-Silva, J.C.; Moreira, C.A.B.; Machado, V.; Alves, T.A.; Siqueira, H.; et al. Dynamic model to predict the association between air quality, COVID-19 cases, and level of lockdown. Environ. Pollut. 2021, 268, 115920. [Google Scholar] [CrossRef] [PubMed]
  27. Chau, P.N.; Zalakeviciute, R.; Thomas, I.; Rybarczyk, Y. Deep Learning Approach for Assessing Air Quality During COVID-19 Lockdown in Quito. Front. Big Data 2022, 5, 1–13. [Google Scholar] [CrossRef]
  28. Grange, S.K.; Carslaw, D.C.; Lewis, A.C.; Boleti, E.; Hueglin, C. Random forest meteorological normalisation models for Swiss PM10 trend analysis. Atmos. Chem. Phys. 2018, 18, 6223–6239. [Google Scholar] [CrossRef] [Green Version]
  29. Grange, S.K.; Carslaw, D.C. Using meteorological normalisation to detect interventions in air quality time series. Sci. Total Environ. 2019, 653, 578–588. [Google Scholar] [CrossRef]
  30. Šimić, I.; Lovrić, M.; Godec, R.; Kröll, M.; Bešlić, I. Applying machine learning methods to better understand, model and estimate mass concentrations of traffic-related pollutants at a typical street canyon. Environ. Pollut. 2020, 263, 114587. [Google Scholar] [CrossRef]
  31. Kuo, C.P.; Fu, J.S. Evaluating the impact of mobility on COVID-19 pandemic with machine learning hybrid predictions. Sci. Total Environ. 2021, 758, 144151. [Google Scholar] [CrossRef]
  32. Brancher, M. Increased ozone pollution alongside reduced nitrogen dioxide concentrations during Vienna’s first COVID-19 lockdown: Significance for air quality management. Environ. Pollut. 2021, 284, 117153. [Google Scholar] [CrossRef]
  33. Etchie, T.O.; Etchie, A.T.; Jauro, A.; Pinker, R.T.; Swaminathan, N. Season, not lockdown, improved air quality during COVID-19 State of Emergency in Nigeria. Sci. Total Environ. 2021, 768, 145187. [Google Scholar] [CrossRef]
  34. Xu, Z.; Cao, R.; Hu, X.; Han, W.; Wang, Y.; Huang, J.; Li, G. The improvement of air quality and associated mortality during the COVID-19 lockdown in one megacity of china: An empirical strategy. Int. J. Environ. Res. Public Health 2021, 18, 8702. [Google Scholar] [CrossRef] [PubMed]
  35. Gope, S.; Dawn, S.; Das, S.S. Effect of COVID-19 pandemic on air quality: A study based on Air Quality Index. Environ. Sci. Pollut. Res. 2021, 28, 35564–35583. [Google Scholar] [CrossRef] [PubMed]
  36. Jakovljević, I.; Pehnec, G.; Vađić, V.; Čačković, M.; Tomašić, V.; Jelinić, J.D. Polycyclic aromatic hydrocarbons in PM10, PM2.5 and PM1 particle fractions in an urban area. Air Qual. Atmos. Health 2018, 11, 843–854. [Google Scholar] [CrossRef]
  37. Pehnec, G.; Jakovljević, I. Carcinogenic Potency of Airborne Polycyclic Aromatic Hydrocarbons in Relation to the Particle Fraction Size. Int. J. Environ. Res. Public Health 2018, 15, 2485. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  38. Godec, R.; Jakovljević, I.; Šega, K.; Čačković, M.; Bešlić, I.; Davila, S.; Pehnec, G. Carbon species in PM10 particle fraction at different monitoring sites. Environ. Pollut. 2016, 216, 700–710. [Google Scholar] [CrossRef] [PubMed]
  39. Air Quality in the Republic of Croatia. Available online: http://iszz.azo.hr/iskzl/ (accessed on 30 May 2022).
  40. WHO. WHO Global Air Quality Guidelines: Particulate Matter (PM2.5 and PM10), Ozone, Nitrogen Dioxide, Sulfur Dioxide and Carbon Monoxide; World Health Organization: Geneva, Switzerland, 2021; ISBN 9789240034228.
  41. Šimić, I.; Lovrić, M.; Bešlić, I.; Pehnec, G.; Godec, R. Particulate matter concentrations (PM1, PM2.5, PM10) since 2009 for a measurement sites in Zagreb, Croatia. Zenodo 2022. [Google Scholar] [CrossRef]
  42. Breiman, L. Random Forests. Mach. Learn. 2001, 45, 5–32. [Google Scholar] [CrossRef] [Green Version]
  43. Ke, G.; Meng, Q.; Finley, T.; Wang, T.; Chen, W.; Ma, W.; Ye, Q.; Liu, T.Y. LightGBM: A highly efficient gradient boosting decision tree. Adv. Neural Inf. Process. Syst. 2017, 2017, 3147–3155. [Google Scholar]
  44. Zhong, J.; Zhang, X.; Gui, K.; Wang, Y.; Che, H.; Shen, X.; Zhang, L.; Zhang, Y.; Sun, J.; Zhang, W. Robust prediction of hourly PM2.5 from meteorological data using LightGBM. Natl. Sci. Rev. 2021, 8, nwaa307. [Google Scholar] [CrossRef]
  45. Lovrić, M.; Meister, R.; Steck, T.; Fadljević, L.; Gerdenitsch, J.; Schuster, S.; Schiefermüller, L.; Lindstaedt, S.; Kern, R. Parasitic resistance as a predictor of faulty anodes in electro galvanizing: A comparison of machine learning, physical and hybrid models. Adv. Model. Simul. Eng. Sci. 2020, 7, 1–16. [Google Scholar] [CrossRef]
  46. Lovrić, M.; Pavlović, K.; Žuvela, P.; Spataru, A.; Lučić, B.; Kern, R.; Wong, M.W. Machine learning in prediction of intrinsic aqueous solubility of drug-like compounds: Generalization, complexity, or predictive ability? J. Chemom. 2021, 35, e3349. [Google Scholar] [CrossRef]
  47. Snoek, J.; Larochelle, H.; Adams, R.P. Practical Bayesian Optimization of Machine Learning Algorithms. In Proceedings of the NIPS 2012, Lake Tahoe, NV, USA, 3–8 December 2012; Volume 4, pp. 2951–2959. [Google Scholar]
  48. About, A. What’s the Government Response Timeline for Croatia? | Ask About. Available online: https://askabout.io/covid-19/ask/what-is-the-government-response-timeline-for-croatia/ (accessed on 30 May 2022).
  49. IUS-INFO: Decisions of the Civil Protection Headquarters of the Republic of Croatia in 2020 and Relevant Regulations in the Conditions of the Coronavirus Epidemic. Available online: https://www.iusinfo.hr/aktualno/u-sredistu/41376#travanj2020 (accessed on 30 May 2022).
  50. Government of the Republic of Croatia—Coronavirus Protection Measures. Available online: https://vlada.gov.hr/coronavirus-protection-measures/28950 (accessed on 30 May 2022).
  51. Stipaničev, D.; Repec, S.; Vucić, M.; Lovrić, M.; Klobučar, G. COVID-19 Lockdowns–Effect on Concentration of Pharmaceuticals and Illicit Drugs in Two Major Croatian Rivers. Toxics 2022, 10, 241. [Google Scholar] [CrossRef] [PubMed]
  52. Federal Office: MeteoSwiss Saharan dust events—MeteoSwiss. Available online: https://www.meteoswiss.admin.ch/home/climate/the-climate-of-switzerland/specialties-of-the-swiss-climate/saharan-dust-events.html (accessed on 31 July 2020).
  53. Davidović, M.; Dmitrašinović, S.; Jovanović, M.; Radonić, J.; Jovašević-stojanović, M. Diurnal, temporal and spatial variations of main air pollutants before and during emergency lockdown in the city of novi sad (Serbia). Appl. Sci. 2021, 11, 1212. [Google Scholar] [CrossRef]
  54. Zoran, M.A.; Savastru, R.S.; Savastru, D.M.; Tautan, M.N. Assessing the relationship between surface levels of PM2.5 and PM10 particulate matter impact on COVID-19 in Milan, Italy. Sci. Total Environ. 2020, 738, 139825. [Google Scholar] [CrossRef] [PubMed]
  55. Collivignarelli, M.C.; Abbà, A.; Bertanza, G.; Pedrazzani, R.; Ricciardi, P.; Carnevale Miino, M. Lockdown for COVID-2019 in Milan: What are the effects on air quality? Sci. Total Environ. 2020, 732, 139280. [Google Scholar] [CrossRef] [PubMed]
  56. Eleftheriadis, K.; Gini, M.I.; Diapouli, E.; Vratolis, S.; Vasilatou, V.; Fetfatzis, P.; Manousakas, M.I. Aerosol microphysics and chemistry reveal the COVID19 lockdown impact on urban air quality. Sci. Rep. 2021, 11, 1–12. [Google Scholar] [CrossRef]
  57. Gorrochategui, E.; Hernandez, I.; Pérez-gabucio, E. Temporal Air Quality (NO2, O3 and PM10) Changes in Urban and Rural Stations in Catalonia during COVID-19 Lockdown: An Association with Human Mobility and Satellite Data. Environ. Sci. Pollut. Res. 2022, 29, 18905–18922. [Google Scholar] [CrossRef]
  58. Tobías, A.; Carnerero, C.; Reche, C.; Massagué, J.; Via, M.; Minguillón, M.C.; Alastuey, A.; Querol, X. Changes in air quality during the lockdown in Barcelona (Spain) one month into the SARS-CoV-2 epidemic. Sci. Total Environ. 2020, 726, 138540. [Google Scholar] [CrossRef]
  59. Higham, J.E.; Ramírez, C.A.; Green, M.A.; Morse, A.P. UK COVID-19 lockdown: 100 days of air pollution reduction? Air Qual. Atmos. Health 2021, 14, 325–332. [Google Scholar] [CrossRef]
  60. Orak, N.H.; Ozdemir, O. The impacts of COVID-19 lockdown on PM10 and SO2 concentrations and association with human mobility across Turkey. Environ. Res. 2021, 197, 111018. [Google Scholar] [CrossRef]
  61. Hashim, B.M.; Al-Naseri, S.K.; Al-Maliki, A.; Al-Ansari, N. Impact of COVID-19 lockdown on NO2, O3, PM2.5 and PM10 concentrations and assessing air quality changes in Baghdad, Iraq. Sci. Total Environ. 2021, 754, 141978. [Google Scholar] [CrossRef] [PubMed]
  62. Al-Hemoud, A.; Al-Khayat, A.; Al-Dashti, H.; Li, J.; Alahmad, B.; Koutrakis, P. PM2.5 and PM10 during COVID-19 lockdown in Kuwait: Mixed effect of dust and meteorological covariates. Environ. Challenges 2021, 5, 100215. [Google Scholar] [CrossRef]
  63. Singh, V.; Singh, S.; Biswal, A.; Kesarkar, A.P.; Mor, S.; Ravindra, K. Diurnal and temporal changes in air pollution during COVID-19 strict lockdown over different regions of India. Environ. Pollut. 2020, 266, 115368. [Google Scholar] [CrossRef] [PubMed]
  64. Meji, M.A.; Dennison, M.S.; Mobisha, M.; Umar, M.M. Impact of COVID-19 Induced Lockdown on the Air Quality in the Southern Key Regions of India. i-Manager’s J. Futur. Eng. Technol. 2021, 16, 11. [Google Scholar] [CrossRef]
  65. Nakada, L.Y.K.; Urban, R.C. COVID-19 pandemic: Impacts on the air quality during the partial lockdown in São Paulo state, Brazil. Sci. Total Environ. 2020, 730, 139087. [Google Scholar] [CrossRef]
  66. Sicard, P.; De Marco, A.; Agathokleous, E.; Feng, Z.; Xu, X.; Paoletti, E.; Rodriguez, J.J.D.; Calatayud, V. Amplified ozone pollution in cities during the COVID-19 lockdown. Sci. Total Environ. 2020, 735, 139542. [Google Scholar] [CrossRef]
  67. Zhang, Z.; Arshad, A.; Zhang, C.; Hussain, S.; Li, W. Unprecedented temporary reduction in global air pollution associated with COVID-19 forced confinement: A continental and city scale analysis. Remote Sens. 2020, 12, 2420. [Google Scholar] [CrossRef]
  68. Rodríguez-Urrego, D.; Rodríguez-Urrego, L. Air quality during the COVID-19: PM2.5 analysis in the 50 most polluted capital cities in the world. Environ. Pollut. 2020, 266, 115042. [Google Scholar] [CrossRef]
  69. Venter, Z.S.; Aunan, K.; Chowdhury, S.; Lelieveld, J. COVID-19 lockdowns cause global air pollution declines. Proc. Natl. Acad. Sci. USA 2020, 117, 18984–18990. [Google Scholar] [CrossRef]
  70. Fu, F.; Purvis-Roberts, K.L.; Williams, B. Impact of the COVID-19 pandemic lockdown on air pollution in 20 major cities around the world. Atmosphere 2020, 11, 1189. [Google Scholar] [CrossRef]
  71. Chauhan, A.; Singh, R.P. Decline in PM2.5 concentrations over major cities around the world associated with COVID-19. Environ. Res. 2020, 187, 109634. [Google Scholar] [CrossRef] [PubMed]
  72. Connerton, P.; de Assunção, J.V.; de Miranda, R.M.; Slovic, A.D.; Pérez-Martínez, P.J.; Ribeiro, H. Air quality during COVID-19 in four megacities: Lessons and challenges for public health. Int. J. Environ. Res. Public Health 2020, 17, 5067. [Google Scholar] [CrossRef] [PubMed]
  73. Ceballos-Santos, S.; González-Pardo, J.; Carslaw, D.C.; Santurtún, A.; Santibáñez, M.; Fernández-Olmo, I. Meteorological Normalisation Using Boosted Regression Trees to Estimate the Impact of COVID-19 Restrictions on Air Quality Levels. Int. J. Environ. Res. Public Health 2021, 18, 13347. [Google Scholar] [CrossRef] [PubMed]
  74. Hoffer, J.G.; Ofner, A.B.; Rohrhofer, F.M.; Lovric, M.; Kern, R.; Lindstaedt, S.; Geiger, B.C. Theory-inspired machine learning—Towards a synergy between knowledge and data. Weld. World 2022, 66, 1291–1304. [Google Scholar] [CrossRef]
Ijerph 19 06937 g001 550
Figure 1. A time series plot of the collected particle mass concentration data from 2009–2020 for PM1, PM2.5 and PM10. For the sake of simplicity, the data is plotted as rolling monthly averages.
Figure 1. A time series plot of the collected particle mass concentration data from 2009–2020 for PM1, PM2.5 and PM10. For the sake of simplicity, the data is plotted as rolling monthly averages.
Ijerph 19 06937 g001
Ijerph 19 06937 g002 550
Figure 2. Datasets for 2020: model validation set-MVS (3 January–15 March); comparison set-CS (January and February); official lockdown (13 March–11 May); lockdown set- LDS (1–30 April); new normal set-NNS (1 June–31 July); construction works (March).
Figure 2. Datasets for 2020: model validation set-MVS (3 January–15 March); comparison set-CS (January and February); official lockdown (13 March–11 May); lockdown set- LDS (1–30 April); new normal set-NNS (1 June–31 July); construction works (March).
Ijerph 19 06937 g002
Ijerph 19 06937 g003 550
Figure 3. Out-of-ordinary events during 2020 shown as ratio between: (a) PM10/PM2.5 and (b) PM10/PM1. The events in March 2020 are assigned to either the Zagreb-earthquake on 22 March, the sand dust event between 24 and 30 March, or construction works at the site in early March, August, and September. These events were excluded from the analyzed timeframes.
Figure 3. Out-of-ordinary events during 2020 shown as ratio between: (a) PM10/PM2.5 and (b) PM10/PM1. The events in March 2020 are assigned to either the Zagreb-earthquake on 22 March, the sand dust event between 24 and 30 March, or construction works at the site in early March, August, and September. These events were excluded from the analyzed timeframes.
Ijerph 19 06937 g003
Ijerph 19 06937 g004 550
Figure 4. An overview of temperature and precipitation in Zagreb, Croatia through the studied timeframe 2009–2020. The data is plotted as a 3-month average and given trend-line by means of a regression line. The grey horizontal lines are the minimal values in of the regression line, showing the deviation from the minimal value.
Figure 4. An overview of temperature and precipitation in Zagreb, Croatia through the studied timeframe 2009–2020. The data is plotted as a 3-month average and given trend-line by means of a regression line. The grey horizontal lines are the minimal values in of the regression line, showing the deviation from the minimal value.
Ijerph 19 06937 g004
Ijerph 19 06937 g005 550
Figure 5. Schematics of the modelling procedure presented in this work.
Figure 5. Schematics of the modelling procedure presented in this work.
Ijerph 19 06937 g005
Ijerph 19 06937 g006 550
Figure 6. (a) Normalized medians for January and February (CS) during 2009–2020, (b) Normalized medians for June and July (NNS) during 2009–2020. The green line is PM10 related to the respective right axis, while blue and orange are PM1 and PM2.5 mass concentrations respectively referred to the left axis. The horizontal lines present the lowest values for the axes respectively.
Figure 6. (a) Normalized medians for January and February (CS) during 2009–2020, (b) Normalized medians for June and July (NNS) during 2009–2020. The green line is PM10 related to the respective right axis, while blue and orange are PM1 and PM2.5 mass concentrations respectively referred to the left axis. The horizontal lines present the lowest values for the axes respectively.
Ijerph 19 06937 g006
Ijerph 19 06937 g007 550
Figure 7. Boxplots of April normalized concentrations through the years for (a) PM1, (b) PM2.5, (c) PM10. The boxes show the quartiles of the dataset while the whiskers show the rest of the distribution, except for points that are determined to be “outliers” (diamond-shaped markers) by being outside of the inter-quartile range.
Figure 7. Boxplots of April normalized concentrations through the years for (a) PM1, (b) PM2.5, (c) PM10. The boxes show the quartiles of the dataset while the whiskers show the rest of the distribution, except for points that are determined to be “outliers” (diamond-shaped markers) by being outside of the inter-quartile range.
Ijerph 19 06937 g007
Table
Table 1. Results of the machine learning models for PM-concentrations shown on the validation set (MVS).
Table 1. Results of the machine learning models for PM-concentrations shown on the validation set (MVS).
PollutantRMSER2 ScoreWinning Algorithm
PM1010.470.77Random Forests
PM2.59.870.78Random Forests
PM16.490.77LightGBM
Table
Table 2. Literature findings on the lockdown’s effects on particulate matter concentration. Abbreviations used: Machine learning (ML), Descriptive statistics (DS), Modelling (MD), Unsupervised methods (UM), Meteorology (Met), Machine learning with normalization (MLN).
Table 2. Literature findings on the lockdown’s effects on particulate matter concentration. Abbreviations used: Machine learning (ML), Descriptive statistics (DS), Modelling (MD), Unsupervised methods (UM), Meteorology (Met), Machine learning with normalization (MLN).
Geographic
Location		Pollutants MethodsData UsedRef.
Zagreb, CroatiaPM10, PM2.5, PM1.0MLNTraining: from 1 January 2019 to 31 December 2019 (114 samples)
Validation: 3 smaller datasets in 2020 (10 samples); Test: 4 May to 13 May 12020 (10 samples)		This study
Zagreb, CroatiaNO2, PM10DSComparison between lockdown period (26 February–7 May 2020) and the same period in 2019[17]
Zagreb, CroatiaNO2, PM1.0, PAHs DSComparison between lockdown period (March–May 2020) and the same period in 2019[18]
Novi Sad, SerbiaPM2.5, NO2, NO, NOx, CO, SO2 + MetDSComparison before and after entering the state of emergency (1 February to 30 April)[53]
Skopje, Bitola, Tetovo, Kumanovo, MacedoniaPM10, PM2.5, NO2, O3, CO, MetDS Comparison of COVID19 period (last week of February 2020 to the end of May 2020) with the same period in 2017–2019 (nonCOVID-19 period)[19]
Milan, ItalyPM10, PM2.5, O3, NO2, SO2, CO, air quality index (AQI) + MetDSComparison between pre-lockdown (January–February 2020) and lockdown period (March–April 2020)[54]
Milan, ItalyPM10, PM2.5, BC, benzene, CO, NO2, O3, NOx + MetDSComparison between periods: CTRL (from 7 February 2020 to February 20), PL (from 9 March 2020 to 22 March 2020), and TL (from 23 March 2020 to 5 April 2020)[55]
Milan, Bologna, Florence, Rome, Naples, and Palermo, ItalyPM10, PM2.5, NO2, O3 + MetDSComparison between 2019-period (25 February–2 May 2019) and 2020-period (24 February–30 April 2020)[20]
Athens, GreecePM2.5, PM1.0, eBC, EC, OC, paricle number size distribution, SO42-, NO3-, Cl-, NH+ + MetDSComparison of reference period (1 January–10 March 2020) the two lockdown periods (11 March–22 March 2020 & 23 March–12 April 2020) with the respective periods in 2018 and 2019 [56]
Barcelona & Catalonia, SpainNO2, O3, PM10—hourly samplesDSComparisons during the before (15 February to 13 March), during (14 March to 21 June) and after lockdown (22 June to 31 August)[57]
Barcelona, SpainPM10, NO2, SO2, O3, BC + MetDSComparison for the periods before (16 February to 13 March) and during the lockdown (14 March to 30 March)[58]
Madrid, Barcelona, SpainNO2—hourly samples + MetDSComparison of March in the years 2018, 2019 and 2020[21]
South East of the UKNO2, PM2.5, PM10, O3 + MetDSComparison between lockdown period (March–May 2020) with the same period in 2015–2019 [22]
UKNO, NO2, NOx, O3, PM10, PM2.5—hourly samplesDSComparison between lockdown period (1 January to 30 June 2020) with the period from 1 January 2015 to 31 December 2019[23]
London, Glasgow, Belfast, Birmingham, Manchester and Liverpool, UKNOx, SO2, PM2.5, O3 + MetDSComparison of 100 days post-lockdown (23 March to 30 June 2020) with the same period from the previous 7 years[59]
TurkeyPM10, SO2, DSComparison of 2020 to the average of the 5-year period (2015–2019)[60]
Baghdad, IraqNO2, O3, PM2.5, PM10, AQIDSComparison of the periods before the lockdown from 16 January to 29 February 2020, and during four periods of partial and total lockdown from (1 March to 24 July 2020)[61]
KuwaitPM10, PM2.5 + MetDSComparison between the lockdown in 2020 with the corresponding periods of the years 2017–2019[62]
IndiaPM2.5, PM10, NO2, O3, CO, SO2 + Met—hourlyDSComparison between lockdown period (25 March–3 May 2020) and the same period in 2017–2019[63]
Southern regions of IndiaPM2.5, PM10, NO, CO, O3DSComparison between lockdown period (1 April to 31 July 2020) and the same periods in 2018 and 2019[64]
Kolkata City, IndiaPM10, PM2.5, O3, SO2, NO2, COUMComparison of lockdown period (25 March to 15 May 2020), with the similar time frame in 2017, 2018 and 2019[24]
Sao Paulo, BrazilNO, NO2, CO, PM2.5, PM10, SO2, O3, NOx DSComparison the partial lockdown periods (25 February 2020 to 23 March 2020 and from 24 March 2020 to 20 April 2020) to the five-year monthly trend (February, March and April of the years 2015, 2016, 2017, 2018 and 2019)[65]
Nice (France), Rome and Turin (Italy), Valencia (Spain) and Wuhan (China)NOx, PM2.5, PM10, O3DSComparison of lock down period (1 January 2017 until 18 April 2020) with the same period over the three previous years (2017–2019)[66]
sixteen selected cities located in South Asia, East Asia, Europe, and North AmericaNO2, CO, PM2.5, O3, SO2DSComparison between from 1 January–15 May for the year of 2015–2019 (defined as baseline period) and 2020 (lockdown)[67]
50 most polluted capital cities in the worldPM2.5, AQIDSComparison between before and during quarantine [68]
34 countriesNO2, O3, PM2.5DSComparison between from 1 January–15 May for the year of 2017–2019 and 2020 (lockdown)[69]
Multiple locations *NO2, SO2, CO, O3, PM10, PM2.5, AQIDSComparison between lockdown period in 2020 to the same period of 2017, 2018 and 2019[70]
New York, Los Angeles, Zaragoza, Rome, Dubai, Delhi, Mumbai, Beijing and ShanghaiPM2.5DSComparison of lockdown period (December 2019–March 2020), and the same period in earlier years 2017–2019[71]
São Paulo in Brazil; Paris in France; and Los Angeles and New York in the USANO2, CO, PM2.5, O3 + meteorologyDSComparison of March in the years 2015–2020[72]
Graz, AustriaNO2, PM10, O3, Ox + MetMLTraining: from 3 January 2014 to 31 December 2019 (daily)
Validation: from 3 January 2020 to 10 March 2020 (daily),
Test:l lockdown set, LD (10 March 2020 to 2 May 2020—daily samples), and a hard lockdown set, HLD (20 March 2020 to 14 April 2020—daily samples) 		[13]
Lombardy, ItalyNO2, PM2.5 + MetMLTraining: from 2012 through 2019
Validation: months from January to April for 2016–2020
Test: from January through early May 2020		[25]
Sao Paulo, BrazilCO, O3, NO2, NO, PM2.5, PM10 + MetML Training: from 1 January to 23 April 2020 (114 samples); Validation: 24 April to 3 May 2020 (10 samples); Test: 4 May to 13 May 2020 (10 samples)[26]
Quito, EcuadorCO, NO2, PM2.5, SO2, O3MLN Training: from 1 January 2016 to 15 January 2020 (2 months before the COVID-19 lockdown)
Test: from 16 January 2020 to 15 March 2020 (the day of the national lockdown).		[27]
Cantabria, SpainNO, NO2, PM10, O3, MetMLNData from 11 stations (2013–2020)
Training data 2013–2019, test set lockdown and new normal 2020		[73]
Vienna, Austria
* Wuhan, Beijing (China), Delhi (India), Tehran (Iran), Istanbul (Turkey) in Asia; Rome (Italy), Madrid (Spain), Paris (France), London (UK), Berlin (Germany) and Moscow (Russia) in Europe; Johannesburg (South Africa) in Africa and Los Angeles, New York City (USA), Mexico city (Mexico), Sao Paulo (Brazil) and Lima (Peru) in North and South America.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Lovrić, M.; Antunović, M.; Šunić, I.; Vuković, M.; Kecorius, S.; Kröll, M.; Bešlić, I.; Godec, R.; Pehnec, G.; Geiger, B.C.; et al. Machine Learning and Meteorological Normalization for Assessment of Particulate Matter Changes during the COVID-19 Lockdown in Zagreb, Croatia. Int. J. Environ. Res. Public Health 2022, 19, 6937. https://doi.org/10.3390/ijerph19116937

AMA Style

Lovrić M, Antunović M, Šunić I, Vuković M, Kecorius S, Kröll M, Bešlić I, Godec R, Pehnec G, Geiger BC, et al. Machine Learning and Meteorological Normalization for Assessment of Particulate Matter Changes during the COVID-19 Lockdown in Zagreb, Croatia. International Journal of Environmental Research and Public Health. 2022; 19(11):6937. https://doi.org/10.3390/ijerph19116937

Chicago/Turabian Style

Lovrić, Mario, Mario Antunović, Iva Šunić, Matej Vuković, Simonas Kecorius, Mark Kröll, Ivan Bešlić, Ranka Godec, Gordana Pehnec, Bernhard C. Geiger, and et al. 2022. "Machine Learning and Meteorological Normalization for Assessment of Particulate Matter Changes during the COVID-19 Lockdown in Zagreb, Croatia" International Journal of Environmental Research and Public Health 19, no. 11: 6937. https://doi.org/10.3390/ijerph19116937

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop