Heatwaves in the Future Warmer Climate of South Africa

: Weather and climate extremes, such as heat waves (HWs), have become more frequent due to climate change, resulting in negative environmental and socioeconomic impacts in many regions of the world. The high vulnerability of South African society to the impacts of warm extreme temperatures makes the study of the e ﬀ ect of climate change on future HWs necessary across the country. We investigated the projected e ﬀ ect of climate change on future of South Africa with a focus on HWs using an ensemble of regional climate model downscalings obtained from the Conformal Cubic Atmospheric Model (CCAM) for the periods 2010–2039, 2040–2069, and 2070–2099, with 1983–2012 as the historical baseline. Simulations were performed under the Representative Concentration Pathway (RCP) 4.5 (moderate greenhouse gas (GHG) concentration) and 8.5 (high GHG concentration) greenhouse gas emission scenarios. We found that the 30-year period average maximum temperatures may rise by up to 6 ◦ C across much of the interior of South Africa by 2070–2099 with respect to 1983–2012, under a high GHG concentration. Simulated HW thresholds for all ensemble members were similar and spatially consistent with observed HW thresholds. Under a high GHG concentration, short lasting HWs (average of 3–4 days) along the coastal areas are expected to increase in frequency in the future climate, however the coasts will continue to experience HWs of relatively shorter duration compared to the interior regions. HWs lasting for shorter duration are expected to be more frequent when compared to HWs of longer durations (over two weeks). The north-western part of South Africa is expected to have the most drastic increase in HWs occurrences across the country. Whilst the central interior is not projected to experience pronounced increases in HW frequency, HWs across this region are expected to last longer under future climate change. Consistent patterns of change are projected for HWs under moderate GHG concentrations, but the changes are smaller in amplitude. Increases in HW frequency and duration across South Africa may have signiﬁcant impacts on human health, economic activities, and livelihoods in vulnerable communities.


Introduction
Characteristics of extreme temperature events are expected to change due to global warming/global mean temperature increases. A recent study indicated that by 2017 the global mean temperature had increased by 1 • C with respect to preindustrial levels [1]. A more recent report [2] indicated that a global mean increase of about 1.5 • C could be seen in the future. Previous studies (e.g., [3,4]) have shown that there is a positive trend in both daily minimum and maximum temperatures across South Africa.

Data and Methods
The CCAM model was used to simulate present-day and future global climate for the period 1961-2099 with surface forcing from six different Global Climate Models (GCMs) forced by two different emission scenarios, namely RCP 4.5 and RCP 8.5. The CCAM was developed by the Commonwealth Scientific and Industrial Research Organization (CSIRO) in Australia [30,31] and is based on a conformal cubic grid with a semi-Lagrangian and semi-implicit dynamical core that solves hydrostatic primitive equations [32]. CCAM is a global model, but it can also be employed as a regional climate model (RCM) using stretched-grid variable-resolution mode [33,34]. CCAM includes a prognostic cloud scheme [35], a cumulus convection scheme with a mass-flux closure [36], the Geophysical Fluid Dynamics Laboratory (GFDL) parameterizations for long-wave and short-wave radiation, and a stability-dependent boundary layer scheme based on Monin Obukhov similarity theory with non-local treatment [37,38].
The CCAM model has here been integrated coupled to the CSIRO Atmosphere-Biosphere Land Exchange model (CABLE). The simulations used in this study were generated on a quasi-uniform global grid of about 50 km resolution in the horizontal and with the RCP 4.5 and RCP 8.5 to account for the uncertainty in the emission scenarios. The RCPs were developed by the Intergovernmental Panel on Climate Change (IPCC) in its Assessment Report Five (AR5), superseding the Special Report on Emissions Scenarios (SRES) [39]. RCP 4.5 represents a moderate GHG concentration future, whilst RCP 8.5 represents a high GHG concentration future. The CCAM model was forced with sea ice and sea surface temperatures from six earth system models (Table 1) under RCP 4.5 and RCP 8.5 forcing conditions and results in a total of 12 ensemble members were produced.
There are several definitions of HWs in different regions across the world. This work adopted the South African Weather Service (SAWS) definition that a HW is "when for at least three consecutive days the maximum temperature of a certain region or grid point is five degrees or higher than the average mean maximum for the hottest month for that particular station or grid point". This definition is similar to the one used in previous studies on HWs (e.g., [40]), while other studies prefer to define HWs as events of at least six consecutive days (e.g., [33]). There are several different thresholds used when defining a HW, based on weather conditions of the area being considered. The threshold used to identify HWs for all the present-day and future periods are based on the threshold that define HWs under present-day climate. Therefore, the definition used does not take into account possible changes in damage and impact due to adaptation. This definition also does not capture relatively extreme events occurring during the colder seasons that may not reach this extreme temperature threshold, but that could result in strong ecological impacts, such as affecting the growing season, and agricultural and ecological development. The definition of HWs used in this study will ensure that shorter-lived events are also detected. The use of a threshold linked to maximum temperature average for each location (grid point in the case of a model simulation) makes it possible for the HWs to vary based on the climatology of each region in South Africa.
Our analyses were split into 30-year averages, where the period 1983-2012 was considered as the present-day baseline, and the periods 2010-2039, 2040-2069, and 2070-2099 were considered as representing the near, mid, and far-future respectively. Simulations obtained from emission scenarios RCP 4.5 and RCP 8.5 are shown separately.
Three different measures of the characteristics of HWs were considered, namely, the frequency, duration, and intensity. Frequency was calculated as total actual occurrences per seasons per 30-year period while duration is the total number of days a HW lasts. The HW intensity was calculated using the average maximum temperature during all the days identified as forming part of the HW event.
The present-day simulations were compared with SAWS observations and also to HWs in future periods. The simulated HWs were identified using thresholds identified for each grid point after interpolation of the simulations to a 0.5 • resolution latitude-longitude grid. Those obtained from the SAWS observations were identified for each station. South Africa has a good network of surface temperature stations compared to other countries in southern Africa, however for this study 24 stations were selected countrywide (Figure 1) based on a 95% data availability criteria during the study period.
Atmosphere 2020, 11, x FOR PEER REVIEW 4 of 19 interpolation of the simulations to a 0.5° resolution latitude-longitude grid. Those obtained from the SAWS observations were identified for each station. South Africa has a good network of surface temperature stations compared to other countries in southern Africa, however for this study 24 stations were selected countrywide (Figure 1) based on a 95% data availability criteria during the study period. RClimDex version 1.0 software [41] was used to quality control the weather station data. RClimDex is an R based statistical tool which is freely downloadable from the Expert Team on Climate Change Detenction, Monitoring and Indices (ETCCDMI) website [42] and outputs of this are shown in Tables 2 and 3. The next section will discuss the results starting with a comparison of simulated present-day events and observations. The projections for the three future periods will subsequently be discussed, followed by a discussion and concluding remarks.  RClimDex version 1.0 software [41] was used to quality control the weather station data. RClimDex is an R based statistical tool which is freely downloadable from the Expert Team on Climate Change Detenction, Monitoring and Indices (ETCCDMI) website [42] and outputs of this are shown in Tables 2  and 3. The next section will discuss the results starting with a comparison of simulated present-day events and observations. The projections for the three future periods will subsequently be discussed, followed by a discussion and concluding remarks. CNRM-CM5 (CNRM-CERFACS, France) [45] GFDL-CM3 (NOAA, USA) [46] MPI-ESM-LR (MPI, Germany) [47] NorESM1M-M (NCC, Norway) [48] Atmosphere 2020, 11, 712 5 of 18  Figure 2a shows the average observed maximum temperature across South Africa during mid-summer, December-February (DJF), for the period 1983 to 2012. Figure 2c shows the 50th percentile of the simulated average temperature for the RCP 4.5 scenario during the same period. The observed and simulated 30-year daily average maximum temperature indicate that South Africa experiences highest maximum temperatures during DJF across the northern parts of the country.

Osberved vs. Simulated Average Daily Maximum Temperature (T X ) and HW Thresholds
The highest values are shown in the dry regions in the northwestern parts, with averages of over 36 • C. The model positions the highest temperature in the northwestern parts which is consistent with the observations. The model is able to capture the observed east-west and south-north temperature gradient, as well as the colder temperatures at higher altitudes across the Drakensberg Mountains in the eastern parts of the country. There are small differences amongst different ensemble members under both RCP 8.5 and RCP 4.5 in the present-day climate. Moreover, the RCP 8.5 and 4.5 emissions are similar until the year 2004, and only then start to diverge, however by the end of 2012 the differences between the two scenarios remain small. Figure 2b,d show the observed and simulated average maximum temperatures for mid-winter, June-August (JJA). Lower daily average maximum temperatures are observed in this season, and the model is also able to capture this seasonal variability. The observed and simulated temperatures average to less than 22 • C in much of the interior and coastal regions for winter during the 1983-2012 period.    The Tn10P and Tx10P columns in Table 3 show that there is a decreasing trend in the number of cold days. The two columns indicate that there is a decrease in the number of days associated with lower minimum and maximum temperatures, respectively. The negative trend is significant across 21 of the 24 selected stations for minimum temperatures, and 14 stations for maximum temperatures. The number of days associated with the highest minimum and maximum temperature were increasing at all the stations. For the maximum temperature, only Cape St Francis has a negative trend, which is statistically significant. These results are in agreement with studies that have shown a decrease in cold extreme events, and an increase in hot extreme events (e.g., [8]).
Daily average maximum temperatures are projected to increase by up to 2 • C during the period 2010 to 2039 as compared to the reference period of 1983 to 2012 during summer across the country's plateau (Figure 3a,b). Most central parts of the subcontinent are expected to experience increases of over 1 • C. Regions with projected temperature increases of more than 1 • C extend into the neighboring countries across the subtropical western parts of Africa for all RCP 8.5 projections. The winter temperatures are also projected to increase, with the median of the RCP 4.5 downscalings projecting increases of less than 1 • C across larger part of the central interior. The coastal areas are projected to have a lower increases as compared to inland.
Atmosphere 2020, 11, x FOR PEER REVIEW 7 of 19  is landlocked in South Africa and characterized by high altitudes. Projected temperature increases over the coastal regions are slower because of the moderating effects of the oceans.
The climate change signal is expected to be largest for the period 2070 to 2099 (Figure 3e,f), because it is the furthest and the anthropogenic greenhouse gas concentrations will be largest especially in the RCP 8.5 emission scenario. The summer temperatures are projected to rise further, with increases going in excess of 6 • C in the RCP 8.5 scenario for the period 2070-2099. The projected increase across the coastal regions is also smaller as with the, near future and medium-term projections. All the ensemble members, both in the RCP 4.5 and RCP 8.5, show an increase in temperature across the whole of South Africa, which indicates that there is a high confidence in the temperature projections. The increase in the mean temperatures are expected to have an effect on the characteristic of HWs in the future climate.

Observed vs. Simulated HWs
The SAWS definition of a HW was used, where all events with an observed or simulated temperature of 5 • C or greater above the daily average maximum temperature of the hottest month, for a minimum of three consecutive days, are identified, for each station or grid point, respectively. The threshold calculation is made based on present-day climate simulations, and the present-day thresholds are applied to future projections. These thresholds are used by the Disaster Risk Reduction (DRR) division of SAWS which issues HW warnings to the public. The HW thresholds are generally higher across the western parts of the country (which are also arid), ranging from 32 to 44 • C, in line with observed average temperatures shown in Figure 2a during DJF (when most HW are experienced). The thresholds are also relatively high in the north-east lowveld, ranging from 32 to 40 • C.
South Africa experiences a considerable spatial variability of HW occurrences seasonally (Figure 4). The country experienced higher HW frequency in the northern parts, with Punda Maria (22.68 • S, 31.02 • E) recording 45 HWs from 1983 to 2012. The north-eastern low-lying regions generally experience relatively high HW occurrences being warm throughout the year.
Atmosphere 2020, 11, x FOR PEER REVIEW 8 of 19 especially in the RCP 8.5 emission scenario. The summer temperatures are projected to rise further, with increases going in excess of 6 °C in the RCP 8.5 scenario for the period 2070-2099. The projected increase across the coastal regions is also smaller as with the, near future and medium-term projections. All the ensemble members, both in the RCP 4.5 and RCP 8.5, show an increase in temperature across the whole of South Africa, which indicates that there is a high confidence in the temperature projections. The increase in the mean temperatures are expected to have an effect on the characteristic of HWs in the future climate.

Observed vs. Simulated HWs
The SAWS definition of a HW was used, where all events with an observed or simulated temperature of 5 °C or greater above the daily average maximum temperature of the hottest month, for a minimum of three consecutive days, are identified, for each station or grid point, respectively. The threshold calculation is made based on present-day climate simulations, and the present-day thresholds are applied to future projections. These thresholds are used by the Disaster Risk Reduction (DRR) division of SAWS which issues HW warnings to the public. The HW thresholds are generally higher across the western parts of the country (which are also arid), ranging from 32 to 44 °C, in line with observed average temperatures shown in Figure 2a during DJF (when most HW are experienced). The thresholds are also relatively high in the north-east lowveld, ranging from 32 to 40 °C.
South Africa experiences a considerable spatial variability of HW occurrences seasonally ( Figure  4). The country experienced higher HW frequency in the northern parts, with Punda Maria (22.68° S, 31.02° E) recording 45 HWs from 1983 to 2012. The north-eastern low-lying regions generally experience relatively high HW occurrences being warm throughout the year.  Observations show that most regions experience HWs during the austral summer, particularly in DJF, with the exception of subtropical western to south-west coastal regions experiencing HWs mainly between May and August. The northern parts of South Africa generally experienced over 12 HWs between September and November (SON) during the period 1983-2012. The western coastal regions of South Africa experienced most HWs in winter when subtropical anticyclones are more equator-ward and when a continental high becomes more established. Descending offshore flow warms dry adiabatically creating a coastal low that precedes the arrival of a cold front during this time. This descending wind leads to berg winds and have been linked to frequent HW occurrences across the south coast of South Africa. Port Nolloth experienced most (71) HW occurrences between 1983 and 2012, 66 of which occurred between April and September. Figure 5a shows the observed average duration of the HWs across the whole of South Africa. The southern parts of the country are observed to experience shorter lived HWs. The extreme northern parts of Limpopo experience HWs that last for six days or more on average. Observations show that most regions experience HWs during the austral summer, particularly in DJF, with the exception of subtropical western to south-west coastal regions experiencing HWs mainly between May and August. The northern parts of South Africa generally experienced over 12 HWs between September and November (SON) during the period 1983-2012. The western coastal regions of South Africa experienced most HWs in winter when subtropical anticyclones are more equator-ward and when a continental high becomes more established. Descending offshore flow warms dry adiabatically creating a coastal low that precedes the arrival of a cold front during this time. This descending wind leads to berg winds and have been linked to frequent HW occurrences across the south coast of South Africa. Port Nolloth experienced most (71) HW occurrences between 1983 and 2012, 66 of which occurred between April and September. Figure 5a shows the observed average duration of the HWs across the whole of South Africa. The southern parts of the country are observed to experience shorter lived HWs. The extreme northern parts of Limpopo experience HWs that last for six days or more on average. Simulated thresholds were calculated for both RCP 4.5 and 8.5 for all ensemble members and were found to exhibit a similar pattern as the observations. Figure 6a,b show examples of these thresholds from one of the members, CCAM forced with CCSM4, for both RCP4.5 and RCP8.5. Simulated HW thresholds are well below 32 °C in much of the eastern parts and coastal regions of the country. The simulated HW thresholds for both RCP 4.5 and 8.5 are spatially consistent with the observed thresholds ( Figure 1). The model shows an east-west and south-north gradient, with the highest thresholds in the north and the west as shown in the observations. The median HW frequency as simulated by CCAM is shown in Figure 6c,d for DJF and JJA, respectively. The number of HWs simulated across northern parts of the country is high for both RCPs, as well as across the interior regions adjacent to the south and east coasts. Figure 4a and b also show that SON and DJF are observed to have more HWs than other seasons. The CCAM model is able to capture the relatively large number of HWs across the Limpopo province. Across the western parts of the country, the largest number of HWs is simulated further west of where the events are observed to peak in frequency. The CCAM model is found to simulate a somewhat larger number of HWs than what is observed in general. It also had challenges simulating winter HWs over the west in the present-day climate. Simulated thresholds were calculated for both RCP 4.5 and 8.5 for all ensemble members and were found to exhibit a similar pattern as the observations. Figure 6a,b show examples of these thresholds from one of the members, CCAM forced with CCSM4, for both RCP4.5 and RCP8.5. Simulated HW thresholds are well below 32 • C in much of the eastern parts and coastal regions of the country. The simulated HW thresholds for both RCP 4.5 and 8.5 are spatially consistent with the observed thresholds ( Figure 1). The model shows an east-west and south-north gradient, with the highest thresholds in the north and the west as shown in the observations. The median HW frequency as simulated by CCAM is shown in Figure 6c,d for DJF and JJA, respectively. The number of HWs simulated across northern parts of the country is high for both RCPs, as well as across the interior regions adjacent to the south and east coasts. Figure 4a and b also show that SON and DJF are observed to have more HWs than other seasons. The CCAM model is able to capture the relatively large number of HWs across the Limpopo province. Across the western parts of the country, the largest number of HWs is simulated further west of where the events are observed to peak in frequency. The CCAM model is found to simulate a somewhat larger number of HWs than what is observed in general. It also had challenges simulating winter HWs over the west in the present-day climate.
Most parts of the country are observed to experience the smallest number of HWs during the winter season with an exception of north-western coast of South Africa where observed HWs are higher. This may be due to prolonged berg winds resulting from the anticyclonic flow by the Atlantic High which is more equator-ward during this time. The eastern part of the country is different, where simulations indicate relatively high number of HWs compared to observations in winter. Most parts of the country are observed to experience the smallest number of HWs during the winter season with an exception of north-western coast of South Africa where observed HWs are higher. This may be due to prolonged berg winds resulting from the anticyclonic flow by the Atlantic High which is more equator-ward during this time. The eastern part of the country is different, where simulations indicate relatively high number of HWs compared to observations in winter.
The CCAM model captures the seasonal differences in the number of HWs well. The smallest number of HWs is simulated to occur in winter, while the highest number is simulated for summer. The number of HWs simulated in spring is higher than those simulated for autumn, however, both autumn and spring HWs are fewer than those that occur in summer, and more than those that occur in winter. The spring simulations are in agreement with observations in that they show the highest number of HWs across the eastern parts of the country. The observations showed a similar seasonal cycle, and hence we can say the CCAM model simulates present-day climate with high skill, and its future projections can be considered with some confidence. The observed average HW duration was shown in Figure 5a. Figure 5b shows the simulated average duration for DJF and JJA for RCP 4.5 during the present-day climate. The south-north gradient is captured in general, with the southern The CCAM model captures the seasonal differences in the number of HWs well. The smallest number of HWs is simulated to occur in winter, while the highest number is simulated for summer. The number of HWs simulated in spring is higher than those simulated for autumn, however, both autumn and spring HWs are fewer than those that occur in summer, and more than those that occur in winter. The spring simulations are in agreement with observations in that they show the highest number of HWs across the eastern parts of the country. The observations showed a similar seasonal cycle, and hence we can say the CCAM model simulates present-day climate with high skill, and its future projections can be considered with some confidence. The observed average HW duration was shown in Figure 5a. Figure 5b shows the simulated average duration for DJF and JJA for RCP 4.5 during the present-day climate. The south-north gradient is captured in general, with the southern parts of the country experiencing shorter lived HWs. However, the model does not capture the location of the longest lived HWs across the northern most parts of the country. It may be noted that the observed duration is based on only 24 stations across the whole of the country, while the peak in the north is a result of interpolating from two stations.
With a warming globe, the expectation is that HW characteristics may not be the same in future climate compared to those occurring in current climate. This is largely due to the fact that average temperatures are increasing, and hence the present-day thresholds for defining a HW may be exceeded easily. We analyzed the projected HWs using the present-day definition, because that is the definition that we experience in our current climate. Very long lasting HWs in the future climate may simply mean the average temperature has gotten higher, and perhaps what is considered as a HW in the current climate will not be considered as a HW in the future climate using future averages. The subsequent sections provide expected HW frequency, duration and intensity in South Africa in the future climates, i.e., 2010 to 2039, 2040 to 2069, and 2070 to 2099.

Period 2010 to 2039
The number of HWs in the immediate future period of 2010 to 2039 is projected to be higher than the HW frequency in 1983-2012. As expected, the number of HW occurrences are higher during DJF (Figure 7a,b) season in 2010-2039, with areas in the east and northern parts of South Africa projected to experience over 80 events over the 30-year period which is about 10 more than the present-day climate frequency. The simulations for 2010-3039 indicate a slight increase in the overall number of HW frequency in South Africa, particularly across the eastern parts of the country during the JJA season (Figure 8a,b). Simulated HWs under the RCP 4.5 and RCP 8.5 scenarios are in agreement with both sets showing an increase in the events in the future climate. The spatial pattern of simulated HW frequency corresponds closely between the two scenarios for the immediate future period. The western parts of South Africa are simulated to experience the smallest number of HWs (less than 10 events in other areas) over this 30-year period during the austral autumn, winter, and spring. These are the parts influenced the most by passing cold fronts which are responsible for the winter rainfall.  The average duration of HWs across South Africa during the 2010-2039 period is projected to increase throughout the year. Regions in the central interior are expected to experience HWs lasting for over a week during SON and DJF. Just like in the present-day-climate, HWs occurring in winter will continue to last for fewer days, i.e., 3-4 days. This is thought to be as a result of the moderating effect by adjacent oceans. HW intensity during this period slightly increases throughout the year compared to the 1983-2012 period. This increase strongly manifests in the 90th percentile under RCP 4.5 and 8.5 in winter with an average increase of about 6 • C/30-year period over the Highveld of Mpumalanga and KwaZulu-Natal (not shown). Moreover, the simulation indicates that HWs will be more intense across the north-western part of South Africa during DJF with a seasonal average ranging from 41 to 44 • C during this period.

Period 2040 to 2069
HW frequencies are projected to be significantly higher during this period than during the present-day climate, particularly across the east. Some regions, such as the high-altitude regions of the Free State and KwaZulu-Natal, which were simulated as not experiencing HWs in winter during present-day climate are beginning to experience such events over the mid-future period. Over 20 events are projected to occur across those areas during austral winter.
Dramatic increases in DJF HW frequency are expected during the 2040-2069 period (Figure 7c,d) when compared to the present climate, particularly along the south and eastern coastal regions (frequency of 30 HW higher over the 30-year period). These regions are expected to experience up to 100 events during summer over the 30-year period; which averages to slightly over three events per summer season per year. However, it has to be noted that these HWs are projected to last for less than a week along the coasts but the number of days per event increases inland. Increases in the number of HWs are also observed in a greater extent across the north-western parts of the country during both summer and autumn and the simulations also indicate that 2040-2069 HWs will have the longest duration, between 9 and 10 days, across this region in summer. This may be associated with increasingly dry conditions across this region, as the region experiences little or no rainfall throughout the year. Drier regions heat up quickly because the heating that reaches the surfaces does not start by evaporating water on the surface, which otherwise result in cooling by latent heat absorption [49] (Figure 9a,b) season as most parts that were experiencing HWs less than 20 events in 1983-2012 in the western half of the country are expected to experience over 40 events which is more than 100% increase.
While projected increases during the 2070-2099 period are slightly higher in all seasons compared to the period 2040-2069, major differences in these two periods are observed most strikingly in the average duration of HWs. HW events across coastal and eastern regions of the country are expected to last for about a week (Figure 9c,d). HWs are projected to become common events across South Africa and much of the country's central interior will experience HWs lasting for over 10 days, particularly during summer. As much as simulations indicate that coastal regions will experience most HWs in summer during this 30-year period, the HWs over the coasts will continue to be the least lasting in the country.
HW intensity is projected to drastically increase from the present-day climate towards the end of the 21st century, in the period 2070-2099. The eastern parts of South Africa are expected to experience HWs of low intensity compared to the western parts during summer seasons of the same period. This may be linked to convective rainfall events in the east during SON-DJF, which may bring rainfall and relief to high atmospheric temperatures. However, it is projected that HWs will be more intense on the east particularly across the northeast reaching an average of over 47 • C/30-year period during the March-May (MAM)-JJA compared to the dry western parts of the country (Figure 10). Atmosphere 2020, 11, x FOR PEER REVIEW 15 of 19

Summary and Conclusions
In this study we have shown that there will be a decrease in the number of cold events, while the hot extreme events are increasing. We used an ensemble of CCAM downscalings of six different GCMs for two different RCPs to simulate present-day climate and project changes through to 2099. The HW definition used in this study does not take into account extremes associated with cooler months and changes in impacts associated with adaptation. The simulations were in agreement with observations that the northern parts of the country experiences the highest temperatures, as indicated by the higher number of HW occurrences. The model was, however, found to simulate a high number of total HW occurrences during the present-day climate than those observed, and did not capture the longest lasting HW across the northern parts of Limpopo province. This study has found that HWs are relatively unusual in the present-day climate but are expected to occur more frequently in the future warmer climate when using both RCP 4.5 and RCP 8.5 emission scenarios, particularly during the far future period, 2070-2099, reaching an over 50% increase in both summer and winter. It is also indicated that HWs are also expected to last longer, particularly across the interior, and become more intense, which is in agreement with earlier findings (e.g., [33]) about the nature of HWs in the 21st century.
It has been indicated that, on the global scale, the magnitudes of the probability of climate extremes occurring are expected to change in both mean and variance in the future. It has previously been suggested that, on average, the world is expected to experience more hot weather events [50]. However, that may not be the case on regional basis. Although several studies have investigated temperature projections across Africa (e.g., [5,26]), few studies (e.g., [51,52]) have investigated the attributes and spatial structure of warm temperature extreme events such as HWs in South Africa. In addition to the nature of HWs in the present climate, this study also investigated temperature trends in South Africa in a future warmer climate. During the 2070-2099 period an average increase of 2 HWs per three-month season per year is expected, which is similar to earlier findings [41] that a median increase of 2.5 events per degree of global warming is expected across southern Africa in the 21st century. Simulations indicated that daily maximum temperatures will continue to rise, increasing by up to 6 • C in most parts of the interior of South Africa throughout the year during the 2070-2099 period. This increases the likelihood of having intense and frequent HW occurrences.
Simulated HW thresholds are consistent with observed HW threshold. It was established in this study that HWs will occur more frequent in future climates, especially during DJF season. The simulations also indicate that HW frequency increases with time and last longer, particularly in the interior of South Africa. Regions that are not prone to HW occurrences in the present-day climate are expected to experience HWs in the future warmer climate.
The country experiences less HW events in mid-winter, and that is also expected to continue in future climates. Winter HW events are expected to extend to the eastern interior and to also increase in intensity reaching over 47 • C/30-year. HWs in South Africa are not only increasing in frequency, but are also expected to last longer and become more intense. The increasing intensity is highly indicated under the RCP 8.5 than 4.5. As in the present-day climate, it is also expected that in future climates HWs will last longer, over two weeks in rare circumstances, across the interior of the country compared to coastal regions. In conclusion, HWs in much of South Africa are expected to occur more frequently, last longer, and become more intense.