Trends in extreme temperature indices in South Africa:

INTERNATIONAL JOURNAL OF CLIMATOLOGY Int. J. Climatol. 33: 661 676 (2013) Published online 12 March 2012 in Wiley Online Library (wileyonlinelibrary.com) DOI: 10.1002/joc.3455 Trends in extreme temperature indices in South Africa: 1962 2009 A. C. Kruger* and S. S. Sekele Climate Service, South African Weather Service, Pretoria, South Africa ABSTRACT: Trends in daily maximum and minimum extreme temperature indices were investigated for 28 weather stations in South Africa, not only for the common period of 1962 2009, but also for longer periods which the individual record lengths of the stations would allow. The utilized weather stations had limited gaps in their time series, did not undergo major moves, or had their exposure compromised during the study period, as to influence the homogeneity of their time series. The indices calculated were forthcoming from those developed by the WMO/CLIVAR Expert Team on Climate Change Detection and Indices (ETCCDI), but only those applicable to the South African climate were selected. The general result is that warm extremes increased and cold extremes decreased for all of the weather stations. The trends however vary on a regional basis, both in magnitude and statistical significance, broadly indicating that the western half, as well as parts of the northeast and east of South Africa, show relatively stronger increases in warm extremes and decreases in cold extremes than elsewhere in the country. These regions coincide to a large degree with the thermal regimes in South Africa which are susceptible to extreme temperatures. The annual absolute maximum and minimum temperatures do not reflect the general trends displayed by the other indices, showing that individual extreme events cannot always be associated with observed long-term climatic trends. The analyses of longer time series than the common period indicate that it is highly likely that warming accelerated since the mid-1960s in South Africa. Copyright 2012 Royal Meteorological Society KEY WORDS temperature trends; temperature extremes; climate change; South Africa Received 30 December 2010; Revised 22 December 2011; Accepted 9 February 2012 1. Introduction According to the latest Intergovernmental Panel on Climate Change (IPCC) report, there is strong scientific evidence that climate change is mostly attributed to human activities (Trenberth et al., 2007). Climate change has resulted in rising temperature trends with associated changes in temperature extremes across the globe (Hansen et al., 2001; Alexander et al., 2006; Brohan et al., 2006; Caesar et al., 2006; Lugina et al., 2005; Smith and Reynolds, 2005). The rate of warming over the last 50 years is almost double that over the last 100 years (0.13 ± 0.03 C vs 0.07 ± 0.02 C per decade). Also, over 74% of the global land area sampled, a significant decrease in the annual occurrence of cold nights (i.e. extremely low minimum temperatures) is shown, while a significant increase in the annual occurrence of warm nights (i.e. extremely high minimum temperatures) took place over 73% of the area. In the TAR (IPCC, 2001), global surface temperature trends were examined for three sub-periods: 1910 1945, 1946 1975 and 1976 2000. The first and third subperiods had rising temperatures, while the second subperiod had relatively stable global mean temperatures * Correspondence to: A. C. Kruger, South African Weather Service, Private Bag X097, Pretoria 0001, South Africa. E-mail: Andries.Kruger@weathersa.co.za (Trenberth et al., 2007). Therefore, since the instrumental period in temperature observations began, the global trend in temperatures was not near-constant, but consisted of periods of stronger and weaker trends. The inclusion of as many regions as possible in assessments of the global trend is essential, as the trends are highly variable, not only on a temporal but also on a regional basis. Examples of the most recent trend studies over parts of southern Africa are by Kruger and Shongwe (2004) who, in their temperature trends study for 1960 2003 in South Africa, found in general positive trends in the annual mean, maximum and minimum temperatures, as well as increases in days and nights with high temperatures and decreases in days and nights with low temperatures. However, the observed warming was not consistent between the analysed weather stations, which indicated regional variations in temperature trends. New et al. (2006), reporting on a regional workshop on temperature trends, also indicated warming trends for most of the southern African subcontinent, with magnitudes similar to those found by Kruger and Shongwe (2004). The aim of this study is to update the state of South African trends of daily maximum and minimum temperatures extremes for the period 1962 2009, and longer where possible, which will be useful in the assessment of climate change in the region. Copyright 2012 Royal Meteorological Society

662 A. C. KRUGER AND S. S. SEKELE 2. Data A common analysis period of 1962 2009 of daily maximum and minimum temperature data from the SAWS (South African Weather Service) climate database was selected, in order to obtain the longest possible period with a reasonable number of available climate stations to cover most regions in South Africa. The data sets of the selected stations were subjected to quality control to remove any values which were possible erroneous. As a first step, the metadata files of the potential climate stations were scrutinized for large movements, inadequate exposure, as well as maintenance problems that could have resulted in inhomogeneities in their data sets. The validity of anomalously high or low values were checked against the reports of the prevailing weather conditions around the time of observation, and also compared to the values measured at available neighbouring stations sharing the same climate regime. Consequently, any suspicious data were removed from the time series. Thereafter, the completeness of the data series was verified to be higher than 90%. A total of 28 weather stations were accepted for analysis. The SAWS climate number, station name, location in terms of latitude and longitude, altitude, mean minimum and maximum temperature, as well as the available period of record, are presented for each station in Table I. The spatial distribution of the weather stations, as well as the provinces of South Africa for reference purposes, is shown in Figure 1. 3. Methodologies 3.1. Trend analysis The trend analysis of the extreme temperature indices were performed using the RClimDex software, which is available from the WMO/CLIVAR/JCOMM ETCCDI website http://cccma.seos.uvic.ca/etccdi. This software was used in 19 regional workshops worldwide, e.g. for southern Africa (New et al., 2006) and most recently Indonesia in 2009 (workshop report available on the ETCCDI website), in an attempt to address the scarcity of climate trend studies over some regions of the world. RClimDex is capable of computing 27 core indices. However, only the indices which could be relevant to South Africa were selected, and are shown in Table II. In particular, the threshold indices were not included, due to the highly variable climate of South Africa, mostly because of the topography, as can be deduced from the altitudes of the weather stations in Table I. Table I. List of weather stations utilized in the study. Climate number Station name Latitude ( S) Longitude ( E) Mean Minimum Temperature ( C) Mean Maximum Temperature ( C) Altitude (m) Period 0 003 020 Cape Agulhas 34.19 22.13 14.3 20.0 8 1911 2009 0 012 221 Mossel Bay 33.97 18.60 14.6 20.7 60 1920 2009 0 021 178 Cape Town 26.95 24.63 11.4 22.0 46 1957 2009 0 035 209 Port Elizabeth 29.97 30.95 13.5 22.3 59 1937 2009 0 050 887 Willowmore 29.02 29.87 9.1 23.2 840 1960 2009 0 059 572 East London 24.98 31.60 14.4 22.9 124 1953 2009 0 106 880 Vredendal 27.27 28.50 11.1 25.8 37 1961 2009 0 127 272 Mthatha 29.10 26.30 10.9 24.1 743 1959 2009 0 169 880 De Aar 27.37 29.88 8.6 24.7 1287 1959 2009 0 193 561 Vanwyksvlei 23.87 29.45 9.7 26.7 962 1939 2009 0 240 808 Durban 33.03 27.83 16.5 25.2 8 1957 2009 0 242 644 Port Nolloth 30.35 21.82 10.8 18.9 4 1960 2009 0 247 668 Pofadder 30.67 24.00 11.4 25.8 989 1941 2009 0 261 516 Bloemfontein 33.30 23.48 7.5 24.4 1354 1962 2009 0 290 468 Kimberley 26.47 20.61 10.9 26.0 1197 1932 2009 0 300 690 Estcourt 29.25 16.87 10.0 23.9 1148 1957 2009 0 317 475 Upington 29.13 19.39 12.5 28.5 841 1952 2009 0 399 894 Bothaville 25.73 28.18 9.7 25.8 1280 1961 2009 0 403 886 Frankfort 28.80 24.77 6.8 23.7 1500 1954 2009 0 406 682 Volksrust 25.50 26.35 8.1 21.7 1652 1954 2009 0 432 237 Vryburg 27.41 26.50 9.5 26.7 1234 1920 2009 0 461 208 Twee Rivieren 28.41 21.26 11.0 29.4 879 1959 2009 0 513 284 Pretoria PUR 34.83 20.01 9.6 25.3 1286 1959 2009 0 546 630 Marico 31.53 28.67 11.7 27.4 1078 1936 2009 0 554 816 Lydenburg 31.67 18.50 9.5 22.9 1434 1957 2009 0 596 179 Skukuza 25.11 30.48 14.3 29.5 263 1951 2009 0 677 802 Polokwane 33.98 25.61 11.7 24.7 1237 1953 2009 0 809 706 Musina 22.27 29.90 15.3 29.6 525 1934 2009 The period indicates the longest period of record, unique for each weather station, for which index trends could be calculated.

TRENDS IN TEMPERATURE EXTREMES IN SOUTH AFRICA 663 Musina Twee Rivieren Polokwane Limpopo Lydenburg Skukuza Marico Pretoria PUR North-West Gauteng Mpumalanga Vryburg Bothaville Frankfort Volksrust Upington Kimberley Port Nolloth Pofadder Bloemfontein Estcourt Northern Cape Durban Van Wyksvlei De Aar Vredendal Mthatha Cape Town Free State Eastern Cape East London Willowmore Port Elizabeth Mossel Bay Cape Agulhas Western Cape KwaZulu-Natal Figure 1. Map of South Africa with provinces, depicting positions of 28 weather stations utilized in the study. Table II. Temperature indices covered in the study. Index Description Units TX90P Annual percentage of days when TX > 90th percentile % TX10P Annual percentage of days when TX < 10th percentile % TXx Annual maximum value of TX C TXn Annual minimum value of TX C WSDI Annual count of days with at least 6 consecutive days when TX > 90th percentile d TNx Annual maximum value of TN C TNn Annual minimum value of TN C TN90P Annual percentage of days when TN > 90th percentile % TN10P Annual percentage of days when TN < 10th percentile % CSDI Annual count of days with at least six consecutive days when TN < 10th percentile d TX indicates daily maximum temperature and TN indicates daily minimum temperature. Percentiles are based on the 1971 2000 period. The base line period of 1971 2000 was used in the analysis, initially for the estimation of threshold values for the identification of anomalous temperature values for quality control purposes, but also the thresholds for the percentile-based indices. Linear trends were calculated with the least-squares method for each index and the correlation factors tested for significance with the t-test at the 95% level of confidence (Wilks, 2006). In addition, the error bars of the trends are included in the results, as provided by the outputs of the RClimDex software, also at the 95% confidence level. 3.2. Cluster analysis Cluster analysis is often used in climatological studies to define regions with similar climatological characteristics, and was performed on the index trend values for this particular purpose. Of the different cluster analysis techniques the most widely applied method is the K-means method, as it is relatively simple to use and also allows reassignment of observations as the analysis proceeds from one number of clusters to the next. The K refers to the number of groups or clusters, which is specified in advance of the analysis. The K-means algorithm usually begins with a random partition of the n data vectors into the pre-specified number of groups. The algorithm proceeds then as follows: 1. Compute the vector means, i.e. x k, k = 1,...,K;for each cluster. 2. Calculate the Euclidian distances between the current data vector x i and each of the K x k s.

664 A. C. KRUGER AND S. S. SEKELE Figure 2. Regions A F from cluster analysis performed on annual mean minimum and maximum temperatures, as well as geographical coordinates of stations. 3. If necessary the x i is reassigned to the group whose mean is closest. Repeat for all x i,i = 1... n. 4. Return to step 1. The algorithm is iterated until a full cycle through all the data vectors produces no reassignments (Wilks, 2006). The K-means method was performed on the trend results for each index, together with the geographical coordinates of the weather stations. The groupings of stations according to the maximum number of clusters that could be resolved are depicted in the maps of the results, to aid in the more objective interpretation of the regional differences thereof. 3.2.1. Identification of thermal regimes As an example of the application of cluster analysis, and to interpret the trend results presented in the following section in the context of the thermal regimes of South Africa, cluster analysis was applied to the annual means of the minimum and maximum temperatures, of which the result is shown in Figure 2. Six clusters, A F, were resolved, with the stations grouped broadly according to the general characteristics of the mean temperatures in South Africa. These clusters can broadly be considered to be homogeneous groups of stations, which exhibit roughly similar annual thermal characteristics or regimes. The clusters or groups of stations identified can be summarized as follows, with reference to Figure 2. Coastal: cluster A represents the western and southern coastal region, which exhibits a generally mild Mediterranean climate, while cluster B represents the southeast and eastern coast, which has a subtropical climate with relatively higher temperatures than A. Interior: cluster C represents the Lowveld, a region in the northeast at low altitude with relatively high mean temperatures throughout the year. Cluster D represents the dry western interior, which is characterized by temperature extremes with relatively high maximum temperatures in summer, usually exceeding 30 C, but low temperatures in winter, when minimum temperatures often drop below freezing. The diurnal range in temperature of cluster D is also higher than elsewhere. The southeastern interior, represented by cluster E, exhibits relatively small diurnal ranges in temperature, most probably due to the moderating influence of the regular influx of oceanic air from the Indian Ocean. However, the northwestern part of the region on the plateau at high altitude, experiences very low temperatures in winter. If it was possible to identify a larger number of clusters, a separate group of stations around the escarpment would probably have been identified. Cluster F in the north has higher temperatures during winter than the larger part of cluster E, and is not as prone to temperature extremes as cluster D. 4. Results and discussions The index trend results are presented in Appendices A and B, for the maximum and minimum temperature indices, respectively. These results, as well as their delineation with cluster analysis, are shown in the maps

TRENDS IN TEMPERATURE EXTREMES IN SOUTH AFRICA 665 1.98* 2.29* 1.81* 1.61* 1.44* 3.4* 2.8* 2.5* 0.41 1.36* 1.82* 0.41 1.79* -0.43 1.82* 1.12 1.9* 0.89* 1.99* 0.52 1.26* 0.43 0.61* 3.94* 2* 0.92* 0.68* 0.83* Figure 3. Trends in annual percentage of days of TX90P/decade. The demarcating lines refer to regions of similar trend magnitude and indicates trends which are statistically significant at the 5% level. -1.11* -0.71* -2.03* -2.47* -1.3* -0.93* 0.08-0.86* -0.59* -1.32* -0.79** -1.1* 0.22-0.24-0.96* -0.76-0.64* -0.22-0.58-0.65* -0.61* -0.5-1.12* -1.6* -1.93* -0.1-0.94-1.09* Figure 4. Trends in annual percentage of days of TX10P/decade. The demarcating lines refer to regions of similar trend magnitude and indicates trends which are statistically significant at the 5% level. of Figures 3 12. The clusters with more noteworthy results, which are highlighted in the discussions either because of their spatial extent or magnitudes of trends, are identified with capital letters in the relevant maps. In the discussions of the results reference is made to the thermal regimes over South Africa, as identified in Figure 2, as it is important to note the extent of how these regimes were affected by the observed trends.

666 A. C. KRUGER AND S. S. SEKELE 0.2 0.18 0.17 0.08-0.02 0.01 0.42* -0.14-0.03 0.14-0.03-0.09 0.8* 0.2 0.26* 0.33* -0.03 0.22 0.25-0.04 0.01 0.22-0.12 0.41-0.15 0.14 0.19-0.07 Figure 5. Trends in TXx/decade. The demarcating lines refer to regions of similar trend magnitude and indicates trends which are statistically significant at the 5% level. 0.19 0.11 0.21* 0.19 0.39* 0.44* 0.33 0.24 0.26 0.28* 0.38* 0.14 0.02 0.45 0.36 0.05 0.43 0.1 0.7* 0.67* 0.19 0.03 0.11 0.21* 0.32* 0.08 0.06 0.01 Figure 6. Trends in TXn/decade. The demarcating lines refers to regions of similar trend magnitude and indicates trends which are statistically significant at the 5% level. 4.1. Maximum temperature-related trends for 1962 2009 The trends in the TX90P index (percentage of days per year when the maximum temperature is greater than the 90th percentile of the 1971 2000 base period) is shown in Figure 3. Most of the weather stations, 22 out of 28 (79%), experienced increases in warm extremes which are statistically significant at the 5% level, with the

TRENDS IN TEMPERATURE EXTREMES IN SOUTH AFRICA 667 1.47* 0.65 6.74* 1.55 0.28 0.21 1.95* 0.10 1.19 0.75-0.10 1.33 0.15 0.71 1.65* 0.78 0.70 1.75 2.64 0.29 0.18 2.84 0.82 0.21 Figure 7. Trends in annual count of days with at least six consecutive days of TX90P/decade. The demarcating lines refer to regions of similar trend magnitude and indicates trends which are statistically significant at the 5% level. denotes stations with no results available, due to too many zeros in the time series. 1.52* 0.78 1.14* -0.43 0.12 1.25* -0.03 0.1 0.95 1.92* 1.1* 0.98 0.66* 1.14* 0.47 0.44 1.98* -1.1* -0.79 1.17* 0.37 0.63 1.03* 1.69* 0.97* 0.4-0.53 0.4 Figure 8. Trends in annual percentage of days of TN90P/decade. The demarcating lines refer to regions of similar trend magnitude and indicates trends which are statistically significant at the 5% level. strongest trends observed along the southern coastline (cluster A with a mean of 2.97% per decade). Other regions which exhibit relatively strong positive trends are the larger part of the Northern Cape province (cluster B with a mean of 2.26% per decade), and parts of the central and northern interior (cluster C with a mean

668 A. C. KRUGER AND S. S. SEKELE -1.38* -1.06* -1.54* 0.21-0.86-3.93* 0.52 0.8-1.14* -2.14*-1.31* -0.73-1.58* -1.47* -1.51* -1.16* -0.51 1.37* -0.16-2.02* -0.29-2.67* -2.18* -1.25* -0.45-0.93* 0.92-0.52 Figure 9. Trends in annual percentage of days of TN10P/decade. The demarcating lines refer to regions of similar trend magnitude and indicates trends which are statistically significant at the 5% level. 0.41* -0.07 0.27 0.05-0.01 0.16-0.04 0.03 0.23 0.18* -0.07-0.34 0.61* 0.3* 0.11 0.18-0.05-0.15-0.03-0.06 0.09 0.04 0.12 0.1 0.03-0.02-0.01-0.07 Figure 10. Trends in TNx/decade. The demarcating lines refer to regions of similar trend magnitude and indicates trends which are statistically significant at the 5% level. of 1.87% per decade). In summary, it is the western half and the northern interior, mostly with relatively warmer thermal regimes, mostly in clusters D and C in Figure 2, which exhibited the strongest trends. In other parts of the country the trends were decidedly weaker. TX10P (percentage of days per year when the maximum temperature is less than the 10th percentile) shows,

TRENDS IN TEMPERATURE EXTREMES IN SOUTH AFRICA 669 0.74* 0.17 0.37 0.05 0.14 0.68* -0.13-0.07 0.62* 0.45* 0.21* 0.08 0.02 0.07 0.22 0.24 0.34* -0.56* -0.07 0.62* 0.14 0.4* 0.4* 0.29 0.05 0.24* -0.09 0.16 Figure 11. Trends in TNn/decade. The demarcating lines refer to regions of similar trend magnitude and indicates trends which are statistically significant at the 5% level. 0.17-0.37 2.15-0.26-0.95* -0.69-0.38-0.3 0.52 0.2-1.65* 0.69-1.8* -4.08* 1.09 1.56* 0.06 0.02 0.28 1.24-0.32 Figure 12. Trends in annual count of days with at least six consecutive days of TN10P/decade. The demarcating lines refer to regions of similar trend magnitude and indicates trend which are statistically significant at the 5% level. denotes stations with no results available, due to too many zeros in the time series. as can be expected, some similarities in the spatial distribution of the results of TX90P but with trends of the opposite sign: 18 of the stations (64%) have pairings of significantly positive and negative trends for TX90P and TX10P, respectively. A total of 19 weather stations (68%) show significantly negative trends for TX10P, as shown

670 A. C. KRUGER AND S. S. SEKELE in Figure 4. The strongest trends are found in the west in cluster A with a mean of 1.57% per decade, which coincides more or less with clusters A and B in Figure 3. Results for TXx (absolute annual maximum of the daily maximum temperature) shown in Figure 5, indicate significantly positive trends for only four weather stations (14%), which are located in the western, northwestern and central interior in clusters A, B and C. The highest trend of 0.80 C per decade is found at Vredendal in the Western Cape province (the only member of cluster A). Clusters A, B and C mostly cover the region in South Africa with the highest summer temperatures (clusters D and F in Figure 2). In general, the results show very small trend values, which indicate that in the larger part of the country no significant trend in TXx is evident. For TXn (absolute annual minimum of the maximum temperature) similar results than TXx was found, in which the number of significant trends, in this case increases at nine stations (32%), and are relatively small compared to the other maximum temperature indices. The strongest trends of 0.70 and 0.67 C per decade are found in the northeast in the Lowveld (cluster C in Figure 2). The other significant trends are widespread, and do not show higher concentrations in particular regions. Warm spell duration index (WSDI annual count of days with at least six consecutive days when maximum temperature is greater than the 90th percentile) shows only four stations (14%) with significant increases, with the highest trend in the northern part of the Northern Cape at 6.74 d per decade, as indicated by the sole member of cluster A in Figure 7. Due to the small number of significant results, no regional trends can be deduced, except that the four stations are all located in the northern half of the country, where extreme temperatures are in general more prevalent than elsewhere in the country (cluster D in Figure 2). For most of the weather stations along the coast, no trends in WSDI could be calculated, due to the small number of years where cases of six or more consecutive days occurred when the daily maximum temperatures were greater than the 90th percentile. 4.2. Minimum temperature-related indices for 1962 2009 TN90P (percentage of days per year when the minimum temperature is greater than the 90th percentile) shows 13 weather stations (46%) with significantly positive trends, as shown in Figure 8. While half of the significant results are found in the southwest and west in cluster A, which has a mean trend of 1.05% per decade, stations with significant trends are also well represented in the northeast in cluster B, with a mean trend of 1.09% per decade. The strongest positive trends are therefore mostly confined to the regions with thermal regimes exhibiting relatively higher temperatures (clusters C and D in Figure 2), but also the western and southern coastal areas, represented by cluster A in Figure 2. Three stations indicate negative trends, situated in North West, Free State and southeastern Eastern Cape provinces, but only statistically significant at Kimberley in cluster C. TN10P (percentage of days per year when the minimum temperature is lesser than the 10th percentile) shows 17 weather stations (61%) with significantly negative trends, as shown in Figure 9. The strongest negative trends are found in the north and east (cluster A with a mean of 2.07% per decade), the southwestern Cape (cluster B with a mean of 1.67% per decade), as well as most of the Northern Cape province (cluster C with a mean of 1.28% per decade). As with the comparison of results of TX90P and TX10P, TN90P and TN10P show similarities in trend results, but with opposite signs. Eleven weather stations (39%) show pairings of significantly positive and negative trends for TN90P and TN10P, respectively. This number is markedly lower than the results for the comparison of TX90P and TX10P (64%), which indicates that the positive shifts in the overall distributions of daily maximum temperatures are more pronounced than the daily minimum temperatures. However, the regions that exhibit significant changes in both the left and right tails of the maximum and minimum temperature distributions coincide more or less, and cover the western half of the country, as well as the northern and northeastern interior. Again, as with most of the previous indices discussed, the regions where these changes occurred coincide with the thermal regimes which exhibit in general higher temperatures. TNx (absolute annual maximum of the daily minimum temperature) shown in Figure 10, shows significantly positive trends for only four weather stations (14%). Two of the weather stations with significant trends are situated in the west of the country in cluster A, while the other two are in the east in clusters B and C, respectively. There are no apparent regional trends as in most of the country the trends are almost non-existent. Trends in TNn (absolute annual minimum of the daily minimum temperature) shown in Figure 11, show ten stations (36%) with statistically significant increases. All five of the stations in cluster A in the northeast indicate significantly positive trends, with a mean value of 0.54 C per decade. Other clusters in which significantly positive trends are found are B and C. In summary, significant positive trends in TNn are evident in the northeast and east, and the southern parts of the western interior. However, while most of the significantly positive results coincide with clusters C and D in Figure 2, other thermal regimes also experienced some significant trends in TNn. CSDI (annual count of days with at least six consecutive days when the minimum temperature is lesser than the 10th percentile), shown in Figure 12, shows five stations (18%) with significantly negative trends. The strongest negative trend of 4.08 d per decade was found in cluster A, in which all three stations showed significantly negative trends with a mean of 3.77 d per decade. It therefore seems that significant decreases in cold spells are confined to the northeast. Cluster B indicates an extensive region with positive trends, although the results were only significant for one station. For many of the climate stations in the south, southeast and east of

TRENDS IN TEMPERATURE EXTREMES IN SOUTH AFRICA 671 the country, the trends could not be calculated due to the insufficient number of available years with occurrences where the minimum temperature was lower than the 10th percentile for six consecutive days. 4.3. Summary and integration of trend results for 1962 2009 To obtain a characterization of the general trends in extreme temperatures in South Africa for the common period of 1962 2009, the trend results of related indices were compared and integrated to obtain a condensed view or summary of the regional trends and their relative differences over South Africa. 4.3.1. Maximum temperature As discussed in the previous section, if there are significant trends in the daily maximum temperature, there will most likely be a general agreement between the trend results of TX90P and TX10P, albeit of opposite signs. Such trends would indicate a long-term shift in the statistical distribution of maximum temperatures, provided that the variance of the annual maximum temperatures stays near-constant throughout the relevant period. The results shown in Figures 3 and 4 confirm this to a large degree, and also indicate a general increase in maximum temperatures over South Africa. Closer inspection of the results of TX90P and TX10P reveal that the western half and the northeastern interior of the country experienced relatively stronger increases in daily maximum temperatures than elsewhere. The stronger increases in warm extremes broadly cover the regions in Figure 2 indicated by clusters A and D in the west, and C in the northeast. The results for TXx and TXn do not broadly agree with those of the previously mentioned indices, with relatively small percentages of weather stations with significant trends, and little spatial coherence between the results. However, it is worth mentioning that most of the significant results are in the western half of the country, which relates to the relatively strong warming there, as indicated by the results for TX90P and TX10P. Trends in the WSDI would indicate possible changes in the number of occurrences of extended periods with high temperatures, such as heat waves, and would not necessarily coincide with those of the other maximum temperature indices, apart from possibly confirming a general increase in maximum temperatures. The results show that the only regions where there are signs of significant increases in heat waves are the extreme northern parts of the western and northeastern interior, which are in thermal regimes prone to extremes in maximum temperatures, especially during the warm summer months. 4.3.2. Minimum temperature The trend results of TN90P and TN10P indicate a general increase in daily minimum temperatures across South Africa, except for parts of the central interior. Relatively stronger increases in daily minimum temperatures were observed in the west, northeast and east of the country, which coincides with the results of TX90P and TX10P, and ultimately the warmer thermal regimes in the country, represented by clusters C and D in Figure 2. The results of TNx and TNn are in some agreement with the regional results of TN90P and TN10P. Trends in CSDI would indicate changes in the frequencies of cold spells, when the daily minimum temperatures are much lower than usual for an extended period of time. The results of CSDI do not necessarily have to reflect the results of the other minimum temperature indices, apart from possibly confirming the general trends in the daily minimum temperature. The regions where significant decreases in cold spells are detected are in the northeastern interior and northern part of the western interior, i.e. parts of larger regions of relatively strong positive trends of TN10P. In summary, the trend results for most of the minimum temperature indices indicate that the western half, as well as the northeastern interior of South Africa experienced relatively stronger increases in minimum temperatures than elsewhere in the country. 4.4. Comparisons between trends of 1962 2009 and longer periods All of the weather stations, except one, have record lengths which are longer than the common analysis period of 1962 2009, as indicated in Table I. Analyses of the index trends for longer periods can provide information on the persistence of the observed trends, especially those that are relatively strong and statistically significant for 1962 2009. It should be considered that the results of trend analysis of variables, with magnitudes which are cyclical in nature, depend heavily on the analysis period. In climate analyses, this is particularly relevant to cyclical behaviour which is near-decadal, where non-existent long-term trends might be inferred if the analysis period spans a small number of decades, from a period with relatively high (low) values to a period with relatively low (high) values in the cycle. In such cases, erroneous or exaggerated long-term trends will be indicated by the trend analysis. An indication of persistence of trend over periods which are much longer than the common 1962 2009 analysis period will increase confidence in the results found so far. Appendixes A and B present the results of trend analysis for both the common period of 1962 2009, as well as for the extended periods P, which are unique for the weather stations utilized in the study. Of interest is a comparison of the trend results for those weather stations, where P is much longer than the 1962 2009 period. These are Cape Agulhas, Mossel Bay and Port Elizabeth on the South Coast, Vanwyksvlei and Kimberley in the Northern Cape, Vryburg and Marico in the North West, and Musina in the Limpopo province. For the maximum temperature indices in Appendix A it is observed that in the south, for Cape Agulhas and Mossel Bay, the increase in warm extremes are noticeably

672 A. C. KRUGER AND S. S. SEKELE Figure 13. RClimDex output of trend analysis of TX90P for Cape Agulhas, for the period 1911 2009. Thin line and circles indicate time series and index values, bold line indicates linear trend and dotted line indicates decadal-scale variations based on Lowess smoother (Cleveland, 1979). weaker over the longer term, with the trends for TX90P and TX10P non-significant for Mossel Bay over the longer period of 1920 2009. Closer inspection of the results reveal an opposite trend from the 1920s to the 1960s, which offsets the increase in warm extremes since the 1960s to 2009. For Cape Agulhas, the differences in trends are less pronounced, but here TX90P shows almost no trend from 1911 to 1965, as shown by the RClimDex output presented in Figure 13. For Port Elizabeth to the east, the differences in trend between the shorter and longer term are much smaller. The stations in the Northern Cape show similar results than those in the south, with weaker trends over the longer analysis periods. For Kimberley, almost no trends were observed for TX90P and TX10P from 1911 to 1975, after which strong increases in warm extremes were observed. Similar weak trends in TX90P and TX10P were observed for Vanwyksvlei, but for the shorter period from 1932 to 1975. In the North West province of Vryburg, warm extremes decreased from 1920 to the mid 1960s, after which a significant increase of said extremes occurred. For Marico, the trend results show the opposite, especially for TX90P where the trend is stronger over the longer period of 1936 2009. A relatively weak increase in warm extremes occurred from the early 1980s to late 1990s, compared to the whole period of 1936 2009. In the extreme north, Musina shows similar results than most of the longer term weather stations, with a small decrease in warm extremes over the period 1934 1965, and a strong increase thereafter. Trends in minimum temperatures over the longer periods, as shown in Appendix B, indicate trends in TN90P along the South Coast comparable to that of 1962 2009, but stronger negative trends in TN10P. In the Northern Cape, the longer period analyses indicate accelerated decreases in cold extremes since the mid- 1960s, similar to the results for the maximum temperature indices, which indicated increases in warm extremes. The weather stations in the North West province show similar trends for the shorter and longer periods. The indices for Musina in the extreme north show, similar to most of the longer term weather stations, an accelerating decrease in cold extremes since the mid-1960s. 5. Summary The study provided an updated analysis of the daily maximum and minimum temperature trends of relevant extreme temperature indices over South Africa, for the period 1962 2009. While the maximum temperature indices show general increases in warm extremes, the minimum temperature indices show general decreases in cold extremes. This indicates that South Africa experienced general warming over the analysis period [also see Kruger et al. (2011) for an update on the trends in annual mean temperatures over South Africa]. Most of the results indicate relatively stronger increases in warm extremes and decreases in cold extremes in the western, northeastern and extreme eastern parts of the country. The results obtained are in general agreement with those of recent temperature trend studies for the region, which show a general warming trend, but with relatively weaker trends in the central parts of South Africa (Kruger and Shongwe, 2004; New et al., 2006). The parts of South Africa that experienced relatively stronger warming can be summarized as in Figure 14. The regions in South Africa with relatively warmer thermal regimes, and which are more prone to hot daily extremes, i.e. the Lowveld in the northeast of the country, the east coast, and the dry western interior as indicated by

TRENDS IN TEMPERATURE EXTREMES IN SOUTH AFRICA 673 Stronger warming Stronger warming Figure 14. Summary of regions of relatively stronger warming in South Africa over the period 1962 2009. clusters B, C and D in Figure 2, experienced the strongest increases in warm extremes. It is envisaged that a persistence in the strong warming observed, particularly in the Northern Cape and parts of the Western Cape, both in the west where the interior can be described as semi-arid with highly variable precipitation, will have a negative effect on the biodiversity, due to habitat loss, and agriculture, due to likely increases in evaporation and consequent heat stress to livestock. Biodiversity has already been affected by rising temperatures in the drier regions of the northern and Western Cape provinces, as evidenced by Foden et al. (2007). This remarkable differential warming over South Africa can most likely be attributed to possible changes in the atmospheric circulation over the subcontinent. Over the western parts it may include possible changes in the strengths of cold fronts moving over the subcontinent from the west, or weaker ridging by the quasi-stationary Atlantic Ocean high pressure system from the south or southeast, especially during the austral summer. In the east, weaker ridging by the Indian Ocean high pressure system might reduce the frequency or strength of the influx of cooler maritime air from the east. It is recommended that the possible changes be investigated with regional model studies and/or the analysis of long-term reanalysis data, both of which falls beyond the scope of the present study. The analyses of longer time series than the common study period of 1962 2009, indicate that for most of the longer term stations which show relatively large differences between trends over the longer term and 1962 2009, the frequencies of warm extremes have accelerated since around the mid-1960s. This finding is in agreement with the mean global temperature trend, where increased warming is evident since the latter part of the 20th century, particularly from the mid-1960s (Hansen et al., 2001; Lugina et al., 2005; Smith and Reynolds, 2005; Brohan et al., 2006). References Alexander LV, Zhang X, Peterson TC, Ceaser J, Gleason B, Klein Tank AMG, Haylock M, Collins D, Trewin B, Rahimzadeh F, Tagipour A, Rupa Kumar K, Revadekar J, Griffiths G, Vincent L, Stephenson DB, Burn J, Aguilar E, Brunet M, Taylor M, New M, Zhai P, Rusticucci M, Vazquez-Aguirre JL. 2006. Global observed changes in daily climate extremes of temperature and precipitation. Journal of Geophysical Research 111: D05109, DOI: 10.1029/2005JD006290. Brohan P, Kennedy JJ, Harris I, Tett SFB, Jones PD. 2006. Uncertainty estimates in regional and global observed temperature changes: a new dataset from 1850. Journal of Geophysical Research 111: D12106, DOI: 10.1029/2005JD006548. Caesar J, Alexander L, Vose R. 2006. Large-scale changes in observed daily maximum and minimum temperatures: creation and analysis of a new gridded data set. Journal of Geophysical Research 111: D05101, DOI: 10.1029/2005JD006280. Cleveland WS. 1979. Robust locally weighted regression and smoothing scatterplots. Journal of the American Statistical Association 74: 829 836. Foden W, Midgley GF, Hughes G, Bond WJ, Thuiller W, Hoffman MT, Kaleme P, Underhill LG. 2007. A changing climate is eroding the geographical range of the Namib Desert tree Aloe through population declines and dispersal lags. Diversity and Distributions 13(5): 645 653, DOI: 10.1111. Hansen J, Easterling D, Imhoff M, Karl T, Lawrence W, Peterson T, Ruedy R, Sato M. 2001. A closer look at United States and global surface temperature change. Journal of Geophysical Research 106: 23947 23963. IPCC. 2001. Climate Change 2001: The Scientific Basis. Contribution of Working Group I to the Third Assessment Report of the Intergovernmental Panel on Climate Change, Houghton JT, Ding Y, Griggs DJ, Noguer M, Van de Linden PJ, Dai X, Maskell K, Johnson CA. (eds). Cambridge University Press: Cambridge, United Kingdom and New York, NY, USA.

674 A. C. KRUGER AND S. S. SEKELE Kruger AC, McBride C, Thiaw WM. 2011. Southern Africa Regional Climate in State of the Climate 2010. Bulletin of the American Meteorological Society 92(6): S161 S163. Kruger AC, Shongwe S. 2004. Temperature trends in South Africa: 1960 2003. International Journal of Climatology 24: 1929 1945. Lugina KM, Groisman PY, Vinnikov KY, Koknaeva VV, Speranskayaet NA. 2005. Monthly surface air temperature time series area-averaged over the 30-degree latitudinal belts of the globe, 1881 2004. In Trends: A Compendium of Data on Global Change. Carbon Dioxide Information Analysis Center, Oak Ridge National Laboratory, US Department of Energy: Oak Ridge, TN, USA. New M, Hewitson B, Stephenson DB, Tsiga A, Kruger A, Manhique A, Gomez B, Coelho CAS, Masisi DN, Kululanga E, Mbambalala E, Adesina F, Saleh H, Kanyanga J, Adosi J, Bulane L, Fortunata L, Mdoka ML, Lajoie R. 2006. Evidence of trends in daily extremes over southern and west Africa. Journal of Geophysical Research 111: D14102, DOI: 10.1029/2005JD006289. Smith TM, Reynolds RW. 2005. A global merged land and sea surface temperature reconstruction based on historical observations (1880 1997). Journal of Climate 18: 2021 2036. Trenberth KE, Jones PD, Ambenje P, Bojariu R, Easterling D, Klein Tank A, Parker D, Rahimzadeh F, Renwick JA, Rusticucci M, Soden B, Zhai P. 2007. Observations: surface and atmospheric climate change. In Climate Change 2007: The Physical Science Basis. Contribution of Working Group I to the Fourth Assessment Report of the Intergovernmental Panel on Climate Change 2001, Solomon S, Qin D, Manning M, Chen Z, Marquis M, Averyt KB, Tignor M, Miller HL(eds). Cambridge University Press: Cambridge, United Kingdom and New York, NY, USA. Wilks DS. 2006. Statistical methods in the Atmospheric Sciences, Elsevier Academic Press: Burlington, MA, USA.

TRENDS IN TEMPERATURE EXTREMES IN SOUTH AFRICA 675 Appendix A. Maximum temperature index trends for the period 1962 2009 Climate No. P Maximum temperature indices TXx TXn TX90P TX10P WSDI TXx (P) TXn (P) TX90P (P) TX10 (P) WSDI (P) 0 003 020 1911 0.041 (0.024) 0.021 (0.010) 0.394 (0.062) 0.160 (0.040) 0.284 (0.152) 0.001 (0.009) 0.021 (0.016) 0.129 (0.019) 0.147 (0.018) 0.078 (0.029) 0 012 221 1920 0.015 (0.028) 0.032 (0.014) 0.200 (0.033) 0.193 (0.040) 0.082 (0.041) 0.005 (0.012) 0.012 (0.005) 0.023 (0.020) 0.001 (0.020) 0.025 (0.029) 0 021 178 1957 0.012 (0.019) 0.011 (0.010) 0.061 (0.024) 0.112 (0.030) 0.003 (0.017) 0.004 (0.009) 0.062 (0.007) 0.120 (0.027) 0 035 209 1937 0.019 (0.021) 0.006 (0.013) 0.068 (0.024) 0.094 (0.024) 0.006 (0.012) 0.011 (0.007) 0.038 (0.014) 0.074 (0.017) 0.020 (0.018) 0 050 887 1960 0.014 (0.012) 0.008 (0.017) 0.092 (0.039) 0.001 (0.023) 0.021 (0.041) 0.004 (0.012) 0.006 (0.016) 0.087 (0.036) 0.000 (0.021) 0.023 (0.038) 0 059 572 1953 0.007 (0.025) 0.001 (0.012) 0.083 (0.026) 0.109 (0.024) 0.003 (0.021) 0.006 (0.010) 0.068 (0.020) 0.102 (0.020) 0 106 880 1961 0.080 (0.015) 0.019 (0.011) 0.144 (0.024) 0.247 (0.030) 0.015 (0.038) 0.076 (0.015) 0.016 (0.011) 0.145 (0.023) 0.244 (0.029) 0.018 (0.037) 0 127 272 1959 0.004 (0.017) 0.010 (0.016) 0.089 (0.033) 0.022 (0.019) 0.029 (0.018) 0.006 (0.015) 0.007 (0.015) 0.080 (0.030) 0.021 (0.017) 0.026 (0.015) 0 169 880 1959 0.003 (0.013) 0.002 (0.017) 0.041 (0.050) 0.022 (0.059) 0.070 (0.097) 0.004 (0.012) 0.003 (0.015) 0.032 (0.044) 0.054 (0.052) 0.052 (0.079) 0 193 561 1939 0.033 (0.010) 0.014 (0.015) 0.250 (0.039) 0.110 (0.035) 0.078 (0.045) 0.014 (0.006) 0.002 (0.008) 0.104 (0.027) 0.051 (0.022) 0.060 (0.022) 0 240 808 1957 0.022 (0.019) 0.003 (0.014) 0.043 (0.029) 0.050 (0.031) 0.003 (0.018) 0.004 (0.012) 0.041 (0.024) 0.029 (0.026) 0 242 644 1960 0.009 (0.026) 0.021 (0.045) 0.181 (0.042) 0.203 (0.053) 0.133 (0.068) 0.011 (0.024) 0.019 (0.010) 0.171 (0.038) 0.231 (0.052) 0.128 (0.062) 0 247 668 1941 0.020 (0.012) 0.028 (0.012) 0.161 (0.033) 0.132 (0.034) 0.071 (0.062) 0.014 (0.007) 0.018 (0.008) 0.104 (0.020) 0.080 (0.022) 0.020 (0.038) 0 261 516 1962 0.025 (0.015) 0.038 (0.018) 0.179 (0.058) 0.079 (0.035) 0.264 (0.187) 0.025 (0.015) 0.038 (0.018) 0.179 (0.058) 0.079 (0.035) 0.264 (0.187) 0 290 468 1932 0.022 (0.012) 0.026 (0.016) 0.182 (0.047) 0.059 (0.028) 0.175 (0.087) 0.014 (0.006) 0.017 (0.011) 0.084 (0.024) 0.022 (0.015) 0 300 690 1957 0.001 (0.012) 0.019 (0.022) 0.126 (0.053) 0.061 (0.025) 0.018 (0.039) 0.00 (0.011) 0.014 (0.019) 0.116 (0.046) 0.055 (0.022) 0.018 (0.039) 0 317 475 1952 0.026 (0.009) 0.024 (0.015) 0.280 (0.035) 0.086 (0.025) 0.165 (0.070) 0.027 (0.006) 0.010 (0.011) 0.225 (0.027) 0.083 (0.021) 0.129 (0.047) 0 399 894 1961 0.042 (0.018) 0.036 (0.027) 0.182 (0.057) 0.096 (0.038) 0.195 (0.071) 0.035 (0.017) 0.028 (0.027) 0.165 (0.055) 0.084 (0.037) 0.173 (0.069) 0 403 886 1954 0.014 (0.018) 0.005 (0.030) 0.112 (0.074) 0.076 (0.042) 0.010 (0.096) 0.014 (0.017) 0.003 (0.029) 0.107 (0.070) 0.073 (0.040) 0.021 (0.092) 0 406 682 1954 0.003 (0.010) 0.043 (0.022) 0.190 (0.062) 0.064 (0.029) 0.119 (0.078) 0.005 (0.009) 0.039 (0.020) 0.159 (0.054) 0.051 (0.025) 0 432 237 1920 0.008 (0.011) 0.044 (0.019) 0.136 (0.059) 0.093 (0.031) 0.155 (0.115) 0.008 (0.006) 0.012 (0.007) 0.001 (0.029) 0.020 (0.015) 0.035 (0.055) 0 461 208 1959 0.017 (0.014) 0.039 (0.016) 0.340 (0.073) 0.130 (0.029) 0.674 (0.121) 0.016 (0.013) 0.033 (0.016) 0.332 (0.069) 0.134 (0.028) 0.651 (0.116) 0 513 284 1959 0.001 (0.015) 0.045 (0.026) 0.043 (0.055) 0.024 (0.037) 0.021 (0.053) 0.002 (0.014) 0.035 (0.023) 0.046 (0.048) 0.029 (0.032) 0.006 (0.049) 0 546 630 1936 0.002 (0.021) 0.033 (0.024) 0.041 (0.086) 0.008 (0.052) 0.028 (0.093) 0.006 (0.011) 0.014 (0.012) 0.080 (0.038) 0.037 (0.028) 0.034 (0.044) 0 554 816 1957 0.014 (0.012) 0.070 (0.024) 0.199 (0.058) 0.058 (0.038) 0.075 (0.048) 0.011 (0.011) 0.062 (0.021) 0.178 (0.049) 0.048 (0.033) 0.061 (0.041) 0 596 179 1951 0.003 (0.014) 0.067 (0.018) 0.052 (0.049) 0.065 (0.031) 0.010 (0.035) 0.004 (0.013) 0.053 (0.017) 0.072 (0.046) 0.068 (0.029) 0 677 802 1953 0.018 (0.011) 0.011 (0.019) 0.229 (0.049) 0.071 (0.031) 0.065 (0.051) 0.011 (0.009) 0.001 (0.015) 0.190 (0.040) 0.072 (0.024) 0.063 (0.038) 0 809 706 1934 0.020 (0.013) 0.019 (0.023) 0.198 (0.057) 0.111 (0.036) 0.147 (0.066) 0.007 (0.007) 0.001 (0.012) 0.074 (0.031) 0.044 (0.021) 0.042 (0.032) The values in brackets indicate the trends calculated for year P to 2009. indicates that the trend is statistically significant at the 5% level. indicates that the trend could not be calculated due to too many zeros in the time series.

676 A. C. KRUGER AND S. S. SEKELE Appendix B. Minimum temperature index trends for the period 1962 2009 Climate Number P Minimum Temperature Indices TNn (P) TNx (P) TN10P (P) TN90P (P) CSDI (P) 0 003 020 1911 0.029 (0.116) 0.010 (0.006) 0.125 (0.036) 0.169 (0.044) 0.032 (0.045) 0.028 (0.007) 0.011 (0.002) 0.244 (0.027) 0.129 (0.013) 0.393 (0.003) 0 012 221 1920 0.005 (0.010) 0.003 (0.008) 0.045 (0.051) 0.097 (0.048) 0.028 (0.109) 0.015 (0.004) 0.005 (0.003) 0.069 (0.018) 0.065 (0.016) 0.013 (0.021) 0 021 178 1957 0.040 (0.008) 0.012 (0.010) 0.218 (0.040) 0.103 (0.033) 0.040 (0.007) 0.011 (0.009) 0.230 (0.036) 0.108 (0.029) 0 035 209 1937 0.009 (0.016) 0.001 (0.009) 0.092 (0.046) 0.053 (0.036) 0.124 (0.200) 0.022 (0.008) 0.010 (0.004) 0.165 (0.029) 0.052 (0.018) 0.028 (0.038) 0 050 887 1960 0.024 (0.011) 0.002 (0.017) 0.093 (0.037) 0.040 (0.041) 0.019 (0.010) 0.003 (0.015) 0.081 (0.034) 0.039 (0.037) 0 059 572 1953 0.016 (0.008) 0.007 (0.009) 0.052 (0.033) 0.040 (0.032) 0.019 (0.010) 0.006 (0.007) 0.078 (0.029) 0.063 (0.025) 0.090 (0.124) 0 106 880 1961 0.023 (0.016) 0.061 (0.026) 0.158 (0.041) 0.066 (0.029) 0.026 (0.029) 0.023 (0.016) 0.057 (0.025) 0.157 (0.039) 0.067 (0.028) 0.021 (0.028) 0 127 272 1959 0.062 (0.012) 0.006 (0.014) 0.202 (0.037) 0.117 (0.004) 0.058 (0.011) 0.002 (0.010) 0.192 (0.034) 0.111 (0.035) 0.033 (0.021) 0 169 880 1959 0.034 (0.012) 0.005 (0.016) 0.051 (0.079) 0.198 (0.049) 0.030 (0.011) 0.013 (0.014) 0.010 (0.070) 0.172 (0.043) 0 193 561 1939 0.024 (0.016) 0.018 (0.015) 0.116 (0.051) 0.044 (0.036) 0.030 (0.036) 0.004 (0.009) 0.003 (0.009) 0.018 (0.027) 0.018 (0.021) 0.007 (0.019) 0 240 808 1957 0.040 (0.013) 0.004 (0.006) 0.267 (0.052) 0.063 (0.040) 0.039 (0.012) 0.006 (0.005) 0.281 (0.045) 0.079 (0.033) 0 242 644 1960 0.008 (0.011) 0.034 (0.04) 0.073 (0.074) 0.098 (0.062) 0.215 (0.118) 0.011 (0.011) 0.032 (0.022) 0.082 (0.068) 0.099 (0.057) 0.187 (0.108) 0 247 668 1941 0.020 (0.022) 0.030 (0.012) 0.147 (0.039) 0.114 (0.032) 0.095 (0.046) 0.015 (0.012) 0.017 (0.007) 0.148 (0.021) 0.099 (0.017) 0.120 (0.032) 0 261 516 1962 0.007 (0.015) 0.003 (0.016) 0.016 (0.050) 0.079 (0.044) 0.006 (0.056) 0.007 (0.015) 0.003 (0.016) 0.016 (0.050) 0.079 (0.044) 0.006 (0.056) 0 290 468 1932 0.056 (0.013) 0.015 (0.013) 0.137 (0.038) 0.110 (0.048) 0.020 (0.031) 0.019 (0.007) 0.005 (0.007) 0.045 (0.019) 0.001 (0.024) 0 300 690 1957 0.014 (0.014) 0.009(0.010) 0.029 (0.036) 0.037 (0.031) 0.002 (0.053) 0.012 (0.012) 0.010 (0.009) 0.047 (0.032) 0.041 (0.028) 0.002 (0.053) 0 317 475 1952 0.022 (0.015) 0.011 (0.013) 0.151 (0.052) 0.047 (0.041) 0.038 (0.052) 0.047 (0.012) 0.011 (0.009) 0.228 (0.043) 0.075 (0.028) 0.163 (0.053) 0 399 894 1961 0.013 (0.024) 0.004 (0.020) 0.052 (0.068) 0.003 (0.037) 0.109 (0.145) 0.013 (0.023) 0.002 (0.019) 0.055 (0.064) 0.007 (0.035) 0.116 (0.138) 0 403 886 1954 0.007 (0.017) 0.003 (0.012) 0.080 (0.050) 0.010 (0.042) 0.156 (0.049) 0.009 (0.016) 0.002 (0.011) 0.057 (0.049) 0.018 (0.043) 0.151 (0.046) 0 406 682 1954 0.062 (0.020) 0.023 (0.015) 0.114 (0.040) 0.095 (0.051) 0.059 (0.016) 0.027 (0.012) 0.120 (0.037) 0.098 (0.041) 0 432 237 1920 0.005 (0.018) 0.005 (0.015) 0.021 (0.038) 0.043 (0.048) 0.052 (0.039) 0.004 (0.007) 0.006 (0.007) 0.037 (0.019) 0.023 (0.017) 0.022 (0.020) 0 461 208 1959 0.037 (0.018) 0.027 (0.015) 0.154 (0.042) 0.114 (0.041) 0.069 (0.041) 0.040 (0.018) 0.029 (0.014) 0.145 (0.041) 0.102 (0.040) 0.059 (0.040) 0 513 284 1959 0.068 (0.015) 0.016 (0.010) 0.393 (0.050) 0.125 (0.030) 0.408 (0.123) 0.052 (0.014) 0.018 (0.009) 0.371 (0.045) 0.092 (0.031) 0.386 (0.110) 0 546 630 1936 0.014 (0.017) 0.001 (0.019) 0.086 (0.061) 0.012 (0.095) 0.069 (0.068) 0.005 (0.010) 0.004 (0.009) 0.076 (0.028) 0.008 (0.042) 0.024 (0.038) 0 554 816 1957 0.045 (0.021) 0.018 (0.008) 0.214 (0.044) 0.192 (0.049) 0.165 (0.062) 0.032 (0.018) 0.015 (0.007) 0.145 (0.042) 0.140 (0.044) 0.091 (0.055) 0 596 179 1951 0.021 (0.013) 0.007 (0.009) 0.131 (0.042 0.110 (0.038) 0.180 (0.073) 0.013 (0.013) 0.008 (0.008) 0.121 (0.039) 0.092 (0.035) 0.057 (0.052) 0 677 802 1953 0.017 (0.012) 0.007 (0.009) 0.106 (0.038) 0.078 (0.043) 0.037 (0.054) 0.013 (0.010) 0.006 (0.007) 0.140 (0.031) 0.088 (0.035) 0.054 (0.044) 0 809 706 1934 0.074 (0.017) 0.041 (0.013) 0.138 (0.044) 0.152 (0.063) 0.017 (0.043) 0.016 (0.011) 0.016 (0.007) 0.080 (0.024) 0.086 (0.032) 0.080 (0.026) The values in brackets indicate the trends calculated for year P to 2009. indicates that the trend is statistically significant at the 5% level. indicates that the trend could not be calculated due to too many zeros in the time series.