The Influence of Socioeconomic Factors on Households’ Vulnerability to Climate Change in Semiarid Towns of Mopani, South Africa

The changing climate and its current rate, frequency, as well as its life-threatening impacts are undoubtedly abnormal and globally worrisome. Its effects are expected to be severely different across segments of the society. It is disposed to leaving no facet of human endeavor immune, particularly in vulnerable cities of developing countries where there is dearth of empirical studies. For the context-specific nature of climate change impacts and place-based character of vulnerability, this study explores the influence of socioeconomic attributes on household vulnerability in Mopani District northeast of South Africa to provide basis for targeting, formulating, evaluating, and monitoring adaptation policies, programs, and projects. The study adopted a multistage random sampling to draw 500 households from six towns in Mopani District, Limpopo Province. Mixed methods approach was used for data collection, while Household Vulnerability Index (HVI) was estimated using principal component analysis and regressed with socioeconomic attributes. The study reveals that climate is changing with high HVI across selected towns. It further depicted that age and marital status have positive and significant relationships with HVI, while gender and educational levels have inverse and significant relationship with HVI in some towns. The study recommends the need for municipalities to partner with private sector to empower household and mainstream micro level coping strategies in urban planning across the district.


Introduction
Climate is changing and the adduced factors to its variability include both natural but substantially human induced [1,2]. The changing climate comes with attendant negative global impacts such as increasing temperature, sea level rise, floods, etc. [3]. These consequences are manifesting in its complex interdictions, uncertainties, and competing interest among the stakeholders as well as the multigenerational nature of its challenges. It is multigenerational because it covers far beyond the lives of people currently attending to its externalities [4]. The current rate at which the climate is changing, particularly, the post 1990, is unprecedented and worrying. The increasing temperature signal is very strong for southern African region with the observed trend indicating an increase of more than 150% than the global rate [5] The "projected warming over 1990-2100 ranges between 1.4 and 2.9 • C for the central emissions scenario (1595a), while sea-level rise ranges between 20 and 86 cm . . . ", particularly for cities of developing countries [1,6] and their vulnerable population, particularly young, old, and women.
Developing countries are highly vulnerable in the region that bears the significant share of the global climate change consequences, despite their insignificant carbon footprint [7]. In the Forth Assessment Report of Intergovernmental Panel on Climate Change's This prompted three fundamental research questions as follows: Is climate changing in Mopani District? Can it be said that households are vulnerable to climate change in Mopani District in South Africa's semiarid (northeast) region? And can it be said that socioeconomic characteristics have influence on the Households' Vulnerability Index (HVI) across the towns and the District? No doubt, the frequency, magnitude, and consequences of climate change differ the world over, in accordance with their peculiarities and the virtue of the people's unique individualities (socioeconomic contents) [34]. It is on this basis that this study aimed at examining the variability and trends of climatic change (rainfall and temperature) in the Mopani District and to assess the HVI as well as analyze the nexus between the HVI and socioeconomic attributes of households in the six selected towns (Tzaneen, Nkowankowa, Hoedspruit, Modjadjiskloof, Phalaborwa, and Giyani). This was intended to provide basis for targeting, formulating, evaluating, and monitoring adaptation policies, programs, and projects.

Description of the Study Area
Mopani District Municipality is a category C municipality. It is located in Limpopo Province, the northern-most province in South Africa ( Figure 1). The district comprises of five local municipalities, namely, Greater Giyani, (district's administrative seat), Maruleng, Greater Letaba, Ba-Phalaborwa, and Greater Tzaneen. The municipality is on longitudes: 29 • 52 E to 31 •   Most part of Mopani District receives about 85% of its rainfall in the summer. The rainfall amount varies from the mountainous zones with about 2000 mm/annum to as low as 400 mm/annum in Kruger National Park, while the district's high mountainous zone experience on the average temperature range of 21-25 • C in the dry low veld areas of Kruger National Park [35].
The population of Mopani District is increasingly becoming susceptible to the increasing incidences of extreme climate events of floods, drought, heat waves, and the spread of climate-sensitive diseases [28][29][30]. These are occasioned by high rates of unemployment and poverty, limited access to resources, and high infrastructure deficit [36,37]. This justifies the need for empirical examination of household vulnerability level and its nexus with the socioeconomic contents across the selected towns in each of the local municipalities.

Methods
Several methods are available in the literature for analyzing vulnerability. The "starting point" known as disaster focused, an approach that assesses vulnerability to hazards such as flooding and heat [38], and "end point," i.e., those approaches that see and assess vulnerability as an aftermath effects of the hazards, are the two most distinct [39,40]. Adopting an isolated method comes with its attendant reservations and shortfalls [40]. Hence, this study integrates the two approaches to eradicate the shortfalls that characterize the usage of an isolated approach because of the indispensability of the two approaches [41,42].
Households' exposure sequels to the establishment of the trend of climate parameters in the study area, whereas households' duration of stay in their various places of abode signifies how long such households' have been exposed to the effects of the trending climate. These values were routed through Likert Scale of Exposure Time spent: not exposed (1) for staying for less than a year; just exposed (2) for 1-5 years of exposure; exposed (3) for >5-10 years of exposure; very exposed (4) for 11-15 years of exposure, and very much exposed (5) for >15 years of exposure.
Households' sensitivity described by the degree to which a system is affected or modified by either internally or externally induced disturbances [43]. Sensitivity was measured from data regarding households' loss of properties, fatalities, i.e., death of a family member and livestock, income that is natural resource-based, or remunerative income ( Table 1). The estimation of both mean and standard deviation was conducted. The Sensitivity Index were further categorized in to three types (not sensitive, sensitive, and very sensitive).
Households' adaptive capacity is taken to be a constituent property of the several indices of livelihood assets. Adopting integrated approach [1,44] household's adaptive capacity was measured using Personal Possession Index (PPI). It is, thus, submitted that for the 0-1 scale, the calculated indicator value is more asymptotic to 1 indicates a higher adaptive capacity, but a calculated indicator asymptotic to zero (0) indicates a weak or low adaptive capacity [1]. However, these indices were further subjected to the Likert Scale of Capability as: very much capable (5) for possession of 81-100% of prescribed assets; very capable (4) for possession of 61-80% of prescribed assets; capable (3) for possession of 41-60% of prescribed assets; not capable (2) for possession of 21-40% of prescribed assets, and not capable at all (1) for possession of less than 21% of prescribed assets.
Estimating vulnerability using indicators-based approach, this study adopted the UN Habitat (2011), which indicates vulnerability as a measure of adaptive capacity, exposure, and sensitivity to climate-prone hazards. This is symbolized as: The Analytical Framework was to establish Household Vulnerability Index. For the population, N made up of n households, i.e., N = hh1, hh2, hh3, . . . , hhn, V is a subset of v households that have some degree of vulnerability, hence having an internal vulnerability.
Thus v ≤ n and v = 0 implies that there are no vulnerable household, and v = n implies that all households are vulnerable.
The first step witnessed the identification of key indicators ( Table 1) that contribute to the vulnerability of a target system. To break down the vulnerability X into m specific dimensions of impact and give a corresponding weight (wi, i = 1, . . . , m) to each indicator, weights was assigned to these indicators rather than assigning equal weights across all indicators or using expert-based judgement [39] because of the criticism for its subjectivity and disagreements among experts [41]. This study adopted principal component analysis (PCA) following [44], as considered appropriate in a work of this nature [45]. Thus, PCA was separately ran for the selected indicators of exposure, sensitivity, and adaptive capacity to determine the magnitude of influence of each variable on the index. The loadings from the first component of the PCA was used as the weights for the indicators.
The study used both quantitative and qualitative data as guided by the hypothesized indicators (Table 1). The quantitatively measured data were, for instance, ages and duration of stay in a locality by the respondents, while housing was measured qualitatively in terms of housing condition. Making comparability of the results possible, we normalized the quantitative data normalized to 0-1 scale following [1], adopting the following equation where X ij represents the standardized indicator value i of the household j; X ij is the indicator value i that corresponds to household j; Max and Min indicate the maximum and minimum scaled indicators' values i, respectively.
In case the standardized value is more asymptotic to 1, it will imply a higher vulnerability, and it is of low vulnerability if the value is more asymptotic to 0.
To provide broad information on the urban Households' Vulnerability Index, the mean and standard deviation of each components (exposure, sensitivity, and adaptive capacity) as well as the vulnerability were computed.
The weighting of the indicators ranges from 0 and 1, signs indicating the path in relation to other indicators for the respective indexes, with dimension values 0 = no impact and 1 = full impact. The magnitude of calculated value of the weights defines the contribution of each indicator to the value of the index.
The sum of the weighted vulnerabilities across all dimensions give the particular household's total vulnerability Vhhi: The summation of the dimensions gives the contributions of each of the identified dimensions to vulnerability.
For the HVI, the sum of the weights is set to To establish the extents of household vulnerability, several authors have adopted the Multinomial Logistic Regression Model (MLRM) but because of its limiting assumptions, that there is no order to the categories of the outcome variable (i.e., the categories are nominal), the information contained in the ordering is lost [46]. The Ordered Logit Regression Model (OLRM) on STRATA 14, an extension of the binary logic model was adopted. The choice of OLRM is informed by its ability to calculate the individual Household Vulnerability Index, the total households, as well as ordered the contributions of each hypothesized indicator [47]. In this study, the ordered dependent variable has three categories: low vulnerability, moderate vulnerability, and high vulnerability. The low vulnerability category was used as the reference category. This estimation will show, for instance, lowly vulnerable households HVI in range from 0 to 47, moderately vulnerable households HVI in range from 47.1 to 63.7, and highly vulnerable households HVI in range from 63.71 to 100.
The study used high-resolution data set from global climatic data set (temperature and rainfall) of monthly climate data based on a grid from latitudes 24.4 • S to 23.3 • S and 30.0 • E and 31.1 • E, during the period of 1958-2017 for temperature and 1958-2016 for rainfall. However, several studies like [48,49] justify this step. Table 1 summarizes the hypothesized indicators that provided guides for the evaluation of HVI.

Climate Change in Mopani District (1958-2017)
The temperature and rainfall data were analyzed by estimating the mean interannual cycle for the trend of rainfall as well as the trend of anomalies of temperature to indicate extreme cases over the District during the period of 1958-2017.
The monthly average result of analysis as contained in Figure 2 illustrates temperature trend between 1958 and 2017. From these figures, the temperature exhibits an increasing trend, rising steadily such that from about the year 2000, they have been mostly above the long-term mean. The implication of this trend is that Mopani District displayed warm temperature over the period of examination. However, despite the long-term trend, there are still seasonal differences from year to year with hotter temperature during drought seasons. It is equally clear that interannual unevenness across seasons in the district as shown in Figure 3 signifies high rainfall variability in the district. These have consequently occasioned a shift in rainfall to a later onset (extending the dry fire season), whereas extreme floods, often due to cyclones and droughts, connected to the El Nino phenomenon, have turned out to be more frequent. The recent incidences of extreme events have continuously put Mopani towns in an untold hardship, particularly, the most vulnerable groups (children, older adults, and women).
This study adopted the temperature projections from existing studies [18,25,28]. The studies have very strong increasing temperature signals for southern African region with more than 150% higher than the global rate. Precipitation is projected to assume about 5%-10% reduction over South Africa in the next five decades [6], with greater uncertainties than those of temperature. The situation calls for pragmatic policies and serious attitudinal change, to forestall devastating human consequences as may likely be escalated by socioeconomic attributes of the people. Table 2 summarizes in no specific order, the households' socioeconomic attributes considered in this study. It is shown that the larger proportion of the household's heads are still in their prime and active ages between the range of 20 and 50 years. The category laid claim to over 80% of the total samples across the district. The number of household heads decrease with increased age, implying that the younger household heads are more than the older ones. The implication of this scenario suggests that any program that is targeted towards economic empowerment of the people that have special focus on gender and young people will significantly be impactful across the selected towns. This will equally be true on programs focusing on climate change adaptation and early warning propagation programs in the district by the local municipalities, targeting this population category. These groups are not only young but also are economically and socially active that, if facilitated, they can be productively engaged to initiate and champion climate change risk sensitization, mitigation, and adaptation initiatives. These findings validate the findings of [66,67] In their respective studies of Livelihood and income diversification strategies among rural farm households in Niger State, (Northern) Nigeria, and that of poverty and rural livelihood diversification among farming households in southwest Nigeria.

Households' Socioeconomic Characteristics
The results further show that about 57% of household heads are married, while about 43% are either divorced, single, widow, or separated. This may suggest additional pressure from extra responsibility or imply an extra burden to adapt during and after climate change extreme events. Table 2, also shows that about one-third of the sampled population is either job seekers or surviving on income below or equivalent to R16.44 per day (R6000, i.e., $312 per annum). Quiet a significant proportion (56%) engaged in activities that earn them between R240,000 ($12,400) and R480,000 ($24,960) per annum, while about 1% earn above R480,000 ($24,960) per annum. Majority of the household heads having financial resources at their disposals (as indicated) signifies a higher household affordability. This is capable of enhancing the households' coping capacity during and after climate change extreme events. These findings validate those of [68,69] in their analysis of the determinants of rural poverty among small holder farmers in South-western, Nigeria and that of the examination of livelihood strategies of indigenous nationalities in Nepal, the case of Chepangs.
Considering the activities from which the incomes are sourced in Mopani District, about 23% of household heads are without any economic engagement, while about 8% accounted for by climate-fed economic activities (farming and lumbering) and 69% are in the category of remunerative income activities. However, it is observed that more than half (51%) of the sampled households have streams of livelihood activities across the town. This is expected to boost households' coping capability with additional income. The findings reinforced the stand of [70] in their study of smallholder farmers' livelihood security options amidst climate variability and change in rural Ghana.
An examination of the educational qualification of the respondents shows a clear indication of average literacy level of households in Mopani District, i.e., about 72% are holders of Matric and various higher qualifications (out of which 34.2% hold degree and above), while 19.8% have qualifications lower than Matric. This may guarantee interest and access to useful climate change-related information that could facilitate early warning and promote propagation of adaptation actions; this finding is contrary to the findings of [1]. More than one in every five of sampled household heads have been exposed to their towns' microclimatic condition by virtue of their home locations since birth while another 51% have been for nothing less than a decade or more. A climatic condition was characterized by an upward temperature trend and erratic precipitation occasioning heat waves, flash floods, and drought. This exposure level might imply high vulnerability especially if low adaptive capacity compliments.
Household heads across the district are dominated by females, about 51% of the entire sample. The results from gender analysis were buttressed as a reflection of the situation in almost all the municipalities in the district according to [71] having more females than males. However, most momentous is Greater Giyani and Greater Letaba municipalities with both Giyani and Modjadjiskloof as administrative headquarters. The scenario was significantly adduced to domestic group's reorganization, low levels of literacy and affluence, worsened by men's mobility for greener pasture away from homes. With some comparable equilibrium in figures flanked by males and females, Phalaborwa, however, has more males than females at working age categories. This may best explain the reason for employments in mining sector by young men than women, coupled with more young women's unwillingness to take on dirty and hard jobs. However, there is a need for a more focus employment balance while gender equity is emphasized in labor.
The number of children per household, in both male and female categories, according to Table 2 depicts that households with no male child accounted for 3.6% while those with no female child are 14.7% across the district. Those categories with one male and female child are 42.5% and 45% 183, respectively, and forms the majority, while the category with three of each of the sexes claims 12.3% (male) and 4.6% (female) accordingly. This implies more male children in households across the district depicting the likelihood to reduce vulnerability. This is not only a reflection of high poverty among women per se but also the apportioned societal gendered roles and restrictions, see [72][73][74].

Household Vulnerability Levels
The results of vulnerability analyses reveal that there exists a high exposure level (92%) to changing climate in Mopani (Table 3). It was also found that Nkowankowa, Giyani, and Phalaborwa towns recorded the highest proportion (96.58%, 91.50%, and 89.28%, respectively) of households in exposed category. Being the town with the least (though still high) in the category, Modjadjiskloof accounted for about 60%. Tzaneen recorded the lowest mean exposure score (7.9%), and Phalaborwa the highest mean score (12.42%). Hoedspruit, Modjadjiskloof, and Giyani have almost the same average exposure score as summarized in Table 3. The high level of exposure across the towns reflects the consistent regional warm temperature, as over 72% have lived and exposed to the prevailed climatic condition continuously for over a decade. However, the variation in exposure levels in the towns was accounted for by the varying households' economic engagements like Hoedspruit with largest proportion of households in climate-fed venture. The towns' central areas were also dominated by impervious surface cover with attendant heatwaves and escalated runoff.
Low level of sensitivity is generally evidenced in the district with 80.24% of the sampled population to be insensitive ( Table 3). The Sensitivity Index across the six towns is equally very low, except in Hoedspruit town where 38.9% of the households are sensitive and a similar proportion is very sensitive. The situation was similar in Giyani town with 21.54% sensitive group and 4.62% very sensitive households. The low levels of sensitivity in the towns explains the low level of recorded casualties (human lives and livestock) in a general term and the low level of reported climate-related ailments. Across the Mopani District, a significant proportion of sampled households (76.2%) has capacity to adapt to climate change impacts. Tzaneen had the highest adaptability rating with an aggregate of 86.5%, while Modjadjiskloof of the six towns has the least (28%) followed by Nkowankowa with 40% (Table 3). With 31.28%, Hoedspruit had the highest mean adaptive capacity, while Phalaborwa followed with 20.51% and Tzaneen with the least. The backlogs of infrastructure and services, coupled with unemployment rate and low income level, have collaboratively compromised households' affordability level to reduce their capacity to adapt to climate change in the towns.  Table 3 is generally high with about 67.14% household across the district vulnerable. The HVI presents a similar scenario in the six sampled towns with exception of Hoedspruit town with about 11% of highly vulnerable households and 33.3% of moderately vulnerable. Tzaneen, Phalaborwa, Giyani, and Nkowankowa accounted for three out of every five highly vulnerable households. With X 2 = 2341 and p = 0.002 as summarized in Table 3, there is statistically significant variation in HVI level across the towns in Mopani District. This implies that households' levels of vulnerability differ from town to town. This may not be unconnected with the variations in socioeconomic attributes that informed households adaptive capacities as the towns were expose to similar climatic conditions. This led the study to investigate the extent of influence of households' socioeconomic factors on the household vulnerability levels.

The Nexus between Socioeconomic and Households' Vulnerability Index
The influence of households' socioeconomic factor on the HVI was examined, using the Ordered Logit Regression Model (OLRM) and the estimated results are summarized and presented in Table 4. In Table 4, the result indicates with respect to the influence of households' age on the households' vulnerability levels that a unit in upward shift in the age of respondent in Mopani District (i.e., going from 0-1) brings about 0.333 increase in the log odd of HVI. It then implies that the older the head of household, the higher the vulnerability level of the household. Likewise, in Nkowankowa town, a positive relationship exists with a unit increase in age of respondent increases the log odd of HVI by 2.345. The finding confirms the submission of [33] in their study of climate change vulnerability in Nordic region. However, R 2 value for all the models are low, i.e., in Nkowankowa (being the highest) where age only explained 19% of HVI variations. Note: Robust standard errors in parentheses; *** means p < 0.01 and significant impact at 1%, ** means p < 0.05 and significant impact at 5%, * means p < 0.1 and significant impact at 10% (Source: Author's Field Data, 2019).
The relationship between genders and HVI as indicated in Table 4 indicates an insignificant relationship existing between gender and HVI in Mopani District, with Wald Chi 2 value of 1.18, as well as the selected towns. This result is contrary to the submissions of [61] in his study of semiarid regions of Africa and Asia. However, the results reflect Wild Chi 2 = (0.78) for Tzaneen, Nkowankowa (0.05), Hoedspruit (1.10), while Modjadjiskloof, Phalaborwa, and Giyani being 0.06, 0.18, and 0.03, respectively. The regression coefficient (R 2 ) across the towns as well as the district are observed to be very low (Table 4), implying that gender only explains for a very small variation in HVI in the models, with Hoedspruit (being the highest) with only 6%. With these findings, it is inferred that since no significant relationship exists between genders of household and HVI levels, the HVI is not a function of household's genders in the towns as well as the district. This finding is contrary to that of [60][61][62].
Analyzing the nexus between marital status and HVI, the variable, marital status was dichotomized to currently married = 1 and not currently married = 0, (this includes every other category of marital status aside being currently married). The result depicts a significant unidirectional (positive) relationship where a unit increase in being married in Mopani District increases the log odd of HVI by 0.535. It is a statistically significant relationship with Wald Chi 2 = 7.95, though with low regression coefficient R 2 = 0.009. The result implies that households with married heads have higher HVI in Mopani District. The result validates the findings of [75] while examining the gendered vulnerabilities to climate change from the semiarid regions of Africa and Asia in 2016. On the contrary, the result of the examination across the selected towns shows that, except for Nkowankowa town, having significant positive relationship, with X 2 =0.618 and all the model in good fits, results for every other town is contrary and statistically insignificant (Table 4). By implication, marital status does not have any significant impact on HVI in these towns.
The ordered logit regression result for the level of educational qualification and its nexus with the HVI shows a statistically insignificant relationship in Hoedspruit, Phalaborwa, and the district. However, results for other towns depict that for Nkowankowa and Modjadjiskloof, the two variables have statistically significant and opposing (negative) relationship with X 2 = 0.253 and 0.788, respectively, while in Tzaneen and Giyani, the educational qualification has unidirectional significant relationship with HVI. This finding corroborates the findings of [76] in their assessment of the effects of educational attainment on climate risk vulnerability.
The number of male children in a household reflected an opposing (negative) relationship with HVI to climate change. Table 4 shows that a unit increase in the number of male children in a household reduces the log odd of HVI in Mopani as a district and Nkowankowa town by 0.204 and 0.363, respectively. This suggests that the higher the number of male children in a household, the lower the HVI. However, for other towns, number of male children does insignificantly relate with HVI. For the number of female children and HVI to climate change, the relationship across the entire Mopani District, as well as the selected towns, is statistically insignificant (Table 4).
Income as another important socioeconomic component was examined against HVI. Table 4 shows that income from primary source has a significant negative relationship with HVI to climate change in Mopani District as a whole, with X 2 =18.56 and p < 0.01. The result indicated that a unit increase in income reduces the log odd of HVI to climate change by 0.147 in Mopani. The result across towns shows that an insignificant relationship occurs between the two variables, with Chi 2 values of 0.22 (Tzaneen), 2.44 (Modjadjiskloof), 1.64 (Phalaborwa), and 0.36 (Giyani) and p > 0.05 in all cases. However, for Nkowankowa town, a unit reduction in income heightens the log odd of HVI to climate change by 0.244, see Table 4. This relationship is found to be statistically significant with X 2 = 5.47 and p = p < 0.01. The submission corroborates that of [40,77].
An investigation into the nexus between duration of stay and HVI shows that Nkowankowa has positive significant relationship between length of stay and HVI, with a unit increase in length of stay occasioning an upward shift in log odd of HVI by 0.244. The result implies that the longer the length of stay, the higher the HVI to climate change in Nkowankowa. However, towns like Modjadjiskloof and Giyani show negative relationships that are statistically significant between length of stay in locality and HVI to climate change. This implies that in these two towns, the longer the length of stay in locality, the lower the HVI to climate change, with the two models significant and in good fit, and X 2 = 11.26 and 7.26 with p = 0.01, respectively, for Modjadjiskloof and Giyani. With this result, one can infer that the longer the duration of stay in locality by household, the lower the HVI to climate change in both Modjadjiskloof and Giyani and vice versa. However, length of stay in locality does not significantly relates with HVI to climate change in Mopani as a district and Tzaneen and Phalaborwa towns.
In case of Vulnerability Index and the Livelihood Diversification Index (LDI) across the district of Mopani, a unit increase in (LDI) reduces HVI log odd by 0.729. This implies that the more diversified the livelihood of a household is, the lower the level of vulnerable to climate change. Similar results were obtained for Tzaneen, Modjadjiskloof, and Giyani. For these three towns, a unit increase in LDI reduces Vulnerability Index by 3.092, 9.238, and 1.29, respectively, in the three towns. On the other hand, contrary results were obtained for Nkowankowa and Phalaborwa, where Livelihood Diversification Index does not significantly affect Vulnerability Index with p = >0.05. The results obtained from the four models depict statistically significant scenario and all in good fit, with X 2 = 0.10, 3.39, 10.91, and 2.80 while the latter two models are not statistically significant. By implication, an increase in the diversification in livelihood among households will reduce the vulnerability level in Tzaneen, Modjadjiskloof, and Giyani towns. This result suggests that any innovative idea that can facilitate diversification in these towns will significantly enhance the households' adaptive capacity and invariably reduces their vulnerability levels.

Discussion
In line with the research questions, three important findings can be drawn from this work. First, climate is changing and the trend reported during the period of examinations revealed an increasing trend in both minimum and maximum temperature, while reducing erratic trend was recorded in rainfalls. The scenario was predicted to have come to stay for some decades to come. Thus, in a business as usual scenario, living and livelihood will continue to be threatened. The development is considered dangerous for households, as continuous subjection to reducing precipitation suggests an intensified risk of water stress [8]. This will sadly aggravate the plights and escalate the several millions of people in Africa that are already battling with acute water trauma, water pressure, and water unavailability. This burden would mostly be confronted by the most vulnerable (women, elderly, and children) class on one hand [11]. On the other hand, the upward trend in both maximum and minimum surface air temperature significantly evidenced local warming with varying degrees of looming costs, ranging from discomfort to heat-related ailments (asthma, malaria, measles, typhoid fever, cholera, diarrhea, and others). Climate changerelated challenges are already constituting serious threat to water resources in South Africa, overwhelming food security, compromising health, destroying infrastructure, as well as endangering the ecosystem services and biodiversity [19]. Bearing in mind the high level of poverty and inequality in South Africa, these critical impacts will be devastating if left unabated. This justifies the need for a pragmatic and integrated interventions to stem the negative impacts of the changing climate phenomenon. The strategic intervention should be eco-friendly and sustainable. Aggressive tree (fruits) planting, among other greening strategies should be embarked upon, which will not only reduce heat (adaptation) and sink carbon (mitigation) but also equally provide food (nutrients) as well as generate additional income for the households.
Second, the four estimated indexes of Household Vulnerability (Exposure index (EI), Sensitivity Index (SI), and Adaptive Capacity Index (ACI) aggregate Household Vulnerability Index (HVI)) suggested: (i) the selected towns across the district are very high in terms of EI, but relatively low in SI and relatively high in ACI, while the estimated aggregate HVI is equally high across the district as well as in the selected towns, except Hoedspruit with relatively low HVI. This proves that urgent attention is required by all stakeholders (particularly government) towards a pragmatic, integrated, and inclusive adaptation policy, programs, and projects to curtail HVI. Although, low level of sensitivity is generally evidenced in the district, nonetheless, efforts must not be surrendered because of the dynamic nature of climate. Increased investment in the basic infrastructure and services development supported by livable housing through Public-Private-Partnership will reduce the households' exposure level and further boost their adaptability of the vulnerable households to the adversities of climate change, thus reducing the vulnerability. The high level of EI in the selected towns was exacerbated by the process of urbanization that is characterized by deforestation and impervious dominance. These towns will require more resources from local and districts municipalities to create a nature-balance development, as they are all prone to becoming more urbanized in the coming decades. This manner of growth is capable of reducing household vulnerability and enhance sustainable development in the context of climate change.
Although commendable attempts by South African national government has been made to develop a National Climate Change Response White Paper, it is unable to actualize the mainstreaming of adaptation policy in daily planning practice as well as medium and long-term planning, at every level (particularly the local municipality) of governance [78]. Therefore, the need to expedite the process of mainstreaming adaptation policy into day-to-day spatial planning activities as well as prioritizing climate change matter in municipality's long-term plan is no time better than now. This require a strong political will because the five-year term of political tenure that was unfavorable to long-term issues as climate change makes them suffer lack of political attention.
Third, the key contributing socioeconomic factors to households' vulnerability were uncovered, which should provide yardstick and position for decision-making in mitigating households' vulnerability in the district and beyond. (a) The manifestation of HVI across the district is influenced decreasingly in the district and Nkowankowa town. This proves that attention is needed to enhance households' income across the district and Nkowankowa, while in other selected towns, the influence of income on HVI is yet not a serious issue. However, with majority of the household heads having financial resources at their disposals signifies a higher household affordability. A further enhancement of economic engagement of households across the district will further boost the households' affordability level and resultantly reduce their vulnerability level. (b) The number of household heads decrease with increased age implies that younger household heads are more than the older ones with majority being females. The aged were more prone to vulnerability to the impacts of climate change in the district as applicable to Nkowankowa, whereas in the cases of Tzaneen, Phalaborwa, and Giyani towns, the influence of age on vulnerability were yet not a serious issue. Thus, extra efforts are required at district municipal level to facilitate the adaptability of the old (aged) with reduction in their exposure level to climatic factors. The scenario of young household head dominating in the selected town suggests that any program that is targeted towards economic empowerment in the district with special focus on young people will impact significantly across the selected towns. This will equally be true of programs that focus on climate change adaptation and early warning propagation in the district, targeting this population category. (c) There is no evidence of statistical significance of gender's influence on the HVI in the study area. However, women accounted for more than half of the entire household in the district. The fact that females outnumber males in every age category, it is imperative that gender issues be taken seriously in climate change-related matters. The men's and gender fora that was established in the district to handle gender-related matters and to facilitate the participation of families, communities, and work place discussions but suffered from budgetary suffocation should be resuscitated, with more creative and involving efforts to bring about the desired turn around impacts. (d) Being married in Mopani District influenced higher vulnerability, but not statistically significant in the selected towns (except in Nkowankowa). With high proportion of household heads in either divorced, single, widow, or separated category may suggest an additional pressure from extra responsibility and burden to adapt for these classes of household during and after climate change extreme events, because they lack supports, as a result of absence of partners. (e) In Tzaneen, Modjadjiskloof, Giyani, and Nkowankowa higher educational qualification reduces HVI. This implies that an improved education among households can help to reduce HVI to climate change [79]. Therefore, further adjustment strategies and more work are required to get more citizen acquire additional qualifications, to enhance their resilience and reduce their vulnerability to climate change. (f) The nonclimate-fed nature of economic activities of significant proportion of households signaled a lower sensitivity. Thus, a collaborative effort by governments at various levels must be mobilized to facilitate jobs creation, through the green economic policy of the national government in renewable energy, this should be championed by the private sector. (g) Moreover, more than one in every two residents had multiple sources of income, thus, the potentials of households' enhanced affordability (with extra incomes) and opportunity to swap among income source options in case one is affected by climate change adversity are higher. This is a pointer to a stronger capacity to cope with climate change hazards and a lowered vulnerability. However, policy and intervention efforts of government should be geared towards income diversification to strengthen the citizens' buffer for cushioning the households' susceptibility during and after climate change extreme events.
The upcoming attempt at national level to scope long-term adaptation scenarios aiming at focusing on creating linkages across sectors in climate change adaptation and responses should pay more attention to active local level participation. Although some notable city-scale and project-based adaptation interventions were already executed, yet the huge institutional challenges still persist. The best time is now to make climate change adaptation responses the responsibility of the local municipality, and they should be household based. This approach will enhance the effectiveness of spatial planning-climate change adaptation designated strategy [32], by conducting assessment on location-specific basis, to cater for the peculiarities and priorities of the local community. Thus, macroplanning policies will be enabled by micro level input [80], which will in a long run mitigate the consistent failure of microlevel-based policy with no sufficient local input. Hence, there is the need to emphasize municipal governments' collaboration with private sector to drive the facilitation of private households' indigenous (microlevel) capacities for adaptation and its mainstreaming in urban planning across the district.