Livelihoods Strategies and Household Resilience to Food Insecurity: A Case Study from Rural Tunisia

Due to the decrease of household incomes, the increase of food prices, and the negative effects of climate change on agricultural production, Tunisia faces a food insecurity challenge, especially in rural and arid areas. The purpose of our research is to understand and explore household resilience to food insecurity in two villages, Selta and Zoghmar, in central Tunisia. A cross-sectional survey of 250 sample households was conducted in the villages. Factor analysis and regression models were employed to analyze the data using SPSS version 21. The results indicate that only around 36% of the households were resilient at different levels. In Selta, 62.8% and in Zoghmar 66.7% of the households were vulnerable. As indicated by the factor loadings and beta coefficients, income and food access, adaptive capacity, and the social safety net were important dimensions of household resilience to food insecurity, being positively correlated with the resilience index. However, asset possession, and climate change negatively affect household resilience. Therefore, interventions must target strategies that address the different levels of resilience reflected by the resilience estimators. These estimators were generated by focusing mainly on building farmers’ knowledge of how to face the different difficulties and challenges.


Introduction
Tunisia, a Mediterranean country bordered by the Sahara, has a variable climate with alternating periods of drought and intense rainfall that generate runoff and inflows of water. These flows are sometimes very important. Tunisia is subjected to an irregular climate and is also characterized by high aridity. In fact, climate projections indicate rising temperatures, changes in rainfall patterns, changes in the frequency of extremes, a rising sea level, and intensification of desertification. Undoubtedly, these changes will increase the country's environmental and socio-economic vulnerability [1]. This is in a country where agriculture is its backbone, economically speaking. Agriculture constitutes 10% of gross domestic product and employs 18% of the active population. Tunisia also suffers from an over-exploitation of ground water. It mobilizes most of its productive natural resources; using 4.3 million ha of its 4.5 million ha of arable land [2]. Agricultural production remains very dependent on the climatic conditions that influence rain-fed agriculture, particularly in the central and Southern regions. This makes water resources the most valuable economic and environmental 'good' for Tunisia. With less than 500 m3/inhabitant/year, the country is already in a situation of water scarcity. Tunisian agriculture is very vulnerable to climate change [3]. As a result, the area under cereal crops will decline because of this. Similarly, the tree farming area for the dried fruit species will decline in terms of its relative importance at the national level-particularly in the central and Southern regions. In addition, the number of livestock (cattle, sheep, and goats) is likely to decrease as well. The largest decreases would be felt mainly in the center and the South of the country because of the restriction of the grazing space, decline of feed production, and diseases caused by climate change-driven new pathogens [4]. Sidi Bouzid governorate, located in central Tunisia, is characterized by a large and rapidly increasing food and feed shortage, highly variable income levels, and limited natural resources-particularly arable land and water. The inhabitants are poor, most are socially disadvantaged, and disfavored in terms of infrastructural and institutional support.
This calls for urgent measures to be taken to challenge the imposed food insecurity situations. Indeed, this requires a global approach in terms of resilience, where all its dimensions should be included. However, empirical research on food insecurity in Tunisia, in general, and in the arid area, in particular, is scarce.
As a concept, references to resilience can be found in more than one discipline-art, literature, law, psychology, and engineering [5] are some examples. This reflects its wide variety of meanings and definitions [6]. It was considered as new in the social sciences [7], but has now become one of the keywords in recent scholarly and policy debates on food security [8]. At the Food and Agriculture Organization (FAO) in 2008, a well-built formulation of the concept for studying food systems and food security was started following the novel contribution of Alinovi et al. [9].
The objective of this study is to identify and analyze the main dimensions of household resilience to food insecurity and compare the results between two villages, Zoghmar and Selta, in Sidi Bouzid. The important contributions identified in this study can serve policy makers to build household resilience to food insecurity in rural and arid areas.

Study Area
Sidi Bouzid governorate is in central Tunisia. It covers an area of 7405 km2 and it has an arid climate with an annual rainfall of between 200 and 300 mm. It is divided into 12 districts, 10 municipalities, 111 sub-districts, and 12 rural councils ( Figure 1). Sidi Bouzid governorate is characterized by low levels of economic activity, a high incidence of drought, and a high concentration of rural population (75%). For many years, it was disfavored in terms of infrastructural and institutional support besides having limited natural resources, particularly arable land and water. A large number of small farmers derive most of their family incomes from a barley/livestock-based system. Sidi Bouzid is typical of the barley-sheep belt in central Tunisia and across North Africa. Barley is a drought tolerant crop and local sheep are very much adapted to the prevailing climate conditions. Sheep fattening practice is quite profitable in the region. According to national statistics, Sidi Bouzid is ranked number one nationwide in terms of collected milk, with a contribution of 293,000 L/day (11 to 15% of the national volume). The dairy cattle population is about 35,000 cows, which are owned by the majority of the small producers (≤6 cows). The daily production volume per farm is between 20 to 60 L [10].

Data Collection and Analysis Procedure
The data used in this empirical application have been collected through two random sample surveys conducted in Zohmar and Selta sub-districts of Sidi Bouzid Governorate. Primary data were obtained through farmer interviews using a structured questionnaire. A socioeconomic questionnaire was also used to collect information required to study characteristics of the communities within the geographical area of interest. From a Rapid Rural Appraisal (RRA), field visits and expert consultations have a simple random sampling approach was used to select a representative sample of 250 households within the two sub-districts (Zoghmar, 105 and Selta, 145).
The questionnaire was designed and developed by the socio-economic teams both at the International Center for Agricultural Research in the Dry Areas (ICARDA) and the Tunisian National Research Agronomic Institute (INRAT). The questionnaire was pre-tested, modified, and implemented within the two regions. Data covers information on household head's background, demographic characteristics, financial, physical, natural and social capitals, household gross income, agricultural production, farmers' practices and inputs used, utilization of farm production, access and management of water resources, livestock production and management, food security, and access to credit. Information on system vulnerability and local coping mechanisms used by households were also collected. The analysis was conducted using SPSS® (IBM) version 21.
The approach adopted was characterized by the steps taken to generate the RI model. These steps included developing a theoretical framework, identifying and developing relevant variables, standardizing to allow for comparison, weighting and aggregating the variables, and conducting uncertainty measures to gauge the robustness of the variables.

Empirical Model
When it comes to household resilience, it has at its core the notion of complex dynamic processes that aptly describes the nature and dynamics of vulnerability and changes [11]. As for its origins, the term "resilience" implies strength and resistance, but more recent applications in ecology, socioecological systems, disaster management, and urban sustainability, emphasize that resilience requires flexibility, learning, and change [12]. Actually, building resilience in the context of a household and/or a community requires full participation at the community level. It is built through community decision making identifying the appropriate interventions and determining what constitutes success at the community level. "A resilient system has the capacity to respond positively to change, maintaining or improving function" [13]. In fact, in order to define resilience, it is necessary

Data Collection and Analysis Procedure
The data used in this empirical application have been collected through two random sample surveys conducted in Zohmar and Selta sub-districts of Sidi Bouzid Governorate. Primary data were obtained through farmer interviews using a structured questionnaire. A socioeconomic questionnaire was also used to collect information required to study characteristics of the communities within the geographical area of interest. From a Rapid Rural Appraisal (RRA), field visits and expert consultations have a simple random sampling approach was used to select a representative sample of 250 households within the two sub-districts (Zoghmar, 105 and Selta, 145).
The questionnaire was designed and developed by the socio-economic teams both at the International Center for Agricultural Research in the Dry Areas (ICARDA) and the Tunisian National Research Agronomic Institute (INRAT). The questionnaire was pre-tested, modified, and implemented within the two regions. Data covers information on household head's background, demographic characteristics, financial, physical, natural and social capitals, household gross income, agricultural production, farmers' practices and inputs used, utilization of farm production, access and management of water resources, livestock production and management, food security, and access to credit. Information on system vulnerability and local coping mechanisms used by households were also collected. The analysis was conducted using SPSS®(IBM) version 21.
The approach adopted was characterized by the steps taken to generate the RI model. These steps included developing a theoretical framework, identifying and developing relevant variables, standardizing to allow for comparison, weighting and aggregating the variables, and conducting uncertainty measures to gauge the robustness of the variables.

Empirical Model
When it comes to household resilience, it has at its core the notion of complex dynamic processes that aptly describes the nature and dynamics of vulnerability and changes [11]. As for its origins, the term "resilience" implies strength and resistance, but more recent applications in ecology, socio-ecological systems, disaster management, and urban sustainability, emphasize that resilience requires flexibility, learning, and change [12]. Actually, building resilience in the context of a household and/or a community requires full participation at the community level. It is built through community decision making identifying the appropriate interventions and determining what constitutes success at the community level. "A resilient system has the capacity to respond positively to change, maintaining or improving function" [13]. In fact, in order to define resilience, it is necessary to clarify conceptual and theoretical relations with several concepts-vulnerability, stability, sustainability, adaptability capacity, robustness, and transformability [14].
The choice of the resilience approach is not arbitrary. Developing and refining models in food security research and studies has resulted in resilient systems that also involve households as components. These components are not passive and simply sub-systems, but also, they are adaptive complex approaches [7]. These have as origins the vulnerability approach and the many flaws and weaknesses such models involve [9]. Resilience thinking is a systemic thinking [14]. Resilience thinking therefore necessarily includes a system thinking, while a system thinking does not imply application of the concepts from resilience thinking. The household is an essential component of the food system that is in continuous interaction with its environment (economical, ecological, cultural, institutional, technological, etc.). In the case of food insecurity, households make important decisions to improve their resilience, thereby reorganizing the functioning of the food system. Improving the food system resilience necessarily depends on how decision-makers decipher the feedback of households to shocks and how they can incorporate it into their intervention measures.
A helpful framework for household resilience to food security has been provided by researchers at the FAO, as mentioned earlier. The analytical framework used in this paper is based on it ( Figure 2). Household resilience to food insecurity is a latent variable observable per se [15]. This latent variable is a function of eight other variables, themselves latent, and estimated from a given number of observed variables. Figure 1 shows the analytical framework that represents a rational measuring of the households' resilience to food insecurity. This model of a household's resilience index can be expressed mathematically as: where RI is the resilience index, AC is the adaptive capacity, AP is asset possession, SSN is the social safety net, IFA is income and food access, S is stability, and CC is climate change. Burton [16] argues that there are two fundamental and important tasks or steps to be taken when it comes to furthering resilience scholarship. The first is the measurement tools in order to better understand the factors contributing to resilience and the second is the selection of the relevant factors of resilience to a particular shock [14]. to clarify conceptual and theoretical relations with several concepts-vulnerability, stability, sustainability, adaptability capacity, robustness, and transformability [14]. The choice of the resilience approach is not arbitrary. Developing and refining models in food security research and studies has resulted in resilient systems that also involve households as components. These components are not passive and simply sub-systems, but also, they are adaptive complex approaches [7]. These have as origins the vulnerability approach and the many flaws and weaknesses such models involve [9]. Resilience thinking is a systemic thinking [14]. Resilience thinking therefore necessarily includes a system thinking, while a system thinking does not imply application of the concepts from resilience thinking. The household is an essential component of the food system that is in continuous interaction with its environment (economical, ecological, cultural, institutional, technological, etc.). In the case of food insecurity, households make important decisions to improve their resilience, thereby reorganizing the functioning of the food system. Improving the food system resilience necessarily depends on how decision-makers decipher the feedback of households to shocks and how they can incorporate it into their intervention measures.
A helpful framework for household resilience to food security has been provided by researchers at the FAO, as mentioned earlier. The analytical framework used in this paper is based on it ( Figure  2). Household resilience to food insecurity is a latent variable observable per se [15]. This latent variable is a function of eight other variables, themselves latent, and estimated from a given number of observed variables. Figure1 shows the analytical framework that represents a rational measuring of the households' resilience to food insecurity. This model of a household's resilience index can be expressed mathematically as: where RI is the resilience index, AC is the adaptive capacity, AP is asset possession, SSN is the social safety net, IFA is income and food access, S is stability, and CC is climate change. Burton [16] argues that there are two fundamental and important tasks or steps to be taken when it comes to furthering resilience scholarship. The first is the measurement tools in order to better understand the factors contributing to resilience and the second is the selection of the relevant factors of resilience to a particular shock [14]. A consensus on how to measure resilience has not been reached among researchers. While some researchers employ qualitative techniques [18], others construct an RI and statistically measure the relative status of the household resilience. In this current paper, quantitative techniques are employed to measure the status of a households' resilience to seasonal food insecurity. The main tool used to perform the analysis is a series of multivariate techniques (principal component analysis, factor analysis, and optimal scaling). Bartlett's [19] scoring methods were used to analyze the data. If the observed variables are measured on a continuous scale, principal component analysis is used to select the parsimonious variables used to estimate the first latent variables. Otherwise the optimal scaling method is used. The initial phase of variable selection was based on the literature and the researchers' experience of the study area.
In addition to its use in variable selection, optimal scaling was used to estimate the factor scores of the latent dimensions from dummy variables. Factor analysis, using principal axis factoring, was used to estimate such scores from the continuous variables. The variance accounted for in each factor in the solution was computed from the eigenvalue used to estimate the weighted scores of each latent variable and the overall RI for each household. This is, per se, also a latent indicator. That is the indices of each household were generated weighting the scores using the following simple formula: where RI n is the resilience index of the nth household, ∂n is the variance explained by each factor, factors1, 2, 3 . . . n are the respective factors generated by the factor analysis representing each latent dimension.

Estimation of the Latent Indicators
Before building the RI, some latent indicators were estimated based on the data survey. As explained in the conceptual framework, identification of the relevant variables is a basic requirement for designing effective future policy and intervention objectives in Sidi Bouzid. For this purpose, two steps of analysis were undertaken. First, relevant multivariate analyses were run (i.e., principal components/factor analysis for continuous variables and optimal scaling for non-continuous ones) using the available indicators of each latent dimension separately [17]. Then, relevant variables were selected based on the factor loadings and other statistical criteria. These criteria included Kaiser-Meyer-Olkin (KMO) statistics of sampling adequacy and Bartlett's test of sphericity, communalities, and variance explained by the factor generated. Secondly, these variables were used to estimate the respective final latent dimensions that are later used to estimate the overall RI. In light of these, the following observed variables were identified, defined, and used to estimate their respective latent dimensions.

Income and Food Access (IFA)
It is important that the latent dimension of the household's resilience is related directly to the household's capacity to access food. Traditionally, this has been measured by income [17]. However, in this current study, our focus was set not only on the existing literature measuring the IFA, but also on other generated indicators that seemed essential to reflect the reality of the farmers facing food insecurity.
The carefully chosen indicators estimating IFA are: • Total income per day per person (TIPDPP). This is a continuous included observed variable measuring the IFA. Its calculation is based on the different sources of income (livestock income, crop income, and other income) • Storage after harvest (SAH). The quantity of grain held in stock (quintals (1 quintal = 100 kg)) • Diversity (D). This measure is based on the number of days/month, on average, that the household consumes certain items, such as wheat, barley, other cereal (rice, sorghum, millet, maize, etc.), roots and tubers (potatoes and sweet potatoes), legumes, sugar or sugar products, groundnuts or cashew nuts, etc. • Index adaptation strategy (IAS). This was estimated from the average number of coping strategies available to the household when it lacks food-purchasing food on credit to overcome daily food insecurity, borrowing food, seeking help from relatives or the community, adults eating less so that children can eat more, consuming seeds saved for next season, etc. Its value lies between 0 and 1 where 0 equates to no coping strategy and 1 equates to a maximum of 11 coping strategies • Food access (FA). This is a continuous observed variable that was estimated as an average of the responses to six questions that sought to understand the household's ability to access food. Its values lie between 0 and 3. The question was concerned mainly with different aspects of the households' abilities to access food. The questions addressed such issues as: • How often was the household worried that it would not have enough food to eat for the month? • How often was the household unable to eat its preferred food because of a lack of income? • How often did the household eat fewer meals per day because of not having enough food?
Due to the high correlation among them, all these IFA estimators can produce a latent variable that fits the common pattern in the data. They all have continuous scales. A factor analysis was run to estimate the IFA using the principal axis factoring method and the scoring method suggested by Bartlett [19], using SPSS version21. Considering the communalities and the initial commonalities before rotation we note that they are all above 0.30, which is good (Table 1). Two factors were extracted and based on the analysis of the eigenvalues. We notice that the total variance explained shows that these two components explain 57.138% of the total variance. The relative size of the factor loading of each variable has an important policy implication (i.e., the higher the loading, the more important it is, and the greater the attention that should be paid to it). These indicators play important roles in estimating the IFA dimension, although they differ in their correlation coefficients. As expected, the factor loadings and correlation coefficients of all variables are different. Only the roles of D and SAH seem to have less correlation with IFA, and therefore, are less representative of it. That may be explained by the farmers not considering diversity but sticking to common Tunisian meals (milk and cereal products) to be food secure (D). Additionally, they rarely consider the harvest savings as a source of food for coping with food insecurity (SAH). Crop harvests are mostly sold to meet household needs. IAS is the factor most correlated with IFA. The households mention using different strategies to overcome food insecurity. The most used strategies were selling small ruminants (67.1%), purchasing food on credit (58%), borrowing money from relatives and neighbors (53.1%), and borrowing food and getting help from relatives and the community (51.7%). TIPDP, FA, and IAS were captured by Factor 1, while SAH and D were captured by Factor 2. The relative size of the factor loading of each variable has, therefore, important policy implication (i.e., the higher the loading, the more important it is, and the more attention that should be given to it). In general, as suggested by the statistical criteria mentioned above, the factor scores can be used to estimate the IFA as follows: The result obtained from this equation was further used as one of the latent indicators in constructing an overall RI for each household.

Adaptive Capacity (AC)
Scholars have analyzed AC from different angles [20].It is considered as an essential dimension of resilience that not only measures the household's ability to adapt to facing shocks [14], but AC is also critical in reducing food insecurity. It refers to the conditions that enable people to [21]:

•
Anticipate and respond to changes • Minimize, cope with, and recover from the consequences of change • Take advantage of new opportunities. All the variable values are based on a counted measurement scale and factor analysis was run, using the principal axis factoring method, to estimate the AC latent variable. Table 2 shows that the measure of the adequacy of the KMO sampling is 0.524 and Bartlett's sphericity test is significant (value p = 0.000, chi squared = 28.798). These results indicate that our sample is suitable for this factor analysis. Two factors are generated from the analysis representing 60.5% of the total variance. The first factor has an eigenvalue of 1.365 that represents more than 34% of the variance and the second has an eigenvalue of 1.055 and represents a variance of 26%. The ID and ACS were captured by the first factor and NOFI and the NHA were captured by the second. We also note that all the variables are positively correlated with the latent variable, AC. The off-farm income and the number of farm and household assets are the ones most correlated with the latent variable, AC. This seems to be explained by rural households being more likely to work in a sector other than agriculture or to sell farm assets to address food insecurity. We also notice that, based on the KMO measure of sampling, almost 52.4% of the variability can be explained by the underlying factors. In other words, it indicates that our sample is suitable for this second factor analysis. Additionally, the Bartlett test is significant (p = 0.000, which is less than 0.05). This means that the variables are sufficiently highly correlated to provide a reasonable basis for factor analysis, as in this case. The initial communalities are above 0.30, which is good. Only two factors were extracted based on the analysis of the eigenvalues. The total variance explained is 60.50% of the total variance. The equation estimating the AC latent variable is: AC = (0.34122 × Factor 1 + 0.26378 × Factor 2)/2 (4)

Social Safety Net (SSN)
The SSN can be considered as the system's capacity to mitigate shocks (in our case, food insecurity). The idea is that the more the households practice these indicators, the stronger they are socially integrated, and the more resilient they become [14]. Accordingly, the following indicators of SSN were identified in the present study, which seeks to estimate the household's resilience index (RI): The variables are not continuous, so an optimal scaling was adequate to estimate the latent variable, SSN, by transforming the observed variables. As expected, the loadings and correlations of all variables with SSN were positive (Table 3). The analysis shows that the total Cronbach's alpha is about 0.8 indicating a high level of internal consistency for our specific sample. Moreover, the first factor generated by this procedure is very satisfying and can be used to represent SSN. This factor alone accounts for about 28% of the variance. By adding the second, both factors account for about 49% of the total variance. In addition, with the exception of the variable MSOA, the factor loadings of all observed variables are greater than or equal to 0.3. BFHRC, SCLR, and BMRN were captured by Factor1, while AFAIO was captured by Factor 2.
Additionally, we note that the correlations of all the variables with SSN are positive. We ought to point out that the variable BMRN is most correlated with the latent variable SSN. The first solution that households favor is to borrow money from relatives or neighbors. The acceptance of food aid is also strongly correlated with SSN, which shows the efforts of non-governmental and international organizations, such as FAO, UNDP, etc., in their poverty alleviation programs. Being an MSOA has little correlation with SSN, which reflects the lack of an associative or political fabric in the rural area.
The SSN score for further estimating overall household resilience was estimated as:

Assets Possession (AP)
AP is a key element of household resilience because the more assets a household has, the more shocks it will be able to weather, thus improving its resilience. Five observed variables were used to estimate the AP components that we considered important for an agricultural household, especially one in a rural environment where sheep husbandry is dominant. The number of livestock is a kind of security guarantee for households for coping with food insecurity or climatic shocks. The households can sell animals to meet their needs (destocking) and the flock size can be rapidly increased when the household financial status improves. The large demand on animal products is also an important driver for households to maintain sheep. The larger the size of the herd, the greater the resilience of the household was supposed to be. Agricultural and non-agricultural assets are also considered means that can be divested when needed. The value of the house is also an important indicator that gives an idea of the standard of living of rural households. The higher the value of the house, the more the opportunities for the household to circumvent food insecurity. Finally, land ownership is also a relevant indicator of a household's ability to cultivate or rent land.
The AP latent variable was estimated from five different observed variables: • Total owned land (ha) (TOwLH): This is an interesting variable, knowing that only 2.9% in Zoghmar, and only 7.6% in Selta are landless farmers • Herd size (HS): was included because it is an also an intriguing observed variable. Given the lack of water resources, the system of production is mostly dependent on raising livestock, and therefore, the practice of sheep herding [10]  Running the factor analysis estimating AP with principal component analysis, promax rotation, and Barlett's scoring method, two indicators were extracted. Table 4 shows that the measure of the adequacy of the KMO sampling is of the order of 0.566 and the test of Bartlett's sphericity is significant (value p = 0.000, chi squared = 114.417). These results indicate that our sample is suitable for this factor analysis. The first factor generated is satisfactory and can be used to represent the latent AP indicator. This factor alone accounts for about 36% of the variance, while the three factors together represent about 58%. Therefore, these two factors were used to estimate our latent variable AP. In addition, the factor loads, and the correlations of each observed variable are positive. The first three variables, TOwLH, HS, and HVal, are captured by Factor 1, while AA, and NAA are captured by Factor 2. The variables HS, AA, and NAA are the most correlated with the latent variable AP. In fact, households are not in favor of selling or renting their home or land. Rather, they are in favor of selling a certain number of animals or a few agricultural or NAA when needed. The following formula is used to estimate the AP: AP = (0.358 × Factor1 + 0.2182 × Factor 2)/2 (6)

Stability (S)
The fifth latent variable in our study to estimate the RI is stability (S). It was estimated from a set of five different variables. In marketing research, survey research, and social and behavioral sciences research, categorical data is dealt with most frequently [22].
The observed variables chosen to reflect the stability latent variable are: All four variables are measured on continuous scales. A component factor analysis was run as the most suitable tool for estimating the latent variable, S. Two factors were extracted. Table 5 shows that the first factor generated was satisfactory in explaining about 44% of the variance. However, the contribution of the second factor was also significant, as together, they accounted for about 65.54% of the total variance. The variables FS and NHHOF are the most correlated with the latent variable, S. Having a large family offers more stability when facing food insecurity, especially if the household has members who work off-farm. In addition, a higher number of females in the household offers greater stability, especially in rural areas where agricultural work is mainly done by women.

Climate Change (CC)
This is one of the most important dimensions of household resilience in our study. Much of the previous literature on resilience has missed or not even considered it as a main latent variable in estimating resilience. Therefore, the CC indicator is a very interesting one to consider-not only the water resources scarcity in Sidi Bouzid, but also the severe weather conditions encountered.
The inter-annual precipitation variability (CV%), which ranges from 38% to 59%, is, on average, 51% [17]. The included variables were chosen using several questions reflecting observed households' shocks and the CC realities encountered and observed (for the previous decade) are: Not all of these variables are measured on a continuous scale. They are all dummy variables that take a value of one or two. They were collected using questions about previously enumerated major shocks that had been encountered by the households in the past 10 years.
As a result, optimal scaling was used to estimate the CC latent variable as it relates to the dummy variables. Table 6 shows that the first factor explains about 33% of the variation. While it is satisfactory and could be used to represent the CC, the contributions of all three factors were significant, accounting for about 61% of the total variation. We also noted that all the variables are positively correlated with the latent variable, CC. The variables ODS, OTFH, OTR, and OHS are the ones most strongly correlated with the latent variable, CC.
All variables have factor loadings that are positive and close to or greater than 0.3. OD is captured by Factor 3, OF is captured by Factor 2, while the rest and most of the observed variables are captured by Factor 1.

Estimating an Overall Resilience Index (RI)
The overall RI for each household is estimated by applying principle component analysis and using the factor method-taking the latent variables estimated earlier, which have been normalized/standardized using factor analysis or optimal scaling (i.e., mean = 0 and variance = 1). The six building blocks (latent factors) were analyzed using Bartlett's method for factor scores saved in the SPSS active data set in order to calculate the overall index. The model shows that three factors together accounted for about 72% of the total variance ( Table 7). The KMO statistic for the model was 0.541 and Bartlett's test of sphericity was significant (p-value = 0.000 with chi-square = 203.538). RI can be successfully written as: RI = 0.29455 × Factor 1 + 0.26124 × Factor 2 + 0.16798 × Factor 3 (9) Consequently, the household RI was built by calculating the weighted sum of the three factors (the weights are proportional to the variance explained by each factor), as suggested by DiStefano [23].
RI, like all composite indices, is a relative measure of the level of household resilience, or vulnerability, to food insecurity. It is a composite index that simply shows the level of resilience among the studied subjects.
In addition to suggesting its own cut-off points to categorize households into different resilience levels, this paper has two main objectives that were successfully accomplished. The first is to identify the important building blocks and determinants of household resilience. The second is to examine the different resilience levels for the two villages, Zoghmar and Selta. The second objective can be used to consider what are the policy implications and the strategies that can be implemented.
The paper uses the following ranges of RI scores, which are randomly proposed and used as a grouping method for four levels of resilience:   Table 9 shows that the households in Zoghmar are more vulnerable than those in Selta. However, in Zoghmar there are more highly resilient households (14.3%) than in Selta (9.7%). The two villages represent different production systems. In Zoghmar, biophysical characteristics facilitate the diversification of agricultural activities. For instance, sheep production is the main source of revenue for more than half the local population, whereas in Selta farmers, they rely more on goats and beehives because of the proximity of Mount Mghila. When it comes to plantations, olive trees are more common in Selta because they require less labor and its soil is stony soil. In Zoghmar farmers rely more on cereals which provide them with usable by-products, such as straw, insuring livestock feed is available throughout the year. Most of the households in both villages are vulnerable-62.8% of our sample's households in Selta and 66.7% in the Zoghmar district. This point out the importance of the political and strategic decisions that should be taken in the two study villages.

Relative Importance of the Latent Dimensions of Household Resilience
It should be acknowledged that the latent dimensions do not equally contribute to household resilience. Some have more significant impacts than others, although the role of each is still crucial in estimating resilience. To identify the effects of each latent variable on household resilience, a simple regression analysis was run using the ordinary least squares (OLS) algorithm. This gives the effects in terms of a standardized (β) coefficient. The household resilience scores were used as dependent variable and five of the latent variables, as explanatory variables. The latent variable S was removed from linear regression to avoid a problem of multicollinearity between the variables.
The following is the empirical model used to predict household resilience (RI).
where RI is household resilience; ∂ is a constant; β1-5 are the coefficients of each variable, is an error term representing the negligible information of the variables used to estimate resilience but excluded from the model above. Our first table of interest (Table 10) shows that the R value, which represents the simple correlation, is 0.966, which indicates a high degree of correlation. In addition, the R 2 value indicates how much of the total variation in the dependent variable, RI, can be explained by the independent variables, IFA, AC, SSN, AP, and CC. In our case, 93% can be explained, which is large. SPSS recommends examining the β coefficients for the relative importance of each variable in determining a dependent variable (in this paper Zoghmar's and Selta's households' resilience to food insecurity). Accordingly, Table 11 shows, that, irrespective of the negative signs, AC was the most important dimension. This contributed more (β = 0.614) to the regression model. It was followed by SSN (β = 0.337), IFA (β = 0.328), AP (β = −0.035), and finally, CC (β = −0.283). A 0.614 unit increase in the household's AC increases the households' resilience by one standard deviation. A 0.283 unit decrease in climate change increases the households' resilience by one standard deviation. The logic behind this is that the return to equilibrium after facing climatic shocks and challenges remains temporary. Neither a long-term solution nor a long-term program has ever existed to allow farmers to adapt the farming systems in the two villages to the changes in terms of expected weather. In addition, AP, like CC, is characterized by a negative β. In fact, a decrease of 0.035 unit in AP increases the resilience of the households in the two sub-districts by one standard deviation. This seems to be explained by two main reasons. The weak AP of smallholders is exacerbated by the under exploitation of farm assets because of poor management. In this direction, the lack of market information is a major handicap for good asset management. The relative importance of this dimension of a household's resilience to seasonal food insecurity has clear policy implications because it minimizes intervention costs. Policy makers and implementers can emphasize the important dimensions of resilience when planning to improve the resilience capacity of households.

Concluding Remarks and Policy Implications
One of the most important conclusions based on this study is the significant need for concrete policy implementations in Selta and Zoghmar. Facing more than one challenge, especially with the accelerated deterioration of natural resources, added to the exhaustion of the development model based on the agricultural sector's [24] resilience, which means that improving governmental policy is crucial to facilitating the development of the two zones. The major focus should be to work on the elaboration and successful implementation of the strategies along with education. The essential goals should be strengthening the households' knowledge of the coping strategies to deal with food shortages and their awareness of the climatic challenges. Different coping strategies and approaches should be adopted taking in account the resilience estimators analyzed in this study. This reflects the important intensity of the call regarding appropriate policy interventions based on the current research. In this regard, AC played the most significant role in household resilience to food insecurity in the study areas. This latent variable is highly correlated with number of off-farm income sources (NOFI) and number of household assets (NHA). Strategies such as creating opportunities for off farm income generating activities especially for youth should be included in any political program and interventions in these areas. SSN and IFA are found to play a moderate role in the household resilience to food insecurity. Regarding the observed variables estimating SSN, we notice the importance of informal solidarity networks among community, especially in terms of borrowing money and food. However, the household heads membership to professional organizations (association, political group, ONG, etc.) was very week. This result suggests the need to focus on developing producers' organizations and household associations, which will contribute to improve the SSN resilience dimension. The best example for those associations are the Mutual Agricultural Services Organization (SMSA) where they could play a strong role in helping poor households through ensuring a simple and attractive credit system, allowing families to purchase indispensable household needs and agricultural inputs, lowering the costs of technologies and services and fostering decent rural employment.
Among strategies that can help the poor actors to do better and improving IFA, is creating local markets close to the poor households, located in in less-favored areas, particularly mountainous areas (i.e., Selta) and characterized by a lack of proper agricultural infrastructure to support farming. Another vital strategy that should be considered is providing micro-finance loans and smart subsidies are proving vital facilities for many farmers in enabling them to manage risk, make investments in water harvesting technologies, diversification of food products, and creating opportunities for farm income generating activities by promoting private-label products (i.e., Sidi Bouzid meat lamb).
As regards AP, the policy recommendation suggests to emphasize on the good management of assets. Strategies related to improving farmers' know-how on land and livestock management are strongly recommended, especially under the dynamic threat of climate change. Furthermore, findings indicate that CC has a significant and negative impact on household's resilience to food insecurity. Autonomous household's adaptation was insufficient to adequately address the threats posed by climate change. Interventions could include, for instance, subsidizing more of what increases households' resilience, such as varieties that are highly drought adaptive, or by implementing programs (drought preparedness plans, soil erosion and water harvesting plans, etc.) that target the farmers' knowledge of how to face climate change difficulties in the best possible ways.
It is also crucial to consolidate public investments to improve the quality of human and animal health services in study area A, One Health approach addressing simultaneously human, animal, and environment health should be an important pillar of the coping strategy to increase resilience of the households and mitigate effects of climate change.
The study shows the relevant role of resilience in synergizing the implementation of the 2030 agenda, especially the first of the 17 Sustainable Development Goals (SDGs) related to the eradication poverty considered as the key challenge to achieve sustainable development, integrating environmental, and social and economic dimensions. A further research including additional resilience indicators related to SDGs in different agro-ecological zones will be useful in improving understanding of food insecurity impacts and vulnerabilities and in formulating, evaluating, and improving response strategies.