Gradient Boosting Machine to Assess the Public Protest Impact on Urban Air Quality

: Political and economic protests build-up due to the ﬁnancial uncertainty and inequality spreading throughout the world. In 2019, Latin America took the main stage in a wave of protests. While the social side of protests is widely explored, the focus of this study is the evolution of gaseous urban air pollutants during and after one of these events. Changes in concentrations of NO 2 , CO, O 3 and SO 2 during and after the strike, were studied in Quito, Ecuador using two approaches: (i) inter-period observational analysis; and (ii) machine learning (ML) gradient boosting machine (GBM) developed business-as-usual (BAU) comparison to the observations. During the strike, both methods showed a large reduction in the concentrations of NO 2 (31.5–32.36%) and CO (15.55–19.85%) and a slight reduction for O 3 and SO 2 . The GBM approach showed an exclusive potential, especially for a lengthier period of predictions, to estimate strike impact on air quality even after the strike was over. This advocates for the use of machine learning techniques to estimate an extended effect of changes in human activities on urban gaseous pollution.


Introduction
Polluted air can contain a substantial amount of chemicals in the form of microparticles and harmful gases. Among the most common atmospheric gaseous pollutants are carbon monoxide (CO), nitrogen oxides (NO and NO 2 ), sulfur dioxide (SO 2 ) and ozone (O 3 ). All of them can cause a variety of respiratory and cardiovascular health problems and are especially dangerous for at-risk populations: the young, the elderly and people with respiratory or cardiovascular problems [1]. The air pollutant concentrations can be affected by social dynamics, such as the trigger of mass protests due to economic troubles, job insecurity and rapidly growing inequality [2]. In 2019, a huge number of public protests and strikes spiked up in developing and developed countries, including Europe, Asia and the Americas [3][4][5][6]. At the end of 2019, Latin America was in the center stage of demonstrations against governments in Mexico [7], Chile [8], Bolivia [9,10], Colombia [11] and Ecuador [12]. Massive protests and riots spread from capital cities overtaking significant parts of the countries, and in a few cases even lasting well into 2020 or 2021 [13,14].
As social dissatisfaction is growing more vocal, it encourages scientific curiosity to question the effects that protests and strikes might have on air quality. During those special events, the protesters often cause restrictions or complete seizure of public and private 2 of 15 transportation, which might improve air quality, however, they might also cause riots, vandalisms, burning and looting [15]. Though protests have previously been investigated from the social science perspective [16,17], the research on air pollution dynamics is still very limited. Previous research analyzing the air quality during strike events has demonstrated a variety of effects on air quality. A study in Hong Kong during large-scale street protests reported a substantial reduction in NO 2 in the central area [5]. An investigation in São Paulo during a truck drivers' strike showed a 50% reduction in CO and NO concentrations, although NO 2 and particulate matter (PM) indicated no changes [18]. In contrast, during public transport strikes in Italy and Spain, CO levels increased due to the greater use of private vehicles [19,20]. Such differences in research findings suggest a need for a broader understanding of the impact of these events on air quality.
In previous studies, pollution change estimates are often based on differences between central and background sites [5] or "trend conditions". The latter are usually calculated using a comparable period (e.g., 3 to 5 years) prior to the protest event [18] or simply comparing pollution pre-strike levels to those during the strike [21]. Depending on the nature of the pollutant, these types of analyses could be subject to biases, as meteorological conditions that influence the production and transport of certain contaminants could vary between the time periods used for comparison [22]. In the case of traffic-borne primary pollutants, however, macroscopic observations before and after events that cause a change in traffic, are useful as differences are directly related to the physical source of pollution. In the case of statistics-based models (e.g., machine learning algorithms), meteorology variables can be included in order to avoid biases and to increase confidence when calculating business-as-usual (BAU) conditions, as previously done in COVID-19 impact studies [23][24][25].
As the number and scale of protests spiral, especially considering the economic burden due to the global COVID-19 pandemic, it is crucial to be able to evaluate the true environmental impact of these events. Hence, the purpose of this study is to find out how local is the impact of these events on urban air quality and how long do the effects last after the protests stop? For that, a new approach is proposed to quantify the effect of such an event on urban pollution. A specific South American case study in Quito, Ecuador, is investigated to gauge changes in air quality during the national strike of 2019 and up to 12 days later. Due to the location of the central government in Quito, protesters, mostly indigenous groups, traveled hundreds of kilometers by foot to Quito from all over the country. Throughout the protests, there was a voluntary (private bus companies) and forced (city transportation lines and private vehicles) lockdown of urban transport circulation. Otherwise, circulating vehicles were often attacked and forced to stop by protesters. These unusual conditions allowed for estimating the impact of the motorized fleet on trace gas levels in this high elevation Andean capital.
In this study, to estimate the protest impact on urban air quality in the South American city, first, the strike impact on concentrations of NO 2 , CO, O 3 and SO 2 was evaluated using an inter-period observational data comparison, considering 12 days before, 12 days during and 12 days after the strike. Then, the evolution of trace gas levels was assessed using machine learning (ML) gradient boosting machine (GBM) modeled business-as-usual (BAU) conditions. The comparison between the two approaches was investigated. Finally, the two model performance evaluations conclude the study.

Study Sites, Data Analysis and Visualization
Geographically divided by the Andes mountains and the equator, Ecuador (17.3 million inhabitants) is the most densely populated South American country (71 people km −2 ) [26,27]. Due to a wide range of elevations influencing ecosystems, this region is known for its extremely high biodiversity. The Ecuadorian capital of Quito is positioned at an average elevation of 2850 m above sea level (m.a.s.l.) on the western branch of the Andes cordillera ( Figure 1). It spreads over an assortment of heights, varying from the slopes of an active Pichincha volcano (4800 m.a.s.l.) towards the inter-Andean valleys (2300 m.a.s.l.). The change in elevation in the capital city is responsible for an array of precipitation-dependent microclimates [28]. Close to 3,000,000 people inhabit the metropolitan district of Quito (DMQ) [29]. This urban center is characterized by a rapidly growing motorized fleet of vehicles powered by poor-quality fossil fuels affecting air quality [30,31].
Appl. Sci. 2021, 11, 12083 3 of 15 age elevation of 2850 m above sea level (m.a.s.l.) on the western branch of the Andes cordillera ( Figure 1). It spreads over an assortment of heights, varying from the slopes of an active Pichincha volcano (4800 m.a.s.l.) towards the inter-Andean valleys (2300 m.a.s.l.). The change in elevation in the capital city is responsible for an array of precipitation-dependent microclimates [28]. Close to 3,000,000 people inhabit the metropolitan district of Quito (DMQ) [29]. This urban center is characterized by a rapidly growing motorized fleet of vehicles powered by poor-quality fossil fuels affecting air quality [30,31].  In this study, air quality changes due to the political protests in Quito were measured at seven sites of the environmental monitoring network of the Secretariat of the Environment of DMQ. The study sites are well distributed across the city (north to south): S1-Carapungo (elev. 2660 m.a.s.l, coord. 78 • (Figure 1).
Air pollution and meteorological variables (Table A1, Appendix A) were sampled using the United States Environmental Protection Agency approved methods [32], described in detail in past publications [30]. Data from the ground-based sites were gathered from 20 September to 25 October 2019. This includes three 12-day periods: before (20 September 2019-1 October 2019), during (2 October 2019-13 October 2019) and after (14 October 2019-25 October 2019) the national strike of Ecuador. We believe that the S2-Cotocollao site had technical issues in measuring the NO 2 concentrations during the study period, thus it was excluded from the NO 2 analysis. The data for atmospheric variables were processed as 24-h averages and accumulation (rain only) for all seven study sites. In addition, 12-day statistics were processed for four main atmospheric gaseous pollutants: NO 2 , CO, O 3 and SO 2 . For each pollutant, the concentration change in each site was computed as the mean percent difference between the 12-days before the strike and values measured during (and after) the strike (e.g., % change = (pre-strike average − strike average)/(pre-strike average) × 100%). The pre-strike period might be presumed as business-as-usual conditions, due to the normal anthropogenic activity. Then, from those, the average values and standard deviations for each pollutant were calculated for the whole of Quito. The accuracy of the inter-period comparative approach was measured through the root mean squared error (RMSE), as follows: where "prediction" is the average concentration of the 12 days before the strike and "observation" stands for the hourly data on the same period of time. The performance of the inter-period comparison method is used as a benchmark to assess the reliability of the machine learning model to predict the pollution change. Air quality is influenced by environmental conditions. Thus, six meteorological parameters (relative humidity, precipitation, temperature, atmospheric pressure, wind speed and wind direction) were included in the analysis. Those were sampled from the same monitoring sites using complete Vaisala WXT536 instrumentation, while a Kipp&Zonnen (Delft, The Netherlands) pyrometer was used to measure global solar radiation. For hourly 24 h and 12-day data analysis and visualization, Microsoft Excel, Igor Pro 6.0 (Wavemetrics Inc., Nimbus, Portland) and QGIS 3.20.2 software were utilized. QGIS 3.20.2 software was utilized for spatial imaging of pollution, using the EPSG.4326 (WGS84) coordinate system (available at https://qgis.org/downloads/, accessed on 20 October 2021). Public access data was used for: country limits from Geographic Military Institute data (retrieved from https://sni.gob.ec/geoservicios-ecuador, accessed on 25 October 2021); for city limits urban growth of Quito (retrieved from http://gobiernoabierto.quito.gob.ec/ ?page_id=1122, accessed on 11 July 2021); and for the monitoring sites (retrieved from http://www.quitoambiente.gob.ec/index.php/generalidades, accessed on 20 October 2021). Roadblock and barricade sites were retrieved from Ecuadorian newscasts (available at https://www.primicias.ec/noticias/sociedad/quito-manifestaciones-vias/, https:// www.elcomercio.com/actualidad/bloqueos-escombros-personas-vias-quito.html, accessed on 11 July 2021) and Satellite XYZ Tiles service from Google (available at http://mt0.google. com/vt/lyrs=s&hl=en&x=\{x\}&y=\{y\}&z=\{z\}, accessed on 11 July 2021). Inverse distance weighting (distance coefficient p = 2; grid resolution 1 m) was used for city data interpolation between the study sites.

Machine Learning Modelling
The impact of the national protest on the concentrations of NO 2 , SO 2 , CO and O 3 was quantified through a machine learning (ML) approach. A model was trained to predict business-as-usual concentrations of the pollutants from hourly measurements of meteorological and temporal variables. Due to the nonlinear relationships between meteorology and air quality, we selected a nonparametric algorithm called gradient boosting machine (GBM). It is a state-of-the-art decision tree-based ensemble method (Friedman, 2001). As explained in our previous study on the strike effect on the concentration of fine particulate matter [15], we applied a meteorology-normalized technique [33] to get a prediction that considers the real climatic conditions during the unusual anthropogenic event.
The meteorology was analysed through seven features: solar radiation, temperature, pressure, relative humidity, precipitation, wind direction and wind speed. The temporal attributes were: hour, weekday, year day, and trend (date index). The dataset was split into a training set (1 January 2015-1 April 2019) and a testing set (2 April 2019-1 October 2019). Therefore, the performance of the models was tested six months prior to the protests. The accuracy of the models was measured through the coefficient of determination (R 2 ) and the root mean squared error (RMSE) (see Equation (1)). In the latter, "prediction" values are provided by the ML model outputs for the 12 days before the strike, while "observation" are observational data during the same period. This RMSE for the ML model will be compared to the RMSE for the inter-period method for each pollutant, in order to estimate the performance of the modeling approach in relation to the baseline (i.e., inter-period accuracy).
GBM builds an ensemble of weak and shallow consecutive trees with each tree learning and improving the previous one. The algorithm sequentially adds new models to the ensemble. At each iteration, a new weak base-learner model is trained from the error of the entire learned ensemble. A model is considered weak if its error rate is just slightly better than random guessing. In this boosting method, a weak model is built at each iterative sequence, in order to improve the remaining errors. The generalization of the algorithm is a stage-wise additive model of n individual regression trees as shown in Equation (2).
For each pollutant, the concentration change was computed from the mean difference between ML-based business-as-usual and the actual value measured during or after the strike (e.g., % change = (BAU − observed)/(BAU) × 100%). Then, from those, the average values and standard deviations for each pollutant were calculated for the whole of Quito. The GBM models were implemented in R, using the H 2 O package. The parameters of the library were tuned for stopping the training (i.e., increment in the number of trees) if ten successive trees did not reduce the prediction error.

Inter-Period Observational Approach
Ground-based observations are reported from seven well-distributed monitoring sites around the city. To better visualize the evolution of spatial observational data, 12-day average data inverse distance weighting images for pre-strike  Figure 2. Numerical analysis of percentage changes is also presented in Table 1. It can be seen that CO (−19.85 ± 6.67%) and NO 2 (−32.36 ± 15.38%) concentrations showed a much higher reduction during the strike when compared to the 12-day average before the strike (see Table 1 and Figure 2a,b,e,f). Meanwhile, O 3 and SO 2 showed much smaller decreases, −5.18 ± 7.32% and −7.20 ± 7.48%, respectively (see Table 1 and Figure 2c,d,g,h). For concentrations of SO 2 , there is a visible reduction in most of the sites (−7.2 ± 7.48%), except for the two sites showing an increase in SO 2 , one of them near the thermoelectric power plant in S7-Chillos (Figure 2h) showing no changes in electricity production. We note that the latter site reaches the highest SO 2 concentrations in the city, which is due to crude oil and poor-quality diesel fuel burned at the plant.
while, O3 and SO2 showed much smaller decreases, −5.18 ± 7.32% and −7.20 ± 7.48%, respectively (see Table 1 and Figure 2c,d,g,h). For concentrations of SO2, there is a visible reduction in most of the sites (−7.2 ± 7.48%), except for the two sites showing an increase in SO2, one of them near the thermoelectric power plant in S7-Chillos (Figure 2h) showing no changes in electricity production. We note that the latter site reaches the highest SO2 concentrations in the city, which is due to crude oil and poor-quality diesel fuel burned at the plant.    Similar results were observed in other studies carried out during strike events involving a suspension of motorized fleet circulation. A significant reduction in CO concentrations was observed in the city of Hyderabad in India. Explicitly, a CO level reduction of −40% was registered in 2009 [34] and 20% in 2010 [21]. In a study in Sao Paulo, Brazil, an even higher reduction, −67%, of CO levels was reported [18]. A reduction of NO 2 concentrations was also observed in Sao Paulo ranging from −7% [18] to −55.5% [35]. As the findings of our study, low reductions in SO 2 (−11%) were reported in another study in Sao Paulo [36]. On the other hand, during the truck strike event in 2009 in Hyderabad [34] a significant reduction in O 3 (−50%) was observed, while in Sao Paulo, the O 3 concentration reduction was lower (−16%) although still superior to our findings [18].

Inter-Period
As expected, after the strike ended, most of the gaseous criteria pollutants returned to the pre-strike levels (Figure 2i-l and Table 1). When NO 2 concentrations still showed a slight reduction (−2.79 ± 5.82%), CO and SO 2 increased a little by +1.62 ± 5.83% and +3.29 ± 13.56%, respectively (Table 1). These overall fluctuations, however, are quite negligible from the pre-strike scenario, thus suggesting a return of human activities to normal levels and a return of anthropogenic emissions to the urban air shed. Some hot spots in the city can be pinpointed. For example, for NO 2 concentrations, the central part of the city flairs back up post-strike (see Figure 2i), almost indistinguishably as before the strike (Figure 2a). This indicates a return of the traffic to the city streets and roads [37,38]. The northwestern part of the city stays less polluted, where most of the urban sprawl to residential, less dense areas is found. A similar effect can be observed for the combustion-related CO concentrations that returned to the industrial part of the city after the strike (Figure 2j).
Compared to the pre-strike levels (Figure 2c), the concentrations of O 3 decreased during (−5.18 ± 7.32%) and after the strike (−25.77 ± 12.01%) (Figure 2g,k). As important as meteorological conditions are to ozone formation and accumulation in the ambient air, a key question is whether the reduction in primary urban emissions, due to activities that were paralyzed by the strike, had an impact on the levels of ozone in Quito. We propose that a shift in emission proportions had an impact in terms of ozone production rates, mainly during those days when solar radiation was more available ( Figure A1, Appendix A). As demonstrated in previous work [39], during the strict COVID-19 lockdown, there was a substantial decrease in traffic emissions in Quito, and thus NO x , comparable to the strike (Figure 2e,f). As a result, the local atmospheric chemistry shifted from a generally NO x -saturated regime into conditions in which free radicals [OH (hydroxyl) and HO 2 (hydroperoxyl)]-the core of the ozone production chemistry-were less depleted to form subspecies such as nitric acid and, instead, produced additional ozone per unit time (ozone production rates, ppbv h −1 ). However, increased rates of ozone production were counterbalanced by meteorological conditions, as the strict lockdown also coincided with a period of increased cloudiness and rain. Such conditions had the effect of cleaning the atmosphere. Therefore, even though the strike did not last very long, it shares similar traits to the COVID-19 lockdown from the perspective that primary emissions decreased and meteorology was similar.
Finally, for concentrations of SO 2 , during the strike (−7.2 ± 7.48%) and post-strike (+3.29 ± 13.56%), the change is much smaller, with sustained peaks in the southern valley of the city, where a thermoelectric plant is located (Figure 2d,h,l). This indicates there were close to no changes in power generation during the strike. However, further away from the power plant, an impact of reduced city traffic circulation can be observed in almost all the sites during the strike.

Machine Learning (ML) Approach
As the straightforward observational data analysis of the inter-period comparison may be skewed by variation of meteorological conditions, it is important to estimate the pollution changes using a business-as-usual scenario modeled using the ML technique [25].  Table A2, Appendix A). However, when the national strike takes place, the observations, as expected, are clearly lower than the BAU concentrations ( Figure 3). Table 2 summarizes the all-site and overall Quito modeled percentage changes for each pollutant due to the strike activities using the machine learning approach. Some sites were more affected by the pollution changes than others, especially the city center, which was the epicenter of the protest. Overall, it can be seen that the highest reduction from the BAU scenario during the strike is for NO 2 (−31.53 ± 17.22%) ( Table 1). SO 2 (−17.65 ± 27.38%) and CO (−15.55 ± 6.64%) also show a significant reduction during the strike. These statistics confirm the importance of the suspension of human activities, especially the circulation of motorized traffic. This has been previously proven globally during the COVID-19 pandemic [25,[40][41][42][43]. While the primary pollutants show a significant reduction due to the reduction in anthropogenic activities during the national strike, O 3 shows only a slight reduction (−6.82 ± 15.11%) in its concentrations. The latter is due to a potential shift from a NO x -saturated regime into a more NO x -limited regime during the event of national transportation protest although more research is needed to better understand ozone.
Furthermore, even up to 10 days after the strike, the pollution levels took time to return to BAU conditions (Figure 3), varying from −21.49 ± 14.72% for O 3 to −12.43 ± 6.25% for NO 2 and −9.71 ± 32.25% for SO 2 as an overall reduction ( Table 2). This is an interesting observation, showing an extension of the effect of the violent strikes on anthropogenic activities even days after the event was over.  The color scale is based on a gradient color in red, yellow, and green. The greener is the cell, the higher is the pollution drop. Furthermore, even up to 10 days after the strike, the pollution levels took time to return to BAU conditions (Figure 3), varying from −21.49 ± 14.72% for O3 to −12.43 ± 6.25% for NO2 and −9.71 ± 32.25% for SO2 as an overall reduction ( Table 2). This is an interesting observation, showing an extension of the effect of the violent strikes on anthropogenic activities even days after the event was over.

Assessment of the Methods
While both methods show an overall reduction in the concentrations of most gaseous pollutants during the social strike in the Ecuadorian capital, some statistics vary. During the strike, both methods show a significant reduction in the concentrations of NO2 (interperiod: −32.36 ± 15.38% vs. ML model: −31.53 ± 17.22%) and CO (inter-period: −19.85 ± 6.67  The color scale is based on a gradient color in red, yellow, and green. The greener is the cell, the higher is the pollution drop.

Assessment of the Methods
While both methods show an overall reduction in the concentrations of most gaseous pollutants during the social strike in the Ecuadorian capital, some statistics vary. During the strike, both methods show a significant reduction in the concentrations of NO 2 (inter-period: −32.36 ± 15.38% vs. ML model: −31.53 ± 17.22%) and CO (inter-period: −19.85 ± 6.67 vs. ML model: −15.55 ± 6.64%). In both cases, the ML method is more conservative. Again, the two methods agree a smaller reduction is observed for O 3 (interperiod: −5.18 ± 7.32% vs. ML model: −6.82 ± 15.11%). These agreements suggest that both approaches are able to estimate equally well the impact of the strike on urban gaseous pollution. This might be due to the fact that the comparison of observations two weeks before the unusual event is within a short enough time to be as good as an ML modeling approach. Meanwhile, SO 2 concentrations show a much larger reduction using the ML model (−17.65 ± 27.38%) when compared to the inter-period analysis (−7.20 ± 7.48%).
After the strike, the only pollutant that both methods fully agree on is a continuous reduction of O 3 concentrations (inter-period: −25.77 ± 12.01% vs. ML model: −21.49 ± 14.72%). Again, the ML model is more conservative, suggesting a possible overestimation of the observational approach. The other three pollutants show a continuous reduction in concentrations using the ML model approach (NO 2 : −12.43 ± 6.25%; CO: −4.29 ± 6.50%; SO 2 : −9.71 ± 32.25%) and nearly no changes using the inter-period analysis (NO 2 : −2.79 ± 5.82%; CO: 1.62 ± 5.83%; SO 2 : 3.29 ± 13.56%). This again might propose the importance of using an ML BAU condition prediction to better estimate the extended effect of the protest on urban air quality.
In order to confirm which of the two approaches is more accurate, a method performance evaluation is presented in Table 3. It can be observed that the ML model, using the GBM method, demonstrates better performance in RMSE indicators in all cases except for a slightly worse performance on SO 2 . For this latter pollutant, the observational inter-period analysis tends to overperform the GBM (except the northern suburban S2-Cotocollao site), even if the difference is minuscule. This result is consistent with previous studies showing a lower performance of the models to predict SO 2 than other pollutants [25,33,44]. It suggests that the ML model should integrate additional features to improve the prediction of this particular contaminant [45]. For instance, gas-particle interaction must be taken into consideration in predicting this pollutant [46]. While the meteorology has been accounted for, the PM levels were not considered in the ML prediction for SO 2 . A previous study [15] showed large spikes in PM 2.5 concentrations during the strike (beyond the BAU conditions) in most of the sites. Some exceptions were observed in the northern sites (S2-Cotocollao and S3-Belisario). This might help explain a poorer ML performance to predict this pollutant, since the unusual PM behavior, due to barricade burning, is not taken into account. Table 3. Comparison of the performance between the two methods: gradient boosting machine (GBM) vs. inter-period analysis (I-P). The best accuracy (i.e., the lowest RMSE) is highlighted in yellow. Opposite to the SO 2 , the RMSE for NO 2 , CO and O 3 computed for the ML modeling approach are considerably smaller than the inter-period method. This outcome suggests that the machine learning method based on weather normalization tends to be more reliable than a simple inter-period comparison to estimate the effect of protests on urban gaseous pollution. However, we do not notice a large difference in the estimate of concentration change between the two methods (see Tables 1 and 2). This apparent contradiction can be explained by the fact that the protest occurred in quite a short period of time, for which the environmental condition (e.g., seasonal variation) did not change significantly. It suggests that a simple averaging approach can provide a fair estimate of pollution change on a few weeks of unchanged meteorology, but might have more complications to account for an event that lasts over several months, such as the COVID-19 lockdown [25].

NO
Ultimately, besides assessing the impact of a protest on urban air quality, the general purpose of this study was to propose a machine learning approach to quantify the effect of an unusual event on urban pollution. This was done by investigating the pollution concentrations during the business-as-usual conditions in an absence of a strike. There are different ways to calculate the BAU conditions. In recent studies focused on the COVID-19 lockdown impact on air pollution, the average values on the weeks [47] or the years [48] preceding the event are used as an approximation of the BAU. A limitation of the latter technique is the fact that it does not consider the meteorological conditions, which can significantly modify the assumed BAU pollution concentrations. In this study, we used the prior period averaging technique, but also provided a more robust estimation of the national protest on urban air pollution by developing a meteorology-normalized ML model [25,43] that takes into account the weather conditions during, and even after, the strike. The performance of the GBM method, compared to the averaging technique, is a more powerful method for predicting the true and lasting effect of an unusual event such as a protest on urban air quality.
This unique investigation offers an example from South America, where the protests often get violent and involve a strategy of burning toxic barricade objects, covering the city in smoke. However, we show that BAU emissions from traffic in this high elevation city are larger than those from protest activities. Therefore, this study can serve as a valuable lesson to help better estimate the future changes in pollution due to changes in human activities, such as pandemics, protests, or strikes anywhere on the planet. In addition, it can serve as a call for urban planning strategies in developing countries to steer the city transportation problems towards greener and more sustainable solutions.

Conclusions
This study successfully employs two different approaches to estimate the effect of a national strike of 2019 in Ecuador on concentrations of urban trace gases (NO 2 , CO, O 3 and SO 2 ): (i) observational inter-period analysis of air quality differences between the pre-strike period versus the periods during and post-strike in the capital city; and (ii) comparing a machine learning (ML) gradient boosting machine (GBM) modeled businessas-usual (BAU) conditions in Quito with the observed concentrations of atmospheric gaseous pollutants. Unexpectedly, the two approaches agree closely for urban NO 2 (ML: −31.5 ± %, inter-period: −32.36 ± 15.39%), CO (ML: −15.55 ± %, inter-period: −19.85 ± 6.67%) and O 3 (ML: −6.82 ± %, inter-period: −5.18 ± 7.32%) concentration changes during the strike. The ML approach overestimates the overall reduction of city SO 2 concentrations (ML: −17.65 ± %, inter-period: −7.20 ± 7.48%) during the national protest, which might show a limitation of this method on predicting SO 2 using selected meteorological parameters. We also estimate the effect of a national strike up to 12-days after the protest is over. While both approaches agree on a continuous decrease in O 3 pollution, suggesting the changes in precursor abundance even post-strike, they could not be more different for the other three gaseous pollutants. The observation-only inter-period approach underestimates the long-lasting effect on NO 2 reduction (ML: −12.43 ± 6.25%, inter-period: −2.79 ± 5.82%) and disagrees on the protest effects on CO (ML: −4.29 ± 6.50%, inter-period: 1.62 ± 5.83%) and SO 2 (ML: −9.71 ± 32.25%, inter-period: 3.29 ± 13.56%) after the protest ended.
We conclude that overall, the GBM model is more robust than the inter-period comparison method for predicting the true effect of an unusual event such as a protest on urban air quality. Nevertheless, some contaminants are more difficult to predict than others. This is the case of SO 2 , for which the prediction accuracy is not different than the baseline (i.e., the inter-period method). For this latter gas, better performance could be expected by considering additional features in the ML model (e.g., PM). This unique study presents an example from frequently violent protests, involving burning toxic barricade objects to block transport circulation, from South America. We show that the emissions from those burnings have less impact on atmospheric pollution at an urban scale than the footprint of "normal" city traffic. It can serve as a valuable tool to help better estimate the future changes in pollution due to changes in human activities, such as pandemics, protests or strikes.