Heat Wave Events over Georgia since 1961: Climatology, Changes and Severity

The Caucasus Region has been affected by an increasing number of heat waves during the last decades, which have had serious impacts on human health, agriculture and natural ecosystems. A dataset of 22 homogenized, daily maximum (Tmax) and minimum (Tmin) air temperature series is developed to quantify climatology and summer heat wave changes for Georgia and Tbilisi station between 1961 and 2010 using the extreme heat factor (EHF) as heat wave index. The EHF is studied with respect to eight heat wave aspects: event number, duration, participating heat wave days, peak and mean magnitude, number of heat wave days, severe and extreme heat wave days. A severity threshold for each station was determined by the climatological distribution of heat wave intensity. Moreover, heat wave series of two indices focusing on the 90th percentile of daily minimum temperature (CTN90p) and the 90th percentile of daily maximum temperature (CTX90p) were compared. The spatial distribution of heat wave characteristics over Georgia showed a concentration of high heat wave amplitudes and mean magnitudes in the Southwest. The longest and most frequently occurring heat wave events were observed in the Southeast of Georgia. Most severe heat wave events were found in both regions. Regarding the monthly distribution of heat waves, the largest proportion of severe events and highest intensities are measured during May. Trends for all Georgia-averaged heat wave aspects demonstrate significant increases in the number, intensity and duration of low-and high-intensity heat waves. However, for the heat wave mean magnitude no change was observed. Heat wave trend 309 magnitudes for Tbilisi mainly exceed the Georgia-averages and its surrounding stations, implying urban heat island (UHI) effects and synergistic interactions between heat waves and UHIs. Comparing heat wave aspects for CTN90p and CTX90p, all trend magnitudes for CTN90p were larger, while the correlation between the annual time-series was very high among all heat wave indices analyzed. This finding reflects the importance of integrating the most suitable heat wave index into a sector-specific impact analysis.


Introduction
Heat waves are among the most threatening meteorological hazards related to global warming posing impacts to society, economy and ecology.Effects of heat waves can include an increase in morbidity and mortality rates (e.g., [1][2][3][4][5]), a rising stress on agricultural resources (e.g., [6,7]) and a strain on infrastructure (e.g., [8]).According to the Intergovernmental Panel on Climate Change (IPCC) anthropogenic influences on climate since the mid-20th century resulted in a change of frequency and intensity of daily temperature extremes and doubled the probability of occurrence of heat waves in some regions of the world [9] and the global averaged frequency, duration and intensity of heat waves are projected to increase [10][11][12][13].The severe heat waves during the years 2001, 2003, 2006 and 2007 in Europe, and Russia in 2010, and their impacts have been investigated in various studies [14][15][16][17][18][19][20][21].The fact that in West Asia heat waves are likely to increase in frequency and/or duration, makes it even more important to investigate observed heat wave changes with special regards to regional impacts, the communities involved, and the climatic fields affected [13].Georgia is located in the Southern Caucasus Region between 41°-44°N and 40°-47°E at the border of Europe and Asia.Although it covers a small area of 69,700 km 2 , the country's physiographic and climatic conditions are very diverse.The relief declines from the Greater Caucasus Range in the North and the Southern Caucasus towards Transcaucasia, which stretches from the Black Sea coast in the West to the Eastern Steppe, near the Caspian Sea in the East.The Surami mountain chain divides Transcaucasia into eastern and western lowlands.While the west of Georgia is characterized by mild winters and hot summers with mean annual air temperatures of 13 to 15 °C and high annual precipitation values (1200-2400 mm), the climate in eastern Georgia is continental with much lower annual precipitation (500-600 mm in the lowlands) and a mean temperature between 10-13 °C.In the mountainous areas mean temperature covers a range of −5 to 10 °C and precipitation varies from 800-1400 mm [22].According to findings by [23] Georgia experienced pronounced summer warming trends for monthly minimum and maximum temperature means and extremes, whereas warm extremes show larger trends than cold extremes.Moreover, the trend for warm spells was observed to be significantly increasing since the 1960s.Heat waves in Georgia are one of the most common natural hazards observed and cause increasing human health impacts and economic losses [24].Most vulnerable groups to severe heat waves can be characterized by poor health, high age and/or pronounced social isolation [25].The rapid increase of the population and urbanization in Georgia and its strong dependence on agricultural production might amplify these negative effects.Tbilisi, Georgia's main capital is by far the most populated city in Georgia, with over one million inhabitants and a rapid rising number of residents.Particularly, urban regions are exposed to more frequent heat wave events, due to the urban heat island effect and synergistic interactions between heat waves and urban heat islands, such as the lack of surface moisture in urban areas, the low wind speed associated with heat waves, the increase in the ambient temperatures and the difference between urban and rural temperatures [26].According to [27][28][29][30], who investigated recent and future heat wave changes in the East Mediterranean and Middle East (EMME) region, the health risk factor due to the projected increase in the number of heat waves, the decrease in precipitation and rising air pollution, will increase, particularly in urban areas.Heat wave impacts are currently under-reported in Georgia and information on morbidity, mortality and economic consequences are difficult to assess.During the last decades Georgia made great efforts in data recovery and monitoring extreme events and their impacts.Currently, 13 meteorological stations and 20 meteorological posts measuring temperature minimum and maximum and precipitation are working on the territory of Georgia [31].In the course of data recovery and digitalization, the National Environmental Agency of Georgia (NEA) kindly contributed metadata to this study giving information about observation periods, and changes of locations of meteorological stations and post.In past studies most heat wave indices are defined by temperature exceeding a fixed or percentile threshold for a given period, consider maximum or minimum and focus on consecutive or single days where defined conditions above the threshold persist [14,19,32,33].A common definition of heat waves in scientific literature does not yet exist.
Following [34,35] the heat wave indicator in the present study is defined by three or more consecutive days above positive Excess Heat Factor conditions.The EHF considers the local geographic acclimatization to temperature, the total heat load, and the recent deviation in temperature from mean temperature and therefore provides a comparative measure of intensity, load, duration and spatial distribution of a heat wave event [25].Moreover, indices based on the 90 th percentile for maximum temperature and the 90 th percentile for minimum temperature have been investigated to compare trends in heat wave characteristics for temperature minimum and maximum series separately.Based on [18,36] multiple heat wave aspects are studied for all three indices: the heat wave number, duration, participating days, the peak and mean magnitudes.The heat wave indices and aspects are recommended by the Commission for Climatology (CCl) Expert Team on Climate Risk and Sector-specific Climate Indices (ET CRSCI) of the World Meteorological Organization (WMO) World Climate Programme.Three new heat wave aspects have been added to the analysis by applying a severity threshold to each time-series: the number of positive, severe and extreme heat wave days, enabling the identification of potential heat-health impacts [37].Heat-health related morbidity and mortality based on the EHF have been investigated in various studies [11,25,[36][37][38][39].The aim of this study is to provide a better understanding of the climatology and changes in the frequency, duration and intensity of summer heat wave events over Georgia and the case study Tbilisi, as the most urbanized city in Georgia.Climatologies and trends of heat wave aspects are presented using a new databank of homogenized daily temperature series from 22 stations well distributed throughout Georgia.Section 2 describes the quality control, homogenization methods of the used observation data, the heat wave indices and aspects investigated and methods used for determining heat wave severity.The annual and seasonal climatology and annual changes in the intensity, duration and frequency of heat waves over Georgia and Tbilisi are presented and discussed in Section 3.Moreover, the five most severe heat waves over Tbilisi during 1961 and 2010 have been identified and illustrated.In Section 4 results are discussed and conclusions are summarized.

Data Quality and Homogeneity Adjustment
Investigating the climatology and trends on observed heat waves over Georgia 87 daily minimum and maximum temperature series over Georgia for the period 1936 to 2010 were provided by the National Environmental Agency of Georgia.The analysis period 1961-2010 was chosen to study changes in heat wave characteristics under anthropogenic influenced climate conditions as well as to optimize spatial coverage and the number of stations available for the trend analysis.During 1988During , 1992During , and 1993 data availability for observation data records was very low and had to be rejected from the analysis of Georgian-averaged climatologies and changes.Climatologies of all heat wave aspects over Georgia have been investigated for the normal period 1961 to 1990.Temperature minimum and maximum series with more than 20% missing values within the analysis period were excluded.Data quality control has been carried out using the computer program RClimDex Software version 1.1 available on http://etccdi.pacificclimate.org.Data quality was tested in order to label potentially wrong values, and to reject them from the analysis.Unphysical values, such as Tmax ≥ 70 °C, Tmin ≤ −50 °C, Tmax ≤ Tmin were identified and set to missing values.Outliers were detected and rejected for daily maximum and minimum temperature exceed ± four standard deviation.After quality control, a final number of 31 minimum and maximum temperature series were left for homogeneity testing and adjustment.During the index calculation process, the following data quality requirements have been applied in order to include as many Georgian temperature series as possible: (1) A summer value is calculated if all months are present (May to September); (2) A month is considered as complete if ≤ 3 days are missing; (3) A station will be rejected from the analysis if more than 5 consecutive months are missing.For threshold indices, a threshold is calculated if at least 70% of data are present.For the analysis period 1961-2010, 31 temperature minimum and maximum series matched all quality criteria.Inhomogeneity of time-series was tested using the software package RHtestV3 and to adjust significant breakpoints.Metadata provided by the National Environmental Agency include information of regarding the station's name, coordinates, altitude, WMO code, observational periods, missing data during an observation period, and the station's relocation date.In order to provide reliable results on the climatology and changes in heat waves, time-series with significant breakpoints not documented in the metadata were excluded from the study.Thirty-one stations were left matching all quality criteria as described above.Seven significant breakpoints were found, whereas not more than one breakpoint per maximum or minimum time-series was detected.For two stations (Tbilisi and Gori) dates of site moves were noted within the metadata, which could be confirmed by the detected breakpoints.In this study the Quantile-Matching (QM) adjustment procedure of RHtestV3 is used to adjust daily temperature series [40,41].Details on the homogeneity testing and adjustment procedure and parameter usage during the QM adjustment procedure are stated in [23].The resulting homogenized daily temperature dataset comprise minimum and maximum temperature series from 22 stations for the analysis period of 1961 to 2010 (see Appendix).

Heat Wave Indices and Characteristics
In 2011 a "core set" of 34 descriptive sector-specific indices has been defined by the World Meteorological Organization (WMO) Commission for Climatology (CCl) Expert Team on Climate Risk and Sector-specific Indices (ET CRSCI) to improve decision-making for planning, operations, risk management and for adaptation to both climate change and variability.These internationally agreed indices were developed in part from the core set of indices that are developed and maintained by the Expert Team on Climate Change Detection and Indices (ETCCDI) monitoring changes in "moderate" extremes.However, the new core set of indices provide nine additional indices, such as five heat wave indices, to facilitate the use of climate information in users' decision-support systems for climate risk management and adaptation strategies [36].All selected indices were calculated using the ClimPACT software, an R-based software provided to calculate the indices from the ET CRSCI website [42].The present study uses the Excess Heat Factor (EHF) as heat wave index defined by [25].Both maximum and minimum temperatures are used in this assessment.The EHF combines a measure of excess heat, the deviation from long term mean temperature, calculated with respect to the period 1961-1990, to take into account local geographic acclimatization, and heat stress, the deviation in temperature from mean temperature for the previous 30 days to measure short-term acclimatization.EHF values are calculated from a three-day mean of forecast temperatures to derive an index of heat wave intensity.Two sub-indices are combined to produce the complete EHF index.The first is a measure of significant excess heat relative to local climatic conditions, the 95th percentile of mean temperature conditions (EHFsig): The second sub-index is a measure of short-term acclimatization to heat, relative to the mean temperature of the previous 30 days (EHFaccl): These two indices are combined to generate the EHF index.The unit of EHF is °C2 : EHF incorporates the effect of humidity on heat tolerance indirectly, by using the mean, rather than the maximum daily temperature, in the calculation.It provides a comparative measure of intensity, load, duration and spatial distribution of a heat wave event and has a strong signal-to-noise ratio.With an increase in Excess Heat (EHFsig) and Heat Stress (EHFaccl), EHF increases as a quadratic response to increasing heat load.Heat waves are defined as a period of at least three days with EHF ≥ 0 and the combined effect of excess heat and heat stress with respect to the local climate [25].The heat wave index EHF was selected with regards to eight heat wave aspects quantifying the intensity, frequency and duration of a heat wave event.This multi-aspect framework is based on that of [18], which was slightly extended by [11] to the following five attributes of heat wave aspects: 1. HWN-the yearly number of heat waves; 2. HWD-the length (in days) of the longest yearly event; 3. HWF-the sum of participating heat wave days per year 4. HWA-the hottest day (amplitude) of the hottest yearly event; 5. HWM-the mean event intensity averaging all participating event days Further to these five heat wave aspects three new aspects have been included in the heat wave analysis based on the Excess Heat Factor: HWday-the number of heat wave days, HWsev-the number of severe heat wave days and HWex-extreme heat wave days.These heat wave aspects enable to differentiate between heat wave days, severe and extreme heat wave days and include the aspect of heat-health considering the exceedance of a station's severe EHF threshold (Section 2.3).The resulting heat wave aspects analyzed are summarized in Table 1.All eight aspects are calculated annually over the summer season, which is defined as a period from May to September (153 days).In order to compare the differences between Georgia-averaged trends in heat wave characteristics based on temperature minimum and maximum series separately the heat wave aspects HWN (the yearly number of heat waves), HWD (the length of the longest yearly event), HWF (the sum of participating heat wave days per year), HWA (the hottest day the hottest yearly event), and HWM (the mean event intensity averaging all participating event days) based on the 90th percentile for maximum temperature (CTX90pct), the 90th percentile for minimum temperature were calculated (CTN90pct).Heat wave indices based on daily Tmin and Tmax are defined as follows: (1) CTN90pct-the calendar day 90th percentile of daily Tmin calculated for a five-day window centered on each calendar day in the base period, (2) CTX90pct-the calendar day 90th percentile of daily Tmax, as described for Tmin.The unit of CTN90pct and CTX90pct is °C.Due to the low amount of measurable events a set of 14 time series for all three heat wave indices was used for this analysis.All trends were calculated by the non-parametric Sen's slope estimator based on Kendall's tau (τ) [43].The annual slopes of trends were converted into slope per decade.The statistical significance has been estimated using the Mann-Kendall test, whereas in the present study a trend was considered to be statistically significant if it was less than or equal to a level of 5% [44,45].Station trends have only been estimated if data requirements are met as described above.Apart from trends for each individual station, trends were also averaged for all Georgian station records.These trends were calculated as the arithmetic average of the summer index values at all stations.

Severe and Extreme Heat Waves
In terms of heat-health, low values of EHF may be considered as uncomfortable and have no or only little impact.Severe heat wave impacts on heat-health are relatively rare, due to the local adaption capacity of the affected population based on the long-term experience to cope with heat waves of low to moderate intensity.In the present study severe heat waves are defined by an event where EHF values exceed a threshold for severity that is specific to the climatology of each location following [25].To identify severe and extreme heat wave events during 1961 and 2010 over Georgia a severe EHF threshold for each station has been detected.Using the example of Tbilisi time series (Figure 1) it is shown that the upper end of the EHF distribution is fat-tailed which allows the identification of a threshold between more-frequent lower-intensity events and less-frequent higher intensity events.The severe EHF threshold is calculated empirically as the 85 th percentile of the distribution of positive EHF values (EHF85) based on the observation record at a given location.This method ensures that all EHF values are truly representative of each site's climatology and avoids potential errors by modeling a distribution.Extreme heat waves are defined as an event where EHF values are well in excess of the severity threshold and result in a wide impact based on a cascade of failing systems [25].In the present study extreme heat waves occur if an EHF value during a heat wave at least doubles the severity threshold (EHF ≥ 2 × EHF85).

Heat Wave Climatology
In this section the climatological Georgia-average during the period 1961-1990 is investigated with respect to selected heat wave aspects.In Table 2, it is shown that the Georgia-average HWN amounts to 1.7 events/year.Highest numbers of annual events range between 1.7 and 1.9 events/year and can be found throughout the whole territory of Georgia.Lowest numbers of HWN are observed at the Western coast (1.3-1.5 events/year), while a general distribution for all heat wave numbers across Georgia does not show clear spatial patterns (Figure 2a).However, spatial patterns can be found for HWD, which measures 5,5 days/year in the Georgia-average.Highest averaged values for HWD are concentrated in the southeastern dry steppe (6-7 days/year).Heat waves with a mean duration of 5 to 6 days per year are predominantly located in the Western part of Georgia.The shortest heat waves (4-5 days/year) are observed in central Georgia (Figure 2b).The Georgia-average for HWF (yearly sum of participating heat wave days) measures 10.4 days/year (Table 2).Eleven to 15 days contribute at most to all summer heat waves in Georgia.Similar to HWD, the highest HWF values can be found in the southeast, lowest values in central Georgia (9-10 days/year) (Figure 2c).HWA means measure 12.2 °C2 /year between 1961 and 1990 (Table 2).Highest station heat wave amplitudes per year measure 17-24 °C2 averaged over Georgia.Strongest heat wave peaks are located in the western part of Georgia with highest HWA means (17-25 °C2 /year) in the southwest (Figure 2d).Similar to HWA, highest station means for HWM (6-9 °C2 /year) are concentrated in the southwestern part of Georgia (Figure 2e) with a Georgia-average of 4.3 °C2 /year (Table 2).Lowest magnitudes (2-4 °C2 /year) are predominantly found in the eastern part of Georgia.Most positive EHF values (11-12 days/year) can be found in the Southeastern part of the study area, while most severe (40-50 days/year) and extreme days (0.4-0.5 days/year) can be observed in the southeast and the southwest (Figure 2f-h).Georgia-averaged climatologies for HWday measure 10.3 days/year.For severe and extreme days an annual mean of 1.3 days/year and 0.3 days/year can be found, respectively (Table 2).The monthly distribution of the mean number and mean intensity (positive EHF values) of summer HWday (Heat Wave Days), HWsev (Severe Heat Wave Days) and HWex (Extreme Heat Wave Days) between 1961 and 2010 is shown in Figure 3. Highest numbers of HWday can be found in August (3.4days/month), lowest in May (1.8 days/month).However, in the Georgia-average most HWsev and HWex occur during May with decreasing numbers toward September.The highest mean intensity for

Changes in Heat Wave Characteristics
In this section heat wave changes in the Georgia-average during the period 1961 to 2010 are investigated.
As Table 3 shows, all trends were found significantly positive, with the exception of the annual trend for HWM (average magnitude of all yearly heat waves), which shows an insignificant trend of approximately zero.For HWN a significant increasing trend of 0.4 events/decade can be observed (Table 3).All single station trends throughout Georgia were found to be positive (Figure 4a).During the period 1961-2010, 12 of 22 stations analyzed indicate significantly positive trends.Five of them show trend magnitudes of 0.5-1 events/decade, while 10 of all stations are approximately zero, indicating that the yearly number of heat waves for these stations did not change during 1961 and 2010.For HWD (the length of the longest yearly heat wave event) a significant magnitude of 0.9 days/decade was found in the Georgia-average (Table 3).Similar to HWN, no decreasing trend could be observed throughout Georgia.Twelve stations show significant trends, of which four stations show trend magnitudes between 1 and 2 days/decade (Figure 4b).The heat wave aspect HWF (yearly sum of participating heat wave days) shows the largest trends of up to 6 days/decade during the period 1961 to 2010 (Figure 4c).Similar to HWN and HWD all station trends are increasing and stations with statistically significant trends for HWN and HWD are also statistically significant for HWF.The Georgia average trend amounts to 2.9 days/decade (Table 3).These findings are in accordance with [36], demonstrating that the high trend magnitude of HWF drives increases in HWN and HWD, as the number of heat wave days represents an influencing factor in the calculation of event length and occurrence.HWA (the averaged peak of the hottest summer day) and HWM (the average magnitude of all summer heat waves) show lower and less significant station trends.Although few negative station trends can be observed for HWA a significant positive trend of 1.1 °C/decade can be found in the Georgia average (Table 3).HWA and HWM (the average magnitude of all summer heat waves) show lower and less significant averaged station trends.For HWM an insignificant averaged trend, approximating zero, can be observed, indicating no change of the mean magnitude of all summer heat waves during 1961 and 2010.However, regarding Georgia-averaged HWday (positive EHF days), HWsev (severe days) and HWex (extreme days) significant increasing trends could be found.For HWday all station trends are increasing and a high proportion of 13 stations show significant trends measuring up to 6 days/decade.Trends with the highest magnitudes are predominantly located in Eastern and Western Georgia (Figure 4f).In the Georgia-average a significant positive trend of 3.3 days/decade for HWday was observed.Severe days and extreme days increased at a significant rate of 0.4 days/decade and 0.05 days/decade, respectively (Table 3).Additionally, for HWsev all station trends are increasing, while only three stations show significant trends.Due to the low number of observed extreme days per year, trends for HWex show either insignificant slopes throughout Georgia or trends approximating zero.There are predominantly positive trends noted, only two station trends were found to be insignificantly negative (Figure 4h).In order to be able to compare trends between temperature minimum and maximum series separately and to investigate which of the series shows more response to trends in heat wave aspects, trends for the 90 th percentile maximum temperature (CTX90pct), the 90th percentile minimum temperature were calculated (CTN90pct).It is shown, that all trends for CTN90pct demonstrate a higher magnitude than for CTX90pct, which corresponds to global heat wave observations [11].
As shown in Table 3, heat wave trends for HWN measure a significant increase of 0.50 heat wave events/decade averaged over Georgia, regarding the heat wave index CTN90pct.However, for CTX90pct a significant trend at a rate of 0.41 events/decade can be observed, suggesting that consecutive days of extreme high temperature maxima are less common over Georgia than consecutive nights of extreme temperature minima.Regarding the heat wave aspect HWD, a significant increase for both indices could be found, while CTN90p shows a higher trend magnitude averaged over Georgia (0.65 days/decade).Similar to results for the heat wave index EHF Georgia-averaged trends for HWF are of higher magnitude compared to HMN and HWD.CTN90pct measures a significant positive trend of 2.8 days/decade.For CTX90pct a lower trend of 1.8 days/decade can be observed.Additionally, for HWA and HWM the index CTN90pct shows a higher trend magnitude throughout Georgia, than CTX90pct.CTN90pct-HWA and CTN90pct-HWM measure significant increasing trends of 0.53 °C/decade and 0.32 °C/decade, respectively.For CTX90pct-HWA and CTX90pct-HWM insignificant positive trends were observed.
Figure 5 shows, that although the quantitative observations of the annual values and trend magnitudes vary for all three indices (Table 4), the annual time-series show qualitatively similar patterns among each of the three indices.A quantitative comparison cannot be conducted, due to the different number of stations used for the analysis on heat wave indices EHF and CTN90pct/CTX90pct.Table 4. Georgia-averaged trends in CTN90pct and CTX90pct between 1961 and 2010 with respective confidence intervals (95%).Trends significant at the 5% level are indicated in bold and highlighted in green.

Severe Heat Waves over Tbilisi
In this section the climatology and changes in heat waves over Tbilisi station have been investigated, considering the city as an urban heat island.Moreover, the five most severe heat waves between 1961 and 2010 have been identified and illustrated.
Tbilisi experiences by far the highest trend magnitudes for all summer temperature variables analyzed compared to the Georgian mean, as well as the surrounding stations (not shown).As shown in Table 5, summertime minimum and maximum temperature averages 16 °C and 27.4 °C, respectively.The annual summer trends measure 0.32 °C/per decade for Tmin and 0.38°/decade for Tmax.Comparing these trends with the surrounding neighboring stations (not shown) or the Georgian-averages for summertime temperature minimum (0.20 °C/decade) and maximum (0.31 °C/decade), Tbilisi shows pronounced high trend magnitudes for Tmin and Tmax and a lower temperature variability.Moreover, for warm nights a much higher trend magnitude (3.4 °C/decade) is observed during summer than for warm days.In addition, compared to its neighboring stations Tbilisi shows a higher trend magnitude (not shown).Regarding Tbilisi's heat wave climatology the city shows a HWN value of 0.7 events/year, the mean heat wave duration measures 6.2 days/year and 11.0 days participate in yearly heat wave events.The peak and mean magnitude for Tbilisi averages 8.2 °C2 /year and 3.2 °C2 /year, respectively.For the number of heat wave days a mean value of 10.9 days/year could be found, while the climatologies for severe heat wave days (extreme heat wave days) measure 1.2 days/year (0.1 days/year).While Tbilisi shows lower means for HWA and HWM, it demonstrates similar mean climatological values for HWN, HWday, HWsev, and highly exceeds the Georgia-averages for the heat wave aspects HWD and HWF (Tables 2 and 3). Figure 6 shows the annual numbers of heat wave days, severe and extreme heat wave days.It is notable that the largest proportion of low-and high-intensity heat wave days can be found in the last 20 decades.An increasing number of heat waves of 0.7 events/decade were observed, the heat wave length increased by 1.0 days/decade and the trend for HWF measured 2.9 days/decade.Additionally, for HWA (0.9 °C2 /decade) and HWday (5.5 days/decade) Tbilisi station shows a higher trend magnitude, than the Georgia averages and its neighboring stations (not shown).However, for the number of severe heat wave days Tbilisi shows a slightly lower trend magnitude (0.32 days/decade) compared to the Georgia-average.For HWM there is a trend measured approximating zero.Due to the low number of extreme heat wave days distributed over 50 summers the trend magnitude for HWex was not denoted.Table 5. Climatology and annual trends for summer Tmin and Tmax, temperature extremes and selected heat wave aspects at Tbilisi, 1961 and 2010.Trends and respective confidence intervals (95%) significant at the 5% level are indicated in bold and highlighted in green.Table 6 presents the five most severe heat wave events measured at Tbilisi station between 1961 and 2010 and their intensity (peak EHF value, °C2 ), EHF event duration (days) and accumulated heat load (°C 2 ).The accumulated heat load is defined as the sum of the consecutive positive EHF values.By means of comparisons of international extreme heat wave events and their impacts, the peak heat load (highest EHF value) over a station's severe EHF threshold represents the most crucial indicator determining the severity of a single station heat wave [1,25].
It is notable that during 1961 and 2010 four of five top-ranked severe heat wave events were found in the second half of the analysis period.The highest peak EHF value observed was measured between 25 May and 3 June 2007.The EHF value of 36.9 °C2 is over six times the magnitude of Tbilisi's severe EHF threshold (7.1 °C2 ), which classifies this heat wave as an extreme event.The high average heat load and EHF event duration of eight severe days result in a huge accumulated heat load (203 °C2 ) and impose high amplitudes of heat stress on the population, and other affected sectors (Figure 7a).The 1971 EHF time series shows Tbilisi's second ranked extreme heat wave.The peak magnitude for this event is almost four times the severe EHF threshold (26.2 °C). Figure 7b shows a short period of positive EHF values with a relatively high average heat load of 12.9 °C2 during the heat wave.The third ranked peak heat load observed (23.7 °C2 ) was measured in 24 May 1995 and is over three times the magnitude of Tbilisi's severe EHF threshold.It was found among six severe heat wave days (Figure 7c).The fourth highest classified heat wave event observed at Tbilisi station occurred in August 1998.The event peak of 20.1 °C is six to seven times the magnitude of Tbilisi's severe EHF threshold (Table 6) and is one of 14 heat wave days and nine consecutive severe EHF days (not shown).The fifth ranked peak EHF event of 17.9 °C2 exceeds the station's severe EHF threshold by circa 2.5 times.

Discussion and Conclusions
This study analyzed climatology and changes in summer heat wave events over Georgia since the 1960s using a new dataset of daily minimum and maximum temperature series.The data was carefully quality controlled and homogeneity was tested using the software RClimDex 1.1.Metadata could be used to identify the date of breakpoints in a time-series.Inhomogeneous time-series could be adjusted using RHVtest 4. For the investigation of heat wave climatology, changes and single severe single heat waves over Georgia and at Tbilisi station the Excess Heat Factor was, chosen a combined indicator of excess heat and heat stress and comparative a measure of intensity, load, duration and spatial distribution of a heat wave event.The EHF was analyzed with respect to eight heat wave aspects: the yearly number of heat waves (HWN), length of the longest yearly event (HWD), the yearly sum of participating heat wave days (HWF), the hottest day of hottest yearly event (HWA), and the average magnitude of all yearly heat waves (HWM).Three new heat wave aspects have been applied to the study: the number of heat wave days (HWday), the number of severe heat wave days (HWsev) and extreme heat wave days (HWex) to allow a differentiated assessment of the spatial and temporal characteristics of lowintensity, severe and extreme heat waves.Severe and extreme heat waves have been detected applying a threshold, the 85th percentile of the distribution of positive EHF values, to each station.This multiaspect heat wave analysis provides a more profound study of heat wave changes in Georgia than carried out in previous scientific studies.Next to EHF, the heat wave indices CTN90p and CTX90p were used to study heat wave aspects regarding temperature minimum and maximum series separately.The study could improve the understanding in the spatial and temporal characteristics and changes in past summer heat wave events over Georgia.Climatologies across the study area show clear patterns for most of the heat wave aspects.Stations with highest climatological values for HWD, HWF and HWday are located in the Southeast.Highest values for mean HWA and HWM are found in the Southwest of the study area.For HWsev and HWex both regions demonstrate highest climatological values.For all aspects lowest heat wave impacts can be expected in central Georgia.The seasonal distribution of HWday, HWsev and HWex between 1961 and 2010 implies the occurrence of highest heat wave intensities in May and September.Highest numbers of high-intensity heat wave days are measured in May, while most low-intensity heat wave days can be found during August.Regarding changes in heat wave events almost all heat wave aspects analyzed based on the indices EHF, CTX90pct and CTN90pct show significant increasing trends in the Georgia-average.However, HWA and HWM show predominantly insignificant and only small changes during the analysis period.HWday and HWF show by far the largest trend magnitudes across all heat wave indices studied (EHF, CTN90pct and CTX90pct), which is consistent with findings of [11,18], stating that occurrence-based heat wave aspects possess larger trend magnitudes.They also found that the high magnitudes of trends for HWF are driving increases in HWN and HWD, due to the fact that the overall number of heat waves (and their duration) will increase when the number of participating days increases.Spatial patterns of trends over Georgia are difficult to assess, which can be attributed to the moderate data quality of several time-series.Nevertheless, it is notable that for the heat wave aspects HWN, HWF, HWD, and HWday the highest proportion of increasing station trends can be observed in the Southeastern and Southwestern region of Georgia.Tbilisi station experienced a pronounced increase in the intensity, frequency and duration of heat waves during 1961 and 2010 compared to the Georgian-average and its neighboring stations.This result corresponds well with findings by [27][28][29][30] projecting an increasing number of heat wave events in the EMME region, particularly in urban areas.Including the strong increase in extreme warm night-time temperatures, these results could be explained by the UHI effect, the enhanced heat absorption and storage in sealed surface during the day and the slow release at night, and combined interactions between heat waves and urban heat islands [26,46].The most severe summer heat waves in Tbilisi have been identified in the years 1971, 1995, 1997, 1998 and 2007.Comparing the indices CTN90pct and CTX90pct, it is shown, that all trends for Tmin (night-time) based events over Georgia increase faster than for days-time events, which corresponds to global heat wave observations [11] and changes for extreme Tmin and Tmax events [16,32,[47][48][49].These findings highlight the importance of selecting the appropriate index related to the most affected sector investigated.Heat wave studies based on minimum temperature rather consider impacts on sectors, such as agriculture, whereas temperature maximum-based heat wave indices are appropriately used to determine impacts on infrastructure.However, for investigations on power consumption and human health impacts a heat wave index based on a combination of minimum and maximum temperature might be the most suitable choice.In order to obtain a more detailed insight in the connections of large-scale drivers on the occurrence of heat waves, a further study on contributions of land-atmosphere coupling and large-scale atmospheric circulation over the Caucasus region is necessary.

Figure 1 .
Figure 1.Determining EHF85 and severe heat wave days for Tbilisi station (1961-2010): Tbilisi cumulative distribution of positive EHF (red line) and the severe EHF threshold at the 85th percentile (black line).

Figure 3 .
Figure 3. Seasonal distribution of (a) the mean number (days) and (b) the mean EHF value (°C) of HWday (green bars), HWsev (yellow bars) and HWex (red bars) averaged over Georgia for the period 1961-2010.

Figure 4 .
Figure 4. Station trends of selected heat wave aspects over Georgia for the period 1961-1990: (a) Heat Wave Number (no. of events); (b) Heat Wave Duration (days); (c) Heat Wave Frequency (days); (d) Heat Wave Amplitude (°C 2 ); (e) Heat Wave Magnitude (°C 2 ); (f) Heat Wave Days (days) and (g) Severe Heat Wave Days (days); and (h) Extreme Heat Wave Days (days).Red triangles indicate warming trends, blue indicate cooling trends.Light blue and red triangles indicate trends not significant at the 5% level.

Figure 5 .
Figure 5.Comparison of annual (a) Heat Wave Duration and (b) Heat Wave Frequency for heat wave indices EHF, CTN90pct and CTX90pct during the period 1961-2010.

Figure 6 .
Figure 6.The annual number of heat wave days (green bars), severe heat wave days (yellow bars) and extreme heat wave days (red bars) at Tbilisi station, 1961-2010.

Figure 7 .
Figure 7.The three most severe summer heat wave events over Tbilisi during 1961 and 2010: Observed daily EHF values and severity threshold (dashed red line) during May 2007 (a); May 1971 (b) and May 1995 (c).Red dots indicate EHF values above the station's severity threshold (EHF = 7.06 °C2 ).The three-day-averaged daily mean temperature (yellow line) is superimposed.

Table 1 .
Heat wave aspects analyzed based on the Excess Heat Factor (EHF) with index names, definitions, and units.
3C2HWM Heat wave meanAverage magnitude of all summer (May-Sep) heat wave days where conditions persist for at least3consecutive days with positive EHF values °C2 HWday Heat wave days The annual number of all summer (May-Sep) heat wave days with positive EHF values days HWsev Heat wave severe days The annual number of all severe summer (May-Sep) heat wave days with positive EHF values above the station's severe EHF threshold at EHF 85 (85th percentile of the station's distribution of positive EHF values) days HWex Heat wave extreme days The annual number of all extreme summer (May-Sep) heat wave days with positive EHF values that at least double the station's severe EHF threshold at EHF 85 (85th percentile of the station's distribution of positive EHF values) days

Table 3 .
Annual averaged trends for investigated heat wave aspects between 1961 and 2010 with respective confidence intervals (95%).Trends significant at the 5% level are indicated in bold and highlighted in green.

Table 6 .
Top ranked severe heat waves measured at Tbilisi station during 1961 and 2010 listed by year: Rank, peak day of event (Date), peak EHF value (°C 2 ), accumulated and average heat load (°C 2 ) and EHF event duration (days).Tbilisi's severe EHF threshold accounts for 7.06 (°C 2 ).