Open Access Article
Lucille Joanna S.
Borlaza
a,
Vy Dinh
Ngoc Thuy
a,
Stuart
Grange
b,
Stéphane
Socquet
c,
Emmanuel
Moussu
c,
Gladys
Mary
c,
Olivier
Favez
d,
Christoph
Hueglin
b,
Jean-Luc
Jaffrezo
a and
Gaëlle
Uzu
*a
aUniversity of Grenoble Alpes, CNRS, IRD, INP-G, IGE (UMR 5001), Grenoble F-38000, France. E-mail: gaelle.uzu@ird.fr
bEmpa, Swiss Federal Laboratories for Materials Science and Technology, 8600 Dübendorf, Switzerland
cAtmo AuRA, Grenoble 38400, France
dINERIS, Parc Technologique Alata, BP 2, 60550 Verneuil-en-Halatte, France
First published on 11th April 2023
In Europe, COVID-19 lockdown restrictions were first implemented in March 2020 to control the spread of the disease from the SARS-CoV-2 virus. Many studies have focused on the influence of the applied measures on pollution levels during this period, but very limited information on the oxidative potential (OP), an emerging metric of particulate matter (PM) exposure. Furthermore, most previous studies also commonly used comparative methods with historical datasets, which may not be estimating the real pollution levels without the lockdown restrictions in place. In this study, the OP of PM collected at urban background (Grenoble, France) and traffic (Bern, Switzerland) sites was assessed using dithiothreitol (DTT) and ascorbic acid (AA) assays. These measurements were also compared with PM and black carbon (BC) mass concentrations, including the wood burning and fossil fuel fractions of BC. To obtain a more realistic pollution level, assuming there were no lockdown restrictions in place, a machine learning technique called the Random Forest (RF) regression model was applied to predict a business-as-usual (BAU) level for OP, PM, and BC in both sites. This model provided a good estimate of the BAU levels, allowing a more realistic assessment of the pollution changes during the lockdown period. The results indicate a clear decrease in OP found in the traffic site, while a more modest change in OP was found at the urban background site, likely due to sustained contributions from wood burning sources for residential heating. Overall, this study confirms the major roles of both of these combustion sources in the OP levels in ambient air.
Environmental significanceCOVID-19 lockdown restrictions have provided a rare opportunity to assess air quality during a window with significantly reduced human-led activities. Compared to other pollutants, very limited work has been done on the oxidative potential (OP) of particulate matter (PM), an emerging metric of PM exposure. Most studies used comparative methods with historical datasets, which may be insufficient to estimate real pollution levels without the lockdown restrictions in place. A machine learning technique called the Random Forest (RF) regression model was applied to predict a more realistic business-as-usual (BAU) level for OP, PM, and BC at urban background and traffic sites. While there is a clear decrease in OP found in the traffic site, more modest changes in OP were found in the urban site, likely due to sustained contributions from wood burning sources for residential heating. Overall, this study confirms the major roles of both combustion sources in the OP levels in ambient air during the lockdown period. |
Atmospheric PM can originate from anthropogenic and natural sources. Anthropogenic sources include vehicular, industrial, residential combustion, and agricultural emissions, while natural sources include sea spray, wildfires, mineral dust, and biogenic sources.9 A bigger percentage of particulate pollution is driven by anthropogenic emissions, especially traffic and residential burning, in urban areas10 and even in background sites.11,12 During the lockdown periods in China and amidst strict restrictions on anthropogenic activities, there were still occurrences of severe PM pollution events, attributed to stagnant meteorological conditions.13–17 Reduction in NOx pollution levels was found in the United States during the lockdown period,18,19 Europe,20,21 and Canada;22 however, this was not reflected in PM levels, as exemplified in France.5 These findings have offered thought-provoking insights about the effect of lockdown restrictions on air quality and urge the need of a more comprehensive assessment that includes a health-based measure of PM exposure.
The OP of PM is commonly defined as the capacity of PM to generate reactive oxygen species (ROS) in biological media23 and is currently gaining attention as an emerging health-based metric of PM exposure.24–27 Acellular (or cell-free) assays are widely used to conveniently measure OP.25,28 Frequently used assays include dithiothreitol (DTT) and ascorbic acid (AA). The DTT assay measures OP based on the ability of PM to transfer electrons from DTT to oxygen, thereby producing ROS in the form of superoxides,29 while the AA assay measures OP based on PM-induced depletion of cellular AA antioxidants.30 Each OP assay appears to be sensitive to different components of PM and has varying levels of association with different health endpoints;25 hence, the use of multiple assays is frequently recommended.26,31
So far, only a few studies have considered OP measurements during the lockdown period. Wang et al.32 (2022) reported a decrease of OP simultaneously observed with decreasing biomass burning emissions in a typical coal-combustion city in China. Paraskevopoulou et al.8 (2022) found that there was no measurable change in the OP of PM2.5, even with the significant decrease in traffic during the COVID-19 lockdown in an urban background site in Athens, Greece. At an urban background site in the metropolitan area of Milan, Italy, the OP of PM2.5 (particularly OP assessed by DTT assay) was only slightly reduced during the lockdown period with much lower change than the other air quality parameters (e.g., chemical species in PM2.5, NO2, and BC).33,34 Overall, the lockdown periods opened a unique opportunity to evaluate OP of PM as a health exposure-relevant metric, especially in areas where no significant reduction in PM levels took place during the lockdown periods.
Due to multiple interplays between several factors (i.e., meteorology, seasonality, natural and anthropogenic emission sources) affecting PM pollution, it is crucial to take all important factors into account when predicting pollution levels. The efforts to reduce PM pollution through implemented air quality policies could also lead to decreasing trend of some emissions from sources, especially from traffic.11 Hence, to obtain a forecast of a business-as-usual (BAU) scenario, the use of simple and robust machine learning models can be useful apart from traditional interpolation using historical data. In fact, the use of stochastic data-driven models for forecasting has been proven useful for various applications in air pollution studies.35–37 One of the less computationally expensive machine learning techniques is called Random Forest (RF).38 This model uses a supervised learning algorithm that uses an ensemble learning method by combining predictions from multiple machine learning algorithms to make a more accurate prediction.
The present study takes advantage of the unique scenario brought about by the COVID-19 pandemic lockdown restrictions, enabling the estimation of the impact of significantly reduced anthropogenic emissions on the OP of atmospheric PM. The BAU levels were estimated using RF regression to more realistically assess the impact of the lockdown restrictions on PM pollution. This study aimed to elucidate the reduction in PM mass concentration and OP, as well as various combustion-related components (total equivalent black carbon (BCtot), wood burning BC (BCwb), and fossil fuel BC (BCff)) in an urban background site (Grenoble, France) and a traffic site (Bern, Switzerland). The findings of this study could provide a baseline impact of the COVID-19 lockdown restrictions on PM redox activity in an urban and traffic setting in Western Europe.
000 inhabitants. This Alpine city is surrounded by several mountain ranges that favours the development of temperature inversions in the atmosphere,39 typical in a valley. The BERN site (traffic) is located in the city centre of Bern, Switzerland, a city with a population of about 145
000 inhabitants. This traffic site is situated within a busy canyon street and close to main train stations (47.0°N, 7.4°E).
The total PM10 mass concentration was measured in the GRE site using tapered element oscillating microbalances equipped with filter dynamics measurement systems (TEOM-FDMS),40 while an optical particle counter (Palas Fidas 200) was used in the BERN site. For both sites, equivalent black carbon (eBC) was measured using an aethalometer (AE33, Magee Scientific, USA).41 For the latter measurements, no further correction was applied to the concentration estimate calculated and provided by the AE33 software. The total equivalent BC concentration obtained in this way (BCtot) was apportioned to its wood burning (BCwb) and fossil fuel (BCff) fractions following the methodologies proposed in previous studies.42–44 Briefly, the contributions from liquid/fossil fuels (e.g., traffic) and solid fuels (typically wood burning) to BCtot were apportioned using the different absorption coefficients at 470 nm and 950 nm by the aethalometer. The absorption Ångström exponent (α) of 0.9 and 1.68 for the BERN site,43 while 1.0 and 2.0 for the GRE site45 were used for BCff and BCwb, respectively. The carbonaceous components, specifically organic carbon (OC) and elemental carbon (EC), were also analysed using the thermo-optical method on a Sunset Lab analyser,46 using the EUSAAR2 temperature protocol.
All of these measurements were conducted following the related EN standards (i.e., EN 12341, EN 14902, EN 16909, and EN 16913) and inter-laboratory comparison exercises for OC and EC within ACTRIS (Aerosols, Clouds and Trace Gases Research InfraStructure Network). Thanks to long-term observations at these given sites, as well as the presence of complementary monitoring stations in both areas, GRE10,47,48 and BERN49 are well characterized in terms of ambient levels and sources of PM pollution.
In the DTT assay, DTT is used as a chemical proxy that mimics the in vivo interaction of PM with biological reducing agents. This interaction represents the consumption of DTT in the assay or the capacity of PM to generate ROS.28 The leftover DTT in the initial reaction with PM is further reacted with 5,5′-dithiobis-(2-nitrobenzoic acid) (DTNB) which produces 5-mercapto-2-nitrobenzoic acid or TNB. The TNB is measured in terms of absorbance at 412 nm wavelength using a plate reader (TECAN spectrophotometer Infinite M200 Pro) and a 96-well multiwall plate (Greiner-Bio). To allow for the consumption of DTT, the absorbance is measured every 10 minutes for a total of 30 minutes of analysis time.
In the AA assay, the AA antioxidant is used in a respiratory tract lining fluid (RTFL)30 which represents the role of antioxidants in preventing the oxidation of lipids and proteins in the lung.51 The PM extract and AA mixture are analysed in terms of absorbance at 265 nm wavelength using a plate reader (TECAN spectrophotometer Infinite M200 Pro) and 96-well UV-transparent multiwall plate (Greiner-Bio) every 4 minutes for a total of 32 minutes of analysis time.
A positive control test is performed in every experiment, specifically a 1,4-naphthoquinone (1,4-NQ) solution for both the DTT (40 μl of 24.7 μM stock solution) and the AA (80 μl of 24.7 μM 1,4-NQ solution) assays. There is a < 3% CV in the positive control tests for both assays.
To improve the numerical stability of the model, both the target variable and the features were scaled using a standard scaling function. Afterwards, the dataset was split into training and testing sets by randomly assigning 70% of the input data as the training set and the rest of the 30% as testing. For the RF models for both OPDTT and OPAA in the BERN site, the test size had to be increased to 0.4 (0.6 for the training set) due to the number of data points (n = 177 and n = 273 for PM2.5 and PM10, respectively). Briefly, the training set is the dataset used to build and fit predictive models, while the test set is a subset of the dataset to assess the likely future performance of a model. A hyper-parameter tuning step was performed using a randomized search using a cross-validation technique based on a fixed number of parameter settings to determine an optimal architecture. The parameters considered and implemented were: the generated random numbers for the random forest, the number of trees in the forest, the minimum number of data points before the sample is split, the minimum number of leaf nodes that are required to be sampled, the maximum depth in a tree, applying a bootstrapping method when building decision trees, and applying a generalization score using out-of-bag samples. Typically, the hyper-parameter tuning step did not change the model performance by a large amount. Finally, the model performance was evaluated by calculating the root mean square error (RMSE) and R-squared value (r2) for both the training and testing sets.36 The optimal models were chosen based on their RMSE (as low as possible) and r2 (as high as possible) on the testing set. The modelling workflow from data collection, data processing, and the machine learning model optimization is presented in Fig. 1.
![]() | ||
| Fig. 1 Scheme of the modelling workflow of the Random Forest (RF) model for predicting business-as-usual (BAU) levels for the COVID-19 lockdown period in 2020. | ||
The percentage median change (%) based on mass concentrations (units: μg m−3 for PM, BCtot, BCwb, and BCff) or OP activity (units: nmol min−1 m−3 for OPDTT and OPAA) by period was calculated using eqn (1):
![]() | (1) |
Generally, it can be inferred that the average PM levels have not varied greatly even with the lockdown restrictions in place in the GRE site. Conversely, there was a noticeable decrease in the PM levels in the BERN site for both PM size fractions (PM10 and PM2.5), except for specific days implicated by a Saharan dust event from March 26 to 30, 2020.57 A comparison of BC mass concentration levels for both sites is also presented in Fig. 2. On the other hand, there was no large reduction in the total BC mass concentrations in the GRE site possibly due to the slight increase in BCwb. On the other hand, there seems to be a noticeable decrease in total BC in the BERN site due to the reduced BC from the fossil fuel fraction (BCff). Interestingly as well, for the OP activities, no drastic decrease is observed at each of the sites investigated in the present study.
In the GRE site, it has been consistently reported that biomass burning is a major contributor to PM10 mass during cold periods.10,58,59 Similar to other sites in Western Europe, increased residential/domestic heating emissions during lockdown due to people being at home likely affected the PM mass concentration, resulting in a less pronounced decrease in total PM levels.60 In fact, an increase in PM pollution was found in Tehran,61 China,62 and Malaysia63 potentially due to local burning activities as well. In Brazil, there is an increase in ozone levels due to the sharp decrease in nitrogen oxides resulting from decreased vehicular emissions.64 In Morocco, the decrease in local emissions was exceeded by regional and/or transboundary pollution resulting in an increase in pollution levels during the lockdown period.65 Overall, the decrease in vehicular emissions and mobility (i.e., travel restrictions) was not enough to offset other pollution sources and, potentially, the impact of secondary aerosol formation could be a more important driver of the increased pollution levels during the lockdown period.6,66
In a reference background site in France, an existing long-term decreasing trend in PM has been observed over the years.11 This has been attributed to the decreasing traffic emissions, likely due to the improvement of regulations restricting vehicular emissions, especially in bigger cities. Hence, it can be challenging to differentiate the potential reduction due to the impact of lockdown restrictions on pollution levels from an already existing decreasing trend. This has motivated an estimation of the concentrations for all variables under the hypothesis of BAU emissions to allow a more realistic assessment of pollution levels during the lockdown period in 2020.
The associations between the observed and predicted values for both the training and test sets are presented in Fig. 3 for OPAA in the GRE site. The same comparisons are provided for the other variables in the ESI (Fig. S4–S17†). The model uncertainties based on the RMSE between the observed and RF predictions are also provided in the same figures.
In general, the model performed well for most variables for both sets for the GRE site (r2 ranging from 0.90 to 0.94 for training and 0.71 to 0.89 for test sets), but some variables were harder to predict in the BERN site (r2 ranging from 0.73 to 0.92 for training and 0.40 to 0.92 for test sets). Better predictions were found in the GRE site possibly due to the higher number of data points available for the model input (Table S1†). It is also important to note that the RF models were built differently, based on the available data in each site. In the GRE site, the available long-term data of OC (a species contributing to about 23 to 28% of total PM10 in the city68) showed high importance score in the RF model, especially for PM10, BCwb, and OPAA predictions. For the BERN site, the RF model relied more on the long-term BC data from other cities (Basel and Zurich) than any of the other considered features. Even the meteorological variables in the BERN site did not significantly increase the performance of the RF model for this site.
Very few studies have estimated BAU levels in terms of OP during the COVID-19 lockdown period.57,69 In the Ciarelli et al.69 (2021) study, a regional air quality model (Comprehensive Air Quality Model with extensions, version 6.50) was used, which relies on emission inventories to model NO2, O3, and PM2.5 in Italy (Milan) and Switzerland (Zurich). The model performance for PM2.5 is associated with r2 equal to 0.41 and 0.73, with a corresponding RMSE of 9.83 and 4.08 for the sites in Italy and Switzerland, respectively. In Lovrić et al.57 (2021), a RF regression model was used to predict “true pollution” levels in multiple traffic-loaded (r2 ranging from 0.42 to 0.61) and less traffic-loaded (r2 ranging from 0.66 to 0.71) sites in Graz, Austria. Overall, the RF models in our study provided a proper generalization based on available long-term data and performed better than the models of the two other studies, allowing a realistic estimation of a BAU scenario of the pollution levels during the COVID-19 lockdown period.
| Variables | Percentage change (%) between observed vs. historical | Percentage change (%) between observed vs. predicted | ||||
|---|---|---|---|---|---|---|
| GRE | BERN | GRE | BERN | |||
| PM10 | PM10 | PM2.5 | PM10 | PM10 | PM2.5 | |
| a The asterisks depict the Mann–Whitney U-test of statistically significant difference between observed and historical/predicted observations at p ≤ 0.05 (*) and p ≤ 0.01 (**). Note: a negative% change represents a decrease in median levels during the lockdown period, while a positive% change presents an increase. | ||||||
| PM | 0.6 | −24.2** | −27.1** | 8.4 | 5.8 | 0.3 |
| OPAA | −35.2* | −17.5 | 30.8 | 12.7 | −24.9 | −39.2** |
| OPDTT | 49.5 | −41.6* | 48.3 | 3.3 | −18.7 | 21.8 |
| BCtot | −31.3* | −31.0** | −14.2 | −32.5** | ||
| BCwb | 24.0 | −17.8** | 47.1 | −16.1** | ||
| BCff | −37.7** | −44.4** | −18.5 | −40.1** | ||
Indeed, this is confirmed for the PM mass concentration in the BERN site with a clearer shift in levels when comparing the observed and historical levels. This shift pertains to a change in median levels of about −24% and −27% for PM10 and PM2.5, respectively, as shown in Table 1, implying decreased levels during the lockdown period. The predicted BAU levels were generally consistently lower than the historical levels of PM mass concentration as a result of an existing decreasing trend. Hence, there were minimal changes between the observed and the predicted BAU levels, implying a marginal effect of the lockdown restrictions on the PM mass concentration in the BERN site. In the GRE site, there was a negligible change in the observed compared to the historical level (<1%), while a slight increase when compared to the predicted BAU level (8%).
In terms of BC mass concentration, a much clearer decrease in levels (both in the BCwb and BCff fractions) was found in the BERN site. A stronger decrease was found in the BCff fraction, as expected, due to a sharp decrease in vehicular emissions during the lockdown period in this traffic site. This refers to a median% change in the observed level of −44% and −40% compared to historical and RF-predicted levels, respectively. In the GRE site, a decrease was also found in the BCff fraction (−38% and −19% for the historical and RF-predicted levels, respectively), while an increase was found in the BCwb fraction. This tends to imply a sustained contribution from wood burning emissions likely due to residential heating emissions during the lockdown period, especially in an urban area such as the GRE site. In fact, the median percentage change for BCwb in the GRE site has increased by 24 and 47% when the observed levels were compared with the historical average and the RF-predicted levels, respectively.
For OPv (both OPAA and OPDTT), the findings were less straight forward. Generally, the distributions of the observed OP were lower than the historical and BAU levels in the PM10 fraction in the BERN site. This suggests that OPv was as much influenced by the decrease in PM mass concentration in the area. This decrease in OP pertains to a median% change of −18% for OPAA and −42% for OPDTT. A similar shift was also found comparing the observed with the predicted BAU levels with −25% for OPAA and −19% for OPDTT. For the OPv of the PM2.5 fraction in the BERN site, only OPAA showed a decrease in the median% change (−39%) after comparing the observed with the predicted BAU level. In other cases, the observed was higher than the historical or the predicted BAU (OPDTT only). For the GRE site, the observed OPDTT was higher than the historical level (50%), but was about the same as the BAU levels. The observed OPAA decreased compared to the historical (−35%), but slightly increased compared to the BAU levels (13%).
The sensitivity of the assays could also play a role as BCff mass concentrations indeed have stronger associations with OPAA (Spearman correlation coefficient, rs of 0.7 and 0.3 for the GRE and BERN sites, respectively) (see Fig. S19–S21 in the ESI†). This is much weaker on OPDTT with rs of 0.3 and 0.2 for the GRE and BERN sites, respectively. The BCwb mass concentrations also showed slightly stronger correlations with OPDTT with rs of 0.5 and 0.6 for the GRE and BERN sites, respectively, than OPAA. Overall, these results present the difference in the impact of lockdown restrictions on OP between two different typologies (traffic vs. urban background sites).
On the other hand, in the BERN site, the reductions in total BC (both in wood burning and fossil fuel BC fractions) were significant, resulting in a clearer impact of the lockdown restrictions on the levels of PM and OP. This site is also close to a train station, which suffered a drastic reduction in traffic during the lockdown. In a traffic area, it is anticipated that road dust resuspension could be a major source of PM10 pollution.26 Road dust is a common source of metals containing highly reactive species leading to increased OP activity in particulates.29 With the decrease in mobility during the lockdown period, this could be the primary cause of that reduction in OP in the traffic site.
Indeed, these findings are consistent with other studies that have reported greater reductions in pollution levels in areas directly impacted by anthropogenic emissions (i.e., bigger cities, traffic areas, and industrial sites).57,66
In an urban background site in Zagreb, Croatia, it was found that a decrease in mobility (e.g., transport to/from work, grocery shopping, etc.) during the lockdown period did not have a significant influence in PM mass concentration in the area.70 A similar case was found in an urban background site in Athens, Greece, revealing that the decrease in traffic did not result in any measurable change on OP (assessed by DTT assay) of PM2.5.8 In fact, this last study suggests unchanged contributions from the biomass burning and secondary aerosols. Interestingly, this is also consistent with the findings in an urban background site in Milan, Italy highlighting that OP (also assessed by DTT assay) has been mostly impacted by biomass burning resulting in weaker reductions than other pollutants during the lockdown period.33 Meanwhile, in a typical coal-combustion city in Linfen, China, a significant decrease in OP was observed attributed to the decrease of several subgroups of primary organic aerosols (OA).32 These authors also found an enhanced formation of secondary OA resulting from high oxidation capacity (i.e., high O3 levels) in the atmosphere during the lockdown period.
Fig. 5 presents the joint distribution of PM and OPAA during the lockdown period in both the GRE and BERN sites. A similar figure is presented for OPDTT in the ESI (Fig. S22†). These figures show a difference in the distributions during the lockdown period compared to the past few years and even on a predicted BAU scenario, particularly the shift towards a lower OPAA even when the PM mass concentrations were relatively similar for both sites. This indicates two important takeaways: first, the lockdown restrictions have potentially resulted in reduced redox activity of PM in the atmosphere and second, the added value of OP as a PM exposure metric. The apportioned BC in both sites has provided some indication of the changes in the contribution of combustion-related sources during the lockdown period. However, a complete chemistry would provide a more comprehensive view of the specific impacts of reduced source emissions to PM-induced oxidative stress.
• The impact of the lockdown restrictions on OP varies based on site types and pollution profiles in the area. In an urban background site, there were more modest changes in the PM mass concentration levels, and a decrease in OP (assessed by AA assay) was found. This decrease was not found in the OPDTT likely due to the sustained contributions from wood burning sources in the urban background site during the lockdown period.
• A clearer decrease in OP was found in the traffic site, highlighting the impact of lockdown restrictions on sites with direct anthropogenic emissions in terms of minimizing PM pollution and the potential health effects upon PM exposure.
• The RF modelling technique provided a good estimate of a BAU scenario, allowing for a more representative assessment of the changes of pollution levels during the lockdown period.
Overall, this study elucidated the differences in the impacts of the COVID-19 lockdown restrictions on OP of PM in two sites in Western Europe, further highlighting the importance of considering the redox activity of PM. This was done by using a more traditional approach through a comparison with historical data and, also, a more innovative approach using a machine learning technique that allows a comparison with a more representative BAU level.
Footnote |
| † Electronic supplementary information (ESI) available. See DOI: https://doi.org/10.1039/d3ea00013c |
| This journal is © The Royal Society of Chemistry 2023 |