Heat Wave Vulnerability Mapping for India
Int. J. Environ. Res. Public Health 2017, 14(4), 357; doi:10.3390/ijerph14040357
The RAND Corporation, 1776 Main Street, Santa Monica, CA 90401, USA
Pardee RAND Graduate School, 1776 Main Street, Santa Monica, CA 90401, USA
Rollins School of Public Health, Emory University, Atlanta, GA 30322, USA
Indian Institute of Public Health, Gandhinagar, Gujarat 382042, India
Public Health Foundation of India, New Delhi 110070, India
Correspondence: Tel.: +1-301-393-0411 (ext. 6056)
Academic Editor: Peng Bi
Received: 7 January 2017 / Accepted: 28 March 2017 / Published: 30 March 2017
Assessing geographic variability in heat wave vulnerability forms the basis for planning appropriate targeted adaptation strategies. Given several recent deadly heatwaves in India, heat is increasingly being recognized as a public health problem. However, to date there has not been a country-wide assessment of heat vulnerability in India. We evaluated demographic, socioeconomic, and environmental vulnerability factors and combined district level data from several sources including the most recent census, health reports, and satellite remote sensing data. We then applied principal component analysis (PCA) on 17 normalized variables for each of the 640 districts to create a composite Heat Vulnerability Index (HVI) for India. Of the total 640 districts, our analysis identified 10 and 97 districts in the very high and high risk categories (> 2SD and 2-1SD HVI) respectively. Mapping showed that the districts with higher heat vulnerability are located in the central parts of the country. On examination, these are less urbanized and have low rates of literacy, access to water and sanitation, and presence of household amenities. Therefore, we concluded that creating and mapping a heat vulnerability index is a useful first step in protecting the public from the health burden of heat. Future work should incorporate heat exposure and health outcome data to validate the index, as well as examine sub-district levels of vulnerability.
heatwave; vulnerability; heat vulnerability index; vulnerability assessment; mapping; India
The Intergovernmental Panel on Climate Change (IPCC) report  highlights the projected increases in heatwave frequency, intensity and duration, and resulting deaths both globally and in India. Heatwave events have caused massive deaths in the past; the most famous among them are the European 2003 and Russian 2010 heatwaves, where tens of thousands died [2,3,4]. India has experienced several heatwaves, and most recently, just in the past two years, thousands have reportedly died . Research has documented an increase in cardiovascular , respiratory , and all-cause  mortality along with increases in ambulance calls and admissions [8,9] resulting from heatwave exposure. While most of the evidence is from North America and Europe, there is an emerging body of evidence from developing countries, including India , where heat wave deaths may currently be underestimated .
At the same time, heat-related deaths are preventable and prevention programs have been shown to be extremely cost effective . Population adaptation  along with preparedness measures have reduced mortality. Indeed, several cities and countries around the world have adopted heatwave preparedness plans . However, in India, this effort is limited to only a few cities . A broader preparedness strategy is particularly important given the large population, difficult local conditions, and limited adaptive capacity.
Health vulnerability can be conceptualized as complex and multidimensional . Vulnerability encompasses individual biophysical characteristics, as well as population-level socio-economic-environmental characteristics. These population measures have typically included measures of age, income, discrimination, social isolation, vegetation, and health characteristics . Incorporating multidimensional data can present a more comprehensive characterization of vulnerability. Given the considerable intra-country variations in these measures that exist, it is prudent to use this characterized vulnerability to identify communities in need of prioritized and focused interventions. Heat Vulnerability Indices (HVI) have been found to be a useful screening tool for targeting heat risk interventions .
Several international studies have explored vulnerabilities at the national, county, and city levels [17,18,19,20], but none have comprehensively examined India. Additionally, most of this work has been performed in the context of urban settings, while the majority of the Indian population resides in rural areas. To our knowledge only one Indian study looked at agricultural vulnerability using census data from 2001 . In fact, a recent review of heat vulnerability indices points out that that the majority of studies have been performed in Europe and the United States and recommends further study in other countries and regions to account for local context . We, therefore, aim to create and map an integrated district level heat vulnerability index for India that can be used to identify the most heat-vulnerable districts in the country.
2. Materials and Methods
2.1. Data Sources
Data was extracted from the Census of India 2011, District Level Household Survey (DLHS)-3 data , and from the Indian Space Research Organization ISRO server Bhuvan. District level census data was downloaded through the census of India portal. We chose districts as the unit for our analysis. While states are too heterogenous, districts are appropriate for planning purposes, and are the smallest unit for which we could get reliable data from multiple sources. We downloaded data for household amenities and the primary census abstract (PCA). Some DLHS variables were extracted from the Annual health survey report from the ministry of health and family welfare website. Satellite data was extracted from the ISRO server’s “Bhuvan” tool where the average vegetation fraction and Normalized Difference Vegetation Index (NDVI) images were layered with a district level India shapefile and these variables (mean, median, maximum, minimum, range, and standard deviation) were calculated for each of the district polygons using GIS software.
2.2. Choice of Variables
Our initial dataset had 140 variables for 640 districts. Based on the existing literature  and removal of duplicative variables that represented similar constructs, we shortlisted variables across the demographic, social, economic, health, and environmental domains. Researchers independently considered variables for inclusion based on knowledge from prior literature. Final variables were then compared and conflicts were resolved through discussion. We were left with 17 variables to be entered for PCA analysis. We did not assign weights to individual variables.
The shortlisted variables plausibly (and have been documented to) affect the heat-health relationship. The final list of variables employed in analysis along with their respective sources are listed in Table 1.
Table 1. Heat-health vulnerability data for 640 districts of India.
Demographic variables included those based on age and gender. For age extremes, we used the percentage of elderly and under five and for gender we used the sex ratio. While studies have documented elderly age to be a risk factor [3,9,23,24,25,26,27,28,29,30,31,32], there is an argument that children [32,33] could also be at a higher risk, given their higher metabolic rate.
Social class was represented through the percentage of people belonging to the scheduled castes (Dalits) and scheduled tribes (Adivasis). These groups are recognized by the Indian constitution as depressed classes and are a target for development and affirmative action programs. We used these as substitutes for race in the Indian context. Though there are no studies which highlight their vulnerability to extreme heat, this is perhaps from an absence of such literature. Studies between ethnic groups and heat-related deaths show equivocal effects [27,34,35]. This may be due to differential ownership of household amenities, chiefly air conditioning.
Socio economic variables were those of literacy, defined as an ability to read and write in any language; and occupational status of a worker, defined as producing goods and services. Education has been seen to be associated with heat mortality [6,30] perhaps from increased awareness and also as a proxy for socio economic status. Similarly, being in a worker status would increase the environmental exposure to heat since the majority of workers work in the agricultural sector as cultivators and agricultural laborers. Income was assessed through the percentage of people in the lowest wealth quintile. Several studies including those from Asia and India have documented the effects of poverty on heat-related deaths [6,25,29,32].
Household amenities were assessed through the presence of household amenities such as the presence of drinking water inside premises, living in a good house, having a mobile phone, radio and TV . While much research has gone into analyzing the effects of air-conditioning [34,36], we did not include this variable because its prevalence is quite low, electricity supply is irregular across the country, and approximately 300 million people have no power connection. However, many of the household amenities assessed require electricity and can serve as a proxy for electricity supply.
Land cover was assessed through Vegetation Fraction (VF) and Normalized Difference Vegetation Index (NDVI). Studies have documented the protective influence of green cover from heat island effects and heat deaths  and these have been included in other indices .
Population health was assessed through the immunization status children (12–23 months) fully immunized (BCG, 3 doses each of DPT, and Polio and Measles) (%) and presence of a health facility (sub-center) within 3 km. Preexisting illnesses have been documented to affect heat mortalities . Because we were unable to obtain prevalence of chronic disease that may be most closely linked to heat vulnerability (e.g., diabetes, cardiovascular diseases) at the district level, we used these factors as a proxy for overall underlying population health status.
2.3. Data Analysis
The data were merged using the census district numbers as unique identifiers. They were then manually checked for any discrepancies such as those between district names and were resolved through crosschecking. Where district level data were missing, we substituted state-wide averages for the districts.
Data were analyzed with STATA ver. 13, 2016 (StataCorp LP, College Station, TX, USA) and qGIS 2.10.1 (Open Source Geospatial Foundation Project). The map was created with R v. 3.3.2 (R Foundation for Statistical Computing, Vienna, Austria). All data used in this study were publicly available data and did not contain any individually identifiable information.
We calculated Spearman’s correlation coefficients between the vulnerability variables. We then employed the PCA technique to reduce the dimensionality of these variables. In order to perform PCA variables need to be on the same scale, we normalized all our variables by calculating their Z scores. All variables were in the same direction i.e., increasing value implies an increasing vulnerability. After application of PCA we tested for unexplained variation and adequacy of analysis by Kaiser-Meyer-Olkin (KMO) test. We retained four factors based on these criteria—Kaiser’s criteria of eigenvalues > 1, break in values in the scree plot test, and the variance explained by the factors. We tried factor rotations to increase the variability among factors. Individual factor scores were predicted for each district. In absence of any information with regards to the shape or interaction of factors, we assumed a linear relationship and calculated the heat vulnerability index by summing up these district level factor scores.
Table 1 demonstrates the sources and descriptive characteristics of the 17 raw variables included. Some of these variables showed interesting variations, for example, on average only 42.3% of households had drinking water inside their premises but it ranged from 2.4% to 93.8% across districts. Similarly TV ownership, immunization status, having nearby health sub-centers and vegetation also showed such marked variations.
Variable correlations showed many of the variables highly correlated with each other at the 0.001 significant level. (not shown here).
Table 2 shows the PCA results with Varimax rotation. The factors have been reduced to four dimensions. These correspond to demographic, socio-economic, vegetation, and health systems. The PCA led to four factors with primary loadings, these appeared to be (1) demographic; (2) socio-economic; (3) environmental; and (4) health factors. Demographic loadings were constituted of extremes of age, socio-economic loadings were driven by household amenities, environmental loadings were contributed by the VF and NDVI scores and health was driven by availability of health facilities nearby.
Table 2. Factor loadings from varimax rotation based on data from 640 districts.
Of the 17 variables that we started with, these four groups of factors were able to account for 78% of the total variation. The scree plot of eigenvalues showed a clean break at four components (Figure 1). This was also in agreement of the Kaiser criteria. The Kaiser–Meyer–Olkin (KMO) test shows adequacy of our analysis (>0.50).
Figure 1. PCA results. (a) Scree plot of Eigenvalues after PCA; (b) component loadings on orthogonal (Varimax) rotation; (c) score variables on orthogonal (Varimax) rotation.
The HVI calculated as a sum of the four individual factors for each district ranged from −11.8 to 9.4, it had a mean close to 0 and SD of 3.5. Figure 2 maps the vulnerability across the country. Spatial clustering of these “hot spots” is observed in central India. These districts have poor socio-economic and development indicators and appear to be high on the heat vulnerability index.
Figure 2. HVI mapping.
We classified these categories based on the SD scores as “very high” (>2SD), “high” (1–2SD), “high normal” (0–1SD), “low normal” (−1–0SD), “low” (−2–1SD), and “very low” (<−2SD). We chose this SD based classification instead of equal categorization to better represent the variation. Table 3 shows the number of districts according to HVI standard deviations.
Table 3. Number (%) of districts by HVI (Heat Vulnerability Index) Categories.
Ten districts had an HVI score of “very high” (>2SD), most of them in central India in the states of Madhya Pradesh and Chhattisgarh (Table 4). Twenty districts had an HVI score of “very low”; most of them were in the relatively developed states of Kerala and Goa and union territories of Lakshadweep.
Table 4. Districts with a “very high” HVI score.
This study provides a relative ranking of heat wave vulnerability for all districts in India. Although much is known about factors that contribute to vulnerability from other settings [16,17,18,19], there has been minimal research conducted within India on heat-related vulnerability. By coupling this knowledge with local context and using methods previously applied in other settings , we created an index that describes relative variation in heat-related vulnerability across all of India. This index can be used by planners, policy makers, and disaster mitigation experts to target climate adaptation efforts.
Similar to the findings of other international studies , our index too identified demographic, socio-economic, environmental, and health system factors. However, there are important differences in the choice of initial variables making this index useful to the Indian and developing country settings.
The high and very high HVI districts were in the central part of the country. With a higher tribal population, these states have been at the lower end of various health, education, economic and population growth indicators. They are referred to as the Empowered Action Group (EAG) states and often targeted for focused interventions. These land-locked, high HVI districts in the North and Central Indian plains are classically known as the “heat belt”.
While the use of air conditioning has been observed to have the greatest impact in reducing heat wave deaths in the US  it is unlikely to be a solution for India at least in the short term because of lack of a reliable and continuous power supply, the high cost and low penetration of air conditioning.
Suitable local adaptation strategies therefore need to be considered. These may include a range of measures, some of which have been discussed in the literature, such as public messaging (Radio, TV), mobile phone-based text messages, automated phone calls, and amber alerts; to others such as traditional adaptation practices of staying indoors, wearing comfortable clothes, and diets. These are often visible in terms of the housing design and construction material used. Simple design features such as having shaded windows and underground water storage tanks can be helpful. Use of insulator housing materials similarly can be an effective method of prevention. Having access to drinking water within housing premises and indoor toilets could be important. We chose several household amenities not just to proxy for income but also for their protective role.
For risk management, it would make sense to observe whether these identified areas of high vulnerability are also the same as those with higher temperatures and humidity. Similarly, had district level heat wave death data been available, we could have used it to validate this index. Our index does show high (>0.70) and significant correlations with literacy rates, low income status, TV ownership, having toilets and drinking water and open defecation practices. These could be seen as starting points framing local adaptation strategies. These correlations highlight the importance of interventions against other associated diseases such as gastrointestinal diseases in children and water-borne illnesses etc. There is a moderate correlation of 0.42 between HVI and average summer land surface temperatures (from satellite data) suggesting a relationship between higher temperatures and heat vulnerability. The index also shows moderate negative correlation (–0.46, p < 0.001) with urbanization signifying a possible greater vulnerability threat in rural areas. Since the majority of Indians reside in rural areas, this could have important implications. Outdoor workers have been identified as being at a greater risk during heatwaves. In rural settings, agricultural practices in different regions of India may also have diverse vulnerability patterns.
Some limitations of this study also arise from availability of data. Cardiovascular and/or respiratory diseases are more closely related to heat vulnerability but prevalence of chronic diseases at the district level was not available. Similarly, there was no pan India district level data on social isolation or electricity. For the three DLHS variables, we had missing data for the state of Nagaland. We used state averages instead. However, since Nagaland is a small state which has not reported heatwave deaths, we believe this substitution is unlikely to have major effects. In calculating HVI, we assumed a linear combination of factors with no weighting as a good first assumption. Inclusion of temperature as an exposure variable could have been helpful but because temperatures vary at country-wide levels in a thermally diverse India, it would serve to bias the index in favor of places with higher normal temperatures. Our approach is in line with established methods  for large areas. Also, district level temperature data is only collected for a small fraction of the total 640 districts. In view of still building evidence base from Indian temperature-mortality studies we have cited western literature and some Asian and Indian studies identifying vulnerability factors. This approach has also been demonstrated previously in the air pollution literature .
In many prior heat vulnerability studies, rural areas have been overlooked, but may have high vulnerability, and this may be especially important in India given the importance of poverty, and agricultural livelihoods in mediating the temperature mortality relationship. However, since the district level data includes both rural and urban areas; by aggregating them we may have missed the differences between these patterns of vulnerability [39,40]. Also, since we only had district level data, if such data was available at the finer block (Taluka) level, we would have a better identification of vulnerable areas. Similarly, intra-city vulnerability patterns would have been interesting to observe given availability of urban ward level data. Given data availability, future work could also identify areas using the Koeppen climate classification assuming that the warm, dry, arid, and humid areas are more vulnerable. As research continues, we may identify more complex relationships and therefore we could conceive of other heat vulnerability indices with non-linear relationships and differential weights.
Despite the above discussed limitations there are several strengths to this work. This is the first study to look at heat vulnerability across India. It provides a preliminary screening to target heat-health and climate adaptation efforts. This methodology can be used to for further investigations into vulnerability.
We developed a heat wave vulnerability index that aggregates indicators across several dimensions for all districts in India. This index can be used in initial efforts to target resources for adaptation efforts. Most heat preparedness plans are designed and implemented at a local level, and this index can help identify metropolitan areas that are at highest risk. Since heat wave vulnerability varies across spatial scales, our methods can be extended for sub-district level analysis and modified to develop urban and rural indices. Further work within this context should include testing the sensitivity of our linear combination assumption and validating this index to health outcomes.
This study was funded by a seed grant from RAND’s Center for Asia Pacific Policy. These included open access publication charges. We thank Ajit Rajiva, Yale University, for his help in generating land surface temperature estimates.
G.A. and J.M. conceived and designed the research; G.A., J.M., and S.S. analyzed the data; P.S.G. and D.M. contributed to the analysis; G.A. wrote the paper; all authors reviewed and commented on the paper at several stages. All authors made significant contributions towards this work.
Conflicts of Interest
The authors declare no conflict of interest. The founding sponsors had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, and in the decision to publish the results.
- Parry, M.L. Climate Change 2007-Impacts, Adaptation and Vulnerability: Working Group II Contribution to the Fourth Assessment Report of the Ipcc; Cambridge University Press: Cambridge, UK, 2007; Volume 4. [Google Scholar]
- Robine, J.-M.; Cheung, S.L.K.; Le Roy, S.; Van Oyen, H.; Griffiths, C.; Michel, J.-P.; Herrmann, F.R. Death toll exceeded 70,000 in Europe during the summer of 2003. Comptes Rendus Biol. 2008, 331, 171–178. [Google Scholar] [CrossRef] [PubMed]
- Fouillet, A.; Rey, G.; Laurent, F.; Pavillon, G.; Bellec, S.; Guihenneuc-Jouyaux, C.; Clavel, J.; Jougla, E.; Hémon, D. Excess mortality related to the august 2003 heat wave in France. Int. Arch. Occup. Environ. Health 2006, 80, 16–24. [Google Scholar] [CrossRef] [PubMed]
- Shaposhnikov, D.; Revich, B.; Bellander, T.; Bedada, G.B.; Bottai, M.; Kharkova, T.; Kvasha, E.; Lezina, E.; Lind, T.; Semutnikova, E. Mortality related to air pollution with the moscow heat wave and wildfire of 2010. Epidemiology2014, 25, 359–364. [Google Scholar] [CrossRef] [PubMed]
- Liberto, T.D. India heat wave kills thousands. In Event Tracker; Climate Government: Washington DC, USA, 2015. [Google Scholar]
- Curriero, F.C.; Heiner, K.S.; Samet, J.M.; Zeger, S.L.; Strug, L.; Patz, J.A. Temperature and mortality in 11 cities of the eastern United States. Am. J. Epidemiol. 2002, 155, 80–87. [Google Scholar] [CrossRef] [PubMed]
- Mastrangelo, G.; Fedeli, U.; Visentin, C.; Milan, G.; Fadda, E.; Spolaore, P. Pattern and determinants of hospitalization during heat waves: An ecologic study. BMC Public Health 2007, 7, 200. [Google Scholar] [CrossRef] [PubMed]
- Hess, J.J.; Saha, S.; Luber, G. Summertime acute heat illness in us emergency departments from 2006 through 2010: Analysis of a nationally representative sample. Environ. Health Perspect. 2014, 122, 1209. [Google Scholar] [PubMed]
- Knowlton, K.; Rotkin-Ellman, M.; King, G.; Margolis, H.G.; Smith, D.; Solomon, G.; Trent, R.; English, P. The 2006 california heat wave: Impacts on hospitalizations and emergency department visits. Environ. Health Perspect. 2009, 117, 61. [Google Scholar] [CrossRef] [PubMed]
- Azhar, G.S.; Mavalankar, D.; Nori-Sarma, A.; Rajiva, A.; Dutta, P.; Jaiswal, A.; Sheffield, P.; Knowlton, K.; Hess, J.J. Heat-related mortality in india: Excess all-cause mortality associated with the 2010 ahmedabad heat wave. PLoS ONE 2014, 9, e91831. [Google Scholar]
- Azhar, G.S.; Rajiva, A.; Dutta, P.; Mavalankar, D. Research and policy disconnect for heat wave deaths in India. Curr. Sci. 2014, 107, 1091–1092. [Google Scholar]
- Ebi, K.L.; Teisberg, T.J.; Kalkstein, L.S.; Robinson, L.; Weiher, R.F. Heat watch/warning systems save lives: Estimated costs and benefits for Philadelphia 1995–1998. Bull. Am. Meteorol. Soc. 2004, 85, 1067–1073. [Google Scholar] [CrossRef]
- Sheridan, S.C.; Dixon, P.G. Spatiotemporal trends in human vulnerability and adaptation to heat across the United States. Anthropocene 2016. [Google Scholar] [CrossRef]
- Knowlton, K.; Kulkarni, S.P.; Azhar, G.S.; Mavalankar, D.; Jaiswal, A.; Connolly, M.; Nori-Sarma, A.; Rajiva, A.; Dutta, P.; Deol, B. Development and implementation of South Asia’s first heat-health action plan in Ahmedabad (Gujarat, India). Int. J. Environ. Res. Public Health 2014, 11, 3473–3492. [Google Scholar] [CrossRef] [PubMed]
- Watts, N.; Adger, W.N.; Agnolucci, P.; Blackstock, J.; Byass, P.; Cai, W.; Chaytor, S.; Colbourn, T.; Collins, M.; Cooper, A. Health and climate change: Policy responses to protect public health. Lancet 2015, 386, 1861–1914. [Google Scholar] [CrossRef]
- Bao, J.; Li, X.; Yu, C. The construction and validation of the heat vulnerability index, a review. Int. J. Environ. Res. Public Health 2015, 12, 7220–7234. [Google Scholar] [CrossRef] [PubMed]
- Reid, C.E.; O’Neill, M.S.; Gronlund, C.J.; Brines, S.J.; Diez-Roux, A.V.; Brown, D.G.; Schwartz, J.D. Mapping community determinants of heat vulnerability. Environ. Health Perspect. 2009, 117, 1730–1736. [Google Scholar] [CrossRef] [PubMed]
- Johnson, D.P.; Stanforth, A.; Lulla, V.; Luber, G. Developing an applied extreme heat vulnerability index utilizing socioeconomic and environmental data. Appl. Geogr. 2012, 35, 23–31. [Google Scholar] [CrossRef]
- Wolf, T.; McGregor, G. The development of a heat wave vulnerability index for london, United Kingdom. Weather Clim. Extrem. 2013, 1, 59–68. [Google Scholar] [CrossRef]
- Madrigano, J.; Ito, K.; Johnson, S.; Kinney, P.L.; Matte, T. A case-only study of vulnerability to heat wave-related mortality in New York city (2000–2011). Environ. Health Perspect. 2015. [Google Scholar] [CrossRef] [PubMed]
- Rama Rao, C.; Raju, B.; Subba Rao, A.; Rao, K.; Rao, V.; Ramachandran, K.; Venkateswarlu, B.; Sikka, A.; Srinivasa Rao, M.; Maheswari, M. A district level assessment of vulnerability of indian agriculture to climate change. Curr. Sci. 2016, 110, 1939. [Google Scholar] [CrossRef]
- International Institute for Population Sciences. District Level Household and Facility Survey 2007–2008; DLHS-3; Ministry of Health and Family Welfare, Government of India: New Delhi, India, 2010.
- Conti, S.; Masocco, M.; Meli, P.; Minelli, G.; Palummeri, E.; Solimini, R.; Toccaceli, V.; Vichi, M. General and specific mortality among the elderly during the 2003 heat wave in Genoa (Italy). Environ. Res. 2007, 103, 267–274. [Google Scholar] [CrossRef] [PubMed]
- Hutter, H.-P.; Moshammer, H.; Wallner, P.; Leitner, B.; Kundi, M. Heatwaves in Vienna: Effects on mortality. Wien. Klin. Wochenschr. 2007, 119, 223–227. [Google Scholar] [CrossRef] [PubMed]
- Naughton, M.P.; Henderson, A.; Mirabelli, M.C.; Kaiser, R.; Wilhelm, J.L.; Kieszak, S.M.; Rubin, C.H.; McGeehin, M.A. Heat-related mortality during a 1999 heat wave in Chicago. Am. J. Prev. Med. 2002, 22, 221–227. [Google Scholar] [CrossRef]
- Stafoggia, M.; Forastiere, F.; Agostini, D.; Biggeri, A.; Bisanti, L.; Cadum, E.; Caranci, N.; de’Donato, F.; De Lisio, S.; De Maria, M. Vulnerability to heat-related mortality: A multicity, population-based, case-crossover analysis. Epidemiology 2006, 17, 315–323. [Google Scholar] [CrossRef] [PubMed]
- Whitman, S.; Good, G.; Donoghue, E.R.; Benbow, N.; Shou, W.; Mou, S. Mortality in Chicago attributed to the July 1995 heat wave. Am. J. Public Health 1997, 87, 1515–1518. [Google Scholar] [CrossRef] [PubMed]
- Centers for Disease Control and Prevention. Heat-related deaths—Los Angeles County, California, 1999–2000, and United States, 1979–1998. MMWR: Morb. Mortal. Wkly. Rep. 2001, 50, 623–626. [Google Scholar]
- Kim, Y.; Joh, S. A vulnerability study of the low-income elderly in the context of high temperature and mortality in Seoul, Korea. Sci. Total Environ. 2006, 371, 82–88. [Google Scholar] [CrossRef] [PubMed]
- Medina-Ramón, M.; Zanobetti, A.; Cavanagh, D.P.; Schwartz, J. Extreme temperatures and mortality: Assessing effect modification by personal characteristics and specific cause of death in a multi-city case-only analysis. Environ. Health Perspect. 2006, 1331–1336. [Google Scholar] [CrossRef]
- Semenza, J.C.; McCullough, J.E.; Flanders, W.D.; McGeehin, M.A.; Lumpkin, J.R. Excess hospital admissions during the July 1995 heat wave in Chicago. Am. J. Prev. Med. 1999, 16, 269–277. [Google Scholar] [CrossRef]
- Tran, K.V.; Azhar, G.S.; Nair, R.; Knowlton, K.; Jaiswal, A.; Sheffield, P.; Mavalankar, D.; Hess, J. A cross-sectional, randomized cluster sample survey of household vulnerability to extreme heat among slum dwellers in Ahmedabad, India. Int. J. Environ. Res. Public Health 2013, 10, 2515–2543. [Google Scholar] [CrossRef] [PubMed]
- Kakkad, K.; Barzaga, M.L.; Wallenstein, S.; Azhar, G.S.; Sheffield, P.E. Neonates in Ahmedabad, India, during the 2010 heat wave: A climate change adaptation study. J. Environ. Public Health 2014, 2014, 946875. [Google Scholar] [CrossRef] [PubMed]
- Braga, A.L.F.; Zanobetti, A.; Schwartz, J. The time course of weather-related deaths. Epidemiology 2001, 12, 662–667. [Google Scholar] [CrossRef] [PubMed]
- Kalkstein, L.S.; Davis, R.E. Weather and human mortality: An evaluation of demographic and interregional responses in the United States. Ann. Assoc. Am. Geogr. 1989, 79, 44–64. [Google Scholar] [CrossRef]
- Barreca, A.; Clay, K.; Deschenes, O.; Greenstone, M.; Shapiro, J.S. Adapting to Climate Change: The Remarkable Decline in the US Temperature-Mortality Relationship over the 20th Century; National Bureau of Economic Research: Cambridge, MA, USA, 2013. [Google Scholar]
- Kilbourne, E.M.; Choi, K.; Jones, T.S.; Thacker, S.B. Risk factors for heatstroke: A case-control study. JAMA 1982, 247, 3332–3336. [Google Scholar] [CrossRef] [PubMed]
- Balakrishnan, K.; Rajarathnam, U. Public Health and Air Pollution in Asia (PAPA): Coordinated Studies of Short-Term Exposure to Air Pollution and Daily Mortality in Two INDIAN Cities; Research Report; Health Effects Institute: Boston, MA, USA, 2011. [Google Scholar]
- Sheridan, S.C.; Dolney, T.J. Heat, mortality, and level of urbanization: Measuring vulnerability across Ohio, USA. Clim. Res. 2003, 24, 255–265. [Google Scholar] [CrossRef]
- Madrigano, J.; Jack, D.; Anderson, G.B.; Bell, M.L.; Kinney, P.L. Temperature, ozone, and mortality in urban and non-urban counties in the northeastern United States. Environ. Health 2015, 14, 3. [Google Scholar] [CrossRef] [PubMed]
© 2017 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license ( http://creativecommons.org/licenses/by/4.0/).
Leave a Reply