Mortality risk during heat waves in the summer 2013-2014 in 18 provinces of Argentina: Ecological study

Risco de mortalidade durante ondas de calor no verão 2013-2014 em 18 províncias da Argentina: Estudo ecológico

Francisco Chesini Natalia Herrera Maria de los Milagros Skansi Carolina González Morinigo Silvia Fontán Francisco Savoy Ernesto de Titto About the authors

Abstract

Increased frequency of heat waves (HWs) is one of the prominent consequences of climate change. Its impact on human health has been mostly reported in the northern hemisphere but has been poorly studied in the southern hemisphere. The aim of this study was to analyze the effects of the HWs waves occurred in the warm season 2013-14 on mortality in the center-north region of Argentina, where 22 million people live. It was carried out an observational study of ecological-type contrasting the mortality occurred during the HWs of the summer 2013-14 with the mortality in the summers 2010-11 to 2012-13, free from HWs. The mortality was analyzed according to the following variables: place of residence, age, sex and cause of death. During the HWs of the summer 2013-14, 1877 (RR=1.23, 95%CI 1.20-1.28) deaths in excess were registered. Moreover, the death risk significantly increased in 13 of the 18 provinces analyzed. The mortality rates by sex revealed heterogeneous behaviour regarding both the time and spatial scale. The death risk increased with age; it was particularly significant in four provinces for the 60-79 years group and in six provinces in people of 80 years and over. The death causes that showed significantly increments were respiratory, cardiovascular, renal diseases and diabetes.

Key words:
Mortality; Heatwave; Climate change; Environmental health; Death risk

Resumo

O aumento da frequência das ondas de calor (OsC) é uma das consequências proeminentes das alterações climáticas. O seu impacto na saúde humana tem sido relatado principalmente no hemisfério norte, mas tem sido mal estudado no hemisfério sul. O objetivo deste estudo foi analisar os efeitos das OsC ocorridas no verão de 2013-14 sobre a mortalidade na região centro-norte da Argentina, onde 22 milhões de pessoas vivem. Foi realizado um estudo observacional de tipo ecológico que contrasta a mortalidade ocorrida durante as OsC do verão 2013-14 com a mortalidade nos verões 2010-11 a 2012-13. A mortalidade foi analisada de acordo com local de residência, idade, sexo e causa de morte. Durante as OsC do verão de 2013-14, registaram-se 1.877 (RR=1.23, 95%IC 1.20-1.28) mortes em excesso. O risco de morte (RM) aumentou significativamente em 13 das 18 províncias analisadas. As taxas de mortalidade por sexo revelaram um comportamento heterogêneo tanto no tempo como na escala espacial. O RM foi aumentado com a idade; este aumento foi particularmente significativo em quatro províncias para o grupo de 60-79 anos e em seis províncias em pessoas com mais de 80 anos. As causas de morte que mostraram significativamente incrementos foram: doenças respiratórias, cardiovasculares, renais e diabetes.

Palavras-chave:
Mortalidade; Onda de calor; Mudanças climáticas; Saúde ambiental; Risco de morte

Introduction

The health sector is paying increasing attention to the close relationship between climate change and human health. Thus, in 2008, the World Health Assembly recognized climate change as an environmental determinant of health11 World Health Organization (WHO). 61ª Asamblea Mundial de la Salud. Cambio climático y salud [Internet]. 2008 [cited 2020 mar 19]. Available from: http://apps.who.int/iris/bitstream/10665/26288/1/A61_R19-sp.pdf.
http://apps.who.int/iris/bitstream/10665...
and in 2009 the Ministers of Health from MERCOSUR established an Action Strategy to protect human health from climate change effects22 Mercado Común del Sur. Resolución MERCOSUR nº 12/2009. Estrategia Cambio Climático y Salud. Montevideo [Internet]. 2009 [cited 2020 mar 19]. Available from: http://200.40.51.218/SAM/GestDoc/PubWeb.nsf/Reunion?OpenAgent&id=CD00D62A0920FFD7032576A50049C42A&lang=ESP.
http://200.40.51.218/SAM/GestDoc/PubWeb....
.

Among the manifestations of climate change with significant impact on human and natural systems stands out an increasing duration and frequency of heat waves (HWs)33 McGregor GR, Bessemoulin P, Ebi K, Menne B. Heatwaves and Health: Guidance on Warning-System Development. Geneva: WHO; 2015., which been documented on several occasions. For example: the HW that occurred in Chicago in July 1995 caused a 147% increase in deaths (during the most intense heat - July 14 through 20 -, there were 485 heat-related deaths and 739 excess deaths)44 Whitman S, Good G, Donoghue ER, Benbow N, Shou W, Mou S. Mortality in Chicago Attributed to the July 1995 Heat Wave. Am J Public Health 1997; 87:1515-1518.; in Europe in 2003, the HW death rates increased 4 to 5 times causing around additional 70.000 deaths55 Robine JM, Cheung SL, Le Roy S, Van Oyen H, Griffiths C, Michel JP, Herrmann FR. Death toll exceeded 70,000 in Europe during the summer of 2003. C R Biol 2008; 331:171-178.,66 Oudin Åström D, Forsberg B, Rocklöv J. Heat wave impact on morbidity and mortality in the elderly population: A review of recent studies. Maturitas 2011; 69:99-105.; and in 2010 the Russian Federation had a HW that lasted 44 days and caused 10.000 excess deaths only in Moscow77 Shaposhnikov D, Revich B, Bellander T, Bedada GB, Bottai M, Kharkova T, Kvasha E, Lezina E, Lind T, Semutnikova E, Pershagen G. Mortality Related to Air Pollution with the Moscow Heat Wave and Wildfire of 2010. Epidemiology 2014; 25:359-364..

Several studies showed the increase in both hospitalizations and mortality during HW days, with risk increases for cardiovascular, cerebrovascular, respiratory and kidney diseases88 D'Ippoliti D, Michelozzi P, Marino C, Donato F, Menne B, Katsouyanni K, Kirchmayer U, Analitis A, Medina-Ramón M, Paldy A, Atkinson R, Kovats S, Bisanti L, Schneider A, Lefranc A, Iñiguez C, Perucci CA. The impact of heat waves on mortality in 9 European cities: results from the EuroHEAT project. Environ Health 2010; 9:37.

9 McGeehin MA, Mirabelli M. The Potential Impacts of Climate Variability and Change on Temperature-Related Morbidity and Mortality in the United States. Environ Health Persp 2001; 109:185-189.
-1010 Smith KR, Woodward A, Campbell-Lendrum D, Chadee DD, Honda Y, Liu Q, Olwoch JM, Revich B, Sauerborn R, Aranda C, Berry H, Butler C, Chafe Z, Cushing L, Ebi KL, Kjellstrom T, Kovats S, Lindsay G, Lipp E, McMichael T, Murray V, Sankoh O, O'Neill M, Shonkoff SB, Sutherland J, Yamamoto S. Human Health: Impacts, Adaptation, and Co-Benefits. In: Field CB, Barros VR, Dokken DJ, Mach KJ, Mastrandrea MD, Bilir TE, Chatterjee M, Ebi KL, Estrada YO, Genova RC, Girma B, Kissel ES, Levy AN, MacCracken S, Mastrandrea PR, White LL, editores. Climate Change 2014: Impacts, Adaptation and Vulnerability. Part A: Global and Sectorial Aspects. Contribution of Working Group II to the Fifth Assessment Report of the Intergovernmental Panel on Climate Change. Cambridge UK: Cambridge University Press; 2014. p. 709-754.. Although all individuals are potentially exposed to a HW, the risk levels can be modified by other factors such as living in large cities, socioeconomic conditions, gender, behaviour and even the consumption of medicines33 McGregor GR, Bessemoulin P, Ebi K, Menne B. Heatwaves and Health: Guidance on Warning-System Development. Geneva: WHO; 2015.,1111 Montero JC, Mirón IJ, Criado JJ, Linares C, Díaz J. Difficulties of defining the term "heat wave" in public health. Int J Environ Health Res 2013; 23:377-379.. The most vulnerable groups are the elderly, the children, the chronically ill and people working outdoors33 McGregor GR, Bessemoulin P, Ebi K, Menne B. Heatwaves and Health: Guidance on Warning-System Development. Geneva: WHO; 2015.,88 D'Ippoliti D, Michelozzi P, Marino C, Donato F, Menne B, Katsouyanni K, Kirchmayer U, Analitis A, Medina-Ramón M, Paldy A, Atkinson R, Kovats S, Bisanti L, Schneider A, Lefranc A, Iñiguez C, Perucci CA. The impact of heat waves on mortality in 9 European cities: results from the EuroHEAT project. Environ Health 2010; 9:37..

Since HW lacks of an internationally agreed definition1111 Montero JC, Mirón IJ, Criado JJ, Linares C, Díaz J. Difficulties of defining the term "heat wave" in public health. Int J Environ Health Res 2013; 23:377-379.,1212 Perkins SE. A review on the scientific understanding of heatwaves. Their measurement, driving mechanisms, and changes at the global scale. Atmospheric Res 2015; 164-165:242-267, the National Meteorological Service of Argentina (NMSA) defines HW as the period in which the maximum and minimum temperatures equal or exceed for at least 3 consecutive days the 90th percentile, calculated from daily data during the months from October to March of the period 1961-20101313 Servicio Meteorológico Nacional (SMN). Argentina. Datos - Clima de Argentina-Olas de calor (Heat waves) [Internet]. [cited 2020 mar 19]. Available from: https://www.smn.gob.ar/estadisticas.
https://www.smn.gob.ar/estadisticas...
. In line with the evolution of global warming, in this period the frequency of HW has increased in the north and east of the country, showing a pronounced decadal variability with the highest numbers in the decade 2001-20101414 Centro de Investigaciones del Mar y la Atmósfera (CIMA). Modelos Climáticos: Resumen Ejecutivo (Climate Models: Executive Summary). Secretaría de Ambiente y Desarrollo Sustentable [Internet]. 2015 [cited 2020 mar 20]. Available from: http://3cn.cima.fcen.uba.ar/3cn_informe.php
http://3cn.cima.fcen.uba.ar/3cn_informe....
,1515 Rusticucci M, Kysely J, Almeira G, Lhotka O. Long-term variability of heat waves in Argentina and recurrence probability of the severe 2008 heat wave in Buenos Aires. Theor Appl Climatol 2015; 124:679-689..

The summer season 2013-2014 was particularly hot. Starting on the month of December 2013, the persistence record for days with extremely high temperatures was exceeded, mainly in the central and northern part of the national territory.

Intraseasonal variability can explain on average at least 32% of summer temperature variance. Moreover, 73% of the HWs in subtropical South America develop in association with an active South Atlantic Convergence Zone (SACZ)1616 Cerne SB, Vera CS. Influence of the intraseasonal variability on heat waves in subtropical South America. Clim Dyn 2011; 36:2265-2277.. The development of the anticyclonic circulation over subtropical South America appears to be strengthened by the subsidence conditions promoted by the active SACZ, which result in temperature rise in the region under relatively dry conditions. As result of these land and atmospheric conditions, a high coupling (water-limited) regime was imposed, promoting the re-amplification of hot spells that resulted in mega HW episodes1717 Geirinhas JL, Russo A, Libonati R, Sousa PM, Miralles DG, Trigo RM. Recent increasing frequency of compound summer drought and heatwaves in Southeast Brazil. Environ Res Lett 2021; 16:034036.. Several studies analyzed strong and long HWs episodes occurred in the austral summer season 2013-2014 over Southern South America1313 Servicio Meteorológico Nacional (SMN). Argentina. Datos - Clima de Argentina-Olas de calor (Heat waves) [Internet]. [cited 2020 mar 19]. Available from: https://www.smn.gob.ar/estadisticas.
https://www.smn.gob.ar/estadisticas...
,1717 Geirinhas JL, Russo A, Libonati R, Sousa PM, Miralles DG, Trigo RM. Recent increasing frequency of compound summer drought and heatwaves in Southeast Brazil. Environ Res Lett 2021; 16:034036.

18 Cerne B, Vera C, Álvarez M. Influencia de la variabilidad intraestacional en el desarrollo de la ola de calor de diciembre de 2013 en el este de Sudamérica. In: Congreso Argentino de Meteorología (CONGREMET) XII. Mar del Plata; 2015 maio 26-29.
-1919 Osman M, Álvarez M. Subseasonal prediction of the heat wave of December 2013 in Southern South America by the POAMA and BCC CPS models. Climate Dynamics 2018; 50(1-2):67-81..

Some studies have been carried out in Argentina on the relationship between meteorological conditions and health. In one of them it was investigated the short-term effects of temperature and atmospheric carbon monoxide changes in daily mortality in Buenos Aires City (CABA)2020 Abrutzky R, Dawidowski L, Matus P, Romero-Lankao P. Health effects on climate and air pollution in Buenos Aires. A first series analysis. J Environ Prot 2012; 3:262-271.. Other study analyzed the relationship between temperature and general mortality in CABA and Rosario, the two most populated cities in the country, finding that the relationship follows a “U” form with increases in deaths at both extremes of temperature (cold and warm)2121 Almeira G, Rusticucci M, Suaya M. Relación entre mortalidad y temperaturas extremas en Buenos Aires y Rosario. Meteorológica 2016; 41:65-79.. Recently, some of us analyzed mortality in CABA during HWs, observing a significant increase in the risk of dying, with predominant impact in the mortality in women and elder age2222 Chesini F, Abrutzky R, de Titto EH. Mortalidad por olas de calor en la Ciudad de Buenos Aires, Argentina (2005-2015). Cad Saude Publica 2019; 35:e00165218..

Given the increasing frequency of HWs, it is necessary to identify the local epidemiological profile and their sanitary burden in order to design proper health promotion strategies to alleviate their impact on the population. For this reason, the aim of this study was to analyze the effects of the HWs occurred in the warm season 2013-14 on mortality in the center-north region of Argentina, where 22 million people live.

Materials and methods

This is an ecological type observational study and the units of analysis are the population of each jurisdiction (17 Provinces and CABA, which for practical purposes was also considered a province); as presented in Figure 1.

Figure 1
Geographical distribution and number of days with high temperature.

Study area

The geographic scope of the present study was based on previous studies that analyzed the long-term variability of HWs1313 Servicio Meteorológico Nacional (SMN). Argentina. Datos - Clima de Argentina-Olas de calor (Heat waves) [Internet]. [cited 2020 mar 19]. Available from: https://www.smn.gob.ar/estadisticas.
https://www.smn.gob.ar/estadisticas...
,1515 Rusticucci M, Kysely J, Almeira G, Lhotka O. Long-term variability of heat waves in Argentina and recurrence probability of the severe 2008 heat wave in Buenos Aires. Theor Appl Climatol 2015; 124:679-689.. Thus, we revised data from the 57 meteorological stations distributed on the provinces located north of the 40th parallel south.

Argentina is a highly urbanized country; according to the last national census (2010) 91% of the population lives in cities2323 Instituto Nacional de Estadística y Censos (INDEC). Argentina. Censo Nacional 2010 de Población, Hogares y Viviendas [Internet]. [cited 2020 mar 18]. Available from: http://www.indec.gov.ar/index.asp.
http://www.indec.gov.ar/index.asp...
. In Argentina the life expectancy is 75.3 years (78.8 for women and 72.8 for men)2424 Dirección Nacional de Estadísticas e Información de la Salud (DEIS). Indicadores Básicos Argentina 2014 [Internet]. Buenos Aires: Ministerio de Salud de la Nación; 2015 [cited 2020 mar 20]. Available from: http://www.deis.gov.ar/publicaciones/archivos/indicadores_2014.pdf.
http://www.deis.gov.ar/publicaciones/arc...
. Since the 1970s, a population aging process has been observed: people over 65 years old represented 7.0% of the national population in 1970 and 10.2% in 2010.

The Patagonian provinces were excluded from this study; since the HWs that can be recorded there are not of sanitary relevance because the threshold values of the minimum temperature are too low. The province of Buenos Aires was also excluded from this study since for administrative reasons mortality cases in this province are charged monthly and there is no way to get access to the actual date of the deaths from secondary data. Moreover, given its extent (approximately 307.000 km²) HWs never cover the whole province.

Analysis unit

The units of analysis are delimited from the political division. Data collection units are composed of the individual cases of death that occurred during the HWs. These individual cases were integrated into units of analysis called provinces, which include the deaths that occurred in each of the 18 provinces (Figure 1).

Analysis

We first identified HWs that occurred in the summer 2013-2014 in the study area, based on the criteria defined by the NMSA1313 Servicio Meteorológico Nacional (SMN). Argentina. Datos - Clima de Argentina-Olas de calor (Heat waves) [Internet]. [cited 2020 mar 19]. Available from: https://www.smn.gob.ar/estadisticas.
https://www.smn.gob.ar/estadisticas...
. Once the HWs were identified, maps were drawn indicating the number of days with HW, to determine the areas affected by the events and the start and end dates of each HW for each jurisdiction. Then, a spatial overlay analysis was carried out in a Geographic Information System (GIS) to be able to indicate the affected departments in each province and to undertake an individualized analysis of them. The layer in “shape” format with the division by departments was provided by the National Geographic Institute.

To analyze mortality we selected natural death cases registered in the categories included between codes A00 and R99 according to the 10th Revision of the International Classification of Diseases (ICD-10), excluding external causes of death88 D'Ippoliti D, Michelozzi P, Marino C, Donato F, Menne B, Katsouyanni K, Kirchmayer U, Analitis A, Medina-Ramón M, Paldy A, Atkinson R, Kovats S, Bisanti L, Schneider A, Lefranc A, Iñiguez C, Perucci CA. The impact of heat waves on mortality in 9 European cities: results from the EuroHEAT project. Environ Health 2010; 9:37.,2525 Culqui DR, Díaz J, Simón F, Linares C. Análisis de los impactos de las olas de calor sobre la mortalidad de la ciudad de Madrid durante el período 1990-2009. Rev Esp Salud Pública 2013; 87:277-282., from the data basis provided by the National Ministry of Health.

The explicative variable (presence of HW) and the outcome variable (mortality) were studied simultaneously2626 Beaglehole R, Bonita R, Kjellström T. Epidemiologia Básica. Washington D.C.: PAHO; 1994.. The study compared the daily mortality that occurred during each HW in the warm semester 2013-2014 (October 2013 to March 2014) with the average of the daily mortality that occurred during the same dates on the three previous warm seasons (2010-2011 to 2012-2013). The difference between the numbers of deaths during both periods was taken as excessive deaths. It must be noted that on selecting the reference period we required that they lacked HWs so that the comparisons were in the same population and social environment2727 Carracedo-Martínez E, Tobías A, Saez M, Taracido M, Figueras A. Fundamentos y aplicaciones del diseño de casos cruzados. Gac Sanit 2009; 23:161-165.. The number of deaths according to the cause of death, age at the time of death, sex and place of residence were taken as response variables.

We calculated the frequencies of cases and rates by age group, sex and cause of death, for HWs and for the reference period. In addition, the relative risks (RR) were calculated as the ratio of mortality rates during each HW and mortality rates in the absence of HW (average of the previous three years). It was also calculated its confidence interval of 95% (95%CI).

Mortality by age was grouped into five categories: 0 to 19; 20 to 39; 40 to 59; 60 to 79 and 80 and more years, as proposed by the specialized bibliography2828 Martín JF. Los factores definitorios de los grandes grupos de edad de la población: tipos, subgrupos y umbrales. Scripta Nova 2005; IX:190..

Information sources

The information analyzed comes from secondary sources. Statistical data on mortality from 2010 to 2014 were provided by the Direction of Health Statistics and Information of the National Ministry of Health. Population data were obtained from the National Institute of Statistics and Census and the meteorological information (maximal and minimal daily temperatures to build the HW dichotomic variable, following NMSA definition) from 1961 to 2010 for the reference period and from the warm season 2013-14 for the study was provided by the NMSA.

Tools to process information were: “Epiinfo” for Windows version 3.5.3, “ArcGis” version 9.3.1 and Surfer 7 for making maps and package Microsoft Office 2007.

Results

Analysis of heat waves in summer 2013-2014

This research analyzed meteorological data from Argentina during the warm semester for the southern hemisphere (October to March) and focused the study in the months of December 2013 and January-February 2014 since there were no HWs in October and November 2013 and March 2014.

The temperature during the summer 2013-2014 was especially high in practically the whole country, as it can be seen on Figure 1, where the number of days within each month with both maximum and minimum temperatures simultaneously exceeding 90th percentile is presented.

The HWs registered in each of the provinces with their start and end dates, their duration and comparison with the HWs that occurred in the 1961-2010 period are presented in Table 1. It is indicated when the HW was higher than the 80th, 85th, 90th or 95th percentile in duration, maximum temperature average (T. max) and/or minimum temperature average (T. min). For example: Entre Rios had a 16 days long HW in December 2013; this duration was greater than 95th percentile, a T. max greater than 90th percentile and a T. min higher than the 80th percentile when compared to HWs in the 1961-2010 period. In the case of Chaco, the HW had duration of 6 days, and it was not above the 80th percentile in any of the mentioned variables as compared to the HWs in the same province in the 1961-2010 period. The dates presented in the table are indicative since some provinces have more than one meteorological station and registered HWs do not necessarily coincide in the start and end dates. In those cases in which the HW did not affect the whole province, the analysis was carried out for the departments covered by the HW. In addition, when an HW was interrupted for a day but then the phenomenon continued, they were analyzed as a single event.

Table 1
Heat waves in all provinces and comparison with the historical (1961-2010) heat waves.

We thus identified 3 major periods under HWs in the warm semester 2013-14 in Argentina, in all cases extended in areas larger than one province. The HW in December 2013 reached 17 provinces with duration from 5 days in several provinces to 16 days in Entre Ríos. In January 2014 provinces under HWs were 16 although they were to some extent shorter, from 3 to 10 days. Different was the case of the HW in February which extended for 11 to 13 days over 4 northeastern provinces, close to the frontiers with Paraguay and Brazil. In consequence we found 37 provincial events of HWs in the time period studied.

It is interesting to note the different evolution of HWs in December and January. The first wave started in the eastern portion of the country, close to the Atlantic Ocean, then went to the west till the Andes, and later it went north coming back finally to the east. Instead, the second wave started by the Andes, then it traversed the continent to the east in direction to the Atlantic Ocean, north to the frontier with Paraguay and west again to the Andes.

When examining provincial results, we found 20 cases (9 in December, 7 in January and 4 in February) with extreme high (>95th percentile) mean duration, 11 cases (5, 4 and 2 respectively) with extreme high (>95th percentile) mean maximum temperature and 14 cases (4, 9 and 1, respectively) with extreme high (>95th percentile) mean minimum temperature than those historically recorded for HWs.

Provinces suffering the most severe impact of HWs (in terms of duration, maximum and minimum temperatures >P95th percentile, when compared with the 1961-2010 period) were La Rioja and Salta in December, La Pampa and Mendoza in January and Misiones in February.

Analysis of General Mortality

During the HWs that occurred in the summer of 2013-14 there were 8.759 deaths, 1.877 deaths in excess over the reference period for the total of the provinces analyzed (RR=1.23, 95%CI 1.20-1.28). Excess numbers of mortal cases were present in 32 of the 37 provincial events identified; only events in Jujuy, La Pampa and Tucumán in December and Chaco and Formosa in January did not present exceeding mortal cases.

General mortality rates and RR for each province under analysis are presented in Table 2. We found RR significantly increased during at least one HW in 13 out of the 18 provinces. In the provinces where the RR was statistically higher during the HWs than in its absence, the increase in the number of deaths accounted for 1.716 (1.28, 95%CI 1.24-1.33) deaths over those registered in the reference period.

Table 2
Provincial mortality rates and Relative Risk of Dying. Argentina, summer 2013-2014.

The HW of December 2013 affected 17 provinces, adding up an excess mortality of 1.046 (RR=1.28, 95%CI 1.23-1.34) deaths. The risk of dying was statistically higher in 7 out of the 15 jurisdictions that had increased mortality. The greatest increase in risk was recorded in Santiago del Estero (RR=1.65, 95%CI 1.21-2.26) where the maximum temperature was higher than 95th percentile, followed by CABA (RR=1.48, 95%CI 1.34-1.65) where the duration was very long (≥95th percentile), the average maximum temperature exceeded 90th percentile and the average minimum temperature exceeded 85th percentile (Table 1).

In January, the HW reached 16 jurisdictions, six of which presented significant increases in the risk of dying. The highest increase in mortality was registered in La Rioja (RR=2.3, 95%CI 1.48-3.58) where the event was extreme in temperatures (>95th percentile) but not in duration, followed by Tucumán (RR=1.36, 95%CI 1.12-1.64).

The HW in February, which was extreme in duration (>95th percentile), covered only the northeast of Argentina, with significant increases in the risk of dying in three out of the four provinces affected: Chaco (RR=1.34, 95%CI 1.18-1.53), Misiones (RR=1.25, 95%CI 1.10-1.43) and Corrientes (RR=1.22, 95%CI 1.08-1.39).

Mortality by sex

The specific mortality rates and the RR by sex were calculated for each province and HW and are presented in Table 3. To present a spatial view of these results, the RR maps for each of the HWs of summer 2013-2014 are presented in Figure 2.

Table 3
Relative risk for sex and age groups during heat waves of summer 2013-2014 in Argentina.

Figure 2
Relative risk maps for women (a, b, c) and men (d, e, f) during heat waves of summer 2013-2014 in Argentina.

A heterogeneous behaviour was observed between the provinces and even within the same province during different events. For example, the risk of dying in CABA during the HW of December was significantly high in both women and men (RR=1.58 and RR=1.34 respectively), while during the HW of January the increase was only significant in males (RR=1.33). Another general observation was that RRs did not paralleled mortality rates on sex separated calculus.

During the HW of December, the risk of dying for women increased significantly in six provinces with Santiago del Estero presenting the highest value (RR=1.67) followed by CABA (RR=1.58) and Entre Ríos (RR=1.53). During the HW of January the mortality of women did not show significantly increased risk for any of the 16 provinces involved, while in the HW of February, out of the four provinces involved only females from the province of Misiones presented a relative risk significantly increased (RR=1.36).

The risk of dying among males was significantly increased in five provinces during HWs in December, another five in January HW and in two out of the four provinces involved in the HW in February. In December, the greatest increases in the risk of dying of men were recorded in Santiago del Estero (RR=1.64), Mendoza (RR=1.36) and CABA (RR=1.34). In the HW of January, the risk of dying for males presented greater increases in La Rioja (RR=2.23), Tucumán (RR=1.57) and Santa Fe (RR=1.40). During the HW that took place in February the highest mortality rates for males and the RR significantly increased were found in the provinces of Chaco (RR=1.42) and Corrientes (RR=1.29).

Mortality by age groups

The groups 0 to 19 and 20 to 39 years old did not show increased risk of dying during HWs in none of the provinces.

Only one event showed a significant increase in the risk of dying in the group of 40 to 59 years old: it was in the province of Santiago del Estero during the December 2013 HW (RR=2.44), where 8 cases in the reference period turned into 21 cases during the HW.

In the age group of 60 to 79 years old the risk of dying was increased in at least one province in each of the three HWs during the summer 2013-14 (Table 3). Particularly, an increase in the risk of dying was found in CABA (RR=1.29) and Santa Fe (RR=1.23) during the HW of December 2013, in Mendoza (RR=1.33) and Santa Fe (RR=1.28) during the HW of January 2014 and in Formosa (RR=1.60) during the HW of February 2014. Spatial distribution of the results with significant increase in RR is presented in Figure 3.

Figure 3
Relative risk maps for age groups 60-79 years old (a, b, c) and over 79 years old (d, e, f) during heat waves of summer 2013-2014 in Argentina.

The age group over 79 years old had a significant increase in the risk of dying during the December HW in several provinces, with the highest RR value registered in Entre Ríos (RR=2.01), CABA (RR=1.67) and Santa Fe (RR=1.47), while in January 2014 HW significantly increased risk of dying for this age group was found in Tucumán (RR=1.58), Santa Fe (RR=1.41) and Córdoba (RR=1.28); and in February HW in Chaco (RR=1.50).

Mortality and causes of death

Globally speaking during the 266 days (adding up all days in all provinces) under HW in the summer 2013-2014 died 8.759 persons, expressing a risk increase of 1.23 (95%CI 1.20-1.28) over the value in absence of HW. Among them, registered causes of death that have epidemiological evidence of being related to high temperature accounted for 3.736 cases, being the more frequent heart diseases (I30-52 in ICD-10), ischemic heart disease (I20-25), cerebrovascular diseases (I60-69) and pneumonia (J12-18), which together accounted for 77,3% of the total.

Provincial mortality rates and RR were calculated to analyze the mortality linked to causes of death that have epidemiological evidence of being related to high temperature. Thus, statistically significant increases for those causes were found in 10 of the 37 provincial events (5 in December’s HW, 4 in January’s and 1 in February’s), in 6 out of the 18 provinces studied, accounting globally for 1.209 deaths, as presented in Table 4.

Table 4
Groups of morbidities with significantly increased risk during Heat Waves in Argentina, summer 2013-2014.

The most frequent causes of death with significantly increased RR were pneumonia and ischemic heart disease. Taking into consideration only those cases, pneumonia was responsible for 393 deaths in seven HW events (CABA, Entre Rios, Mendoza and Santa Fe in December, and CABA, Córdoba and Santa Fe in January), being the highest RR (3.18) the one in Mendoza in December; while ischemic heart disease, accounted for 394 cases in five provincial events with significantly increased RR (CABA and Córdoba in December, Córdoba and Mendoza in January and Misiones in February, this one with the highest RR=1.85).

On the other side, during the December 2013 HW, the 49% (481 cases) of all deaths that occurred in CABA were due to causes with heat linked epidemiological evidence, being lower respiratory disease the cause of death with the higher increase of RR (2.32), followed by renal insufficiency (RR=2.28), and pneumonia (RR=2.20). Similar picture was found in Santa Fe, where the higher increases were for lower respiratory disease (RR=2.02 in December and 2.18 in January), pneumonia (RR=2.03 in December and 1.71 in January) and diabetes mellitus (RR=1.99 in January). In Córdoba significative increases included ischemic heart disease (RR=1.83 in December and 1.84 in January), cerebrovascular diseases (RR=1.54 in December) and pneumonia (RR=1.58 in January). Other significative increases were found in Entre Rios (pneumonia in December RR=1.76), Mendoza (pneumonia in December RR=3.18; ischemic heart disease in January RR=1.82) and in Misiones in February (cerebrovascular diseases: RR=2.88; ischemic heart disease RR=1.85).

It is interesting to note that although “exposure to excessive natural heat” (Code X30 in the ICD-10) is a specific cause of death that can be assigned; we found that it is very rare. In our study only 12 deaths among the 8.759 registered during HWs were recorded as that: five cases in Chaco, three in Santa Fe, two in Entre Ríos, one in Corrientes and one in Córdoba; 50% cases in males, 9 out of them in people older than 60 years.

Discussion

In Argentina the frequency of HWs has increased north of parallel 40° South1515 Rusticucci M, Kysely J, Almeira G, Lhotka O. Long-term variability of heat waves in Argentina and recurrence probability of the severe 2008 heat wave in Buenos Aires. Theor Appl Climatol 2015; 124:679-689.. The summer 2013-2014 was particularly hot with three HWs that covered most of the country, which has been related with a growing number of days with minimal temperature above historical values rather than with changes in maximal temperatures2929 Centro Regional del Clima para el Sur de América del Sur (CRC-SAS). Climatología - Tendencias de Índices de Extremos Climáticos [Internet]. [acessado 2020 mar 20]. Disponível em: https://www.crc-sas.org/es/climatologia_extremos_climaticos.php.
https://www.crc-sas.org/es/climatologia_...
, although regional differences make it difficult to establish a homogeneous pattern.

It is also worth recalling that it has been suggested that the HW that occurred in December 2013 in central and northern Argentina was partially related with anthropogenic forcing; observing that the recurrence time of this event would be 75 years without anthropogenic forcing and 15 years considering such forcing3030 Hannart A, Vera C, Otto FE, Cerne B. Causal influence of anthropogenic forcings on the Argentinean heat wave of December 2013. Special Suppl Bull Am Meteorol Soc 2015; 96:541-545..

Overall mortality during the HWs of the summer of 2013-2014 showed heterogeneous behaviour in provinces, yet increases in absolute and relative values were observed in relation to the three previous summer’s reference period. In absolute terms there were 1.877 excess deaths, as well as significantly increased risk of dying in 11 of the 18 provinces reached by HWs and in 15 out of the 37 provincial events. Interestingly, the four most populated provinces (CABA, Córdoba, Santa Fe and Mendoza) did presented elevated RR during December and January HWs while, on the other hand, five much less populated provinces (Catamarca, Formosa, Jujuy, La Pampa and San Luis) did not increased their RR during any HW; which could be linked to the fact that in small populations, small changes in the number of deaths can cause great variability in rates and RR, leading to a loss of statistical significance. Nevertheless, it is not clear yet why the impact on mortality varies among provinces. For example, provincial events which lasted about the same, with similar severity in terms of temperature intensity, did show excess numbers of mortal cases or not (San Luis vs La Pampa in December; Tucuman vs. Chaco in January). These differences have been observed in similar studies44 Whitman S, Good G, Donoghue ER, Benbow N, Shou W, Mou S. Mortality in Chicago Attributed to the July 1995 Heat Wave. Am J Public Health 1997; 87:1515-1518.,66 Oudin Åström D, Forsberg B, Rocklöv J. Heat wave impact on morbidity and mortality in the elderly population: A review of recent studies. Maturitas 2011; 69:99-105.,77 Shaposhnikov D, Revich B, Bellander T, Bedada GB, Bottai M, Kharkova T, Kvasha E, Lezina E, Lind T, Semutnikova E, Pershagen G. Mortality Related to Air Pollution with the Moscow Heat Wave and Wildfire of 2010. Epidemiology 2014; 25:359-364.,3131 World Meteorological Organization. World Health Organization (WHO). Atlas del Clima y la Salud. Ginebra: WHO; 2012. but the values are difficult to compare due to differences in the scales of phenomena (spatial, temporal and magnitude of HWs). In order to understand these differences new studies at local scale are necessary.

Relative risk of dying could not be undoubtefully linked to any single variable (length of the HW, maximal or minimal temperatures reached, total population, crude number or rate of death in control period). However, higher RR showed a tendency to accompany provinces with bigger cities, as shown in Table 2.

Health impact of a HW could be linked to its intensity, either in terms of temperature and/or the duration, even though differences in baseline reference values do make it difficult to compare phenomena taking place in different locations for a number of reasons. Among them we include differences in local mortality rates, the population socioeconomic profiles - since impoverished people have fewer resources to face hot climate -, importance of activities such as construction and agriculture, based on outdoor tasks, as well as informal work - usually developed in precarious conditions - which naturally generate greater exposure to climatic conditions33 McGregor GR, Bessemoulin P, Ebi K, Menne B. Heatwaves and Health: Guidance on Warning-System Development. Geneva: WHO; 2015.,88 D'Ippoliti D, Michelozzi P, Marino C, Donato F, Menne B, Katsouyanni K, Kirchmayer U, Analitis A, Medina-Ramón M, Paldy A, Atkinson R, Kovats S, Bisanti L, Schneider A, Lefranc A, Iñiguez C, Perucci CA. The impact of heat waves on mortality in 9 European cities: results from the EuroHEAT project. Environ Health 2010; 9:37.

9 McGeehin MA, Mirabelli M. The Potential Impacts of Climate Variability and Change on Temperature-Related Morbidity and Mortality in the United States. Environ Health Persp 2001; 109:185-189.
-1010 Smith KR, Woodward A, Campbell-Lendrum D, Chadee DD, Honda Y, Liu Q, Olwoch JM, Revich B, Sauerborn R, Aranda C, Berry H, Butler C, Chafe Z, Cushing L, Ebi KL, Kjellstrom T, Kovats S, Lindsay G, Lipp E, McMichael T, Murray V, Sankoh O, O'Neill M, Shonkoff SB, Sutherland J, Yamamoto S. Human Health: Impacts, Adaptation, and Co-Benefits. In: Field CB, Barros VR, Dokken DJ, Mach KJ, Mastrandrea MD, Bilir TE, Chatterjee M, Ebi KL, Estrada YO, Genova RC, Girma B, Kissel ES, Levy AN, MacCracken S, Mastrandrea PR, White LL, editores. Climate Change 2014: Impacts, Adaptation and Vulnerability. Part A: Global and Sectorial Aspects. Contribution of Working Group II to the Fifth Assessment Report of the Intergovernmental Panel on Climate Change. Cambridge UK: Cambridge University Press; 2014. p. 709-754.,3232 McGregor GR, Pelling M, Tanja Wolf T, Gosling S. United Kingdom Environment Agency Science Report [Internet]. Bristol; 2007 [cited 2020 mar 20]. Available from: www.environment-agency.gov.uk.
www.environment-agency.gov.uk...

33 Son J-Y, Gouveia N, Bravo MA, Freitas CU, Bell MK. The impact of temperature on mortality in a subtropical city: Effects of cold, heat, and heat waves in São Paulo, Brazil. Int J Biometeorol 2016; 60:113-121.
-3434 Guo Y, Gasparrini A, Armstrong BG, Tawatsupa T, Tobias A, Lavigne E, Coelho MSZS, Pan X, Kim H, Hashizume M, Honda Y, Guo YL, Wu CF, Zanobetti A, Schwartz JD, Bell ML, Scortichini M, Michelozzi P, Punnasiri K, Li S, Tian L, Garcia SDO, Seposo X, Overcenco A, Zeka A, Goodman P, Dang TN, Dung DV, Mayvaneh F, Saldiva PHN, Williams G, Tong S. Heat Wave and Mortality: A Multicountry, Multicommunity Study. Environ Health Perspect 2017; 125(8):087006..

The importance of acclimatization of individuals to the local conditions in which they live3434 Guo Y, Gasparrini A, Armstrong BG, Tawatsupa T, Tobias A, Lavigne E, Coelho MSZS, Pan X, Kim H, Hashizume M, Honda Y, Guo YL, Wu CF, Zanobetti A, Schwartz JD, Bell ML, Scortichini M, Michelozzi P, Punnasiri K, Li S, Tian L, Garcia SDO, Seposo X, Overcenco A, Zeka A, Goodman P, Dang TN, Dung DV, Mayvaneh F, Saldiva PHN, Williams G, Tong S. Heat Wave and Mortality: A Multicountry, Multicommunity Study. Environ Health Perspect 2017; 125(8):087006.

35 Curriero FC, Heiner KS, Samet JM, Zeger SL, Patz JA. Temperature mortality in 11 cities of the Eastern of the United States. Am J Epidemiol 2002; 155(1):80-87.
-3636 Díaz J, Jordán A, García R, López C, Alberdi JC, Hernández E, Otero A. Heat waves in Madrid 1986-1997: effects on the health of the elderly. Int Arch Occup Environ Health 2002; 75(3):163-170. is supported by the fact that the temperature of triggering mortality attributable to heat varies from place to place, and indicates that each person is accustomed to living in a range of daily temperatures that is below some threshold.

Other studies suggest that HWs that occurred early in the warm season may have greater effects on the health of the population due to lack of acclimatization1010 Smith KR, Woodward A, Campbell-Lendrum D, Chadee DD, Honda Y, Liu Q, Olwoch JM, Revich B, Sauerborn R, Aranda C, Berry H, Butler C, Chafe Z, Cushing L, Ebi KL, Kjellstrom T, Kovats S, Lindsay G, Lipp E, McMichael T, Murray V, Sankoh O, O'Neill M, Shonkoff SB, Sutherland J, Yamamoto S. Human Health: Impacts, Adaptation, and Co-Benefits. In: Field CB, Barros VR, Dokken DJ, Mach KJ, Mastrandrea MD, Bilir TE, Chatterjee M, Ebi KL, Estrada YO, Genova RC, Girma B, Kissel ES, Levy AN, MacCracken S, Mastrandrea PR, White LL, editores. Climate Change 2014: Impacts, Adaptation and Vulnerability. Part A: Global and Sectorial Aspects. Contribution of Working Group II to the Fifth Assessment Report of the Intergovernmental Panel on Climate Change. Cambridge UK: Cambridge University Press; 2014. p. 709-754.,3636 Díaz J, Jordán A, García R, López C, Alberdi JC, Hernández E, Otero A. Heat waves in Madrid 1986-1997: effects on the health of the elderly. Int Arch Occup Environ Health 2002; 75(3):163-170.,3737 Anderson GB, Bell ML. Heat Waves in the United States: Mortality Risk during Heat Waves and Effect Modification by Heat Wave Characteristics in 43 U.S. Communities. Environ Health Perspect 2011; 119(2):210-218., since human physical adaptation to high temperatures can take between two and six weeks3838 Guyton AC, Hall JE. Textbook of Medical Physiology. Philadelphia: Elsevier Saunders; 2006.. This was not the case here; even when 55.7% of excess deaths occurred during the December HW (1.046 deaths, mortality rate: 4.65 per 100.000 population), the February HW affecting the NEA caused a total of 196 excess deaths (mortality rate: 4.98 per 100.000 population) and recorded significant increases in the risk of dying in three out of the four provinces reached for the event (RR=1.24, 95%CI 1.12-1.36) an event that occurred late in the summer (January 31 to February 12) and in provinces where high temperatures are frequent.

When death of males and females were analyzed separately, HW impact was heterogeneous on both temporal and spatial scales. In the December HW four provinces (CABA, Córdoba, Santa Fe and Santiago del Estero) showed increased risk of dying for both groups, as well as in two other provinces for females and in one province for males. During the January HW the behaviour was very dissimilar: RR in women did not increase significantly in any province while RR in males showed significant increases in 5 provinces. February HW led to increased RR for women in one province (Misiones) and for men in other two provinces (Chaco and Corrientes) out of the four compromised. As a whole for our study, RR>1 was found for men in 12 provincial events and for women in 7 events. Changes in the risk of dying of women between December and January could be explained by a harvesting effect, although that effect was not verified in a previous study in CABA2222 Chesini F, Abrutzky R, de Titto EH. Mortalidad por olas de calor en la Ciudad de Buenos Aires, Argentina (2005-2015). Cad Saude Publica 2019; 35:e00165218..

Thus, no clear association between sex and mortality rates or RR of dying could be drawn. In some events significant increases included both sexes or just one of them, provinces did not keep any pattern even from one HW to the following one, as previously described for several Spanish3636 Díaz J, Jordán A, García R, López C, Alberdi JC, Hernández E, Otero A. Heat waves in Madrid 1986-1997: effects on the health of the elderly. Int Arch Occup Environ Health 2002; 75(3):163-170. and Latin American39 cities but differently from other studies performed in Europe88 D'Ippoliti D, Michelozzi P, Marino C, Donato F, Menne B, Katsouyanni K, Kirchmayer U, Analitis A, Medina-Ramón M, Paldy A, Atkinson R, Kovats S, Bisanti L, Schneider A, Lefranc A, Iñiguez C, Perucci CA. The impact of heat waves on mortality in 9 European cities: results from the EuroHEAT project. Environ Health 2010; 9:37.,1010 Smith KR, Woodward A, Campbell-Lendrum D, Chadee DD, Honda Y, Liu Q, Olwoch JM, Revich B, Sauerborn R, Aranda C, Berry H, Butler C, Chafe Z, Cushing L, Ebi KL, Kjellstrom T, Kovats S, Lindsay G, Lipp E, McMichael T, Murray V, Sankoh O, O'Neill M, Shonkoff SB, Sutherland J, Yamamoto S. Human Health: Impacts, Adaptation, and Co-Benefits. In: Field CB, Barros VR, Dokken DJ, Mach KJ, Mastrandrea MD, Bilir TE, Chatterjee M, Ebi KL, Estrada YO, Genova RC, Girma B, Kissel ES, Levy AN, MacCracken S, Mastrandrea PR, White LL, editores. Climate Change 2014: Impacts, Adaptation and Vulnerability. Part A: Global and Sectorial Aspects. Contribution of Working Group II to the Fifth Assessment Report of the Intergovernmental Panel on Climate Change. Cambridge UK: Cambridge University Press; 2014. p. 709-754.,4040 Ishigami A, Hajat S, Kovats RS, Bisanti L, Rognoni M, Russo A, Paldy A. An ecological time-series study of heat-related mortality in three European cities. Environ Health 2008; 7:5. and Brazil4141 Geirinhas JL, Trigo RM, Libonati R, Castro LCO, Sousa PM, Coelho CAS, Peres LF, Magalhães FMA. Characterizing the atmospheric conditions during the 2010 heatwave in Rio de Janeiro marked by excessive mortality rates. Sci Total Environ 2019; 650:796-808.,4242 Geirinhas JL, Russo A, Libonati R, Trigo RM, Castro LCO, Peres LF, Magalhães MAFM, Nunes B. Heat-related mortality at the beginning of the twenty-first century in Rio de Janeiro, Brazil. Int J Biometeor 2020; 64(8):1319-1332. suggesting that women are more sensitive than males to HW.

The literature agrees that older adults are at higher risk during HWs44 Whitman S, Good G, Donoghue ER, Benbow N, Shou W, Mou S. Mortality in Chicago Attributed to the July 1995 Heat Wave. Am J Public Health 1997; 87:1515-1518.,1010 Smith KR, Woodward A, Campbell-Lendrum D, Chadee DD, Honda Y, Liu Q, Olwoch JM, Revich B, Sauerborn R, Aranda C, Berry H, Butler C, Chafe Z, Cushing L, Ebi KL, Kjellstrom T, Kovats S, Lindsay G, Lipp E, McMichael T, Murray V, Sankoh O, O'Neill M, Shonkoff SB, Sutherland J, Yamamoto S. Human Health: Impacts, Adaptation, and Co-Benefits. In: Field CB, Barros VR, Dokken DJ, Mach KJ, Mastrandrea MD, Bilir TE, Chatterjee M, Ebi KL, Estrada YO, Genova RC, Girma B, Kissel ES, Levy AN, MacCracken S, Mastrandrea PR, White LL, editores. Climate Change 2014: Impacts, Adaptation and Vulnerability. Part A: Global and Sectorial Aspects. Contribution of Working Group II to the Fifth Assessment Report of the Intergovernmental Panel on Climate Change. Cambridge UK: Cambridge University Press; 2014. p. 709-754.,3232 McGregor GR, Pelling M, Tanja Wolf T, Gosling S. United Kingdom Environment Agency Science Report [Internet]. Bristol; 2007 [cited 2020 mar 20]. Available from: www.environment-agency.gov.uk.
www.environment-agency.gov.uk...
,4343 Conti S, Meli P, Minelli G, Solimini R, Toccaceli V, Vichi M, Beltrano C, Perini L. Epidemiologic study of mortality during the Summer 2003 heat wave in Italy. Environ Res 2005; 98:390-399.

44 Martínez Navarro F, Simón-Soria F, López-Abente G. Valoración del impacto de la ola de calor del verano de 2003 sobre la mortalidad. Gac Sanit 2004; 18(Supl. 1):250-258.
-4545 Linares C, Díaz J. Temperaturas extremadamente elevadas y su impacto sobre la mortalidad diaria según diferentes grupos de edad. Gac Sanit 2008; 22:115-119.. The present study found the risk of dying significantly increased in four jurisdictions for the 60-79 years old group and in six jurisdictions for those 80 years old and over, taking together all HWs of the summer 2013-2014. Only one province presented a significant increase in the risk of dying in the group of 40-59 years old during the December 2013 HW and none province presented differences for lower ages.

All causes of death that had a significant increase in risk during HWs have epidemiological evidence and biological plausibility88 D'Ippoliti D, Michelozzi P, Marino C, Donato F, Menne B, Katsouyanni K, Kirchmayer U, Analitis A, Medina-Ramón M, Paldy A, Atkinson R, Kovats S, Bisanti L, Schneider A, Lefranc A, Iñiguez C, Perucci CA. The impact of heat waves on mortality in 9 European cities: results from the EuroHEAT project. Environ Health 2010; 9:37.

9 McGeehin MA, Mirabelli M. The Potential Impacts of Climate Variability and Change on Temperature-Related Morbidity and Mortality in the United States. Environ Health Persp 2001; 109:185-189.
-1010 Smith KR, Woodward A, Campbell-Lendrum D, Chadee DD, Honda Y, Liu Q, Olwoch JM, Revich B, Sauerborn R, Aranda C, Berry H, Butler C, Chafe Z, Cushing L, Ebi KL, Kjellstrom T, Kovats S, Lindsay G, Lipp E, McMichael T, Murray V, Sankoh O, O'Neill M, Shonkoff SB, Sutherland J, Yamamoto S. Human Health: Impacts, Adaptation, and Co-Benefits. In: Field CB, Barros VR, Dokken DJ, Mach KJ, Mastrandrea MD, Bilir TE, Chatterjee M, Ebi KL, Estrada YO, Genova RC, Girma B, Kissel ES, Levy AN, MacCracken S, Mastrandrea PR, White LL, editores. Climate Change 2014: Impacts, Adaptation and Vulnerability. Part A: Global and Sectorial Aspects. Contribution of Working Group II to the Fifth Assessment Report of the Intergovernmental Panel on Climate Change. Cambridge UK: Cambridge University Press; 2014. p. 709-754.. Exposure to extreme heat conditions for long periods can exacerbate a wide range of pre-existing conditions by triggering the chain of pathological events that lead directly to death.

The highest RR in causes of death was recorded for pneumonia; which also had the highest frequency (seven jurisdictions with significantly increased risk). The accumulated cases of respiratory diseases (pneumonia and lower respiratory diseases) were responsible for 457 deaths, only in those jurisdictions where the risk was significantly increased. The incidence of HWs on respiratory disease mortality is consistent with other studies88 D'Ippoliti D, Michelozzi P, Marino C, Donato F, Menne B, Katsouyanni K, Kirchmayer U, Analitis A, Medina-Ramón M, Paldy A, Atkinson R, Kovats S, Bisanti L, Schneider A, Lefranc A, Iñiguez C, Perucci CA. The impact of heat waves on mortality in 9 European cities: results from the EuroHEAT project. Environ Health 2010; 9:37.. Cardiovascular disease accounted for 48% of the total causes of death at significantly increased risk observed in five jurisdictions. From these analyses, it would seem that public health programs to prevent heat-related mortality would best be directed at the elderly and those persons with cardiovascular and respiratory diseases.

Despite that the literature33 McGregor GR, Bessemoulin P, Ebi K, Menne B. Heatwaves and Health: Guidance on Warning-System Development. Geneva: WHO; 2015.,1111 Montero JC, Mirón IJ, Criado JJ, Linares C, Díaz J. Difficulties of defining the term "heat wave" in public health. Int J Environ Health Res 2013; 23:377-379. agrees that living in large cities and different socioeconomic conditions can modify risk levels we did not observe a clear pattern, at least in this particular summer. Perhaps, the analysis unit selected (provinces) is not sensitive enough to analyze the influence of living in large cities.

In conclusion, this study analyzed a particularly warm summer in which three HWs of long duration and intensity were recorded. The characterization of mortality during HWs revealed increases in both absolute and relative values in overall mortality and mortality by sex, age groups and major causes of death heterogeneously along the country. This information was useful in reviewing and expanding the Early Warning System for Heat and Health Waves (SAT-OCS, for its Spanish denomination) in Argentina, implemented operationally in the whole country on the summer 2017-20184646 Herrera N, Skansi MM, Berón M, Campetella C, Cejas A, Chasco J, Chesini F, de Titto E, Gatto M, Saucedo M, Suaya M. Sistema de Alerta Temprana por Olas de Calor y Salud (SAT-OCS). Technical note SMN 2018-50 [Internet]. 2018 [cited 2020 mar 20]. Available from: http://repositorio.smn.gob.ar/handle/20.500.12160/772.
http://repositorio.smn.gob.ar/handle/20....
. This system operates in an automated manner for 57 locations of the country, issues a daily alert valid for 24 hours and is available online from the first day of October each year to 31 March, with the aim of enabling both the population and civil protection agencies to take appropriate prevention, mitigation and response measures at each alert level.

References

  • 1
    World Health Organization (WHO). 61ª Asamblea Mundial de la Salud. Cambio climático y salud [Internet]. 2008 [cited 2020 mar 19]. Available from: http://apps.who.int/iris/bitstream/10665/26288/1/A61_R19-sp.pdf
    » http://apps.who.int/iris/bitstream/10665/26288/1/A61_R19-sp.pdf
  • 2
    Mercado Común del Sur. Resolución MERCOSUR nº 12/2009. Estrategia Cambio Climático y Salud. Montevideo [Internet]. 2009 [cited 2020 mar 19]. Available from: http://200.40.51.218/SAM/GestDoc/PubWeb.nsf/Reunion?OpenAgent&id=CD00D62A0920FFD7032576A50049C42A&lang=ESP
    » http://200.40.51.218/SAM/GestDoc/PubWeb.nsf/Reunion?OpenAgent&id=CD00D62A0920FFD7032576A50049C42A&lang=ESP
  • 3
    McGregor GR, Bessemoulin P, Ebi K, Menne B. Heatwaves and Health: Guidance on Warning-System Development. Geneva: WHO; 2015.
  • 4
    Whitman S, Good G, Donoghue ER, Benbow N, Shou W, Mou S. Mortality in Chicago Attributed to the July 1995 Heat Wave. Am J Public Health 1997; 87:1515-1518.
  • 5
    Robine JM, Cheung SL, Le Roy S, Van Oyen H, Griffiths C, Michel JP, Herrmann FR. Death toll exceeded 70,000 in Europe during the summer of 2003. C R Biol 2008; 331:171-178.
  • 6
    Oudin Åström D, Forsberg B, Rocklöv J. Heat wave impact on morbidity and mortality in the elderly population: A review of recent studies. Maturitas 2011; 69:99-105.
  • 7
    Shaposhnikov D, Revich B, Bellander T, Bedada GB, Bottai M, Kharkova T, Kvasha E, Lezina E, Lind T, Semutnikova E, Pershagen G. Mortality Related to Air Pollution with the Moscow Heat Wave and Wildfire of 2010. Epidemiology 2014; 25:359-364.
  • 8
    D'Ippoliti D, Michelozzi P, Marino C, Donato F, Menne B, Katsouyanni K, Kirchmayer U, Analitis A, Medina-Ramón M, Paldy A, Atkinson R, Kovats S, Bisanti L, Schneider A, Lefranc A, Iñiguez C, Perucci CA. The impact of heat waves on mortality in 9 European cities: results from the EuroHEAT project. Environ Health 2010; 9:37.
  • 9
    McGeehin MA, Mirabelli M. The Potential Impacts of Climate Variability and Change on Temperature-Related Morbidity and Mortality in the United States. Environ Health Persp 2001; 109:185-189.
  • 10
    Smith KR, Woodward A, Campbell-Lendrum D, Chadee DD, Honda Y, Liu Q, Olwoch JM, Revich B, Sauerborn R, Aranda C, Berry H, Butler C, Chafe Z, Cushing L, Ebi KL, Kjellstrom T, Kovats S, Lindsay G, Lipp E, McMichael T, Murray V, Sankoh O, O'Neill M, Shonkoff SB, Sutherland J, Yamamoto S. Human Health: Impacts, Adaptation, and Co-Benefits. In: Field CB, Barros VR, Dokken DJ, Mach KJ, Mastrandrea MD, Bilir TE, Chatterjee M, Ebi KL, Estrada YO, Genova RC, Girma B, Kissel ES, Levy AN, MacCracken S, Mastrandrea PR, White LL, editores. Climate Change 2014: Impacts, Adaptation and Vulnerability. Part A: Global and Sectorial Aspects. Contribution of Working Group II to the Fifth Assessment Report of the Intergovernmental Panel on Climate Change. Cambridge UK: Cambridge University Press; 2014. p. 709-754.
  • 11
    Montero JC, Mirón IJ, Criado JJ, Linares C, Díaz J. Difficulties of defining the term "heat wave" in public health. Int J Environ Health Res 2013; 23:377-379.
  • 12
    Perkins SE. A review on the scientific understanding of heatwaves. Their measurement, driving mechanisms, and changes at the global scale. Atmospheric Res 2015; 164-165:242-267
  • 13
    Servicio Meteorológico Nacional (SMN). Argentina. Datos - Clima de Argentina-Olas de calor (Heat waves) [Internet]. [cited 2020 mar 19]. Available from: https://www.smn.gob.ar/estadisticas
    » https://www.smn.gob.ar/estadisticas
  • 14
    Centro de Investigaciones del Mar y la Atmósfera (CIMA). Modelos Climáticos: Resumen Ejecutivo (Climate Models: Executive Summary). Secretaría de Ambiente y Desarrollo Sustentable [Internet]. 2015 [cited 2020 mar 20]. Available from: http://3cn.cima.fcen.uba.ar/3cn_informe.php
    » http://3cn.cima.fcen.uba.ar/3cn_informe.php
  • 15
    Rusticucci M, Kysely J, Almeira G, Lhotka O. Long-term variability of heat waves in Argentina and recurrence probability of the severe 2008 heat wave in Buenos Aires. Theor Appl Climatol 2015; 124:679-689.
  • 16
    Cerne SB, Vera CS. Influence of the intraseasonal variability on heat waves in subtropical South America. Clim Dyn 2011; 36:2265-2277.
  • 17
    Geirinhas JL, Russo A, Libonati R, Sousa PM, Miralles DG, Trigo RM. Recent increasing frequency of compound summer drought and heatwaves in Southeast Brazil. Environ Res Lett 2021; 16:034036.
  • 18
    Cerne B, Vera C, Álvarez M. Influencia de la variabilidad intraestacional en el desarrollo de la ola de calor de diciembre de 2013 en el este de Sudamérica. In: Congreso Argentino de Meteorología (CONGREMET) XII. Mar del Plata; 2015 maio 26-29.
  • 19
    Osman M, Álvarez M. Subseasonal prediction of the heat wave of December 2013 in Southern South America by the POAMA and BCC CPS models. Climate Dynamics 2018; 50(1-2):67-81.
  • 20
    Abrutzky R, Dawidowski L, Matus P, Romero-Lankao P. Health effects on climate and air pollution in Buenos Aires. A first series analysis. J Environ Prot 2012; 3:262-271.
  • 21
    Almeira G, Rusticucci M, Suaya M. Relación entre mortalidad y temperaturas extremas en Buenos Aires y Rosario. Meteorológica 2016; 41:65-79.
  • 22
    Chesini F, Abrutzky R, de Titto EH. Mortalidad por olas de calor en la Ciudad de Buenos Aires, Argentina (2005-2015). Cad Saude Publica 2019; 35:e00165218.
  • 23
    Instituto Nacional de Estadística y Censos (INDEC). Argentina. Censo Nacional 2010 de Población, Hogares y Viviendas [Internet]. [cited 2020 mar 18]. Available from: http://www.indec.gov.ar/index.asp
    » http://www.indec.gov.ar/index.asp
  • 24
    Dirección Nacional de Estadísticas e Información de la Salud (DEIS). Indicadores Básicos Argentina 2014 [Internet]. Buenos Aires: Ministerio de Salud de la Nación; 2015 [cited 2020 mar 20]. Available from: http://www.deis.gov.ar/publicaciones/archivos/indicadores_2014.pdf
    » http://www.deis.gov.ar/publicaciones/archivos/indicadores_2014.pdf
  • 25
    Culqui DR, Díaz J, Simón F, Linares C. Análisis de los impactos de las olas de calor sobre la mortalidad de la ciudad de Madrid durante el período 1990-2009. Rev Esp Salud Pública 2013; 87:277-282.
  • 26
    Beaglehole R, Bonita R, Kjellström T. Epidemiologia Básica. Washington D.C.: PAHO; 1994.
  • 27
    Carracedo-Martínez E, Tobías A, Saez M, Taracido M, Figueras A. Fundamentos y aplicaciones del diseño de casos cruzados. Gac Sanit 2009; 23:161-165.
  • 28
    Martín JF. Los factores definitorios de los grandes grupos de edad de la población: tipos, subgrupos y umbrales. Scripta Nova 2005; IX:190.
  • 29
    Centro Regional del Clima para el Sur de América del Sur (CRC-SAS). Climatología - Tendencias de Índices de Extremos Climáticos [Internet]. [acessado 2020 mar 20]. Disponível em: https://www.crc-sas.org/es/climatologia_extremos_climaticos.php
    » https://www.crc-sas.org/es/climatologia_extremos_climaticos.php
  • 30
    Hannart A, Vera C, Otto FE, Cerne B. Causal influence of anthropogenic forcings on the Argentinean heat wave of December 2013. Special Suppl Bull Am Meteorol Soc 2015; 96:541-545.
  • 31
    World Meteorological Organization. World Health Organization (WHO). Atlas del Clima y la Salud. Ginebra: WHO; 2012.
  • 32
    McGregor GR, Pelling M, Tanja Wolf T, Gosling S. United Kingdom Environment Agency Science Report [Internet]. Bristol; 2007 [cited 2020 mar 20]. Available from: www.environment-agency.gov.uk.
    » www.environment-agency.gov.uk
  • 33
    Son J-Y, Gouveia N, Bravo MA, Freitas CU, Bell MK. The impact of temperature on mortality in a subtropical city: Effects of cold, heat, and heat waves in São Paulo, Brazil. Int J Biometeorol 2016; 60:113-121.
  • 34
    Guo Y, Gasparrini A, Armstrong BG, Tawatsupa T, Tobias A, Lavigne E, Coelho MSZS, Pan X, Kim H, Hashizume M, Honda Y, Guo YL, Wu CF, Zanobetti A, Schwartz JD, Bell ML, Scortichini M, Michelozzi P, Punnasiri K, Li S, Tian L, Garcia SDO, Seposo X, Overcenco A, Zeka A, Goodman P, Dang TN, Dung DV, Mayvaneh F, Saldiva PHN, Williams G, Tong S. Heat Wave and Mortality: A Multicountry, Multicommunity Study. Environ Health Perspect 2017; 125(8):087006.
  • 35
    Curriero FC, Heiner KS, Samet JM, Zeger SL, Patz JA. Temperature mortality in 11 cities of the Eastern of the United States. Am J Epidemiol 2002; 155(1):80-87.
  • 36
    Díaz J, Jordán A, García R, López C, Alberdi JC, Hernández E, Otero A. Heat waves in Madrid 1986-1997: effects on the health of the elderly. Int Arch Occup Environ Health 2002; 75(3):163-170.
  • 37
    Anderson GB, Bell ML. Heat Waves in the United States: Mortality Risk during Heat Waves and Effect Modification by Heat Wave Characteristics in 43 U.S. Communities. Environ Health Perspect 2011; 119(2):210-218.
  • 38
    Guyton AC, Hall JE. Textbook of Medical Physiology. Philadelphia: Elsevier Saunders; 2006.
  • 39
    Bell ML, O'Neill MS, Ranjit N, Borja-Aburto VH, Cifuentes LA, Gouveia NC. Vulnerability to heat-related mortality in Latin America: a case-crossover study in São Paulo, Brazil, Santiago, Chile and Mexico City, Mexico. Int J Epidemiol 2008; 37(4):796-804.
  • 40
    Ishigami A, Hajat S, Kovats RS, Bisanti L, Rognoni M, Russo A, Paldy A. An ecological time-series study of heat-related mortality in three European cities. Environ Health 2008; 7:5.
  • 41
    Geirinhas JL, Trigo RM, Libonati R, Castro LCO, Sousa PM, Coelho CAS, Peres LF, Magalhães FMA. Characterizing the atmospheric conditions during the 2010 heatwave in Rio de Janeiro marked by excessive mortality rates. Sci Total Environ 2019; 650:796-808.
  • 42
    Geirinhas JL, Russo A, Libonati R, Trigo RM, Castro LCO, Peres LF, Magalhães MAFM, Nunes B. Heat-related mortality at the beginning of the twenty-first century in Rio de Janeiro, Brazil. Int J Biometeor 2020; 64(8):1319-1332.
  • 43
    Conti S, Meli P, Minelli G, Solimini R, Toccaceli V, Vichi M, Beltrano C, Perini L. Epidemiologic study of mortality during the Summer 2003 heat wave in Italy. Environ Res 2005; 98:390-399.
  • 44
    Martínez Navarro F, Simón-Soria F, López-Abente G. Valoración del impacto de la ola de calor del verano de 2003 sobre la mortalidad. Gac Sanit 2004; 18(Supl. 1):250-258.
  • 45
    Linares C, Díaz J. Temperaturas extremadamente elevadas y su impacto sobre la mortalidad diaria según diferentes grupos de edad. Gac Sanit 2008; 22:115-119.
  • 46
    Herrera N, Skansi MM, Berón M, Campetella C, Cejas A, Chasco J, Chesini F, de Titto E, Gatto M, Saucedo M, Suaya M. Sistema de Alerta Temprana por Olas de Calor y Salud (SAT-OCS). Technical note SMN 2018-50 [Internet]. 2018 [cited 2020 mar 20]. Available from: http://repositorio.smn.gob.ar/handle/20.500.12160/772
    » http://repositorio.smn.gob.ar/handle/20.500.12160/772

Publication Dates

  • Publication in this collection
    04 May 2022
  • Date of issue
    May 2022

History

  • Received
    30 Nov 2020
  • Accepted
    29 July 2021
  • Published
    01 Aug 2021
ABRASCO - Associação Brasileira de Saúde Coletiva Rio de Janeiro - RJ - Brazil
E-mail: revscol@fiocruz.br