- 1Centro de Ciencias de la Atmósfera, Universidad Nacional Autónoma de México, Ciudad de Mexico, Mexico
- 2CONACYT, Centro de Investigación en Ciencias de Información Geoespacial, Querétaro, Mexico
- 3Instituto Nacional de Ecología y Cambio Climático, Coordinación de Contaminación y Salud Ambiental, Ciudad de Mexico, Mexico
Meteorology and long-term trends in air pollutant concentrations may obscure the results from short-term policies implemented to improve air quality. This study presents changes in CO, NO2, O3, SO2, PM10, and PM2.5 based on their anomalies during the COVID-19 partial (Phase 2) and total (Phase 3) lockdowns in Mexico City (MCMA). To minimise the impact of the air pollutant long-term trends, pollutant anomalies were calculated using as baseline truncated Fourier series, fitted with data from 2016 to 2019, and then compared with those from the lockdown. Additionally, days with stagnant conditions and heavy rain were excluded to reduce the impact of extreme weather changes. Satellite observations for NO2 and CO were used to contrast the ground-based derived results. During the lockdown Phase 2, only NO2 exhibited significant decreases (p < 0.05) of between 10 and 23% due to reductions in motor vehicle emissions. By contrast, O3 increased (p < 0.05) between 16 and 40% at the same sites where NO2 decreased. During Phase 3, significant decreases (p < 0.05) were observed for NO2 (43%), PM10 (20%), and PM2.5 (32%) in response to the total lockdown. Although O3 concentrations were lower in Phase 3 than during Phase 2, those did not decrease (p < 0.05) from the baseline at any site despite the total lockdown. SO2 decreased only during Phase 3 in a near-road environment. Satellite observations confirmed that NO2 decreased and CO stabilised during the total lockdown. Air pollutant changes during the lockdown could be overestimated between 2 and 10-fold without accounting for the influences of meteorology and long-term trends in pollutant concentrations. Air quality improved significantly during the lockdown driven by reduced NO2 and PM2.5 emissions despite increases in O3, resulting in health benefits for the MCMA population. A health assessment conducted suggested that around 588 deaths related to air pollution exposure were averted during the lockdown. Our results show that to reduce O3 within the MCMA, policies must focus on reducing VOCs emissions from non-mobile sources. The measures implemented during the COVID-19 lockdowns provide valuable information to reduce air pollution through a range of abatement strategies for emissions other than from motor vehicles.
Introduction
The coronavirus disease (COVID-19) spread rapidly around the world in early 2020, changing anthropogenic activities permanently (1). The first cases of COVID-19 in Mexico were diagnosed on 27 February 2020 (2). To control the COVID-19 outbreak in Mexico, social distancing measures were implemented following a scheme of 2 lockdown phases. Unofficially, lectures and student attendance to schools were the first activities suspended on 17 March, while all academic activities were suspended completely on 24 March after exceeding one thousand COVID-19 confirmed cases in the whole country (Supplementary Figure 1) (3), which led the Federal Government to declare the COVID-19 Phase 2. Followed by Phase 2, on 30 March, partial lockdown measures targeting non-essential activities were introduced for shopping, leisure, public administration services and public gatherings of more than 25 people (Table 1) (4). Finally, with the declaration of COVID-19 Phase 3 on 21 April, the total lockdown prohibited all non-essential activities remaining in action (5). In particular, within the Mexico City Metropolitan Area (MCMA), the total lockdown measures included the suspension of all non-essential industrial production and supplies trading, reduction of the public transport service and the enforcement of “a day without a car” program for all petrol vehicles (Supplementary Figure 1; Table 1).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/642630/fpubh-09-642630-HTML/image_m/fpubh-09-642630-t001.jpg)
Table 1. Summary of the COVID-19 outbreak propagation and lockdown measures implemented within the MCMA.
Numerous studies focusing on evaluating the urban air quality during the COVID-19 lockdown have documented significant reductions in air pollutants concentrations worldwide, associated with reduced emissions from anthropogenic activities, mainly from reduced vehicle activity (6–10) (Supplementary Table 1). The fact that lockdowns improve air quality in urban environments as a result of restrictions targeting reducing anthropogenic activities is expected when compared with prior business as usual periods. Nevertheless, the degree of air pollutants reduction may vary from city-to-city due to local factors such as the severity of lockdown measures, distribution of emission sources, meteorology and trends in pollutant emissions, complicating the precise quantification of changes in air pollutants concentrations, and their possible public health benefits during the COVID-19 lockdown (11, 12). While most studies have focused on comparing pre-lockdown air pollutant concentrations with those under lockdown, the impact of meteorology, and long-term trends on the latter has received less consideration.
The COVID-19 lockdown can be seen as an unplanned experiment to study the effect on air quality of extraordinary reductions in anthropogenic activities. Furthermore, this lockdown may help to identify an air quality baseline to achieve in non-lockdown conditions. Nowadays, most existing public policies aiming to reduce air pollution are focused on the abatement of long-term emissions and increasing energy efficiency standards. Therefore, minimising the influences of meteorology and long-term policies on air pollutant concentrations is required to precisely detect changes in their concentrations ascribed to new interventions (13–15). An even more challenging task is to extract from observations the signal of short-term interventions to reduce extremely high pollution levels since these can be obscured by meteorological fluctuations (16). For instance, air pollutants may accumulate due to air stagnation related to meteorological conditions favouring dry and stable regimes characterised by a lack of ventilation, presence of temperature inversions, and high-pressure systems with influence at synoptic scale (17, 18). Therefore, a scenario in which changes in air pollutant concentrations during the COVID-19 lockdown could be driven by meteorology rather than by reductions in their emissions must be considered.
With a population of more than 21 million people in the MCMA (19), it can be hypothesised that the unprecedented lockdown measures aiming to control the COVID-19 outbreak resulted in significant reductions in air pollutant concentrations over this period. Such reductions were presumably larger in industrial-vehicle environments than in urban background conditions. This study presents changes in CO, NO2, O3, SO2, PM10, and PM2.5 (criteria pollutants) during the COVID-19 lockdown in the MCMA. By contrast with most existing studies, the focus of this study was to (i) minimise the influence of seasonality, long-term trends and rapid weather changes on air pollutant concentrations, (ii) observe differences in pollutant concentrations for particular environments associated with improved air quality during the lockdown, and (iii) calculate public health benefits under lockdown conditions. Observations from the MCMA Air Quality Monitoring Network were obtained from representative sites to capture air pollutants behaviour before the lockdown using Fourier series modelling. Hypothesis tests were used to confirm that our models, integrated by constant deviation, long-term trend and seasonal component, fit better the data than a constant value. Predictions of air pollutants from calibrated models were used to calculate anomalies for the lockdown, and these compared to anomalies during the corresponding periods between 2016 and 2019 (baseline). A 4-year baseline was selected to compare the air pollutant concentrations during the lockdown in order to reduce the impact of inter-annual variability on air pollutant levels.
To more precisely quantify the impact of the two lockdown phases on criteria air pollutants within the MCMA, the obtained anomalies were analysed through statistical tests of probability density functions, diurnal cycles, and overall anomalies. Remote sensing observations were also used to confirm our results. Finally, an air quality health index and changes in the health burden associated to potential outdoor pollution exposure were calculated to identify potential health benefits from reduced air pollution during the lockdown. The methodology reported here can be applied in other cities where air pollutants exhibit significant seasonality, long-term trends or are severely affected by meteorology to better quantify the impact of pollution control strategies on short-term scales.
Materials and Methods
Air Pollutants Data
Criteria air pollutants have been measured continuously within the MCMA since 1986 by the Atmospheric Monitoring System (SIMAT) of the Mexico City government, and are publicly available as 1-h averages after proper validation using US Environmental Protection Agency standards (20). Air pollutants data recorded during lockdown Phases 2 and 3, and corresponding periods from 2016 to 2019, at monitoring sites representative of traffic (TRA), industrial (IND), commercial (COM), residential (RES), and urban background [upwind [UBN] and downwind [UBS]] sites were downloaded from the SIMAT website (Supplementary Table 2) (http://www.aire.cdmx.gob.mx). Henceforth the lockdown phases are referred to only as Phases 2–3. All selected sites have a minimum data capture of 75% 1-h averages for all air pollutants during the studied periods. Supplementary Figure 2 shows the selected monitoring sites location within the MCMA. Additional details of the monitoring sites location, description and characteristics can be found elsewhere (20).
Calculation of Anomalies in Air Pollutant Concentrations
To minimise the influence of seasonality and long-term trends, we used an approach based on anomalies in air pollutant concentrations during the lockdown in comparison with a 4-year baseline. The air pollutant anomalies during lockdowns were obtained by subtracting predicted data using truncated Fourier series from observed data as follows. Firstly, a modelled fitted function was computed for all air pollutants using truncated Fourier series with a frequency of k = 2 as shown in Equation (1):
where t is the normalised time, that is to say, the number of days in a given calendar year divided by 365 to fit a full sine/cosine cycle. The model captures the long-term trend component with the first two terms (c0 + c1), while the seasonal component is captured by the aggregation of all terms. The Fourier model is a natural form of modelling seasonal trends, it should be noted that a large number of frequencies could lead to overfitting, therefore, we use only two frequencies for achieving a statistically significant model that visually fits the air pollutants trends. This reduces the possibility of meteorological and seasonal bias in contrast with existing studies (10, 21). The Fourier model was calibrated using daily averages from 2016 to 2019 for all air pollutants to reduce high-frequency variability. For the lockdown periods, the modelled pollutants data were obtained from the model predictions, using P(t) with the corresponding constants c0, c1, a1, a2, b1, and b2. Finally, the hourly anomalies were calculated as shown in Equation (2):
where A(t) are the hourly anomalies, D(t) are the hourly observations and P(t) are the modelled pollutant daily data for the day t. Air pollutants modelling was made using the R software (22).
Meteorological Analyses
To discard the effects of air pollutants accumulation, days with stagnant conditions in which meteorology might have a significant influence were excluded from the analysis. Stagnant events for the lockdown and baseline periods were obtained using the air stagnation index (ASI), which is a binary index based on daily precipitation thresholds together with surface and upper wind speed (WS) (13, 14). The following daily mean thresholds must be met together to define a day as stagnant: surface WS <3.2 m s−1, WS at 500 mb <13 m s−1, and daily total precipitation <1 mm. The purpose of calculating the ASI in this work was to discard days with significant meteorological influence on air pollutant concentrations (18), rather than suggesting a cause-effect relationship, an issue that is beyond the scope of this study. The effect of heavy rain on air pollutant concentrations was minimised also by excluding days with heavy precipitation from the analysis (>10 mm as daily average) (23). The removed days account for <30% of the studied period in each year. Hourly surface meteorological data for temperature (T), relative humidity (RH), WS and wind direction (WD) were obtained from the SIMAT website. Data of WS at 500 mb were obtained from the upper-air soundings of the Mexico City International Airport and are readily available at the University of Wyoming data repository (http://weather.uwyo.edu/upperair/sounding.html). Daily precipitation data for the area encompassing the MCMA were obtained from the Level 3 product of the Global Precipitation Measurement (GPM) at 10 km resolution (24). The GPM is an international network of satellites that provide the next-generation global observations of rain and snow.
Satellite-Based Observations
Changes in surface NO2 and CO during lockdown phases 2–3 were compared with those determined in the corresponding period of previous years using satellite-based observations. As for surface observations, days with stagnant conditions and heavy rain were excluded (see section Meteorological Analyses) from the NO2 and CO satellite-based analyses. For NO2, we used the standard product V003 (25) from the Ozone Monitoring Instrument (OMI) onboard the Aura satellite of the National Aeronautics and Space Administration (NASA), downloaded from the Earth Data Website (https://earthdata.nasa.gov/). From the OMI NO2 data product (26), the tropospheric column was used and filtered with a cloud fraction lower than 20% and quality flag 0. NO2 distribution maps were constructed using data averaged for the 2016–2019 lockdown corresponding period and compared with 2020 data. Specific details of the methodology for the construction of maps can be found in (27).
The averaged distribution of the CO column over the MCMA during the lockdown phases was reconstructed from the TROPOMI satellite instrument. The TROPOspheric Monitoring Instrument is, as in Nadir geometry, a measuring instrument on-board the Sentinel 5P satellite, which has a sun synchronised orbit and a daily equatorial crossing time at 13:30 LST and a similar overflight time for the MCMA. The swath of the instrument provides daily coverage over the MCMA and data from the level 2 product were downloaded (28, 29). The distribution maps were reconstructed using the oversampling method reported in UNIATMOS (30) with a defined mesh of 50 × 50 grid points. The calculated CO concentration at each grid point represents the average of all measurements with a pixel centre closer than 5 km to each grid point. Five kilometre was empirically determined as a trade off between the mesh grid (ca. 2 km), the noise and density of the data for this period, the footprint of the TROPOMI-measurement and the inhomogeneity expected for the MCMA.
Mobility Analyses
Associations between air pollutants and mobility were identified from correlations computed between all mobility types reported in the Google COVID-19 Community Mobility Report and concentrations recorded at the selected monitoring sites. Mobility data for the MCMA were obtained from the Google COVID-19 Community Mobility Report (31), which shows changes in visits and length of stay at different places compared to a baseline. The mobility categories reported are retail, grocery, parks, transport stations, residential and workplaces. Further details of how the mobility trends are calculated and descriptions of the place categories can be found on the Google COVID-19 Community Mobility Reports website (31). The subset from 17 March to 21 April 2020 of Google mobility data was selected considering the beginning and end of major mobility changes within the MCMA as shown in Supplementary Figure 1. Significant correlations that passed a hypothesis test with a null hypothesis of 0 correlation are reported with a significance value of α = 0.05, i.e., a high positive significant correlation means that the decrease of an air pollutant is highly correlated with the decrease in mobility.
Calculation of Potential Health Benefits
The potential health benefits for the MCMA population ascribed to changes in air pollution during the COVID-19 lockdown were addressed using (i) an air quality health index that can reflect the overall improvement in air quality and (ii) calculating changes in the excess risk (ER) of premature mortality. Here, we use the AQHI of Canada to provide insights into the effects of combined changes in NO2, O3, and PM2.5 during the lockdown (32). It measures the air quality on a scale from 1 to 10 which is assigned to a category that describes the level of health risk associated with the index calculated from Low to Very High Health Risk (>10). The AQHI is calculated on a city basis (one or more monitoring sites) and is based on 3-h moving average concentrations as follows (33):
i) The average concentrations of O3, NO2, and PM2.5 are calculated at each site for the preceding 3-h with a threshold of 2 out of 3-h at each site, obtaining three pollutant averages for each site. If more than 1 of the preceding 3-h is missing, the site average was not calculated.
ii) The 3-h city average for each parameter is calculated from the 3-h pollutant averages at the available sites, obtaining 3 city parameter averages. If no sites were available for a parameter, such parameter was not calculated.
iii) If the three city parameter averages are available, the city AQHI was calculated. Equation (3) shows the AQHI formulation:
where all pollutants are entered as 3-h moving average concentrations in ppb for O3 and NO2 and in μg m−3 for PM2.5.
Changes in ER-associated with variations in air pollutant concentrations during the lockdown were calculated as follows. Firstly, the relative risk (RR) associated with short-term exposure to criteria air pollutants was calculated using Equation (4) (34):
where RRi corresponds to the RR of the i-pollutant, βi is the exposure-response constant indicating the additional health risk per unit of the i-pollutant after exceeding a threshold, and ci is the daily average for pollutant i. The World Health Organization (WHO) reported that adverse health effects can occur at any concentration of NO2, O3, PM, and SO2, apart from CO (35). Therefore, we assumed a value of c0 = 0 for all air pollutants, which implies that pollutant concentrations equal to this value are associated with no excess risk (i.e., RR = 1), while for CO c0 was equal to 2 mg m−3 as reported in (36). The β values used here were 0.038, 0.032, 0.081, 0.13, and 0.048% per μg m−3 of PM2.5, PM10, SO2, NO2, and O3, respectively, while for CO a value of β = 3.7% per mg m−3 was used (34). Finally, the ER for the i-pollutant was calculated using Equation (5) (37).
Additionally, the potential health benefits in the MCMA during the lockdown were quantified by estimated changes in the health burden (ΔHB) [number of premature deaths; Equation (6)] due to outdoor air pollution exposure following the methodology reported in (38, 39). The averted health burden (ΔHB) during Phase 2 and 3 for the i-pollutant was estimated using Equations (6–8).
where BM is the baseline mortality per 100,000 inhabitants of all age groups, obtained from the Global Burden of Disease study of 2017 (40). ExpPop is the exposed population calculated by multiplying the MCMA population by a 75% factor obtained from the Historical Analysis of Population Health Benefits associated with Air Quality in Mexico City during 1990–2015 (41). Finally, AF represents the attributable fraction for the RR associated with the i-pollutant load during Phases 2–3 and the corresponding reference periods (38).
Results
Meteorology in the MCMA
Supplementary Figure 3 shows the monthly anomalies in T, RH and WS within the MCMA for the lockdown Phases 2–3 within the MCMA (March-May 2020). The temperature anomalies suggest a slightly warmer than the average period during March and April between 1 and ~1.5°C, while RH anomalies suggest a decrease only in May. WS anomalies showed similar differences from March to May with the long-term average. Nevertheless, the analysis with the non-parametric Wilcoxon test showed that the observed differences in T, RH and WS were not statistically significant at α = 0.05. Similarly, the WD anomalies did not show significant differences (p < 0.05) for the prevailing northerly component between the baseline and lockdown period for March, April and May (dominant WD ~2, ~350, and ~351°, respectively). This suggests that, on average, T, RH, WS and WD were similar during the COVID-19 lockdown in comparison with their long-term averages.
Modelling of Air Pollutants Time-Series
Air pollution within the MCMA peaked during the early 1990s and has decreased ever since as a result of the introduction of emission control policies (42). Although the largest decreases in air pollutant concentrations occurred during the 1990s, some air pollutants still exhibit monotonic trends. Figure 1 shows daily averages for all air pollutants at the representative sites within the MCMA from 2016 to 2020, the adjusted curve using Equation (1) and the linear trend of the time-series. An F-test about the model used in this study suggests that the regression computed using truncated Fourier series performs better than the linear model solely (Supplementary Table 3). Furthermore, the residual standard error was <10% for all air pollutant time-series, apart from PM10 which showed errors <17% at all sites. This suggests that our model captures the air pollutants behaviour significantly well. It should be noted that the residual standard error is proportional to the magnitudes of observations, and for all air pollutants shows a good fitting.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/642630/fpubh-09-642630-HTML/image_m/fpubh-09-642630-g001.jpg)
Figure 1. Time-series of criteria air pollutants from 2016 to 2020 at the six selected monitoring sites within the MCMA. The continuous lines show the Fourier fitting and modelled linear trend.
Overall, all air pollutants exhibited clear seasonal cycles during 2016–2020, while apparent monotonic trends were only observed for CO, NO2 and O3 (Table 2). All CO decreasing rates were significant (p < 0.001) and of similar magnitude (0.1 ppm yr−1; 11–15% yr−1), while NO2 decreased by 1.1–1.2 ppb yr−1 (3–6% yr−1) at all sites (p < 0.05), apart from UBN, where it decreased by 0.3 ppb yr−1 (1% yr−1; p < 0.1). In contrast, O3 has increased (p < 0.1) between 0.9 and 1.9 ppb yr−1 (4–9% yr−1), apart from the background sites (UBS and UBN) where no apparent trends were detected likely because of its location on the periphery of the MCMA. PM10 increased significantly only at the downwind UBS site (1% yr−1), whereas PM2.5 decreased between 0.5 μg m−3 at TRA and 0.8 μg m−3 at COM and IND (~2.5% yr−1). The existence of monotonic trends highlight that robust analyses are required to determine whether air pollutants changed within the MCMA during the COVID-19 lockdown.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/642630/fpubh-09-642630-HTML/image_m/fpubh-09-642630-t002.jpg)
Table 2. Trends for criteria air pollutants expressed in units of concentration during 2016–2020 at the selected monitoring sites within the MCMA.
The Air Pollutant Anomalies
The analysis of anomalies can provide robust information about odd air pollutants data (43), like those expected during the COVID-19 lockdown. To better isolate the COVID-19 lockdown effects in the MCMA, pollutant hourly anomalies for the Phases 2–3 calculated as described in section Calculation of Anomalies in Air Pollutant Concentrations were compared with those during the 2016–2019 corresponding periods (4-year reference), instead of year-to-year, which can reduce inter-annual variation caused by extraordinary local up to synoptic events (44, 45). Additionally, days with stagnant and heavy rain conditions were excluded from the analysis to minimise rapid weather changes (16, 18, 46). Supplementary Table 4 summarises the minimum, average (mean ± SD) and maximum hourly anomalies during 2016–2019 and Phases 2–3. On average, no significant differences (p > 0.05) in the minimum values of anomalies were observed in Phase 2 at all sites for all air pollutants, while the averages were significantly lower, apart from CO and O3 that showed increases.
A marked reduction in the maximum anomalies was detected during Phase 2 compared with those during 2016–2019 for all air pollutants, apart from O3 and PM10. Figure 2 shows the probability density function (PDF) for hourly anomalies calculated during Phase 2. Maximum PDF values for CO, NO2, and SO2 were higher in Phase 2 than during 2016–2019 due to lower dispersion of anomalies in the upper quartiles. O3 anomalies showed a lower maximum PDF at upwind and central sites in Phase 2 than during 2016–2019, while at RES and UBS O3 exhibited PDFs consistent with those for CO, NO2, and SO2. In contrast, PM10 and PM2.5 PDFs showed less variation at background sites (UBN and UBS) but for central sites, the maximum PDF values occurred during Phase 2.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/642630/fpubh-09-642630-HTML/image_m/fpubh-09-642630-g002.jpg)
Figure 2. Probability density function plots for hourly anomalies for CO, NO2, O3, PM2.5, PM10, and SO2 during the COVID-19 lockdown Phase 2 and baseline within the MCMA.
NO2, PM10, PM2.5, and SO2 showed higher minimum and lower average anomalies during Phase 3 than in 2016–2019 in good agreement with PDFs values observed in Phase 2, while CO averages were significantly lower in Phase 3 only at UBN and RES (Figure 3). In contrast to Phase 2, O3 maximum anomalies decreased in Phase 3 but the average and minimum anomalies were higher than during 2016–2019. PDF showed higher maximums during Phase 3 than from 2016 to 2019 for all pollutants and sites, apart from O3. No marked differences in PDFs were observed for all air pollutant anomalies at the UBN in comparison with previous years, suggesting that on average, anomalies were similar in background conditions. However, this was not observed for UBS that showed higher PDFs for negative anomalies in Phase 3 excluding O3. Indeed, higher densities of negative anomalies for PM10, PM2.5, NO2, and SO2 were observed in Phase 3 at most of the sites. This suggests that clear reductions in emissions of these air pollutants occurred within the MCMA only in Phase 3 under the total restriction of all non-essential activities. Furthermore, the differences in PDF distributions between Phase 2 and 3 suggest marked spatial variation in air pollutants concentrations within the different environments of MCMA.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/642630/fpubh-09-642630-HTML/image_m/fpubh-09-642630-g003.jpg)
Figure 3. Probability density function plots for hourly anomalies for CO, NO2, O3, PM2.5, PM10, and SO2 during the COVID-19 lockdown Phase 3 and baseline within the MCMA.
Air Pollutants Diurnal Cycles
The air pollutant diurnal cycles arise from changes in meteorology and pollutant emissions, mostly from fossil fuel combustion (47–49). The lockdown measures within the MCMA included restrictions to population mobility, industrial production and office and leisure activities which may have reduced air pollutant emissions. Here, the effect of restrictions on normal anthropogenic activities in the MCMA context is analysed by comparing de-trended air pollutant diurnal cycles constructed from hourly anomalies during Phases 2–3 with the corresponding anomalies during 2016–2019. Changes in the amplitude values (AVs) of diurnal cycles during the lockdown may reflect local changes in emission sources (50). AVs of diurnal cycles for all criteria air pollutants were calculated as the difference trough-to-peak and the significance of changes in AVs was evaluated through Welch's t-tests.
Figure 4 shows that during Phase 2, changes in AVs were evident mostly at sites in the MCMA centre, while background sites on the MCMA periphery showed less marked changes. SO2 AVs exhibited significant (p < 0.05) decreases at all sites between 58 and 72% probably due to the low concentrations typical within the MCMA. O3 AVs did not change at any site in Phase 2 despite decreases in NO2 of 38 and 24% at RES and TRA, respectively. RES was the only site that showed significant decreases in AV for all air pollutants, apart from O3 which did not change, and only UBS showed increases both in PM10 and PM2.5 AVs. During Phase 3, most air pollutants AVs showed decreases within the MCMA (Figure 5), apart from O3 which did not change. As in Phase 2, only RES showed significant (p < 0.05) decreases in AVs for almost all air pollutants and it was the only site where CO decreased (58%) in Phase 3. NO2, PM10, and PM2.5 AVs decreased at all sites in the MCMA centre between 16 and 32, 29–45, and 31–47%, respectively. SO2 AVs also decreased in Phase 3, apart from IND, between 21 and 48%, but increased at UBS (8%). The observed changes in AVs suggest that reductions in air pollutant emissions were generalised and significant within the MCMA only after the restriction of all non-essential activities.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/642630/fpubh-09-642630-HTML/image_m/fpubh-09-642630-g004.jpg)
Figure 4. Average diurnal cycles within the MCMA during 2016–2019 and lockdown Phase 2. De-trended diurnal cycles were constructed by subtracting daily averages modelled with Equation (1) from hourly averages to remove the impact of long-term trends. The shading shows the 95% confidence intervals of the average, calculated through bootstrap re-sampling for 1,000 iterations. The numbers show significant changes in diurnal amplitudes, through-to peak, expressed as percentage with 2016–2019 cycles as baseline. *Level of significance p < 0.05. **Level of significance p < 0.01. ***Level of significance p < 0.001.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/642630/fpubh-09-642630-HTML/image_m/fpubh-09-642630-g005.jpg)
Figure 5. Average diurnal cycles within the MCMA during 2016–2019 and lockdown Phase 3. De-trended diurnal cycles were constructed by subtracting daily averages modelled with Equation (1) from hourly averages to remove the impact of long-term trends. The shading shows the 95% confidence intervals of the average, calculated through bootstrap re-sampling for 1,000 iterations. The numbers show significant changes in diurnal amplitudes, through-to peak, expressed as percentage with 2016–2019 cycles as baseline. *Level of significance p < 0.05. **Level of significance p < 0.01. ***Level of significance p < 0.001.
Net Changes in Air Pollutants During the Lockdown
The overall changes in air pollutant concentrations during the COVID-19 lockdown within the MCMA were determined by comparing hourly anomalies calculated for the lockdown period with those from the 2016–2019 corresponding period. The significance of the obtained changes was calculated with the non-parametric Wilcoxon test. Figure 6 shows that during Phase 2, only NO2 decreased significantly (p < 0.05) between 3 ppb (10%) and 8 ppb (23%) at UBN and TRA, respectively, with an average decrease of 4 ppb (around 20%). By contrast, significant increases (p < 0.05) were observed for O3 ranging from 7 ppb (16%) at COM to 11 ppb (40%) at TRA and UBS, likely because of reduced O3 titration by NO2. CO only increased at UBS (0.1 ppm, 4%), whereas SO2, PM10, and PM2.5 did not exhibit significant (p > 0.05) changes.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/642630/fpubh-09-642630-HTML/image_m/fpubh-09-642630-g006.jpg)
Figure 6. Net changes in anomalies of criteria air pollutants during lockdown Phases 2–3 in comparison to 2016–2019 anomalies as baseline within the MCMA. The stars show significant changes during the COVID-19 lockdown calculated with the non-parametric Wilcoxon test, while the vertical bars show the 95% confidence intervals. *Level of significance p < 0.05. **Level of significance p < 0.01. ***Level of significance p < 0.001.
During Phase 3, NO2 decreased (p < 0.05) significantly between 7 and 10 ppb (average decrease of 43%) at most sites, which are larger decreases than in Phase 2. Similarly, significant decreases (p < 0.05) were observed at all sites both for PM10 and PM2.5 during Phase 3, apart from UBN. Nevertheless, while PM2.5 decreased consistently by 7–8 μg m−3 (20%), large spatial variability was observed for PM10 which decreased between 10 and 18 μg m−3, with an average decrease of 32%. SO2 only decreased significantly at TRA (2 ppb, 31%) likely ascribed to reduced emissions from heavy-duty vehicles and buses. O3 showed large spatial heterogeneity in Phase 3, with significant increases (p < 0.05) seen only at TRA (4 ppb, 14%) and UBS (7 ppb, 19%). Our results allow us to hypothesise that the decreases observed in concentrations of primary pollutants during Phase 3 within the MCMA can be ascribed to the strict lockdown focused on controlling the COVID-19 spread.
To better observe the effects of long-term trends and rapid weather events on air pollution, we compared changes in air pollutants calculated from raw observations with those from anomalies including and excluding stagnant conditions and heavy rain events during Phase 2 (Supplementary Figure 4) and Phase 3 (Supplementary Figure 5). The largest differences were observed clearly for CO, with overestimated changes > 10 times in both Phase 2 and 3, due to the influence of long-term trends. For the other air pollutants, most of the changes during Phase 2 could have been overestimated between 0.5 and 2 times without accounting for the rapid weather changes. Similarly, an overestimation generalised of around 2 times was observed for all air pollutants, excluding CO, during Phase 3. In particular, the overestimation of O3 increases may mislead the design of future control policies by making it difficult to scale emission inventories and affecting photochemical models performance.
Changes in the NO2 and CO Columns Over the MCMA
Figure 7 shows averaged NO2 tropospheric column distribution maps constructed using satellite-based data above the MCMA for the Phases 2–3 and corresponding periods during 2016–2019. Reductions in the amount of NO2 molecules present in the atmospheric column were observed for both lockdown Phases, with the largest decreases observed in Phase 3. Spatially, the largest decrease was seen in the north of MCMA where the maximum NO2 column values decreased from ~1.1 × 1016 to ~0.8 × 1016 molecules (molec.) cm−2 from 2016–2019 to Phase 2 (−30%), and from ~0.8 × 1016 to ~0.4 × 1015 molec. cm−2 during 2016–2019 to Phase 3 (−50%). The MCMA central and southern regions exhibited decreases of between 10 and 40% in the NO2 column during Phase 3, which aligns well with the decrease calculated from NO2 surface observations. The difference in NO2 decreases between Phase 2 and 3 can be explained by a gradual reduction in motor vehicle circulation during Phase 2, while large reductions in NOX emissions from motor vehicles and other combustion sources occurred in Phase 3 because of the total lockdown.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/642630/fpubh-09-642630-HTML/image_m/fpubh-09-642630-g007.jpg)
Figure 7. NO2 tropospheric column distribution maps constructed using satellite-based data from OMI over the MCMA during 2016–2019 (baseline) and lockdown Phases 2–3.
Figure 8 shows averaged CO column distribution maps from satellite-based data over the MCMA during the lockdown Phases 2–3 and the corresponding periods in 2019. Differences between the CO columns of 2019 and Phase 2 were noticed for the central and northern MCMA regions, which corresponded to decreases of around 20 and 33%, respectively. Such differences were less evident for the southern MCMA region. No marked changes in the CO column were seen during Phase 3 compared to 2019 over most of the MCMA. Only over Phase 3, the satellite-based data are consistent with surface observations that suggest no significant changes in CO. Stremme et al. (51) reported that the CO column over MCMA is occasionally influenced by biomass burning events, which could have biased the satellite data during Phase 2 because only data for 2019 was available to use as a baseline. Although the CO is emitted mostly by motor vehicles within the MCMA (48), it is not clear why no reductions were observed during the lockdown despite reduced road traffic. A plausible explanation could be an increase in CO emissions during the lockdown from enhanced domestic liquid petroleum and natural gas burning because of the stay-at-home order, however, further investigation is required to clarify it.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/642630/fpubh-09-642630-HTML/image_m/fpubh-09-642630-g008.jpg)
Figure 8. Averaged CO column distribution maps constructed from satellite-based data from TROPOMI over the MCMA during lockdown Phases 2–3 of 2020 (right) with respect to the previous year (left).
Changes in AQHI and Excess Risk
The AQHI was used here because it can provide information about the additive effects of O3, NO2, and PM2.5 which are known to harm human health. Supplementary Figure 6 shows the boxplot of maximum daily AQHI values during 2016–2020 and Table 3 lists the calculated changes. During Phase 2, all sites showed similar maximum AQHI daily values of 15–16. Nevertheless, only the decreases (ΔAQHI = −3) at RES, IND and COM were significant (p < 0.05). The maximum AQHI daily values also decreased significantly (p < 0.05) at all sites during Phase 3 from 2 to 6 units (average = 4). The decrease in the maximum daily AQHI values during the lockdown suggest that, despite increases in O3, the reductions in NO2, and PM2.5 could represent a potential benefit for the MCMA air quality.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/642630/fpubh-09-642630-HTML/image_m/fpubh-09-642630-t003.jpg)
Table 3. Averages and changes in maximum daily AQHI values during 2016–2019 and lockdown Phases 2–3 by monitoring site within the MCMA.
We calculated the ER of premature mortality for each criteria air pollutant to evaluate possible health benefits within the MCMA during the lockdown as a result of changes in pollutant concentrations (Table 4; Supplementary Figure 7). A significant decrease (p < 0.01) in the ER of premature mortality was observed during Phase 2 only for NO2 (27%) while for O3 it increased by 25%. During Phase 3, the ER of premature mortality showed significant decreases (p < 0.01) for NO2, PM10, PM2.5, and SO2, with the largest decrease observed for NO2 (36%) and the lowest for PM10 (24%). Finally, the health assessment conducted suggested that around 588 deaths were averted during the lockdown, explained by decreases in the health burden from the reference period to the lockdown. This decrease in health burden can be ascribed to reduced population exposure to outdoor air pollution during the lockdown. Nevertheless, while around 152 deaths were averted during Phase 2, an estimate of averted deaths of 436 was found for Phase 3. This is in line with a decline in the ER of premature death larger during Phase 3 than for Phase 2. Together the AQHI and ER of premature mortality analyses can confirm that a significant decrease in potential outdoor exposure to air pollutants occurred in the MCMA as a result of the lockdown measures and health benefits that occurred because of these measures.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/642630/fpubh-09-642630-HTML/image_m/fpubh-09-642630-t004.jpg)
Table 4. Average excess risk within the MCMA for each criteria air pollutant during 2016–2019 (Reference) and COVID-19 lockdown Phases 2–3.
Discussion
The Origin of Changes in Air Pollutants During the Lockdown
Urban mobility showed significant reductions worldwide during the lockdowns (Supplementary Table 1). Within the MCMA, mobility in all categories decreased between 50 and 65% during Phase 2 (Supplementary Figure 1), apart from the residential which increased by 25%, and all stabilised in Phase 3. This suggests that during Phase 3, the population spent longer periods in residential environments because of the restrictions on mobility and labour activities (31). A generalised correlation was observed between mobility and the decreases for NO2 (p < 0.05), contrasting with O3 which showed negative correlations (Supplementary Table 5). This may confirm that the reductions observed in NO2 of up to 60% within the MCMA during Phase 3 were ascribed to the overall reduction of motor vehicles circulation (21, 52).
Lower emissions from fossil fuel combustion by motor vehicles can be confirmed by decreases in average petrol and diesel sales within the MCMA during the lockdown phases compared to 2016–2019 sales (Supplementary Figure 8) (53). Furthermore, the largest decrease in petrol consumption was observed during Phase 3, thus suggesting a decrease in VOCs emissions but only from incomplete combustion. Although CO also showed correlations with mobility at some sites, the changes in CO were non-significant, despite a reported 28% reduction in CO emissions (52). This could be due to increased domestic burning of natural and liquified petroleum gas because of the time spent at home, as suggested by the mobility analysis, however, further analyses are required to clarify this. SO2, PM10 and PM2.5 showed less marked correlations despite the reduction in diesel consumption as motor vehicles are not their major source of emissions.
O3 formation in the MCMA central-southern regions is VOC-limited but NOX-limited on the outskirts (54, 55). The reductions in NOX emissions from motor vehicles could thus explain the observed increases in O3 during Phase 2; as a result of reduced O3 titration (NO+ O3 = NO2 + O2). VOC emissions did not present significant reductions since some activities related to area emissions such as residential sources did not decrease or even increase during Phase 2. Marked decreases in petrol and diesel sales within the MCMA during the lockdown suggest a decline in VOC emissions from incomplete combustion accentuated in Phase 3 (Supplementary Figure 8). The generalised stabilisation of O3 as shown by the averaged anomalies during Phase 3, suggests that combined reductions of NOX and VOC emissions occurred when vehicle emissions troughed and most of the stationary and area emission sources other than domestic were out of operation. This hypothesis can be supported by studies addressing the O3 weekend effect (56, 57) and the power plant influence located north of the MCMA (58), which have reported that in VOC-limited conditions, reductions in NOX emissions even combined with reductions in VOCs, would lead to increments in the O3 levels.
The MCMA exhibits currently low SO2 levels because since 2009 the diesel sold in this is ultra-low sulphur (42, 48), which may explain that SO2 showed the lower changes among all air pollutants during the lockdown. However, the decrease in SO2 seen at TRA may suggest that additional reductions in SO2 emissions from heavy-duty vehicles, which currently contribute with around 21% of total emissions, could be obtained by increasing the availability of ultra-low sulphur diesel in other areas of central Mexico. Indeed, it is required for the operation of less polluting heavy-duty technologies such as EPA−10 and Euro VI. On the other hand, reductions in the operation of stationary and area sources appeared to have no significant effect on SO2 levels despite contributing to 33 and 46% of total SO2 emissions (48). PM10 and PM2.5 only decreased significantly in Phase 3 when their major emitters, stationary and area sources, were out of operation thus suggesting that to reduce their levels air quality policies must include undoubtedly these sources.
Air Pollutants Changes in the MCMA in the International Context
Marked decreases in NO2 emissions of up to 30% were observed worldwide due to reductions in transport activities implemented to control the COVID-19 pandemic in urban centres (59). In the MCMA, NO2 concentrations decreased between 10 and 50% during the lockdown phases in good agreement with reduced energy and fuel demand as observed in other cities, confirming that motor vehicles represent the major share of NO2 urban emissions (48). Increases in O3 observed during the lockdown within the MCMA are consistent with those in major urban areas around the world where O3 production has been determined as VOC-limited. For instance, this behaviour was observed in Barcelona, Spain (10), São Paulo, Brazil (60), and Delhi, India (61), and overall was ascribed to large NOX emission reductions in VOC-limited environments. The O3 observations during the COVID-19 pandemic within the MCMA suggest that measures focused only on reducing NOX or VOCs–NOX combined emissions can lead to increases in O3. Nevertheless, qualitative and quantitative data of ambient VOC levels are required to gain a thorough understanding of the effect of precursor emission reductions on O3 levels.
Area sources represent the major share of PM10 (58%) and PM2.5 (47%) total emissions within the MCMA, followed by motor vehicle sources with 29 and 36%, respectively. These contributions aggregated, account for 87 and 83% of PM10 and PM2.5 total emissions, which highlights that both pollutants decreased at all sites only after the combined total restriction of non-essential commercial, industrial and mobility activities occurred during Phase 3 (62). This is consistent with the change in PM10 levels observed in Rio de Janeiro by Dantas et al. (63) who reported that PM10 decreased only when most of the lockdown measures were in force. That said, Wang et al. (64) reported that pollution events occurred in 9 Chinese cities due to unfavourable meteorology, despite decreases in PM2.5 emissions during the lockdown. Thus, the observed reductions both in PM10 and PM2.5 within the MCMA suggest that only combined strategies targeting different emission sources can be effective for reducing PM. Only TRA showed a decrease in Phase 3 in SO2, highlighting the fact that MCMA is a low SO2 city with the major contribution of heavy-duty vehicles. An opposite scenario was experienced in northern India where upwind emissions from power plants increased the SO2 levels due to the enhanced domestic energy demand because of the lockdown (45).
Several studies have addressed changes in population mortality due to air pollution exposure during the lockdown. Here, we used city average relative risks (RR) for each pollutant by lockdown phase with the aim to include an estimation of averted deaths relative to the lockdown measures applicable. In line with the estimated reductions in health burden in five Indian cities reported by Kumar et al. (39), we can remark that air pollution-related deaths declined during the lockdown period. We estimated a higher number of averted deaths at city scale than that of Kumar et al. (39) likely because they only studied the impact of PM2.5. Similarly, Sharma et al. (45) reported that the average decrease in the ER premature mortality for all air pollutants in India observed under lockdown could avert around 0.65 million deaths in a year if pollutant concentrations decrease. For the MCMA, around 3,385 deaths per year could be averted if the ER of premature mortality observed during Phase 3 was maintained, while this estimated would decrease to 1,852 averted deaths with the ER of premature mortality decrease observed during Phase 2. In conclusion, these estimates suggest that future air quality policies can be oriented using the information gathered during the COVID-19 lockdown to protect the health of the MCMA inhabitants.
Contribution and Limitations
The dry-hot season between March and May in the MCMA is characterised by recurrent pollution episodes of O3, PM10, and PM2.5 (49), associated with calm winds, low RH and low-temperature fluctuations at synoptic scale (16). For instance, during 10–17 May 2019, pollutant emissions from wildfires in central and southern Mexico affected severely the air quality within the MCMA during the occurrence of a high-pressure system (65). In an effort to reduce the O3, PM10, and PM2.5 levels, the MCMA local authorities implemented extraordinary measures by the beginning of the episode such as reducing motor vehicle circulation, restricting the maintenance and operation of LP gas plants, suspension of highly polluting industries and all cleaning and maintenance activities. Despite such measures, the episode lasted for 5 more days after their implementation. This event highlighted how extreme weather conditions in the MCMA may obscure the benefits of extraordinary measures introduced to reduce air pollutant emissions and should be taken into account when conducting public policies evaluation.
Indeed, air pollution episode action-plans may yield significant benefits in terms of primary and secondary pollutants load in the airshed during these episodes. However, the potential benefits of those plans may not be properly assessed during pollution episodes because, in most cases, these were triggered by external drivers such as strong high-pressure systems or forest fires. The air pollutant emission reductions during the lockdown Phase 3 in the MCMA were larger than any achievable reduction resulting from the implementation of extraordinary measures, such as those implemented during the 2019 pollution episode (65). Thus, the COVID-19 lockdown allowed us to assess the potential benefits of potential public policies on air quality under non-episodic conditions.
We also aimed to minimise the influence of long-term trends and inter-annual variability on air pollutant concentrations during the COVID-19 lockdown since some pollutants exhibited monotonic trends. The comparison between changes from raw data with those from data including and excluding stagnant and heavy rain events during Phases 2–3 is shown in Supplementary Figures 4, 5, respectively. CO exhibited the largest differences between estimated changes with an overestimation of changes >10 times for raw data during both phases. The other air pollutants showed less variability in Phase 2, with most overestimations ranging between 0.5 and 2 times for raw data. During Phase 3, a generalised overestimation of around 2-fold was observed for all air pollutant changes, apart from CO. On the other hand, a comparison between changes in air pollutants calculated using the 4-year baseline with those year-to-year showed differences <10% for all pollutant and sites, apart from O3, PM10, and PM2.5 at UBN (Supplementary Table 6). This is explained by large inter-annual pollutant variability at UBN. Therefore, the 4-year baseline used here may both minimise the impact of inter-annual variability and long-term trends (42, 44).
Since the photochemistry is a non-linear process, it requires a set of model sensitivity analyses and field experiments to quantify the impacts on emission reductions on air pollutant concentrations. In particular, the overestimation of O3 increases may mislead the design of future control policies by making it difficult to scale emission inventories and affecting the performance of photochemical models. In this sense, the COVID-19 lockdown measures reduced anthropogenic activities providing an opportunity to observe a real-life experiment for emission reductions and their effects on the air quality in the MCMA. By precisely quantifying the reductions and air quality changes during the lockdowns, it is possible to calibrate photochemical models and scale the emissions inventory to evaluate other policy strategies for improving the air quality in the studied area. We suggest that those strategies can be applied in other regions around the world with similar air quality problems.
Finally, it is important to bear in mind that our findings can be limited to the analysis of data available in the urban environments selected. While the monitoring sites selected in the MCMA may represent typical environments existing in most cities across the world, future analyses must include a larger number of monitoring sites to better capture the air pollution spatial and temporal variation. Also, the public health benefits related to reduced air pollution cannot be extrapolated to other cities and can be subject to uncertainty due to the variables used in the calculation process. Nevertheless, such analyses may contribute to a better understanding of the impact of extraordinary public policies on air pollution and public health.
Conclusions
The MCMA megacity has a long history of air pollution problems. MCMA local authorities have regularly implemented and updated measures to reduce pollutant emissions since the early 1990s, however, the COVID-19 lockdown represents a unique opportunity to evaluate a what-if scenario. While most of the published studies about air pollution changes during the lockdown periods have based their results on raw observations, we presented results based on air pollutant anomalies which are more resistant to meteorological influences and long-term trends influences. Hypothesis tests at a confidence level of α = 0.05 confirmed that our models fitted more accurately the air pollutant data than means. Therefore, the calculated differences in air pollutant anomalies provided confident evidence of the major drivers of air pollutant changes during the lockdown. Monitoring sites representative of different environments were selected based on air pollutant data captured during the lockdown within the MCMA. Although the selected sites showed consistent results for all air pollutants and with those derived from remote sensing data, additional data for the CO satellite column could help to confirm the non-significant changes.
Our results suggest that while NO2 decreased significantly because of reduced motor vehicle emissions, other criteria pollutants require more stringent emission control policies for their abatement within the MCMA. These results confirm that to reduce O3 levels within the MCMA, future strategies must be introduced to reduce VOC emissions from other sources rather than just from motor vehicles. Overall, air quality improved during the lockdown in response to reduced NO2 and PM2.5 emissions despite the increases in O3 levels. A health assessment conducted suggested an estimate of around 588 averted deaths as a result of reduced air pollution during the lockdown. The analysis presented here could be extended to a regional scale to assess the influence of air pollutant emission sources on the MCMA outskirts and surrounding states. Finally, the measures implemented during the COVID-19 lockdown provide valuable information to reduce air pollution through a range of abatement strategies for emission other than from motor vehicles, which have received the most attention, particularly in economically developing countries.
Data Availability Statement
Publicly available datasets were analyzed in this study. This data can be found at: http://www.aire.cdmx.gob.mx/default.php.
Author Contributions
IH-P and MG conceived the study, participated in its design, and coordination. IH-P, SV, VA, CR-C, WS, LR-S, AG-R, and MG collected the data, entered study data, assisted in data analysis, and interpretation of study results. IH-P, SV, VA, CR-C, WS, MG, and LR-S performed final analyses and co-drafted the manuscript. All authors read, revised, and approved the final manuscript.
Funding
The authors received financial support for the research from Centro de Ciencias de la Atmósfera-Universidad Nacional Autónoma de México.
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Acknowledgments
Grateful acknowledgements are made to the Secretariat for the Environment (SEDEMA) of Mexico City for the public domain records.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fpubh.2021.642630/full#supplementary-material
Abbreviations
AF, Attributable fraction; AQHI, Air quality health index; ASI, Air Stagnation Index; AV, Amplitude value of diurnal cycle; CO, Carbon monoxide; COM, Commercial site; COVID-19, Coronavirus disease; ER, Excess risk; EPA, US Environmental Protection Agency; GPM, Global Precipitation Measurement; HB, Health burden; IND, Industrial site; MCMA, Mexico City Metropolitan Area; Molec, molecules; NO2, Nitrogen dioxide; NOX, Nitrogen oxides; O3, Ozone; OMI, Ozone Monitoring Instrument; PDF, Probability density function; PM2.5, Particulate matter with aerodynamic diameter <2.5 μm; PM10, Particulate matter with aerodynamic diameter <10 μm; RES, Residential site; RH, Relative humidity; RR, Relative risk; SD, Standard deviation; SIMAT, Atmospheric Monitoring System of the Mexico City government; SO2, Sulphur dioxide; T, Temperature; TRA, Traffic site; TROPOMI, Tropospheric Monitoring Instrument; UBN, Northern background site; UBS, Southern background site; VOC, Volatile organic compounds; WHO, World Health Organization; WD, Wind direction; WS, Wind speed.
References
1. Kharroubi S, Saleh F. Are lockdown measures effective against COVID-19? Front Public Health. (2020) 8:610. doi: 10.3389/fpubh.2020.549692
2. MxGob (Mexican Government) 2020 Covid-19 México Casos diarios por Estado + Nacional. Available online at: https://coronavirus.gob.mx/datos/#DownZCSV (accessed August 22, 2020).
3. SEP (Public Education Secretariat). 2020 ACUERDO Número 02/03/20 por el que se Suspenden las Clases en las Escuelas de Educación Preescolar Primaria Secundaria Normal y Demás Para la Formación de Maestros de Educación Básica del Sistema Educativo Nacional así Como Aquellas de los Tipos Medio Superior y Superior Dependientes de la Secretaría de Educación Pública. Available online at: https://www.dof.gob.mx/nota_detalle.php?codigo=5589479&fecha=16/03/2020 (accessed August 18, 2020).
4. CSG (CONSEJO DE SALUBRIDAD GENERAL). ACUERDO por el que se Declara Como Emergencia Sanitaria por Causa de Fuerza Mayor a la Epidemia de Enfermedad Generada por el Virus SARS-CoV2 (COVID-19). Diario Oficial de la Federación Lunes 30 de Marzo de 2020. (2020). Available online at: http://dof.gob.mx/2020/CSG/CSG_300320_VES.pdf (accessed August 22, 2020).
5. MxGob (Mexican Government) 2020 Comunicado Salud Inicia la fase 3 por COVID-19. Available online at: https://coronavirus.gob.mx/2020/04/21/inicia-la-fase-3-por-covid-19/ (accessed August 20, 2020).
6. Bao R, Zhang A. Does lockdown reduce air pollution? Evidence from 44 cities in northern China. Sci Total Environ. (2020) 31:139052. doi: 10.1016/j.scitotenv.2020.139052
7. Karuppasamy MB, Seshachalam S, Natesan U, Ayyamperumal R, Karuppannan S, Gopalakrishnan G, et al. Air pollution improvement and mortality rate during COVID-19 pandemic in India: global intersectional study. Air Qual Atmos Health. (2020) 13:1375–84. doi: 10.1007/s11869-020-00892-w
8. Mendez-Espinosa JF, Rojas NY, Vargas J, Pachón JE, Belalcazar LC, Ramírez O. Air quality variations in Northern South America during the COVID-19 lockdown. Sci Total Environ. (2020) 749:141621. doi: 10.1016/j.scitotenv.2020.141621
9. Shehzad K, Sarfraz M, Shah SGM. The impact of COVID-19 as a necessary evil on air pollution in India during the lockdown. Environ Poll. (2020) 266:115080. doi: 10.1016/j.envpol.2020.115080
10. Tobías A, Carnerero C, Reche C, Massagué J, Via M, Minguillón MC, Alastuey A, et al. Changes in air quality during the lockdown in Barcelona (Spain) one month into the SARS-CoV-2 epidemic. Sci Total Environ. (2020) 726:138540. doi: 10.1016/j.scitotenv.2020.138540
11. Cameletti M. The Effect of Corona Virus Lockdown on Air Pollution: evidence from the City of Brescia in Lombardia Region (Italy). Atmos Environ. (2020) 239:117794. doi: 10.1016/j.atmosenv.2020.117794
12. Petetin H, Bowdalo D, Soret A, Guevara M, Jorba O, Serradell K, et al. Meteorology-normalized impact of COVID-19 lockdown upon NO2 pollution in Spain. Atmos Chem Phys Discuss. (2020) 20:11119–41. doi: 10.5194/acp-2020-446
13. Rao ST, Zurbenko IG. Detecting and tracking changes in ozone air quality. Air Waste. (1994) 44:1089–92. doi: 10.1080/10473289.1994.10467303
14. Carslaw DC, Williams ML, Barratt B. A short-term intervention study—Impact of airport closure due to the eruption of Eyjafjallajökull on near-field air quality. Atmos Environ. (2012) 54:328–36. doi: 10.1016/j.atmosenv.2012.02.020
15. Grange SK, Carslaw DC. Using meteorological normalisation to detect interventions in air quality time series. Sci Total Environ. (2019) 653:578–88. doi: 10.1016/j.scitotenv.2018.10.344
16. Silva-Quiroz R, Rivera AL, Ordoñez P, Gay-Garcia C, Frank A. Atmospheric blockages as trigger of environmental contingencies in Mexico City. Heliyon. (2019) 5:e02099. doi: 10.1016/j.heliyon.2019.e02099
17. Horton DE, Diffenbaugh NS. Response of air stagnation frequency to anthropogenically enhanced radiative forcing. Environ Res Lett. (2012) 7:044034. doi: 10.1088/1748-9326/7/4/044034
18. Kerr GH, Waugh DW. Connections between summer air pollution and stagnation. Environ Res Lett. (2018) 13:084001. doi: 10.1088/1748-9326/aad2e2
19. INEGI (National Institute of Statistics and Geography). Población. (2015). Available online at: https://www.inegi.org.mx/temas/estructura/ (accessed August 10, 2020).
20. SEDEMA. Monitoreo de Contaminantes Criterio. Secretaría del Medio Ambiente de la Ciudad de México. (2020). Available online at: http://www.aire.cdmx.gob.mx/ (accessed August 18, 2020).
21. Kerimray A, Baimatova N, Ibragimova OP, Bukenov B, Kenessov B, Plotitsyn P, et al. Assessing air quality changes in large cities during COVID-19 lockdown: the impacts of traffic-free urban conditions in Almaty Kazakhstan. Sci Total Environ. (2020) 730:139179. doi: 10.1016/j.scitotenv.2020.139179
22. R: A language and environment for statistical computing. R Foundation for Statistical Computing. Vienna. Available online at: https://www.R-project.org/ (accessed January 15, 2021).
23. Sun Y, Solomon S, Dai A, Portmann RW. How often does it rain? J Climate. (2006) 19:916–34. doi: 10.1175/JCLI3672.1
24. Huffman GJEF, Stocker DT, Bolvin EJ, Nelkin-Jackson T. GPM IMERG Late Precipitation L3 1 Day 0.1 Degree x 0.1 Degree V06 Edited by Andrey Savtchenko Greenbelt MD Goddard Earth Sciences Data and Information Services Center (GES DISC) (2019). Available online at: https://disc.gsfc.nasa.gov/datasets/GPM_3IMERGDF_06/summary
25. Krotkov NA, Lamsal LN, Celarier EA, Swartz WH, Marchenko SV, Bucsela EJ, et al. The version 3 OMI NO2 standard product. Atmos Meas Tech. (2017) 10:3133–49. doi: 10.5194/amt-10-3133-2017
26. Krotkov NA, Lamsal LN, Marchenko SV, Celarier EA, Bucsela EJ, Swartz WH, et al. OMI/Aura Nitrogen Dioxide (NO2) Total and Tropospheric Column 1-orbit L2 Swath 13x24 km V003, Greenbelt MD USA, Goddard Earth Sciences Data and Information Services Center (GES DISC). (2019). Available online at: https://earthdata.nasa.gov/ (accessed August 15, 2020).
27. Rivera C, Stremme W, Grutter M. Nitrogen dioxide DOAS measurements from ground and space: comparison of zenith scattered sunlight ground-based measurements and OMI data in Central Mexico. Atmósfera. (2013) 26:401–14. doi: 10.1016/S0187-6236(13)71085-3
28. Borsdorff T, Hu H, Hasekamp O, Sussmann R, Rettinger M, Hase F, et al. Mapping carbon monoxide pollution from space down to city scales with daily global coverage. Atmos Meas Tech. (2018) 11:5507–18. doi: 10.5194/amt-11-5507-2018
29. Borsdorff T, García Reynoso A, Maldonado G, Mar-Morales B, Stremme W, Grutter M, et al. Monitoring CO emissions of the metropolis Mexico City using TROPOMI CO observations. Atmos Chem Phys. (2020) 20:15761–74. doi: 10.5194/acp-20-15761-2020
30. UNIATMOS (Repositorio Institucional Centro de Ciencias de la Atmósfera UNAM). Technical Information. (2020). Available online at: http://ri.atmosfera.unam.mx:8586/geonetwork/srv/spa/catalog.search#/metadata/adafdad6-44fc-4182-8dca-0d5d919e3d03 (accessed August 16, 2020).
31. Google. Google COVID-19 Community Mobility Report. (2020). Available online at: https://www.google.com/covid19/mobility/ (accessed May 20, 2020).
32. Gov-Canada (Government of Canada). About the Air Quality Health Index. (2020). Available online at: https://www.canada.ca/en/environment-climate-change/services/air-quality-health-index/about.html (accessed August 15, 2020).
33. Stieb DM, Burnett RT, Smith-Doiron M, Brion O, Shin HH, Economou V. A new multipollutant no-threshold air quality health index based on short-term associations observed in daily time-series analyses. J Air Waste Manag Assoc. (2008) 58:435–50. doi: 10.3155/1047-3289.58.3.435
34. Hu J, Ying Q, Wang Y, Zhang H. Characterizing multi-pollutant air pollution in China: Comparison of three air quality indices. Environ Int. (2015) 84:17–25. doi: 10.1016/j.envint.2015.06.014
35. World Health Organization. WHO Air Quality Guidelines for Particulate Matter, Ozone, Nitrogen Dioxide and Sulfur Dioxide: Global Update 2005: Summary of Risk Assessment. Occupational and Environmental Health Team (2006). Available online at: https://apps.who.int/iris/handle/10665/69477 (accessed January 20, 2021).
36. Sahu SK, Kota SH. Significance of PM2.5 air quality at the Indian capital. Aerosol Air Qual Res. (2017) 17:588–97. doi: 10.4209/aaqr.2016.06.0262
37. Cairncross EK, John J, Zunckel M. A novel air pollution index based on the relative risk of daily mortality associated with short-term exposure to common air pollutants. Atmos Environ. (2007) 41:8442–54. doi: 10.1016/j.atmosenv.2007.07.003
38. Chen Y, Wild O, Conibear L, Ran L, He J, Wang L, et al. Local characteristics of and exposure to fine particulate matter (PM2.5) in four indian megacities. Atmos Environ. (2020) 5:100052. doi: 10.1016/j.aeaoa.2019.100052
39. Kumar P, Hama S, Omidvarborna H, Sharma A, Sahani J, Abhijith KV, et al. A Temporary reduction in fine particulate matter due to ‘anthropogenic emissions switch-off’during COVID-19 lockdown in Indian cities. Sustain Cities Soc. (2020) 62:102382. doi: 10.1016/j.scs.2020.102382
40. GBD. Global Burden of Disease Study 2017 Data Resources: GHDx. (2017). Available online at: http://ghdx.healthdata.org/gbd-2017 (accessed January 20, 2021).
41. SEDEMA. Historical Analysis of Population Health Benefits Associated With Air Quality in Mexico City During 1990-2015. Secretaría del Medio Ambiente de la Ciudad de México; Ciudad de México (2018). Available online at: https://www.google.com/url?sa=t&rct=j&q=&esrc=s&source=web&cd=&ved=2ahUKEwjaiozcq8PuAhVKEawKHQ8TAdwQFjAAegQICBAC&url=http%3A%2F%2Fwww.data.sedema.cdmx.gob.mx%2Fbeneficios-en-salud-por-la-mejora-de-la-calidad-del-aire%2Fdescargas%2Fanalisis-ingles.pdf&usg=AOvVaw3n0lDfsTjYEQjKcQhPuZpd (accessed January 21, 2021).
42. Molina LT, Velasco E, Retama A, Zavala M. Experience from integrated air quality management in the Mexico City metropolitan area and Singapore. Atmosphere. (2019) 10:512. doi: 10.3390/atmos10090512
43. Anandharajan TRV, Vignajeth KK, Hariharan GA, Jijendiran R. Identification of outliers in pollution concentration levels using anomaly detection. In: 2016 International Conference on Computational Techniques in Information and Communication Technologies (ICCTICT). New Delhi (2016). p. 433–8. doi: 10.1109/ICCTICT.2016.7514620
44. Venter ZS, Aunan K, Chowdhury S, Lelieveld J. COVID-19 lockdowns cause global air pollution declines. Proc Natl Acad Sci. (2020) 17:18984–90. doi: 10.1073/pnas.2006853117
45. Sharma S, Zhang M, Gao J, Zhang H, Kota SH. Effect of restricted emissions during COVID-19 on air quality in India. Sci Total Environ. (2020) 728:138878. doi: 10.1016/j.scitotenv.2020.138878
46. Wang X, Zhang L, Moran MD. Bulk or modal parameterizations for below-cloud scavenging of fine, coarse, and giant particles by both rain and snow. J Adv Model Earth Syst. (2014) 6:1301–10. doi: 10.1002/2014MS000392
47. Hossain KM, Easa SM. Pollutant dispersion characteristics in Dhaka city Bangladesh. Asia Pac J Atmos Sci. (2012) 48:35–41. doi: 10.1007/s13143-012-0004-8
48. SEDEMA. Inventario de Emisiones de la Ciudad de México 2016 Contaminantes Criterio Tóxicos y Compuestos de Efecto Invernadero. Secretaría del Medio Ambiente de la Ciudad de México (2018). Available online at: https://www.dof.gob.mx/nota_detalle.php?codigo=5589479&fecha=16/03/2020 (accessed August 18, 2020).
49. SEDEMA. Calidad del Aire en la Ciudad de México Informe 2017. Secretaría del Medio Ambiente de la Ciudad de México; Ciudad de México (2018). Available online at: http://www.aire.cdmx.gob.mx/descargas/publicaciones/flippingbook/inventario-emisiones2016/mobile/ (accessed August 19, 2020).
50. Hernández-Paniagua IY, Clemitshaw KC, Mendoza A. Observed trends in ground-level O3 in Monterrey Mexico during 1993–2014: comparison with Mexico City and Guadalajara. Atmos Chem Phys. (2017) 17:9163. doi: 10.5194/acp-17-9163-2017
51. Stremme W, Grutter M, Rivera C, Bezanilla A, Garcia AR, Ortega I, et al. Top-down estimation of carbon monoxide emissions from the Mexico Megacity based on FTIR measurements from ground and space. Atmos Chem Phys. (2013) 13:1357–76. doi: 10.5194/acp-13-1357-2013
52. SEDEMA. Calidad del Aire en la ZMVM y Pandemia COVID-19. Secretaría del Medio Ambiente de la Ciudad de México (2020). Available online at: https://www.sedema.cdmx.gob.mx/comunicacion/nota/calidad-del-aire-en-la-zmvm-y-pandemia-covid-19 (accessed August 20, 2020).
53. SIE. Sistema de Información Energética, Secretaría de Energía SENER. (2020). Available online at: http://sie.energia.gob.mx/bdiController.do?action=cuadro&cvecua=PMXE2C03 (accessed August 20, 2020).
54. Song J, Lei W, Bei N, Zavala M, De Foy B, Volkamer R, et al. Ozone response to emission changes: a modeling study during the MCMA-2006/MILAGRO Campaign. Atmos Chem Phys. (2010) 10:3827–46. doi: 10.5194/acp-10-3827-2010
55. Kanda I, Basaldud R, Magaña M, Retama A, Kubo R, Wakamatsu S. Comparison of ozone production regimes between two Mexican cities: Guadalajara and Mexico City. Atmosphere. (2016) 7:91. doi: 10.3390/atmos7070091
56. García-Reynoso A, Jazcilevich A, Ruiz-Suárez LG, Torres-Jardón R, Suárez-Lastra M, Reséndiz Juárez NA. Ozone weekend effect analysis in México City. Atmósfera. (2009) 22:281–97. Available online at: https://www.revistascca.unam.mx/atm/index.php/atm/article/view/8627
57. Hernández-Paniagua IY, Lopez-Farias R, Piña-Mondragón JJ, Pichardo-Corpus JA, Delgadillo-Ruiz O, Flores-Torres A, et al. Increasing weekend effect in ground-level O3 in metropolitan areas of mexico during 1988–2016. Sustainability. (2018) 10:3330. doi: 10.3390/su10093330
58. García-Escalante JS, García-Reynoso JA, Jazcilevich-Diamant A, Ruiz-Suárez LG. The influence of the Tula Hidalgo complex on the air quality of the Mexico City Metropolitan Area. Atmósfera. (2014) 27:215–25. doi: 10.1016/S0187-6236(14)71111-7
59. Muhammad S, Long X, Salman M. COVID-19 pandemic and environmental pollution: a blessing in disguise? Sci Total Environ. (2020) 728:138820. doi: 10.1016/j.scitotenv.2020.138820
60. Nakada LYK, Urban RC. COVID-19 pandemic: impacts on the air quality during the partial lockdown in São Paulo state Brazil. Sci Total Environ. (2020) 139087. doi: 10.1016/j.scitotenv.2020.139087
61. Mahato S, Pal S, Ghosh KG. Effect of lockdown amid COVID-19 pandemic on air quality of the megacity Delhi India. Sci Total Environ. (2020) 730:139086. doi: 10.1016/j.scitotenv.2020.139086
62. CDMX-Gob. New Restriction Measures During the COVID-19 Phase 3. Mexico City Local Government (2020). Available online at: https://www.cdmx.gob.mx/portal/articulo/nuevas-medidas-para-la-fase-3-por-covid-9 (accessed August 20, 2020).
63. Dantas G, Siciliano B, França BB, da Silva CM, Arbilla G. The impact of COVID-19 partial lockdown on the air quality of the city of Rio de Janeiro Brazil. Sci Total Environ. (2020) 729:139085. doi: 10.1016/j.scitotenv.2020.139085
64. Wang P, Chen K, Zhu S, Wang P, Zhang H. Severe air pollution events not avoided by reduced anthropogenic activities during COVID-19 outbreak. Resour Conserva Recycl. (2020) 158:104814. doi: 10.1016/j.resconrec.2020.104814
65. SEDEMA. Se Activa Contingencia Ambiental Atmosférica Extraordinaria por Partículas (PM2.5) y Ozono en la ZMVM. Secretaría del Medio Ambiente de la Ciudad de México (2019). Available online at: http://189.240.89.22/boletines-contingencia/14-05-2019/comunicado2.pdf (accessed August 10, 2020).
Keywords: COVID-19, Fourier series, meteorology, remote sensing, air pollutant emissions
Citation: Hernández-Paniagua IY, Valdez SI, Almanza V, Rivera-Cárdenas C, Grutter M, Stremme W, García-Reynoso A and Ruiz-Suárez LG (2021) Impact of the COVID-19 Lockdown on Air Quality and Resulting Public Health Benefits in the Mexico City Metropolitan Area. Front. Public Health 9:642630. doi: 10.3389/fpubh.2021.642630
Received: 16 December 2020; Accepted: 25 February 2021;
Published: 25 March 2021.
Edited by:
Jennifer Ann Salmond, The University of Auckland, New ZealandReviewed by:
Gulnihal Ozbay, Delaware State University, United StatesBehzad Heibati, University of Oulu, Finland
Copyright © 2021 Hernández-Paniagua, Valdez, Almanza, Rivera-Cárdenas, Grutter, Stremme, García-Reynoso and Ruiz-Suárez. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Iván Y. Hernández-Paniagua, ivan.hernandez@atmosfera.unam.mx