Climate Change and Livelihood Vulnerability in Mixed Crop–Livestock Areas: The Case of Province Punjab, Pakistan

: The mixed crop–livestock system is a primary source of livelihood in developing countries. Erratic climate changes are severely a ﬀ ecting the livelihoods of people who depend upon mixed crop–livestock production. By employing the livelihood vulnerability index (LVI), the Intergovernmental Panel on Climate Change LVI (LVI IPCC ), and the livelihood e ﬀ ect index (LEI), this study evaluated livelihood vulnerability in southern Punjab, Pakistan. The study provides a range of indicators for national and local policy makers to improve resilience in the face of livelihood vulnerability. By incorporating more major components and subcomponents, this study identiﬁes more speciﬁc challenges of livelihood vulnerability for future policy directions. It is interesting to ﬁnd that credit and cash used for crop inputs are critical ﬁnancial constraints for farmers. From the estimated indicators, this study also provides some speciﬁc policy recommendations for the four study districts of Punjab Province. These results are helpful in identifying and highlighting vulnerability determinants and indicators. Initiating and promoting better adaptive capacity and starting resilience projects for households are urgent actions required by donors and governments to reduce the livelihood vulnerability of mixed crop–livestock households in arid and semiarid areas.


Introduction
Climate change has exerted an adverse impact on mixed crop-livestock production systems and has become a global challenge, increasing the vulnerability of people who are involved in these systems. The intensity of heat waves has increased since 1950 [1], and will become more devastating in the coming decades [2][3][4][5]. Estimates show that approximately 341 climate-related disasters worldwide were recorded each year from 2000 to 2015-an increase of 50% from the previous 15 years' records-and badly affected people's livelihoods, particularly in the agriculture sector [6]. Mostly, developing countries are under the threat of erratic climate changes, which result in natural disasters including floods, droughts, and heavy rainfalls, with negative impacts on their livelihoods, leading to environmental vulnerability for farmers. Given its production characteristics, the agriculture sector is the most vulnerable and sensitive to climate change [7][8][9]. This is even worse in developing countries, particularly poor and agrarian communities, where people's livelihoods depend upon mixed crop-livestock production [10,11]. This is mainly due to the limited adaptation measures in production resources [8,[12][13][14][15][16]. As a result, climate change poses more serious challenges to local economic, social, and ecological systems [17,18]. households complemented with secondary data of climate variability such as rain and temperature. They categorized indices under exposure (climate variability, early warnings, and monetary loss due to climate events and natural disasters), sensitivity (health, water, food), and adaptive capacity (livelihood strategies, social networks, socio-demographic factors), used to construct indicators of vulnerability. Such an index-based approach established precedence over other frameworks because it is dependent on secondary data and better planned to provide a local perspective, with a context-based vision for local needs and adaptation response. Therefore, many scholars and researchers followed this approach using different contexts of natural hazards [47,[57][58][59]72,[86][87][88][89][90][91][92]. Studies that are based on primary or secondary data or both provide better vision and a platform to identify and determine the drivers of vulnerability by considering vulnerability as a socially constructed subject [93]. This study also uses household and climate data (temperature and precipitation, 2001-2010) to measure the livelihood vulnerability level through the livelihood vulnerability index (LVI), Intergovernmental Panel on Climate Change LVI (LVI IPCC ) [61], and livelihood effect index (LEI) [94][95][96][97][98].

Methods
Measurements were performed under 3 types of sets: (i) calculation of livelihood vulnerability index (LVI) through balanced weight average, (ii) estimation of LVI based on the Intergovernmental Panel on Climate Change [3], and (iii) calculation of livelihood effect index (LEI).
Hahn et al. [61] and Shah et al. [72] developed a methodology to measure the level of vulnerability caused by climate change and variability and also suggested replication for future studies. Thus, we modified and added 6 more major components (housing, land and livestock, natural resources, finance and income, knowledge and skills, and infrastructure) to better understand the sensitivity and vulnerability of households. We took into account 13 major components (socio-demographic, livelihood strategy, health, social networks, food, water, natural disasters and climate variability, housing, land and livestock, knowledge and skills, natural resources, finance and income, and physical infrastructure) and 79 subcomponents (see Appendix B). The combinations of the 13 major components with various subcomponents have equivalent influence in overall calculations [99].

LVI Indicator
We assigned equal weights to all major components to measure LVI [61]. Then, the subcomponents were calculated on a different scale to standardize their comparability and the mathematical values by using the Human Development Index [100] as follows: where S d is a subcomponent of district d, and S max and S min are maximum and minimum values, reflecting high and low vulnerability respectively. This applies to each subcomponent of all 4 districts. For instance, the average distance to travel for access to a source of water subcomponent ranged from 0 to 4 km in 4 districts. The minimum (S min = 0) and maximum (S max = 4) values were used to transform this index S d into a standardized index for the water component of the LVI (maximum and minimum values are actual data, which are not shown in Appendix C; the same for frequencies, average livestock diversity, crop diversity, etc.). For the variable that measures frequency, for example, percent of reported households with orphans, the maximum value was set as 100% and minimum as 0, whereas the subcomponent average livestock diversity index was set as the inverse of the crude indicator by assuming that household vulnerability decreases with increased number of livestock species. We created an index by taking the inverse and allocated greater vulnerability to households with fewer livestock species. We expressed the livestock diversity index by using the formula [1/(number of livestock species + 1)]. Following this logic, minimum and maximum values in Equation (1) were transformed to standardize other subcomponents. Then, we created an index to estimate the vulnerability of each major component by averaging the standardized subcomponents: where M d is each of the 13 major components for district d, index sdi shows the subcomponents for i making up each major component, and n is the number of subcomponents under each major component. After calculating the vulnerability value of all major components from each district, they were averaged by using Equation (3) to calculate the LVI at the district level: where LVI d represents the livelihood vulnerability index for district d, which is equal to the average weight of all 13 major components, and W Mi is the weight of each major component, and is estimated from the sum of all subcomponents that contain each major component and contributes equally to the overall livelihood vulnerability index [99]. We defined the following: socio-demographic (SD), livelihood strategy (L), health (H), social networks (SN), food (F), water (W), natural disasters and climate variability (NDCV), housing (H), land (L), knowledge and skills (KS), natural resources (NR), finance and income (FI), and physical infrastructure (PI). Equation (3) can be specified in detail as: where the LVI values range from 0 (least vulnerable) to 0.5 (most vulnerable). A comprehensive example for LVI estimation of major component health for Bhakkar (BHK) district is illustrated in Appendix D.

LVI IPCC Indicator
We used an alternative approach to calculate LVI based on the IPCC tools of adaptive capacity, sensitivity, and exposure [61,72]. Adaptive capacity is framed under socio-demographic, livelihood strategy, social network, knowledge and skills, finance and income, natural resources, and infrastructure. Natural disaster and climate variability are framed under exposure, and health, land and livestock, food, water, and housing are framed under sensitivity (Table 1). Exposure in this study includes a number of natural disasters events in the last 5 years. Climate variability is measured by the average number of hot months with an average temperature above 30 • C and standard deviation of monthly precipitation and maximum monthly temperature over the period 2001-2010. Equations (1)-(3) are used to estimate LVI IPCC . Then Equation (5) integrates the major components into LVI IPCC contributing factors: where CF d is contributing factors of adaptive capacity, sensitivity, and exposure for district d, M di is major components for district d, indexed by i, W Mi is the weight of each major component, and n is the number of major components in each contributing factor. Then, we estimated adaptive capacity, sensitivity, and exposure by combining the contributing factors using the following equation: where LVI IPCC d is the LVI for district d defined by the IPCC vulnerability approach and CF e d , CF ac d , and CF s d are the estimated contributing factors of exposure, adaptive capacity, and sensitivity for district d, respectively. Then, we ranked LVI IPCC from −1 (least vulnerable) to 1 (most vulnerable). Appendix E provides a comprehensive example of LVIIPCC estimation for BHK district. The LEI was constructed to identify the livelihood vulnerability of households with different types of capital or assets (human, natural, financial, social, and physical) following the sustainable livelihood approach. Based on this capital, we grouped the 13 major components in to 5 types of livelihood capital for LEI estimation (see Appendix B). We used Equations (1)-(3) to calculate the weighted average of all major components that make up above livelihood capital. Then these weighted averages were used in Equation (7) to estimate the vulnerability value of each type of livelihood capital or asset: where LEI is the vulnerability index for 5 types of livelihood capital, equivalent to the weighted average of major components that make up livelihood capital, and W Mi is the weight of each major component, and is computed by the number of subcomponents in each major component. We scaled the LEI from 0 (least vulnerable) to 1 (most vulnerable).

Study Area
Punjab is the most populated province (both people and animals) in Pakistan [29]. The geographic position of Punjab is approximately 70,000 • E and 30,000 • N in the semiarid l zone [101]. The province plays an important role in the country's economy due to fertile agricultural land with extensive irrigation networks. Punjab accounts for 56% of the agricultural land area, contributing 57% to overall agricultural GDP and approximately 74% to cereal production in the country [29,102]. Over 1970-2001, the annual mean minimum temperature of the province ranged from 16.3 to 18.2 • C and annual mean maximum temperature ranged from 29.3 to 31.2 • C. Rainfall in Punjab is linked with monsoon winds, and the province receives 50-75% of annual rainfall during the monsoon season [25].
The main sources of local livelihoods are crop farming and animal raising, but they are underdeveloped due to traditional farming practices in remote areas, low labor quality, and lack of income-generating opportunities. Approximately 67.5% of the people of Pakistan reside in rural areas and are involved directly or indirectly in agriculture. Punjab is the most populated province, accounting for 56% of the total population. The province consists of northern and southern Punjab. The northern area has better developed infrastructure and agricultural techniques than southern Punjab, while the southern area is the poorest region due to large household sizes, higher dependency and lower literacy rate, unemployment, poor infrastructure, and lack of access to markets, especially in remote areas [103]. Approximately 36% of the rural population is poor, and most of the poorest are found in southern Punjab [104,105]. Previous studies stated that farmers in southern Punjab had observed changes in weather patterns and experienced extreme temperatures and delayed precipitation during summer monsoons [25,58,106]. Therefore, this study focuses on the 4 districts of southern Punjab, where the majority of households raise livestock combined with subsistence crop production in the mixed crop-livestock production area (see Figure 1). Rahim Yar Khan (RYK) district covers an area of 11,880 km 2 , including plain and very fertile land distributed among 3 types of areas: riverine, canal irrigated, and the Cholistan Desert. The Indus River flows from the northwest side of the district and is flooded every year. The climate of the district is dry and hot with temperature reaching 49 °C and average annual rainfall of 165 mm [25]. The major crops grown in the district are wheat, cotton, rice, and sugarcane. Total livestock (buffaloes, cattle, sheep, and goats) total 4.8 million in the district [107].
Multan (MLT) district spreads over an area of 3720 km 2 , covering plain and very fertile land. This district is very close to the Chanab River and is mostly flooded due to heavy monsoon rains. Multan has an extreme climate with maximum temperature up to 49 °C and average annual rainfall of 127 mm [108]. Wheat, cotton, rice, and sugarcane are major crops in this district. Total livestock (buffaloes, cattle, sheep, and goats) number 3.9 million in this district [108].
Bhakkar (BHK) district covers an area of 8153 km 2 , including plains and desert. The western side of the district boundary touches the Indus River. Summer temperature has been recorded up to 50 °C. The main sources of livelihood are major crops such as moong, bajra, cotton, gram, sugarcane, and wheat, with total livestock (buffaloes, cattle, sheep, and goats) numbering 2.1 million in the district [109].
Dera Ghazi Khan (DGK) district includes an area of 11,922 km 2 , embracing plain and mountains. It is located in a strip between the Suleman Mountains and the Indus River and is washed out by floods and hill torrents every year. It has an almost uniform climate except in the hilly portions. The major crops grown in the district are wheat, cotton, rice, and sugarcane. Total livestock raised (buffaloes, cattle, sheep, and goats) number 1.8 million in this district [110].

Data Source
The data used in this study were taken from the Pakistan Rural Household Survey (PRHS-2012). The survey was conducted in rural areas of Pakistan, covering 3 provinces, Punjab, Sindh, and Khyber Pakhtunkhawa. Then, we selected Punjab and 4 districts, Bhakkar, Dera Ghazi Khan, Rahim Yar Khan, and Multan, by random sampling with a total of 438 households engaging in mixed croplivestock farming (Table 2, Figure 2). We gathered detailed information about major components and subcomponents (see Appendix B). Furthermore, we added secondary data on climate change such as maximum daily temperature and precipitation over the period 2001-2010 from the Pakistan Metrological Department (PMD). Rahim Yar Khan (RYK) district covers an area of 11,880 km 2 , including plain and very fertile land distributed among 3 types of areas: riverine, canal irrigated, and the Cholistan Desert. The Indus River flows from the northwest side of the district and is flooded every year. The climate of the district is dry and hot with temperature reaching 49 • C and average annual rainfall of 165 mm [25]. The major crops grown in the district are wheat, cotton, rice, and sugarcane. Total livestock (buffaloes, cattle, sheep, and goats) total 4.8 million in the district [107].
Multan (MLT) district spreads over an area of 3720 km 2 , covering plain and very fertile land. This district is very close to the Chanab River and is mostly flooded due to heavy monsoon rains. Multan has an extreme climate with maximum temperature up to 49 • C and average annual rainfall of 127 mm [108]. Wheat, cotton, rice, and sugarcane are major crops in this district. Total livestock (buffaloes, cattle, sheep, and goats) number 3.9 million in this district [108].
Bhakkar (BHK) district covers an area of 8153 km 2 , including plains and desert. The western side of the district boundary touches the Indus River. Summer temperature has been recorded up to 50 • C. The main sources of livelihood are major crops such as moong, bajra, cotton, gram, sugarcane, and wheat, with total livestock (buffaloes, cattle, sheep, and goats) numbering 2.1 million in the district [109].
Dera Ghazi Khan (DGK) district includes an area of 11,922 km 2 , embracing plain and mountains. It is located in a strip between the Suleman Mountains and the Indus River and is washed out by floods and hill torrents every year. It has an almost uniform climate except in the hilly portions. The major crops grown in the district are wheat, cotton, rice, and sugarcane. Total livestock raised (buffaloes, cattle, sheep, and goats) number 1.8 million in this district [110].

Data Source
The data used in this study were taken from the Pakistan Rural Household Survey (PRHS-2012). The survey was conducted in rural areas of Pakistan, covering 3 provinces, Punjab, Sindh, and Khyber Pakhtunkhawa. Then, we selected Punjab and 4 districts, Bhakkar, Dera Ghazi Khan, Rahim Yar Khan, and Multan, by random sampling with a total of 438 households engaging in mixed crop-livestock farming (Table 2, Figure 2). We gathered detailed information about major components and subcomponents (see Appendix B). Furthermore, we added secondary data on climate change such as maximum daily temperature and precipitation over the period 2001-2010 from the Pakistan Metrological Department (PMD).

Results
According to the findings, BHK and MLT were the most vulnerable, followed by RYK and DGK according to all three livelihood vulnerability approaches. On the other hand, RYK and DGK were the most vulnerable in the most important and sensitive subcomponents of water and health components. The spider diagram of vulnerability in Figure 3 shows LVI values for all 13 major components, which are estimated from 79 subcomponents (see Appendix C for subcomponent results) and scaled from 0 (least vulnerable) to 0.900 (most vulnerable). Table 3 illustrates index values of all subcomponents and LVI results of major components. Overall, the LVI results indicate that BHK (0.378) and MLT (0.376) have more vulnerability than DGK (0.364) and RYK (0.363). BHK is more vulnerable in livelihood strategy (0.373), land and livestock (0.337), natural resources (0.551), and food (0.523); DGK is more vulnerable in housing (0.442), knowledge and skills (0.839), and infrastructure (0.576); MLT is more vulnerable in finance and income (0.612), water (0.516), and social networks (0.360); and RYK is more vulnerable in health (0.322), natural disasters and climate variability (0.463), and socio-demographic (0.244) ( Table 3 and Figure 3).

Results
According to the findings, BHK and MLT were the most vulnerable, followed by RYK and DGK according to all three livelihood vulnerability approaches. On the other hand, RYK and DGK were the most vulnerable in the most important and sensitive subcomponents of water and health components. The spider diagram of vulnerability in Figure 3 shows LVI values for all 13 major components, which are estimated from 79 subcomponents (see Appendix C for subcomponent results) and scaled from 0 (least vulnerable) to 0.900 (most vulnerable). Table 3 illustrates index values of all subcomponents and LVI results of major components. Overall, the LVI results indicate that BHK (0.378) and MLT (0.376) have more vulnerability than DGK (0.364) and RYK (0.363). BHK is more vulnerable in livelihood strategy (0.373), land and livestock (0.337), natural resources (0.551), and food (0.523); DGK is more vulnerable in housing (0.442), knowledge and skills (0.839), and infrastructure (0.576); MLT is more vulnerable in finance and income (0.612), water (0.516), and social networks (0.360); and RYK is more vulnerable in health (0.322), natural disasters and climate variability (0.463), and socio-demographic (0.244) ( Table 3 and Figure 3).    Table 4 and Figure 4 show LVIIPCC indicator results of sensitivity, exposure, and adaptive capacity that enable mixed crop-livestock households to cope with climate change. The LVIIPCC scores range from −1 (least vulnerable) to 1 (most vulnerable). Overall, results indicate that BHK has the highest livelihood vulnerability index (0.028), followed by DGK (0.001), RYK (−0.021), and MLT (−0.013) districts. BHK also shows the highest level of sensitivity (0.369) and exposure (0.464) with less adaptive capacity (0.388), after DGK (0.378) and RYK (0.375).   Table 4 and Figure 4 show LVI IPCC indicator results of sensitivity, exposure, and adaptive capacity that enable mixed crop-livestock households to cope with climate change. The LVI IPCC scores range from −1 (least vulnerable) to 1 (most vulnerable). Overall, results indicate that BHK has the highest livelihood vulnerability index (0.028), followed by DGK (0.001), RYK (−0.021), and MLT (−0.013) districts. BHK also shows the highest level of sensitivity (0.369) and exposure (0.464) with less adaptive capacity (0.388), after DGK (0.378) and RYK (0.375).   Table 5 and Figure 5 show LEI results, which indicate that BHK (0.412) and MLT (0.403) are more vulnerable than RYK and DGK. BHK is most vulnerable in human and physical capital; in contrast, MLT is most vulnerable in social, natural, and finance capital.   Table 5 and Figure 5 show LEI results, which indicate that BHK (0.412) and MLT (0.403) are more vulnerable than RYK and DGK. BHK is most vulnerable in human and physical capital; in contrast, MLT is most vulnerable in social, natural, and finance capital.

Discussion
BHK and MLT districts are found to be the most vulnerable, followed by RYK and DGK, in all three livelihood vulnerability approaches. Furthermore, BHK households (HHs) tend to specialize in

Discussion
BHK and MLT districts are found to be the most vulnerable, followed by RYK and DGK, in all three livelihood vulnerability approaches. Furthermore, BHK households (HHs) tend to specialize in a singular livelihood (e.g., mixed crop-livestock production), usually the one that has the most income potential and/or relies on the available natural resources of the area. In fact, they are also more vulnerable to climate changes, which have a negative impact on their livelihood. However, BHK HHs are more sensitive and lack climate change adaptive measures (Table 3, Figure 3), due to lower education levels and inappropriate local government efforts and attention to this issue. The results also show that more than 89% of households in BHK district reported that they were not satisfied with current efforts by the local government to share knowledge of climate change adaptation and livelihood measures, and approximately 98% of households were not familiar with technical training programs regarding their livelihood vulnerability (Appendix C, M 51 ).
Food safety increases households' flexibility to external stress in the face of severe climate events [111], and this makes BHK HHs more vulnerable in the food component (0.523) ( Table 3) compared with other districts. The reason is that BHK is more dependent on natural resources and more sensitive to climate change; 78% of HHs there are limited to mixed crop-livestock production as major source of livelihood, while only 22% are involved in off-farm employment (Appendix C, M 91 ). Furthermore, fewer livelihood strategies and vulnerability in terms of human capital (0.501) (Table 5, Figure 5) lead to food insecurity in BHK. Following these results, policy makers, governments, and donors should take urgent steps toward economic transformation of BHK to decrease food insecurity in the future, with a focus on helping farmers transition from monoculture (crops and livestock) to economic diversification, such as growing more crops and raising more species of livestock, with more off-farm employment opportunities. This is because households in this district have less crop diversification (0.342) and fewer animal species (0.227) (Appendix C, M 93 , M 94 ) than other districts.
In contrast, MLT HHs are more dependent on off-farm work as a major source of livelihood and are less active in mixed crop-livestock production. We found that the effect of the head of the household really matters here, because 19% of heads of households are women in MLT (Appendix C, M 32 ), which is more than in the other three districts. Male heads of households bring authority and greater control over household resources, including family labor, whereas, women heads of households are linked with certain disadvantages, particularly in a country like Pakistan, due to cultural norms. It is not easy to for women without adult male support to manage farming activities alone, particularly when adult men work in non-farm sectors during peak agricultural seasons. Therefore, when men are not available to head households and work outside the community, gender issues could pose greater institutional difficulties and compound the constraints on women during farming activities. Our results also show that women are less involved in performing agricultural activities, especially in Multan district; therefore, MLT is found to be more vulnerable in terms of social network, which makes it more vulnerable in terms of climate change adaptation, social capital, natural capital, and finance and income due to women's role as head of household (Table 5, Figure 5). In fact, 71% of household members work in communities far from their homes (Appendix C, M 91 ). In this case, mostly male household heads travel for work and return home after several months, and the control of the household shifts to women. In fact, women do not have access to sources of communication media such as social networking, which results in vulnerability in social networks and livelihood [112]. Social networking is a good asset (e.g., innovations, dreams, good relations, and shared values with financial exchange), which makes a stronger information network for farmers to communicate, and in particular to deal with unpleasant and emergency events [113][114][115][116]. This finding suggests that in a country like Pakistan, where women are less active in farming activities, they should be supported by the community through female activism, such as farmers' cooperatives sharing information, and exchanging help and goods and technical support, to help farmers better understand emerging and current issues and develop better strategies and planning [117]. Women do not have authority and do not participate in social network programs in rural Pakistan due to cultural norms. In a nutshell, as MLT is more vulnerable in social networks, water, and finance and income components, we can conclude that female headship is the main reason in this district, because men are considered to play the main role in all three components in rural Pakistan due to cultural norms. To reduce the vulnerability, women should be given priority with regard to participating in farm and non-farm activities. For this purpose, local governments must create opportunities for women to access well-paid work in villages, such as supporting women-owned cottage industries, by providing training and access to credit and markets. In addition, women should be formally integrated into the value chain and efforts must be made to reduce the wage gap between men and women.
A lack of basic facilities such as health care, shelter, and clean water are noticed in RYK and MLT, which need urgent attention by policy makers to reduce this vulnerability in a timely manner. The results show that health facilities for both animals and people are located at long distances in terms of access. Because of the long distance, most poor households are not aware of their health care due to financial constraints and start treatment on their own without doctors' prescriptions when they are sick. Particularly, it is very difficult when women are pregnant and give birth on their way to the hospital before even checking in. This is an emerging issue that should be reviewed with regard to how to fix it and make sure basic health facilities are available in existing health units. Local governments should communicate with and mobilize farmers so that they do not use backward practices to treat diseases and strengthen the village health care system. For this purpose, governments should provide mobile health facilities in remote areas and announce their schedules.
The other factor is the absence of clean drinking water sources in all study areas, which makes most people vulnerable in terms of health. More important, households in RYK (16%) and DGK (5%) are drinking from natural sources (river, canals, ponds, and rain), which makes most of them vulnerable in terms of health (Appendix C, M 6 ). On the other hand, water plants or hand pumps should be installed in MLT communities in order to reduce this vulnerability by reducing the need to travel to fetch water outside the community, so that in the absence of adult men in households, women can fetch water themselves. On the other hand, for canal water availability, the irrigation department should check about allocated times per acre to overcome the water theft issue through farmers' cooperatives and make sure water is fully allocated to distributary canals to reduce the vulnerability. In a nutshell, the development of infrastructure can remove this vulnerability, particularly in remote areas, which have fewer facilities for health care and clean drinking water. This needs serious attention by policy makers and donor organizations to start clean water projects and make sure there is access to clean water, particularly in DGK and RYK.
Overall, the education and knowledge levels of heads of household are very low in all study districts (Appendix C, M 5 ). Hence, a lack of education could lead to fewer income opportunities and livelihood diversification strategies, which leads to vulnerability. As farmers are adults and it is impossible for them to go to school, as children do, local governments should be involved in training farmers in formal and informal ways. It could be possible to disseminate information about livelihood diversification to improve income by arranging skits, plays, floating activities, theaters, and corner meetings at the household level to help support livelihood resources. In fact, as the farming experience of less educated farmers increases, they become more conservative and do not adopt new farming practices. These farmers should be aware of new farming technologies with updated knowledge and climate-relevant information, with the active involvement of public and private advisory services. However, better education can help farmers to improve their ability to deal with hardship and find possible solutions under environmental vulnerability [47,59,61,118]. Rural household livelihood options indicate that the role of education and productivity in diversifying livelihood income into off-farm and non-farm activities to cope with the diverse challenges and risks of climate change can improve their livelihood in a sustainable manner by adopting higher return and sustainable livelihood strategies [119,120]. BHK and DGK HHs have a narrow perspective on the knowledge of livelihood skills and diversification (Table 3, Figure 3), with less opportunity to improve their living conditions, such as by adopting new agricultural techniques. These results are consistent with [121][122][123][124]. It is essential to share knowledge and information with farmers to cope with climatic stressors, which could be helpful in protecting their livelihood from natural disasters [25,125]. This suggests that immediate action should provide farmers with livelihood means by arranging training and knowledge-sharing activities at the farm level to improve their adaptive capacity and reduce their vulnerability due to climate change. Governments should set pathways for modifying education policies, not only for school-age children but also for adults, as mentioned above. For instance, most farmers are illiterate and rely more on informal sources for agricultural advisory services than on public or private sources [126].
The tremendous impact of shocks and the mixed crop-livestock production system in erratic climatic change events have attracted attention from policy makers and academic researchers regarding how to sustain the livelihoods of these people. Agriculture is the most severely affected sector and the most vulnerable to such disasters and climate change [16,[127][128][129], and Pakistan is listed among the countries that are extremely susceptible to climate change and natural disasters. Our results show BHK is the most vulnerable in terms of higher exposure (0.464) and sensitivity (0.369) due to climate change and natural disasters (Table 4, Figure 4, Appendix C, M 9 ). These vulnerabilities and variations are consistent with previous studies concluding that southern Punjab is more vulnerable to climate change [58,106]. Erratic climatic change and natural disasters such as floods, droughts, and heavy rainfalls badly affect people's livelihoods and increase their vulnerability. In particular, people in areas that are affected by erratic climate change and disasters need help so that that they can be rescued from falling into poverty traps, because they generally lack alternative income sources [128]. If they do not recover well, small farms in particular become more vulnerable [25,45,46], and further, they can be closed down. For example, Ahmad and Ma [130] empirically found that small farms were ignored during restoration of livelihoods by agricultural assistance (seeds and fertilizers) after super flood disasters in Pakistan in 2010, and as a result, they closed down their farms and changed their occupations in the very next season of farming. Therefore, it is very important to sustain those affected after they lose their shelter, crops, and livelihood sources to natural hazards [131,132], in order to reduce their livelihood vulnerability.
It is important to strengthen housing for low-income and poor households to improve their living conditions and build up their resilience to escape from poverty [13,59,133]. Therefore, it is suggested that governments should provide professional and technical assistance in term of housing construction, particularly in disaster-prone areas, and develop policies in cooperation with brick companies, particularly in DGK and RYK districts, as they are found to be more vulnerable to climate change in terms of maximum temperature and number of hot months (Appendix C, M 7 ). On the other hand, governments and donors should introduce heat-resistant crop varieties and animal species, along with tree plantation and other relevant information. Agriculture and livestock insurance policies could also be initiated against heavy losses due to adverse climate events [81,134,135]. For farmers' technical capacity building, governments should take the initiative regarding veterinary camps and periodic agricultural services to reduce the vulnerability of livestock and agriculture [88]. Livelihood diversification can improve income and is more effective in adapting to natural hazards and risks in the future [136][137][138][139]. Meanwhile, governments and donors should launch welfare projects to mobilize people to build welfare amenities, such as farmers with larger farms renting out their lands to landless farmers under low-payment contracts, particularly in MLT district, where 79% of households are landless (Appendix C, M 111 ). In addition, providing high-output animal breeds to farmers could be a good initiative, because less educated farmers follow their neighbors and believe blindly when neighbors have practical experience of things.
Finance and income vulnerability occur due to less income diversification and a lack of education and knowledge. However, our results suggest that another factor contributing to financial constraints is crop inputs used on credit from dealers or money lenders. For example, farmers who use crop inputs on credit have a heavy markup imposed and are bound to sell their crops to input dealers or lenders at lower than market price. Consequently, sometimes if the crop is not good and the farmer cannot pay the loan on time, the loan goes for the next crop season or year with a heavy markup imposed.
In MLT (32%) and RYK (20%), farmers have to pay previous crops' input loans if they did not pay due to crop losses or low output, and have a markup imposed until they pay the loans (Appendix C, M 131 ). Therefore, financial loan institutions can target these areas of vulnerability and provide loans with low markup to start small business and promote agricultural activities. Usually, agricultural loans are provided to agricultural landowners only, and livestock keepers are not considered for loan schemes. Thus, policy makers should develop loan policy schemes for livestock keepers so that they can diversify their livelihood and income by raising more species of animals to reduce their livelihood vulnerability.
Financial, human, natural, and physical capital have more impact on the livelihoods of mixed crop-livestock households, and these types of capital increase livelihood vulnerability [140]. Social capital shows less vulnerability in all districts despite comprising major socio-demographic and social network components (Table 5, Figure 5). The reason may be a qualitative estimation of LEI values, which indicates whether individual factors affect someone but without measuring the effect [57,92,98,[141][142][143].

Conclusions and Implications
By employing the LVI, LVI IPCC , and LEI, this study evaluated livelihood vulnerability in southern Punjab, Pakistan, using major components and subcomponents. Therefore, this study provides a range of different indicators for national and local policy makers to more accurately target ways to improve resilience in the face of livelihood vulnerability due to climate change. In addition, this research framework can be used in other countries and sectors in arid and semiarid areas. Moreover, this work proposes additional components that help in the determination of livelihood vulnerability for mixed crop-livestock households, an element of fundamental importance in semiarid regions.
The LVI approach can provide a benchmark for policy makers to evaluate livelihood vulnerability from different aspects and develop policy measures at both the macro and micro level to help allocate resources for adaptation and mitigation in the most vulnerable regions. The LVI is a useful assessment tool for critical indicators by applying equal weight to all indicators and provides spatial comparisons across regions at the household level. In addition, focusing attention on and contributing further studies and applications of this methodology is fundamental at present, and this should be included in environmental protection plans and laws worldwide.
By incorporating more major components and subcomponents, this study identified more specific challenges of livelihood vulnerability for future policy directions. From the estimated indicators, these policy targets include clean water projects, farm advisory services, locating residences away from rivers (particularly in flood-prone areas), and crop production inputs. More importantly, this study also found that the credit and cash mostly used for crop inputs are critical financial constraints for farmers, which force them to sell their crops to input dealers or lenders at lower than market price. Therefore, financial constraints should be used to identify livelihood vulnerability more precisely.
Regarding the four study districts, this study provides the following specific policy recommendations. First, in BHK, the priority is to diversify livelihoods and transform the economy from monoculture of crops and livestock and to provide more off-farm employment opportunities. Second, MLT needs to promote agricultural activities, with better access to social networks through community practices and cooperatives to obtain more professional support, and to install water plants or boreholes in communities. Third, RYK needs to initiate clean water projects, particularly in natural water utilization areas (river, ponds, canals, and rain); provide better health care services and provide mobile health unit facilities for HHs; introduce heat-resistance crop varieties and animal species; and initiate better tree plantation to cope with climate variability. Fourth, DGK requires professional and technical assistance in the construction of climate-resilient houses.
This study has scope for future research in livelihood transition and trade-offs between crops and livestock. As the cultivated land area is shrinking annually due to climate change and water shortages, cereal crop production has been severely affected in countries like Pakistan. Therefore, the livestock sector may be expected to see faster growth than crops in the coming years. In addition, it would be interesting to conduct a review to compare the methods, measurements, and results of previous studies at the country or regional level. Vulnerability refers to characteristics of human or social-ecological systems exposed to hazardous climatic (droughts, floods, etc.) or non-climatic (increasing temperature, sea level rise) events and trends.

IPCC [144], 2012a
Vulnerability is dynamic and context specific, determined by human behavior and societal organization, and influences the susceptibility of people (e.g., by marginalization) and their ability to cope with and adapt to hazards. Defining key vulnerabilities in the context of particular societal groups or ecosystem services also takes into account the conditions that make these population groups or ecosystems highly vulnerable, such as processes of social marginalization or the degradation of ecosystems.
IPCC [144], 2012a; Birkmann et al. [154], 2013a; Warner et al. [155], 2012 Societies, communities, and social-ecological systems have a limited ability to cope with the adverse consequences of climate-related hazards and to build adaptive capacities to reduce or limit their effects. Coping and adaptive capacities are part of the formula that determines vulnerability. Severe limits of coping and adaptation provide criteria for defining vulnerability, as they are core factors that increase vulnerability to climatic hazards.
Renaud et al. [156], 2010 Vulnerabilities are considered key when they are persistent and difficult to alter when susceptibility is high and coping and adaptive capacities are very low. For example, irreversible degradation of ecosystems (e.g., warm water coral reefs), chronic poverty and marginalization, and insecure land tenure arrangements are drivers of vulnerability in combination with climatic hazards.
Menkhaus [157] Conditions that make societies highly susceptible to cumulative stressors in complex and multiple-interacting systems also lead to vulnerabilities; for example, conditions of social-ecological systems that are highly susceptible to the effects of additional climatic hazards. Also, the critical dependence of societies on highly interdependent infrastructures (e.g., energy/power supply, transport, and health care) leads to key vulnerabilities regarding multiple interacting systems where capacity to cope with or adapt to their failure is low.
Fussel [67], 2007 There are three ways to understand vulnerability: (i) socioeconomic dynamics to respond any shock, (ii) risk vulnerability that consider risk experience of exposure to particular hazards, and (iii) an integrated approach that combines the two.
Turner et al. [68], 2003 The three models of vulnerability are grouped into risk hazard, pressure and release, and expanded model considering the synergies between biophysical and human systems.
Ford and Smit [73], 2004; Deressa et al. [74], 2011 Vulnerability assessment is the first step in adapting to and reducing the risk of climate change by planning programs and policies.

Social networks Maintained
Percentage of HHs that received cash aid in the past 12 months (private, government, NGOs, friends) (yes = 1, 0 otherwise) Modified Percent Aid improves recovery and reduces vulnerability, which leads to better adaptive capacity Percentage of HHs that had contact with local government office/officials for help in the past 12 months (yes = 1, 0 otherwise) Modified Percent More contact with local officials for help will highlight issues of particular areas; this increases sharing knowledge of coping strategies and strengthens adaptive capacity, with less vulnerability Percentage of HHs that borrowed or lent money (yes = 1, 0 otherwise) Maintained Percent High borrowing indicates financial stress and less adaptive capacity Percentage of HHs that contacted community leader for help in the past 12 months (yes = 1, 0 otherwise) New Percent More contact with community leaders leads to influencing individuals to deal with different events Percentage of HHs not members of any organizations (yes = 1, 0 otherwise) New Percent Information sharing and group insurance indicate less vulnerability and increased adaptive capacity Percentage of HHs that have no TV/radio/telephone at home (yes = 1, 0 otherwise) Modified Percent Access to communication media means more awareness of natural hazards and preparation Percentage of HHs that have easy access to community cooperative leader, political and government officials (yes = 1, 0 otherwise) New Percent Physical support and information sharing and experience increase adaptive capacity and reduce vulnerability

Food Maintained
Percentage of HHs that do not save grain crops (yes = 1, 0 otherwise) New Percent Higher level indicates higher sensitivity to disasters Percentage of HHs that save crop seeds for next season (yes = 1, 0 otherwise) New Percent Lower level means higher sensitivity to disasters Percentage of HHs that use agriculture production for sale of products only (yes = 1, 0 otherwise) New Percent Commercial sale of agricultural products as income contributes to less vulnerability and increases adaptive capacity Percentage of HHs that use domestic animal products as food (milk, butter, meat, eggs, etc.) (yes = 1, 0 otherwise) New Percent Using domestic animal products indicates less vulnerability and sensitivity Percentage of HHs who struggled and had food shortage in last 30 days (yes = 1, 0 otherwise) New Percent Food shortage increases vulnerability and sensitivity

Knowledge and skills Modified
Percentage of HHs not satisfied with local government efforts in sharing knowledge of climate change (yes = 1, 0 otherwise) New Percent Less local government effort to share knowledge increases vulnerability and reduces adaptability to natural hazards Percentage of HH members who have not taken any kind of vocational training (yes = 1, 0 otherwise) New Percent Training makes people better at coping with adaptive strategies HH head illiterate (yes = 1, 0 otherwise) Maintained Percent Lacking both reading and writing reduces adaptive capacity

Water Modified
HHs utilizing hand pumps for drinking water (yes = 1, 0 otherwise) New Percent Less access to good quality water increases sensitivity Average distance to water source (km) New km Shorter time reduces vulnerability and sensitivity Percentage of HHs that store water (yes = 1, 0 otherwise) New Percent Less access to fresh drinking water means higher sensitivity Percentage of HHs that have no access to canal water for irrigation (yes = 1, 0 otherwise) New Percent Less access to diverse irrigation sources means higher sensitivity Percentage of HHs utilizing water from natural resources (river, canal, wells, ponds, rain) (yes = 1, 0 otherwise) Maintained Percent Less access to fresh drinking water increases vulnerability to illness and sensitivity Percentage of HHs receiving water through public water system (water supply) (yes = 1, 0 otherwise) Maintained Percent Consistent water supply reduces vulnerability and sensitivity  Calculating steps for indices of subcomponents and major components as follows: Step 1: Repeat for all subcomponent indicators (refer to Appendix C): Step 2: Repeat step 1 for subcomponents of other major components (refer to Table 3) and then use Equation (2) to calculate indicators for all major components (refer to Table 3 Step 3: Repeat for all other major components in step 2 for LVI (refer to Table 3); for example, BHK district: Appendix E --Natural disasters and calamity variability (M 7 ) 0.464 6 Notes: Calculating steps of contributing factors of adaptive capacity, sensitivity, and exposure and LVI IPCC are as follows: Step 1: Calculate index of subcomponent indicators and major components as in step 1 in Appendix D and take the inverse of the adaptive capacity subcomponent indicators.