Skip to main content

Air pollution, neighbourhood and maternal-level factors modify the effect of smoking on birth weight: a multilevel analysis in British Columbia, Canada



Maternal smoking during pregnancy negatively impacts fetal growth, but the effect is not homogenous across the population. We sought to determine how the relationship between cigarette use and fetal growth is modified by the social and physical environment.


Birth records with covariates were obtained from the BC Perinatal Database Registry (N = 232,291). Maternal smoking status was self-reported as the number of cigarettes smoked per day usually at the first prenatal care visit. Census dissemination areas (DAs) were used as neighbourhood-level units and linked to individual births using residential postal codes to assign exposure to particulate air pollution (PM2.5) and neighbourhood-level attributes such as socioeconomic status (SES), proportion of post-secondary education, immigrant density and living in a rural place. Random coefficient models were used with cigarettes/day modeled with a random slope to estimate its between-DA variability and test cross-level interactions with the neighbourhood-level variables on continuous birth weight.


A significant negative and non-linear association was found between maternal smoking and birth weight. There was significant between-DA intercept variability in birth weight as well as between-DA slope variability of maternal smoking on birth weight of which 68 and 30 % respectively was explained with the inclusion of DA-level variables and their cross-level interactions. High DA-level SES had a strong positive association with birth weight but the effect was moderated with increased cigarettes/day. Conversely, heavy smokers showed the largest increases in birth weight with rising neighbourhood education levels. Increased levels of PM2.5 and immigrant density were negatively associated with birth weight, but showed positive interactions with increased levels of smoking. Older maternal age and suspected drug or alcohol use both had negative interactions with increased levels of maternal smoking.


Maternal smoking had a negative and non-linear dose-response association with birth weight which was highly variable between neighbourhoods and evidence of effect modification with neighbourhood-level factors. These results suggest that focusing exclusively on individual behaviours may have limited success in improving outcomes without addressing the contextual influences at the neighbourhood-level. Further studies are needed to corroborate our findings and to understand how neighbourhood-level attributes interact with smoking to affect birth outcomes.

Peer Review reports


Smoking during pregnancy is a modifiable risk factor associated with adverse birth outcomes and may impart long-term health consequences [13]. This relationship however is confounded by the presence of many other risk factors, including maternal age, education, alcohol or drug use [46]. Furthermore, it’s been shown that these individual-level risk factors have a dose-response association with the level of smoking, with a distinction between heavy smokers (greater than 10 cigarettes per day) and moderate or light smokers [4]. For example, while the prevalence of smoking during pregnancy decreases with increasing maternal age, the level of smoking is heavier among the older mothers who do smoke. As a result, the effect of smoking on birth weight has been shown to be modified by maternal age or other correlated risk factors [7, 8]. Similarly, neighbourhood-level factors might directly or indirectly modify the effect of smoking on birth weight such as neighbourhood deprivation or levels of particulate air pollution [911].

Exposure to the fine fraction of particulate matter (PM2.5, particles with aerodynamic diameter ≤ 2.5 μm) has shown to be a consistent risk factor associated with reduced birth weight [12]. The complex mixture of PM2.5 includes elemental and organic carbon compounds, metals and gases that stem predominantly from vehicle exhaust, residential heating and industrial emissions [13]. The mechanisms by which PM2.5 and its constituents adversely affect the reproductive system are not fully understood; however, evidence supports the potential for a shared mode of developmental toxicity with tobacco smoke exposure [1416]. With similar chemical components, both PM2.5 and tobacco smoke penetrate deep into pulmonary alveolar tissues and translocate to extrapulmonary tissues causing systemic cardiovascular and immunological alterations, including platelet activation, coagulation, endothelial dysfunction, DNA damage and mutagenesis [13, 16, 17].

Low SES remains one of the most robust predictors of adverse pregnancy outcomes such as fetal growth restriction despite universal health care programs in Canada and Europe [9, 18, 19]. The society-level determinants such as poverty, poor education, income inequality and social discrimination and marginalization act indirectly on the placenta and fetus through the promotion of ‘downstream’ or mediating exposures, stresses and behaviours [20, 21]. In studies of cardiovascular disease, neighbourhood-level factors were associated with increased levels of smoking and other risk factors such as obesity, lack of exercise, lower health knowledge and lower positive behaviour changes [11, 22]. These epidemiological observations have been shown with the use of multilevel statistical models capable of separating the individual-level effects from the context of their social and physical environments [23]. The use of multilevel models in perinatal epidemiology has uncovered neighbourhood-level factors that interact with maternal-level risk factors to either buffer or mediate adverse birth outcomes [21, 24].

We present a multilevel cross-sectional analysis of birth registry data in British Columbia, Canada (population 4.6 million) to investigate neighbourhood-level differences in the effect of cigarette smoking during pregnancy on birth weight and to quantify the degree to which individual and neighbourhood-level variables explain any observed differences. Specifically, we sought to determine whether exposure to PM2.5 and living in low SES neighbourhoods explain between-neighbourhood differences in the effect of maternal smoking on birth weight. We also examine whether these neighbourhood-level factors modify the direct effect of maternal smoking on birth weight. Birth weight is among the most important factors affecting neonatal mortality and is a significant determinant of post-neonatal infant mortality and childhood morbidity [25]. Understanding the underlying individual and interactive effects of exposures on birth weight is crucial for effective community planning and strategic interventions to improving reproductive health outcomes.


This was a population-based cross-sectional study of singleton births in British Columbia from 2001 to 2006 (N = 237,470). Data from the BC Perinatal Database Registry were provided by Perinatal Services BC (PSBC), and included information on individual-level maternal-infant health status and outcomes, reproductive history, socio-demographics, risk factors, and residential postal codes. The Registry accounts for nearly 100 % of births and stillbirths in BC of at least 20 weeks gestation or at least 500 g birth weight. Research data access is provided by a Partnership Accord /Memorandum of Agreement between all BC Health Authorities and PSBC through the Freedom of Information and Privacy Protection Act [26]. Research ethics board approval was granted by the University of Victoria (protocol #11-043).

The outcome variable was continuous birth weight of singleton births. Included were all births (stillbirth and live) for gestational ages of 20 to 42 weeks. Excluded birth records included: out-of-province and invalid postal codes (n = 1096), non-viable births prior to 20 weeks gestation or less than 500 g (n = 15), and the list-wise deletion of births missing important data including: cigarettes smoked per day (cigarettes/day, n = 2510), PM2.5 (n = 1512), birth weight (n = 46). Table 1 provides the full list of covariates used along with their summary statistics. All continuous independent variables, except cigarettes/day, were grand-mean centred and standardized to ease interpretation and aid model convergence. Thus, a value of zero represents the transformed variable’s mean and reference value and has a standard deviation equal to one. The variable cigarettes/day was kept un-transformed since the value zero (i.e. non-smokers) was the desired reference level. Smoking levels were capped at 20 cigarettes/day with higher values assigned a value of 20 to stabilize the distribution tail (n = 245, min.21 max.80). Two variables indicating the use of alcohol or drugs (prescription, non-prescription, illicit) to be a risk factor in pregnancy as identified by a physician were combined into a single dichotomous variable.

Table 1 Descriptive statistics# for individual (Level-1) and DA (Level-2) covariates on term birth weight

Birth records were geocoded based on the latitude-longitude coordinate of the mother’s residential postal code at the time of delivery using GeoRef by DMTI [27]. Birth records were then linked to their corresponding census dissemination area (DA) by performing a point-in-polygon spatial join procedure in ArcGIS 10.2 [28]. DAs represent the smallest geographical unit for which census data are available with a spatial coverage ranging between 200 and 800 people depending on the level of urban development. While DAs do not necessarily represent existing neighbourhood communities [29], they can act as proxies for a general catchment area of personal home-life activities [21, 30]. Birth records were identified as being either rural or urban using the Statistics Canada Metropolitan Influence Zone (MIZ) codes which are based on commuting flows of small towns into larger cities and metropolitan areas [31].

PM2.5 exposure was estimated using a national land-use regression (LUR) model developed to estimate PM2.5 at the census street block-face level [32]. The model used a number of predictors including satellite measures, proximity to major roads and industry to account for 46 % of the variability in measured annual PM2.5 concentrations. Individual birth records were related to the block-face point estimates using a nearest-point procedure in ArcGIS10.2. Street block-face point estimates were related to individual birth records using a nearest-point procedure in ArcGIS10.2 and then aggregated to their DA-level mean to represent an area-level air pollution variable on individual births.

Three related but independent datasets all based on the 2006 Statistics Canada national census were used to represent the DA-level SES and demographic data. The first was a Canadian SES index (SESi) developed by Chan et al. which provides a measure of overall socioeconomic neighbourhood well-being [33]. The second was the proportion of population over 15 with any post-secondary education, including college, trades, or university representing higher DA-level education attainment levels. The third was the proportion of continental Asian immigrants by DA. It’s been shown in BC and elsewhere that healthy babies from Asian and South Asian backgrounds are constitutionally smaller compared to Caucasian babies [34, 35]. Asian and South Asian ethnicities are well-represented throughout BC but particularly in concentrated pockets throughout the major urban center of Metro Vancouver where levels of PM2.5 are also high and could therefore confound any PM2.5 effect. Furthermore, concentrated ethnic communities may impart buffering mechanisms through enhanced social interactions and support networks [21, 24]. A sequential regression technique was used to remove the collinearity between sets of DA-level variables [36]. Here, immigrant density was regressed against SESi and PM2.5 with the saved residuals representing the uncorrelated and independent contribution of immigrant density on birth weight freed from its collinearity with SESi and PM2.5 (r = −0.62 and 0.53 respectively). This method was repeated for SESi and education (r = 0.26) creating a residual immigrant density and residual education variable. The education and immigrant data were obtained by access to ABACUS via the Data Liberation Initiative [37].

Imputation for missing SES, education and immigrant density values was performed in order to avoid data loss of rural DAs with low population counts. Taking advantage of the nested hierarchical structure of the administrative census and health boundaries, the mean SESi value for a larger encompassing census subdivision (CSD) or local health area (LHA) was imputed for a nested DA with a missing value. There were 1441 values imputed in 52 DAs for SESi (0.6 % of final N, 0.8 % of DAs), and 3170 values imputed in 108 DAs for both education and immigrant density (1.4 % of final N, 1.7 % of DAs). Sensitivity analyses were performed using only the non-imputed data.

Hierarchical (multilevel) linear regression models were used to test our research questions, thereby accounting for the clustering, or non-independence, of individuals (level-1) belonging to a given DA neighbourhood (level-2). The multilevel model allows the intercept and slope to act as random parameters having between-area (DA) variability from an overall (BC-wide) mean intercept and slope. Therefore each DA has its own intercept and slope in which their variability from the overall mean intercept and slope can be investigated with the addition of individual (level-1) and DA-level (level-2) variables and their interactions [38]. We followed a bottom-up approach to model building to quantify the explained proportional change in variance (PCV) with the addition of sets of variables, the multilevel model equivalent to an R 2 [23]. We started with the empty (Null) random intercept model without any independent variables in which birth weight is only a function of the mother’s residential DA. The presence of significant random intercept variance indicates there are unexplained differences between neighbourhood means of birth weight. The proportion of the total variance in birth weight that arises due to neighbourhood differences can be quantified by computing the intra-class correlation (ICC) which represents the degree of clustering of individual birth weight within neighbourhoods [23].

The Null model was followed by Model that included the individual-level covariates as well as the addition of a random slope for the continuous variable of maternal smoking (cigarettes/day, self-reported at the first prenatal visit). By allowing cigarettes/day to be random, the mean within-DA effect of maternal smoking is allowed to differ between DAs. The presence of a significant random slope indicates that its effect on birth weight is not constant (or equal) for all DAs; that is, there are important unexplained differences between the within-DA group effects of maternal smoking on birth weight. Subsequent models included the DA-level variables along with cross-level interactions to assess their fixed effects on birth weight but to also determine if their inclusion addresses any unexplained slope variance. Several models were tested using the Akaike Information Criterion (AIC) to evaluate model performance. We report the results of three models to compare the degree of change between the level-1 and level-2 homogeneous (non-interaction) models and a model with effect-measure variation. All statistical analyses were conducted in Stata 13IC [39].


After exclusions, the final dataset included 232,291 singleton (live and stillborn) births located in 6338 neighbourhood DAs (min. = 1, max. = 782, avg. = 37). Table 1 summarizes the untransformed individual and neighbourhood covariates (non-centered, non-standardized). The prevalence of maternal smoking in this population was 10.3 % (n = 23,836) with an average of 7.5 cigarettes/day among smokers. Table 2 reports the adjusted coefficients for the individual and DA-level covariate fixed effects on continuous birth weight (Model 1 and 2). Model 1 was a level-1 model that included only the maternal-level covariates. The relationship between birth weight and cigarettes/day was found to be non-linear and was best modeled using a quadratic term indicating a subdued dose-response with increasing exposure (Fig. 1). Model 2 added the DA-level variables. Their fixed effects show that DAs with higher SES and higher proportion of post-secondary education were significantly associated with higher birth weights; whereas DAs with increased levels of PM2.5, higher Asian immigrant density and rural DAs were all significantly associated with lower birth weights. Season of birth (fall or winter) was also significantly association with reduced birth weight. The results in Table 2 represent the fixed effects from homogeneous models (i.e. those without any modeled heterogeneity of the effect measure for maternal smoking).

Table 2 Adjusted fixed effects for level-1 and level-2 covariates on continuous term birth weight
Fig. 1
figure 1

Adjusted Predicted Effects of Maternal Smoking on Birth Weight. Predicted effects of maternal smoking (cigarettes/day) on birth weight with 95 % confidence intervals are conditional on model covariates included in Model 3. Black vertical lines represent the frequency distribution of cigarettes/day (non-smokers, 0 cigarettes/day, have been omitted for display purposes)

Model 3 tested interactions with cigarettes/day including cross-level (level-1 by level2) and level-1 by level-1 interactions to explain the between-DA random intercept and random slope variability. The model results are presented in Table 3 including the main effects as well as the interaction effects with cigarettes/day. The degree of heterogeneity across levels of maternal smoking modified by the DA-level contextual factors is graphically presented in Fig. 2. The five graphs show the predicted conditional fixed effects of SESi, education, PM2.5, Asian immigrant density and rural residence on birth weight and their interactions with specified levels of maternal smoking (Fig. 2ae respectively). For example, Table 3 and Fig. 2a show that higher SES has a significant positive association with birth weight but is less pronounced with increased levels of maternal smoking whereby very heavy smokers (≥20 cigarettes/day) do not incur any benefit of higher SES. Conversely, very heavy smokers showed the greatest gains in birth weight in DAs with higher proportions of post-secondary educated people (Fig. 2b). Recall that the higher education variable was an uncorrelated residual variable independent of SESi, and therefore these observed associations are in addition to the education-related effect captured by SESi.

Table 3 Adjusted individual and DA-level fixed effects on continuous birth weight and their modification by maternal smoking (Model 3)
Fig. 2
figure 2

Adjusted Predicted Effects of Maternal Smoking on Birth Weight across DA-level Factors. a Socioeconomic Status Index (SESi) b Proportion of Population with Post-secondary Education c Particulate Matter ≤ 2.5 μm (PM2.5) d Asian Immigrant Density e Rural Residence. Predicted effects on birth weight with 95 % confidence intervals are conditional on model covariates included in Model 3. Black vertical lines represent the frequency distribution of the variable on the x-axis (except Fig. 2e which shows the frequency distribution of cigarettes/day)

Increasing PM2.5 levels had a significant non-linear association with reduced birth weight; however, it showed a positive interaction with maternal smoking such that the effect of increased smoking on birth weight was attenuated in DAs with higher levels of PM2.5 (Fig. 2c). Similarly, higher Asian immigrant density was significantly associated with lower birth weights but had a positive interaction with increased cigarette use demonstrating a protective effect of higher immigrant density DAs (Fig. 2d). Rural DAs had a significant negative interaction with maternal smoking indicating a further reduction in birth weight with increased cigarette use among rural residents (Fig. 2e).

Two level-1 interactions with maternal smoking were significant, maternal age and suspected drug or alcohol use. The predicted conditional marginal effects of these two interactions are show in Fig. 3a and 3b respectively indicating that the reduction of birth weight among heavier smokers is exasperated by older maternal age and those suspected of drug or alcohol use. A variable for neighbourhood-level smoking (DA-average cigarettes/day) was created and tested in models along with a cross-level interaction with maternal-level cigarettes/day but neither parameters were significant nor explained any additional variability.

Fig. 3
figure 3

Adjusted Predicted Effects of Maternal Smoking on Birth Weight across Maternal-level Factors. a Maternal age b Suspected Drug or Alcohol Use. Predicted effects on birth weight with 95 % confidence intervals are conditional on model covariates included in Model 3. Black vertical lines represent the frequency distribution of the variable on the x-axis (except Fig. 3b which shows the frequency distribution of cigarettes/day)

The random effects, the explained proportional change in variance (PVC), and model diagnostics are presented in Table 4. The unadjusted ICC for the Null random-intercept model was 0.019, indicating that 1.9 % of the total residual differences in birth weight are attributable to DA-level contextual factors. However, the ICC increased to 2.2 % for Model 1 with the inclusion of the level-1 covariates and random slope for cigarettes/day. This was due to the reduction in the level-1 residual variance (560.5 to 554.7) relative to the increase in the level-2 random intercept variance (78.7 to 83.4). The (now adjusted) ICCadj is conditional for the individual composition of the DAs including the random slope for cigarettes/day held constant at 0 (i.e. non-smokers). The addition of DA-level variables in Model 2 removed a lot of the DA-level variance reducing the ICCadj to 0.6 %.

Table 4 Random Effects and Model Diagnostics

The level-2 random intercept variance term (reported as standard deviations in Table 4) indicates that the mean birth weight for every DA has a degree of variability from the overall (BC-wide) mean birth weight. For the Null model, the overall birth weight intercept is 3434.3 g with a standard deviation of 78.7 giving an 8.6 % difference in range between 95 % of the DAs (3434.3 ± (1.96 × 78.7) = 3280.0 and 3588.6 g). The quadratic form of the random slope for cigarettes/day in Model 1 prevents a similar calculation to be performed, but Fig. 4 gives an indication of the large between-DA slope variability which shows the DA-specific slopes of maternal smoking on birth weight. The intercept-slope correlation shown in Table 4 indicates the presence of DA-level heterogeneity signifying that DAs with higher average birth weights from non-smoking mothers have a lower within-DA effect of smoking (i.e. higher average DA intercepts tend to have lower average slopes for smoking) [23, 38].

Fig. 4
figure 4

Neighbourhood-specific slopes of maternal smoking on birth weight. Empirical Bayes predictions of DA-specific regression lines for Model 1

The level-1 and level-2 explained PCV (L1-PCV & L2-PCV) summarizes the relative degree of explained variance at the different levels between the different models (Table 4). Using the Null model as the reference, the Model 1 resulted in an L1-PCV of 2.0 %, and the L2-PCV in the random intercept was−12.3 %. The negative L2-PVC is a result of the larger level-2 intercept variance relative to the Null model. The addition of the DA-level variables in Model 2 explained 68.2 % of the DA-level variance compared to the Null model. Model 3 accounted for an additional 0.3 % of the L2-PCV.

Sensitivity analyses using only the non-imputed DAs (N1 = 229,067 in 6230 DAs) showed very minor differences in magnitude of significant variables. Most of the DAs that were missing data were in rural areas with small population numbers, the likely reason why their data were suppressed from the census tables. While some parameters were slightly attenuated, many of the interaction terms increased in magnitude. The random-intercept standard deviation was also slightly smaller than that of the same model using the full dataset while random-slope standard deviation showed no difference. In a second sensitivity analysis, we restricted the sample to only term births excluding stillbirths and congenital anomalies. As expected, there was a large reduction in the random slope variability (L2 slope (sd) = 5.7) and a small decrease in the random intercept variability (L2 intercept (sd) = 43.0) due to using only term births. Changes in the coefficients for the DA-level variables as well as their cross-level interactions with cigarettes/day were minor, within their 95 % confidence intervals reported in Table 3. The exception was the main effect of education which was no longer significant (p = 0.151), but its interaction with cigarettes/day did remain significant (p = 0.025). The maternal-level variables were attenuated but remained significant with the exception of the interaction between drug/alcohol flag and cigarettes/day which was no longer significant (p = 0.106).


This study employed multilevel random coefficient models to assess whether neighbourhood-level contextual factors can modify the effect of maternal smoking on birth weight. Our results show that the effect of maternal smoking on birth weight, self-reported as the number of cigarettes smoked per day, is modified by both individual-level and neighbourhood-level variables. However, the observed direction of the effect modification was not always as expected. Furthermore, through the use of random-slope models we show that the average effect of maternal smoking on birth weight can vary considerably between neighbourhoods which was only partially explained by the cross-level interactions. After adjusting for individual-level covariates and DA-level socio-economic, socio-demographic and air quality variables, there was a significant non-linear effect between cigarettes/day and birth weight in BC for singleton births from 2001 to 2006. This association was robust to the exclusion of stillbirths and congenital anomalies as well as the use of only term births demonstrating that selection bias does not likely affect the observed results.

The observed non-linear association between cigarettes/day and birth weight shown in Fig. 1 suggests that the largest potential effects are seen at the low to middle range of smoking levels. England et al. [40] found a very similar non-linear association of maternal smoking on term birth weight using self-reported cigarettes/day as well as using urine cotinine concentrations. Therefore, efforts to reduce the number of cigarettes smoked during pregnancy may have limited results for moderate and heavy smokers without substantial reductions or full cessation [41]. Interestingly, we found a similar curvilinear relationship with increasing levels of modeled PM2.5 and birth weight (Fig. 2c), a dose-response phenomenon observed in other exposure-disease contexts [42].

Beyond the non-linear association between cigarettes/day and birth weight, other factors were able to modify this relationship both positively and negatively. Our analysis confirm previously shown modification of the smoking-birth weight relationship by maternal risk factors [7, 8]; however, to our knowledge this is the first study to show that neighbourhood-level factors are able to modify this relationship. We found a significant negative interaction between cigarettes/day and neighbourhood-level SESi that resulted in the attenuation of the beneficial role of rising neighbourhood-level SES on birth weight with increased levels of maternal smoking (Fig. 2a). The predicted effects presented in Fig. 2a suggests is that maternal smoking may have little relevance in affecting birth weight in very low SES neighbourhoods, but becomes more prominent as neighbourhood-level SES increases and perhaps other stressors negatively impacting birth weight are reduced. Hence interventions focusing exclusively on individual behaviours may have limited success without addressing the contextual influences at the neighbourhood-level [9, 4345].

Conversely, the small but significant positive interaction between higher proportions of neighbourhood-level post-secondary education and cigarettes/day found that heavy smokers may benefit the most by living in higher educated neighbourhoods (Fig. 2b). This type of cross-level effect has been observed in other epidemiological scenarios where higher risk individuals have better outcomes than would be expected due to some beneficial capacity of the neighbourhood context [11, 22]. The mechanisms by which neighbourhood-level factors affect individual health is indirectly exerted through individual-level processes, such as behaviours, adaptations and attitudes which may be transmitted between people [46, 47]. Meng et al. found that low education neighbourhoods exert an impact on low birth weight and preterm birth through unhealthy behaviours, psycho-social stress (i.e. sense of control) and SES-related support [21]. Therefore it could be that smoking cessation rates in pregnancy are higher in better educated neighbourhoods where healthier behaviours are more common [48, 49]. Fig. 2b suggests that living in higher educated neighbourhoods may encourage moderate and heavy smokers to reduce their smoking frequency to less than five cigarettes/day.

Neighbourhood social supports and transmission of behaviours could also explain the observed interactions with higher immigrant density and rural address, albeit in opposite directions. The positive interaction between higher immigrant density and maternal smoking (Fig. 2d) may reflect the buffering effect of strong community cohesiveness and beneficial cultural practices [21, 43, 47]. Conversely, the observed negative interaction between rural address and cigarettes/day (Fig. 2e) could be due to the transmission of negative behaviours due to such behaviours being more common [50], and where less support for cessation may lead to smoking throughout pregnancy [51]. The dichotomized definition used to represent rural residential addresses may obscure mechanisms which can be modified by maternal factors such as education [52].

The buffering effect of PM2.5 with increased levels of maternal smoking (Fig. 2c) is curious but could provide evidence for a protective pre-conditioning stress that activates an adaptive response and increases biological resistance to cigarette-induced harms [53, 54]. We found a similar positive interaction between suspected alcohol and drug use and PM2.5 in a different analysis [55]. The suspicion of survival bias due to competing risks was partly mitigated by using a near full population sample that included stillbirths, congenital anomalies and preterm births, although we were not able to control for fetal loss prior to 20 weeks gestation. Other explanations require further scrutiny as evidence of the opposite (negative and synergistic) effect between smoking and air pollutants has been shown [10, 16].

We have shown in an earlier paper that women who reported smoking 10 or more cigarettes/day at their first prenatal visit were significantly more likely to have other maternal risk factors, such as lower education, suspected drug or alcohol use, and fewer prenatal care visits [4]. Our current results compliment the previous study by showing that the cumulative impact of multiple risk factors can have more than an additive effect on birth weight reduction. The negative association between older maternal age and birth weight was markedly greater with increased levels of maternal smoking, particularly among the heaviest smoking group (Fig. 3a). Similarly, those who reported higher levels of smoking who were also suspected of drug or alcohol use showed a pronounced effect compared to those who reported to not smoke (Fig. 3b). These results corroborate the established literature showing similar synergistic interactions between both maternal alcohol use and smoking on lower birth weights [8, 56], as well as between maternal smoking and older maternal age on birth weight [7, 57].

While the application of multilevel models in perinatal epidemiology have become more common [58], most have been random-intercept models with very few including a random-slope parameter. Permitting the slope for the maternal cigarettes/day exposure to be random provides information on how its effect on birth weight differs between neighbourhoods and enables the search for neighbourhood-level variables to help explain the between-neighbourhood variance [38]. For example, the random-slope standard deviation presented in Table 4 drops from 10.7 in Model 1 to 9.0 in the fully adjusted Model 3. This represents a 30 % change in explained random-slope variance (10.72–9.02/10.72). Furthermore, the addition of the level-2 variables explained 68.5 % of the random-intercept variance compared to the Null (empty) model. However in light of these findings, significant inter-DA variance remained for both the random intercept and slope.

This study used self-reports of smoking (cigarettes/day) recorded at the first prenatal visit; however, there were no data on exposure to environmental tobacco smoke or whether smoking reduction or cessation occurred during the pregnancy. The self-reporting bias of cigarette consumption can lead to the attenuation of the true effect of smoking on birth weight [59], and may therefore alter observed interactions. Studies of smoking misclassification in the United States has estimated non-disclosure to be around 20 % [60, 61]. The demographic predictors of non-disclosure include former smokers and younger maternal age which could partially explain the observed interaction between maternal age and cigarettes/day [61]. Similarly, recall bias and perceived stigma may result in under-reporting of actual consumption habits. This could account for the observed curvilinear effect on birth weight if women smoking 10 cigarettes/day report only smoking 5 per day, although England et al. observed a similar slope using urine-cotinine concentrations [40]. While relatively small, sensitivity analyses regarding the list-wise deletion of observations with missing smoking data was to exclude potentially at-risk pregnancies and could therefore alter coefficient estimates (n = 2501, 1.1 % of sample).

Another limitation includes potential measurement error and misclassification bias in the PM2.5 exposure assessment which could affect its estimates. First, the LUR PM2.5 concentrations may be underestimated with less variability compared to compiled monitoring data which could potentially underestimate its association with birth weight in certain areas [62]. Also the PM2.5 LUR model is cross-sectional based on 2006 air quality monitoring data, and we therefore assume that the study population was exposed to the same levels of PM2.5 across 6 year study period based on their residential DA. Finally, our analysis was based on maternal place of residence at delivery, and therefore intra-urban commuting and potential inter-urban relocation within the pregnancy period was not accounted for. Time-activity patterns show that pregnant women spend more time at home in the later stages of pregnancy, but mobility patterns may differ by age, parity and SES [63, 64].

A main strength of this study is the quality of the perinatal registry data [65]. The near 100 % ascertainment of birth records for the province of BC and quality control measures used in database management practices produces highly reliable data on maternal and newborn health outcomes, co-morbidities and exposures. However, the inability to control for individual-level SES, particularly maternal education, may influence the neighbourhood-level effect estimates and interactions. Maternal education is a variable provided in the PSBC Perinatal Registry, but was only available for 10 % of our population cohort. The adjustment for socially-patterned behavioural risk factors such as maternal smoking and suspected drug or alcohol use will control for some individual-level SES differences [4]. Notwithstanding, our results suggest that reported number of cigarettes smoked correlates with a substantial reduction in birth weight and is modified by socioeconomic, demographic and environmental risk factors suggesting the information as provided will help identify those at highest risk.


The effect of maternal smoking on birth weight is not constant across geography, but rather is context specific given the social and physical environment. The use of random coefficient models revealed neighbourhood-level differences in how maternal smoking negatively impacted birth weight demonstrating effect modification by neighbourhood and maternal-level factors. The inclusion of the DA-level SES, demographic and PM2.5 variables explained 68.5 % of the random intercept variability in DA-mean birth weight. However, the random slope variability was only partially explained by the cross-level interactions suggesting other contextual factors are involved in determining the magnitude of maternal smoking on birth weight. Further studies are needed to corroborate our findings and to understand how neighbourhood-level attributes interact with smoking to affect birth outcomes.


  1. U.S. Department of Health and Human Services. Reproductive and Developmental Effects. How Tobacco Smoke Causes Disease: The Biology and Behavioural Basis for Smoking Attributable Disease. Atlanta: A Report of the Surgeon General; 2010. p. 521–701.

    Google Scholar 

  2. Burstyn I, Kuhle S, Allen AC, Veugelers P. The role of maternal smoking in effect of fetal growth restriction on poor scholastic achievement in elementary school. Int J Environ Res Public Health. 2012;9(2):408–20. Molecular Diversity Preservation International.

    Article  PubMed  PubMed Central  Google Scholar 

  3. Behl M, Rao D, Aagaard K, Davidson TL, Levin ED, Slotkin TA, et al. Evaluation of the association between maternal smoking, childhood obesity, and metabolic disorders: a national toxicology program workshop review. Environ Health Perspect. 2013;121(2):170–81.

    PubMed  Google Scholar 

  4. Erickson AC, Arbour LT. Heavy smoking during pregnancy as a marker for other risk factors of adverse birth outcomes: a population-based study in British Columbia Canada. BMC Public Health. 2012;12(1):102.

    Article  PubMed  PubMed Central  Google Scholar 

  5. Mohsin M, Bauman AE. Socio-demographic factors associated with smoking and smoking cessation among 426,344 pregnant women in New South Wales Australia. BMC Public Health. 2005;5:138.

    Article  PubMed  PubMed Central  Google Scholar 

  6. Al-Sahab B, Saqib M, Hauser G, Tamim H. Prevalence of smoking during pregnancy and associated risk factors among Canadian women: a national survey. BMC Pregnancy Childbirth. 2010;10:24.

    Article  PubMed  PubMed Central  Google Scholar 

  7. Cnattingius S. Does age potentiate the smoking-related risk of fetal growth retardation? Early Hum Dev. 1989;20(3–4):203–11.

    Article  CAS  PubMed  Google Scholar 

  8. Aliyu MH, Wilson RE, Zoorob R, Brown K, Alio AP, Clayton H, et al. Prenatal alcohol consumption and fetal growth restriction: potentiation effect by concomitant smoking. Nicotine Tob Res. 2009;11(1):36–43.

    Article  PubMed  Google Scholar 

  9. Blumenshine P, Egerter S, Barclay CJ, Cubbin C, Braveman PA. Socioeconomic disparities in adverse birth outcomes: a systematic review. Am J Prev Med. 2010;39(3):263–72.

    Article  PubMed  Google Scholar 

  10. Currie J, Neidell M, Schmieder JF. Air pollution and infant health: lessons from New Jersey. J Health Econ. 2009;28:688–703.

    Article  PubMed  PubMed Central  Google Scholar 

  11. Cubbin C, Sundquist K, Ahlén H, Johansson S-E, Winkleby MA, Sundquist J. Neighborhood deprivation and cardiovascular disease risk factors: protective and harmful effects. Scand J Public Health. 2006;34(3):228–37.

    Article  PubMed  Google Scholar 

  12. Stieb DM, Chen L, Eshoul M, Judek S. Ambient air pollution, birth weight and preterm birth: a systematic review and meta-analysis. Environ Res. 2012;117:100–11.

    Article  CAS  PubMed  Google Scholar 

  13. HEI Review Panel on Ultrafine Particles. Understanding the Health Effects of Ambient Ultrafine Particles. HEI Perspectives 3. Boston, MA: Health Effects Institute; 2013.

  14. Erickson AC, Arbour L. The shared pathoetiological effects of particulate air pollution and the social environment on fetal-placental development. J Environ Public Health. 2014;2014:901017.

    Article  PubMed  PubMed Central  Google Scholar 

  15. Al-Gubory KH. Environmental pollutants and lifestyle factors induce oxidative stress and poor prenatal development. Reprod Biomed Online. 2014;29(1):17–31.

    Article  CAS  PubMed  Google Scholar 

  16. Valavanidis A, Vlachogianni T, Fiotakis K. Tobacco smoke: involvement of reactive oxygen species and stable free radicals in mechanisms of oxidative damage, carcinogenesis and synergistic effects with other respirable particles. Int J Environ Res Public Health. 2009;6(2):445–62.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  17. Jauniaux E, Burton GJ. Morphological and biological effects of maternal exposure to tobacco smoke on the feto-placental unit. Early Hum Dev. 2007;83(11):699–706.

    Article  CAS  PubMed  Google Scholar 

  18. Joseph KS, Liston RM, Dodds L, Dahlgren L, Allen AC. Socioeconomic status and perinatal outcomes in a setting with universal access to essential health care services. CMAJ. 2007;177(6):583–90.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  19. Mortensen LH, Lauridsen JT, Diderichsen F, Kaplan GA, Gissler M, Andersen A-MN. Income-related and educational inequality in small-for-gestational age and preterm birth in Denmark and Finland 1987–2003. Scand J Public Health. 2010;38(1):40–5.

    Article  PubMed  Google Scholar 

  20. Kramer MS, Seguin L, Lydon J, Goulet L. Socio-economic disparities in pregnancy outcome: why do the poor fare so poorly? Paediatr Perinat Epidemiol. 2000;14(3):194–210.

    Article  CAS  PubMed  Google Scholar 

  21. Meng G, Thompson ME, Hall GB. Pathways of neighbourhood-level socio-economic determinants of adverse birth outcomes. Int J Health Geogr. 2013;12(1):32.

    Article  PubMed  PubMed Central  Google Scholar 

  22. Cubbin C, Winkleby MA. Protective and harmful effects of neighborhood-level deprivation on individual-level health knowledge, behavior changes, and risk of coronary heart disease. Am J Epidemiol. 2005;162(6):559–68.

    Article  PubMed  Google Scholar 

  23. Merlo J, Yang M, Chaix B, Lynch J, Råstam L. A brief conceptual tutorial on multilevel analysis in social epidemiology: investigating contextual phenomena in different groups of people. J Epidemiol Community Health. 2005;59(9):729–36.

    Article  PubMed  PubMed Central  Google Scholar 

  24. Auger N, Giraud J, Daniel M. The joint influence of area income, income inequality, and immigrant density on adverse birth outcomes: a population-based study. BMC Public Health. 2009;9:237.

    Article  PubMed  PubMed Central  Google Scholar 

  25. Kramer MS. Determinants of low birth weight: methodological assessment and meta-analysis. Bull World Health Organ. 1987;65(5):663–737.

    CAS  PubMed  PubMed Central  Google Scholar 

  26. Perinatal Services BC. British Columbia Perinatal Data Registry [Internet]. Vancouver; [cited 2015 Dec 7]. Available from: Accessed 7 Dec 2015.

  27. DMTI Spatial. GeoPinpoint 6.4. Markham, Ont; 2009.

  28. ESRI. ArcGIS 10.2. Redlands: Environmental Systems Research Institute; 2010.

    Google Scholar 

  29. Gauvin L, Robitaille E, Riva M, McLaren L, Dassa C, Potvin L. Conceptualizing and operationalizing neighbourhoods: the conundrum of identifying territorial units. Can J public Heal. 2007;98 Suppl 1:S18–26.

    Google Scholar 

  30. Holowaty EJ, Norwood T a, Wanigaratne S, Abellan JJ, Beale L. Feasibility and utility of mapping disease risk at the neighbourhood level within a Canadian public health unit: an ecological study. Int J Health Geogr. 2010;9:21.

    Article  PubMed  PubMed Central  Google Scholar 

  31. du Plessis V, Beshiri R, Bollman RD, Clemenson H. Definitions of “Rural”. Agriculture and Rural Working Paper Series, No.: 61. Ottawa ON: Statistics Canada, Agriculture Division; 2002.

  32. Hystad P, Setton E, Cervantes A, Poplawski K, Deschenes S, Brauer M, et al. Creating national air pollution models for population exposure assessment in Canada. Env Heal Perspect. 2011;119(8):1123–9.

    Article  CAS  Google Scholar 

  33. Chan E, Serrano J, Chen L, Stieb DM, Jerrett M, Osornio-Vargas A. Development of a Canadian socioeconomic status index for the study of health outcomes related to environmental pollution. BMC Public Health. 2015;15(1):714.

    Article  PubMed  PubMed Central  Google Scholar 

  34. Shiono PH, Klebanoff MA, Graubard BI, Berendes HW, Rhoads GG. Birth weight among women of different ethnic groups. JAMA. 1986;255(1):48–52.

    Article  CAS  PubMed  Google Scholar 

  35. Janssen P, Thiessen P, Klein M, Whitfield MF, Macnab YC, Cullus-Kuhl SC. Standards for the measurement of birth weight, length and head circumfrence at term in neonates of European, Chinese and south Asian ancestry. Open Med. 2007;1(2):E74–88.

    PubMed  PubMed Central  Google Scholar 

  36. Graham MH. Confronting multicollinearity in ecological multiple regression. Ecology. 2003;84(11):2809–15.

    Article  Google Scholar 

  37. Statistics Canada. Census of Canada. Topic-based Tabulations, 2006. Ottawa: Statistics Canada; 2006.

  38. Snijders TAB, Bosker RJ. Multilevel analysis: an introduction to basic and advanced multilevel modeling. 2nd ed. London: Sage; 2012.

    Google Scholar 

  39. StataCorp LP. Stata 13. College Station: StataCorp LP; 2014.

    Google Scholar 

  40. England LJ, Kendrick JS, Gargiullo PM, Zahniser C, Hannon WH. Measures of maternal tobacco exposure and infant birth weight at term. Am J Epidemiol. 2001;153(10):954–60.

    Article  CAS  PubMed  Google Scholar 

  41. Secker-Walker RH, Vacek PM, Flynn BS, Mead PB. Estimated gains in birth weight associated with reductions in smoking during pregnancy. J Reprod Med. 1998;43(11):967–74.

    CAS  PubMed  Google Scholar 

  42. Wigle DT, Lanphear BP. Human health risks from low-level environmental exposures: no apparent safety thresholds. PLoS Med. 2005;2(12), e350.

    Article  PubMed  PubMed Central  Google Scholar 

  43. Buka S, Brennan RT, Rich-Edwards JW, Raudenbush SW, Earls F. Neighborhood support and the birthweight of urban infants. Am J Epidemiol. 2002;157(1):1–8.

    Article  Google Scholar 

  44. McCormick MC, Brooks-Gunn J, Shorter T, Holmes JH, Wallace CY, Heagarty MC. Factors associated with smoking in low-income pregnant women: relationship to birth weight, stressful life events, social support, health behaviors and mental distress. J Clin Epidemiol. 1990;43(5):441–8.

    Article  CAS  PubMed  Google Scholar 

  45. Lumley J, Chamberlain C, Dowswell T, Oliver S, Oakley L, Watson L. Interventions for promoting smoking cessation during pregnancy. Cochrane Database Syst Rev. 2009;(3). Art.No.:CD001055. doi:10.1002/14651858.

  46. Diez Roux AV. Investigating neighborhood and area effects on health. Am J Public Health. 2001;91(11):1783–9.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  47. Mitchell RE, Trickett EJ. Task force report: social networks as mediators of social support. Community Ment Health J. 1980;16(1):27–44.

    Article  CAS  PubMed  Google Scholar 

  48. Connor SK, McIntyre L. The sociodemographic predictors of smoking cessation among pregnant women in Canada. Can J public Heal. 1999;90(5):352–5.

    CAS  Google Scholar 

  49. Lu Y, Tong S, Oldenburg B. Determinants of smoking and cessation during and after pregnancy. Health Promot Int. 2001;16(4):355–67.

    Article  CAS  PubMed  Google Scholar 

  50. Lisonkova S, Sheps SB, Janssen PA, Lee SK, Dahlgren L, Macnab YC. Birth outcomes among older mothers in rural versus urban areas: a residence-based approach. J Rural Health. 2011;27(2):211–9.

    Article  PubMed  Google Scholar 

  51. Borland T, Babayan A, Irfan S, Schwartz R. Exploring the adequacy of smoking cessation support for pregnant and postpartum women. BMC Public Health. 2013;13:472.

    Article  PubMed  PubMed Central  Google Scholar 

  52. Auger N, Authier M-A, Martinez J, Daniel M. The association between rural-urban continuum, maternal education and adverse birth outcomes in Québec, Canada. J Rural Heal. 2009;25(4):342–51.

    Article  Google Scholar 

  53. Sanders CL, Scott BR. Smoking and hormesis as confounding factors in radiation pulmonary carcinogenesis. Dose-response. 2008;6(1):53–79.

    Article  CAS  Google Scholar 

  54. Calabrese EJ, Bachmann K a, Bailer a J, Bolger PM, Borak J, Cai L, et al. Biological stress response terminology: integrating the concepts of adaptive response and preconditioning stress within a hormetic dose-response framework. Toxicol Appl Pharmacol. 2007;222(1):122–8.

    Article  CAS  PubMed  Google Scholar 

  55. Erickson AC, Ostry A, Chan LHM, Arbour L. The reduction of birth weight by fine particulate matter and its modification by maternal and neighbourhood-level factors: a multilevel analysis in British Columbia, Canada. Environmental Health; 2016;15:51.

  56. Windham GC, Fenster L, Hopkins B, Swan SH. The association of moderate maternal and paternal alcohol consumption with birthweight and gestational age. Epidemiology. 1995;6(6):591–7.

    Article  CAS  PubMed  Google Scholar 

  57. Fox SH, Koepsell TD, Daling JR. Birth weight and smoking during pregnancy—effect modification by maternal age. Am J Epidemiol. 1994;139(10):1008–15.

    CAS  PubMed  Google Scholar 

  58. Metcalfe A, Lail P, Ghali WA, Sauve RS. The association between neighbourhoods and adverse birth outcomes: a systematic review and meta-analysis of multi-level studies. Paediatr Perinat Epidemiol. 2011;25(3):236–45.

    Article  PubMed  Google Scholar 

  59. England LJ, Grauman A, Qian C, Wilkins DG, Schisterman EF, Yu KF, et al. Misclassification of maternal smoking status and its effects on an epidemiologic study of pregnancy outcomes. Nicotine Tob Res. 2007;9(10):1005–13.

    Article  PubMed  Google Scholar 

  60. Dietz PM, Homa D, England LJ, Burley K, Tong VT, Dube SR, et al. Estimates of nondisclosure of cigarette smoking among pregnant and nonpregnant women of reproductive Age in the united states. Am J Epidemiol. 2011;173(3):355–9.

    Article  PubMed  Google Scholar 

  61. Shipton D, Tappin DM, Vadiveloo T, Crossley JA, Aitken DA, Chalmers J. Reliability of self reported smoking status by pregnant women for estimating smoking prevalence: a retrospective, cross sectional study. BMJ. 2009;339(1):b4347.

    Article  PubMed  PubMed Central  Google Scholar 

  62. Canadian Council of Ministers of the Environment. Canada-wide standards for particulate matter and ozone: five year report: 2000–2005. Report No.: PN 1374. Winnipeg; 2006.

  63. Nethery E, Brauer M, Janssen P. Time-activity patterns of pregnant women and changes during the course of pregnancy. J Expo Sci Env Epidemiol. 2009;19(3):317–24.

    Article  Google Scholar 

  64. Bell ML, Belanger K. Review of research on residential mobility during pregnancy: consequences for assessment of prenatal environmental exposures. J Expo Sci Environ Epidemiol. 2012;22(5):429–38.

    Article  PubMed  PubMed Central  Google Scholar 

  65. Frosst G, Hutcheon J, Joseph KS, Kinniburgh B, Johnson C, Lee L. Validating the British Columbia perinatal data registry: a chart re-abstraction study. BMC Pregnancy Childbirth. 2015;15(1):123.

    Article  PubMed  PubMed Central  Google Scholar 

Download references


We would like to acknowledge Dr. Scott Venners for his contributions in obtaining the CIHR Operational grant. We acknowledge the staff at Perinatal Services BC for their support regarding data access and review of manuscript. All inferences, opinions, and conclusions drawn in this publication are those of the authors and do not reflect the opinions or policies of Perinatal Services BC.


This publication was funded in part by the Canadian Institute of Health Research (CIHR) Operational Grant (protocol #: 200903-202069). Funders provided no role in the design, conduct, data collection, management, analysis, interpretation, preparation, review and approval of this study.

Availability of data and materials

The Data Sharing Agreement signed with Perinatal Services BC (PSBC) does not allow for the data to be shared; however, permissions can be granted under a formal request to and approval from PSBC.

Authors’ contributions

AE was primarily responsible for the study design, acquisition of data, analysis and interpretation of results, drafting the article and making revisions. LA critically assessed the concept, study design, analysis and interpretation of results. LA provided critical edits of the content. Final approval of submitted manuscript was provided by both authors. LC and AO provided useful edits on draft manuscripts and input into the manuscript framework. All authors give final approval of the version to be published.

Competing interests

The authors declare they have no competing financial interests.

Ethics approval and consent to participate

Research ethics board approval was granted by the University of Victoria (protocol #11-043).

Author information

Authors and Affiliations


Corresponding author

Correspondence to Laura Arbour.

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (, which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Erickson, A.C., Ostry, A., Chan, H.M. et al. Air pollution, neighbourhood and maternal-level factors modify the effect of smoking on birth weight: a multilevel analysis in British Columbia, Canada. BMC Public Health 16, 585 (2016).

Download citation

  • Received:

  • Revised:

  • Accepted:

  • Published:

  • DOI: