Understanding the Pathways from Prenatal and Post-Birth PM2.5 Exposure to Infant Death: An Observational Analysis Using US Vital Records (2011-2013).

We studied the relationship of prenatal and post-birth exposure to particulate matter < 2.5 μm in diameter (PM2.5) with infant mortality for all births between 2011 and 2013 in the conterminous United States. Prenatal exposure was defined separately for each trimester, post-birth exposure was defined in the 12 months following the prenatal period, and infant mortality was defined as death in the first year of life. For the analysis, we merged over 10 million cohort-linked live birth-infant death records with daily, county-level PM2.5 concentration data and then fit a Structural Equation Model controlling for several individual- and county-level confounders. We estimated direct paths from the two exposures to infant death as well as indirect paths from the prenatal exposure to the outcome through preterm birth and low birth weight. Prenatal PM2.5 exposure was positively associated with infant death across all trimesters, although the relationship was strongest in the third trimester. The direct pathway from the prenatal exposure to the outcome accounted for most of this association. Estimates for the post-birth PM2.5-infant death relationship were less precise. The results from our study add to a growing literature that provides evidence in favor of the potential harmful effects on human health of low levels of air pollution.

1. Introduction

Although several studies have investigated the relationship between early life ambient air pollution exposure and infant mortality, the evidence about the importance of prenatal versus post-birth exposure is mixed [1,2]. Some studies have found that exposure to higher levels of post-birth air pollution increases the risk of infant death but increased prenatal air pollution exposure does not [3,4]. In contrast, other studies have found that higher levels of prenatal air pollution exposure increases the risk of infant death, but post-birth air pollution exposure does not [5,6]. Furthermore, some studies have also suggested that both prenatal and post-birth air pollution exposure may increase the risk of infant mortality [7,8].

There are plausible biological mechanisms to link both prenatal and post-birth air pollution exposure to infant death. For instance, there is very strong evidence linking prenatal air pollution exposure with increased risk of adverse birth outcomes such as preterm birth and low birth weight [9]. Such a relationship would suggest that prenatal air pollution may be linked to infant death because prematurity and low birth weight are two of the most important drivers of infant mortality in the US [10]. At the same time, there are biologically plausible direct pathways from prenatal air pollution exposure to infant mortality as well: for example, exposure to fine particulate matter can cause imbalances in an individual’s autonomic nervous system and result in oxidative stress, both of which when experienced by pregnant people can affect maternal and fetal health, which in turn can have implications for infant health [11,12]. Fine particulate matter can also traverse the placental barrier, thus directly affecting the health of the developing fetus, which in turn could lead to worse infant health outcomes and, subsequently, infant death [13]. In terms of post-birth exposure, infants who are exposed to high levels of air pollution may demonstrate similar pathophysiological responses as adults which, in combination with their more immature immune and lung systems, could increase their risk of morbidity and death as well [14,15].

Given the biological plausibility, the mixed nature of the evidence on the relationship of prenatal and post-birth air pollution exposure with infant death may be due to several reasons. First, studies are often conducted in different geographic settings and over different time periods, which have important implications for the type and level of pollution exposure. Second, the exposure may be measured differently across different studies, with some using modeled estimates, while others use distance-weighted estimates of pollution directly from monitors. Finally, analytic strategies may be different across different studies, which could potentially affect the results. For example, variables on the causal pathway, such as gestational length and birth weight, may be adjusted for without accounting for potential collider stratification bias [1,9,16,17,18].

In this study, we revisit the relationship between early life ambient air pollution exposure and infant mortality in the US with regard to particulate matter less than 2.5 μm in diameter (PM2.5). Specifically, we aim to estimate direct pathways from prenatal and post-birth PM2.5 to infant death as well as indirect pathways from the prenatal exposure to the outcome through preterm birth and low birth weight. Investigating the pathways in such a manner not only allows us to compare the relative associations of prenatal versus post-birth exposure but also helps us understand the relative importance of the indirect versus direct pathway in terms of the prenatal PM2.5–infant death relationship. PM2.5 is an important pollutant to study in this context because of its known association with various indicators of fetal health such as fetal growth and organ development as well as its relationship with infant morbidity, especially respiratory diseases [19,20,21].

2. Materials and Methods

2.1. Data Sources

To define our exposure variables, we used modeled estimates of daily, population-weighted mean PM2.5 concentration at the county level within the conterminous US between 2009 and 2014. These publicly available estimates are provided by the Centers for Disease Control and Prevention and are based on the Environmental Protection Agency’s (EPA) Downscaler model [22,23]. The Downscalar model fuses together modeled estimates of air pollution concentration from the EPA’s Community Multi-Scale Air Quality model with data drawn directly from air pollution monitors [23,24].

To define our outcome and mediators, we used restricted access, cohort-linked birth–infant death data between 2011 and 2013. These data were provided by the National Center for Health Statistics and represent the universe of live births and infant deaths for children born between 2011 and 2013. These data also contain information on the pregnant person’s county of residence at delivery as well as characteristics related to parental demographics, pregnancy, delivery, and infant death.

We collected information on county-level confounders from several publicly available data sources. We extracted information on monthly average temperature and rainfall from the National Oceanic and Atmospheric Administration [25]. From the US Census Bureau, we extracted information on annual county-level demographic and economic variables including racial composition, total population, poverty rate, and median household income [26,27]. Similarly, from the Bureau of Labor Statistics, we extracted information on monthly unstandardized unemployment rate at the county level [28]. Finally, from the Centers for Medicare & Medicaid Services, we extracted information on the annual number of physicians in a county [29].

2.2. Outcome Definition

Our primary outcome was infant death, which we defined as an indicator variable that takes the value one if a baby dies in the first year of life for any reason and zero otherwise.

2.3. Exposure Definition

We defined prenatal air pollution exposure as the average PM2.5 concentration in each trimester in the pregnant person’s county of residence at delivery. Disaggregating the prenatal exposure by trimester reflected results from the literature which suggest that the prenatal air pollution–adverse birth outcome relationship varies by pregnancy trimester [30]. The first trimester was defined as the first three months following the date of conception. The second and third trimesters were defined as the three months following the first and second trimesters, respectively. As such, the prenatal exposure period spanned a total of nine months from the date of conception. We defined the prenatal exposure over a nine-month period as opposed to the actual length of gestation because the latter is a function of prenatal air pollution [9].

We similarly defined post-birth exposure as the average PM2.5 concentration in the pregnant person’s county of residence in the 12-month period following the end of the nine-month prenatal period. We defined post-birth exposure over a 12-month period since the number of days alive in the first year of life may be a function of air pollution exposure as well. We also defined the post-birth exposure from the end of the nine-month prenatal period as opposed to the date of birth to ensure that the prenatal and post-birth periods did not overlap temporally.

2.4. Mediator Definition

We used preterm birth and low birth weight as the two mediators in our analysis. Preterm birth was coded as an indicator variable that equaled one if the obstetric/clinical estimate of gestational length was less than 37 weeks and zero otherwise. Low birth weight was coded as an indicator variable that equaled one if birth weight was reported to be less than 2500 g and zero otherwise.

2.5. Constructing the Analytic Dataset

To construct our analytic sample, we applied the following inclusion criteria to the cohort-linked birth–infant death records: (1) birth records were collected using the latest revision of the birth certificate (i.e., the 2003 birth certificate revision, which had revised a number of instruments for collecting data, including maternal and paternal race as well as cigarette smoking during pregnancy); and (2) pregnant person’s county of residence had daily PM2.5 concentration information available.

To merge the birth–infant death records with the exposure data, we began by estimating the date of conception using information on the year, month, day of week of birth, and the obstetric/clinical estimate of gestational length. Specifically, using the day of the week, month, and year of birth information, we randomly assigned each birth to a date of birth by assuming a uniform probability distribution over the day of week within any given month–year. Then, we subtracted the length of gestation from the estimated date of birth to get our estimate of the date of last menstrual period (LMP). Finally, we assigned the date of conception for each live birth observation by adding two weeks to the LMP date under the assumption that conception occurs, on average, two weeks after the LMP.

Then, we merged the birth–infant death records with the PM2.5 data based on each pregnant person’s reported county of residence and the date of conception. Similarly, we merged the monthly temperature, rainfall, and unemployment data with the live birth–infant death records based on the county of residence and month and year of conception. Finally, we merged data on annual population, poverty rate, housing value, and healthcare access based on the county of residence and year of conception.

2.6. Structural Equation Model

We applied a Structural Equation Modeling (SEM) framework to estimate the direct and indirect pathways from prenatal air pollution exposure to infant death as well as the direct path from post-birth air pollution exposure to the outcome. We also modeled a direct path from prematurity to low birth weight, since birth weight is a function of gestational length [31]. An SEM framework is appropriate for this analysis because it allows us to estimate all pathways simultaneously, which improves statistical power and allows us to estimate total direct and indirect associations easily [32]. A graphical representation of our SEM is presented in Figure 1.

Figure 1

Graphical representation of the Structural Equation Model.

To identify the direct path from prenatal PM2.5 exposure to adverse birth outcomes, we controlled for a variety of individual, pregnancy, and county-level covariates. At the individual level, we controlled for the pregnant person’s age, race, highest level of education, marital status, parity, pre-pregnancy smoking behavior, and average PM2.5 exposure in the nine months prior to conception. We also controlled for father’s age, race, and highest level of education. In terms of pregnancy characteristics, we controlled for the method of payment for delivery, multiple birth, and child’s sex. Finally, at the county level, we controlled for average temperature, precipitation, and unemployment rate over the prenatal period. We also controlled for annual county racial composition, poverty rate, median housing value, and number of physicians per 1000 individuals. To identify the direct path from post-birth PM2.5 exposure to infant mortality, we controlled for these same variables.

To identify the direct path from adverse birth outcomes to infant mortality as well as the path from preterm birth to low birth weight, we additionally controlled for cigarette smoking during pregnancy. We did not control for variables such as number of prenatal care visits because they may plausibly fall on the direct pathway from prenatal air pollution exposure to infant death. As such, controlling for these variables would lead to potentially biased estimates of the direct association between prenatal PM2.5 exposure and infant death.

Finally, to account for unobserved time-invariant county-specific sources of confounding, we allowed for county fixed effects in all models by applying the fixed effects transformation (i.e., demeaning the data at the county level) [33]. In addition, we flexibly modeled time trends in the outcome in all models by including month-of-year fixed effects and a linear, county-specific time trend over the study period. Table A1 presents detailed definitions of all covariates used in the SEM.

Table A1

Covariate definition and data source.

Covariate Definition and Categories	Source of Data
Individual-level variables
Mother’s age (19 years or under; 20–24 years; 25–29 years; 30–34 years; 35–39 years; 40–44 years; 45 years and older)	NCHS linked birth-infant death records
Mother’s race (non-Hispanic white; non-Hispanic black; non-Hispanic other; Hispanic)
Mother’s marital status (Married; Unmarried)
Mother’s education (Without a high school degree, With a high school degree only, Any tertiary degree)
Mother smoked before pregnancy (Yes; No)
Mother smoked during pregnancy (Yes; No)
Parity (First birth, Second birth, Three of more births)
Multiple birth (Singleton; Multiple)
Child’s sex (Female; Male)
Delivery payment source (Medicaid; Private insurance; Self-pay; Other)
Father’s age (19 years or under; 20–24 years; 25–29 years; 30–34 years; 35–39 years; 40–44 years; 45 years and older)
Father’s education (Without a high school degree; With a high school degree only; Any tertiary degree)
Father’s race (non-Hispanic white; non-Hispanic black; non-Hispanic other; Hispanic)Conception month and year
County-of-residence-level variables
Average temperature during pregnancy (°F)	NOAA
Average precipitation during pregnancy (in)
Annual racial composition of county (Proportion non-Hispanic white; proportion non-Hispanic black; proportion non-Hispanic other; proportion Hispanic)	US Census Bureau
Annual population
Annual poverty rate
Annual median household income
Monthly unemployment rate	Bureau of Labor Statistics
Number of physicians	CMS

Notes: NCHS = National Center for Health Statistics; NOAA = National Oceanic and Atmospheric Administration; °F = degrees Fahrenheit; in = inches; CMS = Centers for Medicare & Medicaid Services.

2.7. Statistical Analysis

Before fitting the SEM, we estimated summary statistics of the exposures, mediators, covariates, and outcome. For ease of illustrating summary statistics, we used average PM2.5 concentration over the entire nine-month prenatal period as opposed to each trimester.

To prepare our data for the SEM, we imputed missing values in the analytic sample five times by assuming that the data were missing at random and that the observed and unobserved data followed a multivariate normal distribution. We estimated the SEM in each imputed dataset using the standard maximum likelihood method with an identity link function in all equations. Using the identity link function allowed us to easily estimate direct, indirect, and total associations. We corrected our standard errors by using Huber’s cluster-robust variance estimator at the county-level, which allowed us to account for correlated outcomes within a county [34]. We combined the estimates from the SEM across all imputed datasets by using Rubin’s rules [35].

The literature on air pollution and infant mortality suggests that the exposure–outcome relationship may be non-linear [36]. To account for such potential non-linearities, we categorized both the prenatal and post-birth exposure as <8 μg/m3; 8–10 μg/m3; 10–12 μg/m3; and, ≥12 μg/m3. These cutoffs reflected the distribution of the exposure variables in our dataset (Appendix A Figure A1) and captured the annual air quality guidance of the World Health Organization (WHO; threshold = 10 μg/m3) and the EPA (threshold = 12 μg/m3) [37,38].

Figure A1

Distribution of average prenatal and post-birth PM2.5 concentration; Notes: Panel (a) shows the distribution of average prenatal PM2.5 exposure (defined over 9 months from the date of conception) while panel (b) illustrates the distribution of PM2.5 exposure in the post-birth period (defined over 12 months from the end of the 9-month prenatal period). The red dashed lines on each histogram reflect the threshold values used to create the categorical exposure variable used in the Structural Equation Model. The first red line is at 8 μg/m3, the second is at 10 μg/m3, and the third is at 12 μg/m3.

2.8. Robustness Checks

We conducted two robustness checks. For both robustness checks, we redefined the post-birth exposure variable and re-fit the SEM. Specifically, we redefined the post-birth exposure as average PM2.5 concentration in the one month and two months following the end of the nine-month prenatal period. This reflected the fact that over three-quarters of all infant deaths in the US occurred within the first two months of life (Table A2).

Table A2

Frequency table of number of months alive conditional on dying in the first year of life over the study periods.

Months Alive	Frequency	Proportion
1	40,304	68.4%
2	4833	8.2%
3	3509	6.0%
4	2586	4.4%
5	1966	3.3%
6	1525	2.6%
7	1158	2.0%
8	850	1.4%
9	691	1.2%
10	555	0.9%
11	483	0.8%
12	453	0.8%
Total	58,913

2.9. Software

We used Stata/IC 15 to clean the data, create the analytic sample, and estimate the SEM model [39]. We used RStudio to estimate values of the exposure in the preconception, prenatal, and post-birth periods, create all the figures, and impute the data (package: Amelia II) [40,41].

2.10. Ethical Statement

This study was exempted from human subjects review by the Institutional Review Board at the Harvard T. H. Chan School of Public Health as per regulations found at 45 CFR 46.104(d) (4).

3. Results

Our analytic sample consisted of 10,017,357 live births and 58,913 infant deaths in 3053 counties across 48 states in the conterminous US. In addition to Alaska and Hawaii, Washington DC was also excluded from our analysis due to a lack of air pollution data. Appendix A Figure A2 shows that the primary reason for excluding live birth and infant death observations from our analytic sample was births not being recorded using the 2003 revision of the birth certificate. Missing observations were relatively rare (<10%) except in the case of the father’s age, race, and education, where missingness was between 12% and 15% (Table A3).

Figure A2

Constructing the analytic sample.

Table A3

Frequency and proportion of missing data among variables in the analytic sample.

	Number of Missing Observations	Total Number of Observations	Proportion Missing
Father’s education	1,533,908	10,017,357	15.31
Father’s race	1,404,016	10,017,357	14.02
Father’s age	1,287,549	10,017,357	12.85
Mother smoked cigarettes during pregnancy	608,210	10,017,357	6.07
Mother smoked cigarettes pre-pregnancy	606,380	10,017,357	6.05
Payment source for delivery	144,634	10,017,357	1.44
Post-birth PM2.5 exposure	130,515	10,017,357	1.3
Mother’s education	120,953	10,017,357	1.21
Parity	113,745	10,017,357	1.14
Mother’s race	71,191	10,017,357	0.71
Birth weight	30	10,017,357	~0
Prenatal PM2.5 exposure	0	10,017,357	0
Preconception PM2.5 concentration	0	10,017,357	0
Gestational length	0	10,017,357	0
Infant death	0	10,017,357	0
Mother is married	0	10,017,357	0
Mother’s age	0	10,017,357	0
Singleton delivery	0	10,017,357	0
Child born female	0	10,017,357	0
Average temperature during pregnancy	0	10,017,357	0
Average precipitation during pregnancy	0	10,017,357	0
Average unemployment during pregnancy	0	10,017,357	0
Proportion of Non-Hispanic whites in county	0	10,017,357	0
Proportion of Non-Hispanic blacks in county	0	10,017,357	0
Proportion of Hispanics	0	10,017,357	0
Proportion of Non-Hispanic other races	0	10,017,357	0
Average poverty rate	0	10,017,357	0
Median household income	0	10,017,357	0
Physicians per 1000 individuals	0	10,017,357	0

Notes: The variables are ordered in descending order of missingness. They represent the variables that were used in the Structural Equation Model, either as the exposure, outcome, mediators, or confounders.

Figure 2 shows the annual average PM2.5 concentration in counties where a conception occurred between 2010 and 2013, which were the earliest and latest conception years in our analytic sample. The average annual PM2.5 concentration decreased over the study period from 9.20 μg/m3 in 2010 to 8.34 μg/m3 in 2013. PM2.5 concentration also varied substantially by region: counties in the Interior Midwest, the South, and the Southwest experienced the highest levels of air pollution, while counties in the Great Plains experienced the lowest levels of air pollution.

Figure 2

Average annual PM2.5 concentration by conception year in counties in which at least one conception occurred. Note: Each panel represents average PM2.5 concentration in counties in the conterminous US in the indicated year. Gray-colored counties are counties where no birth occurred in the respective year.

In Figure 3, panels (a) and (b) respectively show unadjusted averages of infant mortality by the four categories of prenatal and post-birth PM2.5 exposure. Prenatal exposure here is defined as the average PM2.5 concentration over the entire nine-month prenatal period as opposed to the average concentration over each trimester. Infant mortality was increasing with prenatal air pollution over the study period; however, in terms of the post-birth exposure, infant mortality increased over the first three air pollution categories and decreased in the final category. Panels (c,d), which present unadjusted proportions of preterm birth and low birth weight by categories of the prenatal exposure variable, follow a similar pattern with the proportion of these two adverse birth outcomes increasing over the first three air pollution categories and then decreasing.

Figure 3

Unadjusted averages of infant mortality, preterm birth, and low birth weight by categories of the prenatal and post-birth PM2.5 exposure. Notes: Prenatal exposure is defined as average PM2.5 concentration in the pregnant person’s county of residence over the entire nine-month prenatal period. We define prenatal exposure over the entire nine-month period for ease of presentation. The categories of both prenatal and post-birth exposure are <8 μg/m3; 8–10 μg/m3; 10–12 μg/m3; and, ≥12 μg/m3. Panel (a) shows unadjusted infant deaths per 1000 live births by categories of the prenatal exposure variable. Panel (b) shows unadjusted infant deaths per 1000 live births by categories of the post-birth exposure variable. Panels (c,d) show unadjusted averages of preterm birth and low birth weight by categories of the prenatal exposure variable.

Approximately three-quarters of the pregnant people in our sample were between 20 and 34 years, a majority were non-Hispanic white, and more than 80% had a high school degree or higher (Table 1). In addition, approximately 48% and 43% of mothers paid for their deliveries using private insurance and Medicaid, respectively. When disaggregating these covariates by categories of the prenatal and post-birth exposure, we find that relative to the lowest air pollution category, the highest category had a lower proportion of individuals with a high school degree or more, higher proportion of deliveries paid for using Medicaid, and a higher average poverty rate (Table A4 and Table A5). Additionally, in terms of the prenatal exposure, the highest air pollution category had a higher proportion of Non-Hispanic Blacks and Hispanics relative to the lowest category (Table A4).

Table 1

Distribution of individual, delivery, and county-level covariates in the analytic sample.

	Mean (SD)
	Overall
Individual-level variables
Preconception PM2.5 concentration (μg/m3)	9.67 (1.76)
Mother’s age
≤19 years	7.83 (26.86)
20–24 years	23.23 (42.23)
25–29 years	28.65 (45.21)
30–34 years	25.59 (43.64)
35–39 years	11.81 (32.27)
40–44 years	2.7 (16.21)
≥45 years	0.19 (4.36)
Mother’s race
Non-Hispanic white	55.33 (49.72)
Non-Hispanic black	14.42 (35.13)
Non-Hispanic other	6.75 (25.09)
Hispanic	23.5 (42.4)
Mother’s education
No high school	16.8 (37.38)
High school/some college	46.47 (49.87)
College or more	36.74 (48.21)
Mother is married	59.54 (49.08)
Mother smoked cigarettes pre-pregnancy	11.71 (32.16)
Parity
First child	32.97 (47.01)
Second child	28.44 (45.11)
Third or more child	38.59 (48.68)
Payment source for delivery
Medicaid	43.16 (49.53)
Private insurance	47.59 (49.94)
Self-pay	4.25 (20.16)
Other	5 (21.8)
Child born female	48.82 (49.99)
Singleton delivery	96.56 (18.22)
Father’s age
≤19 years	3.05 (17.2)
20–24 years	15.3 (36)
25–29 years	25.55 (43.61)
30–34 years	28.41 (45.1)
35–39 years	17 (37.56)
40–44 years	7.29 (26)
≥45 years	3.41 (18.14)
Father’s race
Non-Hispanic white	56.08 (49.63)
Non-Hispanic black	12.75 (33.36)
Non-Hispanic other	7.44 (26.24)
Hispanic	23.73 (42.54)
Father’s education
No high school	15.96 (36.62)
High school/some college	48.62 (49.98)
College or more	35.43 (47.83)
County-level variables
Average temperature during pregnancy	58.05 (9.54)
Average precipitation during pregnancy	3.06 (1.54)
Average unemployment during pregnancy	8.51 (2.48)
County racial composition
Non-Hispanic white	61.99 (22.11)
Non-Hispanic black	12.3 (12.62)
Non-Hispanic other	5.6 (6.14)
Hispanic	18.26 (18.22)
Average poverty rate	16.1 (5.47)
Median household income (USD)	52,311.3 (13,252.36)
Physicians per 1000 individuals	0.35 (0.85)

Notes: SD = standard deviation. USD = United States dollars.

Table A4

Distribution of individual, delivery, and county-level covariates in the analytic sample by categories of prenatal PM2.5 exposure.

	Mean (SD)
	Prenatal PM2.5 concentration
	Category 1	Category 2	Category 3	Category 4
	[3.63 μg/m3, 8.00 μg/m3)	[8.00 μg/m3, 10.00 μg/m3)	[10.00 μg/m3, 12.00 μg/m3)	[12.00 μg/m3, 19.16 μg/m3)
Individual-level variables
Preconception PM2.5 concentration (μg/m3)	7.46 (1.19)	9.43 (1.15)	10.84 (1.15)	11.75 (1.39)
Mother’s age
≤19 years	7.54 (26.4)	7.72 (26.69)	8.01 (27.14)	8.47 (27.84)
20–24 years	23.95 (42.68)	23.36 (42.31)	22.79 (41.95)	22.47 (41.74)
25–29 years	29.29 (45.51)	28.78 (45.28)	28.33 (45.06)	27.61 (44.7)
30–34 years	24.82 (43.2)	25.62 (43.65)	25.95 (43.83)	25.89 (43.8)
35–39 years	11.54 (31.95)	11.7 (32.14)	11.98 (32.48)	12.42 (32.98)
40–44 years	2.68 (16.14)	2.64 (16.04)	2.74 (16.32)	2.93 (16.87)
≥45 years	0.18 (4.27)	0.19 (4.31)	0.2 (4.42)	0.21 (4.6)
Mother’s race
Non-Hispanic white	57.33 (49.46)	58.74 (49.23)	52.16 (49.95)	44.66 (49.71)
Non-Hispanic black	9.8 (29.73)	13.95 (34.65)	16.86 (37.44)	18.43 (38.77)
Non-Hispanic other	6.4 (24.48)	6.67 (24.94)	6.93 (25.4)	7.42 (26.21)
Hispanic	26.47 (44.12)	20.64 (40.47)	24.04 (42.73)	29.49 (45.6)
Mother’s education
No high school	15.94 (36.6)	16.07 (36.73)	17.68 (38.15)	19.31 (39.47)
High school/some college	49.08 (49.99)	46.22 (49.86)	45.53 (49.8)	44.96 (49.75)
College or more	34.98 (47.69)	37.71 (48.47)	36.79 (48.22)	35.73 (47.92)
Mother is married	59.98 (48.99)	60.95 (48.79)	58.37 (49.3)	55.66 (49.68)
Mother smoked cigarettes pre-pregnancy	11.22 (31.56)	12.44 (33)	11.52 (31.92)	9.62 (29.49)
Parity
First child	32.25 (46.75)	33.01 (47.02)	33.17 (47.08)	34.01 (47.37)
Second child	28.28 (45.03)	28.64 (45.21)	28.38 (45.08)	27.89 (44.85)
Third or more child	39.47 (48.88)	38.35 (48.62)	38.45 (48.65)	38.1 (48.56)
Payment source for delivery
Medicaid	44.39 (49.68)	41.88 (49.34)	43.13 (49.53)	47.92 (49.96)
Private insurance	44.7 (49.72)	49.02 (49.99)	48.12 (49.96)	44.39 (49.68)
Self-pay	5.17 (22.14)	4.01 (19.63)	4.16 (19.97)	3.28 (17.81)
Other	5.74 (23.25)	5.09 (21.98)	4.59 (20.93)	4.41 (20.54)
Child born female	48.74 (49.98)	48.85 (49.99)	48.83 (49.99)	48.88 (49.99)
Singleton delivery	96.78 (17.66)	96.5 (18.38)	96.5 (18.38)	96.64 (18.02)
Father’s age
≤19 years	2.96 (16.95)	2.95 (16.92)	3.15 (17.48)	3.43 (18.19)
20–24 years	15.82 (36.49)	15.36 (36.06)	14.98 (35.69)	14.96 (35.67)
25–29 years	26.45 (44.11)	25.69 (43.69)	25.04 (43.32)	24.51 (43.02)
30–34 years	27.55 (44.68)	28.56 (45.17)	28.74 (45.26)	28.31 (45.05)
35–39 years	16.41 (37.04)	16.98 (37.55)	17.27 (37.8)	17.54 (38.03)
40–44 years	7.23 (25.9)	7.16 (25.79)	7.4 (26.18)	7.65 (26.59)
≥45 years	3.57 (18.56)	3.3 (17.86)	3.41 (18.14)	3.6 (18.62)
Father’s race
Non-Hispanic white	57.48 (49.44)	59.3 (49.13)	53.23 (49.9)	46.43 (49.87)
Non-Hispanic black	9.45 (29.26)	12.64 (33.23)	14.42 (35.13)	14.83 (35.54)
Non-Hispanic other	6.87 (25.29)	7.39 (26.17)	7.72 (26.69)	8.04 (27.2)
Hispanic	26.2 (43.97)	20.67 (40.49)	24.63 (43.08)	30.7 (46.13)
Father’s education
No high school	15.45 (36.14)	15.2 (35.91)	16.7 (37.3)	18.37 (38.73)
High school/some college	51.55 (49.98)	48.73 (49.98)	47.25 (49.92)	46.06 (49.84)
College or more	33.01 (47.02)	36.07 (48.02)	36.06 (48.02)	35.57 (47.87)
County-level variables
Average temperature during pregnancy	58.25 (12.4)	58 (8.86)	58.37 (8.79)	55.82 (6.52)
Average precipitation during pregnancy	3.12 (1.97)	3.24 (1.47)	2.89 (1.32)	2.67 (1.46)
Average unemployment during pregnancy	8.41 (2.58)	8.01 (2.4)	8.8 (2.3)	10.37 (2.53)
County racial composition
Non-Hispanic white	64.46 (22.13)	64.9 (20.92)	59 (22.59)	52.06 (21.82)
Non-Hispanic black	7.33 (8.52)	12.07 (13.08)	14.69 (12.96)	16.21 (13.56)
Non-Hispanic other	5.07 (7.49)	5.46 (6.4)	5.89 (5.08)	6.51 (4.82)
Hispanic	21.21 (19.51)	15.66 (17.43)	18.65 (17.45)	23.52 (20.36)
Average poverty rate	15.97 (4.98)	15.76 (6.03)	16.26 (5.07)	17.86 (4.91)
Median household income	50,023.7 (11434.19)	53,576.41 (15402.29)	52,584.86 (11941.65)	49,845.36 (8493.12)
Physicians per 1000 individuals	0.36 (1.15)	0.35 (0.93)	0.33 (0.54)	0.39 (0.44)

Notes: SD = standard deviation. USD = United States dollars.

Table A5

Distribution of individual, delivery, and county-level covariates in the analytic sample by categories of post-birth PM2.5 exposure.

	Mean (SD)
	Post-birth PM2.5 concentration
	Category 1	Category 2	Category 3	Category 4
	[2.70 μg/m3, 8.00 μg/m3)	[8.00 μg/m3, 10.00 μg/m3)	[10.00 μg/m3, 12.00 μg/m3)	[12.00 μg/m3, 17.19 μg/m3)
Individual-level variables
Preconception PM2.5 concentration (μg/m3)	7.4 (1.14)	9.38 (1.1)	10.95 (1.16)	11.82 (1.3)
Mother’s age
≤19 years	7.78 (26.78)	7.89 (26.95)	7.83 (26.86)	7.8 (26.82)
20–24 years	24.39 (42.94)	23.55 (42.43)	22.49 (41.75)	21.43 (41.03)
25–29 years	29.38 (45.55)	28.88 (45.32)	28.33 (45.06)	26.49 (44.13)
30–34 years	24.46 (42.98)	25.39 (43.52)	26.3 (44.03)	26.48 (44.12)
35–39 years	11.23 (31.57)	11.52 (31.93)	12.11 (32.62)	14.02 (34.72)
40–44 years	2.59 (15.89)	2.6 (15.91)	2.75 (16.34)	3.5 (18.37)
≥45 years	0.18 (4.18)	0.18 (4.2)	0.2 (4.48)	0.27 (5.21)
Mother’s race
Non-Hispanic white	57.47 (49.44)	59.92 (49.01)	52.31 (49.95)	30.43 (46.01)
Non-Hispanic black	10.32 (30.43)	14.09 (34.79)	17.29 (37.82)	10.91 (31.18)
Non-Hispanic other	5.75 (23.28)	6.65 (24.92)	6.87 (25.3)	10.8 (31.04)
Hispanic	26.46 (44.11)	19.33 (39.49)	23.52 (42.41)	47.86 (49.95)
Mother’s education
No high school	16.06 (36.71)	16.2 (36.84)	17.26 (37.79)	21.85 (41.32)
High school/some college	49.39 (50)	46.15 (49.85)	45.28 (49.78)	45.52 (49.8)
College or more	34.55 (47.55)	37.65 (48.45)	37.46 (48.4)	32.63 (46.89)
Mother is married	59.21 (49.14)	61.2 (48.73)	58.52 (49.27)	55.08 (49.74)
Mother smoked cigarettes pre-pregnancy	11.47 (31.86)	12.98 (33.61)	11.42 (31.8)	4.14 (19.91)
Mother smoked cigarettes during pregnancy	9.06 (28.7)	9.99 (29.98)	8.64 (28.09)	3.07 (17.24)
Parity
First child	32.29 (46.76)	33.1 (47.06)	33.26 (47.11)	33.46 (47.19)
Second child	28.23 (45.01)	28.64 (45.21)	28.38 (45.08)	28.12 (44.96)
Third or more child	39.48 (48.88)	38.26 (48.6)	38.37 (48.63)	38.41 (48.64)
Payment source for delivery
Medicaid	44.82 (49.73)	42.06 (49.37)	42.27 (49.4)	52.07 (49.96)
Private insurance	43.72 (49.6)	49.04 (49.99)	48.91 (49.99)	41.23 (49.23)
Self-pay	5.5 (22.8)	3.84 (19.21)	4.14 (19.91)	3.57 (18.56)
Other	5.96 (23.67)	5.07 (21.93)	4.68 (21.12)	3.13 (17.41)
Child born female	48.79 (49.99)	48.84 (49.99)	48.87 (49.99)	48.76 (49.98)
Singleton delivery	96.88 (17.39)	96.59 (18.14)	96.55 (18.26)	96.86 (17.43)
Father’s age
≤19 years	3.05 (17.2)	2.98 (17)	3.08 (17.29)	3.51 (18.4)
20–24 years	16.12 (36.77)	15.5 (36.19)	14.69 (35.4)	14.85 (35.56)
25–29 years	26.68 (44.23)	25.88 (43.8)	24.88 (43.23)	23.38 (42.32)
30–34 years	27.37 (44.58)	28.48 (45.13)	29.03 (45.39)	27.37 (44.59)
35–39 years	16.15 (36.79)	16.82 (37.41)	17.47 (37.98)	18.28 (38.65)
40–44 years	7.09 (25.66)	7.1 (25.68)	7.44 (26.24)	8.43 (27.78)
≥45 years	3.55 (18.51)	3.23 (17.69)	3.4 (18.13)	4.18 (20)
Father’s race
Non-Hispanic white	57.75 (49.4)	60.47 (48.89)	53.57 (49.87)	31.59 (46.49)
Non-Hispanic black	9.86 (29.81)	12.67 (33.26)	14.75 (35.46)	9.39 (29.16)
Non-Hispanic other	6.27 (24.24)	7.44 (26.24)	7.6 (26.49)	10.97 (31.25)
Hispanic	26.12 (43.93)	19.42 (39.56)	24.08 (42.76)	48.05 (49.96)
Father’s education
No high school	15.36 (36.05)	15.24 (35.94)	16.22 (36.86)	22.69 (41.88)
High school/some college	52.04 (49.96)	48.64 (49.98)	46.98 (49.91)	46.25 (49.86)
College or more	32.61 (46.88)	36.13 (48.04)	36.8 (48.23)	31.06 (46.27)
County-level variables
Average temperature during pregnancy	58.67 (12.81)	57.64 (9.08)	58.09 (8.25)	57.56 (6.51)
Average precipitation during pregnancy	3.1 (1.92)	3.26 (1.44)	2.97 (1.37)	1.81 (1.29)
Average unemployment during pregnancy	8.45 (2.57)	8.06 (2.42)	8.76 (2.2)	11.13 (2.78)
County racial composition
Non-Hispanic white	64.64 (22.69)	65.9 (20.65)	59.16 (21.91)	38.94 (15.62)
Non-Hispanic black	7.57 (8.41)	12.32 (13.56)	15.03 (13.03)	10.21 (8.06)
Non-Hispanic other	4.49 (6.8)	5.37 (6.54)	5.87 (5.06)	10.04 (5.04)
Hispanic	21.46 (20.41)	14.47 (16.69)	18.16 (16.68)	39.04 (16.16)
Average poverty rate	16.08 (5.03)	15.72 (6.05)	16.14 (4.95)	19.23 (4.49)
Median household income	48923.01 (10438.92)	53399.45 (15823.55)	52909.86 (11474.82)	51726.89 (8044.06)
Physicians per 1000 individuals	0.37 (1.18)	0.35 (0.96)	0.33 (0.46)	0.32 (0.25)

Notes: SD = standard deviation. USD = United States dollars.

Table 2 presents the primary results from our analysis. The average Standardized Root Mean Square Residual (SRMR) of the overall model across all five imputed datasets was less than 0.01, which suggests good fit with the data. Panel D presents associations of trimester-wise prenatal PM2.5 exposure and post-birth PM2.5 exposure with infant death. There are two key take-aways from the results presented in this panel: first, the association between prenatal PM2.5 exposure and infant death is positive in all three trimesters but appears to be strongest in the third trimester. For instance, in the third trimester, relative to the reference exposure category of less than 8 μg/m3, being exposed to, on average, 8–10 μg/m3, 10–12 μg/m3, and ≥12 μg/m3 of PM2.5 concentration was associated with a 0.03 percentage point (95 percent confidence interval: 0.02, 0.05), 0.07 percentage point (0.05, 0.09), and 0.1 percentage point (0.06, 0.13) increase in the risk of infant death. These point estimates are generally larger in magnitude than those estimated for the first and second trimesters, although there is substantial overlap in the 95% confidence intervals. The second take-away is that the risk differences associated with post-birth PM2.5 exposure are less precisely estimated: for instance, relative to the reference category of less than 8 μg/m3, risk difference estimates related to being exposed to 10–12 μg/m3 of PM2.5 concentration over the study period were compatible with anything between a 0.03 percentage point decrease and a 0.05 percentage point increase in the risk of infant death.

Table 2

Estimates of the association between prenatal PM2.5 exposure, post-birth PM2.5 exposure, preterm birth, low birth weight, and infant mortality from the Structural Equation Model.

	Percentage Point Change	95% Confidence Interval
Panel A: Direct association between prenatal PM2.5 exposure and preterm birth
First trimester
Ref: <8 μg/m3
[8.00–10.00 μg/m3)	0.15	[0.03, 0.26]
[10.00–12.00 μg/m3)	0.08	[−0.11, 0.27]
[12.00–19.16 μg/m3]	−0.08	[−0.33, 0.18]
Second trimester
Ref: <8 μg/m3
[8.00–10.00 μg/m3)	0.37	[0.25, 0.49]
[10.00–12.00 μg/m3)	0.62	[0.48, 0.77]
[12.00–19.16 μg/m3]	0.74	[0.56, 0.92]
Third trimester
Ref: <8 μg/m3
[8.00–10.00 μg/m3)	0.38	[0.26, 0.49]
[10.00–12.00 μg/m3)	0.71	[0.55, 0.87]
[12.00–19.16 μg/m3]	1.02	[0.82, 1.21]
Panel B: Direct association between prenatal PM2.5 exposure and low birth weight
First trimester
Ref: <8 μg/m3
[8.00–10.00 μg/m3)	0.07	[0.01, 0.12]
[10.00–12.00 μg/m3)	0.08	[0, 0.15]
[12.00–19.16 μg/m3]	0	[−0.09, 0.08]
Second trimester
Ref: <8 μg/m3
[8.00–10.00 μg/m3)	0.11	[0.06, 0.17]
[10.00–12.00 μg/m3)	0.15	[0.08, 0.23]
[12.00–19.16 μg/m3]	0.23	[0.14, 0.31]
Third trimester
Ref: <8 μg/m3
[8.00–10.00 μg/m3)	0.13	[0.07, 0.18]
[10.00–12.00 μg/m3)	0.22	[0.15, 0.29]
[12.00–19.16 μg/m3]	0.35	[0.26, 0.43]
Panel C: Direct association between preterm birth and low birth weight
Preterm birth	49.67	[49.29, 50.05]
Panel D: Direct association of prenatal and post-birth PM2.5 exposure with infant death
Prenatal exposure
First trimester
Ref: <8 μg/m3
[8.00–10.00 μg/m3)	0.03	[0.01, 0.05]
[10.00–12.00 μg/m3)	0.04	[0.01, 0.06]
[12.00–19.16 μg/m3]	0.03	[0, 0.06]
Second trimester
Ref: <8 μg/m3
[8.00–10.00 μg/m3)	0.04	[0.02, 0.05]
[10.00–12.00 μg/m3)	0.04	[0.01, 0.06]
[12.00–19.16 μg/m3]	0.07	[0.04, 0.09]
Third trimester
Ref: <8 μg/m3
[8.00–10.00 μg/m3)	0.03	[0.02, 0.05]
[10.00–12.00 μg/m3)	0.07	[0.05, 0.09]
[12.00–19.16 μg/m3]	0.1	[0.06, 0.13]
Post-birth exposure
Ref: <8 μg/m3
[8.00–10.00 μg/m3)	0.04	[0.01, 0.07]
[10.00–12.00 μg/m3)	0	[−0.03, 0.04]
[12.00–17.19 μg/m3]	−0.01	[−0.07, 0.05]
Panel E: Direct association of preterm birth and low birth weight with infant death
Preterm birth	2	[1.93, 2.07]
Low birth weight	3.64	[3.55, 3.73]
Number of observations	10,017,357
Average SRMR	0

Notes: All coefficients are expressed as percentage point changes in the respective outcomes. 95% confidence intervals were estimated using standard errors that were clustered at the county-level. The post-birth PM2.5 exposure is estimated over a 12-month period following the end of the nine-month prenatal period. SRMR = Standardized Root Mean Square Residual. The average SRMR was calculated as the average of the SRMR of the Structural Equation Model fit in each of the five imputed datasets.

Panels A and B of Table 2 show the association of prenatal PM2.5 exposure by trimester with preterm birth and low birth weight, respectively. The results presented in these panels suggest that increased PM2.5 exposure in the second and third trimesters was especially strongly associated with increased risk of experiencing the two adverse birth outcomes. For example, the second trimester results in Panel A suggest that being exposed to, on average, 8–10 μg/m3, 10–12 μg/m3, and ≥12 μg/m3 of PM2.5 concentration was associated with a 0.37 percentage point (0.25, 0.49), 0.62 percentage point (0.48, 0.77), and 0.74 percentage point (0.56, 0.92) increase in the risk of prematurity relative to the reference exposure category of less than 8 μg/m3. Panel C presents the association between preterm birth and low birth weight, while Panel E presents the association of the two adverse birth outcomes with infant death. Panel C suggests that preterm birth and low birth weight are strongly, positively associated; similarly, Panel E shows that both adverse birth outcomes are strongly associated with increased risk of infant mortality.

Table 3 shows the direct and indirect associations of prenatal PM2.5 exposure with infant death (Panel A) as well as the direct association of post-birth PM2.5 exposure (Panel B) with our primary outcome. The total association between prenatal air pollution and infant death is positive across all three trimesters; however, as shown in Table 2, the magnitude of the exposure–outcome relationship appears to be strongest in the third trimester. For instance, relative to the lowest exposure category, being exposed to ≥12 μg/m3 of PM2.5 concentration in the third trimester was associated with a 0.15 percentage point (0.11, 0.18) increase in the risk of infant death; in the first trimester, the estimated risk difference was compatible with any value between −0.01 percentage points and 0.07 percentage points. The proportion mediated across all three trimesters is less than 50%, which suggests that the majority of the association between prenatal PM2.5 and infant death is being driven through the direct pathway between the exposure and the outcome.

Table 3

Direct and indirect associations of prenatal and post-birth PM2.5 exposure (defined over 12 months) with infant mortality.

	Direct Association		Indirect Association		Total Association		Proportion Mediated (%)
	Percentage Point Change	95% Confidence Interval	Percentage Point Change	95% Confidence Interval	Percentage Point Change	95% Confidence Interval
Panel A: Direct and indirect associations of the prenatal PM2.5 exposure
First trimester
Ref: <8 μg/m3
[8.00–10.00 μg/m3)	0.03	[0.01, 0.05]	0.01	[0, 0.02]	0.04	[0.01, 0.06]	25%
[10.00–12.00 μg/m3)	0.04	[0.01, 0.06]	0.01	[0, 0.01]	0.04	[0.01, 0.07]	25%
[12.00–19.16 μg/m3]	0.03	[0, 0.06]	0	[−0.02, 0.01]	0.03	[−0.01, 0.07]	0%
Second trimester
Ref: <8 μg/m3
[8.00–10.00 μg/m3)	0.04	[0.02, 0.05]	0.02	[−0.01, 0.05]	0.05	[0.03, 0.07]	40%
[10.00–12.00 μg/m3)	0.04	[0.01, 0.06]	0.03	[−0.02, 0.08]	0.07	[0.04, 0.09]	43%
[12.00–19.16 μg/m3]	0.07	[0.04, 0.09]	0.04	[−0.03, 0.1]	0.1	[0.07, 0.13]	40%
Third trimester
Ref: <8 μg/m3
[8.00–10.00 μg/m3)	0.03	[0.02, 0.05]	0.02	[−0.01, 0.05]	0.05	[0.03, 0.07]	40%
[10.00–12.00 μg/m3)	0.07	[0.05, 0.09]	0.04	[−0.03, 0.1]	0.1	[0.08, 0.13]	40%
[12.00–19.16 μg/m3]	0.1	[0.06, 0.13]	0.05	[−0.04, 0.14]	0.15	[0.11, 0.18]	33%
Panel B: Direct and indirect associations of the postnatal PM2.5 exposure
Ref: <8 μg/m3
[8.00–10.00 μg/m3)	0.04	[0.01, 0.07]	-	-	0.04	[0.01, 0.07]	-
[10.00–12.00 μg/m3)	0	[−0.03, 0.04]	-	-	0	[−0.03, 0.04]	-
[12.00–19.16 μg/m3]	−0.01	[−0.07, 0.05]	-	-	−0.01	[−0.07, 0.05]	-

Notes: The indirect association from the prenatal PM2.5 exposure to infant mortality reflects two paths: prenatal exposure → preterm birth → infant mortality; and, prenatal exposure → preterm birth → low birth weight → infant mortality. The post-birth PM2.5 exposure is estimated over a 12-month period following the end of the nine-month prenatal period. Since there is only a direct path from the post-birth exposure to infant mortality, there are no results for the indirect association between the post-birth exposure and the outcome.

Table A6 and Table A7 present results from our robustness checks in which we redefined the post-birth exposure over a one- and two-month period, respectively. In contrast to our primary results, when we defined the post-birth exposure over a one-month period following the end of the nine-month prenatal period, we were able to estimate a positive association more precisely between post-birth air pollution and infant death across all levels of the exposure. However, as in our primary results, the post-birth air pollution–infant death relationship was less precisely estimated when we defined the post-birth exposure over a two-month period. Table A8 and Table A9 agree with our primary results in that they suggest that a majority of the association between the prenatal air pollution exposure and infant death is driven by the direct path from the exposure to the outcome and not the path through the two adverse birth outcome mediators.

Table A6

Estimates of the association between prenatal PM2.5 exposure, post-birth PM2.5 exposure (defined over 1 month), preterm birth, low birth weight, and infant mortality from the Structural Equation Model.

	Percentage Point Change	95% Confidence Interval
Panel A: Direct association between prenatal PM2.5 exposure and preterm birth
First trimester
Ref: <8 μg/m3
[8.00–10.00 μg/m3)	0.15	[0.03, 0.26]
[10.00–12.00 μg/m3)	0.08	[−0.11, 0.27]
[12.00–19.16 μg/m3]	−0.08	[−0.33, 0.18]
Second trimester
Ref: <8 μg/m3
[8.00–10.00 μg/m3)	0.37	[0.25, 0.49]
[10.00–12.00 μg/m3)	0.62	[0.48, 0.77]
[12.00–19.16 μg/m3]	0.74	[0.56, 0.92]
Third trimester
Ref: <8 μg/m3
[8.00–10.00 μg/m3)	0.38	[0.26, 0.49]
[10.00–12.00 μg/m3)	0.71	[0.55, 0.87]
[12.00–19.16 μg/m3]	1.02	[0.82, 1.21]
Panel B: Direct association between prenatal PM2.5 exposure and low birth weight
First trimester
Ref: <8 μg/m3
[8.00–10.00 μg/m3)	0.07	[0.01, 0.12]
[10.00–12.00 μg/m3)	0.08	[0, 0.15]
[12.00–19.16 μg/m3]	0	[−0.09, 0.08]
Second trimester
Ref: <8 μg/m3
[8.00–10.00 μg/m3)	0.11	[0.06, 0.17]
[10.00–12.00 μg/m3)	0.15	[0.08, 0.23]
[12.00–19.16 μg/m3]	0.23	[0.14, 0.31]
Third trimester
Ref: <8 μg/m3
[8.00–10.00 μg/m3)	0.13	[0.07, 0.18]
[10.00–12.00 μg/m3)	0.22	[0.15, 0.29]
[12.00–19.16 μg/m3]	0.35	[0.26, 0.43]
Panel C: Direct association between preterm birth and low birth weight
Preterm birth	49.67	[49.29, 50.05]
Panel D: Direct association of prenatal and post-birth PM2.5 exposure with infant death
Prenatal exposure
First trimester
Ref: <8 μg/m3
[8.00–10.00 μg/m3)	0.03	[0.01, 0.05]
[10.00–12.00 μg/m3)	0.04	[0.02, 0.07]
[12.00–19.16 μg/m3]	0.04	[0.01, 0.07]
Second trimester
Ref: <8 μg/m3
[8.00–10.00 μg/m3)	0.04	[0.02, 0.06]
[10.00–12.00 μg/m3)	0.04	[0.02, 0.07]
[12.00–19.16 μg/m3]	0.08	[0.05, 0.11]
Third trimester
Ref: <8 μg/m3
[8.00–10.00 μg/m3)	0.03	[0.01, 0.05]
[10.00–12.00 μg/m3)	0.07	[0.05, 0.09]
[12.00–19.16 μg/m3]	0.1	[0.07, 0.13]
Post-birth exposure
Ref: <8 μg/m3
[8.00–10.00 μg/m3)	0.03	[0.02, 0.05]
[10.00–12.00 μg/m3)	0.05	[0.03, 0.07]
[12.00–17.19 μg/m3]	0.05	[0.02, 0.08]
Panel E: Direct association of preterm birth and low birth weight with infant death
Preterm birth	2	[1.93, 2.07]
Low birth weight	3.64	[3.55, 3.73]
Number of observations	10,017,357
Average SRMR	0

Notes: All coefficients are expressed as percentage point changes in the respective outcomes. 95% confidence intervals were estimated using standard errors that were clustered at the county level. The post-birth PM2.5 exposure is estimated over a one-month period following the end of the nine-month prenatal period. SRMR = Standardized Root Mean Square Residual. The average SRMR was calculated as the average of the SRMR of the Structural Equation Model fit in each of the five imputed datasets.

Table A7

Estimates of the association between prenatal PM2.5 exposure, post-birth PM2.5 exposure (defined over 2 months), preterm birth, low birth weight, and infant mortality from the Structural Equation Model.

	Percentage Point Change	95% Confidence Interval
Panel A: Direct association between prenatal PM2.5 exposure and preterm birth
First trimester
Ref: <8 μg/m3
[8.00–10.00 μg/m3)	0.15	[0.03, 0.26]
[10.00–12.00 μg/m3)	0.08	[−0.11, 0.27]
[12.00–19.16 μg/m3]	−0.08	[−0.33, 0.18]
Second trimester
Ref: <8 μg/m3
[8.00–10.00 μg/m3)	0.37	[0.25, 0.49]
[10.00–12.00 μg/m3)	0.62	[0.48, 0.77]
[12.00–19.16 μg/m3]	0.74	[0.56, 0.92]
Third trimester
Ref: <8 μg/m3
[8.00–10.00 μg/m3)	0.38	[0.26, 0.49]
[10.00–12.00 μg/m3)	0.71	[0.55, 0.87]
[12.00–19.16 μg/m3]	1.02	[0.82, 1.21]
Panel B: Direct association between prenatal PM2.5 exposure and low birth weight
First trimester
Ref: <8 μg/m3
[8.00–10.00 μg/m3)	0.07	[0.01, 0.12]
[10.00–12.00 μg/m3)	0.08	[0, 0.15]
[12.00–19.16 μg/m3]	0	[−0.09, 0.08]
Second trimester
Ref: <8 μg/m3
[8.00–10.00 μg/m3)	0.11	[0.06, 0.17]
[10.00–12.00 μg/m3)	0.15	[0.08, 0.23]
[12.00–19.16 μg/m3]	0.23	[0.14, 0.31]
Third trimester
Ref: <8 μg/m3
[8.00–10.00 μg/m3)	0.13	[0.07, 0.18]
[10.00–12.00 μg/m3)	0.22	[0.15, 0.29]
[12.00–19.16 μg/m3]	0.35	[0.26, 0.43]
Panel C: Direct association between preterm birth and low birth weight
Preterm birth	49.67	[49.29, 50.05]
Panel D: Direct association of prenatal and post-birth PM2.5 exposure with infant death
Prenatal exposure
First trimester
Ref: <8 μg/m3
[8.00–10.00 μg/m3)	0.01	[−0.01, 0.03]
[10.00–12.00 μg/m3)	0.02	[−0.01, 0.05]
[12.00–19.16 μg/m3]	0.03	[0, 0.07]
Second trimester
Ref: <8 μg/m3
[8.00–10.00 μg/m3)	0.04	[0.02, 0.05]
[10.00–12.00 μg/m3)	0.04	[0.01, 0.06]
[12.00–19.16 μg/m3]	0.07	[0.04, 0.09]
Third trimester
Ref: <8 μg/m3
[8.00–10.00 μg/m3)	0.03	[0.01, 0.05]
[10.00–12.00 μg/m3)	0.07	[0.05, 0.09]
[12.00–19.16 μg/m3]	0.1	[0.07, 0.13]
Post-birth exposure
Ref: <8 μg/m3
[8.00–10.00 μg/m3)	0.03	[0.01, 0.05]
[10.00–12.00 μg/m3)	0.02	[0, 0.05]
[12.00–17.19 μg/m3]	0	[−0.03, 0.03]
Panel E: Direct association of preterm birth and low birth weight with infant death
Preterm birth	2	[1.93, 2.07]
Low birth weight	3.64	[3.55, 3.73]
Number of observations	10,017,357
Average SRMR	0

Notes: All coefficients are expressed as percentage point changes in the respective outcomes. 95% confidence intervals were estimated using standard errors that were clustered at the county level. The post-birth PM2.5 exposure is estimated over a two-month period following the end of the nine-month prenatal period. SRMR = Standardized Root Mean Square Residual. The average SRMR was calculated as the average of the SRMR of the Structural Equation Model fit in each of the five imputed datasets.

Table A8

Direct and indirect associations of prenatal and post-birth PM2.5 exposure (defined over 1 month) with infant mortality.

	Direct Association		Indirect Association		Total Association		Proportion Mediated (%)
	Percentage Point Change	95% Confidence Interval	Percentage Point Change	95% Confidence Interval	Percentage Point Change	95% Confidence Interval
Panel A: Direct and indirect associations of the prenatal PM2.5 exposure
First trimester
Ref: <8 μg/m3
[8.00–10.00 μg/m3)	0.03	[0.01, 0.05]	0.01	[0, 0.02]	0.04	[0.02, 0.07]	25%
[10.00–12.00 μg/m3)	0.04	[0.02, 0.07]	0.01	[0, 0.01]	0.05	[0.02, 0.08]	20%
[12.00–19.16 μg/m3]	0.04	[0.01, 0.07]	0	[−0.02, 0.01]	0.04	[0, 0.07]	0%
Second trimester
Ref: <8 μg/m3
[8.00–10.00 μg/m3)	0.04	[0.02, 0.06]	0.02	[−0.01, 0.05]	0.06	[0.04, 0.08]	33%
[10.00–12.00 μg/m3)	0.04	[0.02, 0.07]	0.03	[−0.02, 0.08]	0.07	[0.05, 0.1]	43%
[12.00–19.16 μg/m3]	0.08	[0.05, 0.11]	0.04	[−0.03, 0.1]	0.11	[0.08, 0.15]	36%
Third trimester
Ref: <8 μg/m3
[8.00–10.00 μg/m3)	0.03	[0.01, 0.05]	0.02	[−0.01, 0.05]	0.05	[0.03, 0.07]	40%
[10.00–12.00 μg/m3)	0.07	[0.05, 0.09]	0.04	[−0.03, 0.1]	0.1	[0.08, 0.13]	40%
[12.00–19.16 μg/m3]	0.1	[0.07, 0.13]	0.05	[−0.04, 0.14]	0.15	[0.11, 0.18]	33%
Panel B: Direct and indirect associations of the postnatal PM2.5 exposure
Ref: <8 μg/m3
[8.00–10.00 μg/m3)	0.03	[0.02, 0.05]	-	-	0.03	[0.02, 0.05]	-
[10.00–12.00 μg/m3)	0.05	[0.03, 0.07]	-	-	0.05	[0.03, 0.07]	-
[12.00–19.16 μg/m3]	0.05	[0.02, 0.08]	-	-	0.05	[0.02, 0.08]	-

Notes: Direct and indirect association estimates were estimated using the Structural Equation Model presented in Table A6. The indirect association from the prenatal PM2.5 exposure to infant mortality reflects two paths: prenatal exposure → preterm birth → infant mortality; and, prenatal exposure → preterm birth → low birth weight → infant mortality. The post-birth PM2.5 exposure is estimated over a one-month period following the end of the nine-month prenatal period. Since there is only a direct path from the post-birth exposure to infant mortality, there are no results for the indirect association between the post-birth exposure and the outcome.

Table A9

Direct and indirect associations of prenatal and post-birth PM2.5 exposure (defined over 2 months) with infant mortality.

	Direct Association		Indirect Association		Total Association		Proportion Mediated (%)
	Percentage Point Change	95% Confidence Interval	Percentage Point Change	95% Confidence Interval	Percentage Point Change	95% Confidence Interval
Panel A: Direct and indirect associations of the prenatal PM2.5 exposure
First trimester
Ref: <8 μg/m3
[8.00–10.00 μg/m3)	0.01	[−0.01, 0.03]	0.01	[0, 0.02]	0.02	[0, 0.04]	50%
[10.00–12.00 μg/m3)	0.02	[−0.01, 0.05]	0.01	[0, 0.01]	0.03	[−0.01, 0.06]	33%
[12.00–19.16 μg/m3]	0.03	[0, 0.07]	0	[−0.02, 0.01]	0.03	[−0.01, 0.07]	0%
Second trimester
Ref: <8 μg/m3
[8.00–10.00 μg/m3)	0.04	[0.02, 0.05]	0.02	[−0.01, 0.05]	0.05	[0.03, 0.07]	40%
[10.00–12.00 μg/m3)	0.04	[0.01, 0.06]	0.03	[−0.02, 0.08]	0.07	[0.04, 0.09]	43%
[12.00–19.16 μg/m3]	0.07	[0.04, 0.09]	0.04	[−0.03, 0.1]	0.1	[0.07, 0.13]	40%
Third trimester
Ref: <8 μg/m3
[8.00–10.00 μg/m3)	0.03	[0.01, 0.05]	0.02	[−0.01, 0.05]	0.05	[0.03, 0.07]	40%
[10.00–12.00 μg/m3)	0.07	[0.05, 0.09]	0.04	[−0.03, 0.1]	0.1	[0.08, 0.13]	40%
[12.00–19.16 μg/m3]	0.1	[0.07, 0.13]	0.05	[−0.04, 0.14]	0.15	[0.11, 0.18]	33%
Panel B: Direct and indirect associations of the postnatal PM2.5 exposure
Ref: <8 μg/m3
[8.00–10.00 μg/m3)	0.03	[0.01, 0.05]	-	-	0.03	[0.01, 0.05]	-
[10.00–12.00 μg/m3)	0.02	[0, 0.05]	-	-	0.02	[0, 0.05]	-
[12.00–19.16 μg/m3]	0	[−0.03, 0.03]	-	-	0	[−0.03, 0.03]	-

Notes: Direct and indirect association estimates were estimated using the Structural Equation Model presented in Table A7. The indirect association from the prenatal PM2.5 exposure to infant mortality reflects two paths: prenatal exposure → preterm birth → infant mortality; and, prenatal exposure → preterm birth → low birth weight → infant mortality. The post-birth PM2.5 exposure is estimated over a two-month period following the end of the nine-month prenatal period. Since there is only a direct path from the post-birth exposure to infant mortality, there are no results for the indirect association between the post-birth exposure and the outcome.

4. Discussion

We studied the associations of prenatal and post-birth exposure to PM2.5 with the risk of infant death for all births that took place in the US between 2011 and 2013. We found that increased prenatal exposure to PM2.5 across all three trimesters was associated with an increased risk of infant death, although associations were strongest for air pollution exposure in the third trimester. We also found that the majority of the association between prenatal PM2.5 and infant death was driven by the direct path from the exposure to the outcome rather than the indirect paths through preterm birth and low birth weight. We were able to estimate the association less precisely between post-birth PM2.5 and infant death in our primary analysis; however, results from our robustness checks suggest the possibility of a positive association between the post-birth air pollution exposure and the infant death outcome.

Overall, the results from our study provide evidence in favor of increased air pollution exposure during gestation and possibly during the first year of life being harmful in terms of infant health, even at levels below the threshold set by the WHO and the EPA. This implication is consistent with recent findings from Di et al. (2018), which suggested that increased PM2.5 exposure below the EPA standard was associated with an increased risk of death among older Americans [42]. It is also nominally consistent with the results from Chay and Greenstone (2003), who demonstrated a positive association between total suspended particulates and infant death at levels below the EPA mandated threshold, although their analysis used data from the early 1980s when the EPA threshold was different from what it is currently [8].

Additionally, our result that over 50% of the association between prenatal PM2.5 exposure and infant death can be explained by the direct pathway suggests that there are mechanisms other than preterm birth and low birth weight by which in utero PM2.5 exposure may affect the risk of infant death. Previous studies have suggested that prenatal air pollution exposure may be linked with intrauterine growth retardation or congenital heart defects, all of which could affect the risk of infant mortality [43,44]. Future research should consider explicitly characterizing these pathways.

Finally, the discrepancy between the results from our primary analysis and robustness checks for the post-birth exposure highlights the complexities in defining the exposure in studies investigating health outcomes in early life. A strength of our analysis relative to the literature is that we define our exposure variables independent of the actual length of gestation or time alive. While this allows us to avoid defining the exposure using metrics that themselves may be a function of the exposure, it also introduces error into the exposure variables, which could affect both the point estimate and its associated standard error. Future research should consider determining the most appropriate methods of defining in utero and post-birth exposures.

Our study improves on the existing literature in several ways: first, we used high quality air pollution and infant death data with wide geographical coverage. Second, following lessons from the causal mediation literature, we carefully controlled for exposure–outcome, exposure–mediator, and mediator–outcome confounders [45]. We included county fixed effects in all our models, which allowed us to control for any time-invariant sources of confounding at the county level. We also modeled the time trend flexibly at the month-of-year level and the county level, which allowed us to net out any trends in the outcome. Third, we used an SEM framework to measure the different pathways from prenatal and post-birth PM2.5 exposure to infant death. To the best of our knowledge, our paper is one of the first to use SEM to understand this system of relationships.

Nevertheless, our analysis is still subject to several limitations. First, we were unable to disaggregate the pollution measure by type of PM2.5 pollutant. Previous studies have shown that there is a high degree of heterogeneity in terms of particulate matter composition in the US by region and season [46]. Furthermore, studies have also suggested that early life exposure to carbonaceous PM2.5 changes the risk of infant death, but sea salt or mineral dust does not [47]. Our use of county fixed effects somewhat addresses this issue by ensuring that we only use within-county variation in the exposure for our analysis.

Second, we did not include other pollutants in our model due to a lack of high-quality, daily data on them. Pollutants such as ozone, carbon monoxide, and nitrogen dioxide have been shown to be associated with infant death, while other studies have indicated that these pollutants may be correlated with PM2.5 levels as well [48,49,50,51,52,53,54,55,56]. Our inability to control for these pollutants means that we cannot rule out that the associations we present reflect the relationship between general pollution levels and infant death.

Third, the fact that air pollution data were only available at the county level may have introduced measurement error into our exposure variable. However, this issue is somewhat mitigated by the fact that the air pollution exposure is weighted based on the county’s population.

Fourth, in defining the exposure variables, we only have the pregnant individual’s county of residence at the time of delivery. Therefore, we are not able to account for changes in county of residence during pregnancy or in the first year of the child’s life. In addition, we are also not able to account for movement across counties for work or visits that might affect an individual’s exposure levels. Although there are relatively few studies examining the movement of pregnant individual’s during pregnancy, the existing literature suggests that when pregnant people do move, the median distance travelled is under 10 km [57,58]. This suggests that any issues due to the movement of individuals during pregnancy may be quite small.

Fifth, we were unable to account for indoor air pollution in our analysis. Despite being implicated in increasing the risk of several morbidities such as respiratory ailments and heart disease, we have relatively low understanding of the level of indoor air pollution in the US [59,60]. Future studies should consider studying the interaction between the level and type of indoor and outdoor air pollutants on infant morbidity and mortality.

Finally, a small but growing literature suggests that increased in utero air pollution exposure may increase the risk of pregnancy loss [61,62]. Since we use vital statistics data, our study effectively conditions on live births to analyze the relationship between prenatal air pollution exposure and infant death. As such, our estimates may be subject to a form of selection bias known as live birth bias [63,64,65]. Future studies should further investigate the degree to which differential fetal loss by levels of prenatal air pollution exposure impacts investigations of the relationship between early life air pollution exposure and infant health outcomes.

5. Conclusions

We showed that prenatal PM2.5 exposure was positively associated with infant mortality, with much of the association being driven by the direct path from the exposure to the outcome. Estimates for the relationship between post-birth PM2.5 exposure and infant mortality were positive but less precisely estimated. Our study motivates the need to understand if specific PM2.5 pollutants are responsible for these observed relationships or whether the associations we present represent general air pollution effects. Importantly, our study also adds to the literature that provides evidence in favor of low levels of air pollution being potentially harmful for human health. As more high-quality evidence on this topic accumulates, it may motivate the need to rethink the existing air pollution standards set by the WHO and the EPA.