Identifying exposure pathways mediating adverse birth outcomes near active surface mines in Central Appalachia

doi:10.21203/rs.3.rs-948330/v2

PDF

Research

Identifying exposure pathways mediating adverse birth outcomes near active surface mines in Central Appalachia

https://doi.org/10.21203/rs.3.rs-948330/v2

This work is licensed under a CC BY 4.0 License

Version 2

posted 17 Apr, 2022

You are reading this latest preprint version

Background: Previous work has determined an association between proximity to active surface mining within Central Appalachia and an increased risk of preterm birth (PTB) and low birthweight (LBW). Multiple potential exposure pathways may exist, however, including inhalation of particulate matter (airshed exposure), or exposure to impacted surface waters (watershed exposure). We hypothesize that this relationship is mediated by exposure to contaminants along one or both of these pathways.

Methods: We geolocated 194,084 birth records through health departments in WV, KY, VA and TN between 1990 and 2015. We performed a mediation analysis, iteratively including within our models: a) the percent of active surface mining within 5 km of maternal residence during gestation; b) the cumulative surface mining airshed trajectories experienced during gestation, and; c) the percent of active surface mining occurring within the watershed of residency during gestation.

Results: Our baseline models found that active surface mining was associated with an increased odds of PTB (1.09, 1.05 – 1.13) and LBW (1.06, 1.02 – 1.11), controlling for individual-level predictors. When mediators were added to the baseline model, the association between active mining and birth outcomes became non-significant (PTB: 0.48, 0.14 – 1.58; LBW 0.78, 0.19 – 3.00), while the association between PTB and LBW remained significant by airshed exposure (PTB: 1.14, 1.11 – 1.18; LBW: 1.06, 1.03 – 1.10).

Conclusions: Our results found that surface mining airsheds at least partially explained the association between active mining and adverse birth outcomes, consistent with a hypothesis of mediation, while mediation via the watershed pathway was less evident.

Environmental Policy

Mediation

surface mining

maternal health

environmental health

pre-term birth

low birthweight

coal mining

What this study adds: Previous work published in Environmental Epidemiology has determined an association between living near surface mining activities within Central Appalachia and an increased risk of adverse birth outcomes. The exposure pathways underlying this association, however, remain unclear, including inhalation of particulate matter versus exposure to contaminated water. This work aims to identify the potential exposure pathways underpinning this relationship, which may be mediated by one or both pathways. Our results suggest that the association between surface mining and adverse birth outcomes is driven in part by exposure to airsheds near active surface mining activities, while mediation via watersheds was less clear.

Communities within the Central Appalachia region, containing areas within Virginia, West Virginia, Kentucky and Tennessee, experience some of the most severe economic and health disparities across the United States.¹ While poverty has generally decreased since the mid-1900’s, the region has not experienced corresponding expected health gains,² and continues to lag behind national poverty estimates, with increased rates of unemployment, lower levels of post-secondary education, and lower life expectancy.³ Surface mining became more widespread in the 1990s, and represented over 60% of coal mined within 2019.⁴ Within Central Appalachia, coal production has declined in the 2010s and there were nearly double the amount of surface mines as underground mines in 2019, despite relatively similar levels of coal production.⁴ Surface mining has disrupted the local environment, where forest is cleared, rocks and topsoil are loosened through explosives to expose the underlying coal seams, and the resultant rubble is dumped into nearby valleys and headstreams.^5,6

These surface mining activities have been previously associated with adverse health outcomes.^1,7,8 For example, a county-level cross-sectional analysis found higher rates of respiratory disease, hypertension and kidney disease among populations living in counties with greater coal production, as compared to those living further away from coal mining activities.⁹ However, the relative importance of the exposure pathways underlying these associations remain less clear. Previous research outlining the public health impacts of surface mining activities has proposed two possible biological mechanisms to explain these relationships.⁸ Firstly, there is evidence that inhalation of particulate matter is increased near surface mining activities.^10,11 Within Central Appalachia, surface mining activities result in increased emissions of acidic aerosols such as sulfur dioxide and nitrous dioxide,¹² as well as particulate matter smaller than 10 microns in diameter (PM₁₀) as a result of trucks hauling coal.¹³ Studies have found that particulate matter collected from nearby surface mines results in increased inflammatory responses and vascular tissue damage in a rodent model ^10,11 and increased carcinogenic markers in human lung cells,¹⁴ which is consistent with the higher rates of cardiovascular and lung disease observed in the region.¹¹ Exposure to air pollution is known to result in a range of adverse health outcomes across life stages, including adverse birth outcomes.^15,16

Secondly, exposure to heavy metals in drinking water contaminated by surface mining activities has also been proposed as a causal pathway linking adverse health outcomes and proximity to surface mining activities.⁸ However, little research is available that quantifies the in situ exposure to heavy metals, water quality and sources of drinking water used within Appalachia.¹ While associations between metals in water and increased rates of cancer and chronic disease have been noted, researchers have called for a more thorough investigation of these associations linked specifically to exposure to waterborne contaminants and dose of exposure, while controlling for behavioral risk factors such as tobacco use.⁸ Despite this lack of evidence, other studies within the field of aquatic ecology have noted elevated concentrations of heavy metals and nonmetallic ions downstream of active and historic surface mining activities, resulting in loss of aquatic biodiversity within Appalachian headwater streams.⁶ Given that these streams are inevitably hydrologically connected to source waters for regional drinking water supplies, degradation of aquatic ecological health potentially indicates important concurrent public health risks.

Surface mining activities have more specifically been shown to correlate with adverse birth outcomes, such as birth defects, low birthweight, and premature birth. Among nearly 2 million births across four Central Appalachian states, Ahern et al. suggested that rates of birth defects were higher in areas with mountaintop removal as compared to areas with no mountaintop removal, after controlling for covariates.¹⁷ Other studies have further suggested that air pollution containing fine particulate matter, such as particulates resulting from surface mining activities, can reach the fetal placenta and trigger inflammatory responses,¹⁸ which may in turn explain associations between air pollution and adverse birth outcomes such as low birthweight and pre-term birth.^19,20 Within the Central Appalachia region, Ahern et al. have documented an independent risk of low birthweight among women living within an area with coal mining activities in West Virginia.²¹ More recently, Buttling et al. suggested an increased risk of low birthweight and pre-term birth among women living within 5km of an active surface mine.²²

Here, this work builds upon the findings outlined in Buttling et al., to examine the potential underlying pathways mediating the documented association between active surface mining activities and adverse birth outcomes. More specifically, we aim to explore airshed and watershed exposure pathways as potential mediators of the association observed between proximity to active surface mining activities and adverse birth outcomes.²²We undertake a mediation analysis within a logistic regression modelling framework to examine the odds of preterm birth, low birthweight, and term low birthweight outcomes as predicted by maternal proximity to surface mining activities during gestation from 1990 – 2015, overlaying surface mine airshed and watershed boundaries. To undertake these analyses, we use geolocated birth records spanning four states in the Central Appalachia region, representing nearly 200,000 births over 2.5 decades. Specifically, we hypothesize that increased odds of preterm birth and low birthweight associated with proximity to surface mining activities during gestation may be mediated by airshed or watershed pathways, or both, and examine the potential additive effects of these mediators.

Data

We obtained birth records for four states across Central Appalachia, provided by departments of health in Kentucky (KY), Tennessee (TN), Virginia (VA), and West Virginia (WV) (Figure 1). This dataset is described in detail in Buttling et al.;²² briefly, a total of 409,394 birth records were obtained from departments of health with street address of reported maternal address. Records were removed due to missingness in street addresses with many records reporting only mailing addresses (such as a rural route or P.O. box), resulting in a final dataset comprising 194,084 births between 1990 and 2015 with a geocoded maternal residential address.. Our previous analysis shows differences in maternal characteristics of births from the original and final dataset are minimal and an additional analysis using a zip-code level exposure variable that allowed inclusion of most of the original records resulted in similar effect estimates.²² The exposure variable of interest within this dataset was the proportion of land designated as ‘active surface mining’ within a 5km radius of the maternal address during the majority gestation year. Pre-mining areas were defined as areas that were classified as actively mined in future years of the study period. Mined areas <40 acres in size were removed from the analysis, as the Office of Surface Mining Reclamation and Enforcement reported that economically viable mines are generally at least 40 acres in size.²³

More specifically, yearly active surface mining activities were defined using a combination of land cover change and mining permits, obtained from Marston and Kolivras for the years 1986 through 2015 using 30 m resolution Landsat remote sensing imagery.²⁴ Specifically, this dataset used remote sensing technology and satellite imagery to quantify changes in the extent of surface mining areas greater than 40-acres over the corresponding years among areas where mining permits were requested. For each year, barren land on which vegetation had been removed was identified, and pixels that were tied to other types of land disturbance (such as clear cutting of timber) were manually excluded. Remaining barren pixels lying within areas permitted for mining within the USGS-defined Appalachian coalfields were designated as places where active mining was likely taking place. Active surface mining was delineated from post-mining areas through classification of vegetation, where actively mined areas tend to have major land disruption and degradation, while post-mining areas tend to show some revegetation through reclamation efforts.²⁵

In order to assess the accuracy of the mined area dataset, Marston et al. applied a standard validation process within remote sensing to the identified mined layer. Specifically, 2,250 points were placed within eight randomly selected counties across the Appalachian coalfield region of Central Appalachia. Data from these counties had not previously been used for the initial identification (training) of mined areas. The classification of these points (mined or not) was compared to aerial imagery from the National Agricultural Imagery Program (NAIP, 1 meter resolution), land cover data from the National Land Cover Dataset, and Landsat images that were not used for initial classification. The overall classification accuracy was 0.88, or 88%, and the kappa coefficient to measure agreement was also 0.88, indicating “strong” agreement.²⁶ When errors were identified, they tended to be along the edges of mined areas. If a pixel at the edge of a cleared area is half barren and half forested, the classification process will classify that pixel as one or the other, with the potential for error to occur.

Figure 1

Airsheds for each of the active surface mines in each year of analysis were estimated using the HYbrid Single-Particle Lagrangian Integrated Trajectory (HYSPLIT4) atmospheric trajectory model,²⁷ which uses meteorological data to compute atmospheric trajectories, particle dispersion, and air concentrations. Detailed methods using HYSPLIT4 to characterize airsheds of surface mines are described in further detail in McKnight et al.²⁸ Briefly, we modeled individual airsheds of each active surface mine in the study area. The results of this process were raster data sets (e.g., a pixelated or gridded surface where each pixel or grid square corresponds to a specific geographical location and has an associated value) comprised of frequency values relating to air movement from surface mines per gridded cell. To quantify total exposure in pregnancy from all mines, we summed the frequency values of all airsheds extracted at each maternal address location for the cumulative amount of air from surface mines experienced at the maternal residential address. Watersheds were classified using 10-digit hydrologic unit codes (HUC10) within the United States Geological Survey’s Watershed Boundary Dataset, which represent the areal extent of surface water drainage using an aggregate collection of hydrologic unit data and amount of active surface mining within watersheds for each year of the study period was calculated.²⁹

Analysis

We employed mediation analyses to explore the potential exposure pathways that could explain associations between living in close proximity to active surface mining activities and increased odds of preterm birth and low birthweight described in Buttling et al., accounting for pre-mining differences and additional individual-level covariates available on birth records.²² Mediation analyses are helpful in exploring the underlying mechanisms underpinning a known relationship between an exposure and outcome.³⁰ Generally, mediation is suggested when four criteria are met, as outlined in Table 1.³¹

Table 1

Our independent variable of interest was the percent of land within a 5km radius of the maternal residential address that was designated as an active surface mine during the majority year of gestation (n = 23,733). Participant characteristics by exposure status are reported in Buttling et al. 2021, and are included in Supplementary Information (Table S5). Our dependent variables of interest included 1) preterm birth (PTB), defined as birth before 37 weeks of gestation; 2) low birthweight (LBW), defined as birthweight less than 2,500 grams; and 3) term low birthweight (TLBW), defined as birth occurring at ≥ 37 weeks gestation and birthweight less than 2,500 grams. Using these criteria, we tested whether mediation between proximity to active surface mining activities during the year containing the majority of the pregnancy (majority gestation year) and these adverse birth outcomes occurred via the airshed pathway or the watershed pathway. We quantified these potential mediators (Figure 2), respectively, as: 1) the cumulative potential exposure to air pollutants via the airshed experienced at the maternal residential address during the majority gestation year, and 2) the percent of land experiencing active surface mining within the watershed of residency during the majority gestation year. The cumulative potential exposure to airborne pollutants via the airshed is outlined in further detail within Mcknight et al.;²⁸ briefly, these values represent the cumulative frequency of air originating from active surface mines, as modeled via the HYSPLIT4 atmospheric trajectory model. Residential addresses associated with higher values have higher potential exposure to surface mining air pollution, including fine particulate matter.

Figure 2

We firstly employed a Sobel test using the bda package in R software to quantify whether mediators significantly influenced the relationship between the independent and dependent variables.^32,33 We further performed a logistic regression analysis using the stats package within the base R software. Covariates within our model included categorical variables found on birth records including for maternal age(‘18 - 35 years’,‘<18 years’,‘> 35 years’), highest education attained by the mother at the child’s birth (‘8th grade or less’, ‘9th - 12th grade (includes high school graduates)’, ‘Post high-school education (with or without degree)’), race (‘White’, ’Black’, ’Other’), ethnicity (‘Hispanic’, ‘Not Hispanic’), self-reported tobacco use during pregnancy (‘Yes’,’No’), sex of the child (‘Male’, ‘Female’), payment type for birth medical services (‘Medicaid’, ‘Private Insurance’, ‘Self-Pay’, ‘Other’), state (‘Kentucky’, ‘Tennessee’, ‘Virginia’, and ‘West Virginia’), and continuous percent of land within 5 km of maternal residence that was not actively experiencing mining activities during the majority gestation year, but would be subsequently mined in later years (referred to as ‘pre-mining’ activities), to account for any temporal baseline difference prior to active mining. Because mining activities tend to show spatial autocorrelation (e.g., active mining tends to move progressively across the landscape and is therefore closely correlated with pre-mining landcover), we further included an interaction term between the amount of pre-mined land and surface mining land within a 5 km buffer of maternal residence. Lastly, to allow for non-linear temporal trends observed within the data and account for serial autocorrelation, we included a spline with 4 degrees of freedom in the year covariate using the splines package within R software. This approach is a mixed effects modeling approach, similar to including random effects in the model, where fixed effects are allowed to vary at inflection (or spline) points, estimated by maximizing a penalized likelihood function.³⁴ This model builds on the model detailed in Buttling et al.,²² which can generally be specified as:

where logit(y_i) represents the odds of an adverse outcome (PTB, LBW, and TLBW) occurring for a given woman’s birth, i; M_irepresents the % of pre-mining land area within a 5km buffer of the maternal residence of individual i; T_irepresents the % of active surface mining area within a 5km buffer of the maternal residence of individual i; represents the interaction effect of having land within 5 km of the maternal residence of individual, i, that is both likely actively being mined (active mining) and will subsequently be mined (pre-mining); bs(t) represents a non-linear spline with 4 degrees of freedom to fit temporal trends in the data; and represents the suite of fixed covariates K as outlined above.²²

To explore hypothesized mediation pathways, we tested three model types which iteratively built on each other (Figure 3).³³ Our base adjusted model (Model 0) included the suite of fixed effects described above, including the independent variable of interest (e.g., % surface mining activities in a 5 km buffer). Model 1 included the same suite of variables as Model 0, with the addition of the airshed mediator variable alone, while Model 2 included the watershed mediator alone. Lastly, Model 3 included the addition of both the airshed and watershed mediators. All models included an interaction term between pre-mining surface area and active surface area as described above, plus an interaction term between the mediator(s) and exposure variable. We report model performance metrics, including McFadden’s pseudo-R² statistic, Akaike Information Criterion (AIC), and mean adjusted error (MAE) and root-mean squared error (RMSE) measures. AIC represents a measure of model predictive power as a trade-off with model complexity, with lower values generally representing models with better fit.³⁵ MAE and RMSE represent model precision and bias, with values closer to zero representing better model fit, while pseudo-R² values represent variance explained by the model, with values between 0.2 – 0.4 representing reasonable model fit across a range of applications.^36–39

Figure 3

Sobel test

Results of the Sobel test suggesting evidence of whether mediation is occurring are presented in Table 2. For each dependent variable of interest, we conducted an independent Sobel test exploring whether mediation is occurring via the airshed pathway (M = air) as well as the watershed pathway (M = water), with the independent variable of interest defined as proximity to surface mining. Evidence suggested that mediation was occurring via the airshed pathway for preterm birth and low birthweight outcomes (p < 0.001), while evidence was less clear for mediation of term low birthweight. There was further evidence suggesting that mediation may be occurring via the watershed pathway for preterm birth (p = 0.012) but mediation was less clear for other birth outcomes.

Table 2

Unadjusted model results

We firstly explored unadjusted model results for each dependent variable of interest, as outlined in Table 3. Specifically, the results of these unadjusted models are used to test whether the independent variable was significantly associated with the mediator, a key criterion indicating mediation. Towards this, we found that proximity to surface mining strongly influenced both hypothesized mediator pathways (p < 0.001). We further conducted unadjusted models on the dependent variables of interest before conducting adjusted models. Across all outcomes, both mediation pathways were generally strongly significant, with the exception of term low birthweight. Of note, results tended to be more strongly significant among the airshed mediation pathways, as opposed to the watershed mediation pathway.

Table 3

Adjusted model results

To test the remaining criteria (Table 1) for whether mediation occurred between adverse birth outcomes and proximity to surface mining, we iteratively conducted a logistic regression with the addition of mediating variables. Full model results are listed in supplemental tables S1 – S3. Overall, in our baseline models (e.g., without the addition of mediation pathways), we found that sociodemographic characteristics such as older age (> 35 years), self-reported tobacco use during pregnancy, having a female child, race, and type of payment used to cover birth services (e.g., Medicaid and ‘Other’) were associated with increased odds of preterm birth, low birthweight, and term low birthweight outcomes. We further found significant differences in the odds of preterm birth across states, with women giving birth in Central Appalachian counties within Tennessee at higher odds of having a preterm birth and women giving birth in Virginia at lower odds, as compared to women giving birth in Kentucky. Across all outcomes, self-reported tobacco use during pregnancy resulted in the highest odds ratios of adverse birth outcomes as compared to other sociodemographic characteristics within the model, particularly for term low birthweight (OR: 3.06, 95% CI: 2.80 – 3.35).

To explore mediation pathways, we firstly examined whether the independent variable was significantly associated with the dependent variables in absence of the mediators (mediation criterion #2), denoted as Model 0. We found proximity to surface mining was associated with preterm birth, such that a 1% increase in the amount of land within a 5 km buffer around the maternal address being mined was associated with a 9% increase in the odds of having preterm birth (OR: 1.09; 1.05 – 1.13) and 6% increase in the odds of low birthweight (OR: 1.06, 1.02 – 1.11). In practical terms, the highest percentage of surface mining within a 5 km buffer surrounding a maternal address was observed to be just over 15%, with a mean of under 1% (0.13%).

To test whether mediators significantly and uniquely associated with the dependent variable (mediation criterion #3), we added the airshed pathway variable (Model 1), watershed pathway variable (Model 2), and both (Model 3). Results of the adjusted models are shown in Table 4. The airshed pathway was strongly significant in the preterm birth model (p < 0.001) and low birthweight model (p = 0.004), even after the addition of the watershed pathway in Model 3. Evidence of watershed mediation was less clear among all outcomes, and among term low birthweight, evidence of mediation via either airshed or watershed pathways was lacking. Table 4 and Figure 4 show whether the effect of the independent variable was reduced after the addition of the mediating variables, outlining the odds ratio and 95% confidence intervals of proximity to surface mining in predicting each birth outcome. Beige lines represent the base model (Model 0) with no mediation variables, while pink represents the addition of the airshed variable alone (Model 1), maroon represents the addition of the watershed variable alone (Model 2), and brown represents the addition of both air and watershed variables together (Model 3). Across all outcomes, the odds of preterm birth, low birthweight, and term low birthweight as explained by proximity to surface mining were reduced with the addition of the airshed mediator, either alone (Model 1) or combined (Model 3). Interestingly, the association between surface mining and birth outcomes became non-significant among preterm birth and low birthweight in the models including the airshed mediator (pink/brown), but remained significant among models with the watershed mediator alone (maroon), with the exception of term low birthweight, which was not significant regardless of the addition of mediators.

Table 4

Figure 4

Previous studies have suggested proximity to surface mining activities is associated with an increased odds of adverse health outcomes,⁷ including adverse birth outcomes,²² and birth defects.^17,21 The direct and indirect causal pathways underlying these associations, however, remain less clear.⁴⁰ Here, we aimed to clarify these associations by exploring hypothesized mediating pathways linking exposure to coal mining activities and resultant birth outcomes. Overall, we found evidence that living within airsheds of active surface mines is associated with higher odds of preterm birth and low birthweight outcomes, with airsheds mediating the effect of living in close proximity to surface mines on these birth outcomes. Specifically, the airshed pathway met all four mediation criteria (Table 1): 1) the airshed pathway was significantly associated with active surface mining (Table 3); 2) proximity to active surface mining was significantly associated with both preterm birth and low birthweight outcomes in the absence of the airshed pathway (Tables S1 and S2, Model 0); the airshed pathway significantly and uniquely was associated with preterm birth and low birthweight outcomes (Tables S1 and S2, Model 1); and, 4) the association between proximity to active surface mining sites and both preterm birth and low birthweight outcomes was reduced with the addition of the airshed pathway (Table 4). Evidence of mediation was less clear for term low birthweight. While criteria 1 and 4 were met, the effect of proximity to active surface mining on term low birthweight was not significant to suggest an association. While mediation was not clear, this outcome was strongly explained by maternal age at birth, race, self-reported tobacco use, payment type for birth services and child sex, with the highest odds of term low birthweight among those reporting tobacco use (OR: 3.06, 2.80 – 3.35).

Evidence for mediation along the watershed pathway across outcomes was less clear. While the watershed pathway was significantly associated with active surface mining, and the effect of surface mining on the dependent variables were generally reduced with its addition, the watershed pathway was not significantly and uniquely associated with the adverse birth outcomes, as evidenced in Table 4 (Model 3). These findings are in line with results from the Sobel test (Table 2), which suggested mediation along the airshed pathway for preterm birth and low birthweight, but not along the watershed pathway.

These findings suggest that the most likely exposure pathway underlying the association between active surface mining and preterm birth and low birthweight outcomes likely occurs via air, potentially influenced by particulate matter emissions from coal mining activities. Evidence suggested mediation via the watershed pathway was less likely, which is unsurprising given the pathway between water pollution and drinking water is more complex than exposure to air pollution. For example, the majority of residents in this region are served by municipal drinking water, which is treated prior to distribution and consumption;⁴¹ adverse changes to source water quality will only have an impact if these treatment systems fail. In addition, it is well-documented that changes in household water quality can alter consumption source patterns: specific to this region, a residential survey conducted in West Virginia found that perceptions of poor water quality within the area resulted in behavior changes among the respondents, such as consuming bottled water over tap water.⁴² Exposure to developmental toxicants via contamination of surface or groundwaters therefore represents a potential indirect exposure pathway, which might be alleviated by water treatment prior to reaching the household tap or alternative drinking water sources (e.g., bottled water).

These analyses should be considered within the context of their associated limitations. Firstly, we did not directly measure air pollutants or water pollutants, but instead used proxy atmospheric models and geological boundaries. For example, we used satellite imagery to designate temporal surface mining activities based on mine permits and barren land cover; while this classification likely represents active surface mining activities because the change in landcover due to surface mining is stark, it does not capture on-the-ground information about mining activities and air or water emissions. Using these proxies therefore does not quantify an individual’s actual exposure to air and water pollutants, heavy metals, etc. and future research could aim to better characterize individual exposure via prospective studies. Further, the time scale of the main exposure variable (proximity to active surface mining) is annual due to limitations in obtaining reliable imagery of surface mining activity on a monthly timescale.²⁴ We calculated an airshed for each mine for every year of the study time period and matched the relevant year’s cumulative airshed frequency value to a birth record’s gestational majority year. Given that we had one mine boundary for each year, we were not able to calculate airsheds at a finer time scale than one year (e.g., at the gestational trimester level). While we don’t expect that mine boundaries shift significantly at a finer time scale, exploring exposure at the gestational trimester represents a promising opportunity for future study.

Secondly, our analyses are limited to women with an identifiable street address (e.g., excluding post office boxes, etc.). It was necessary, however, to restrict our analyses at this level in order to link individual births to surface mining airsheds and watersheds. Our previous analyses using ZIP code level exposures that allowed for inclusion of most birth records in the dataset showed that ZIP code and street-level analyses provided similar results, ²² increasing our confidence that the current analysis is minimally biased due to this limitation. Thirdly, our analyses rely on secondary data sources and are therefore limited to socioeconomic and demographic characteristics existing on birth records. Namely, tobacco use within these records was self-reported, which may not represent actual tobacco use in pregnancy. Further, other important characteristics such as illicit substance use, household income, etc. were not included in our analyses, as these covariates did not exist on birth records, and could be a source of potential bias in the present analysis. However, while these non-observed maternal characteristics could be associated both with mothers' place of residence and birth outcomes, active surface mining is estimated yearly, such that birth outcomes from the same area prior to or after active surface mining are compared to births in which gestation occurs during active mining. Future research could incorporate additional area-level sociodemographic characteristics within a hierarchical framework.

Surface mining activities have dramatically altered the Central Appalachia landscape, and the environmental and public health consequences are still being characterized. While previous literature links proximity to surface mining activities and resultant adverse birth outcomes, the exposure pathways underlying these associations remain poorly understood. Our study represents a novel mediation analysis using nearly 200,000 geolocated birth records spanning four states in Central Appalachia from 1990 – 2015 to explore airshed and watershed exposure pathways as mediators of the relationship between surface mining activities and adverse birth outcomes. Our findings suggest that the airshed pathway mediates the association between surface mining activities and preterm birth and low birthweight outcomes, but evidence of mediation via the watershed pathway was less clear.

Acknowledgements

The authors would like to acknowledge the contributions made to compiling the datasets used in these analyses by the following individuals: M Marston, L Marr, and E Smith. The authors would also like to thank A Mendrinos and B Ramesh for their feedback while these analyses were being conducted.

Data Sharing

Data use agreements with state health agencies prohibit authors from sharing birth record data used in the analysis. Satellite-derived datasets are publicly available, or available on request. Code required to replicate results are available upon request.

Table 1. Mediation criteria, adapted from ³¹
Criteria Assessed		Statistical Test(s) Used
1	Independent variable significantly influences the mediator	Unadjusted regression model	Sobel test
2	Independent variable significantly influences the dependent variable (in absence of the mediator)	Adjusted regression model
3	Mediator significantly and uniquely influences the dependent variable
4	The effect of the independent variable on the dependent variable shrinks with addition of the mediator

Table 2. Sobel test exploring mediation via airshed and watershed pathways from surface mining operations on birth outcomes
	Z-value	Significance
M = air, X = surface mining
Preterm birth	14.20	9.78e-46 ***
Low birthweight	5.66	1.53e-08 ***
Term low birthweight	1.12	0.261
M = water, X = surface mining
Preterm birth	2.522	0.012 *
Low birthweight	1.158	0.247
Term low birthweight	1.857	0.063
Signif. codes: 0 ‘*’ 0.001 ‘’ 0.01 ‘*’

Table 3. Unadjusted model exploring birth outcomes by exposure to surface mining via airshed and watershed mediation
Preterm Birth
Covariates	β	SE	Significance

PTB ~ air	0.138	0.009	< 2e-16 ***
PTB ~ water	0.061	0.015	2.59e-05 ***
Low Birthweight
LBW ~ air	0.069	0.011	7.04e-11 ***
LBW ~ water	0.047	0.017	0.00473 **
Term Low Birthweight
TLBW ~ air	0.028	0.017	0.0928
TLBW ~ water	0.069	0.026	0.00865 **
Mediators
surface mining ~ air	0.228	0.002	<2e-16 ***
surface mining ~ water	0.729	0.002	<2e-16 ***
Signif. codes: 0 ‘*’ 0.001 ‘’ 0.01 ‘*’

Table 4. Selected adjusted model results exploring birth outcomes by exposure to surface mining via airshed and watershed mediation
Preterm Birth
Covariates	Model 0 OR (95% CI)	Model 1 OR (95% CI)	Model 2 OR (95% CI)	Model 3 OR (95% CI)
% Active mining	1.09 (1.05,1.13)	0.69 (0.38,1.22)	1.10 (1.04,1.16)	0.48 (0.14,1.58)
Airshed exposure		1.12 (1.09,1.15)		1.14 (1.11,1.18)
Watershed exposure			1.07 (1.01,1.12)	1.38 (0.45,4.17)
Low Birthweight
% Active mining	1.06 (1.02,1.11)	0.81 (0.42,1.53)	1.06 (1.00,1.13)	0.78 (0.20,3.07)
Airshed exposure		1.05 (1.01,1.08)		1.06 (1.03,1.10)
Watershed exposure			0.99 (0.93,1.06)	0.63 (0.16,2.35)
Term Low Birthweight
% Active mining	1.00 (0.93,1.06)	0.87 (0.28,2.64)	0.98 (0.88,1.09)	0.32 (0.03,3.27)
Airshed exposure		1.02 (0.97,1.07)		1.00 (0.94,1.05)
Watershed exposure			1.04 (0.94,1.14)	0.77 (0.10,5.66)

1. Krometis LA, Gohlke J, Kolivras K, Satterwhite E, Marmagas SW, Marr LC. Environmental health disparities in the Central Appalachian region of the United States. Reviews on Environmental Health. 2017;32(3):253-266. doi:10.1515/reveh-2017-0012

2. Gohlke JM. Human Health in Coalfield Communities of Appalachia. In: Zipper CE, Skousen J, eds. Appalachia’s Coal-Mined Landscapes: Resources and Communities in a New Energy Era. Springer International Publishing; 2021:311-336. doi:10.1007/978-3-030-57780-3_12

3. PDA, INC, APPALACHIAN REGIONAL COMMISSION, THE CECIL G. SHEPS CENTER FOR HEALTH SERVICES RESEARCH. Health Disparities in Appalachia. APPALACHIAN REGIONAL COMMISSION; 2017. Accessed July 28, 2021. https://www.arc.gov/wp-content/uploads/2020/06/Health_Disparities_in_Appalachia_August_2017.pdf

4. US Energy Information Administration. Annual Coal Report 2019. US Department of Energy; 2019:87.

5. U.S. Energy Information Administration (EIA). Coal and the environment. Accessed September 8, 2021. https://www.eia.gov/energyexplained/coal/coal-and-the-environment.php

6. Palmer MA, Bernhardt ES, Schlesinger WH, et al. Mountaintop Mining Consequences. Science. 2010;327(5962):148-149. doi:10.1126/science.1180543

7. Boyles AL, Blain RB, Rochester JR, et al. Systematic review of community health impacts of mountaintop removal mining. Environ Int. 2017;107:163-172. doi:10.1016/j.envint.2017.07.002

8. Hendryx M. The public health impacts of surface coal mining. The Extractive Industries and Society. 2015;2(4):820-826. doi:10.1016/j.exis.2015.08.006

9. Hendryx M, Ahern MM. Relations Between Health Indicators and Residential Proximity to Coal Mining in West Virginia. Am J Public Health. 2008;98(4):669-671. doi:10.2105/AJPH.2007.113472

10. Nichols CE, Baseler WA, Thapa D, et al. Mountain-top mining particulate matter exposure increases markers of mitochondrially-driven apoptosis in rat cardiac tissue. The FASEB Journal. 2012;26(S1):1036.15-1036.15. doi:10.1096/fasebj.26.1_supplement.1036.15

11. Knuckles TL, Stapleton PA, Minarchick VC, et al. Air Pollution Particulate Matter Collected from an Appalachian Mountaintop Mining Site Induces Microvascular Dysfunction. Microcirculation. 2013;20(2):158-169. doi:10.1111/micc.12014

12. Spengler JD, Koutrakis P, Dockery DW, Raizenne M, Speizer FE. Health effects of acid aerosols on North American children: air pollution exposures. Environ Health Perspect. 1996;104(5):492-499. doi:10.1289/ehp.96104492

13. Aneja VP, Isherwood A, Morgan P. Characterization of particulate matter (PM10) related to surface coal mining operations in Appalachia. Atmospheric Environment. 2012;54:496-501. doi:10.1016/j.atmosenv.2012.02.063

14. Luanpitpong S, Chen M, Knuckles T, et al. Appalachian Mountaintop Mining Particulate Matter Induces Neoplastic Transformation of Human Bronchial Epithelial Cells and Promotes Tumor Formation. Environ Sci Technol. 2014;48(21):12912-12919. doi:10.1021/es504263u

15. Dadvand P, Parker J, Bell ML, et al. Maternal Exposure to Particulate Air Pollution and Term Birth Weight: A Multi-Country Evaluation of Effect and Heterogeneity. Environmental Health Perspectives. 2013;121(3):267-373. doi:10.1289/ehp.1205575

16. Gray SC, Edwards SE, Schultz BD, Miranda ML. Assessing the impact of race, social factors and air pollution on birth outcomes: a population-based study. Environmental Health. 2014;13:4. doi:10.1186/1476-069X-13-4

17. Ahern M, Hendryx M, Conley J, Fedorko E, Ducatman A, Zullig KJ. The association between mountaintop mining and birth defects among live births in central Appalachia, 1996-2003. Environ Res. 2011;111(6):838-846. doi:10.1016/j.envres.2011.05.019

18. Luyten LJ, Saenen ND, Janssen BG, et al. Air pollution and the fetal origin of disease: A systematic review of the molecular signatures of air pollution exposure in human placenta. Environmental Research. 2018;166:310-323. doi:10.1016/j.envres.2018.03.025

19. Wang X, Ding H, Ryan L, Xu X. Association between air pollution and low birth weight: a community-based study. Environmental Health Perspectives. 1997;105(5):514-520. doi:10.1289/ehp.97105514

20. Bobak M. Outdoor air pollution, low birth weight, and prematurity. Environmental Health Perspectives. 2000;108(2):173-176. doi:10.1289/ehp.00108173

21. Ahern M, Mullett M, MacKay K, Hamilton C. Residence in Coal-Mining Areas and Low-Birth-Weight Outcomes. Matern Child Health J. 2011;15(7):974-979. doi:10.1007/s10995-009-0555-1

22. Buttling LG, McKnight MX, Kolivras KN, Ranganathan S, Gohlke JM. Maternal proximity to Central Appalachia surface mining and birth outcomes. Environ Epidemiol. 2021;5(1). doi:10.1097/EE9.0000000000000128

23. United States Office of Surface Mining Reclamation and Enforcement. Valid Existing Rights - Proposed Revisions to the Permanent Program Regulations Implementing Section 522(E) of the Surface Mining Control and Reclamation Act of 1977 and Proposed Rulemaking Clarifying the Applicability of Section 522(E) to Subsidence from Underground Mining : Environmental Impact Statement.; 1999. Accessed September 3, 2021. https://onlinebooks.library.upenn.edu/webbin/book/lookupid?key=ha100977853

24. Marston ML, Kolivras KN. Identifying Surface Mine Extent Across Central Appalachia Using Time Series Analysis, 1984-2015. International Journal of Applied Geospatial Research (IJAGR). 2021;12(1):38-52. doi:10.4018/IJAGR.2021010103

25. RATHORE† CS, WRIGHT R. Monitoring environmental impacts of surface coal mining. International Journal of Remote Sensing. 1993;14(6):1021-1042. doi:10.1080/01431169308904394

26. McHugh ML. Interrater reliability: the kappa statistic. Biochemia Medica. 2012;22(3):276-282.

27. Draxler RR, Hess GD. Description of the HYSPLIT4 Modeling System. NOAA Air Resources Laboratory; 1997:31.

28. McKnight MX, Kolivras KN, Buttling LG, et al. Associations Between Surface Mining Airsheds and Birth Outcomes in Central Appalachia at Multiple Spatial Scales. Spatial and Spatio-temporal Epidemiology. Published online in review.

29. United States Department of Agriculture-Natural Resources Conservation Service (USDA-NRCS), United States Geological Survey (USGS), Environmental Protection Agency (EPA). HUC10: USGS Watershed Boundary Dataset of Watersheds. Published April 22, 2017. Accessed May 10, 2021. https://developers.google.com/earth-engine/datasets/catalog/USGS_WBD_2017_HUC10

30. VanderWeele TJ. Mediation Analysis: A Practitioner’s Guide. Annu Rev Public Health. 2016;37(1):17-32. doi:10.1146/annurev-publhealth-032315-021402

31. Preacher KJ, Hayes AF. Asymptotic and resampling strategies for assessing and comparing indirect effects in multiple mediator models. Behavior Research Methods. 2008;40(3):879-891. doi:10.3758/BRM.40.3.879

32. Baron RM, Kenny DA. The moderator-mediator variable distinction in social psychological research: conceptual, strategic, and statistical considerations. J Pers Soc Psychol. 1986;51(6):1173-1182. doi:10.1037//0022-3514.51.6.1173

33. Mackinnon DP, Dwyer JH. Estimating Mediated Effects in Prevention Studies. Eval Rev. 1993;17(2):144-158. doi:10.1177/0193841X9301700202

34. Wang Y. Mixed effects smoothing spline analysis of variance. Journal of the Royal Statistical Society: Series B (Statistical Methodology). 1998;60(1):159-174. doi:10.1111/1467-9868.00115

35. Bozdogan H. Model selection and Akaike’s Information Criterion (AIC): The general theory and its analytical extensions. Psychometrika. 1987;52(3):345-370. doi:10.1007/BF02294361

36. Bosco C, Alegana V, Bird T, et al. Exploring the high-resolution mapping of gender-disaggregated development indicators. Journal of The Royal Society Interface. 2017;14(129):20160825. doi:10.1098/rsif.2016.0825

37. Ruktanonchai CW, Nilsen K, Alegana VA, et al. Temporal trends in spatial inequalities of maternal and newborn health services among four east African countries, 1999–2015. BMC Public Health. 2018;18(1):1339. doi:10.1186/s12889-018-6241-8

38. Smith TJ, McKenna CM. A Comparison of Logistic Regression Pseudo R2 Indices. Multiple Linear Regression Viewpoints. 2013;39(2):10.

39. McFadden D. Conditional logit analysis of qualitative choice behavior. In: Frontiers in Econometrics. P. Zarembka. Academic Press; 1974:105-142. Accessed August 18, 2021. https://eml.berkeley.edu/reprints/mcfadden/zarembka.pdf

40. Lamm SH, Li J, Robbins SA, Dissen E, Chen R, Feinleib M. Are residents of mountain-top mining counties more likely to have infants with birth defects? The west virginia experience. Birth Defects Research Part A: Clinical and Molecular Teratology. 2015;103(2):76-84. doi:10.1002/bdra.23322

41. Hughes J, Whisnant R, Weller L, et al. Drinking Water and Wastewater Infrastructure in Appalachia. Appalachian Regional Commission; :150.

42. McSpirit S, Reid C. Residents’ Perceptions of Tap Water and Decisions to Purchase Bottled Water: A Survey Analysis from the Appalachian, Big Sandy Coal Mining Region of West Virginia. Society & Natural Resources. 2011;24(5):511-520. doi:10.1080/08941920903401432

43. US Energy Information Administration. Annual Coal Report 2015. US Department of Energy; 2016. Accessed September 7, 2021. https://www.eia.gov/coal/annual/archive/05842015.pdf

44. US Energy Information Administration. Coal Industry Annual 1993. US Department of Energy; 1994:287.

AIC Akaike Information Criterion

CI Confidence interval

HYSPLIT4 HYbrid Single-Particle Lagrangian Integrated Trajectory

HUC Hydrologic unit code

KY Kentucky

LBW Low birthweight

MAE Mean adjusted error

MGY Majority gestation year

OR Odds ratio

PTB Pre-term birth

RMSE Root mean square error

TLBW Term low birthweight

TN Tennessee

VA Virginia

WV West Virginia

SupplementalInformationR1.docx

PDF

Version 2

posted 17 Apr, 2022

You are reading this latest preprint version

Identifying exposure pathways mediating adverse birth outcomes near active surface mines in Central Appalachia

Status:

Version 2

Abstract

Figures

Highlights

Introduction

Methods

Results

Discussion

Conclusions

Declarations

Tables

References

Abbreviations

Supplementary Files

Status:

Version 2