Skip to main content

Increase in fertility following coal and oil power plant retirements in California



Few studies have explored the relationship between air pollution and fertility. We used a natural experiment in California when coal and oil power plants retired to estimate associations with nearby fertility rates.


We used a difference-in-differences negative binomial model on the incident rate ratio scale to analyze the change in annual fertility rates among California mothers living within 0-5 km and 5-10 km of 8 retired power plants between 2001 and 2011. The difference-in-differences method isolates the portion of the pre- versus post-retirement contrast in the 0-5 km and 5-10 km bins, respectively, that is due to retirement rather than secular trends. We controlled for secular trends with mothers living 10-20 km away. Adjusted models included fixed effects for power plant, proportion Hispanic, Black, high school educated, and aged > 30 years mothers, and neighborhood poverty and educational attainment.


Analyses included 58,909 live births. In adjusted models, we estimated that after power plant retirement annual fertility rates per 1000 women aged 15–44 years increased by 8 births within 5 km and 2 births within 5-10 km of power plants, corresponding to incident rate ratios of 1.2 (95% CI: 1.1–1.4) and 1.1 (95% CI: 1.0–1.2), respectively. We implemented a negative exposure control by randomly selecting power plants that did not retire and repeating our analysis with those locations using the retirement dates from original 8 power plants. There was no association, suggesting that statewide temporal trends may not account for results.


Fertility rates among nearby populations appeared to increase after coal and oil power plant retirements. Our study design limited the possibility that our findings resulted from temporal trends or changes in population composition. These results require confirmation in other populations, given known methodological limitations of ecologic study designs.

Peer Review reports


Oil and coal power plants can emit significant air pollution, including particulate matter (PM), sulfur dioxide (SO2), benzene, lead, mercury, and other known health hazards [1]. The amount of emitted pollution varies, however, based on fuel usage and type, control technologies (e.g., scrubbers), stack heights, and meteorological conditions [2]. While populations living closest to power plants experience higher levels of primary air pollutants, substantial exposure to secondary particulate matter can occur up to 500 km away [3].

Residential proximity to power plants has been linked to health outcomes such as cardiovascular and respiratory disease [4] and adverse birth outcomes [5, 6]. Coal and oil plants, especially older and less efficient ones, typically emit more pollution than those run on natural gas [7]. Between 2007 and 2016 coal-fired electricity generation dropped by 40% in the U.S [8]. Future coal and oil power plant retirements–driven by energy efficiency, abundant and cheap natural gas, and legislation related to climate change–may benefit the health of those living nearby.

Another possible result of power plant retirement is improved fertility among nearby populations. Some evidence suggests links between environmental factors and infertility [9], but few studies have investigated the role of air pollution on reproductive outcomes such as live birth rates, infertility, and miscarriage [10]. Three studies have reported associations between traffic-related air pollution and longer time-to-pregnancy, infertility, and spontaneous abortion [11,12,13]. Fewer than expected miscarriages and stillbirths were reported after a coal plant retirement in Croatia [14].

Only Nieuwenhuijsen et al. (2014) have examined the relationship between air pollution and fertility rates in a general population. They reported an association between higher levels of traffic-related pollution and coarse fraction particulate matter (PM2.5–10) in Barcelona and reduced fertility [15]. The components of air pollution differ between that generated by traffic and power plants, however, and the authors noted the potential for residual confounding due to their ecological design.

Research suggests that residential proximity to industrial facilities, including power plants, is also associated with adverse mental health outcomes, including stress and depression, due to perceptions of neighborhood disorder and feelings of personal powerlessness [16, 17]. This may have implications for fertility based on evidence linking higher levels of stress with difficulty conceiving [18].

We sought to expand this limited research by capitalizing on a natural experiment when 8 coal and oil power plants retired in California between 2001 and 2011. Natural experiments enable researchers to investigate the impact of modifiable factors or interventions on health outcomes by exploiting changes in exposure that are difficult to manipulate experimentally [19]. Natural experiments effectively randomize exposure, reducing the threat of residual confounding inherent in observational studies [20]. We implemented a quasi-experimental difference-in-differences approach to estimate the association between coal and oil power plant retirements and nearby fertility rates.


Study population

We obtained birth certificate data from January 1, 2001 to December 31, 2011 from the California Department of Public Health. We abstracted infant sex and birth date, maternal address, age, date of last menstrual period (LMP), race/ethnicity, and educational attainment from birth certificates. After exclusions (Additional file 1: Figure S1, available in Supplement), 558,308 live births took place within 20 km of a retiring coal or oil power plant during the study period.

Power plant exposure

We downloaded data on coal and oil power plants in California from the U.S. Energy Information Agency (EIA) [8], the U.S. Environmental Protection Agency (EPA) Air Markets Program [21], and the California EPA Air Resources Board (CARB) [22]. These data contained plant latitude and longitude, fuel type, start date, retirement date, and monthly fuel consumption. CARB also provided annual NOx emissions data on power plants that retired after 2007. We obtained NOx emissions data on the Hunters Point plant from the EPA Air Markets Program database but were unable to locate emissions data on the two smallest power plants that retired in 2002. We included power plants as exposure areas in the analysis if they used coal or oil as a primary fuel (e.g., bituminous coal, distillate fuel oil, and residual fuel oil) and if they retired between 2001 and 2011. Eight California power plants met these criteria (Fig. 1), 6 oil and 2 coal. Two power plants were eventually converted to burn biomass, but this transition required additional permitting and took over a year, meaning they were not active during the study period. For simplicity, we refer to the date of retirement from coal or oil use as “retirement,” even if the plant later transitioned to a new fuel.

Fig. 1
figure 1

Location and names of 8 coal and oil power plants in California with year of retirement

Based on maternal address at the time of delivery, we used QGIS (, QGIS Development Team) to identify births that took place within a 20 km radius of one of the eight power plants. In order to assess a gradient of exposure, we assigned each birth to one of three area bins within the larger 20 km-radius circular area: 0-5 km; 5-10 km; and 10-20 km from a power plant. We used these area bins–0-5 km; 5-10 km; and 10-20 km–and the population of women aged 15–44 years living within them as the units of analysis. Next, we applied temporal restriction criteria. We defined births to mothers with an LMP in the year following power plant retirement as unexposed and births to mothers with an LMP between 1 and 2 years prior to power plant retirement as exposed (Fig. 2). These definitions had the advantage of providing a one-year washout period and of matching exposed and unexposed periods to account for fixed cohort bias, seasonal trends in fertility, and power plant emissions (which peak during the summer months). We anticipated that fertility trends among women living 10-20 km from power plants would be similar to those located ≤10 km in the absence of power plant retirements. Compared to a ≤ 10 km group, the 10-20 km group would experience a smaller change in power plant-related exposure after retirement. Therefore, this 10-20 km group allowed us to control for temporal trends in fertility and isolate changes due to plant retirement.

Fig. 2
figure 2

Identification of exposed and unexposed births based on last menstrual period among California mothers, 2001–2011

Fertility rate

For each of the 8 plants, we estimated the number of women aged 15–44 years living within 20 km before and after retirement using block group level data from the closest, non-overlapping censuses: 2000 U.S. census and 2005–2009, 2006–2010, 2008–2012, 2009–2013, and 2011–2015 American Community Survey (ACS) [23]. For example, the Kings Beach plant retired in December 2002 and we estimated the population before retirement with 2000 census data and the population after retirement with 2005–2009 ACS data. When block groups did not fall entirely in specific area bins, we used geographic-area weighting to allocate the population. We operationalized fertility rate as the number of live births per 1000 women aged 15–44 years within 5 km, 5-10 km, and 10-20 km of power plants before and after retirement.

Neighborhood data

For our primary analysis, we downloaded block group level data on the number of individuals living below the federal poverty threshold and who attained less than a high school education from the 2000 U.S. census and the ACS [23]. Based on mother’s address at the time of birth we linked block group data. In sensitivity analyses we used additional years of ACS data as described below. We also obtained address-level 2005–2011 annual foreclosure data across California from CoreLogic (formerly DataQuick), which we assigned based on mother’s census block group and infant’s birth year.

Statistical analysis

We used a difference-in-differences design [24] to compare fertility rates within 20 km of power plants before and after retirement. The difference-in-differences method isolates the portion of the pre- versus post-retirement contrast in the 0-5 km and 5-10 km bins, respectively, that is due to retirement rather than secular trends [24]. We used a negative binomial model to estimate the association between power plant retirements and annual fertility rates on the incident rate ratio scale [25]. Under the assumption that secular trends are parallel in both groups and that the negative binomial was correctly specified, the resulting difference-in-differences estimator corresponds to log-fertility incidence rate ratio. We obtained the difference-in-differences estimates from the product interaction term of two indicator variables: one for the area bin (i.e., 0-5 km or 5-10 km, with 10-20 km as the reference) and one for exposure (i.e., post-retirement, with pre-retirement as the reference). We hypothesized that we would observe a difference-in-differences estimators greater than 1, representing increased fertility rates in the communities 0-5 km and 5-10 km from the power plants after retirement. We identified predictors of fertility a priori [26,27,28] and controlled for them to increase precision of our model [29, 30]. Adjusted models included fixed effects for power plant, and the proportion births to mothers ≥30 years, mothers of Black or Hispanic race/ethnicity, and mothers with ≤ a high school degree in each plant-area bin. They also included census block group poverty and educational attainment to account for neighborhood characteristics. Finally, we used the estimated parameters of the adjusted negative binomial model to calculate the number of births per 1000 women at the mean values of all covariates. Regression analyses were implemented in Stata version 13.1 (StataCorp LP, College Station, TX) and R Statistical Software, version 3.3.2 (, R Core Team).

Sensitivity analyses

We performed several sensitivity analyses. First, to assess if small numbers were driving the observed associations, we repeated the analysis after omitting Kings Beach and Portola plants, the plants with the fewest women of childbearing age living within 20 km. Second, we stratified the overall analyses by primary fuel type: coal or oil. Third, we implemented a negative exposure control [31] by randomly selecting 8 operational (i.e., non-retiring) power plants in California. We repeated the analyses using the geographic location of these operating plants and by randomly assigning to them the retirement dates of the original 8 retired plants. An association with the negative exposure control would have suggested that observed temporal changes in fertility were unrelated to power plant retirement.

Demographic shifts, particularly those related to the Great Recession, a period of substantial economic downturn, and the housing foreclosure crisis that peaked in California in the summer of 2008 [32,33,34], could have influenced our results. Because block group level ACS data only exists in 5-year windows, we included the annual number of block group level foreclosures as an indicator of neighborhood socioeconomic deprivation and re-ran analyses for the years 2005–2011. We also used block group level data [23] prior to and after power plant retirements to explore changes in the number and composition of people living near plants after they retired. We estimated the change in the number of total, non-Hispanic white residents, and women aged 15–44 years, the percent of population living below the federal poverty threshold, and median household income.

Finally, we linked daily PM2.5 data from the U.S. Environmental Protection Agency Community Multiscale Air Quality Model (CMAQ) to mothers’ block groups to assess changes in air pollution exposures near retiring power plants. We compared average annual PM2.5 levels in the 1–2 years prior to the year after power plant retirements.


Eight California power plants–2 coal and 6 oil–retired between 2001 and 2011 (Additional file 1: Table S1 available in supplement). The U.S. EIA reported that the amount of coal used for electricity generation in California declined over the study period from nearly 1 million short tons per year in 2000 to 539,439 short tons in 2012 [35]. The amount of distillate fuel oil and residual fuel oil also declined (Additional file 1: eFigure S2 available in supplement). Average nitrogen oxide (NOx) emissions from the 8 power plants in our analysis fell from 177 tons in the year prior to retirements to 4 tons per year in the year post-retirement [21, 22]. CMAQ data also suggested larger reductions in PM2.5 after power plant retirement in the < 5 km area bin (median pre vs. post: 12.8 μg/m3 vs. 10.8 μg/m3) compared to the 10-20 km area bin (median pre vs. post: 10.8 μg/m3 vs. 10.2 μg/m3) (Additional file 1: eTable S2).

The main analysis included 58,909 live births. We observed an average annual fertility rate of 45 births per 1000 women aged 15–44 years (SD = 24). Temporal trends in maternal and neighborhood characteristics appeared similar by area bin (Table 1). We did note, however, a greater increase in the proportion of births to Hispanic and non-Hispanic Black mothers in the < 5 km bin compared to the 5-10 km and 10-20 km bins. The number of foreclosures increased across all three exposure buffers in the unexposed period, reflecting temporal changes in the housing market.

Table 1 Maternal and neighborhood characteristics of births occurring in California before and after power plant retirement by area bin around power plants between 2001 and 2011

Unadjusted fertility rates varied by power plant and area bin (Fig. 3). We expected minimal confounding due to the difference-in-differences design, and a comparison of point estimates between the unadjusted and adjusted models suggested that the included covariates did not confound the association (Table 2). In the < 5 km area bin, the incident rate ratio (IRR) from the adjusted model was 1.2 (95% CI: 1.1–1.4); this corresponds to an annual fertility rate of 51 births per 1000 women aged 15–44 years post-retirement and 43 births per 1000 women aged 15–44 years pre-retirement. In the 5-10 km buffer, the IRR was 1.1 (95% CI: 1.0–1.2), corresponding to an annual fertility rate of 47 births per 1000 women aged 15–44 years post-retirement and 45 births per 1000 women aged 15–44 years pre-retirement. We observed little change in annual fertility rates in the 10-20 km area bin; the annual fertility rate in the pre-period was 41 (95% CI: 26–57) compared to 42 (95% CI: 25–59) live births per 1000 women aged 15–44 years in the post-period.

Fig. 3
figure 3

Unadjusted annual fertility rate (number of live births per 1000 women aged 15–44 years) by area bin and power plant before and after power plant retirement. Purple bars represent the closest area bin (< 5 km), orange bars the 5-10 km bin, and green bars the 10-20 km bin. Within area bin, the more saturated bars denote the pre-period and the less saturated bars the post-period. The power plants are ordered by the average number of women aged 15–44 years of age living within 20 km, ranging from 820 women near the Portola plant to 464,599 within 20 km of the Hunters Point plant

Table 2 Change in the annual fertility rate (live births per 1000 women aged 15–44 years) from 2001 to 2011 in California after 8 coal and oil power plant retirements

We completed several robustness checks. When we omitted two power plants with the fewest women of childbearing age living nearby results remained similar (IRR = 1.2 [95% CI: 1.0–1.4], IRR = 1.1 [95% CI: 1.0–1.3] in the 0-5 km and 5-10 km bins, respectively). This suggests small, variable rates did not drive the observed association. We observed little change in adjusted associations when stratifying the overall analysis by primary fuel type, except we no longer observed an association between power plant retirement and fertility rates in the 5-10 km bin for coal power plants (Additional file 1: Table S3). While the 10-20 km bins controlled for secular trends local to each plant, in order to further assess the possible effect of statewide trends on our results we implemented negative exposure control analysis, using power plants that did not retire during the study period. This included biomass, geothermal, hydroelectric, natural gas, oil, and solar plants. When combining the location of these plants and their surrounding births with the retirement dates of the original eight power plants we found no significant changes in fertility rates (IRR = 1.1 [95% CI: 0.8–1.3] and IRR = 1.1 [95% CI: 0.9–1.3] in the 0-5 km and 5-10 km bins, respectively, Additional file 1: Table S4). This result suggested that statewide temporal trends may not have accounted for our results. When we integrated foreclosures occurring in the year of birth in block groups within each bin, estimates remained similar (7 additional births per 1000 women in the 0-5 km bin and 1 additional birth per 1000 women in the 5-10 km bin, Additional file 1: Table S5). Finally, we did not identify any major population shifts between the unexposed and exposed time periods that could explain results (Additional file 1: Table S6).


Our study leveraged a natural experiment of 8 coal and oil power plants retirements and found an association with increased fertility rates nearby. Annual fertility rates increased by 8 births per 1000 women aged 15–44 years in the immediate 5 km radius. We observed a gradient, fertility increased the most within 5 km and to a lesser extent within 5-10 km. Because exposed periods occurred prior to unexposed periods and several plant retirements occurred proximal to the housing foreclosure crisis [32,33,34] temporal trends could have influenced our results. Therefore, we conducted several sensitivity analyses to assess demographic changes over time among residents living near power plants. We found no evidence that these compositional changes explained our findings.

Prior literature has assessed the association between power plant proximity and birth outcomes [5, 6, 36, 37]. We extended this literature, finding an association between power plant retirements and increased fertility. Combustion of petroleum products releases numerous air pollutants, such as SO2, NOx, PM, and trace metals like mercury [38] previously linked to reduced human fertility [11, 15]. Stress and mental health among individuals living near power plants may also have played a role [18].

Several mechanisms could explain our findings of increased fertility in the year following plant retirements. Fecundability (i.e., the probability of conception) may have improved. Even exposure to relatively low levels of PM2.5, PM10, and ozone has been associated with changes in sperm morphology and increases in sperm containing increased DNA fragmentation [39]. In a Czech study of 1916 couples, a 10 μg/m3 increase in PM2.5 level was associated with a 22% reduction in fecundability [40]. Air pollution may lead to oxidative stress or endocrine disruption that could explain observed changes in sperm and conception more broadly. In vitro fertilization studies have also demonstrated associations between some–NO2, PM2.5, and PM10–but not all, air pollutants and disrupted fertilization and implantation [41,42,43].

Prior literature has reported changes of a similar magnitude in fertility associated with air pollution exposures. While traffic and power plant air pollution differ [44], traffic studies provide information about potential mechanisms of action related to increased fertility. Nieuwenhuiisen (2014) found a 13% reduction in census tract fertility rates with IQR increases in PM2.5–10 from traffic in Barcelona, Spain during 2011–2012 [15]. A study using data from 116,430 nurses on self-reported infertility (defined as attempting to become pregnant for one year without success) found an 11% increase in the hazard of infertility among women living < 200 m (versus ≥200 m) from a major roadway [11]. Self-reported data indicated that ovulatory disorder and male infertility may have driven these associations [11]. In Texas and Michigan, among 393 couples attempting to conceive, couples were 3% more likely to become pregnant for each additional 200 m they lived away from a major roadway [12]. These changes in fertility may result from endocrine disruption, oxidative stress, or DNA modifications caused by pollutants [10].

A reduction in spontaneous abortion may also help explain our findings. As many as 10–20% of clinically recognized pregnancies end in a pregnancy loss [45], and research has shown that these rates increase following period of stress, such as natural disasters [46]. Our results are consistent with evidence that reports associations between pollutant exposures and increased fetal loss. One study found reduced miscarriage and stillbirth following the temporary closure of a coal plant in Croatia [14]. Two observational studies (350,000 and 750,000 live births, respectively) in New Jersey and Ohio found elevated odds of stillbirth (versus live birth) with increased exposure to ambient SO2, carbon monoxide, and PM2.5 in the third trimester [47, 48]. In Mongolia, which experiences high levels of air pollution from coal combustion in winter months, researchers found a correlation between ambient PM2.5, NO2 and SO2 levels and fetal death [49].

In addition to changes in air quality, the closure of oil and coal power plants may have changed economic conditions for individuals living nearby and thus affected fertility. For instance, if residential proximity was related to employment at the power plants, closure may have resulted in more job loss among individuals living closer to the power plants (i.e., < 5 km). We do not have data to assess this pathway, but job loss would have likely reduced fertility [50], biasing our results towards the null. In addition, other aspects of power plants could have affected fertility rates. Some studies suggest a link between industrial activities and stress, which have been associated with reduced fertility [16, 18].

We assessed whether compositional changes and economic improvement around power plants after retirement might have explained our results. First, we restricted the unexposed period to the single year after power plant retirement, limiting compositional changes that may have occurred over longer periods of time. Second, using ACS data we found minimal changes in population composition, such as small changes in neighborhood poverty following plant retirements. We also assessed whether changes in the housing market due to the foreclosure crisis affected fertility, unrelated to contemporaneous power plant retirements. Models adjusted for the number of annual foreclosures by area bin did not differ from our original results. This adjustment, however, may have been insufficient if temporal changes in foreclosures differed by area bin pre- and post-retirement [51]. Finally, we used communities located 10-20 km from power plants as controls under the assumption that they might have had similar composition to communities located closer but would experience smaller changes in air pollution due to power plant retirements. The choice of control group may have produced conservative estimates, since pollutants from coal plants can travel over 10 km [1].

Assessment of spatial and temporal exposures has been noted as a methodological concern in studies of air pollution and reproductive health [52]. Our design, utilizing a natural experiment of plant closures, addresses these issues through specification of exposure location and timing. As fertility did not cause the closing of plants, staggered closing dates allowed us to approximate a randomized design using difference-in-differences methods.

Our study had several limitations. Because fertility is governed by both biological and behavioral factors [28, 53], our outcome variable (i.e., the fertility rate) lacks specificity to identify pure biological effects. The estimated benefit of power plant retirement for area-level fertility rate may not apply to individual couples attempting to conceive. In addition, the fertility rate denominator estimate included error. The ACS provides population estimates at the block group level in 5-year windows, not annually, and block groups did not perfectly overlap with our selected area bins. Furthermore, we assumed a constant population density of women 15–44 years across block groups, which also added uncertainty to fertility rate estimates. Several of the plants had a relatively small number of births occur within 20 km, which could have resulted in variable rates and spurious associations. To assess this issue, we repeated analysis after removing the two plants with the fewest women of childbearing age nearby and noted no qualitative change in results. We also lacked air quality measurements at mothers’ homes or places of work meaning we could not account for varying levels of background air pollution or identify a single pollutant as the potential causative agent. However, CMAQ PM2.5 data did suggest more reduction in air pollution closer to power plants after retirement. Future analyses could explore whether other sources of air pollution, for example biomass plants or traffic, contribute to changes in fertility. Two of the eight power plants eventually transitioned to burning biomass, which release air pollutants. Such transitions take time. For example, the Port of Stockton coal plant retired in January of 2011 and did not reopen to burn biomass until February of 2014 [54]. We treated these transitions, therefore, as retirements for the purpose of this study where we compared the year after retirement to the period 1–2 years prior. In addition, we assigned exposure based on the residential address on the birth certificate and could not account for exposure accrued from work, travel, or residential mobility during pregnancy. We were unable to assess the association of plant closures with neighborhood perceptions and maternal stress. Finally, an assumption of our difference-in-differences analysis is that secular trends in fertility in the 10-20 km area bin mirror that of the < 5 km and 5-10 km area bins, conditional on measured covariates. Although the demographics of the distance bins may have exhibited some differences over the study period [51], unadjusted and adjusted estimates were similar, suggesting that any differences in the measured covariates over time did not confound our results.


Our results indicate potential fertility benefits associated with coal and oil power plant retirements in California. These results require confirmation in other populations, given known methodological limitations of ecologic study designs. The two coal-fired generators retired in this study had capacities–54 MW and 24 MW–about 1 order of magnitude lower than the nationwide average coal plant. Therefore, we might expect stronger associations in other locations. The U.S. EIA has identified nearly 7000 MW of coal and oil capacity in 8 states slated for retirement in 2018 alone. The Clean Power Plan [55]–the federal government’s strategy to curb carbon dioxide emissions from electricity generation–and similar climate change policies will also likely lead to more power plant retirements as states strive to meet emissions standards. Our study design limited the possibility that the estimated increase in fertility rates after coal and oil power plant retirements resulted from temporal trends or changes in population composition. Ecologic designs have known methodological limitations, however, so our results require confirmation in other populations. Future plant retirements will afford opportunities to further assess their potential reproductive health benefits.



American Community Survey


California Air Resources Board


Community Multiscale Air Quality Model


incident rate ratio


last menstrual period

NOx :

nitrogen oxide


particulate matter

SO2 :

sulfur dioxide


U.S. Energy Information Agency


U.S. Environmental Protection Agency


  1. Toxic Air: The Case for Cleaning Up Coal-Fired Power Plants []. Accessed 15 Feb 2016.

  2. Levy JI, Baxter LK, Schwartz J. Uncertainty and variability in health-related damages from coal-fired power plants in the United States. Risk Anal. 2009;29:1000–14.

    Article  Google Scholar 

  3. Levy JI, Spengler JD, Hlinka D, Sullivan D, Moon D. Using CALPUFF to evaluate the impacts of power plant emissions in Illinois: model sensitivity and implications. Atmos Environ. 2002;36:1063–75.

    Article  CAS  Google Scholar 

  4. Hoek G, Krishnan RM, Beelen R, Peters A, Ostro B, Brunekreef B, Kaufman JD. Long-term air pollution exposure and cardio- respiratory mortality: a review. Environ Health. 2013;12:43.

    Article  CAS  Google Scholar 

  5. Ha S, Hu H, Roth J, Kan H, Xu X. Associations between residential proximity to power plants and adverse birth outcomes. Am J Epidemiol. 2015;182:215–24.

    Article  Google Scholar 

  6. Mohorovic L. First two months of pregnancy—critical time for preterm delivery and low birthweight caused by adverse effects of coal combustion toxics. Early Hum Dev. 2004;80:115–23.

    Article  CAS  Google Scholar 

  7. Jaramillo P, Griffin WM, Matthews HS. Comparative life-cycle air emissions of coal, domestic natural gas, LNG, and SNG for electricity generation. Environ Sci Technol. 2007;41:6290–6.

    Article  CAS  Google Scholar 

  8. U.S. Energy Information Administration. Electric Power Monthly. Washington, DC. U.S. Department of Energy; 2017.

    Google Scholar 

  9. Woodruff TJ, Janssen SJ, Guillette LJ Jr, Giudice LC. Environmental Impacts on Reproductive Health and Fertility. Cambridge: University Press; 2010. p. 125–44.

    Book  Google Scholar 

  10. Carré J, Gatimel N, Moreau J, Parinaud J, Léandri R. Does air pollution play a role in infertility?: a systematic review. Environ Health. 2017;16:82.

    Article  Google Scholar 

  11. Mahalingaiah S, Hart JE, Laden F, Farland LV, Hewlett MM, Chavarro J, Aschengrau A, Missmer SA. Adult air pollution exposure and risk of infertility in the Nurses' health study II. Hum Reprod. 2016;31:638–47.

    Article  CAS  Google Scholar 

  12. Mendola P, Sundaram R, Louis GMB, Sun L, Wallace ME, Smarr MM, Sherman S, Zhu Y, Ying Q, Liu D. Proximity to major roadways and prospectively-measured time-to-pregnancy and infertility. Sci Total Environ. 2017;576:172–7.

    Article  CAS  Google Scholar 

  13. Green RS, Malig B, Windham GC, Fenster L, Ostro B, Swan S. Residential exposure to traffic and spontaneous abortion. Environ Health Perspect. 2009;117:1939–44.

    Article  CAS  Google Scholar 

  14. Mohorovic L, Petrovic O, Haller H, Micovic V. Pregnancy loss and maternal methemoglobin levels: an indirect explanation of the association of environmental toxics and their adverse effects on the mother and the fetus. Int J Environ Res Public Health. 2010;7:4203–12.

    Article  Google Scholar 

  15. Nieuwenhuijsen MJ, Basagana X, Dadvand P, Martinez D, Cirach M, Beelen R, Jacquemin B. Air pollution and human fertility rates. Environ Int. 2014;70:9–14.

    Article  CAS  Google Scholar 

  16. Boardman JD, Downey L, Jackson JS, Merrill JB, Saint Onge JM, Williams DR. Proximate industrial activity and psychological distress. Popul Environ. 2008;30:3–25.

    Article  Google Scholar 

  17. Downey L, Van Willigen M. Environmental stressors: the mental health impacts of living near industrial activity. J Health Soc Behav. 2005;46:289–305.

    Article  Google Scholar 

  18. Akhter S, Marcus M, Kerber RA, Kong M, Taylor KC. The impact of periconceptional maternal stress on fecundability. Ann Epidemiol. 2016;26:710–6. e717

    Article  Google Scholar 

  19. Craig P, Cooper C, Gunnell D, Haw S, Lawson K, Macintyre S, Ogilvie D, Petticrew M, Reeves B, Sutton M: Using natural experiments to evaluate population health interventions: new Medical Research Council guidance. J Epidemiol Community Health 2012:jech-2011–200375.

  20. Basu S, Meghani A, Siddiqi A. Evaluating the health impact of large-scale public policy changes: classical and novel approaches. Annu Rev Public Health. 2017;38:351–70.

    Article  Google Scholar 

  21. Air Markets Program Data []. Accessed 6 Feb 2017.

  22. Facility Search Engine []. Accessed 12 Jan 2017.

  23. Steven Manson JS. David van riper, and Steven Ruggles: IPUMS National Historical Geographic Information System: version 12.0 [database]. Minneapolis: University of Minnesota; 2017.

    Google Scholar 

  24. Angrist JD, Krueger AB. Empirical strategies in labor economics. Handbook of Labor Economics. 1999;3:1277–366.

    Article  Google Scholar 

  25. Hilbe JM. Negative Binomial Regression. New York: Cambridge University Press; 2011.

  26. Finer LB, Zolna MR. Shifts in intended and unintended pregnancies in the United States, 2001-2008. Am J Public Health. 2014;104(Suppl 1):S43–8.

    Article  Google Scholar 

  27. Sweeney MM, Raley RK. Race, ethnicity, and the changing context of childbearing in the United States. Annu Rev Sociol. 2014;40:539–58.

    Article  Google Scholar 

  28. Bongaarts J. A framework for analyzing the proximate determinants of fertility. Popul Dev Rev. 1978:105–32.

  29. Colantuoni E, Rosenblum M. Leveraging prognostic baseline variables to gain precision in randomized trials. Stat Med. 2015;34:2602–17.

    Article  Google Scholar 

  30. Zhang M, Tsiatis AA, Davidian M. Improving efficiency of inferences in randomized clinical trials using auxiliary covariates. Biometrics. 2008;64:707–15.

    Article  Google Scholar 

  31. Lipsitch M, Tchetgen Tchetgen E, Cohen T. Negative controls: a tool for detecting confounding and bias in observational studies. Epidemiology. 2010;21:383–8.

    Article  Google Scholar 

  32. Sobotka T, Skirbekk V, Philipov D. Economic recession and fertility in the developed world. Popul Dev Rev. 2011;37:267–306.

    Article  Google Scholar 

  33. Currie J, Tekin E. Is there a link between foreclosure and health? Am Econ J Econ Policy. 2015;7:63–94.

    Article  Google Scholar 

  34. Downing J. The health effects of the foreclosure crisis and unaffordable housing: a systematic review and explanation of evidence. Soc Sci Med. 2016;162:88–96.

    Article  Google Scholar 

  35. U.S. Energy Information Administration. Coal Data Browser. Accessed 29 Dec 2017.

  36. Severnini E. Impacts of nuclear plant shutdown on coal-fired power generation and infant health in the Tennessee Valley in the 1980s. Nat Energy. 2017;2:17051.

    Article  Google Scholar 

  37. Yang M, Bhatta RA, Chou SY, Hsieh CI. The impact of prenatal exposure to power plant emissions on birth weight: evidence from a Pennsylvania power plant located upwind of New Jersey. J Policy Anal Manage. 2017;36:557–83.

    Article  Google Scholar 

  38. Pirrone N, Cinnirella S, Feng X, Finkelman R, Friedli H, Leaner J, Mason R, Mukherjee A, Stracher G, Streets D. Global mercury emissions to the atmosphere from anthropogenic and natural sources. Atmos Chem Phys. 2010;10:5951–64.

    Article  CAS  Google Scholar 

  39. Lafuente R, García-Blàquez N, Jacquemin B, Checa MA. Outdoor air pollution and sperm quality. Fertil Steril. 2016;106:880–96.

    Article  CAS  Google Scholar 

  40. Slama R, Bottagisi S, Solansky I, Lepeule J, Giorgis-Allemand L, Sram R. Short-term impact of atmospheric pollution on fecundability. Epidemiology. 2013;24:871–9.

    Article  Google Scholar 

  41. Legro RS, Sauer MV, Mottla GL, Richter KS, Li X, Dodson WC, Liao D. Effect of air quality on assisted human reproduction. Hum Reprod. 2010;25:1317–24.

    Article  Google Scholar 

  42. Perin PM, Maluf M, Czeresnia CE, Januário DANF, Saldiva PHN. Effects of exposure to high levels of particulate air pollution during the follicular phase of the conception cycle on pregnancy outcome in couples undergoing in vitro fertilization and embryo transfer. Fertil Steril. 2010;93:301–3.

    Article  Google Scholar 

  43. Perin PM, Maluf M, Czeresnia CE, Januario DA, Saldiva PH. Impact of short-term preconceptional exposure to particulate air pollution on treatment outcome in couples undergoing in vitro fertilization and embryo transfer (IVF/ET). J Assist Reprod Genet. 2010;27:371–82.

    Article  Google Scholar 

  44. Brunekreef B, Holgate ST. Air pollution and health. Lancet. 2002;360:1233–42.

    Article  CAS  Google Scholar 

  45. Ammon Avalos L, Galindo C, Li DK. A systematic review to calculate background miscarriage rates using life table analysis. Birth Defects Res A Clin Mol Teratol. 2012;94:417–23.

    Article  CAS  Google Scholar 

  46. Evans RW, Hu Y, Zhao Z. the fertility effect of catastrophe: US hurricane births. J Popul Econ. 2010;23:1–36.

    Article  Google Scholar 

  47. Faiz AS, Rhoads GG, Demissie K, Kruse L, Lin Y, Rich DQ. Ambient air pollution and the risk of stillbirth. Am J Epidemiol. 2012;176:308–16.

    Article  Google Scholar 

  48. DeFranco E, Hall E, Hossain M, Chen A, Haynes EN, Jones D, Ren S, Lu L, Muglia L. Air pollution and stillbirth risk: exposure to airborne particulate matter during pregnancy is associated with fetal death. PLoS One. 2015;10:e0120594.

    Article  Google Scholar 

  49. Enkhmaa D, Warburton N, Javzandulam B, Uyanga J, Khishigsuren Y, Lodoysamba S, Enkhtur S, Warburton D. Seasonal ambient air pollution correlates strongly with spontaneous abortion in Mongolia. BMC Pregnancy Childbirth. 2014;14:146.

    Article  Google Scholar 

  50. Del Bono E, Weber A, Winter-Ebmer R. Clash of career and family: fertility decisions after job displacement. J Eur Econ Assoc. 2012;10:659–83.

    Article  Google Scholar 

  51. Buckley JP, Doherty BT, Keil AP, Engel SM. Statistical approaches for estimating sex-specific effects in endocrine disruptors research. Environ Health Perspect. 2017;125:067013.

    Article  Google Scholar 

  52. Woodruff TJ, Parker JD, Darrow LA, Slama R, Bell ML, Choi H, Glinianaia S, Hoggatt KJ, Karr CJ, Lobdell DT, Wilhelm M. Methodological issues in studies of air pollution and reproductive health. Environ Res. 2009;109:311–20.

    Article  CAS  Google Scholar 

  53. Bongaarts J. The fertility-inhibiting effects of the intermediate fertility variables. Stud Fam Plan. 1982;13:179–89.

    Article  CAS  Google Scholar 

  54. DTE Energy's Northern California biomass plant begins operations []. Accessed 29 Dec 2017.

  55. Complying with President Trump's Executive Order on Energy Independence []. Accessed 29 Dec 2017.

Download references


We thank Elena Krieger and Kristina Dang for their assistance in assembling environmental data and Ralph Catalano for his advice regarding the analysis.


This work was supported in part by the University of California, San Francisco Preterm Birth Initiative. Additional support was provided by the National Institute of Environmental Health Sciences through grant K99ES027023 (Dr. Casey) and grants P01ES022841 and R01ES027051 (Dr. Morello-Frosch), and the grant U.S. EPA RD-83543301 (Dr. Morello-Frosch). This publication was made possible in part by support from the Berkeley Research Impact Initiative (BRII) sponsored by the UC Berkeley Library.

Availability of data and materials

The power plant data used in this article is available from the U.S. EPA at their Air Markets Program website and the U.S. Energy Information Administration through form 923. The birth data used in this study are available from the California Department of Public Health (CDPH). Data are available from the authors upon reasonable request and with permission of the CDPH.

Author information

Authors and Affiliations



JAC, ELO, and RMF designed the study. JAC analyzed the data, drafted the manuscript, and revised drafts. AG and DK contributed to writing the Introduction and Discussion. AG, ELO, and DEG provided statistical support. All authors guided the analyses and reviewed and revised the manuscript. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Joan A. Casey.

Ethics declarations

Ethics approval and consent to participate

The State of California, Committee for the Protection of Humans Subjects and the Committee for Protection of Humans Subjects at the University of California, Berkeley approved this study (Protocol # 2013–10-5693).

Competing interests

The authors declare that they have no competing interests.

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Additional file

Additional file 1:

Figure S1. Flow chart of study population assembly. Figure S2. Changes in coal and oil used for electricity generation in California between 2000 and 2012. Table S1. Characteristics of the 8 coal and oil power plants that retired between 2001 and 2011. Table S2. Change in median and interquartile range PM2.5 concentrations by area bin. Table S3. Analysis stratified by fuel type for the change in the annual fertility rate between 2001 and 2011in California after 8 coal and oil power plant retirements. Table S4. Negative control analysis for the change in the annual fertility rate between 2001 and 2011in California after 8 coal and oil power plant retirements. Table S5. Change in the annual fertility rate between 2001–2011in California after 8 coal and oil power plant retirements, further adjusted for annual housing foreclosures. Table S6. Demographic changes comparing before and after coal and oil power plant retirements in California, between 2001–2011. (DOCX 379 kb)

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

Casey, J.A., Gemmill, A., Karasek, D. et al. Increase in fertility following coal and oil power plant retirements in California. Environ Health 17, 44 (2018).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: