- Open Open Peer Review
Associations between ambient air temperature, low birth weight and small for gestational age in term neonates in southern Israel
Environmental Healthvolume 17, Article number: 76 (2018)
The increase in ambient temperatures (Ta) and emissions of greenhouse gases over the last century has focused attention on the effects of ambient temperatures on health outcomes. We aimed to investigate the association between Ta and the clinical measures of term low birth weight (tLBW) and small for gestational age (SGA) in singleton term infants using a decade of regional hospital data in southern Israel.
We linked all births in Soroka University Medical Center in the southern district of Israel insured by Clalit Health Services with pregnancy Ta estimated by our novel hybrid spatio-temporally resolved prediction model. Logistic regression generalized additive models and general linear models were used, with either tLBW or SGA as the dependent variable, modeling entire pregnancy and trimester-specific Ta adjusting for seasonality, time trend, particulate matter, maternal age, gravidity, parity, ethnicity, sex, poverty index and population density.
The study population included 56,141 singleton term newborns, with 1716 (3.1%) cases of tLBW and 8634 (15.4%) cases of SGA. The average and the median Ta across the entire pregnancy were 19.9 (SD: 1.77, range: 14.6–24.9) degrees centigrade. The lowest Ta quartile (Ta = < 18.5) was associated with higher risk of tLBW (odds ratio = 1.33, 95%CI 1.11–1.58) while the highest Ta quartile (Ta > =21.3) was not significantly associated with tLBW (odds ratio = 1.17, 95%CI 0.99–1.38), in comparison to the two intermediate quartiles. When analyzing SGA as the dependent variable, the lowest Ta quartile was associated with significantly higher risk of SGA (odds ratio = 1.18, 95%CI 1.09–1.29) while the highest quartile was associated with significantly lower risk of SGA (odds ratio = 0.91, 95%CI 0.84–0.99) in comparison to the two intermediate quartiles.
Our findings suggest that lower pregnancy Ta may increase the risk of tLBW and SGA, and higher pregnancy Ta may decrease the risk of SGA in singleton term infants in southern Israel.
The rise in temperatures over the last century and emissions of greenhouse gases has focused attention on the effects of increasing heat . The effects of temperatures and extreme events on health outcomes have been extensively studied [2,3,4,5,6].
Many studies have shown that heat is related to various health outcomes including an increase in total and cause specific mortality [4, 7,8,9,10,11,12] as well hospital admissions, emergency room visits and ambulance calls [13, 14]. The vulnerability of human populations to extreme temperatures is a function of their sensitivity to the exposure, of the character, magnitude and rate of the climate extreme, and of the adaptation measures and actions in place. The association between temperature and mortality has been identified and described as a non-linear U-, J- or V-shaped function, with the lowest mortality rates recorded at moderate temperatures, rising progressively as temperatures increase or decrease [7, 10, 13]. Finally, these studies have shown that the effects of temperature vary geographically [15, 16] and depend on local climatic conditions and population characteristics such as demography, socioeconomic conditions and health status .
Among many susceptible groups, pregnant women, fetuses and children are particularly vulnerable to the physiologic stress of extreme temperatures and sudden temperature variation [17, 18]. Yet only in recent years, studies have started to examine the association between maternal exposure to extreme ambient air temperatures (Ta) and adverse pregnancy outcomes. Basu and colleagues  examined apparent temperature (a combination of temperature and relative humidity) with full-term low birth weight (tLBW) comparing 43,629 tLBW infants and 2,032,601 normal weight infants in California from 1999 to 2013. They presented 13% (95% confidence interval (CI): 4.1–22.7%) increased risk of tLBW per 10 °F with entire pregnancy mean apparent temperature above 55 °F (12.8 C).
Li and colleagues , examined the effects of ambient temperature at three trimesters of pregnancy on preterm birth and stillbirth, and evaluate the effect changes during 1994–2013. Their study presented results showing that both low and high temperatures at the 3rd trimester of pregnancy significantly increased the risk of preterm birth, with similar hazard ratios (HR) for low (HR = 1.21,95% CI: 1.16–1.27) and high (HR = 1.2, 95% CI: 1.16–1.26) temperatures in comparison with minimum-prevalence temperatures for each trimester. A number of studies have looked at seasonal patterns in pregnancy outcomes, and recently some studies have focused on the relationship between ambient temperature exposure at different time-windows during pregnancy and outcomes including preterm birth and birth weight [18, 21,22,23].
Extremes of both high [24,25,26,27,28,29] and low  temperatures, mainly estimated during the last month or days of pregnancy, have been associated with an increased risk of preterm birth. Similarly, lower birth weight has been associated with both higher [31,32,33,34] and lower  ambient temperature during specific trimesters of pregnancy. In contrast, a large study of more than 480,000 births found no association between preterm birth and a variety of factors including temperature, humidity and barometric pressure .
These conflicting results published to date on the relationship of Ta with birth outcomes may be due to several reasons. There could be differences between local climatic and geographical conditions of the study areas as well as different levels of acclimatization among the pregnant women across the populations . One established reason could be attributed to differences in temperature exposure assessment [11, 32]. In addition, there is limited evidence for populations living in suburban and rural areas due to the limited temperature data available. Most of these studies have been carried out in urban areas primarily due to data availability in terms of exposure and the number of outcome counts, large enough to ensure statistical power. This may result in lower variability of the exposure and possibly bias in the estimates related to the exposure considered due to uncertainty in assessment of individual exposure levels .
While preterm birth and multiple births pregnancy are the major reasons for low birth weight (LBW), approximately 90% of newborns are singleton term infants. tLBW and small for gestational age (SGA) are two main clinical measure that define fetal growth restriction in term infants, a phenomenon that predicts various complications in early childhood and beyond, including perinatal mortality, impaired immune system, respiratory morbidity and neurobehavioral disorders [38, 39]. In this study we make use of novel Ta estimation from our spatio-temporally resolved satellite-based hybrid models in Israel  to study the association between Ta and clinical measures of birth weight in term infants - tLBW and SGA, using regional hospital data from 2004 to 2013. Using these models, we can overcome the common limitation of insufficiently detailed climate data – in space and time – for estimating robust exposure to temperature and examine its associations with tLBW and SGA. In addition, we study this question in unique geo-climatic conditions that vary widely across southern of Israel, including 4 known climate zones (Mediterranean climate, semi-arid climate, arid climate and extreme arid climate) with maximum temperature which can exceed 45 °C in extreme days. These climate conditions and high Ta are important for insights into the effects of temperature on birth outcomes in a changing climate with extreme weather events.
Study domain and population
The research area is the Southern district of Israel (see Fig. 1). The research area has unique geo-climatic conditions that vary widely across Southern of Israel, including 4 known climate zones (Mediterranean climate, semi-arid climate, arid climate and extreme arid climate) defined by the Israel Meteorological Service  and based on the Kopen classification . All these zones are characterized by long, hot, rainless summers and relatively short, cool, rainy winters. The rain period in the region takes place across the months of October–May . The climate conditions within each zone of the study region vary locally by altitude, latitude, and the proximity to the Mediterranean. Summers are very humid along the Mediterranean coast but dry in the Rift Valley and the Negev Desert. The terrain also varies widely from the deepest point in the Jordan Rift Valley (− 430 m below sea level).
The study population included all singleton term births in Soroka University Medical Center (SUMC) during October 2004 – December 2013, residing in the southern district of Israel and insured by Clalit Health Services, the largest health care provider in Israel, covering approximately 70% of a population of 730,000 residents in the area [44, 45]. The rest of the population in the area is covered similarly by one of other three healthcare organizations acting under the national health act, and each resident is free to choose his/her provider among these four providers and move from one to another, with very few limitations. SUMC is a tertiary hospital in southern Israel, which receives most birth referrals in its area (Fig. 1).
Maternal and birth data were routinely collected in the Admission-Transfer-Discharge SUMC computerized database. The study population contained 56,184 singleton term births, from which we excluded subjects with missing Ta (N = 23) or birth weight (N = 9) and subjects with very extreme Ta (> 3 standard deviations, N = 11), leaving 56,141 subjects for analysis. tLBW was defined as birth weight < 2500 g, and SGA was defined as birth weight < 10th percentile (in Israel), according to newborn sex and gestational age at delivery .
Air temperature: Predicted 24 h mean Air temperature (Ta) exposure data were generated by a novel hybrid spatio-temporally resolved prediction model . Briefly, in these prediction model we used mixed effects models to first calibrate earth observation satellite surface temperature (Ts) data to daily mean Ta measurements from ground monitors, regressing Ta measurements against day-specific random intercepts, fixed and random Ts slopes and several spatial and temporal predictors (normalized difference vegetation index, percent urban and elevation). Then to capture the ability of neighboring cells to fill in the cells with missing Ts values, we regressed the Ta predicted from the first stage model against the mean of the Ta measurements on that day from interpolated inverse distance weighted models, discretely for each grid cell. Our model performance was excellent for both days with available Ts and days without Ts observations (cross-validated R2 results: 0.966–0.986). To estimate 24 h mean Ta exposure at each mother’s residence during pregnancy, we spatially joined each mother’s residence to the corresponding Ta grid cell centroid (1 km*1 km resolution) into which it fell (Fig. 1). Ta Exposure history was calculated by averaging mean daily Ta over the entire pregnancy (calculated by exact birth date and gestational age at birth for each pregnancy) and for each trimester. In general, studies in environmental epidemiology looking at acute events (such as heat waves) and the associations with health outcomes focus on acute exposure windows. This is partly due to the difficulty of detangling these acute effects from the long-term effect. In line with such previous studies , we decided to focus on the chronic/long term exposure. It should be noted that acute heat wave events are still included in chronic exposure moving averages and thus accounted for in the model.
Air pollution: Daily average particulate matter < 2.5 μm (PM2.5) concentrations were assessed using a hybrid satellite based model incorporating daily satellite remote sensing data at 1 × 1 km spatial resolution . Briefly, we used an algorithm developed by NASA Multi-Angle Implementation to Atmospheric Correction, which provides aerosol optical depth (AOD) data in a high resolution. Using mixed models, we regressed daily PM2.5 mass concentration from the Ministry of Environmental Protection monitoring sites against AOD, traditional land use regression temporal and spatial predictors to generate daily high resolution PM2.5 estimations. More in depth details can be found in Kloog 2015 .
For initial analyses, we used logistic regression generalized additive models with either tLBW or SGA as the dependent variable, modeling entire pregnancy Ta using a thin plate regression spline (allowing for non-linear associations). The analysis was adjusted for seasonality and time trend using categorical variables for calendar month and year, PM2.5, maternal age at birth (continuous, squared), gravidity, parity, ethnicity (Jewish / non-Jewish), newborn sex, census-level poverty index (a categorical variable with 20 levels, representing socio-economic level) and population density (based on the population residing in the same small-statistical-area). The two latter variables were defined and calculated by the Israel Central Bureau of Statistics .
Further analyses used logistic regression generalized linear models with either tLBW or SGA as the dependent variable, modeling entire pregnancy Ta using a categorical variable and each of the two extreme Ta quartiles with the two middle quartiles, adjusting for the same covariates as described above. Possible effect modification by ethnicity or urbanity were tested by adding interaction terms to these models. Trimester-specific Ta exposures were modeled similarly, in a mutually adjusted model that included the exposure in each trimester as a separate variable, in order to avoid bias from auto-correlation among Ta trimester values .
In a sensitivity analysis, we fit a model that included only those newborns for whom addresses were geocoded at the house or street level (N = 18,412). Statistical analyses were performed in R statistical software version 3.3.2 (R Foundation for Statistical Computing, Vienna, Austria) with general additive models using package mgcv . A priory level of significance was 0.05 (two-sided).
The study population is described in Table 1. Of the 56,141 singleton term births included in our final dataset for analysis, 51% were males and 40% were Jewish. The mean birth weight was 3272 g (SD: 430 g), with 1716 (3.1%) cases of tLBW and 8634 (15.4%) cases of SGA. The average and the median daily Ta across the entire pregnancy were 19.9 (SD: 1.77, range: 14.6–24.9) degrees centigrade. Spearman correlation among trimester-specific Ta were: r = 0.048, r = − 0.9 and r = − 0.013 for pairs of trimesters 1–2, 2–3 and 1–3 respectively.
Generalized additive logistic regression model revealed a significant (p < 0.001), nonlinear (estimated degrees of freedom = 3.2), non-monotonic association between entire pregnancy average Ta and tLBW, with both low and high temperatures associated with higher risk in comparison to intermediate temperatures, although only on low Ta was the association significant (Fig. 2a). A similar model for SGA as the dependent variable resulted in a decreasing, almost linear (estimated degrees of freedom = 2.2) curve between Ta and SGA (p < 0.001) (Fig. 2b).
We further categorized the entire pregnancy average Ta into quartiles and compared the lowest Ta quartile (Ta = < 18.5) and the highest Ta quartile (Ta > =21.3) to the two intermediate quartiles (18.5 < Ta < 21.3) using logistic regression models adjusted for the same covariates as above. The lowest Ta quartile was associated with higher risk of tLBW (odds ratio [OR] = 1.33, 95%CI 1.11–1.58) while the highest Ta quartile was not significantly associated with tLBW (OR = 1.17, 95%CI 0.99–1.38), in comparison to the two intermediate quartiles. When analyzing SGA as the dependent variable, the lowest Ta quartile was associated with significantly higher risk of SGA (OR = 1.18, 95%CI 1.09–1.29) while the highest quartile was associated with significantly lower risk of SGA (OR = 0.91, 95%CI 0.84–0.99) in comparison to the two intermediate quartiles (Table 2).
When examining associations among tLBW, SGA and trimester-specific Ta exposures, we found decreased risk of SGA with high Ta quartile during the 1st trimester (OR = 0.82, 95%CI 0.73–0.92), and increased risk for SGA with low Ta quartile during the 3rd trimester (OR = 1.17, 95%CI 1.05–1.31), while all other associations were not statistically significant (Table 3).
In a sensitivity analysis for the entire pregnancy average Ta, including only those newborns whose addresses were geocoded at the building or the street level (N = 18,412), the lowest Ta quartile showed a stronger association (OR = 1.48, 95%CI 1.06–2.06), and the highest Ta quartile was still not associated (OR = 1.01, 95%CI 0.74–1.39) with tLBW, in comparison with the two intermediate quartiles. When examining Ta quartiles association with SGA in this sub-population, the associations were similar to the primary analyses but not statistically significant (OR = 1.09, 95%CI 0.94–1.28 and OR = 0.93, 95%CI 0.80–1.08 for the lowest and the highest Ta quartiles, respectively). When adding an interaction term to the main models to test for possible effect modification by ethnicity, we did not find any significant interaction between Ta and ethnicity (p > 0.45 and p > 0.2 for LBW and SGA, respectively). Similarly, no significant interaction with urbanity was found (p > 0.2 for all models).
In the presented study, we examined the association of Ta with tLBW and SGA in a study of term singleton births in Southern Israel between 2004 and 2013. Using our novel exposure model, we were able to assign exposure to almost all subjects with substantially reduced exposure misclassification compared to using monitored data. We found that lower Ta was associated with increased risk of tLBW and SGA, and higher Ta was associated with decreased risk of SGA in this population. When examining trimester-specific Ta, the association with higher Ta was restricted to the 1st trimester, while the association with lower Ta was restricted to the 3rd trimester. In other words, higher than usual Ta in early phases of the pregnancy seems to be protective against SGA in term infants, and lower than usual Ta in the late phases of pregnancy may be a risk factor of SGA in our study population.
Biological mechanisms that could explain a possible effect of temperature on tLBW and SGA have been rarely investigated, but several potential hypotheses exist. Some studies suggest that variation in temperature and humidity, especially at high temperature and humidity levels, can increase stress on cardiovascular function , which, in pregnancy, is already highly solicited. On the other hand, our results suggest higher risk of SGA in term infants with lower than average Ta, but not with higher Ta. The second and third trimesters are important times for fetal growth . As extreme temperature may affect uterine blood flow and placental exchange necessary for fetal growth , disruption to the mechanism needed for proper growth during these time windows would have the greatest impact.
The results of our study regarding tLBW are quite consistent with those of a recent study from California, with very similar U-shape curves describing associations between apparent temperature and risk of tLBW . Another recent large American study of participants from the consortium on safe labor, with almost 200,000 term infants, found both extreme high and extreme low Ta to be associated with risk of tLBW . This study, however, used more extreme Ta cutoffs (<5th and > 95th percentiles) than we did. Entire pregnancy Ta was not significantly associated with SGA (without limiting to term infants) in that study, and colder Ta during the 3rd trimester were associated with lower risk of SGA, in contrast to our findings. These inconsistencies among the studies may be a result of: 1) including preterm infants when analyzing SGA; 2) differences in the populations studied; and 3) the different range of Ta in Southern Israel in comparison the colder environments inspected in other studies.
Another factor that may explain differences among study results is differences in exposure misclassification, introduced by different Ta exposure assessment methods. Stronger non-differential misclassification, as expected when using lower spatial resolution in exposure assessment, is expected to bias the estimate towards the null. Since temperature can vary greatly both spatially and temporally, the use of simple models or even temperature from sparse monitor networks and relatively large hospital capture areas can introduce considerable measurement error. Another point to note is that this lack of robust spatially resolved daily Ta exposures restricts these studies to populations surrounding monitoring stations, which may not be representative of the population as a whole. On the other hand, urbanity did not modify our effect estimates significantly, suggesting that excluding rural populations may limit sample size, but not necessarily bias effect estimates.
Our study also detected decreased risk of SGA in term infants exposed to higher Ta levels during the 1st trimester. This finding is novel and requires inspection and reproduction in other cohorts as well. We are currently not aware of biological mechanisms that may explain this finding, but it may be related to other exposures, which are a result of pregnant women changing their behavior in response to high ambient temperatures, such as differences in physical activity or time spent indoors. Previous studies that examined 1st trimester exposures did not result in a similar finding, but their published results were also not adjusted for exposures during the other trimesters [19, 31]. This adjustment is important in order to avoid bias resulting from correlation of the exposure among various exposure windows, as demonstrated using simulations in the case of air pollution .
Comparison with the results of other studies of Ta and birth weight is more complicated: Most studies took a different analytical approach, examining the associations with birth weight as a continuous variable, unlike the approach we took, concentrating on clinical categories (tLBW and SGA) which are known predictors of health and disease. From the studies that did examine LBW we found 3 ecological studies [54,55,56] and one retrospective cohort  . From these, only one of the ecological studies was restricted to term births like the current study, and found no significant associations of LBW with or cold waves . However, the differences in study design (ecological vs cohort) and the limited sample size in that study (approximately 1500 births) can account for the difference in the results.
Another thing to note is that there are clear and significant changes in average Ta in Israel in recent decades. Since the 1970s, average annual temperature increased by close to 1 °C . These trends, if persist, have serious environmental and health implications in the future. However, we believe that they did not impact our study results, which are based on births during 2004–2013 and compare stronger Ta exposure contrasts (e.g. 14.6–18.5 and 18.5–21.3 for the two lowest comparison groups).
There are several limitations in the present study. First, the spatial resolution of the exposure was 1 × 1 km and exposures were assigned based on best-available residential address geolocation. We are planning to improve our hybrid models using new available satellites and machine learning methodologies to generate higher spatial resolution data reaching 100 m, which will further reduce exposure error. Other limitations with the SUMC database include the lack of some health and personal level data such as maternal weight and height, chronic hypertension, preeclampsia or gestational hypertension, maternal smoking and physical activity, which were not available. On the other hand, we do not have a reason to think that these characteristics are a common cause for both ambient Ta and tLBW or SGA, and therefore they should not be a source of concerns of residual confounding.
Our findings suggest that lower pregnancy Ta may increase the risk of tLBW and SGA, and higher pregnancy Ta may decrease the risk of SGA in singleton term infants in southern Israel. Future studies may examine associations of Ta with birth weight in preterm infants and examine possible mechanisms of effects of Ta on fetal growth.
Aerosol optical depth
Low birth weight
- PM2.5 :
Particulate matter < 2.5 μm
Small for gestational age
Soroka University Medical Center
Ambient air temperature
Term low birth weight
Crowley TJ. Causes of climate change over the past 1000 years. Science. 2000;289:270–7.
Basu R, Samet JM. Relation between elevated ambient temperature and mortality: a review of the epidemiologic evidence. Epidemiol Rev. 2002;24:190–202.
Benmarhnia T, Deguen S, Kaufman JS, Smargiassi A. Vulnerability to heat-related mortality: a systematic review, meta-analysis, and meta-regression analysis. Epidemiology. 2015;26:781–93.
Nordio F, Zanobetti A, Colicino E, Kloog I, Schwartz J. Changing Patterns of the Temperature-mortality Association By Time and Location in the Us, and Implications for Climate Change. 2015;81:80–86.
Robine J-M, Cheung SLK, Le Roy S, Van Oyen H, Griffiths C, Michel J-P, et al. Death toll exceeded 70,000 in Europe during the summer of 2003. C R Biol. 2008;331:171–8.
Xu Z, FitzGerald G, Guo Y, Jalaludin B, Tong S. Impact of heatwave on mortality under different heatwave definitions: a systematic review and meta-analysis. Environ Int. 2016;89:193–203.
Baccini M, Biggeri A, Accetta G, Kosatsky T, Katsouyanni K, Analitis A, et al. Heat effects on mortality in 15 European cities. Epidemiology. 2008;19:711–9.
Basu R. High ambient temperature and mortality: a review of epidemiologic studies from 2001 to 2008. Environ Health. 2009;8:40.
Bunker A, Wildenhain J, Vandenbergh A, Henschke N, Rocklöv J, Hajat S, et al. Effects of air temperature on climate-sensitive mortality and morbidity outcomes in the elderly; a systematic review and meta-analysis of epidemiological evidence. EBioMedicine. 2016;6:258–68.
Curriero FC, Heiner KS, Samet JM, Zeger SL, Strug L, Patz JA. Temperature and mortality in 11 cities of the eastern United States. Am J Epidemiol. 2002;155:80–7.
Shi L, Kloog I, Zanobetti A, Liu P, Schwartz J. Impacts of Temperature and Its Variability on Mortality in New England. 2015;
Von Klot S, Peters A, Aalto P, Bellander T, Berglind N, D’Ippoliti D, et al. Ambient Air Pollution Is Associated With Increased Risk of Hospital Cardiac Readmissions of Myocardial Infarction Survivors in Five European Cities, vol. 112; 2005. p. 3073.
Schwartz J, Samet JM, Patz JA. Hospital admissions for heart disease: the effects of temperature and humidity. Epidemiology. 2004;15:755–61.
Williams S, Nitschke M, Weinstein P, Pisaniello DL, Parton KA, Bi P. The impact of summer temperatures and heatwaves on mortality and morbidity in Perth, Australia 1994–2008. Environ Int. 2012;40:33–8.
Guo Y, Gasparrini A, Armstrong B, Li S, Tawatsupa B, Tobias A, et al. Global variation in the effects of ambient temperature on mortality: a systematic evaluation. Epidemiol Camb Mass. 2014;25:781.
Zanobetti A, Coull BA, Gryparis A, Kloog I, Sparrow D, Vokonas PS, et al. Associations Between Arrhythmia Episodes and Temporally and Spatially Resolved Black Carbon and Particulate Matter in Elderly Patients. 2013;71:201–207.
Balbus JM, Malina C. Identifying vulnerable subpopulations for climate change health effects in the United States. J Occup Environ Med. 2009;51:33.
Strand LB, Barnett AG, Tong S. The Influence of Season and Ambient Temperature on Birth Outcomes: a Review of the Epidemiological Literature. 2011;111:451–462.
Basu R, Rau R, Pearson D, Malig B. Temperature and Term Low Birth Weight in California. Am J Epidemiol [Internet]. 2018 [cited 2018 Jul 8]; Available from: https://academic.oup.com/aje/advance-article/doi/10.1093/aje/kwy116/5036486
Li S, Chen G, Jaakkola JJK, Williams G, Guo Y. Temporal change in the impacts of ambient temperature on preterm birth and stillbirth: Brisbane, 1994–2013. Sci Total Environ. 2018;634:579–85.
Auger N, Naimi AI, Smargiassi A, Lo E, Kosatsky T. Extreme heat and risk of early delivery among preterm and term pregnancies. Epidemiology. 2014;25:344–50.
Beltran AJ, Wu J, Laurent O. Associations of meteorology with adverse pregnancy outcomes: a systematic review of preeclampsia, preterm birth and birth weight. Int J Environ Res Public Health. 2013;11:91–172.
Ha S, Liu D, Zhu Y, Kim SS, Sherman S, Mendola P. Ambient temperature and early delivery of singleton pregnancies. Environ Health Perspect. 2017;125:453–9.
Basu R, Malig B, Ostro B. High ambient temperature and the risk of preterm delivery. Am J Epidemiol. 2010;172:1108–17.
Dadvand P, Basagaña X, Sartini C, Figueras F, Vrijheid M, De Nazelle A, et al. Climate extremes and the length of gestation. Environ Health Perspect. 2011;119:1449.
Schifano P, Lallo A, Asta F, De Sario M, Davoli M, Michelozzi P. Effect of Ambient Temperature and Air Pollutants on the Risk of Preterm Birth, Rome 2001–2010. 2013;61:77–87.
Strand LB, Barnett AG, Tong S. Maternal Exposure To Ambient Temperature and the Risks of Preterm Birth and Stillbirth in Brisbane, Australia. 2012;175:99–107.
Vicedo-Cabrera AM, Iñíguez C, Barona C, Ballester F. Exposure to elevated temperatures and risk of preterm birth in Valencia, Spain. Environ Res. 2014;134:210–7.
Wang MZ, Zheng S, He SL, Li B, Teng HJ, Wang SG, et al. The association between diurnal temperature range and emergency room admissions for cardiovascular, respiratory, digestive and genitourinary disease among the elderly: a time series study. Sci Total Environ. 2013;456–457:370–5.
Bruckner TA, Modin B. V\a agerö D. cold ambient temperature in utero and birth outcomes in Uppsala, Sweden, 1915–1929. Ann Epidemiol. 2014;24:116–21.
Ha S, Zhu Y, Liu D, Sherman S, Mendola P. Ambient temperature and air quality in relation to small for gestational age and term low birthweight. Environ Res. 2017;155:394–400.
Kloog I, Melly SJ, Coull BA, Nordio F, Schwartz JD. Using Satellite-Based Spatiotemporal Resolved Air Temperature Exposure To Study the Association Between Ambient Air Temperature and Birth Outcomes in Massachusetts. 2015; Available from: https://doi.org/10.1289/ehp.1308075
Lawlor DA, Leon DA, Smith GD. The association of ambient outdoor temperature throughout pregnancy and offspring birthweight: findings from the Aberdeen children of the 1950s cohort. BJOG Int J Obstet Gynaecol. 2005;112:647–57.
Murray LJ, O’reilly DPJ, Betts N, Patterson CC, Smith GD, Evans AE. Season and outdoor ambient temperature: effects on birth weight. Obstet Gynecol. 2000;96:689–95.
Hartig T, Catalano R. Cold Summer Weather, Constrained Restoration, and Very Low Birth Weight in Sweden, vol. 22; 2013. p. 68–74.
Lee SJ, Hajat S, Steer PJ, Filippi V. A time-series analysis of any short-term effects of meteorological and air pollution factors on preterm births in London, UK. Environ Res. 2008;106:185–94.
Spiegelman D. Approaches to uncertainty in exposure assessment in environmental epidemiology. Annu Rev Public Health. 2010;31:149–63.
Arcangeli T, Thilaganathan B, Hooper R, Khan KS, Bhide A. Neurodevelopmental delay in small babies at term: a systematic review. Ultrasound Obstet Gynecol. 2012;40:267–75.
Pallotto EK, Kilbride HW. Perinatal outcome and later implications of intrauterine growth restriction. Clin Obstet Gynecol. 2006;49:257–69.
Rosenfeld A, Dorman M, Schwartz J, Novack V, Just AC, Kloog I. Estimating daily minimum, maximum, and mean near surface air temperature using hybrid satellite models across Israel. Environ Res. 2017;159:297–312.
Ziv B, Saaroni H, Pargament R, Harpaz T, Alpert P. Trends in rainfall regime over Israel, 1975–2010, and their relationship to large-scale variability. Reg Environ Chang. 2014;14:1751–64.
Peel MC, Finlayson BL, McMahon TA. Updated world map of the Köppen-Geiger climate classification. Hydrol Earth Syst Sci Discuss. 2007;4:439–73.
Goldreich Y. The climate of Israel: observation. Res Appl Springer N Y. 2003.
Sade MY, Kloog I, Liberty IF, Katra I, Novack L, Novack V. Air pollution and serum glucose levels. A Population Based Study. 2015.
Sade MY, Novack V, Ifergane G, Horev A, Kloog I. Air pollution and ischemic stroke among young adults. Stroke. 2015;46:3348–53.
Dollberg S, Haklai Z, Mimouni FB, Gorfein I, Gordon E-S. Birth weight standards in the live-born population in Israel. Isr Med Assoc J. 2005;7:311–4.
Kloog I, Sorek-Hamer M, Lyapustin A, Coull B, Wang Y, Just AC, et al. Estimating daily PM 2.5 and PM 10 across the complex geo-climate region of Israel using MAIAC satellite-based AOD data. Atmos Environ. 2015;122:409–16.
ICBS. Isr Cent Bur Stat [Internet]. 2017; Available from: http://www.cbs.gov.il
Wilson A, Chiu Y-HM, Hsu H-HL, Wright RO, Wright RJ, Coull BA. Potential for Bias when estimating critical windows for air pollution in Children’s health. Am J Epidemiol. 2017;186:1281–9.
Wood SN. Fast stable restricted maximum likelihood and marginal likelihood estimation of semiparametric generalized linear models. J R Stat Soc B. 2011;73:3–36.
McGeehin MA, Mirabelli M. The potential impacts of climate variability and change on temperature-related morbidity and mortality in the United States. Environ Health Perspect. 2001;109:185.
Louis GMB, Grewal J, Albert PS, Sciscione A, Wing DA, Grobman WA, et al. Racial/ethnic standards for fetal growth: the NICHD fetal growth studies. Am J Obstet Gynecol. 2015;213:449–e1.
Browne VA, Julian CG, Toledo-Jaldin L, Cioffi-Ragan D, Vargas E, Moore LG. Uterine artery blood flow, fetal hypoxia and fetal growth. Phil Trans R Soc B. 2015;370:20140068.
Wolf J, Armstrong B. The association of season and temperature with adverse pregnancy outcome in two German states, a time-series analysis. PLoS One. 2012;7:e40228.
Arroyo V, Díaz J, Carmona R, Ortiz C, Linares C. Impact of air pollution and temperature on adverse birth outcomes: Madrid, 2001–2009. Environ Pollut. 2016;218:1154–61.
Díaz J, Arroyo V, Ortiz C, Carmona R, Linares C. Effect of environmental factors on low weight in non-premature births: a time series analysis. PLoS One. 2016;11:e0164741.
Yosef Y, Aguilar E, Alpert P. Detecting and adjusting artificial biases of long-term temperature records in Israel. Int J Climatol. 2018;38:3273–89.
This study was funded by the BSF Grant No. 2017277, EHF grant RPGA 1401 and Marie Currie grant FP7-PEOPLE-2013-CIG- STRPM.
Availability of data and materials
The datasets generated and analyzed during this study are not publically available due to IRB restrictions. Non-identifiable data are however available from the authors upon reasonable request and with permission from the IRB at Soroka University Medical Center.
Ethics approval and consent to participate
Approval to the study was granted by Soroka University Medical Center IRB (Ref: 0233–16). No consent from participants was required by the IRB.
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.