Changes in Thornthwaite Moisture Index and Reactive Soil Movements under Current and Future Climate Scenarios—A Case Study

: Expansive soils go through signiﬁcant volume changes due to seasonal moisture variations resulting in ground movements. The ground movement related problems are likely to worsen in the future due to climate change. It is important to understand and incorporate likely future changes in design to ensure the resilience of structures built on such soils. However, there has been a limited amount of work quantifying the effect of climate change on expansive soils movement and related behaviour of structures. The Thornthwaite Moisture Index (TMI) is one of the commonly used climate classiﬁers in quantifying the effect of atmospheric boundary on soil behaviour. Using the long-term weather data and predicted future changes under different emission scenarios, a series of TMI maps are developed for South Australia. Potential changes in ground movement are then estimated for a selected area using a simpliﬁed methodology where the effect of future climate is captured through changes in TMI. Results indicate that South Australia is likely to face a signiﬁcant reduction in TMI under all emission scenarios considered in this study. The changes in TMI will lead to a considerable increase in potential ground movement which will inﬂuence the behaviour of structures built on them and in some areas may lead to premature failure if not considered in the design.


Introduction
South Australia is the driest state on the driest inhabited continent in the world. The climate here has been going through quantifiable changes since the 1960s. Similar to other parts of the world, climate change here is likely to have a significant impact on the ecosystem, agriculture, and the built environment [1]. Climate change has already been attributed to reduced rainfall and increased temperatures in South-Eastern Australia [2][3][4]. Cai and Cowan [2] indicated climate change is responsible for a 20% reduction in rainfall during winter since the 1960s and with the possibility of increased greenhouse gases resulting in even greater rates of rainfall reduction in the future [5]. A decrease in cool seasonal rainfall across many regions of South Australia, more intense but shorter duration heavy rainfall across Australia, and increased temperatures are going to be the features of Australian climate [1, 6,7].
Climate change is already having an impact on the behaviour of geotechnical structures in many different countries [8][9][10][11][12][13][14][15][16]. Depending on the type of structure, local soil type, local climate and expected future changes, the impact can be different in different scenarios. Many contemporary construction or ground improvement methods use temperature and moisture-sensitive materials [17][18][19][20] and these changes need to be carefully considered to avoid unexpected losses in efficiencies of these techniques.
Climate change and its influence on soil behaviour have been an active area of research over the past two decades. However, the majority of the work has been focused on understanding the interactions involved in the soil-atmospheric boundary investigating soil-moisture dynamics and related changes in effective stress and stress-strain behaviour [16,[21][22][23][24][25][26][27][28]. One of the problems that is likely to worsen due to climate change and has been paid limited attention in the literature is reactive soil movement and its effects on shallow lightweight structures such as pavements, pipelines, or footings for residential buildings [29][30][31][32][33]. Approximately 30% of Australia's surface soil can be classified as reactive (also known as expansive soils). These soils go through significant volume change due to changes in moisture content resulting in ground movement, which can influence the safety and serviceability of the structures built on them.
Leao [34] investigated the effect of changes in ground movement due to climate change on the residential footing for the state of Victoria, Australia and assessed the associated risks. One of the conclusions from their study was that most urban growth areas were susceptible to medium to high risk for damage. However, the focus of the study was limited to changes in climate from the years 2001 to 2011. Mitchell [31] discussed the importance of considering climate change in the design of structures on expansive soils and used the Thornthwaite Moisture Index (TMI) as a climate classifier in his study. The scope of his study, however, was limited to four capital cities of Australia (Adelaide, Melbourne, Sydney, and Perth). Yue [35] investigated different methods of calculating TMI and compared three different methods [36][37][38] and indicated that when using a TMI-based approach in predicting climate-related unsaturated soil parameters, caution should be exercised as large variabilities among correlations/results from different methods were observed. Karunarathne [33] and Karunarathne [30] investigated the effect of the atmospheric boundary interaction, how it may affect ground movement, and how the estimation procedure may influence the footing design. Sun et al. [32] assessed the effect of climate change on expansive soil movement and site classification. However, their study mainly focused on the state of Victoria and the five major cities of Australia. Several past studies [39][40][41] have also investigated the influence of atmospheric boundary interaction and climate change on the behaviour of lightweight structures such as road pavements on expansive soils. However, the findings may not directly translate to the Australian climate, more specifically the South Australian climate, as the studies were conducted for a different part of the world. There has been very limited research on the effect of climate change on expansive soils in South Australia. To the author's best knowledge, no predictive TMI map exists for South Australia based on the most recent climate predictors that designers can use as a reference tool.
AS2870-2011 [42], the Australian standard relating to the design of footings for lightweight structures on expansive soils, provides a guideline for the estimation of the depth of design soil suction change (H s ) based on the TMI. H s is the depth below which the short-term climatic variation on a built-up greenfield site does not affect the soil moisture content. Therefore, the soil below H s does not contribute to ground movement. With knowledge of the soil profile and the rate of vertical strain for a unit change in suction for the various soils throughout the soil profile, the design soil suction changes for the site are used to estimate the potential differential ground movement, or design characteristic ground movement (y s ) [43]. This parameter influences the designed stiffness of the near-surface footing, which is required to maintain the serviceability of the supporting structure. Since TMI affects the depth of design suction changes, it indirectly affects footing design for footings built in/on expansive soils. TMI may have similar influences on other geotechnical structures. Therefore, prediction of the potential TMI changes in the coming future is necessary to ensure that infrastructure built today can last into the future [44,45]. Given that changes to temperature and rainfall caused by climate change will lead to changes in TMI, geotechnical designers must be adequately aware and prepare for such effects [46].
In this paper, the literature for three commonly used methods for calculating TMI is briefly discussed. Sixty-six years of weather data from a weather station in South Australia were used to calculate TMI using the methods discussed and compared. Based on the comparison, a selected method is used to calculate a set of "current TMI" values for 70 locations across South Australia and plotted in a GIS map. The "current TMI" calculation is based on the most recent 30 years of weather data from weather stations at the selected locations. Based on the forecasted changes in climate for the years 2050 and 2080 in different emission scenarios, a series of predictive TMI maps are also developed. A simplified methodology is developed to estimate the future changes in ground movement and implemented for a selected area of South Australia. A series of ground movement maps are produced representing current and different future climate scenarios. The consequence of these changes is also discussed using a simplified footing design example.

Methods for Estimating TMI
TMI, calculated as a function of rainfall and potential evapotranspiration, was first introduced in 1931 by Thornthwaite [47] and has been widely accepted and employed within the scientific community in many applications. Thornthwaite [36] initially applied TMI as a tool to map the climate zones of the United States of America, which was further extended in applications in the agricultural field as an indicator of suitable climates for various crop types [48]. TMI has also been used in other areas, including ecological community, prediction of water catchment capacities and forest growth, and estimation of potential soil suction beneath covered areas and pavements [49][50][51]. In Australia, TMI has been adopted to understand general differences in climate across the country [49,52] and estimate the depth of design soil suction changes H s, which is needed to determine y s [42].
Over several decades, various methods/equations have been developed to estimate TMI. The original method, developed by Thornthwaite [36], calculated TMI as, where, I h = humidity index = 100 R PET ; I a = aridity index = 100 D PET ; R = moisture runoff or quantity of moisture that cannot be absorbed in an already wet site; D = moisture deficit; PET = potential evapotranspiration or the possible loss of moisture from the ground surface and vegetation. This method will be referred to as Method 1 from this point onward.
According to Method 1, I h and I a quantitatively captures the soil's ability to store water between the extreme cases of completely dry and fully saturated states and could be used to distinguish extremes in seasonal moisture conditions, such as droughts and periods of prolonged rain. Method 1 also attempted to account for the effects of vegetation on the moisture balance. The 0.6 multiplier to the aridity index was introduced into the method to account for an assumption that water can enter the soil at a faster rate than it can be extracted from the soil profile [5,36,53]. This method requires estimation of (i) monthly moisture surplus-R, (ii) monthly moisture deficit-D, and (iii) monthly adjusted potential evapotranspiration-PET. The first two parameters (i.e., R and D) can be derived using a water balance approach and the calculation procedure has been described by McKeen and Johnson [54]. Empirical methods for calculating PET have evolved over the past decades; however, the common approach in the literature is by using the adjusted potential evapotranspiration for the month, i (PET i ) and can be determined by using the Thornthwaite [36] evapotranspiration equation as presented in Equation (2) below. The equation simply relies on precipitation and temperature data which are readily available from most of the weather stations. here, D i is the day length correction factor for the month i; N i is the number of days in the month i; pe i is the non-adjusted potential evapotranspiration for the month i and can be calculated as: where, (t i ) is the mean monthly temperature in • C of the study location for a given month (i) and can be calculated as the average of the minimum and maximum temperature. H y is the yearly heat index and can be determined using t i as: where, y denotes the given year. The power term a in Equation (3) is a coefficient that ranges between 0 < a < 4.25, and can be calculated as: It is important to note that a positive value of TMI indicates a humid or wet climate with higher soil moisture, while a negative value of TMI represents an arid climate with less moisture in the soil. Zero TMI indicates there is no net soil moisture change for the study period. In other words, the soil moisture is in balance with environmental conditions. Thornthwaite and Mather [37] modified the TMI equation by removing the coefficient of 0.6 for the aridity index, as shown in Equation 6, and this made the calculation process considerably simpler.
here, P is the precipitation. Equation (6) will be referred to as Method 2 in this paper. Witzec et al. [38] further modified Method 2 as below, This method will be referred to as Method 3 in this paper.

Comparison of Methods 1, 2, and 3
A comparative study was conducted using the three methods discussed in the previous section. A total of 66 years of weather data (1955 to 2020) were collected from a weather station in Adelaide (Adelaide Airport, Lat. −34.95 degree, Long. 138.52 degree). Collected weather data included monthly rainfall, mean maximum, and mean minimum temperature. TMIs were calculated using the collected weather data and Methods 1, 2, and 3. The calculated TMIs from the different methods are presented in Figure 1 below. Significant differences could be observed between Method 1 and Method 2. Calculations from Method 1 and Method 3, however, were very close to each other. Provided that AS2870-2011 [42] uses Method 1 for TMI estimation and Method 3 is significantly simpler to apply yet produces results very close to Method 1, Method 3 was used for all further calculations of TMI in this study.

Future Climate of South Australia
Climate Change Australia [6] presents a series of future climate scenarios for different regions (clusters) of Australia. The state of South Australia falls into three of those climate clusters, namely Rangelands, Southern and South-Western Flatlands, and Murray Basin. The climate clusters are shown in Figure 2. Tables 1 and 2 summarize the projected changes in temperature and rainfall respectively for different emission scenarios for 2050 and 2080. Three different emission scenarios, i.e., representative concentration pathways RCP2.6, RCP4.5, and RCP8.5 are summarized here. RCP2.6 assumes the greenhouse gas emission peaking around the year 2020 and then declining rapidly. The CO2 concentration in the atmosphere is expected to reach 440 ppm by 2040 and then reduce to 420 by the end of the 21st century. RCP4.5 assumes the emission to peak around 2040 and the CO2 concentration to reach 540 ppm by the year 2100. RCP8.5 is the most pessimistic scenario where minor changes in the current emission trend are assumed and CO2 rises to 660 ppm by the end of this century [55].
Under different emission scenarios, the temperature of different regions of South Australia has been predicted to rise by 0.8 to 1.9 °C by 2050 and by 0.95 to 4.2 °C by 2080. In most areas of South Australia, the rainfall is expected to decrease in all seasons of the year, with the maximum reduction expected in winter and spring. It is to be noted that Climate Change Australia [6] provides the probabilistic changes in rainfall and temperature pattern; the numbers presented here are deduced average values.

Future Climate of South Australia
Climate Change Australia [6] presents a series of future climate scenarios for different regions (clusters) of Australia. The state of South Australia falls into three of those climate clusters, namely Rangelands, Southern and South-Western Flatlands, and Murray Basin. The climate clusters are shown in Figure 2. Tables 1 and 2 summarize the projected changes in temperature and rainfall respectively for different emission scenarios for 2050 and 2080. Three different emission scenarios, i.e., representative concentration pathways RCP2.6, RCP4.5, and RCP8.5 are summarized here. RCP2.6 assumes the greenhouse gas emission peaking around the year 2020 and then declining rapidly. The CO 2 concentration in the atmosphere is expected to reach 440 ppm by 2040 and then reduce to 420 by the end of the 21st century. RCP4.5 assumes the emission to peak around 2040 and the CO 2 concentration to reach 540 ppm by the year 2100. RCP8.5 is the most pessimistic scenario where minor changes in the current emission trend are assumed and CO 2 rises to 660 ppm by the end of this century [55].

TMI Map of South Australia for Current and Future Climate Scenarios
TMI is an annual index and is calculated from monthly historical data. It can be calculated for a particular location for any specific year, as shown in Section 3 of this paper, and this approach is generally known as the year-by-year method. An alternative to the year-by-year method is to assemble monthly average temperatures and rainfalls using all existing data and calculate TMI for a single average year. This method is known as average year analysis. For simplifying the calculation, an average year approach was used in this paper. The calculation steps are outlined below,  Year 2050 Year 2080 Under different emission scenarios, the temperature of different regions of South Australia has been predicted to rise by 0.8 to 1.9 • C by 2050 and by 0.95 to 4.2 • C by 2080. In most areas of South Australia, the rainfall is expected to decrease in all seasons of the year, with the maximum reduction expected in winter and spring. It is to be noted that Climate Change Australia [6] provides the probabilistic changes in rainfall and temperature pattern; the numbers presented here are deduced average values.

TMI Map of South Australia for Current and Future Climate Scenarios
TMI is an annual index and is calculated from monthly historical data. It can be calculated for a particular location for any specific year, as shown in Section 3 of this paper, and this approach is generally known as the year-by-year method. An alternative to the year-by-year method is to assemble monthly average temperatures and rainfalls using all existing data and calculate TMI for a single average year. This method is known as average year analysis. For simplifying the calculation, an average year approach was used in this paper. The calculation steps are outlined below, (1). The mean temperature for each month of the year was determined by averaging the maximum and minimum monthly temperatures over 30 years for that particular month. (2). The heat index for each month was then calculated using Equations (4) and (5). (3). The next step involved the calculation of unadjusted potential evapotranspiration using Equation (3). (4). Calculated unadjusted values of potential evapotranspiration were then adjusted for day length using Equation (2). The day length correction factors for relevant latitudes are presented in Table 3. (5). TMI was calculated using Equation (7). The past 30 years of weather data (1991 to 2020) from 70 weather stations scattered across South Australia were used for this exercise, and the calculated set of TMI values were tagged as "Current TMI". Using GIS software, the "Current TMI" data were plotted in a contour map with contour intervals of 10. The map is presented in Figure 3a with the city of Adelaide and surrounding areas highlighted. The past 30 years of weather data (1991 to 2020) from 70 weather stations scattered across South Australia were used for this exercise, and the calculated set of TMI values were tagged as "Current TMI". Using GIS software, the "Current TMI" data were plotted in a contour map with contour intervals of 10. The map is presented in Figure 3a with the city of Adelaide and surrounding areas highlighted.
Four other sets of TMI values were calculated for 2050 low, 2050 high, 2080 low, and 2080 high emission scenarios. The TMI map for the 2050 low emission scenario was developed based on an average temperature increase of 0.95 °C and a 4% reduction in rainfall (yearly average across the three clusters for RCP 2.6) applied to the past 30 years of average monthly temperature and rainfall data, respectively. That is, an average monthly temperature increase of 0.95 °C was added and a 4% rainfall was deducted from the past 30 years monthly average data.
Similarly, the 2050 high emission scenario was based on a 1.72 °C rise in temperature accompanied by a 5.9% reduction in rainfall. The 2080 high and low scenario TMIs were calculated based on 0.98 and 3.75 °C rises in temperature and 4.75 and 8.7% reduction in rainfall, respectively, related to RCP2.6 and RCP8.5 climate change scenarios. The predictive TMI maps for 2050 and 2080 are presented in Figure 3b-e. Calculated TMI values for all weather station locations are presented in Table A1 in Appendix A at the end of this paper.   Table 4 below. One interesting observation could be made from Table 4 and Appendix A that the areas with relatively wetter climates (higher TMI) are affected more significantly (in future climate scenarios) compared to areas with lower TMI under the current climate.
It is difficult to quantify the changes with a high degree of certainty as climate change is a function of a range of different factors. However, the developed maps are indicative of future changes and engineers working on design projects should accommodate these changes in their current design if the structure is expected to last over its design life, for example, 50 or 100 years.

Changes in Ground Movement in a Local Government area of South Australia
The characteristic expansive soil movement was calculated using a simplified method. According to AS 2870-2011 [42], for Adelaide and surrounding areas, the changes in suction at the soil surface between driest and wettest points on an allotment over the structure's design life (Δu) can be taken as 1.2 pF, and the depth of cracking (Hc) can be taken as 3 m. A triangular distribution of soil suction changes then can be assumed, which extends to the depth Hs as shown in Figure 4 below. AS2870-2011 [42] provides guidelines that can be used for estimating Hs based on TMI of a particular area. The characteristic ground movement, ys can be calculated as, Four other sets of TMI values were calculated for 2050 low, 2050 high, 2080 low, and 2080 high emission scenarios. The TMI map for the 2050 low emission scenario was developed based on an average temperature increase of 0.95 • C and a 4% reduction in rainfall (yearly average across the three clusters for RCP 2.6) applied to the past 30 years of average monthly temperature and rainfall data, respectively. That is, an average monthly temperature increase of 0.95 • C was added and a 4% rainfall was deducted from the past 30 years monthly average data.
Similarly, the 2050 high emission scenario was based on a 1.72 • C rise in temperature accompanied by a 5.9% reduction in rainfall. The 2080 high and low scenario TMIs were calculated based on 0.98 and 3.75 • C rises in temperature and 4.75 and 8.7% reduction in rainfall, respectively, related to RCP2.6 and RCP8.5 climate change scenarios. The predictive TMI maps for 2050 and 2080 are presented in Figure 3b-e. Calculated TMI values for all weather station locations are presented in Table 1 in Appendix A at the end of this paper.
The comparison exercise (see Figure 3a-e) indicates that the South Australian climate is likely to get drier under all four future climate scenarios considered here. An average (average values from all 70 weather stations) TMI change of −3.4 was observed for the 2050 low emission scenario, which increased to −5.5 in the high emission scenario for the same year. For 2080, the average changes were −3.6 and −10.4 for high and low emission scenarios, respectively. The average TMI values for South Australia along with average changes and minimum and maximum values are summarized in Table 4 below. One interesting observation could be made from Table 4 and Appendix A that the areas with relatively wetter climates (higher TMI) are affected more significantly (in future climate scenarios) compared to areas with lower TMI under the current climate.
It is difficult to quantify the changes with a high degree of certainty as climate change is a function of a range of different factors. However, the developed maps are indicative of future changes and engineers working on design projects should accommodate these changes in their current design if the structure is expected to last over its design life, for example, 50 or 100 years.

Changes in Ground Movement in a Local Government area of South Australia
The characteristic expansive soil movement was calculated using a simplified method. According to AS 2870-2011 [42], for Adelaide and surrounding areas, the changes in suction at the soil surface between driest and wettest points on an allotment over the structure's design life (∆u) can be taken as 1.2 pF, and the depth of cracking (H c ) can be taken as 3 m. A triangular distribution of soil suction changes then can be assumed, which extends to the depth H s as shown in Figure 4 below. AS2870-2011 [42] provides guidelines that can be used for estimating H s based on TMI of a particular area. The characteristic ground movement, y s can be calculated as, y s = ∆z∆u αI ps (8) were, α = 2 − z/5 above H c , and α = 1 below H c ; z is the depth, I ps is the soil reactivity. α takes into account the degree of lateral confinement and to some extent overburden pressure.
Energies 2021, 14, x FOR PEER REVIEW 10 of 18 were, = 2 − 5 ⁄ above , and = 1 below ; z is the depth, is the soil reactivity. α takes into account the degree of lateral confinement and to some extent overburden pressure. For this exercise, a town centre and surrounding area in South Australia (not named here) was selected based on data availability. A soil reactivity database of 135 locations within the selected area was created based on a published soil movement map [56]. Hs was estimated for each of the 135 data points based on the corresponding location's TMI for the "current" climate scenario using a correlation established from AS2870 [42]. The same exercise was repeated for 2050 low and 2080 high emission future climate scenarios. Based on the estimated Hs and using Equation (8), ys was calculated for each of the 135 locations for each of the different climate scenarios.
The soil movement map developed based on Hs estimated from current TMI values is presented in Figure 5a. It is to be noted that only the effect of changes in climate on Hs was considered in this exercise and the changes in Δu was ignored to avoid complexity and lack of reliable methods to calculate them. Increased temperature and reduced rainfall are likely to increase Δu for many soil types and thus ignoring this may have contributed to some underestimation of ys in the current exercise.  For this exercise, a town centre and surrounding area in South Australia (not named here) was selected based on data availability. A soil reactivity database of 135 locations within the selected area was created based on a published soil movement map [56]. H s was estimated for each of the 135 data points based on the corresponding location's TMI for the "current" climate scenario using a correlation established from AS2870 [42]. The same exercise was repeated for 2050 low and 2080 high emission future climate scenarios. Based on the estimated H s and using Equation (8), y s was calculated for each of the 135 locations for each of the different climate scenarios.
The soil movement map developed based on H s estimated from current TMI values is presented in Figure 5a. It is to be noted that only the effect of changes in climate on H s was considered in this exercise and the changes in ∆u was ignored to avoid complexity and lack of reliable methods to calculate them. Increased temperature and reduced rainfall are likely to increase ∆u for many soil types and thus ignoring this may have contributed to some underestimation of y s in the current exercise. The soil movement map developed based on Hs estimated from current TMI values is presented in Figure 5a. It is to be noted that only the effect of changes in climate on Hs was considered in this exercise and the changes in Δu was ignored to avoid complexity and lack of reliable methods to calculate them. Increased temperature and reduced rainfall are likely to increase Δu for many soil types and thus ignoring this may have contributed to some underestimation of ys in the current exercise.  Figure 5a shows that majority of the study area under the current climate has ys less than 40 mm. According to AS2870 [42], these sites can be classified as slightly to moderately reactive clay or silt sites with deep-seated moisture changes. A small proportion of the area has ys > 40 mm with an even smaller proportion of the area with ys > 50 mm.  Figure 5a shows that majority of the study area under the current climate has y s less than 40 mm. According to AS2870 [42], these sites can be classified as slightly to moderately reactive clay or silt sites with deep-seated moisture changes. A small proportion of the area has y s > 40 mm with an even smaller proportion of the area with y s > 50 mm. AS2870 [42] considers these sites to be experiencing high to very high ground movement.
In Figure 5b,c, the y s for the same area under 2050 low emission and 2080 high emission scenarios are presented. y s values have been calculated using a similar approach as Figure 5a but with H s calculated based on changed TMI. The average H s increased by 0.17 m and 0.51 m respectively for 2050 low and 2080 high emission scenarios resulting in significant changes in y s in the two forecast movement maps. The area with y s < 40 mm reduced significantly with an accompanying increase in areas with larger movement potential. The area with y s > 50 mm which was limited to less than approximately 1% of the investigated area under the current climate scenario, increased by approximately 10-fold in the 2080 high emission scenario.
On average the y s values increased by 4.6% and 13.4% with maximum changes of 4.23 and 12.2 mm for 2050 low and 2080 high emission scenarios respectively. Increased ground movement potential will mean structures built on such sites will have to cope with additional stresses exerted to them by the soil due to seasonal variations in moisture in soils and may fail if the expected future changes are not incorporated in their design.
To highlight the possible implication of the changes in ground movement, a hypothetical case is discussed here. A rectangular stiffened raft footing (15 m × 10 m) with three and four beams in the long and short directions, respectively) is designed for a lightweight residential building on an expansive soil site characterized by y s of 50 mm under the current climate. The design exercise is repeated for a possible future scenario where a 10 mm increase in y s is expected. It is to be noted that in the design of shallow footings for lightweight structures on reactive soil sites, the conventional method based on bearing capacity principles do not apply. The design is dictated by stresses applied to footing structures by expansive soil movement and works on the principle of limiting the differential movement of the footing structure within the tolerance of the superstructure. The design inputs are summarized in Table 5 below. A computer program CORD [57] which incorporates the Walsh method [58,59] for estimating the mound shape was used for the design.
A stiffened slab designed under the current climate scenario would require 300 mm beam depth for an articulated masonry veneer type superstructure. For a 10 mm increase in ground movement, the design for the same superstructure will require a beam depth of 380 mm and a slightly higher reinforcement area.
The analysis indicates, for geotechnical structures where soil-atmospheric boundary interaction significantly affects the stability, the changes in climate need to be taken into account. Using the methodology outlined here, tools can be developed to estimate expected future soil movements and incorporated into the design rather than using soil's current movement potential. The maps presented here can also assist in the establishment of policies and strategies towards the development of resilient infrastructure.
It should be noted that the extrapolation used in estimating the ground movement may not be very accurate as they are based on gross approximation and can be argued upon and certainly requires further research and refinement. To estimate the effect of future climate, the changes in temperature and rainfall have simply been added to the past 30 years' average values. This can be looked at as an oversimplification of the interaction that occurs at the soil-vegetation-atmosphere boundary. However, the change in suction depth is a seasonal process and is contributed to by incremental changes due to the soil-atmospheric boundary interactions over months (often 6 months, i.e., the difference in time between the end of summer to the end of winter). The current methodology cannot capture the effect of individual weather events explicitly especially considering more intense but shorter duration heavy rainfall will be a feature of the future climate of Australia. However, considering the suction depth changes occurs over a long time, the assumptions can be treated as reasonable. Nonetheless, it highlights the importance of incorporating/accommodating climate change-related effects on the designs of geotechnical structures today, especially those located at shallow depths such as pipelines, pavements, and so forth.

Conclusions
The followings are the conclusions from this study, • An increase in temperature and reduction in rainfall will be the feature of the future climate in South Australia with a 0.8 to 4.2 • C rise in temperature and up to 16

Conflicts of Interest:
The authors declare no conflict of interest.

Disclaimer:
The opinion expressed in this paper are of the authors and the assumptions must be understood and engineering judgement should be used for its use in any design or other purpose.