==== Front Int J Environ Res Public Health Int J Environ Res Public Health ijerph International Journal of Environmental Research and Public Health 1661-7827 1660-4601 MDPI 33260804 10.3390/ijerph17238805 ijerph-17-08805 Article Estimating the Combined Effects of Natural and Built Environmental Exposures on Birthweight among Urban Residents in Massachusetts https://orcid.org/0000-0003-4453-0401Yitshak-Sade Maayan 1* https://orcid.org/0000-0002-1658-3349Fabian M. Patricia 2 Lane Kevin J. 2 Hart Jaime E. 34 Schwartz Joel D. 34 Laden Francine 34 James Peter 35 https://orcid.org/0000-0002-8659-6794Fong Kelvin C. 36 Kloog Itai 7 Zanobetti Antonella 3 1 Department of Environmental Medicine and Public Health, Icahn School of Medicine at Mount Sinai, New York, NY 10029, USA 2 Department of Environmental Health, Boston University School of Public Health, Boston, MA 02118, USA; pfabian@bu.edu (M.P.F.); klane@bu.edu (K.J.L.) 3 Harvard T.H. Chan School of Public Health, Boston, MA 02115, USA; rejch@channing.harvard.edu (J.E.H.); joel@hsph.harvard.edu (J.D.S.); Francine.laden@channing.harvard.edu (F.L.); pjames@hsph.harvard.edu (P.J.); kcf502@mail.harvard.edu (K.C.F.); azanobet@hsph.harvard.edu (A.Z.) 4 Channing Division of Network Medicine, Department of Medicine, Brigham and Women’s Hospital and Harvard Medical School, Boston, MA 02115, USA 5 Department of Population Medicine, Harvard Medical School and Harvard Pilgrim Health Care Institute, Boston, MA 02215, USA 6 School of the Environment, Yale University, New Haven, MA 06511, USA 7 Department of Geography and Environmental Development, Faculty of Humanities and Social Sciences, Ben-Gurion University, Beer-Sheva 84105, Israel; ikloog@bgu.ac.il * Correspondence: maayan.yitshak-sade@mssm.edu 27 11 2020 12 2020 17 23 880513 10 2020 21 11 2020 © 2020 by the authors.2020Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (http://creativecommons.org/licenses/by/4.0/).Intrauterine growth has health implications both in childhood and adulthood. Birthweight is partially determined by prenatal environmental exposures. We aim to identify important predictors of birthweight out of a set of environmental, built environment exposures, and socioeconomic environment variables during pregnancy (i.e., fine particulate matter (PM2.5), temperature, greenness, walkability, noise, and economic indices). We included all singleton live births of mothers who resided in urban census block-groups and delivered in Massachusetts between 2001 and 2011 (n = 640,659). We used an elastic-net model to select important predictors of birthweight and constructed a multivariate model including the selected predictors, with adjustment for confounders. We additionally used a weighted quantile sum regression to assess the contribution of each exposure to differences in birthweight. All exposures were selected as important predictors of birthweight. In the multivariate model, lower birthweight was significantly associated with lower greenness and with higher temperature, walkability, noise, and segregation of the “high income” group. Treating the exposures individually, nighttime noise had the highest weight in its contribution to lower birthweight. In conclusion, after accounting for individual confounders, maternal environmental exposures, built environment exposures, and socioeconomic environment during pregnancy were important predictors of birthweight, emphasizing the role of these exposures in fetal growth and development. birthweightexposomeair pollutiontemperaturenoisegreenness ==== Body 1. Introduction Low birthweight is a strong predictor of perinatal morbidity and mortality [1]. Intrauterine growth has implications on cardiac health, bone development, and other health consequences later in life [2]. Studies suggest that up to 80% of the variation in birthweight is determined by genetic factors [3]. The residual variance is determined partially by maternal health and exposures during pregnancy [4,5,6,7,8,9,10]. Environmental exposures are estimated to contribute about 25% of the variation in birthweight [3]. The contribution of environmental exposures to birthweight reduction might not be large enough to result in a shift from normal to low birthweight (i.e., birthweight < 2500 g. However, given the multiple exposures, and the extent of the population exposed on a daily basis throughout the pregnancy, birthweight differences within the normal range may have meaningful implications on public health on the population level. There is extensive literature documenting the association between multiple environmental exposures (e.g., air pollution, temperature, and noise) [11,12,13,14,15] or neighborhood built environment characteristics (e.g., greenness) [16] and birthweight. However, most studies assessed the effect of a single environmental exposure, either without considering other exposures or with the consideration of one additional exposure (usually temperature or fine particulate matter (PM2.5)) [17,18]. Given the wide range of environmental exposures that can potentially affect birthweight, it has been challenging to isolate important environmental predictors of birthweight and to identify the exposures that are more highly weighted in their contribution to differences in birthweight. There is a gap in studies that consider the impact of multiple environmental exposures on birthweight while accounting for the correlations between individual and environmental risk factors. For decades, the research focus of the effects of external environmental exposures has been on the effects of single pollutants. In recent years, the exposome concept has emphasized the need for epidemiological studies that consider many exposures [19]. Only a few recent studies have taken such an exposomic approach when investigating the environmental predictors of birthweight. Hystad et al., for example, have assessed the effect of greenness exposure on birthweight while accounting for multiple environmental exposures (e.g., PM2.5, black carbon, NO2, noise, walkability, and park proximity) [20]. Another recent study of six European cohorts have estimated associations between low birthweight and multiple exposures defined as the urban exposome including the built environment, air pollution, road traffic noise, meteorology, natural space, and road traffic [21]. Finally, Agier et al. have examined over 100 predictors of birthweight and found lead exposure as a significant risk factor for lower birthweight. After correcting for multiple testing, no other exposure was significantly associated with birthweight [22]. Still, the ability to account for multiple exposures is often limited by sample size and statistical power or by the availability of exposure data. In addition, with multiple spatially and temporally correlated exposures, a priori selection of predictors can be challenging. In the current study, we apply a combination of variable selection, traditional, and novel regression methods. Using a large cohort of births in Massachusetts, and a wide range of prenatal exposures, we identified important predictors of birthweight out of a set of environmental exposures, built environment characteristics and economic indices, and assessed the relative contribution of each of the selected exposures to differences in birthweight. 2. Materials and Methods 2.1. Study Population Our study population were all singleton live births for the years 2001–2011 (n = 791,324) extracted from the Massachusetts birth registry data. For each birth, the data included the newborn’s birthweight in grams and sociodemographic and clinical characteristics: date of birth, government support for prenatal care (yes/no—prenatal care paid by the government), race (black/white/other), maternal age at birth, parity (first-born, not first-born), maternal smoking prior to or during pregnancy (yes/no), maternal highest level of education (less than high school, high school, some college, college, advanced degree beyond college), diabetes mellitus (yes/no), chronic hypertension (yes/no), gestational age (in weeks), and the newborn sex. We excluded birthweight < 500 g (n = 76) due to higher probabilities of lower birthweight related to high-risk pregnancies and prematurity. We additionally excluded records of gestational age > 42 weeks (n = 6362), birthweight over 6 kg (n = 133), and records that had missing residence information or missing covariate data (n = 116,175) to avoid misclassification due to inaccurate or missing data. Finally, we excluded records of mothers who resided in rural areas (n = 38,501) as most of our exposures of interest had little or no variability in non-urban areas. We used the 2010 Census urban and rural classification to define urban/rural areas. This classification is based on population density, and nearby geographic areas containing non-urban land use [23]. In total, we excluded 19% (n = 150,665) of all births. The excluded population was similar except for a higher proportion of white women, college graduates, and a higher gestational age at birth (Table S1). For each birth, the maternal residence address was geocoded by the Massachusetts Department of Public Health. This study was approved by the Massachusetts Department of Public Health (approval number 249812) and the human subjects committee at the Harvard T. H. Chan School of Public Health (approval number 20222). 2.2. Outcome Our outcome was birthweight in grams. Birthweight is commonly treated as a marker of the intrauterine environment and is an important predictor of normal growth, healthy development, and survival [2,24]. 2.3. Exposures We a priori selected potential environmental exposures based on associations with birthweight in previous studies [7,18,20,25,26]. We examined the following exposures. PM2.5: We obtained daily average PM2.5 predictions from a model generating predictions at a 1 × 1 km spatial resolution. Out-of-sample “tenfold” cross-validation of the prediction model showed excellent model performance (mean out-of-sample R2 = 0.88). For a more in-depth description, refer to Kloog et al. [27]. Temperature: We obtained daily average temperature predictions from a model that generated predictions at a 1 × 1 km spatial resolution. Out-of-sample tenfold cross-validation showed excellent model performance (mean out-of-sample R2 = 0.94). For more in-depth information refer to Kloog et al. [28]. For both the PM and Air temperature modeling we implemented the same rigorous cross validation scheme. We used a “tenfold out of sample” cross-validation scheme to validate our model prediction versus Environmental Protection Agency (EPA) and National Climatic Data Center (NCDC) ground monitoring stations. We randomly divided our data into 90 and 10 percent splits ten times. We predicted values for the 10% data sets using the model fitted from the remaining 90% of the data. To test our results for bias we regressed the measured PM or temperature value for a given site and day against the corresponding predicted value. We estimated the model prediction precision by taking the square root of the mean squared prediction errors (RMSPE). In addition, we calculated prediction errors from a model that contained the spatial components only and calculated temporal R2 by regressing Delta measurements of monitoring stations against Delta predicted values where: Delta is the difference between the actual PM/Temperature in place i at time j and the annual mean PM/Temperature at that location, and Delta predicted is defined similarly for the predicted values generated from the model. Spatial R2 was calculated by regressing the site-specific annual means in observed PM/Temperature versus the same annual means for predicted PM/Temperature. Models were validated globally but also per region (i.e., Massachussetts) [27,28]. Greenness: We used monthly remote sensing data from the MODIS satellites, and averaged the values within each season to calculate the seasonal normalized difference vegetation index (NDVI). NDVI incorporate the plants’ absorption and reflection of light. Its values range from −1.0 to 1.0 with higher values indicate higher levels of vegetative density. Based on previous studies of neighborhood greenness, we used a 250 m buffer around each address to approximate the measure of residential and walking greenness exposure [29,30,31]. Walkability index: We used a walkability index which incorporates three census block group level components obtained from the Environmental Protection Agency (EPA) Smart Location Database. [32,33] The index combines z-scores of population density (people/acre) on unprotected land; street intersection density (weighted, with auto-oriented intersections eliminated); and land use diversity based on the mix of retail, office, service, industrial, entertainment, education, healthcare, and public administration employment in the census block group. A higher walkability index indicates a more walkable neighborhood [32]. Nighttime noise: We defined median nighttime noise as the exposure of interest. We obtained noise exposure from a non-time varying geospatial sound model that provides nighttime noise level at a 270 m × 270 m spatial resolution [34]. The model incorporates geospatial features (i.e., climate, anthropogenic activity, and empirical acoustical data) sampled at 492 sites across the U.S. between 2000 and 2014. The method utilized random forest, a tree-based machine learning algorithm to predict noise based on geospatial features. The noise model was constructed using a training set containing all available observations, and the fit showed excellent correlation with the empirical data (R2 = 0.98). In addition, the model was cross-validated using an exhaustive, leave-one-out process to evaluate the accuracy of national scale projections. Each site was omitted, a new model was constructed, and the model output was tested against the observations for the omitted site. The cross-validation procedure showed a root mean squared error of 4.91 dB. As rare conditions are more prone to error, the main systematic source of error was overestimation of sound levels at the quietest sites and underestimation of sound levels at the loudest sites [34,35]. Socioeconomic environment: we used two measures of economic segregation and dissimilarity to assess the neighborhood’s socioeconomic status. We used the formulas published by Krieger et al. [36] to calculate measures of economic segregation and dissimilarity at the block group level. These indices have previously been used to assess how the geographic distribution of economic groups affects human health [37,38]. a. Economic residential segregation (ERS): measures the uniformity of the proportion of people living in high income households (≥$10,000/year) versus the proportion of people living in low income households (<$25,000/year). The ERS ranges between −1 and 1 and grouped into five categories, where the middle category indicates equal distribution of the two extreme economic groups, positive values closer to 1 indicate higher segregation of the “high income” group, and negative values closer to −1 indicate higher segregation of the “low income” group. b. Index of economic dissimilarity (IED): compares the distribution of the proportion of people living in households with high and low income in a smaller area (e.g., block groups) versus the distribution in a larger area (e.g., census tracts) (Table S2). We linked all the exposures to the residential address at birth. Specifically, we assigned temperature, PM2.5, noise, and greenness exposures to each mother based on proximity to the residential address at birth. We assigned walkability and economic segregation and dissimilarity indices to each mother based on the block group of residential address. For walkability, noise and economic segregation and dissimilarity indices, we assumed fixed exposures over time. PM2.5, temperature, and greenness predicted exposures did vary with time; we therefore calculated the average exposure for each trimester of pregnancy. 2.4. Statistical Analysis To assess the relative contribution of multiple environmental exposures during pregnancy on birthweight, we followed several steps. We used Spearman correlation tests to understand the correlations between all exposures. Then, to describe the direction of the relationship between each exposure and birthweight, we fitted separate linear regression models for each exposure. Before doing so, we assessed the linearity of the relationship between the exposures and birthweight by creating a multivariate model including penalized spline functions for all the exposures. Both the shape dose response curve and the potential biological plausibility were considered in our decision whether to include nonlinear terms. Due to the complexity of the analysis, we tended to include linear terms. Although there is evidence of nonlinear effects between temperature exposure and birth outcomes [4,12,39], we found a linear dose response curve for temperature, and therefore treated the exposure as linear. For ERS, we found a nonlinear relationship, and therefore created a categorical variable as described in the Exposure section. The association with nighttime noise was not linear, showing a protective effect for levels higher than 46 dB. However, as the majority of the noise distribution was lower than 46 dB (80%), and as there is probably no reasonable explanation for the protective effect other than a measurement error, we decided to treat noise exposure as linear. For all other exposures, we used linear terms (Figure S1). We adjusted each model for the covariates available in the birth records shown in prior studies as potential confounders of the associations tested in our study [9,39]: season and year of birth, government support for prenatal care, race, age at birth, parity, maternal smoking before or during pregnancy, maternal highest level of education attained, diabetes mellitus, chronic hypertension, and gestational age at birth [5,9,29]. We defined the season of birth as follows; Winter—December to February, Spring—March to May, Summer—June to August, and Fall—September to November. Multiple linear regression is the classical approach to examine the association between maternal or neonatal risk factors and birthweight. However, in cases of multicollinearity (when the exposures are highly correlated) effect estimates obtained from a standard multivariate regression model may be unreliable. In cases of over fitting (where the number of predictors in the model approaches the number of observations in the analyzed population), the variance can be inflated, providing insufficient statistical power to detect small effects. To reduce the impact of collinearity and model overfitting, we used elastic net regression modeling to select the exposures (PM2.5, temperature and NDVI in each trimester, walkability, noise, ERS, and IED) to be included in the final multivariate model [40,41,42]. The method applies a weighted penalty to the regression coefficients to achieve asymptotic normality and consistent selection [40,41,42]. In this model, to fully adjust for all potential confounders, we included all the fixed covariates mentioned above, and we applied the penalty to the standardized exposures but not the confounders. We selected the model that minimizes the cross-validation error, with the alpha value that minimized the mean standardized error (0.1), which yields a model result in which some coefficients are shrunk exactly to zero and some are non-zero, and we selected the exposures with non-zero coefficients. As the elastic net model leaves the redundant covariates out of the model, we refer to the selected exposures throughout the manuscript as important predictors of birthweight. Finally, we ran a multivariate linear regression model, including all the selected exposures, to obtain the unpenalized direct effect estimates of each exposure on birthweight. We present the results as the difference in birthweight (with 95% confidence intervals (CI) for an IQR rise in exposures). As a secondary analysis, to assess the relative contribution of each of the selected exposures to birthweight at term, we used a Weighted Quantile Sum (WQS) regression. The WQS regression summarizes all the exposures into one index while taking the association with the outcome into account, and the contribution of each exposure is weighted based on its relevance to the overall association with the outcome. The weights were assigned to quartile-scored exposures within the composite index and were constrained to be between zero and one and summing up to one. Therefore, the exposures with the highest weights are those which contribute most to the outcome, taking the correlation between the exposures into account [43]. We dedicated 60% of the data for training and 40% for validation. As WQS can only incorporate linear predictors, we treated ERS as a confounder, rather than an exposure, for this analysis. As a WQS model requires all the associations to follow a consistent direction, we could only include exposures that are negatively associated with birth weight. For the WQS analysis we therefore transformed the greenness exposures, which were positively associated with birthweight, into new variables where higher values indicate lower greenness. 2.5. Sensitivity Analyses As exposures are clustered within block-groups, we repeated the multivariate regression model adding a random intercept for each block-group of residence. Second, to make sure the associations found in our study are not confounded by socioeconomic status (SES), we repeated the models with adjustment for census block-group median household income and percent poverty. In addition, we have restricted our study population to term-births, to make sure our findings are robust. Finally, we added two WQS sensitivity analyses calculating the weighted contribution of each of the selected exposures to a negative difference in birthweight with weights assigned to tertile-scored and quintile-scored exposures. 3. Results Table 1 shows the maternal and neonatal characteristics of the 640,659 births. The mean birthweight was 3379.3 g, 48.8% of the newborns were female, and the mean gestational age at birth was 39.0 weeks. The mothers were 30.1 years old on average, 69.9% were white, and 33.1% received governmental support. Figures S2–S8 show the maps with the distributions of the exposures across urban census block-groups. The mean (standard deviation) of the exposures during pregnancy were 10.4 (1.9) µg/m3 for PM2.5, 11.1 (4.7 °C) for temperature, and 0.5 (0.2) for NDVI. For the walkability index, where higher values indicate neighborhoods that encourage walking, the values ranged between −2.8 and 16.9 and the mean index was 1.1 (1.9). The mean nighttime noise level was 43.5 (2.9 dB). Thirty-nine percent of the mothers resided in block-groups with similar distribution of the high (≥$100 K/year) and low (<$25 K/year) income groups. The mean index of economic dissimilarity was 16.9%, indicating that in general, the distribution of the “high income” and “low income” groups were similar at the block-group and the census tract level (Table 2). We observed large negative correlations between the average temperature in the first and third trimesters (r = −0.9) and correlations around 0.5 between PM2.5 exposure in the first and third trimester and between temperature and NDVI exposures (Figure 1). Table 3 shows the multivariable adjusted difference in birthweight for an IQR rise in each exposure, obtained from single exposure models. We observed lower birthweight associated with higher PM2.5 exposure in the three trimesters of pregnancy; temperature exposure in the 2nd and 3rd trimester; and with higher walkability, noise, and IED. Conversely, we found that higher greenspace in the three trimesters of pregnancy and residing in block groups with more residents of the ‘high income’ group were associated with higher birthweight. In the multivariate model including all exposures, temperature exposure in all three trimesters was significantly associated with birthweight. Although attenuated due to adjustment to the other covariates, the associations with noise, ERS, walkability, and NDVI in the 1st and 3rd trimesters remained significantly associated with birthweight. IED, NDVI in the 2nd trimester, and PM2.5 were no longer associated with birthweight (Table 3). The associations with the full set of model covariates are presented in Table S3. All the effect estimates for the associations with birthweight were obtained from linear regressions. The single exposure models tested the association with each exposure separately, with adjustment for maternal and neonatal covariates. The multivariate model includes all the exposures selected by the elastic net model, with adjustment for maternal and neonatal covariates. Each model was adjusted for the following covariates: season and year of birth, government support for prenatal care, the maternal highest level of education attained, race, maternal age, parity, maternal smoking before or during pregnancy, diabetes mellitus, chronic hypertension, and gestational age at birth. We present the results as the difference in birthweight (in grams) for an interquartile (IQR increase in the exposures). We conducted several sensitivity analyses to make sure the associations found in our study were robust. The addition of random effects to the multivariate regression yielded larger standard errors, as expected, but the significance and the magnitude of the effects remained similar for all the exposures except for temperature. With the inclusion of a random intercept, temperature in the first and second trimesters was no longer significantly associated with birthweight (Table S4a). The inclusion of neighborhood SES adjustment (Table S4b) and the restriction of the study population to term births (Table S4c) did not change the results. Figure 2 shows results from the WQS showing the weighted contribution associated with rises in each of the exposures and birthweight. Figure 2 shows the weighted contribution of each exposure to a negative difference in birth weight. Results were obtained from a WQS regression. The WQS regression summarizes all the exposures into one index while taking the association with the outcome into account, and the contribution of each exposure is weighted based on its relevance to the overall association with the outcome. The weights were assigned to quartile-scored exposures within the composite index and were constrained to be between zero and one and summing up to one. The combined index of exposures was significantly associated with lower birthweight. We found nighttime noise to be the largest contributor to lower birthweight, accounting for 18% of the weights. This result was consistent across the sensitivity analyses using different quantiles. Temperature in the third trimester was equally weighted to noise, but results differed in the sensitivity analyses. IED, temperature in the first trimester, and PM2.5 exposures across pregnancy had the smallest contribution to lower birthweight. The relative contribution of NDVI exposures largely differed between the analyses. However, when summing all three trimesters of exposure, NDVI accounted for 34% of the wights in the main analysis and 40% of the weights in the sensitivity analyses (Figures S9 and S10). 4. Discussion We found environmental exposures, built environment exposures, and economic indices during pregnancy to be important predictors of birthweight, even after accounting for known individual maternal and neonatal risk factors. We observed lower birthweights in infants born to women with higher temperature exposures during pregnancy, living in areas with less greenness and higher noise, living in more walkable areas, and in areas with more of the “low income” population. Treating the exposures individually, nighttime noise had the highest weight in its contribution to lower birthweight (graphical abstract). Summing trimestral exposures together, NDVI was the most highly weighed exposure. Our findings are similar to previous studies where associations with birthweight were observed with maternal exposure to various environmental exposures, among them PM2.5 [11,17,25,44,45], temperature [18,46,47], noise [6,7], greenness [5,9], and walkability [20]. Although, most of these studies have assessed the effect of a single exposure at a time, some did try to distinguish between these spatially correlated exposures and studied the effects of multiple exposures simultaneously. For example, Smith et al. found an increased risk for low birthweight associated with PM2.5 exposure, independent of road traffic noise. Traffic noise, however, was no longer associated with the outcome once PM2.5 was controlled for [6]. Hystad et al. explored the pathways related to the protective effect of greenness on birthweight and found that the beneficial effect of greenness was robust to adjustment for PM2.5, black carbon, NO2, noise, walkability, and park proximity [20]. Although we observed significant associations of PM2.5 with birthweight in the single exposure models, the effects were attenuated in the multi-exposure model. Moreover, although PM2.5 exposures were selected by the elastic net, the effect sizes were very small. It is possible that our analysis was underpowered to detect a significant effect of PM2.5 exposure in a multivariate model, while adjusting for other environmental exposures. Similar to our findings, Hystad et al. found no association between PM2.5 and birthweight when including NDVI in the model [20]. Nieuwenhuijsen et al., however, did find a significant association between PM2.5 and birthweight when including other urban exposome exposures in the model. Like our study, which identified neighborhood greenness as one of the most important predictors of differences in birthweight, this study also found the most consistent statistically significant associations between green space exposure and birthweight [21]. Low socioeconomic status (SES) has been linked to poor obstetric outcomes in many studies [48,49,50]. Most studies that assessed the effect of segregation on birth outcomes focused mainly on the effect of racial segregation in the neighborhood rather than economic segregation [51,52,53]. One study, by Vinikoor et al., found that being economically privileged is protective against low birthweight only in predominantly black neighborhoods. The Vinikoor study, however, used a different economic segregation measure and assessed the effect of maternal residence in census tracts with a higher household income than would be expected based on the mother’s education and marital status [54]. Here, we show that residing in a census block-groups with “low income” economic segregation (i.e., high proportion of people living in low income versus high income households) poses a risk for low birthweight deliveries, independent of the mother’s race. Economic dissimilarity was not associated with birthweight in the multivariate model. Unlike the established association between birthweight and exposure to air pollution and SES, evidence on noise, greenness, walkability, and temperature exposures is scarce and results are inconsistent. We found birthweight to be inversely associated with temperature and nighttime noise and positively associated with greenness. Most previous studies found birthweight decreases associated with hot temperatures [18,46], some with cold temperatures [55], and others observed no association at all [8,56]. A recent systematic review found evidence for an association between noise exposure during pregnancy and lower birthweight; however, it concluded that most current studies are of low quality and require better exposure assessment and adjustment for SES and environmental confounders [7]. Several studies have demonstrated health benefits of residence in greener neighborhoods [57]. Although a recent review of residential greenness and birth outcomes found consistent positive associations between neighborhood greenness and birthweight, the associations were weak [16]. We found increases in birthweight associated with greenness exposure in the first and third trimesters of pregnancy, among women who reside in urban neighborhoods. We also found lower birthweight associated with higher walkability. Similar to our findings, a study in Connecticut found that urbanicity was negatively associated with birthweight [5]. The exposures tested in our study can affect birthweight through different mechanisms. Exposure to high temperatures during early stages of pregnancy may cause a reduction in placental weight and umbilical cord flow [4]. High-temperature exposures during the final stages of pregnancy may also affect birthweight, but the mechanism for this association is yet to be clear [4]. Exposure to higher air pollution levels may cause oxidative stress, impairing the transport of oxygen and nutrient to the fetus and affect the intrauterine growth [58,59]. Exposure to high noise levels can alter the mother’s levels of stress hormones, reduce ovarian and uterine blood flow and inhibit fetal growth [60]. Residence in a greener neighborhood can have a beneficial effect on fetal growth through mitigation of harmful environmental exposures (i.e., filter air pollution and provide shade), or through the opportunity for social interactions and physical activity, contributing to better maternal health and lower stress levels [9]. Higher neighborhood walkability can be associated with higher physical activity and therefore a healthier maternal lifestyle and higher birthweight. On the one hand, it can be associated with higher urbanicity, traffic, and pollution levels, and therefore associated with lower birthweight [20]. To the best of our knowledge, there are no studies that assessed the effect of economic residential segregation on birthweight. However, there is a study that examined the effect of racial segregation on birthweight among black mothers, which found a higher probability of low birthweight in more segregated metropolitan areas [61]. Residence in neighborhoods with lower SES can increase the risk for lower birthweight through limited access to housing, healthcare services, healthy nutrition, quality education and more [48,49]. Similarly, we hypothesize that residence in block groups with higher economic segregation, similar to racial segregation, can have adverse effects on birthweight due to higher levels of poverty and crime, higher exposure to stress, fewer options for purchasing healthy food, and limited access to health care [61,62]. The major strength of our study is the inclusion of a large sample of births, and a wide range of environmental predictors of birthweight. Another strength is the combination of two statistical approaches. The elastic net regression, which allowed us to select important predictors of birthweight while accounting for the complex correlations between these multiple exposures. The analysis of health outcomes that are affected by a complex net of exposures requires tools that will help distinguish the contribution and weight of each exposure while accounting for other exposures. As the elastic net model has the ability to distinguish between the redundant exposures, and to identify exposures that are important predictors of birthweight, this method is an important tool in environmental research, where health outcomes are affected by multiple exposures, and where exposures are often correlated. Within the pool of important predictors, we also used WQS regression, and identified the exposures that had higher weight in their contribution to a negative difference in birthweight. These results can be of use to expecting families in their choice of place of residence. No place of residence is free of environmental exposures, and there is always a tradeoff of exposures. This tool, however, can provide insights on exposures that contribute more to differences in birthweight, and help focus the efforts to reduce exposure to nighttime noise and temperature, and increase exposure to greenness. Our study had several limitations. First, we excluded 12% of the records due to missing data. However, this problem is common in studies that use routinely collected data. In addition, as the excluded population was mostly similar to the population included in the analysis, we do not expect it to cause a selection bias. Second, as we assigned exposure based on maternal residence reported at the time of birth, we might have had misclassified exposure for women who changed addresses as conception. It is unlikely that women based their decision regarding the new place of residence on air pollution or temperature exposure. Additionally, it is unlikely that a pregnancy will result in a move to a new residence with a very different socioeconomic environment. We therefore expect the measurement error for these exposures to be non-differential and therefore we do not expect it to bias our results. It is possible, however, that pregnant women will relocate to quieter, greener, and/or less walkable areas. Therefore, for these exposures, the measurement error could potentially bias the results toward the null. Third, we did not have information on other maternal factors that may confound the association between the tested exposures and birth weight. Maternal weight for example was reported to be associated with birthweight [63], as well as with air pollution, neighborhood greenness, and noise exposure [64,65]. Factors such as genetic susceptibility [3] and substance-using mothers [66], although related to birthweight, are not likely to be related to environmental exposures and therefore are not likely to confound the associations in our study. Finally, we could not evaluate trimester-specific exposures to walkability and noise, which were only calculated once during the study period. However, as these land use characteristics likely vary very slowly over time, we do not expect it to bias our results. In addition to the aforementioned limitations, our study is subjected to the limitations of the methods used in our analysis. The elastic net model selects variables with non-zero coefficients, regardless of the effect size and clinical meaning of this effect. In addition, although the elastic net model is capable of identifying important predictors of birthweight, it may have a tendency to incorrectly select a large number of predictors in cases of strong associations with the outcome [67]. However, as we estimated small environmental effect sizes, we do not expect this limitation to affect the variable selection process in our analysis. Lastly, a multi-phase approach in a statistical analysis may result in model overfitting. However, as the aim of the first stage in our statistical approach (i.e., single exposure models) was descriptive and was not used to select variables, this stage was not a part of the selection process and therefore unlikely to result in overfitting. Furthermore, as the elastic net model found all the covariates as predictors of the outcome, overfitting of the final multivariate regression, which included the exact same covariates, is unlikely in our analysis. A major limitation of the WQS regression is that it can only include variables for which the effect on the outcome is in the same direction. We therefore could not include economic residential segregation as an exposure in this analysis. In addition, the WQS regression may lose information on the continuous range of exposure due to the use of quartiles. We therefore added sensitivity analyses, assigning different quantiles to exposures to create the weights. Although the results differed depending on the quantiles used, nighttime noise was consistently the most highly weighted exposure. Last, WQS requires strong assumptions of linear effects and lack of interactions between exposures [43]. 5. Conclusions In conclusion, even after accounting for maternal and neonatal characteristics, we found multiple environmental exposures, built environment exposures, and neighborhood economic indices during pregnancy to be associated with differences in birthweight, emphasizing the complex role of multiple exposures in fetal growth and development. The inclusion of multiple exposures in the model provides avenues for intervention across different exposure dimensions. Acknowledgments We would like to thank Daniel J. Mennitt and Kurt M. Fristrup from the National Park Service, CO, USA, who contributed greatly to this work by providing the noise exposure model used in this analysis. Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations. Supplementary Materials The following are available online at https://www.mdpi.com/1660-4601/17/23/8805/s1, Table S1: Maternal and neonatal characteristics of women who were excluded from the analysis. Table S2: Census area economic segregation and dissimilarity indices. Table S3: A multi-exposure model for the association between birthweight and each environmental exposure, built environment characteristic and economic indices during pregnancy. Table S4: Sensitivity analyses: the association between birthweight and each environmental exposure, built environment characteristic, and economic indexes during pregnancy in a multi-exposure model. Figure S1: The associations between the exposures and birthweight using penalized splines to allow for nonlinear associations: results of a multivariate regression. Figure S2: The mean PM2.5 exposure during pregnancy per census block group, among women who gave birth in Massachusetts between 2001 and 2011. Figure S3: The mean temperature exposure during pregnancy per census block group, among women who gave birth in Massachusetts between 2001 and 2011. Figure S4: The mean greenness exposure during pregnancy per census block group, among women who gave birth in Massachusetts between 2001 and 2011. Figure S5: The mean noise exposure during pregnancy per census block group, among women who gave birth in Massachusetts between 2001 and 2011. Figure S6: The mean walkability index during pregnancy per census block group, among women who gave birth in Massachusetts between 2001 and 2011. Figure S7: The mean economic residential segregation (ERS) exposure during pregnancy per census block group, among women who gave birth in Massachusetts between 2001 and 2011. Figure S8: The mean index of economic dissimilarity (IED) exposure during pregnancy per census block group, among women who gave birth in Massachusetts between 2001 and 2011. Figure S9: The weighted contribution of each of the selected exposures to a negative difference in birthweight: results of a Weighted Quantile Sum (WQS) regression with weights assigned to tertile-scored exposures. Figure S10: The weighted contribution of each of the selected exposures to a negative difference in birthweight: results of a Weighted Quantile Sum (WQS) regression with weights assigned to quintile-scored exposures. Click here for additional data file. Author Contributions Conceptualization, A.Z., M.P.F., J.E.H., and F.L.; Methodology, A.Z. and M.Y.-S.; Formal Analysis, M.Y.-S.; Resources, J.D.S., P.J., K.J.L., K.C.F., and I.K.; Writing—Original Draft Preparation, M.Y.-S.; Writing—Review and Editing, M.P.F., K.J.L., J.E.H., J.D.S., F.L., P.J., K.C.F., I.K., and A.Z.; Supervision, A.Z., M.P.F., and F.L.; Funding Acquisition, A.Z., M.P.F., J.D.S., and F.L. All authors have read and agreed to the published version of the manuscript. Funding This study was supported by the National Institute of Minority Health and Health Disparities grant (P50MD010428); National Institution of Health grants (P30 ES000002, R01 ES024332, and R01 HL150119); and Environmental Protection Agency grants (RD83615601 and RD83587201). Its contents are solely the responsibility of the grantee and do not necessarily represent the official views of the USEPA. Further, USEPA does not endorse the purchase of any commercial products or services mentioned in the publication. Conflicts of Interest The authors declare no conflict of interest. Figure 1 The correlations between the environmental exposures, the built environment characteristics, and the economic indices during pregnancy. Figure 1 shows the correlations between all the exposures. Higher positive correlations are marked with larger blue circles, and higher negative correlations are marked with larger red circles. ERS = economic residential segregation; IED = index of economic dissimilarity; NDVI 1 to 3 = normalized difference vegetation index in the 1st, 2nd, and 3rd trimesters; PM 1 to 3 = fine particulate matter in the 1st, 2nd, and 3rd trimesters; TMP 1 to 3 = temperature in the 1st, 2nd, and 3rd trimesters. Figure 2 The weighted contribution of each of the selected exposures to a negative difference in birthweight: results of a Weighted Quantile Sum (WQS) regression. ijerph-17-08805-t001_Table 1Table 1 Maternal and neonatal characteristics of women who gave birth in Massachusetts between 2001 and 2011. N 640,659 Maternal characteristics Age, mean ± SD 30.1 (6.8) Race, n (%) White 446,053 (69.9%) Black 53,359 (8.5%) Other 141,247 (21.6%) Education, n (%) Less than high school 72,606 (11.3%) High school 161,774 (25.3%) Some college 137,743 (21.5%) College 166,484 (26.0%) More than college 102,052 (15.9%) Government support 213,462 (33.3%) Smoking, n (%) 91,050 (14.2%) Parity > 1, n (%) 292,604 (54.7%) Diabetes mellitus, n (%) 6095 (1.0%) Chronic hypertension, n (%) 8180 (1.3%) Neonatal characteristics Birthweight, mean ± SD (g) 3379.3 ± 533.2 Sex, n (%) Female 312,706 (48.8%) Male 327,953 (51.2%) Gestational age, mean ± SD 39.0 ± 1.6 Census-block SES characteristics Percent poverty 2.5% ± 3.2% Median household income 70,627 ± 35,693 Table 1 shows the maternal, neonatal, and census block group socioeconomic characteristics of the study population. Maternal diabetes and hypertension refer to chronic conditions, not induced by gestation. Gestational age at birth is calculated from the last menstrual period. ijerph-17-08805-t002_Table 2Table 2 Summary statistics of the environmental exposures, built environment characteristics, and economic indices during pregnancy of women who gave birth in Massachusetts between 2001 and 2011 (n = 640,695). Exposure Summary Statistics Range (Min; Max) 25th Percentile 50th Percentile 75th Percentile PM2.5 (µg/m3, Mean ± SD) 1st trimester 10.5 ± 1.9 (3.9; 22.9) 9.0 10.3 11.8 2nd trimester 10.4 ± 1.7 (4.3; 20.5) 9.1 10.2 11.5 3rd trimester 10.4 ± 2.0 (3.9; 23.3) 8.8 10.2 11.7 Temperature (°C, Mean ± SD) 1st trimester 10.9 ± 4.8 (−1.6; 29.7) 6.5 10.8 15.4 2nd trimester 11.2 ± 4.5 (−1.2; 33.2) 7.1 11.2 15.4 3rd trimester 11.4 ± 4.8 (−4.1; 29.6) 7.1 11.7 15.8 Greenness (NDVI, Mean ± SD) 1st trimester 0.5 ± 0.2 (−0.2; 0.9) 0.3 0.5 0.7 2nd trimester 0.5 ± 0.2 (−0.2; 0.9) 0.3 0.5 0.7 3rd trimester 0.5 ± 0.2 (−0.2; 0.9) 0.3 0.5 0.6 Walkability,( Mean ± SD) 1.1 ± 1.9 (−2.8; 16.9) −0.1 0.7 2.0 ERS, n (%) −1 to −0.6 16,729 (2.6) −0.6 to −0.2 121,562 (19.0) −0.2 to 0.2 251,704 (39.3) 0.2 to 0.6 214,642 (33.5) More than 0.6 36,022 (5.6) IED (%, Mean ± SD) 16.9 ± 11.2 (0.0; 51.3) 8.4 14.8 23.3 Median Nighttime Noise (dB, Mean ± SD) 43.5 ± 2.9 (29.7; 53.5) 41.5 42.9 45.6 ijerph-17-08805-t003_Table 3Table 3 The association between birthweight and each environmental exposure, built environment exposure, and economic indices during pregnancy in single, and multi-exposure models (n = 640,695). IQR Difference in Weight (g per IQR Increase (95% CI)) Exposure Single Exposure Models Multi-Exposure Models PM2.5, 1st trimester (µg/m3) 2.8 −5.4 (−7.53; −3.28) * 0.96 (−1.65; 3.56) PM2.5, 2nd trimester (µg/m3) 2.4 −9.4 (−11.59; −7.21) * −0.83 (−3.34; 1.69) PM2.5, 3rd trimester (µg/m3) 2.8 −7.76 (−10.07; −5.46) * −0.44 (−3.16; 2.29) Temperature, 1st trimester (°C) 8.8 0.59 (−3.3; 4.48) −10.27 (−17.68; −2.87) * Temperature, 2nd trimester (°C) 8.3 −10.6 (−14.58; −6.62) * −4.85 (−9.08; −0.63) * Temperature, 3rd trimester (°C) 8.6 −15.75 (−19.84; −11.66) * −19.3 (−27.04; −11.57) * NDVI, 1st trimester 0.3 20.89 (18.55; 23.23) * 5.71 (2.97; 8.45) * NDVI, 2nd trimester 0.3 17.94 (15.7; 20.19) * 2.53 (−0.17; 5.23) NDVI, 3rd trimester 0.3 20.03 (17.79; 22.26) * 5.87 (3.22; 8.53) * Walkability 2.2 −5.4 (−7.53; −3.28) * −5.63 (−7.24; −4.02) * ERS −1 to −0.6 Reference Reference −0.6 to −0.2 21.39 (14.28; 28.5) * 15.99 (8.85; 23.14) * −0.2 to 0.2 42.68 (35.68; 49.67) * 28.64 (21.52; 35.77) * 0.2 to 0.6 53.39 (46.19; 60.6) * 32.12 (24.61; 39.63) * 0.6+ 44.88 (36.48; 53.27) * 19.46 (10.58; 28.34) * IED 14.9 −2.14 (−3.58; −0.71) * −0.32 (−1.83; 1.18) Noise levels, (dB) 4.1 −16.88 (−18.49; −15.27) * −5.63 (−7.52; −3.73) * ERS = Economic Residential Segregation; IED = Index of Economic Dissimilarity; NDVI = normalized difference vegetation index. * p < 0.05. ==== Refs References 1. Goldenberg R.L. Culhane J.F. Low birth weight in the United States Am. J. Clin. Nutr. 2007 85 584S 590S 10.1093/ajcn/85.2.584S 17284760 2. Belbasis L. Savvidou M.D. Kanu C. Evangelou E. Tzoulaki I. Birth weight in relation to health and disease in later life: An umbrella review of systematic reviews and meta-analyses BMC Med. 2016 14 147 10.1186/s12916-016-0692-5 27677312 3. Mann N. Birth weight symposium Arch. Dis. Child. 2002 86 F2 10.1136/fn.86.1.F2 4. Wells J.C.K. Thermal environment and human birth weight J. Theor. Biol. 2002 214 413 425 10.1006/jtbi.2001.2465 11846599 5. Ebisu K. Holford T.R. Bell M.L. Association between greenness, urbanicity, and birth weight Sci. Total Environ. 2016 542 750 756 10.1016/j.scitotenv.2015.10.111 26546769 6. Smith R.B. Fecht D. Gulliver J. Beevers S.D. Dajnak D. Blangiardo M. Ghosh R.E. Hansell A.L. Kelly F.J. Anderson H.R. Impact of London’s road traffic air and noise pollution on birth weight: Retrospective population based cohort study BMJ Br. Med. J. 2017 359 j5299 10.1136/bmj.j5299 29208602 7. Nieuwenhuijsen M.J. Ristovska G. Dadvand P. WHO Environmental Noise Guidelines for the European Region: A Systematic Review on Environmental Noise and Adverse Birth Outcomes Int. J. Environ. Res. Public Health 2017 14 1252 10.3390/ijerph14101252 29048350 8. Diaz J. Arroyo V. Ortiz C. Carmona R.O. Linares C. Effect of Environmental Factors on Low Weight in Non-Premature Births: A Time Series Analysis PLoS ONE 2016 11 e0164741 27788159 9. Fong K.C. Kloog I. Coull B.A. Koutrakis P. Laden F. Schwartz J.D. James P. Residential Greenness and Birthweight in the State of Massachusetts, USA Int. J. Environ. Res. Public Health 2018 15 1248 10.3390/ijerph15061248 10. Yitshak-Sade M. Novack L. Landau D. Kloog I. Sarov B. Hershkovitz R. Karakis I. Relationship of ambient air pollutants and hazardous household factors with birth weight among Bedouin-Arabs Chemosphere 2016 160 314 322 10.1016/j.chemosphere.2016.06.104 27391054 11. Li X.Y. Huang S.Q. Jiao A.Q. Yang X.H. Yun J.F. Wang Y.X. Xue X.W. Chu Y.Y. Liu F.F. Liu Y.S. Association between ambient fine particulate matter and preterm birth or term low birth weight: An updated systematic review and meta-analysis Environ. Pollut. 2017 227 596 605 10.1016/j.envpol.2017.03.055 28457735 12. Kloog I. Novack L. Erez O. Just A.C. Raz R. Associations between ambient air temperature, low birth weight and small for gestational age in term neonates in southern Israel Environ. Health 2018 17 9 10.1186/s12940-018-0420-z 29347948 13. Gehring U. Tamburic L. Sbihi H. Davies H.W. Brauer M. Impact of Noise and Air Pollution on Pregnancy Outcomes Epidemiology 2014 25 351 358 10.1097/EDE.0000000000000073 24595395 14. Li S. Wang J. Xu Z. Wang X. Xu G. Zhang J. Shen X. Tong S. Exploring associations of maternal exposure to ambient temperature with duration of gestation and birth weight: A prospective study BMC Pregnancy Childbirth 2018 18 513 10.1186/s12884-018-2100-y 30594173 15. Sun S. Spangler K.R. Weinberger K.R. Yanosky J.D. Braun J.M. Wellenius G.A. Ambient Temperature and Markers of Fetal Growth: A Retrospective Observational Study of 29 Million U.S. Singleton Births Envion. Health Perspect. 2019 127 67005 10.1289/EHP4648 16. Banay R.F. Bezold C.P. James P. Hart J.E. Laden F. Residential greenness: Current perspectives on its impact on maternal health and pregnancy outcomes Int. J. Womens Health 2017 9 133 144 10.2147/IJWH.S125358 28280395 17. Zhu X.X. Liu Y. Chen Y.Y. Yao C.J. Che Z. Cao J.Y. Maternal exposure to fine particulate matter (PM 2.5) and pregnancy outcomes: A meta-analysis Environ. Sci. Pollut. Res. 2015 22 3383 3396 10.1007/s11356-014-3458-7 18. Zhang Y.Q. Yu C.H. Wang L. Temperature exposure during pregnancy and birth outcomes: An updated systematic review of epidemiological evidence Environ. Pollut. 2017 225 700 712 10.1016/j.envpol.2017.02.066 28284544 19. Wild C.P. The exposome: From concept to utility Int. J. Epidemiol. 2012 41 24 32 10.1093/ije/dyr236 22296988 20. Hystad P. Davies H.W. Frank L. Van Loon J. Gehring U. Tamburic L. Brauer M. Residential Greenness and Birth Outcomes: Evaluating the Influence of Spatially Correlated Built-Environment Factors Environ. Health Perspect. 2014 122 1095 1102 10.1289/ehp.1308049 25014041 21. Nieuwenhuijsen M.J. Agier L. Basagana X. Urquiza J. Tamayo-Uria I. Giorgis-Allemand L. Robinson O. Siroux V. Maitre L. de Castro M. Influence of the Urban Exposome on Birth Weight Environ. Health Perspect. 2019 127 047007 10.1289/EHP3971 31009264 22. Agier L. Basagaña X. Hernandez-Ferrer C. Maitre L. Tamayo Uria I. Urquiza J. Andrusaityte S. Casas M. de Castro M. Cequier E. Association between the pregnancy exposome and fetal growth Int. J. Epidemiol. 2020 49 572 586 10.1093/ije/dyaa017 32167557 23. United States Central Bureau of Statistics (USCBS) Census United States Department of Commerce Washington, DC, USA 2010 Available online: https://www.census.gov/en.html (accessed on 15 March 2019) 24. Ali M.K. Quarashi M.S.A. Sultana S. Rahman M.Z. Observation of Birth Weight of Babies in relation on maternal age, parity and gestational age in Tertiary Level Hospital Bangladesh J. Med Sci. 2020 19 291 295 10.3329/bjms.v19i2.45010 25. Dadvand P. Ostro B. Figueras F. Foraster M. Basagana X. Valentin A. Martinez D. Beelen R. Cirach M. Hoek G. Residential Proximity to Major Roads and Term Low Birth Weight the Roles of Air Pollution, Heat, Noise, and Road-Adjacent Trees Epidemiology 2014 25 518 525 10.1097/EDE.0000000000000107 24787556 26. Sapkota A. Chelikowsky A.P. Nachman K.E. Cohen A.J. Ritz B. Exposure to particulate matter and adverse birth outcomes: A comprehensive review and meta-analysis Air Qual. Atmos. Health 2012 5 369 381 10.1007/s11869-010-0106-3 27. Kloog I. Chudnovsky A.A. Just A.C. Nordio F. Koutrakis P. Coull B.A. Lyapustin A. Wang Y.J. Schwartz J. A new hybrid spatio-temporal model for estimating daily multi-year PM2.5 concentrations across northeastern USA using high resolution aerosol optical depth data Atmos. Environ. 2014 95 581 590 10.1016/j.atmosenv.2014.07.014 28. Kloog I. Nordio F. Coull B.A. Schwartz J. Predicting spatiotemporal mean air temperature using MODIS satellite surface temperature measurements across the Northeastern USA Remote Sens. Environ. 2014 150 132 139 10.1016/j.rse.2014.04.024 29. Casey J.A. James P. Rudolph K.E. Wu C.D. Schwartz B.S. Greenness and Birth Outcomes in a Range of Pennsylvania Communities Int. J. Environ. Res. Public Health 2016 13 311 10.3390/ijerph13030311 30. Lane K.J. Stokes E.C. Seto K.C. Thanikachalam S. Thanikachalam M. Bell M.L. Associations between Greenness, Impervious Surface Area, and Nighttime Lights on Biomarkers of Vascular Aging in Chennai, India Environ. Health Perspect. 2017 125 087003 10.1289/EHP541 28886599 31. James P. Hart J.E. Banay R.F. Laden F. Exposure to Greenness and Mortality in a Nationwide Prospective Cohort Study of Women Environ. Health Perspect. 2016 124 1344 1352 10.1289/ehp.1510363 27074702 32. Environmental Protection Agency (EPA) Smart Location Database. Technical Documentation and User Guide 2014 Available online: https://www.epa.gov/smartgrowth/smart-location-mapping (accessed on 10 April 2015) 33. James P. Kioumourtzoglou M.A. Hart J.E. Banay R.E. Kloog I. Laden F. Interrelationships between Walkability, Air Pollution, Greenness, and Body Mass Index Epidemiology 2017 28 780 788 10.1097/EDE.0000000000000724 28767514 34. Mennitt D.J. Fristrup K.M. Influential factors and spatiotemporal patterns of environmental sound levels in the contiguous United States Noise Control Eng. J. 2016 64 342 353 10.3397/1/376384 35. Casey J.A. Morello-Frosch R. Mennitt D.J. Fristrup K. Ogburn E.L. James P. Race/Ethnicity, Socioeconomic Status, Residential Segregation, and Spatial Variation in Noise Exposure in the Contiguous United States Environ. Health Perspect. 2017 125 077017 10.1289/EHP898 28749369 36. Krieger N. Waterman P.D. Batra N. Murphy J. Dooley D. Shah S. Measures of Local Segregation for Monitoring Health Inequities by Local Health Departments AJPH Methods 2017 107 903 906 10.2105/AJPH.2017.303713 37. Massay D.S. The prodigal paradigm returns: Ecology comes back to sociology Does It Take a Village? Community Effects on Children, Adolescents, and Families Lawrence Erlbaum Associates Mahwah, NJ, USA 2001 41 48 38. Massey D.S. Reflections on the Dimensions of Segregation Soc. Forces 2012 91 39 43 10.1093/sf/sos118 24920864 39. Kloog I. Air pollution, ambient temperature, green space and preterm birth Curr. Opin. Pediatr. 2019 31 237 243 10.1097/MOP.0000000000000736 30640892 40. Tibshirani R. Regression shrinkage and selection via the lasso: A retrospective J. R. Stat. Soc. Ser. B Stat. Methodol. 2011 73 273 282 10.1111/j.1467-9868.2011.00771.x 41. Fan J.Q. Li R.Z. Variable selection via nonconcave penalized likelihood and its oracle properties J. Am. Stat. Assoc. 2001 96 1348 1360 10.1198/016214501753382273 42. Zou H. The adaptive lasso and its oracle properties J. Am. Stat. Assoc. 2006 101 1418 1429 10.1198/016214506000000735 43. Carrico C. Gennings C. Wheeler D.C. Factor-Litvak P. Characterization of Weighted Quantile Sum Regression for Highly Correlated Data in a Risk Analysis Setting J. Agric. Biol. Environ. Stat. 2015 20 100 120 10.1007/s13253-014-0180-3 30505142 44. Fleischer N.L. Merialdi M. van Donkelaar A. Vadillo-Ortega F. Martin R.V. Betran A.P. Souza J.P. O’Neill M.S. Outdoor Air Pollution, Preterm Birth, and Low Birth Weight: Analysis of the World Health Organization Global Survey on Maternal and Perinatal Health Environ. Health Perspect. 2014 122 425 430 10.1289/ehp.1306837 24508912 45. Hyder A. Lee H.J. Ebisu K. Koutrakis P. Belanger K. Bell M.L. PM2.5 Exposure and Birth Outcomes Use of Satellite- and Monitor-Based Data Epidemiology 2014 25 58 67 10.1097/EDE.0000000000000027 24240652 46. Kloog I. Melly S.J. Coull B.A. Nordio F. Schwartz J.D. Using Satellite-Based Spatiotemporal Resolved Air Temperature Exposure to Study the Association between Ambient Air Temperature and Birth Outcomes in Massachusetts Environ. Health Perspect. 2015 123 1053 1058 10.1289/ehp.1308075 25850104 47. Kuehn L. McCormick S. Heat Exposure and Maternal Health in the Face of Climate Change Int. J. Environ. Res. Public Health 2017 14 853 10.3390/ijerph14080853 48. Morello-Frosch R. Lopez R. The riskscape and the color line: Examining the role of segregation in environmental health disparities Environ. Res. 2006 102 181 196 10.1016/j.envres.2006.05.007 16828737 49. Shmool J.L.C. Kubzansky L.D. Newman O.D. Spengler J. Shepard P. Clougherty J.E. Social stressors and air pollution across New York City communities: A spatial approach for assessing correlations among multiple exposures Environ. Health 2014 13 91 10.1186/1476-069X-13-91 25374310 50. Gavin A.R. Thompson E. Rue T. Guo Y.Q. Maternal Early Life Risk Factors for Offspring Birth Weight: Findings from the Add Health Study Prev. Sci. 2012 13 162 172 10.1007/s11121-011-0253-2 21986991 51. Debbink M.P. Bader M.D.M. Racial Residential Segregation and Low Birth Weight in Michigan’s Metropolitan Areas Am. J. Public Health 2011 101 1714 1720 10.2105/AJPH.2011.300152 21778487 52. Roberts E.M. Neighborhood social environments and the distribution of low birthweight in Chicago Am. J. Public Health 1997 87 597 603 10.2105/AJPH.87.4.597 9146438 53. Bell J.F. Zimmerman F.J. Almgren G.R. Mayer J.D. Huebner C.E. Birth outcomes among urban African-American women: A multilevel analysis of the role of racial residential segregation Soc. Sci. Med. 2006 63 3030 3045 10.1016/j.socscimed.2006.08.011 16997438 54. Vinikoor L.C. Kaufman J.S. MacLehose R.F. Laraia B.A. Effects of racial density and income incongruity on pregnancy outcomes in less segregated communities Soc. Sci. Med. 2008 66 255 259 10.1016/j.socscimed.2007.08.016 17920176 55. Elter K. Ay E. Uyar E. Kavak Z.N. Exposure to low outdoor temperature in the midtrimester is associated with low birth weight Aust. N. Z. J. Obstet. Gynaecol. 2004 44 553 557 10.1111/j.1479-828X.2004.00314.x 15598296 56. Tustin K. Gross J. Hayne H. Maternal exposure to first-trimester sunshine is associated with increased birth weight in human infants Dev. Psychobiol. 2004 45 221 230 10.1002/dev.20030 15549686 57. Kondo M.C. Fluehr J.M. McKeon T. Branas C.C. Urban Green Space and Its Impact on Human Health Int. J. Environ. Res. Public Health 2018 15 445 10.3390/ijerph15030445 29510520 58. Proietti E. Roosli M. Frey U. Latzin P. Air Pollution during Pregnancy and Neonatal Outcome: A Review J. Aerosol Med. Pulm. Drug Deliv. 2013 26 9 23 10.1089/jamp.2011.0932 22856675 59. Schlesinger R.B. Kunzli N. Hidy G.M. Gotschi T. Jerrett M. The health relevance of ambient particulate matter characteristics: Coherence of toxicological and epidemiological inferences Inhal. Toxicol. 2006 18 95 125 10.1080/08958370500306016 16393926 60. Ristovska G. Laszlo H.E. Hansell A.L. Reproductive Outcomes Associated with Noise Exposure—A Systematic Review of the Literature Int. J. Environ. Res. Public Health 2014 11 7931 7952 10.3390/ijerph110807931 25101773 61. Ellen Gould I. Is Segregation Bad for Your Health? The Case of Low Birth Weight Papers on Urban Affairs Press BI Washington, DC, USA 1999 203 229 62. Kramer M.R. Hogue C.R. Is Segregation Bad for Your Health? Epidemiol. Rev. 2009 31 178 194 10.1093/epirev/mxp001 19465747 63. Rode L. Hegaard H.K. Kjtaergaard H. Moller L.F. Tabor A. Ottesen B. Association between maternal weight gain and birth weight Obstet. Gynecol. 2007 109 1309 1315 10.1097/01.AOG.0000266556.69952.de 17540802 64. Sallis J.F. Floyd M.F. Rodriguez D.A. Saelens B.E. Role of Built Environments in Physical Activity, Obesity, and Cardiovascular Disease Circulation 2012 125 729 737 10.1161/CIRCULATIONAHA.110.969022 22311885 65. Lachowycz K. Jones A.P. Greenspace and obesity: A systematic review of the evidence Obes. Rev. 2011 12 e183 e189 21348919 66. Oni H.T. Khan M.N. Abdel-Latif M. Buultjens M. Islam M.M. Short-term health outcomes of newborn infants of substance-using mothers in Australia and New Zealand: A systematic review J. Obstet. Gynaecol. Res. 2019 45 1783 1795 10.1111/jog.14051 31313404 67. Czarnota J. Gennings C. Wheeler D.C. Assessment of Weighted Quantile Sum Regression for Modeling Chemical Mixtures and Cancer Risk Cancer Inform. 2015 14 159 171 26005323