{'fraction_non_alphanumeric': 0.033220782049885635, 'fraction_numerical': 0.029789783248012198, 'mean_word_length': 4.767116834170854, 'pattern_counts': {'":': 0, '<': 2, '<?xml version=': 0, '>': 2, 'https://': 0, 'lorem ipsum': 0, 'www.': 0, 'xml': 0}, 'pii_count': 0, 'substrings_counts': 0, 'word_list_counts': {'cursed_substrings.json': 0, 'profanity_word_list.json': 0, 'sexual_word_list.json': 0, 'zh_pornsignals.json': 0}}
Air pollution exposure and COVID-19: A look at mortality in Mexico City using individual-level data

We use individual-level data to estimate the effects of long- and short-term exposure to air pollution (PM2.5) on the probability of dying from COVID-19. To the best of our knowledge, our study is the first to look at this relationship using individual-level data. We find that for Mexico City there is evidence of a positive relationship between pollution and mortality that significantly grows with age and that appears to be mostly driven by long- rather than short-term exposure. By using a rich set of individual- and municipal-level covariates we are able to isolate the effect of exposure to pollution from other crucial factors, thus alleviating endogeneity concerns related to selection. Our results provide yet another reason for the need to implement environmental strategies that will reduce the exposure to air pollution: it is a key element to improve the general population's health. In addition, and considering that at this moment we do not know when the pandemic will stop or if SARS-CoV-2 will become a recurrent threat, the relationship that we uncovered suggests that financial resources should be allocated to improve medical services in those areas where PM2.5 concentrations tend to be high.


Introduction
The COVID-19 pandemic has had devastating effects both in terms of human lives and health, and of economic and social costs. Given this challenging context, understanding the factors that might increase mortality risks from COVID-19 is essential so that policy makers can design and implement effective policies to protect individuals with high-risk characteristics and living in contexts that might enhance their vulnerability. While the global pandemic has spread across all regions, the death toll it has brought along has varied enormously between countries, plausibly due to a combination of population characteristics (i.e. age structure, comorbidity prevalence), public policies, and environmental factors. In this study, we focus on the impact of one important environmental factor, air pollution, which has been explored recently in academic studies as a factor contributing to COVID-19 mortality (Conticini et al., 2020;Copat et al., 2020). We analyze the case of Mexico which, as of October 28th, ranks fourth in the world in terms of total deaths (>90,000) and is among the 10 countries with the highest number of deaths per 100 thousand inhabitants (>70) (JH-CSSE, 2020). COVID-19 deaths are largely concentrated in the Mexico City Metropolitan Area (MCMA), where air pollution levels, in spite of no longer being at the historically high levels observed during the end of the last century, remain very high.
Polluted air gains access to the body through the respiratory tract but it has systemic effects that can damage many other organs (Schraufnagel et al., 2019). Therefore, air pollution is linked to many health problems and can lead to premature death in children and adults (Brunekreef and Holgate, 2002;Cohen et al., 2017;Burnett et al., 2018). In particular, air pollutants can cause a persistent inflammatory response and thus increase the risk of infection from viruses that target the respiratory tract (Brunekreef and Holgate, 2002;Travaglio et al., 2020). Hence, it has been argued that prolonged exposure to air pollution makes people more vulnerable to COVID-19 (Yao et al., 2020;Wu et al., 2020). Furthermore, although the role of air pollution in the airborne transmission of the virus is still uncertain, preliminary evidence confirms that SARS-CoV-2 might be present on particulate matter (Setti et al., 2020). If that is the case, current exposure to air pollution could contribute to the spread of the virus as it is potentially carried deep into the lungs on the surface of particulate matter (Schraufnagel et al., 2019). Results for China seem to confirm this as they show a positive relationship between ambient air pollutants and COVID-19 infection (Zhu et al., 2020).
Air pollution is a pressing issue for many megacities, in particular for those located in developing countries (Marlier et al., 2016). The MCMA, the fifth largest in the world (United Nations, 2018), is no exception as its pollution levels frequently exceed the limits suggested by the World Health Organization (WHO) (SEDEMA, 2017). Furthermore, the MCMA has been hit hard by the COVID-19 pandemic: as of October 28th, there have been more than 21,000 COVID-19 deaths in the MCMA-only 13 countries in the world have had more deaths.
Recent research has used aggregate data (e.g., county, municipal, etc.) to look at the relationship between air pollution and COVID-19 (e.g., Cole et al., 2020;Coker et al., 2020, andWu et al., 2020). To the best of our knowledge, our study is the first to look at the relationship between air pollution and COVID-19 mortality at the individual level. By doing so, we are able to isolate the effect of exposure to air pollution from other crucial factors like age, sex, obesity, diabetes, and hypertension, thus partially alleviating endogeneity concerns. A second contribution of this paper is the use of both short-and long-term exposure to air pollution to disentangle the effects of each type of exposure on the probability of dying from COVID-19. Finally, the context of the countries analyzed in the available studies, with the exception of Vasquez-Apestegui et al. (2020), is very different from that of Mexico and other developing countries, in particular those in Latin America, where health systems have a limited capacity and where the share of the population without formal employment is high (López-Feldman et al., 2020). Our results can contribute to the design of public policies that are better focused to the realities of countries with such characteristics.

Data
Data for confirmed COVID-19 cases and deaths come from the official statistics released daily by Mexico's Ministry of Health (SEDESA, 2020; SSA, 2020). 1 Data for each individual is collected by medical personnel. The information collected and reported includes basic demographic characteristics, health variables (such as being diagnosed with diabetes, hypertension or obesity), smoking status, date when symptoms began, and whether the individual died. Additionally, individuals' municipality of residence is reported, thus allowing us to match each individual to municipal level covariates as well as different measures of pollution exposure.
The municipal level variables that we use in our analysis are population size, population density, hospital beds per capita, poverty measures ("health-service access deficiency" and "food insecurity"), and percentage of jobs that can be done at home. All the information comes directly from official sources (Coneval, 2017;INEGI, 2015) except the percentage of jobs that can be done from home which was estimated based on information from (INEGI, 2015) following Dingel and Neiman (2020).
We rely on two sources of information for air pollution, one for longterm exposure and the other for short-term exposure. The average annual concentrations of particulate matter finer than 2.5 micrograms per cubic meter (PM 2.5 μg/m 3 ) were obtained from Hammer et al. (2020). Their global dataset, estimated from NASA's satellite products combined with ground level observations from the WHO, has been widely used to analyze the impact of air pollution on health outcomes. Wu et al. (2020) and Knittel and Knittel and Ozaltun (2020) use their data to study the relationship between PM 2.5 exposure and COVID-19 deaths in the US. Our analogous municipality-averaged measure was obtained after averaging 1.1 × 1.1 km 2 grid-cells (0.01 × 0.01 degrees) within each municipality each year. The mean exposure between 2000 and 2018 is shown in Fig. 1; although all 76 municipalities comprised in the MCMA (which encompasses Mexico City) exceed the WHO's standard of annual 10 μg/m 3 , variability in exposure across the territory is comparable to other recent studies focusing in small countries or specific regions within a country (e.g., Cole et al., 2020;Coker et al., 2020). In the MCMA, the highest exposure is observed in central and northeastern municipalities close to the boundaries between Mexico City and the State of Mexico where most economic activity 1 We use the data update available on October 7th, 2020. For all our estimations we restrict the sample to those individuals for which the symptoms started at least 14 days before that date. occurs, and most people inhabit. Fig. 1 suggest that such areas are also experiencing high numbers of COVID-19 deaths. Our contemporaneous PM 2.5 exposure variables are estimated using information published by Mexico City's Secretary of Environment from ground-based air quality monitors (SEDEMA, 2020). Their Automatic Network of Atmospheric Monitoring comprises 34 monitors that report hourly measures of concentrations of different pollutants throughout the MCMA, but that are more densely distributed within Mexico City and bordering municipalities from the State of Mexico. In order to assign more accurately pollution concentrations to each municipality, we restrict our analysis to those municipalities whose centroid is at most 7 km away from a monitoring station. Therefore, given data availability and the geographic location of the air quality monitors, we can measure short-term exposure for 14 out of 16 of the MCMA municipalities in Mexico City using 18 out of the 21 stations that measure PM 2.5 . 2 With the information available, we calculate the weighted average of daily concentrations of PM 2.5 based on the readings from all monitors located within a 7 km radius from each municipality centroid (the average municipality had 2 monitoring stations located within this radius). 3 Weights are given by the inverse of the distance between the municipality centroid and each monitor. Using this information, we calculate alternative versions of our short-term exposure measure that we use as robustness checks.
We also calculate an alternative version of the long-term exposure variable based on ground-based pollution concentrations during 2019. Temporally, this variable could be considered a more accurate measure of the actual exposure of the average inhabitant of a given municipality if there is a concern of people changing residential locations from one municipality to another between 2018 and 2020. Spatially, however, the satellite-based measures are constructed from finer-grid cells within a municipality and thus could be considered to more accurately represent exposure of inhabitants in a given municipality.
It is worth noting that, in addition to the substantial geographic variation in PM 2.5 concentrations within the MCMA that is shown in Fig. 1, we also observe ample variation in the exposure to air pollution across individuals with different characteristics and comorbidities. In particular, Fig. 2 shows variation in PM 2.5 long-term exposure across age, diagnoses for diabetes, hypertension, obesity, and smoking status for individuals living in the MCMA. These characteristics and comorbidities have been found to be important modifiers of the risk of dying from COVID-19; it is therefore particularly important to have variation in exposure to air pollution across them.
As Table 1 shows, the MCMA and the Mexico City sub-sample are very similar in terms of individual characteristics, with slightly less than half of the cases being female, a mean age of 45, and comorbidity rates between 10% and 18%. In terms of pollution exposure and socioeconomic characteristics, the Mexico City sub-sample lives in slightly denser and more polluted areas but has larger access to hospital beds, health care, food security, and telework. Importantly, the death rate among confirmed cases is respectively 10% and 11%. Worth mentioning is that in Mexico the number of total COVID-19 tests performed per capita is very low (Roser et al., 2020). This implies that the number of COVID-19 deaths is likely to be underestimated. The low level of testing is the result of an explicit decision made by the federal government to not perform widespread testing. This policy is being applied in the whole country and there is no evidence that suggests that testing levels vary inside the MCMA. 4 Therefore, we argue that underestimation is not correlated to municipal pollution levels and hence should not bias our estimations in any meaningful way. 5 As is the case in many other countries, there are instances where deaths are not immediately attributed to COVID-19 due to lack of testing. These deaths remain classified as suspicious until a committee decides if they should be attributed to COVID-19. Our estimations use the most up-to-date information, which includes deaths that were originally classified as suspicious and later on were attributed to COVID-19. There is no evidence to suggest that in the MCMA the number of deaths classified as suspicious varies with pollution levels.

Econometric specification
We specify the following probit model to estimate the effect that pollution has on the probability of dying from COVID-19: where Φ is the standard normal cumulative distribution function, Death im is a dummy variable equal to one if individual i from municipality m who has been infected by COVID-19 dies and zero otherwise, pollution m is a measure of air pollution at municipality m, and X im is a vector of individual-and municipal-level variables that control for potential confounding factors. 6 Based on the literature on air pollution and health, and on the results available for the relationship between exposure to PM 2.5 and COVID-19, we expect β to be positive. In addition to testing this hypothesis, we are interested in estimating the effect that an increase of 1 μg/m 3 in PM 2.5 has on the probability that an individual will die after having contracted COVID-19. To do so, we take the derivative of the right-hand side of Eq. (1) with respect to pollution and estimate the average marginal effect of PM 2.5 pollution on the probability of dying (for more on the average marginal effect of the probit model see Wooldridge (2002), p. 467).
We are also interested in the effects that short-term exposure might have on the probability of dying from COVID-19. We use the air pollution data that is available for Mexico City to analyze this. This allows us to estimate a modified version of Eq. (1) where, in addition to longterm exposure, we include additional measures of exposure to PM 2.5 . In particular, in addition to our long-term exposure variable, we include a variable that captures the cumulative average concentration of PM 2.5 during the 14 days before the onset of symptoms (PM 2.5 -14). Additionally, in some of the models we use the variable that measures average concentrations in Mexico City during 2019 (PM 2.5 -2019).
Addressing endogeneity concerns is crucial in order to establish a causal relationship between exposure to pollution and health. In particular, the evidence suggests that environmental quality is an important dimension for choosing residential location, which is likely to lead to a correlation between higher pollution levels and both observable and unobservable characteristics correlated with health-potentially leading to omitted variable bias. Additionally, avoidance behavior may similarly 2 Incorporating the two remaining municipalities in Mexico City, representing 6% percent of the city's population, required the maximum distance from a municipality to a monitor to be increased to 15 km. This distance is particularly unlikely to provide reliable pollution exposure measures for these two municipalities given that the weightedaverage would be estimated based on a single and relatively remote monitoring station. 3 We also restricted our measure to consider only readings from 7:00 am-7:00 pm which we consider to better represent daily exposure due to economic and social activities that mostly takes place during those hours. Nonetheless, all the estimations presented in the paper were also done using 24-hour readings, the results (available upon request) are very similar to those obtained using 12-hour readings. 4 This is different to what happens in other countries, like Italy, where testing policies vary among regions and therefore measurement error is non-randomly distributed (Coker et al., 2020). 5 In order to empirically test this, we regress the share of positive COVID-19 tests at the municipality-month level on a linear time trend and municipal-level covariates (population size, population density, hospital beds per capita, poverty measures, and percentage of jobs that can be done from home) and obtain residuals. We then regress the residuals on the mean PM 2.5 municipal variable and find no statistically significant associations (95% CI: −0.006, 0.002) between residuals and air pollution, suggesting that undermeasurement is indeed uncorrelated to PM 2.5 levels. Intuitively, this shows that being above or below the general trend for MCMA in terms of the share of positive tests is uncorrelated with pollution levels conditional on municipal-level covariates. This relation can be seen graphically in Figure A1 of the Appendix. 6 In some of the estimations for Mexico City we take advantage of the fact that for that sub-sample we have information at the locality level (a smaller administrative area than a municipality) and hence use pollution concentrations at that level.
A. López-Feldman, D. Heres and F. Marquez-Padilla Science of the Total Environment 756 (2021) 143929 lead to nonrandom assignment of pollution exposure (Currie et al., 2014). In order to partially address these concerns, we include a rich set of municipal-and individual-level covariates that we use to control for factors that may be correlated to pollution exposure and simultaneously affect mortality. In our preferred specification we include population size, population density, hospital beds per capita, poverty measures ("health-service access deficiency" and "food insecurity"), and percentage of jobs that can be done at home, all measured at the municipal level, as covariates. We also include individual indicators for sex, obesity, smoking status, diabetes, hypertension, age, age squared, and the day in which symptoms started to control for time effects. The inclusion of these covariates is likely to alleviate endogeneity concerns -Cole et al. (2020) and Wu et al. (2020)  COVID-19. Our standard errors are clustered at the municipality level in order to correct for spatial correlation of the individuals' unobserved variation that may remain after including covariates.

Results
To begin the econometric analysis, we use the full sample for the MCMA and estimate a version of Eq. (1) that only includes long-term exposure to PM 2.5 as the explanatory variable. Our results (Model 1, Table 2) show that there is a positive and statistically significant relationship between exposure and the probability of dying from COVID-19. The relationship remains valid when we include municipal level covariates in Model 2. The last column (Model 3) shows the results of including individual level covariates in our analysis and is our preferred specification. The positive relationship between long-term exposure to PM 2.5 and the probability of dying from COVID-19 is robust to the inclusion of covariates at both the individual and municipal levels. 7 In order to better interpret our results, we estimate marginal effects for the three versions of the model, which are shown in Fig. 3. For the version without covariates (Model 1, Table 2) we estimate an average marginal effect of 0.0076, which implies that for an increase of 1 μg/m 3 in PM 2.5 the probability of dying from COVID-19 increases on average by 0.76 percentage points. For the model with municipal-level covariates (Model 2) the marginal effect is 0.0129. Our preferred specification is the one that takes into account age, obesity, diabetes, hypertension, smoking, and municipal level covariates, as in Model 3. For that model we estimate that on average the probability of dying from COVID-19 increases by 0.77 percentage points with an increase of 1 μg/m 3 in PM 2.5 . Given that the average mortality rate for our sample is of approximately 10.4%, our estimates suggest that an increase of 1 μg/m 3 in PM 2.5 may increase mortality risk by approximately 7.4%. Fig. 4 presents the average marginal effects of pollution across age. It shows that the effects of long-term exposure on the probability of dying are heterogeneous and increasing with age. Long-term exposure to air pollution has a very small effect on children (0.02 percentage points from one additional μg/m 3 of PM 2.5 for 5-year-old children). At first, the magnitude of the effect increases slowly-individuals with 20 years of age have a marginal effect of 0.09 percentage points-but starts increasing faster gradually-the marginal effect is 0.45 percentage points for individuals that are 40 years old and 1.23 percentages points for those that are 60 years old, until it plateaus around 2 percentage points for individuals aged 80 years or older. 8 Table 3 shows the results of using the data from air pollution monitors that is available for Mexico City. With these data we are able to, in Table 2 The relationship between long-term exposure to pollution and the probability of dying from COVID-19 in the MCMA. The relationship is also robust to the way in which the probability function is modeled; our results remain valid when we estimate a linear probability model instead of a probit model (Table A1 of the Appendix). 8 We do not find evidence that the magnitude of the effect of air pollution on COVID-19 mortality risk varies according to comorbidities (hypertension, type 2 diabetes, obesity) or smoking status-although these factors do affect the probability of death from COVID-19 directly ( Figure A2 of the Appendix). Sources: Own estimation with data from Coneval (2017), SEDESA (2020), SSA (2020), Dingel and Neiman (2020), INEGI (2015), SEDEMA (2020) and Hammer et al. (2020). The MCMA is comprised of 76 municipalities in the states of Hidalgo, Mexico, and Mexico City. All municipalities (16) from Mexico City are part of the MCMA and are included in the MCMA sample. The Mexico City sample only considers the 14 municipalities for which we have PM 2.5 for 2019 and 2020. PM 2.5 -2019 and PM 2.5 -14 are calculated from the weighted average of the daily concentrations from 7:00 am to 7:00 pm based on all monitors located within a 7 km radius from each municipality centroid, with weights given by the inverse of the distance between the municipality centroid and each monitor.
A. López-Feldman, D. Heres and F. Marquez-Padilla Science of the Total Environment 756 (2021) 143929 addition to our long-term exposure variable, include short-term PM 2.5 cumulative average concentration during the 14 days before the onset of symptoms. Model 1 shows that the coefficient of short-term pollution is significant only at the 10% level. In Model 2, we use PM 2.5 concentrations during 2019 instead of using the long-term variable. We only find weak evidence to support that there is an effect of short-term exposure on the probability of dying (Model 2, Table 3). The effect of long-term exposure seems to be partially captured with the variable PM 2.5 -2019, which is weakly statistically significant. This is reassuring; both variables (PM 2.5 -2019 and PM 2.5 ) seem to be in fact capturing the same phenomenon. As a way to test that this is in fact the case, the last column shows the results of estimating the model with the three pollution variables simultaneously (Model 3). The two long-term exposure variables appear to be substitutes of each other, while once again we only find a slightly statistically significant effect for short-term exposure. 9 As a robustness check, we also estimate all the models presented in Table 3 using cumulative average concentrations during 2, 7, 5 and 10 9 All models presented in Table 3 were also estimated using linear probability models.
The results, shown in Table A2 of the Appendix, are very similar to those presented in Table 3. As an additional robustness check, we also included the average temperature that the individual experienced in the 14 days before the onset of symptoms as an additional variable in all the models presented in Table 3. Results for PM 2.5 and PM 2.5-2019 remain virtually identical, while PM 2.5 -14 ceases to be statistically significant (Table A3 of the  Appendix).  [2000][2001][2002][2003][2004][2005][2006][2007][2008][2009][2010][2011][2012][2013][2014][2015][2016][2017][2018] shown with 95% clusterrobust confidence intervals. The models behind each estimation are those described in Table 2.  Table 2) includes municipal-level covariates (population density, population, density of hospital beds, percentage of population without access to health care, percentage of population with moderate or severe food insecurity, and percentage of labor force with jobs that can be done from home), and individual-level covariates (gender, age and age squared, obesity, diabetes, hypertension, smoking status, and day in which symptoms started).

Table 3
The relationship between long-and short-term exposure to pollution and the probability on dying from COVID-19 in Mexico City. Notes: The dependent variable is a dummy equal to one if an individual diagnosed with COVID-19 dies and zero otherwise. Estimations are done using a probit model. Cluster robust standard errors at the municipal-level shown in brackets. Municipal-level covariates are: Population density, population, density of hospital beds, percentage of population without access to health care, percentage of population with moderate or severe food insecurity, and percentage of labor force with jobs that can be done from home. Individuallevel covariates are: gender, age and age squared, obesity, diabetes, hypertension, smoking status, and day in which symptoms started. The Mexico City sample only considers the 14 municipalities for which, in addition to long-term exposure to PM 2.5 , we know PM 2.5 for 2019 and 2020. ⁎ p < 0.1. ⁎⁎ p < 0.05.
A. López-Feldman, D. Heres and F. Marquez-Padilla Science of the Total Environment 756 (2021) 143929 days instead of 14 days (results are presented in Table A4 of the Appendix). While the results for long-term exposure are similar to those presented in Table 3, results for short-term exposure are not statistically significant for any of the alternative exposure windows. Finally, as shown in Fig. 5, the marginal effects for changes in long-term PM2.5 exposure estimated for Mexico City are very similar to those estimated for the MCMA (Fig. 3). Marginal effects for short-term exposure are very small and only weakly significant (90% CI: 0.00006, 0.00099).

Discussion
Given the magnitude of the health crisis caused by SARS-CoV-2 and the amount of research that it has spurred across many fields, many researchers have looked at the relationship between exposure to particulate matter and morbi-mortality by COVID-19, perhaps unsurprisingly. Many of these studies have analyzed the relationship between exposure and infection and have shown that there is in fact a positive association between high levels of PM 10 and PM 2.5 and confirmed cases (Coccia, 2020;Zhu et al., 2020). Consistent with results for the SARS 2003 outbreak (Cui et al., 2003), the research that looks at the effects on mortality has found, as we do, that there is a positive relationship between exposure to PM 2.5 and COVID-19 mortality.
Our results cannot be directly compared to other results for COVID-19 as ours is, to the best of our knowledge, the first study to estimate the effects of air pollution on the probability of dying at the individual level. 10 Nonetheless, a natural point of comparison are the results from Wu et al. (2020) for the U.S., as we use similar covariates in our analysis and the long-term exposure variable that we use comes from the same source, although they look at exposure from 2000 to 2016 while we use 2000 to 2018. Their reported estimate that a 1 μg/m 3 increase in PM 2.5 is associated to an 8% increase in the death rate, is consistent with our estimation of an increase in the COVID-19 mortality rate of 0.77 percentage points, which is equivalent to an increase in the probability of dying from COVID-19 of 7.4%. Our results are in line with those published by Petroni et al. (2020) who find that a 1 μg/m 3 increase in PM 2.5 is correlated with a 7% increase in the mortality rate ratio in the U.S. Cole et al. (2020) and Coker et al. (2020) also look at long-term exposure to PM 2.5 and mortality. They use municipal-level data for Netherlands and Italy, respectively, and find evidence of a positive relationship. Coker et al. (2020) use excess mortality as dependent variable and estimate that a one-unit increase in PM 2.5 concentration (μg/m 3 ) is associated with a 9% increase in mortality. Meanwhile, Cole et al. (2020) find that the same increase in PM 2.5 will lead to 2.3 more deaths at the municipality level. Although not directly comparable, our results are similar in magnitude to those form Coker et al. (2020) and higher than those from Cole et al. (2020). This is consistent with the fact that Mexico's per-capita COVID-19 mortality is similar to Italy's and much higher to that of the Netherlands.
The work of Yao et al. (2020) for 49 cities in China is also related to ours as it uses short-term and long-term exposure indicators, average PM concentrations from January 15, 2020 to February 29, 2020 and average concentrations of PM annual mean levels in years 2015-2019. Using data at the city-level they show, like we do, that the association between short-and long-term exposure and fatality is positive.
Contrary to all these studies that use data at an aggregate level, by using individual-level information we are able to estimate how exposure to PM 2.5 can have a differential effect on the probability of dying according to the age of the infected individual. We show that the effect of PM 2.5 is monotonically increasing with age and reaches a maximum at around 80 years of age, a finding that is consistent with the available evidence that clearly shows that COVID-19 fatality rates vary with age (Roser et al., 2020). Furthermore, although the effects are small, according to our results exposure to PM 2.5 can actually increase the probability that infants and small children die from COVID-19. This is consistent with evidence showing that exposure to PM 2.5 can lead to systemic inflammation and suppress early immune responses to infections even in young people (Conticini et al., 2020;Wu et al., 2020).
This study is subject to some limitations. First, although we include a rich set of municipal-and individual-level covariates, available information does not allow us to directly control for individuals' residential choice or avoidance behavior. Both factors could lead to nonrandom assignment of pollution exposure and hence represent a threat to establishing a causal relationship. Second, we are implicitly assuming that  [2000][2001][2002][2003][2004][2005][2006][2007][2008][2009][2010][2011][2012][2013][2014][2015][2016][2017][2018] shown with 95% clusterrobust confidence intervals. The models behind each estimation are those described in Table 3. 10 Travaglio et al. (2021) use individual-level data for England but they estimate the probability of being infected, not the probability of dying.
A. López-Feldman, D. Heres and F. Marquez-Padilla Science of the Total Environment 756 (2021) 143929 air pollution is evenly distributed within each municipality (locality). If that is not the case and the distribution of PM 2.5 inside a municipality (locality) is correlated to some individual characteristics that also affect the probability of dying, then our results could be biased. Third, although we know the residence of all the individuals in our sample, we do not know for how long they have been living there or how much of their time they actually spend at home. In our estimations, individuals' long-term exposure to PM 2.5 is based on residence. Nonetheless, if an individual recently moved to a municipality, their true exposure to air pollution would be different to the one that we are attributing. Although we partially control for this by using exposure during 2019 as an alternative measure to exposure during 2000-2018, it is still possible for our results to be biased because of this. Fourth, we do not find robust evidence of an effect of short-term PM 2.5 exposure on the probability of dying from COVID-19. We had expected to find a statistically significant relationship because pollution exposure could be expected to generate an acute deterioration of the immunological system and to facilitate transportation of the virus. However, due to heterogeneous changes in daily activities and mobility across the population during the pandemic, measured PM 2.5 concentrations may be less indicative of actual exposure and thus become a bad predictor of the probability of dying from COVID-19. In spite of these limitations, our study makes relevant contributions to the literature.

Concluding remarks
We find evidence of a positive relationship between PM 2.5 air pollution and the probability that an individual will die after contracting COVID-19; this relationship increases with age, especially for individuals that are 40 years old or more. Although our results suggest that this is mainly driven by long-term exposure, we cannot rule out that short-term exposure also has an effect. These results are robust to model specification and inclusion of confounders at the municipaland individual-level. The last is the most important contribution of this paper as, in contrast to the previous literature, we are able to explicitly include in our estimations not only the age and gender of the person affected by COVID-19 but also if she has underlying health conditions. By doing this, we reduce the likelihood that unobservables could explain away the effects that we report. In spite of these efforts, our robust statistical evidence cannot be fully interpreted as causal as we do not have experimental or quasi-experimental data. Furthermore, although in our estimations we use cluster-robust standard errors, there might be spatial autocorrelation issues that cannot be fully accounted for by our empirical strategy that is based on individual-level data.
Mexico is still in the midst of the pandemic; health conditions are continuously evolving, and data related to COVID-19 cases and mortality is being updated regularly. In the coming months we should have better data and a cleaner picture of the full extent of the effects of the pandemic both in the MCMA and in the country. Meanwhile, our results provide yet another reason for the need to implement more stringent controls to air pollution. In addition, and considering that at this moment we do not know when the pandemic will stop or if SARS-CoV-2 will become a recurrent threat, the relationship that we uncovered suggests that financial resources should be allocated to improve medical services in those areas where PM 2.5 concentrations tend to be high; long-term exposure has made people living in those areas more vulnerable to the disease. Additionally, policy makers and governments should rethink and reinforce environmental strategies that will reduce levels of air pollution as a key element of public health policies.
In our analysis we focus on the MCMA because of its high levels of air pollution and high rate of COVID-19 mortality. Nonetheless, it is not clear if our results are valid for other big metropolitan areas in the world or even for other parts of Mexico. The fact that our results are in general consistent with other studies on the relation between air pollution and COVID-19 mortality gives us greater confidence that our findings are externally valid-at least to some extent. Finally, our findings suggest that there are groups that, given their exposure to PM 2.5 , could be particularly vulnerable to COVID-19 in developing countries with characteristics similar to Mexico. Notorious among these groups are users of public transportation, street vendors, and households that use solid fuels for all or part of their cooking or heating needs. Understanding how COVID-19 affects these groups will be essential for the design of public policies that can successfully decrease the toll that this pandemic is exerting on the population.