Abstract

With temperatures expected to increase because of climate change, it is essential to study the health outcomes of elevated temperature in vulnerable populations, such as expectant mothers. In this study, the authors estimated the association between heat and humidity, as measured by apparent temperature, and preterm delivery. They conducted a case-crossover analysis of almost 60,000 births spanning 16 counties in California that occurred from 1999 to 2006 between May and September. The authors identified cases of preterm birth from a state registry of births, which were combined with meteorologic and air pollution monitoring data based on residential zip code. High ambient temperature was significantly associated with preterm birth for all mothers, regardless of maternal racial/ethnic group, maternal age, maternal education, or sex of the infant. Results indicated that an 8.6% increase (95% confidence interval: 6.0, 11.3) in preterm delivery was associated with a 10°F (5.6°C) increase in weekly average (lag06) apparent temperature. Greater associations were observed for younger mothers, blacks, and Asians. These associations were independent of air pollutants. Given the significant associations for apparent temperature and preterm delivery found in this study, more large-scale studies of temperature and preterm delivery are warranted.

Preterm delivery of an infant before 37 gestational weeks occurred in 12.8% of all livebirths in the United States in 2006 (1). Preterm infants are more likely to experience adverse health outcomes, such as death (2), respiratory distress during hospitalization at birth, and hospital readmission after birth (3). Families of preterm infants also bear a substantial economic burden, including inpatient hospital costs (4) and long-term health care costs during the first year (5) and through early childhood (6, 7). Infants born preterm can suffer lifelong impacts including decreased motor and cognitive functioning, increased behavioral disorders, and lower educational attainment in adulthood when compared with those born at full term (8, 9). Thus, researchers need to identify risk factors for preterm delivery and work toward prevention.

The etiology of the majority of preterm births remains unknown. Maternal risk factors that have been identified include chronic infections and hypertension (10). More recently, preterm delivery has been linked to environmental factors such as air pollution and traffic exposure (11–16). In a recent investigation of preterm delivery, a peak was observed in August and September (17). Only a few studies have focused on the associations between temperature or other meteorologic variables and preterm delivery, and the findings from these studies have been inconclusive (18, 19).

We studied the association between elevated ambient temperature and preterm delivery during the warm season in California from 1999 to 2006. In addition, several demographic subgroups were examined.

MATERIALS AND METHODS

Exposure classification

Meteorologic data consisting of daily mean, maximum, and minimum temperatures and relative humidity were provided by the California Irrigation Management Information System (20) and the US Environmental Protection Agency Air Quality System (21) for the warm season, designated as May 1 to September 30, from 1999 to 2006. Temperatures in California generally rise sharply in May and June, level off in July and August, and start to decline in September. Apparent temperature was calculated using the following formula: −2.653 + (0.994 × temperature in °C) + 0.0153 × (dew-point temperature in °C)2. Each mother was assigned the resulting apparent temperature value from the weather monitor closest to the reported residential zip code centroid at the time she gave birth. Only those cases residing in zip codes with centroids located within 10 km of a weather monitor, as identified by Hawth's Tools for ArcGIS 9.3 (22), were eligible for this study.

Since previous investigators have reported an association between preterm delivery and several air pollutants, such as particulate matter less than 2.5 μm in aerodynamic diameter (PM2.5) (15), sulfur dioxide (13), and carbon monoxide (12), potential confounding or effect modification by these pollutants was considered. Nitrogen dioxide was also examined, since traffic measures were found to be associated with preterm delivery in some studies (16, 23), and nitrogen dioxide serves as a proxy for traffic (16). Potential confounding or effect modification by ozone was assessed, since ozone levels are generally elevated during the warm season in California, and ozone has been found to be associated with intrauterine growth retardation (24).

Mean daily data from air pollution monitors were provided by the California Air Resources Board (25). The pollution monitor closest to each mother's residential zip code centroid was used to assign exposure, restricting cases to those residing within 10 km of an air pollution monitor. The daily 1-hour maximum was used for all pollutants except PM2.5, which was based on a 24-hour average.

Outcome measure

California's Office of Vital Records supplied information on infants born between 1999 and 2006 (26). To focus on heat exposures, births were restricted to those occurring from May 1 to September 30. Length of gestation was abstracted from birth certificate data, with preterm delivery defined as births from 20 to 36 gestational weeks. Other variables collected from the birth certificate included date of infant's birth, infant's sex, mother's age, mother's years of completed education, mother's race/ethnicity, and mother's residential zip code at the time she gave birth. Only singleton births in counties with at least 20 preterm births for which data on the demographic variables of interest were available were included in the analyses. Deliveries induced prematurely because of pregnancy complications were excluded. Births at an improbable gestational age and/or birth weight were also excluded, using cutoffs as described previously (27).

Statistical analysis

The case-crossover method, a statistical technique ideally suited to examine short-term exposures with acute outcomes, was used for the analysis. The case-crossover study is a modification of the matched case-control study, where each person serves as his or her own control so that known and unknown time-invariant confounders, such as body mass index, are inherently adjusted for by study design. In this study, exposures for up to 1 week before a birth were compared with exposures for the same mother at other times in the infant's birth month and year. Confining the control periods to the month of birth reduces the influence of long-term time trends. Control periods were also limited to the same day of the week as for the case to adjust for day of the week. Consequently, there could be a maximum of 4 control periods per case occurring a minimum of 7 days and a maximum of 28 days before or after the case period.

Two sampling schemes using single-day lags are demonstrated in Figure 1. If a case of preterm delivery occurred in the middle of the month, control periods would have occurred both before and after the case period, as shown in the first example. In the second scenario, a case of preterm delivery occurred at the end of the month, so that all control periods occurred prior to the case. Similarly, if a case of preterm delivery occurred at the beginning of the month, all control periods would be selected after the case. If a case occurred on the first of the month, the lag period could be selected from the last week of the previous month. This sampling scheme is known as the time-stratified design, the most current methodology used in studies with the case-crossover design; bidirectional and unidirectional sampling of control days may introduce bias (28). For the time-stratified design to be valid, it is essential for control periods to be selected at random with respect to when the case of preterm delivery occurred. As mentioned above, this design may use postevent control periods, even though a pregnant woman is theoretically no longer at risk of giving birth. However, any bias from using control periods after the risk period is small if the event is rare (29). In our birth cohort, giving birth on a specific gestational day was a rare occurrence, never exceeding a 4% probability.

Figure 1.

Schematic diagram of a time-stratified case-crossover study design. C1, C2, C3, C4: control periods 1, 2, 3, and 4; Case: case period; T: time that the case (preterm birth) occurred; T − 7, T − 14, T − 21, T − 28, T + 7, T + 14, T + 21: time that the control periods occurred (7, 14, 21, and 28 days prior to the case period and 7, 14, and 21 days following the case period, respectively).

Figure 1.

Schematic diagram of a time-stratified case-crossover study design. C1, C2, C3, C4: control periods 1, 2, 3, and 4; Case: case period; T: time that the case (preterm birth) occurred; T − 7, T − 14, T − 21, T − 28, T + 7, T + 14, T + 21: time that the control periods occurred (7, 14, 21, and 28 days prior to the case period and 7, 14, and 21 days following the case period, respectively).

A linear term for apparent temperature was included in a conditional logistic regression model, and the log(odds) of preterm delivery (yes/no) served as the outcome measure. Nonlinear associations at higher temperatures were explored (i.e., adding a term for the 90th percentile for apparent temperature to the model), but such terms were not significant and were excluded from subsequent analyses. All analyses were performed in 2 steps: first, we calculated the county-level estimate based on residential zip codes; second, we combined the county-level estimates to produce an overall estimate using meta-analytical techniques (30). In this paper, all estimates are reported as percent change per 10°F (5.6°C) increase in apparent temperature with corresponding 95% confidence intervals. Estimates with 95% confidence intervals that do not include zero are considered statistically significant. PROC PHREG in SAS version 9.1 software (31) was used to conduct the first stage of the analysis, and version 2.6.2 of the R program (32) was used for the meta-analyses.

Several lag periods of mean, maximum, and minimum apparent temperature were examined to determine which metric demonstrated the best model fit according to t tests. The exposure periods tested included same-day lag (lag 0), 6 single-day lags (lag 1, lag 2, lag 3, lag 4, lag 5, and lag 6), and 3 cumulative average lag periods (lag01, lag03, lag06). The single-day lags referred to apparent temperature on the previous day (lag 1), 2 days prior (lag 2), and so forth, whereas the cumulative lag periods were averages of several days, such as the exposure on the same day and the previous 3 days (lag03).

We used the apparent temperature metric and lag period with the best model fit to examine potentially vulnerable subgroups of high ambient temperature and preterm delivery. In separate models, we stratified by maternal racial/ethnic group (non-Hispanic white, non-Hispanic black, Hispanic, non-Hispanic Asian), maternal age (<20, <25, 25–34, ≥35 years), infant's sex (male, female), and highest level of maternal education completed at the time of giving birth (high school graduation or less vs. some college). We also examined whether effect estimates differed between earlier (20–33 gestational weeks) and later (34–36 gestational weeks) preterm births, for all full-term births (37–44 gestational weeks), and specifically for 37 and 38 gestational weeks. Differential impacts among demographic groups were evaluated using t tests (33).

We assessed the potential confounding or effect modification of the apparent temperature and preterm delivery association by air pollutants. Bivariate models with a weekly average of each air pollutant considered in addition to the weekly average apparent temperature term were used to assess potential confounding. An interaction term incorporating apparent temperature and each air pollutant was added when assessing potential effect modification.

Effect estimates were compared by month and by year to examine whether any month or year influenced the overall estimate. For the same reason, we also compared effect estimates between inland (Fresno, Kern, Merced, Riverside, Sacramento, San Bernardino, San Joaquin, Tulare) and coastal (Alameda, Contra Costa, Los Angeles, Orange, San Diego, Santa Clara, Solano, Ventura) areas and between northern and southern (designated as Kern and San Bernardino counties and below) California. Finally, we restricted the control periods to 2 weeks instead of 1 month since baseline risk for preterm delivery changes by gestational week.

Before this study was conducted, the research protocol was approved by the California Department of Public Health Internal Review Board.

RESULTS

The study population included 58,681 preterm births throughout 16 counties in California from May to September, 1999–2006 (Table 1). Mothers who were Hispanic (55%), who were between 25 and 34 years of age (47%), or whose educational attainment had not exceeded high school (62%) constituted most of the study population. In Table 2, the overall study mean, 5th percentile, and 95th percentile values for case mean apparent temperature (overall mean: 70.4°F (1°F = 0.56°C)), case maximum apparent temperature (88.7°F), case minimum temperature (57.6°F), carbon monoxide (1.0 ppm), nitrogen dioxide (37 ppb), PM2.5 (15.3 μg/m3), ozone (63 ppb), and sulfur dioxide (5 ppm) are listed by county. The ranges for mean and maximum apparent temperature were larger than those for minimum apparent temperature. Although the probability of birth rises with increasing gestation, the design is valid if the probability of giving birth does not change throughout the month at the aggregate level. We found that approximately the same number of births occurred every week (21% the first week, 22% the second week, 22% the third week, and 25% the fourth week).

Table 1.

Demographics (%) of the Study Population, by California County, May–September, 1999–2006

County No. of Preterm Births Maternal Race/Ethnicity
 
Maternal Age, Years
 
Maternal Education
 
Infant Sex Male 
White Black Hispanic Asian <20 <25 25–34 ≥35 ≤High School Some College 
Alameda 2,812 19 23 32 24 10 29 50 21 53 47 55 
Contra Costa 1,107 40 35 14 27 49 24 49 51 54 
Fresno 3,100 24 10 54 15 44 43 14 67 33 53 
Kern 2,570 24 63 18 46 42 12 72 28 54 
Los Angeles 20,900 13 10 65 11 12 33 47 19 65 35 54 
Merced 422 23 61 19 51 37 12 73 27 55 
Orange 3,560 20 58 18 10 30 49 21 61 39 53 
Riverside 2,828 27 61 13 39 46 16 67 33 54 
Sacramento 3,165 43 15 26 12 11 36 47 17 57 43 54 
San Bernardino 4,869 21 15 58 14 42 44 14 70 30 55 
San Diego 6,006 32 41 13 10 33 49 19 53 47 55 
San Joaquin 2,240 22 11 45 18 16 42 43 15 67 33 55 
Santa Clara 1,231 30 36 29 23 53 24 35 65 56 
Solano 474 31 19 33 14 13 36 46 17 57 43 52 
Tulare 998 24 67 20 49 40 10 71 29 56 
Ventura 2,399 31 55 11 31 48 21 55 45 54 
Total 58,681 22 10 55 11 12 35 47 18 62 38 52 
County No. of Preterm Births Maternal Race/Ethnicity
 
Maternal Age, Years
 
Maternal Education
 
Infant Sex Male 
White Black Hispanic Asian <20 <25 25–34 ≥35 ≤High School Some College 
Alameda 2,812 19 23 32 24 10 29 50 21 53 47 55 
Contra Costa 1,107 40 35 14 27 49 24 49 51 54 
Fresno 3,100 24 10 54 15 44 43 14 67 33 53 
Kern 2,570 24 63 18 46 42 12 72 28 54 
Los Angeles 20,900 13 10 65 11 12 33 47 19 65 35 54 
Merced 422 23 61 19 51 37 12 73 27 55 
Orange 3,560 20 58 18 10 30 49 21 61 39 53 
Riverside 2,828 27 61 13 39 46 16 67 33 54 
Sacramento 3,165 43 15 26 12 11 36 47 17 57 43 54 
San Bernardino 4,869 21 15 58 14 42 44 14 70 30 55 
San Diego 6,006 32 41 13 10 33 49 19 53 47 55 
San Joaquin 2,240 22 11 45 18 16 42 43 15 67 33 55 
Santa Clara 1,231 30 36 29 23 53 24 35 65 56 
Solano 474 31 19 33 14 13 36 46 17 57 43 52 
Tulare 998 24 67 20 49 40 10 71 29 56 
Ventura 2,399 31 55 11 31 48 21 55 45 54 
Total 58,681 22 10 55 11 12 35 47 18 62 38 52 
Table 2.

Exposure Metrics for Weekly Average (Lag06) Apparent Temperature and Air Pollutants, by California County, May–September, 1999–2006

County Case Mean Apparent Temperature, °Fa (5th, 95th Percentile) Case-Control Difference in Mean Apparent Temperature, °F Case Maximum Apparent Temperature, °F (5th, 95th Percentile) Case Minimum Apparent Temperature, °F (5th, 95th Percentile) Case Mean O3 Level, ppb (5th, 95th Percentile) Case Mean CO Level, ppm (5th, 95th Percentile) Case Mean NO2 Level, ppb (5th, 95th Percentile) Case Mean SO2 Level, ppb (5th, 95th Percentile) Case Mean PM2.5 Level, μg/m3 (5th, 95th Percentile) 
Alameda 61.3 (53.2, 68.0) 2.5 77.9 (64.7, 91.7) 50.2 (43.6, 55.5) 35 (19, 52) 0.7 (0.4, 1.3) 23 (12, 36) 3 (1, 6)  
Contra Costa 64.0 (54.2, 74.3) 2.9 84.8 (68.1, 100.5) 49.5 (42.5, 55.0) 51 (37, 71) 0.5 (0.3, 0.8) 17 (10, 28) 5 (1, 12)  
Fresno 75.7 (62.5, 87.5) 4.4 95.9 (78.4, 109.1) 57.9 (47.6, 68.4) 80 (57, 104) 0.6 (0.2, 1.1) 28 (16, 44) 3 (1, 5) 10.6 (6.3, 16.4) 
Kern 76.9 (62.3, 89.3) 4.4 93.8 (75.6, 107.7) 61.1 (48.6, 73.6) 79 (60, 98) 1 (0.5, 1.6) 40 (19, 58) 5 (3, 7) 12.1 (8.3, 16.2) 
Los Angeles 70.3 (59.7, 80.6) 2.9 87.6 (70.7, 104.3) 59.0 (49.9, 66.6) 60 (40, 85) 1.2 (0.5, 2.1) 46 (27, 73) 6 (1, 14) 18.1 (10, 27.9) 
Merced 69.9 (58.5, 81.4) 3.3 95.0 (79.2, 107.7) 49.6 (42.2, 58.3) 79 (57, 99)  23 (13, 36)   
Orange 70.3 (59.6, 80.2) 2.5 86.2 (70.6, 99.6) 59.4 (50.4, 68.4) 62 (46, 78) 0.7 (0.3, 1.4) 27 (14, 44) 3 (1, 6) 14.4 (9, 21.5) 
Riverside 77.1 (61.3, 94.4) 3.9 98.9 (78.8, 113.7) 59.9 (48.0, 76.0) 82 (58, 105) 1.1 (0.3, 2.3) 33 (16, 56) 4 (0, 11) 24.9 (14.2, 38.3) 
Sacramento 71.6 (59.9, 82.0) 4.1 93.6 (75.7, 107.4) 55.3 (47.1, 63.0) 64 (46, 86) 0.7 (0.3, 1.5) 24 (12, 43) 4 (1, 9) 8.5 (4.9, 14.6) 
San Bernardino 73.9 (59.7, 85.7) 3.8 96.2 (75.4, 111.6) 57.7 (47.9, 67.5) 84 (57, 109) 1 (0.4, 1.7) 50 (32, 70) 3 (1, 7) 23.7 (13.6, 37.7) 
San Diego 68.4 (58.3, 78.0) 2.4 83.4 (67.3, 101.5) 58.4 (48.2, 66.9) 52 (38, 67) 1 (0.5, 1.7) 25 (15, 37) 6 (2, 12) 12.9 (7.6, 18.1) 
San Joaquin 70.7 (59.0, 81.2) 3.8 91.9 (75.2, 104.3) 54.6 (45.1, 64.0) 58 (43, 78) 0.7 (0.4, 1.2) 27 (17, 43)   
Santa Clara 63.9 (53.6, 71.5) 2.7 83.0 (67.4, 96.5) 50.7 (41.5, 56.6) 43 (26, 66) 1.1 (0.6, 1.6) 33 (19, 50)   
Solano 65.1 (55.7, 74.9) 2.9 86.4 (70.7, 100.7) 49.8 (41.8, 55.3) 46 (34, 61) 0.6 (0.3, 0.9) 17 (10, 26) 3 (1, 6)  
Tulare 74.4 (60.7, 85.5) 4.1 96.8 (77.7, 110.1) 55.6 (45.2, 65.4) 80 (59, 100) 0.7 (0.3, 1.2) 31 (20, 48)   
Ventura 65.0 (56.6, 73.4) 2.5 77.7 (65.6, 91.7) 54.6 (46.8, 62.1) 59 (43, 84) 0.6 (0.3, 1.3) 21 (13, 38) 3 (0, 11)  
Overall 70.4 (58.1, 83.1) 3.1 88.7 (70.2, 106.7) 57.6 (47.1, 67.2) 63 (37, 95) 1.0 (0.3, 1.9) 37 (15, 65) 5 (1, 12) 15.3 (7.5, 26.7) 
County Case Mean Apparent Temperature, °Fa (5th, 95th Percentile) Case-Control Difference in Mean Apparent Temperature, °F Case Maximum Apparent Temperature, °F (5th, 95th Percentile) Case Minimum Apparent Temperature, °F (5th, 95th Percentile) Case Mean O3 Level, ppb (5th, 95th Percentile) Case Mean CO Level, ppm (5th, 95th Percentile) Case Mean NO2 Level, ppb (5th, 95th Percentile) Case Mean SO2 Level, ppb (5th, 95th Percentile) Case Mean PM2.5 Level, μg/m3 (5th, 95th Percentile) 
Alameda 61.3 (53.2, 68.0) 2.5 77.9 (64.7, 91.7) 50.2 (43.6, 55.5) 35 (19, 52) 0.7 (0.4, 1.3) 23 (12, 36) 3 (1, 6)  
Contra Costa 64.0 (54.2, 74.3) 2.9 84.8 (68.1, 100.5) 49.5 (42.5, 55.0) 51 (37, 71) 0.5 (0.3, 0.8) 17 (10, 28) 5 (1, 12)  
Fresno 75.7 (62.5, 87.5) 4.4 95.9 (78.4, 109.1) 57.9 (47.6, 68.4) 80 (57, 104) 0.6 (0.2, 1.1) 28 (16, 44) 3 (1, 5) 10.6 (6.3, 16.4) 
Kern 76.9 (62.3, 89.3) 4.4 93.8 (75.6, 107.7) 61.1 (48.6, 73.6) 79 (60, 98) 1 (0.5, 1.6) 40 (19, 58) 5 (3, 7) 12.1 (8.3, 16.2) 
Los Angeles 70.3 (59.7, 80.6) 2.9 87.6 (70.7, 104.3) 59.0 (49.9, 66.6) 60 (40, 85) 1.2 (0.5, 2.1) 46 (27, 73) 6 (1, 14) 18.1 (10, 27.9) 
Merced 69.9 (58.5, 81.4) 3.3 95.0 (79.2, 107.7) 49.6 (42.2, 58.3) 79 (57, 99)  23 (13, 36)   
Orange 70.3 (59.6, 80.2) 2.5 86.2 (70.6, 99.6) 59.4 (50.4, 68.4) 62 (46, 78) 0.7 (0.3, 1.4) 27 (14, 44) 3 (1, 6) 14.4 (9, 21.5) 
Riverside 77.1 (61.3, 94.4) 3.9 98.9 (78.8, 113.7) 59.9 (48.0, 76.0) 82 (58, 105) 1.1 (0.3, 2.3) 33 (16, 56) 4 (0, 11) 24.9 (14.2, 38.3) 
Sacramento 71.6 (59.9, 82.0) 4.1 93.6 (75.7, 107.4) 55.3 (47.1, 63.0) 64 (46, 86) 0.7 (0.3, 1.5) 24 (12, 43) 4 (1, 9) 8.5 (4.9, 14.6) 
San Bernardino 73.9 (59.7, 85.7) 3.8 96.2 (75.4, 111.6) 57.7 (47.9, 67.5) 84 (57, 109) 1 (0.4, 1.7) 50 (32, 70) 3 (1, 7) 23.7 (13.6, 37.7) 
San Diego 68.4 (58.3, 78.0) 2.4 83.4 (67.3, 101.5) 58.4 (48.2, 66.9) 52 (38, 67) 1 (0.5, 1.7) 25 (15, 37) 6 (2, 12) 12.9 (7.6, 18.1) 
San Joaquin 70.7 (59.0, 81.2) 3.8 91.9 (75.2, 104.3) 54.6 (45.1, 64.0) 58 (43, 78) 0.7 (0.4, 1.2) 27 (17, 43)   
Santa Clara 63.9 (53.6, 71.5) 2.7 83.0 (67.4, 96.5) 50.7 (41.5, 56.6) 43 (26, 66) 1.1 (0.6, 1.6) 33 (19, 50)   
Solano 65.1 (55.7, 74.9) 2.9 86.4 (70.7, 100.7) 49.8 (41.8, 55.3) 46 (34, 61) 0.6 (0.3, 0.9) 17 (10, 26) 3 (1, 6)  
Tulare 74.4 (60.7, 85.5) 4.1 96.8 (77.7, 110.1) 55.6 (45.2, 65.4) 80 (59, 100) 0.7 (0.3, 1.2) 31 (20, 48)   
Ventura 65.0 (56.6, 73.4) 2.5 77.7 (65.6, 91.7) 54.6 (46.8, 62.1) 59 (43, 84) 0.6 (0.3, 1.3) 21 (13, 38) 3 (0, 11)  
Overall 70.4 (58.1, 83.1) 3.1 88.7 (70.2, 106.7) 57.6 (47.1, 67.2) 63 (37, 95) 1.0 (0.3, 1.9) 37 (15, 65) 5 (1, 12) 15.3 (7.5, 26.7) 

Abbreviations: CO, carbon monoxide; NO2, nitrogen dioxide; O3, ozone; PM2.5, particulate matter less than 2.5 μm in aerodynamic diameter; SO2, sulfur dioxide.

a

1°F = 0.56°C.

The cumulative average weekly lag (lag 06) provided the best model fit for mean, minimum, and maximum apparent temperatures (Figure 2), although all temperature metrics were associated with significant elevated risks. Mean and minimum apparent temperature had higher per-unit effect estimates than maximum apparent temperature, although ranges for these metrics varied. Scaling by the difference between the 5th and 95th percentiles in each metric yielded virtually identical excess risks. Because mean apparent temperature consistently supplied the best overall model fit, this temperature metric with lag06 was used in subsequent analyses.

Figure 2.

Map of California showing county-level estimates for the estimated percent change associated with a 10°F (5.6°C) increase in weekly average apparent temperature (lag06) and preterm birth, May–September, 1999–2006. For each county, 95% confidence intervals are given in parentheses.

Figure 2.

Map of California showing county-level estimates for the estimated percent change associated with a 10°F (5.6°C) increase in weekly average apparent temperature (lag06) and preterm birth, May–September, 1999–2006. For each county, 95% confidence intervals are given in parentheses.

Figure 3 is a map of California with county-level estimates for the 16 counties studied. As shown, most counties had elevated risks, except for Merced and Solano, the 2 smallest, which had nonsignificant negative risks. In meta-analyses, northern California (12.1%, 95% confidence interval (CI): 7.5, 16.8) had a higher association than southern California (6.9%, 95% CI: 3.7, 10.1) (P = 0.07 for the difference), whereas coastal (7.9%, 95% CI: 4.1, 11.9) and inland (8.7%, 95% CI: 5.5, 13.1) areas had similar associations (P = 0.72). When considering associations by month of exposure, we found that the months May–July were associated with elevated risks, whereas August and September were not. We did not find a significant differential impact by year.

Figure 3.

Estimated percent change associated with a 10°F (5.6°C) increase in mean, minimum, and maximum apparent temperature and preterm birth by lag time for 16 counties in California, May–September, 1999–2006. CI, confidence interval.

Figure 3.

Estimated percent change associated with a 10°F (5.6°C) increase in mean, minimum, and maximum apparent temperature and preterm birth by lag time for 16 counties in California, May–September, 1999–2006. CI, confidence interval.

As demonstrated in Figure 4, all racial/ethnic groups analyzed had significantly elevated risks for weekly apparent temperature and preterm delivery. Black mothers had the highest risk, followed by Asians, Hispanics, and whites. Differences between the estimates for any racial/ethnic group were not significant (P = 0.17 between blacks and whites).

Figure 4.

Estimated percent change associated with a 10°F (5.6°C) increase in weekly average apparent temperature and preterm birth by maternal race/ethnicity for 16 counties in California, May–September, 1999–2006. CI, confidence interval.

Figure 4.

Estimated percent change associated with a 10°F (5.6°C) increase in weekly average apparent temperature and preterm birth by maternal race/ethnicity for 16 counties in California, May–September, 1999–2006. CI, confidence interval.

Elevated risk was found across all maternal age group categories (Figure 5). A slightly higher risk was found for the younger age groups, specifically for mothers less than 20 years of age (P = 0.04 compared with mothers at least 35 years of age). A significant risk was also found for mothers between the ages of 25 and 34 years, whereas mothers who were at least 35 years of age had a slightly elevated but nonsignificant risk.

Figure 5.

Estimated percent change associated with a 10°F (5.6°C) increase in weekly average apparent temperature and preterm birth by maternal age (years) for 16 counties in California, May–September, 1999–2006. CI, confidence interval.

Figure 5.

Estimated percent change associated with a 10°F (5.6°C) increase in weekly average apparent temperature and preterm birth by maternal age (years) for 16 counties in California, May–September, 1999–2006. CI, confidence interval.

No significant differences were found by either infant's sex or maternal education. Female infants (10.1%, 95% CI: 5.9, 14.4 per 10°F increase in apparent temperature) had a slightly more elevated risk, but it was not significantly different from that for male infants (7.8%, 95% CI: 4.2, 11.6; P = 0.44). Mothers who completed at most a high school education yielded slightly higher estimates (9.7%, 95% CI: 6.2, 13.2) than mothers who had some college education (7.0%, 95% CI: 2.7, 11.6), but this difference was not significant (P = 0.36).

For a sufficient number of cases, nitrogen dioxide and ozone data were available for all 16 counties, as in the primary analysis. However, some counties did not have enough cases with carbon monoxide, sulfur dioxide, or PM2.5 data for analysis, and these counties were excluded (refer to Table 2 for counties with missing data). No significant confounding was found for any pollutant examined; the estimates for apparent temperature did not decrease with any air pollutant in the model (not shown). In the same models, no pollutant demonstrated an independently significant positive association with preterm delivery. There was also no consistent significant interaction between apparent temperature and any air pollutant.

In the analysis of preterm delivery by gestational week, the strongest impacts were observed for preterm births from 34 to 36 gestational weeks (13.5%, 95% CI: 10.2, 16.9). For earlier preterm births, however, the association was significantly negative (−5.2%, 95% CI: −9.9, −0.4). When the data were limited to only full-term births, a relatively small but significant association (1.8%, 95% CI: 0.8, 2.7) was found for the same-day lag of mean daily apparent temperature. No association was found for full-term births when weekly apparent temperature was considered (0.0%, 95% CI: −0.9, 0.9). A significantly negative association was found for both gestational weeks 37 and 38.

DISCUSSION

A significant positive association was found between apparent temperature and preterm delivery during the warm season in California. Mean, maximum, and minimum apparent temperatures all had significantly elevated associations for lag days up to 1 week. When we considered vulnerable subgroups, all studied showed increased risks regardless of maternal age, racial/ethnic group, level of educational attainment, or infant's sex. A stronger association observed earlier in the warm season could be explained by different conception patterns. Temperatures earlier in the season may also have greater impacts than those later in the season because of acclimatization.

No significant difference by infant's sex or maternal education was observed. Although a distinction by infant's sex was not expected, maternal education was used as a proxy for socioeconomic status. Maternal race/ethnicity and maternal age may have served as better markers of socioeconomic status, since younger mothers had greater risks. Younger mothers are generally of lower socioeconomic status, have completed less education, and exhibit more risky behavioral patterns (34), particularly mothers of preterm infants (35–37). Our findings were significant across all racial/ethnic groups, with the strongest associations for black and Asian mothers. Previous investigators have also indicated disparities in birth outcomes by race/ethnicity for black (38–40), Hispanic (39), and Native American (41) infants.

Among preterm births, the effect estimates were all significantly elevated for up to 6 days, with a weekly average of apparent temperature being the best predictor. The same magnitude of association and statistical significance was not found for full-term births at 37–44 gestational weeks. The estimates tapered off after the first day for the full-term births and became significantly negative after a few days, suggestive of a harvesting pattern not observed in the preterm delivery analysis. Harvesting refers to the phenomenon by which an exposure merely accelerates the occurrence of an outcome by a few days. In other analyses, we found that preterm births occurring later in the pregnancy had greater associations with apparent temperature than those occurring earlier. However, there may have been some misclassification of gestational age, since it was based on the date of the last menstrual period, which may be unknown for some women. A broader definition for preterm delivery such as 20–36 weeks, as used in our main analyses, is subject to less error.

Much of the prior evidence for preterm delivery from environmental exposures has been based on air pollutants (11–16). In this study, temperature was found to be associated with preterm delivery independent of air pollutants. Because this analysis was restricted to the warm season and short-term associations of up to 1 week, any association with air pollutants may have been masked by ambient heat exposure. Air pollutants, on the other hand, may have had a larger impact if longer-term exposure windows, such as trimester averages, were considered (14, 42).

To our knowledge, only 2 previous investigations have focused on preterm labor or delivery and ambient temperature (18, 19). In a large time-series study recently conducted in London, United Kingdom, no associations were reported between preterm delivery and short-term meteorologic or air pollution variables (18). In a preliminary analysis in New York examining heat-humidity index, an association with preterm labor was found, but not with preterm delivery (19). With only a few cases in the study, it is possible that potentially small associations may have been missed. Investigators from both studies found that a weekly lag, as used in our study, was most relevant in terms of preterm delivery. Single dates could result in misclassification bias since labor may last more than 1 day, or a mother may not be admitted until a day following the high heat-humidity index (19). These 2 previous studies did not analyze specific subgroups. Health effects from ambient temperature are based primarily on studies of mortality (43, 44).

Since preterm delivery has multiple etiologies, a clear biologic mechanism or cause is unknown. One possible explanation may be increased dehydration with heat exposure, which could decrease uterine blood flow and increase pituitary secretion of antidiuretic hormone and oxytocin to induce labor (45). Pregnant women may not be able to thermoregulate efficiently. When body temperatures rise, the body generally shifts blood flow from the vital organs to the skin's surface in an effort to cool down (46). Thus, thermoregulation may be inadequate when too much blood is diverted from the vital organs of the mother and the developing fetus (46). Increased blood viscosity, elevated cholesterol levels associated with higher temperatures, and a higher sweating threshold have also been reported in susceptible subgroups (47).

There are several limitations to this study. We relied on ecologically monitored data rather than on individual monitoring to assign exposure. By limiting the study to mothers whose residential zip codes were within 10 km of a meteorologic monitor, we minimized exposure misclassification, but it still could have been a problem, especially in coastal areas. We had to rely on data provided on the birth certificates, and we did not have information on several characteristics, such as individual income level, to adjust for socioeconomic status more efficiently than relying on maternal age and racial/ethnic group. We also did not have information on occupational address or air conditioning status, so we could not account for these variables in our analysis. For most cases, utilizing the zip code for the residential address reported on the birth certificate should be sufficient to characterize exposure for the week prior to birth.

This is the first known large-scale, population-based epidemiologic study of ambient temperature and preterm delivery in the United States covering an extensive time period. It is also the first study known to find significant positive associations and to examine potentially vulnerable subgroups. All demographic subgroups examined had elevated risks, although some differences were found between groups. This study focused on high ambient temperature, and the association during a heat wave is likely to be substantially greater, as was found in a recent mortality study (48). The Intergovernmental Panel on Climate Change states that hot days and nights, as well as heat waves, have become more frequent in recent years, and their duration and intensity are likely to increase in the future (49). Given the significant associations for apparent temperature and preterm delivery found in this study, more large-scale studies of temperature and preterm delivery are warranted. Biologic acclimatization and susceptible subgroups should be considered.

Abbreviations

    Abbreviations
  • CI

    confidence interval

  • PM2.5

    particulate matter less than 2.5 μm in aerodynamic diameter

Author affiliation: California Office of Environmental Health Hazard Assessment, Air Pollution Epidemiology Section, Oakland, California (Rupa Basu, Brian Malig, Bart Ostro).

The opinions expressed in this article are solely those of the authors and do not represent the policy or position of the state of California or the California Environmental Protection Agency.

Conflict of interest: none declared.

References

1.
Martin
JA
Hamilton
BE
Sutton
PD
, et al.  . 
Births: final data for 2006
Natl Vital Stat Rep.
 , 
2009
, vol. 
57
 
7
(pg. 
1
-
102
)
2.
Mathews
TJ
MacDorman
MF
Infant mortality statistics from the 2004 period linked birth/infant death data set
Natl Vital Stat Rep.
 , 
2007
, vol. 
55
 
14
(pg. 
1
-
32
)
3.
Escobar
GJ
Clark
RH
Greene
JD
Short-term outcomes of infants born at 35 and 36 weeks gestation: we need to ask more questions
Semin Perinatol.
 , 
2006
, vol. 
30
 
1
(pg. 
28
-
33
)
4.
Mangham
LJ
Petrou
S
Doyle
LW
, et al.  . 
The cost of preterm birth throughout childhood in England and Wales
Pediatrics.
 , 
2009
, vol. 
123
 
2
(pg. 
e312
-
e327
)
5.
Cuevas
KD
Silver
DR
Brooten
D
, et al.  . 
The cost of prematurity: hospital charges at birth and frequency of rehospitalizations and acute care visits over the first year of life: a comparison by gestational age and birth weight
Am J Nurs.
 , 
2005
, vol. 
105
 
7
(pg. 
56
-
64
)
6.
Petrou
S
The economic consequences of preterm birth during the first 10 years of life
BJOG.
 , 
2005
, vol. 
112
 
suppl 1
(pg. 
S10
-
S15
)
7.
Petrou
S
Mehta
Z
Hockley
C
, et al.  . 
The impact of preterm birth on hospital inpatient admissions and costs during the first 5 years of life
Pediatrics.
 , 
2003
, vol. 
112
 
6 pt 1
(pg. 
1290
-
1297
)
8.
Marlow
N
Wolke
D
Bracewell
MA
, et al.  . 
Neurologic and developmental disability at six years of age after extremely preterm birth
N Engl J Med.
 , 
2005
, vol. 
352
 
1
(pg. 
9
-
19
)
9.
Anderson
P
Doyle
LW
Neurobehavioral outcomes of school-age children born extremely low birth weight or very preterm in the 1990s
JAMA.
 , 
2003
, vol. 
289
 
24
(pg. 
3264
-
3272
)
10.
Mattison
DR
Wilson
S
Coussens
C
, et al.  . 
The Role of Environmental Hazards in Premature Birth
 , 
2003
Washington, DC
National Academies Press
11.
Ritz
B
Wilhelm
M
Hoggatt
KJ
, et al.  . 
Ambient air pollution and preterm birth in the environment and pregnancy outcomes study at the University of California, Los Angeles
Am J Epidemiol.
 , 
2007
, vol. 
166
 
9
(pg. 
1045
-
1052
)
12.
Wilhelm
M
Ritz
B
Local variations in CO and particulate air pollution and adverse birth outcomes in Los Angeles County, California, USA
Environ Health Perspect.
 , 
2005
, vol. 
113
 
9
(pg. 
1212
-
1221
)
13.
Sagiv
SK
Mendola
P
Loomis
D
, et al.  . 
A time-series analysis of air pollution and preterm birth in Pennsylvania, 1997–2001
Environ Health Perspect.
 , 
2005
, vol. 
113
 
5
(pg. 
602
-
606
)
14.
Huynh
M
Woodruff
TJ
Parker
JD
, et al.  . 
Relationships between air pollution and preterm birth in California
Paediatr Perinat Epidemiol.
 , 
2006
, vol. 
20
 
6
(pg. 
454
-
461
)
15.
Brauer
M
Lencar
C
Tamburic
L
, et al.  . 
A cohort study of traffic-related air pollution impacts on birth outcomes
Environ Health Perspect.
 , 
2008
, vol. 
116
 
5
(pg. 
680
-
686
)
16.
Studnicka
M
Hackl
E
Pischinger
J
, et al.  . 
Traffic-related NO2 and the prevalence of asthma and respiratory symptoms in seven year olds
Eur Respir J.
 , 
1997
, vol. 
10
 
10
(pg. 
2275
-
2278
)
17.
Darrow
LA
Strickland
MJ
Klein
M
, et al.  . 
Seasonality of birth and implications for temporal studies of preterm birth
Epidemiology.
 , 
2009
, vol. 
20
 
5
(pg. 
699
-
706
)
18.
Lee
SJ
Hajat
S
Steer
PJ
, et al.  . 
A time-series analysis of any short-term effects of meteorological and air pollution factors on preterm births in London, UK
Environ Res.
 , 
2008
, vol. 
106
 
2
(pg. 
185
-
194
)
19.
Lajinian
S
Hudson
S
Applewhite
L
, et al.  . 
An association between the heat-humidity index and preterm labor and delivery: a preliminary analysis
Am J Public Health.
 , 
1997
, vol. 
87
 
7
(pg. 
1205
-
1207
)
20.
California Irrigation Management Information System (CIMIS)
 
Davis, CA: Office of Water Use Efficiency, Department of Water Resources; 2009. (http://www.cimis.water.ca.gov/cimis/welcome.jsp). (Accessed January 29, 2009)
21.
US Environmental Protection Agency
Air Quality System Data Mart [Internet database]
 , 
2009
Washington, DC
Environmental Protection Agency
 
(http://www.epa.gov/ttn/airs/aqsdatamart/). (Accessed August 18, 2009)
22.
Beyer
HL
Hawth's Analysis Tools for ArcGIS
 , vol. 
2004
  
(http://www.spatialecology.com/htools). (Accessed September 13, 2010)
23.
Wilhelm
M
Ritz
B
Residential proximity to traffic and adverse birth outcomes in Los Angeles County, California, 1994–1996
Environ Health Perspect.
 , 
2003
, vol. 
111
 
2
(pg. 
207
-
216
)
24.
Salam
MT
Millstein
J
Li
YF
, et al.  . 
Birth outcomes and prenatal exposure to ozone, carbon monoxide, and particulate matter: results from the Children's Health Study
Environ Health Perspect.
 , 
2005
, vol. 
113
 
11
(pg. 
1638
-
1644
)
25.
California Air Resources Board
Air Quality Data Section, Planning and Technical Support Division. California Air Quality Data [Internet database]
 , 
2006
Sacramento, CA
California Air Resources Board
 
(http://www.arb.ca.gov/aqd/aqdcd/aqdcddld.htm). (Accessed August 18, 2008)
26.
California Department of Health Services
Birth Statistical Master Files, 1999–2006
 
Sacramento, CA
California Department of Health Services
27.
Alexander
GR
Himes
JH
Kaufman
RB
, et al.  . 
A United States national reference for fetal growth
Obstet Gynecol.
 , 
1996
, vol. 
87
 
2
(pg. 
163
-
168
)
28.
Levy
D
Lumley
T
Sheppard
L
, et al.  . 
Referent selection in case-crossover analyses of acute health effects of air pollution
Epidemiology.
 , 
2001
, vol. 
12
 
2
(pg. 
186
-
192
)
29.
Lumley
T
Levy
D
Bias in the case-crossover design: implications for studies of air pollution
Environmetrics.
 , 
2000
, vol. 
11
 (pg. 
689
-
704
)
30.
DerSimonian
R
Laird
N
Meta-analysis in clinical trials
Control Clin Trials.
 , 
1986
, vol. 
7
 
3
(pg. 
177
-
188
)
31.
SAS Statistical Software [computer program]. Version 9.1
 , 
2008
Cary, NC
SAS Institute, Inc
32.
R: A Language and Environment for Statistical Computing, version 2.6.2
 , 
2008
Vienna, Austria
R Foundation for Statistical Computing, R Development Core Team
 
(http://www.R-project.org). (Accessed September 13, 2010)
33.
Wilson
WE
Mar
TF
Koenig
JQ
Influence of exposure error and effect modification by socioeconomic status on the association of acute cardiovascular mortality with particulate matter in Phoenix
J Expo Sci Environ Epidemiol.
 , 
2007
, vol. 
17
 
suppl 2
(pg. 
S11
-
S19
)
34.
McDonald
JA
Suellentrop
K
Paulozzi
LJ
, et al.  . 
Reproductive health of the rapidly growing Hispanic population: data from the Pregnancy Risk Assessment Monitoring System, 2002
Matern Child Health J.
 , 
2008
, vol. 
12
 
3
(pg. 
342
-
356
)
35.
Cho
J
Holditch-Davis
D
Belyea
M
Gender, ethnicity, and the interactions of prematurely born children and their mothers
J Pediatr Nurs.
 , 
2004
, vol. 
19
 
3
(pg. 
163
-
175
)
36.
Holditch-Davis
D
Miles
MS
Weaver
MA
, et al.  . 
Patterns of distress in African-American mothers of preterm infants
J Dev Behav Pediatr.
 , 
2009
, vol. 
30
 
3
(pg. 
193
-
205
)
37.
Holditch-Davis
D
Schwartz
T
Black
B
, et al.  . 
Correlates of mother-premature infant interactions
Res Nurs Health.
 , 
2007
, vol. 
30
 
3
(pg. 
333
-
346
)
38.
Whitley
R
Kirmayer
LJ
Perceived stigmatisation of young mothers: an exploratory study of psychological and social experience
Soc Sci Med.
 , 
2008
, vol. 
66
 
2
(pg. 
339
-
348
)
39.
Shone
LP
Dick
AW
Brach
C
, et al.  . 
The role of race and ethnicity in the State Children's Health Insurance Program (SCHIP) in four states: are there baseline disparities, and what do they mean for SCHIP? Pediatrics
2003
, vol. 
112
 
6
pg. 
e521
  
(doi:10.1542/peds.112.6.SE1.521)
40.
Buescher
PA
Mittal
M
Racial disparities in birth outcomes increase with maternal age: recent data from North Carolina
N C Med J.
 , 
2006
, vol. 
67
 
1
(pg. 
16
-
20
)
41.
Sparks
PJ
Do biological, sociodemographic, and behavioral characteristics explain racial/ethnic disparities in preterm births?
Soc Sci Med.
 , 
2009
, vol. 
68
 
9
(pg. 
1667
-
1675
)
42.
Parker
JD
Woodruff
TJ
Basu
R
, et al.  . 
Air pollution and birth weight among term infants in California
Pediatrics.
 , 
2005
, vol. 
115
 
1
(pg. 
121
-
128
)
43.
Basu
R
High ambient temperature and mortality: a review of epidemiologic studies from 2001 to 2008
Environ Health.
 , 
2009
, vol. 
8
 pg. 
40
  
(doi: 10.1186/1476-069X-8-40)
44.
Basu
R
Samet
JM
Relation between elevated ambient temperature and mortality: a review of the epidemiologic evidence
Epidemiol Rev.
 , 
2002
, vol. 
24
 
2
(pg. 
190
-
202
)
45.
Stan
C
Boulvain
M
Hirsbrunner-Amagbaly
P
, et al.  . 
Hydration for treatment of preterm labour
Cochrane Database Syst Rev.
 , 
2002
2
pg. 
CD003096
  
(doi: 10.1002/14651858.CD003096)
46.
Bouchama
A
Knochel
JP
Heat stroke
N Engl J Med.
 , 
2002
, vol. 
346
 
25
(pg. 
1978
-
1988
)
47.
Astrand
PO
Rodahl
K
Dahl
HA
, et al.  . 
Textbook of Work Physiology: Physiological Bases of Exercise
 , 
2003
4th ed.
Windsor, Canada
Human Kinetics
48.
Ostro
BD
Roth
LA
Green
RS
, et al.  . 
Estimating the mortality effect of the July 2006 California heat wave
Environ Res.
 , 
2009
, vol. 
109
 
5
(pg. 
614
-
619
)
49.
Intergovernmental Panel on Climate Change
Climate Change 2007: Synthesis Report
  
(Intergovernmental Panel on Climate Change Fourth Assessment Report). Geneva, Switzerland: World Meteorological Organization; 2007. (http://www.ipcc.ch/pdf/assessment-report/ar4/syr/ar4_syr.pdf). (Accessed August 11, 2009)

Author notes

Editor's note: An invited commentary on this article appears on page 1118, and the authors’ response is published on page 1121.