Original Research
Frontiers in Environmental Economics ·
Volume 5 - 2026 | https://doi.org/10.3389/frevc.2026.1781478
Understanding climate vulnerability and economic impacts on rural livelihoods: lessons from Narowal, Punjab
- Laila Shahzad1*corresponding author
- Maham Saleem1
- Faiza Sharif1
- Arshad Ali Shedayi2
- Muhammad Umar Hayyat1
- Mahammad Farhan1
- Maryam Dogar3
- 1Sustainable Development Study Center, GC University, Lahore, Pakistan
- 2Karakorum International University, Gilgit, Pakistan
- 3National University of Science and Technology, Islamabad, Pakistan
Abstract
Climate change is directly affecting the agriculture sector, and the underlying vulnerability of farmers is often overlooked. Despite agriculture's vital role in economy of Pakistan, localized assessments of farmers' climate vulnerability and adaptation remain limited. This study assesses the vulnerability of rice–wheat farmers in District Narowal, Punjab, to multi-hazard climatic variability using a climate vulnerability index (CVI). Data was collected from 300 households through door-to-door surveys and focus group discussions across two Tehsils of District Narowal. Overall, IPCC-CVI was accounted as 0.26 and 0.91 for Tehsil Narowal and Tehsil Shakargarh respectively. These findings reveal that extreme temperatures, erratic rainfall, and recurrent floods significantly disrupt crop cycles, reduce yields, and exacerbate food insecurity. Biological hazards, including pest outbreaks and resilient weeds, further compound risks, while limited financial, human, and physical capital restrict adaptive capacity. Socio-economic factors such as age, education, land size, and ownership rights emerged as key determinants of vulnerability, making Tehsil Narowal slightly more vulnerable than Tehsil Shakargarh. Farmers exhibit moderate adaptation through crop diversification and altered sowing dates; however, sociocultural barriers, particularly for women, hinder technology adoption. The study underscores the urgent need to mainstream climate resilience into agricultural development projects and provincial policies to safeguard rural livelihoods and stabilize Pakistan's agricultural economy.
Keywords
1 Introduction
Climate change and its increasing variability have become major environmental and social challenges, driving a rapid rise in climate–induced hazards and their destructive impacts worldwide (Rehman et al., 2024). Worldwide, 1.57 billion acreages of land are covered by crop production (Fuglie et al., 2024). The agricultural system provides livelihood to 3.83 billion people worldwide (FAO, 2023) and shares 26% of gross domestic production in the economy of many developing countries. In many developing countries, where poverty is widespread and agriculture remains the primary source of livelihood, climate variability poses a critical challenge and has become a
major focus of research and policy attention (Huong et al., 2019; Shahzad et al., 2019a). In agriculture sector, among all the natural events, climate related natural events like floods and droughts have the most devastating effects (Tran et al., 2022). In the last 30 years, an estimated USD 3.8 trillion in crops and livestock production has been lost due to disasters (FAO, 2023a). This shows that on average disasters have caused loss of USD 123 billion annually, which makes up of 5% of the annual global agricultural gross domestic product. According to recent global studies, climate change is a threat to human wellbeing and planetary health, with higher risk and projected adverse impacts due to incremental global warming (IPCC, 2023; Khan et al., 2022; Shahzad, 2025). The intense droughts and high temperature are expected to lower water availability by 19%, further causing up to 30% decrease in food crops by 2050 (Aryal et al., 2020).
Climate change combined with economic shocks has become essential factor that extremely affect the food unavailability issues (Shahzad, 2025). The State of Food Security and Nutrition in the World Report highlights that 2.4 billion people experienced food insecurity of moderate or severe levels while 900 million people faced severe levels of food insecurity (FAO, 2023b). Due to changing climate, the earth's temperature continues to increase, thus increasing the likelihood of disasters like floods, heat and cold waves (Shahzad, 2024). (Hultgren et al. (2025) show that climate change continues to significantly reduce global agricultural productivity, with cross–continental analysis underscores the urgency for crop–level and region–specific adaptation evidence (Rehman et al., 2015). It is recognized that low production of food crops is due to un-accessibility and non- availability of technical and latest knowledge required to adapt to cope with climate variability Dorji and Morrison-Saunders, 2026). In addition, climate variability also affects other sources of livelihoods such as animals raised on farms, rangelands, fruits, vegetables and horticulture. Vulnerability in terms of land productivity is reflected through indicators such as farm intensity, crop yield, soil fertility, and farmland location (Field and Barros, 2014). These factors are further influenced by alterations in weather patterns, including sudden flooding events, temperature fluctuations, and erratic rainfall (Teshome, 2016; Shahzad et al., 2019b).
Across South Asia, rising temperatures, erratic rainfall, and CO2 emissions are undermining food security and agricultural output (Behera, 2024; Shahzad et al., 2024). A recent regional assessment projects a 7.5% drop in agricultural output by 2050 and highlights key adaptation priorities, including irrigation modernization and expanded digital climate–information services for farmers (Gbadegesin et al., 2026).
A recent south Asian climatic review highlighted 60% reduction in wheat yield in Bangladesh; 20%−25% cotton loss in India, and situating Pakistan within a broader region of high agricultural vulnerability (Shahzad et al., 2024; Thakur and Bajagain, 2019). Pakistan is experiencing pronounced impacts of climate change, evidenced by rising trends in hot days and nights, increased frequency of cold nights, and more frequent extreme precipitation events (Nida et al., 2025).
According to Pakistan Economic Survey (2023–24), agriculture sector is crucial to the country's economic growth, food security, generation of employment and particularly poverty reduction at farm level. It contributes 24% to the gross domestic production and provides 37.4% employment to the labor force. Moreover, 65%−70% of the population depends on agriculture for their income source (Pakistan Economic Survey, 2024). Agriculture sector in Pakistan is dominated by the production of cash crops like, rice, wheat, maize and sugarcane, which majorly contributes to the overall performance of agriculture sector. According to Government of Pakistan 2019, the country's average temperature has increased by 0.6 °C during past decade (GOP, 2019). Pakistan has experienced over USD 3.7 billion loss in agriculture and livestock sector due to floods of 2022 (World Bank, 2022). Additionally, due to water shortages, agriculture sector is adversely affected in terms of maintaining growth (0.85%), leading to poverty and food unavailability in farm households (Shahzad et al., 2024; Pakistan Bureau of Statistics, 2025; Riaz et al., 2026; OCHA, 2019). Climate–driven declines in fodder quality and quantity are also undermining livestock productivity, weakening rural livelihoods, and posing a growing threat to national food security (Dogar et al., 2025). National Nutritional Survey reveals that 36.9% of households in Pakistan faced “food insecurity,” with 18.3% experiencing severe levels of it (UNICEF, 2018). In Pakistan, like other developing countries, climate simulation models have forecasted the decline of rice production from 15.2 to 17.2% by the year 2040–2069 (Ahmad et al., 2015).
Understanding food security under climate change requires an examination of household-level vulnerability, as these factors determine the capacity to cope with water stress and agricultural risks (Anwar et al., 2020). Vulnerability assessments in agricultural households are commonly conducted using indicator-based approaches (Tate, 2012; Williams et al., 2018). While widely applied, these methods often fail to capture variations in vulnerability (Vincent, 2007; Bennett et al., 2016). Vulnerability is a multidimensional concept that varies across spatial and temporal scales, influenced by socio-economic, environmental, and geographic factors (Füssel and Klein, 2006). Indicator-based approaches help identify and prioritize vulnerable regions and populations, but challenges remain regarding weighting, data accuracy, and conceptual framing (Birkmann, 2006; Shukla et al., 2021). Frameworks such as the Sustainable Livelihood Approach (SLA), which considers human, physical, financial, natural, and social capital, have been used to develop climate vulnerability indices (Shahzad et al., 2019c; Ullah et al., 2021). In Pakistan, limited studies assess vulnerability and adaptation at the household level, with existing research primarily focusing on climate impacts on agricultural productivity (Ahmad and Ma, 2020; Khan et al., 2021, 2024; Faisal et al., 2021).
This study is novel because it addresses a clear research gap: despite the central role of Punjab's rice–wheat system in national food security, farm–level vulnerability assessments that integrate socioeconomic, livelihood, and climate-related factors remain scarce. By examining how these determinants jointly shape household exposure, sensitivity, and adaptive capacity, the study makes an important contribution to the empirical literature, offering a more nuanced and context-specific understanding of vulnerability than existing, more aggregated analyses. Its findings carry strong significance, as they provide actionable insights into how agricultural households navigate climate variability, thereby informing more targeted, equitable, and effective adaptation interventions in one of Pakistan's most climate-sensitive agricultural regions.
To achieve this, we employ the Intergovernmental Panel on Climate Change (IPCC) framework—encompassing exposure, sensitivity, and adaptive capacity—within a Climate Vulnerability Index (CVI) model, integrated with the Analytical Hierarchy Process (AHP). This combined approach enables a comprehensive assessment of agricultural households' vulnerability and adaptive responses. The primary objectives of this study are, first, to explicitly assess the extent to which changing meteorological patterns—such as shifts in temperature, rainfall, and extreme weather events—have affected the livelihoods, production systems, and overall wellbeing of farm households; and second, to systematically examine the range of adaptation strategies adopted by rice–wheat farmers in Punjab, Pakistan, in response to these evolving climatic conditions.
1.1 Conceptual framework of climate vulnerability assessment and adaptation measures in agricultural households
The conceptual framework for climate vulnerability assessment and adaptation in agricultural households of the current study is grounded in the Intergovernmental Panel on Climate Change (IPCC) approach, which identifies three key components of vulnerability: exposure, sensitivity, and adaptive capacity (IPCC, 2007) (Figure 1, Table 1). Vulnerability, as defined by the IPCC Climate Change (2014), refers to the degree to which a system is unable to cope with the adverse effects of climate variability and extremes. It is shaped by multiple socioeconomic and environmental factors, including fragile economic growth, high dependence on agriculture, rapid population growth, urbanization, weak policy decisions, low literacy rates, gender inequality, water scarcity, land degradation, and declining forest cover (Ahmed et al., 2021; Mekonen and Berlie, 2021; Abid et al., 2016). In this study, vulnerability is assessed across climatic, biological, and biophysical dimensions (Hinkel, 2011). Climatic factors include temperature and rainfall variability; biological factors encompass pest infestations and plant diseases; and biophysical factors involve riverbank erosion, soil fertility loss, floods, and droughts. Sensitivity reflects the extent to which farming systems are affected by climate stimuli, either directly—such as changes in crop yield due to temperature and rainfall variation—or indirectly, through impacts like soil degradation and resource scarcity (IPCC, 2018). Farm-level sensitivity includes lack of equipment, limited land for cereal cultivation, soil fertility loss, seed shortages, and inadequate fertilizer availability, as well as livestock-related constraints. Adaptation, on the other hand, refers to the capacity to minimize yield losses and capitalize on favorable climatic conditions (Adams, 1999), and depends on regional, local, and community-level capacities (Maddison, 2007; Gbetibouo, 2009; Gbetibouo et al., 2010). Adaptive capacity is evaluated through five capital assets: natural (land availability and soil fertility), physical (market access, livestock ownership, irrigation type), financial (credit access, farm insurance, income sources), human (education, farming experience, health access), and social (community support and organizational membership). This integrated framework allows for the inclusion of both qualitative and quantitative indicators across multiple dimensions—climatic, biophysical, social, financial, human, and natural—providing a comprehensive basis for assessing resilience and adaptation strategies (Datta and Behera, 2022; Fahad et al., 2020; Naqvi et al., 2020).

Table 1. Criteria indicators of climate change risks in Tehsil Narowal and Shakargarh.
| Components | Sub components | Variables | Units | Variable description and reported question | References |
|---|---|---|---|---|---|
| Exposure | Climatic exposure | Floods intensity | Percentage | 3 components of climatic exposures like flood intensity, temperature, and rainfall variability affect the livelihood of many farmers. | Ahmed et al., 2021 |
| Change in temperature | Celsius | Mehmood et al., 2021 | |||
| Change in rainfall patterns | Millimeter | Ullah et al., 2021 | |||
| Natural exposure | Ground water depletion | Percentage | The respondents were asked whether the erosion of the riverbank and degradation of ground water quality affect their farms. | Khan et al., 2021 | |
| Erosion of riverbank | Percentage | Shahzad et al., 2019c | |||
| Biological exposure | Severe crop pest | Percentage | This section described how biological hazards like crop and disease attacks affect agricultural activities. | Hahn et al., 2009 | |
| Sensitivity | Livelihood sensitivity | Livestock sensitivity | Percentage | The respondents were asked about their sensitivity regarding livestock ownership, access to veterinary services, and whether they had basic health and water services. | Hahn et al., 2009 |
| Veterinary sensitivity | Percentage | Khan et al., 2021 | |||
| Water sensitivity | Percentage | Shahzad et al., 2021 | |||
| Health sensitivity | Percentage | Shahzad et al., 2021 | |||
| Social sensitivity | Displacement Sensitivity | Percentage | This section described the social sensitivity of the farmers, as most of the respondents lost their farms due to uneven temperature and rainfall patterns and migrates to the urban areas. | ||
| Farm sensitivity | Soil fertility loss | Percentage | This section described the 2 components of farm sensitivity, as loss of soil fertility and unavailability of fertilizers affected the farm productivity. The unavailability of fertilizers affects the farm yield negatively. | ||
| Unavailability of fertilizers | Percentage | Young et al., 2009 | |||
| Adaptive capacity | Natural capital (NC) | Land ownership rights | Percentage | Greater land ownership helps to increase adaptive capacity due to its correlation with the use of natural resources. Better soil quality and land size help farmers grow different crops and overall it increases the financial capital of the region. | |
| Soil quality | Percentage | Shahzad et al., 2021 | |||
| Land size | Percentage | ||||
| Physical capital (PC) | Type of irrigation | Percentage | Different types of irrigation will help the farmers swap between them according to climatic changes. Livestock ownership increases the adaptive capacity of households as they can be sold for income. Access to the main road helps to increase the farmer's adaptive capacity | ||
| Livestock ownership | Percentage | Naqvi et al., 2020 | |||
| Main road accessibility | Percentage | Naqvi et al., 2020 | |||
| Financial capital (FC) | Household income sources | Percentage | A greater number of income sources increases the overall financial capital of households. Farm insurance and access to loans help farmers to increase their adaptive capacity as it helps farmers when they do not have enough money to grow crops | Fahad et al., 2020 | |
| Farm Insurance | Percentage | Duffy et al., 2021 | |||
| Access to loans | Percentage | Hahn et al., 2009 | |||
| Human capital (HC) | Farming experience | Percentage | Greater farming experience and farmer education increase the adaptive capacity as they know very well how to tackle difficult situations during different hazards. Less distance increases the farmer's adaptive capacity as it helps farmers to maintain their health without wasting time and resources on commute. | ||
| Farmer education | Percentage | Shahzad et al., 2021 | |||
| Distance to the basic health facility | Percentage | Fahad et al., 2020 | |||
| Social capital (SC) | Institutional support provided during climatic hazards | Percentage | The support provided by the local institution at right time during hazards enhances the adaptive capacity of rural farmers. In addition, attending any training or educational program helps to improve their adaptive capacity toward climate change. A greater number of household members increase the adaptive capacity as they help in farms to increase yield. | ||
| Joined any community organization | Percentage | Naqvi et al., 2020 | |||
| Household members | Percentage |
2 Methodology
This study assesses farmers' perceptions of vulnerability and adaptation at the household level using a combination of statistical techniques and Geographic Information System (GIS). The methodology includes: (i) study area and site selection, (ii) identification of risk perception and adaptation indicators, (iii) data collection, (iv) construction of vulnerability and adaptation indices, and (v) analysis of indicator impacts.
2.1 Study area and site selection
Punjab is Pakistan's leading agricultural province, accounting for nearly half of the country's arable land and producing over 70% of its grain crops (GOP, 2019). This study focuses on District Narowal (32°05′58″N, 74°52′29″E) in Gujranwala Division, a major rice–wheat agro-ecological zone located on fertile alluvial plains with an average elevation of 236 meters (Figure 2). The district covers 4,623 sq. km and comprises three tehsils—Narowal, Shakargarh, and Zafarwal—with a population of approximately 1.7 million, of which 85% resides in rural areas and depends largely on agriculture, livestock, and related activities for livelihoods (Pakistan Bureau of Statistics, 2025). The region experiences a humid subtropical climate, with mean temperatures rising in recent decades, summer peaks exceeding 42 °C, winter averages around 17 °C, and annual rainfall of about 900 mm, more than 70% of which occurs during the monsoon season, often in intense and erratic patterns (GOP (Government of Pakistan), 2020). Narowal was selected due to its classification as highly vulnerable under the National Monsoon Contingency Plan (National Disaster Management Authority (NDMA), 2018), having experienced major flood events in 2010, 2012, and 2022 that inundated over 65 villages and large areas of cropland, particularly affecting rice and livestock systems. In addition to flooding, the district faces droughts, pest outbreaks, waterlogging, and soil degradation, creating a multi-hazard risk environment (Mehmood, 2022). Within the district, Shakargarh and Narowal tehsils were purposively selected to capture contrasting vulnerability conditions: Shakargarh is highly exposed to riverine flooding due to its proximity to the Ravi River, while Narowal faces comparatively lower flood intensity but greater challenges from rainfall variability and resource constraints. This combination of climatic variability, agro-ecological sensitivity, and socio-economic dependence on agriculture makes Narowal an appropriate case for analyzing multi-hazard climate vulnerability and adaptation strategies.

2.2 Data sources and collection tools
A mixed-method approach combining quantitative and qualitative techniques was employed. Primary data were collected through structured questionnaires, focus group discussions (FGDs), and field observations, while secondary data on temperature and precipitation trends (2002–2021) were obtained from the Pakistan Meteorological Department (PMD).
The study population comprised rice–wheat farming households located in flood-prone riverbank areas of District Narowal. A total sample of 300 farm households was selected using a 95% confidence level and 6% margin of error (Fahad et al., 2018). A multistage purposive sampling approach was applied: two tehsils were selected, followed by three union councils per tehsil and four villages per council, prioritizing locations proximate to riverbanks. From each village, 10–12 farmers engaged in rice–wheat farming were interviewed. Data were collected between December 2022 and March 2023 in Punjabi and Urdu. Farmers' perceptions of vulnerability and adaptive capacity were assessed using a Likert scale (none, low, moderate, high).
2.3 Climate change vulnerability and index measurement
Farmers' vulnerability perceptions were assessed using the Climate Change Risk Perception Score (CCRPS). Respondents rated 13 climate-related risks on a four-point Likert scale: 0 (none), 1 (low), 2 (moderate), and 3 (high). The CCRPS was calculated as a weighted aggregate:
CCRPS=(0×CCRPn)+(1×CCRPl)+(2×CCRPm)+(3×CCRPh)
where CCRP*n*, CCRP*l*, CCRP*m*, and CCRP*h*represent the number of respondents reporting no, low, moderate, and high perception levels, respectively.
To ensure comparability across risks, the CCRPS was standardized into the Standardized Climate Change Risk Perception Index (SCCRPI), scaled from 0 to 100:
SCCRPI=(CCRPSMaximum Possible CCRPS)×100
where the maximum possible CCRPS equals the total number of respondents multiplied by the highest Likert score (3). The SCCRPI was calculated separately for each risk and then analyzed collectively to assess overall exposure patterns.
2.4 Assessment of farmers' adaptation strategies
Adaptation was evaluated using the Sustainable Livelihood Framework (SLF) across five capital assets: natural, physical, financial, human, and social. Indicator weights were derived using the Analytic Hierarchy Process (AHP) based on expert judgement (Saaty, 2008). Pairwise comparison matrices were normalized, and consistency was verified using the Consistency Ratio (CR), with CR < 0.10 considered acceptable.
The Adaptive Capacity Index (ACI) was calculated as a weighted sum of indicators:
ACI=∑(wi×Xi)
where wi represents the AHP-derived weight and Xi the corresponding indicator value. All indicators were normalized using the min–max method to a 0–1 scale prior to aggregation.
2.5 Climate Vulnerability Index (IPCC-CVI)
The study applied the IPCC vulnerability framework, where vulnerability is a function of exposure (E), sensitivity (S), and adaptive capacity (AC):
CVI=(E-AC)×S
Also shown in detail components in Table 1. All three components were normalized to a 0–1 scale to ensure comparability. Exposure and sensitivity indices were derived from SCCRPI and related indicators, while adaptive capacity was calculated using the ACI. Equal weighting was applied to the three components to avoid bias, following standard IPCC-based approaches.
The final CVI ranged from −1 (least vulnerable) to +1 (most vulnerable), enabling comparative assessment across tehsils.
2.6 Climate data analysis and change detection
Historical climate data, including temperature and rainfall records for the period 2002–2021, were obtained from the Pakistan Meteorological Department (PMD) to assess trends and variability. Simple Linear Regression (SLR) was applied to detect changes in mean minimum and maximum temperatures, as well as rainfall patterns over time.
Rainfall data were further classified into seasonal categories, distinguishing between winter and summer periods. The summer season was subdivided into pre-monsoon, mid-monsoon, and post-monsoon phases to capture intra-seasonal variability.
Statistical analyses, including descriptive statistics and index calculations, were conducted using SPSS (Version 21). Spatial distribution maps of the study area were developed using ArcGIS to support the visualization of climatic patterns and trends.
3 Results
3.1 Demographic and climatic trends in District Narowal
Most respondents in both tehsils were middle-aged (36–50 years), followed by young (36.4 and 29.4%) and older groups (30.7 and 28.8%). Household heads were predominantly male (90.7 and 82.5%), with only 9.3 and 17.5% female-headed households. Education levels were low, as 37.1 and 43.8% had primary education, while 29.3% were illiterate. Average family size was six, with most households having 6–10 members. Over half of respondents (53.6 and 40.6%) owned medium-sized farms (2.5–6.5 acres), and most belonged to medium-income groups (≈ PKR 51,000/month). Farming experience averaged 47–51 years, with 41%−44% having 21–30 years of experience. Cemented houses were common (62.1 and 53.1%).
Farmers' perceptions of climate variability were validated using 20 years of meteorological data (2002–2021) on temperature and rainfall (Figure 3). Simple Linear Regression (SLR) revealed a declining trend in mean maximum annual temperature (y = −0.03 + 29.88) and an increasing trend in mean minimum temperature (y = 0.08 + 16.67). The highest mean maximum temperature (30.42 °C) occurred in 2009, while the highest mean minimum temperature (26.03 °C) was recorded in 2018. Rainfall patterns also showed significant variability, with an increasing trend (y = 0.91 + 71.39). Extreme rainfall events were observed in 2003 (133 mm), 2013 (116 mm), and 2014 (110 mm), contributing to severe floods that damaged crops and livestock. These findings align with farmers' reports of erratic rainfall and rising temperatures, particularly during the Kharif season (May–October), when monthly maximum temperatures ranged between 31 °C and 38 °C.

3.2 Climate vulnerability analysis and farmers' response
This section presents the results of climate vulnerability analysis based on farmers' responses. First, the Standardized Climate Change Risk Perception Index (SCCRPI) was calculated along with exposure and sensitivity components, which were ranked to reflect farmers' perceptions of climatic events. Second, farmers' adaptive capacity was assessed and integrated into the IPCC-CVI framework to determine overall vulnerability. Tehsil Narowal exhibited higher exposure (index = 0.83) compared to Shakargarh (index = 0.62) as shown in Tables 2 and 3 respectively. CCRPS values ranged from 269–351 in Narowal and 210–326 in Shakargarh, indicating moderate to high exposure in Narowal and moderate to low exposure in Shakargarh. Flood intensity ranked highest in both tehsils (SCCRPI: Narowal = 95.48; Shakargarh = 94.01), followed by rainfall and temperature fluctuations. SCCRPI values ranged from 95.48 to 73.17 in Narowal and 94.01 to 60.74 in Shakargarh, showing heterogeneous exposure patterns. Both tehsils demonstrated moderate to high sensitivity, with indices of 0.76 (Narowal) and 0.71 (Shakargarh). Sensitivity was influenced by climatic variability and associated biological hazards, affecting farm assets such as fertilizer availability, soil fertility, and access to water and health services. SCCRPI values for sensitivity ranged from 83.47 to 62.19 in Narowal and 85.33 to 70 in Shakargarh.
Table 2. Results of CCRPS and SCCRPI of exposure and sensitivity in tehsil Narowal.
| Sr. No. | Climatic events | No | Low | Medium | High | CCRPS | SCCRPI | Rank |
|---|---|---|---|---|---|---|---|---|
| 1 | Flood intensity | 0 | 17 | 34 | 89 | 351 | 95.48 | 1 |
| 2 | Ground water depletion | 5 | 37 | 55 | 43 | 276 | 75.08 | 9 |
| 3 | Erosion of river bank | 5 | 29 | 48 | 58 | 299 | 81.38 | 6 |
| 4 | Change in temperature | 4 | 31 | 27 | 78 | 319 | 86.79 | 4 |
| 5 | Change in rainfall patterns | 7 | 18 | 23 | 92 | 340 | 92.49 | 2 |
| 6 | Severe crop pest | 17 | 31 | 38 | 54 | 269 | 73.17 | 10 |
| 7 | Soil fertility loss | 12 | 25 | 52 | 51 | 282 | 62.19 | 14 |
| 8 | Livestock ownership | 4 | 50 | 38 | 48 | 270 | 73.49 | 11 |
| 9 | Access to veterinary services | 16 | 37 | 48 | 39 | 250 | 68.08 | 13 |
| 10 | Basic health services | 15 | 94 | 15 | 16 | 172 | 76.13 | 7 |
| 11 | Drinking water scarcity | 10 | 29 | 58 | 43 | 274 | 74.53 | 8 |
| 12 | Migrate due to financial crises | 11 | 14 | 53 | 62 | 306 | 82.24 | 5 |
| 13 | Outbreak of animal diseases | 19 | 32 | 31 | 58 | 268 | 72.05 | 12 |
| 14 | Unavailability of fertilizers | 4 | 16 | 31 | 89 | 329 | 83.47 | 3 |
Table 3. Results of CCRPS and SCCRPI in tehsil Shakargarh.
| Sr. No. | Climatic events | No | Low | Medium | High | CCRPS | SCCRPI | Rank |
|---|---|---|---|---|---|---|---|---|
| 1 | Flood intensity | 7 | 102 | 22 | 29 | 233 | 94.01 | 1 |
| 2 | Ground water depletion | 8 | 80 | 59 | 13 | 237 | 68.56 | 12 |
| 3 | Erosion of river bank | 21 | 89 | 29 | 21 | 210 | 60.74 | 14 |
| 4 | Change in temperature | 6 | 59 | 38 | 57 | 306 | 88.16 | 3 |
| 5 | Change in rainfall patterns | 13 | 37 | 41 | 69 | 326 | 82.52 | 6 |
| 6 | Severe crop pest | 10 | 62 | 42 | 46 | 284 | 67.99 | 13 |
| 7 | Soil fertility loss | 11 | 75 | 55 | 19 | 242 | 70 | 11 |
| 8 | Livestock ownership | 8 | 53 | 56 | 43 | 294 | 85.04 | 5 |
| 9 | Access to veterinary services | 17 | 64 | 57 | 22 | 244 | 70.81 | 10 |
| 10 | Basic health services | 5 | 74 | 56 | 25 | 261 | 75.89 | 8 |
| 11 | Drinking water scarcity | 7 | 78 | 51 | 24 | 252 | 72.95 | 9 |
| 12 | Migrate due to financial crises | 20 | 30 | 67 | 43 | 293 | 84.55 | 4 |
| 13 | Outbreak of animal diseases | 4 | 77 | 39 | 40 | 275 | 79.48 | 7 |
| 14 | Unavailability of fertilizers | 8 | 44 | 43 | 65 | 325 | 85.33 | 2 |
3.3 Adaptive capacity assessment using AHP
Farmers' adaptive capacity was evaluated using the Analytical Hierarchy Process (AHP) to assign weights to indicators under five capital assets: natural, physical, financial, human, and social. Pairwise comparisons produced a Consistency Ratio (CR) of 0.09 for Narowal and 0.05 for Shakargarh, indicating acceptable consistency in judgments (Tables 4, 5). Under natural capital, land ownership had the highest weight (0.13 in Narowal; 0.15 in Shakargarh). For physical capital, irrigation sources were critical (0.06 and 0.07), reflecting reliance on natural water systems. Financial capital indicators such as farm insurance scored 0.10 and 0.08, highlighting limited financial resilience. Social capital emphasized support during hazards (0.10 and 0.09), pointing to gaps in institutional assistance (Tables 6).
Table 4. Normalized comparison matrix and criteria weight of Tehsil Narowal.
| Indicators | NC1 | NC2 | NC3 | PC1 | PC2 | PC3 | FC1 | FC2 | FC3 | HC1 | HC2 | HC3 | SC1 | SC2 | SC3 | Criteria weight |
|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|
| NC1 | 0.05 | 0.1 | 0.1 | 0.02 | 0.14 | 0.05 | 0.07 | 0.03 | 0.02 | 0.02 | 0.1 | 0.04 | 0.04 | 0.07 | 0.08 | 0.13 |
| NC2 | 0.03 | 0.05 | 0.02 | 0.09 | 0.03 | 0.05 | 0.11 | 0.05 | 0.08 | 0.08 | 0.07 | 0.04 | 0.02 | 0.07 | 0.08 | 0.06 |
| NC3 | 0.02 | 0.1 | 0.05 | 0.09 | 0.03 | 0.05 | 0.07 | 0.05 | 0.08 | 0.02 | 0.07 | 0.04 | 0.02 | 0.07 | 0.08 | 0.05 |
| PC1 | 0.1 | 0.02 | 0.02 | 0.04 | 0.13 | 0.05 | 0.07 | 0.05 | 0.02 | 0.08 | 0.07 | 0.02 | 0.02 | 0.07 | 0.13 | 0.06 |
| PC2 | 0.02 | 0.1 | 0.1 | 0.02 | 0.07 | 0.1 | 0.07 | 0.05 | 0.08 | 0.12 | 0.07 | 0.16 | 0.04 | 0.07 | 0.08 | 0.07 |
| PC3 | 0.1 | 0.1 | 0.1 | 0.09 | 0.07 | 0.1 | 0.07 | 0.21 | 0.08 | 0.08 | 0.02 | 0.16 | 0.08 | 0.07 | 0.08 | 0.09 |
| FC1 | 0.02 | 0.01 | 0.02 | 0.02 | 0.03 | 0.05 | 0.04 | 0.05 | 0.08 | 0.08 | 0.07 | 0.02 | 0.02 | 0.02 | 0.02 | 0.03 |
| FC2 | 0.15 | 0.1 | 0.1 | 0.09 | 0.14 | 0.05 | 0.07 | 0.10 | 0.08 | 0.08 | 0.07 | 0.16 | 0.17 | 0.07 | 0.08 | 0.10 |
| FC3 | 0.1 | 0.02 | 0.1 | 0.09 | 0.03 | 0.05 | 0.02 | 0.05 | 0.04 | 0.08 | 0.07 | 0.04 | 0.02 | 0.07 | 0.02 | 0.05 |
| HC1 | 0.1 | 0.02 | 0.02 | 0.02 | 0.02 | 0.05 | 0.02 | 0.05 | 0.02 | 0.04 | 0.07 | 0.04 | 0.02 | 0.07 | 0.02 | 0.04 |
| HC2 | 0.01 | 0.02 | 0.1 | 0.02 | 0.03 | 0.19 | 0.02 | 0.05 | 0.02 | 0.02 | 0.03 | 0.04 | 0.04 | 0.02 | 0.01 | 0.03 |
| HC3 | 0.1 | 0.1 | 0.1 | 0.14 | 0.03 | 0.05 | 0.11 | 0.05 | 0.08 | 0.08 | 0.07 | 0.08 | 0.25 | 0.07 | 0.08 | 0.09 |
| SC1 | 0.1 | 0.15 | 0.15 | 0.19 | 0.14 | 0.1 | 0.11 | 0.05 | 0.13 | 0.12 | 0.07 | 0.02 | 0.08 | 0.11 | 0.04 | 0.10 |
| SC2 | 0.02 | 0.02 | 0.02 | 0.02 | 0.03 | 0.05 | 0.07 | 0.05 | 0.02 | 0.02 | 0.07 | 0.04 | 0.02 | 0.03 | 0.08 | 0.04 |
| SC3 | 0.02 | 0.02 | 0.02 | 0.01 | 0.03 | 0.05 | 0.07 | 0.05 | 0.08 | 0.08 | 0.1 | 0.04 | 0.08 | 0.02 | 0.04 | 0.05 |
| λmaxax = 17.17484127 CI= 0.15 RI= 1.59 CR=CI/RI C.R= 0.09 |
NC1-NC3, Natural capital; PC1-PC3, physical capital; FC1-FC3, financial capital; HC1-HC3, human capital; SC1-SC3, social capital of Tehsil Narowal.
Table 5. Normalized comparison matrix and criteria weight of tehsil Shakargarh.
| Indicators | NC1 | NC2 | NC3 | PC1 | PC2 | PC3 | FC1 | FC2 | FC3 | HC1 | HC2 | HC3 | SC1 | SC2 | SC3 | Criteria weight |
|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|
| NC1 | 0.04 | 0.1 | 0.03 | 0.02 | 0.1 | 0.02 | 0.08 | 0.02 | 0.02 | 0.07 | 0.07 | 0.04 | 0.04 | 0.07 | 0.08 | 0.15 |
| NC2 | 0.02 | 0.05 | 0.03 | 0.11 | 0.02 | 0.04 | 0.02 | 0.03 | 0.13 | 0.07 | 0.07 | 0.04 | 0.16 | 0.07 | 0.08 | 0.06 |
| NC3 | 0.09 | 0.1 | 0.07 | 0.11 | 0.1 | 0.17 | 0.08 | 0.14 | 0.03 | 0.07 | 0.07 | 0.04 | 0.33 | 0.07 | 0.08 | 0.08 |
| PC1 | 0.09 | 0.5 | 0.03 | 0.05 | 0.1 | 0.04 | 0.08 | 0.03 | 0.13 | 0.07 | 0.07 | 0.02 | 0.16 | 0.07 | 0.02 | 0.07 |
| PC2 | 0.02 | 0.1 | 0.03 | 0.02 | 0.05 | 0.17 | 0.08 | 0.03 | 0.03 | 0.07 | 0.07 | 0.04 | 0.04 | 0.01 | 0.08 | 0.06 |
| PC3 | 0.13 | 0.1 | 0.03 | 0.11 | 0.02 | 0.08 | 0.13 | 0.14 | 0.13 | 0.07 | 0.07 | 0.08 | 0.16 | 0.07 | 0.08 | 0.09 |
| FC1 | 0.02 | 0.1 | 0.03 | 0.02 | 0.02 | 0.02 | 0.04 | 0.5 | 0.03 | 0.07 | 0.07 | 0.04 | 0.04 | 0.07 | 0.08 | 0.05 |
| FC2 | 0.13 | 0.1 | 0.03 | 0.11 | 0.1 | 0.04 | 0.08 | 0.07 | 0.13 | 0.07 | 0.07 | 0.04 | 0.04 | 0.07 | 0.08 | 0.08 |
| FC3 | 0.13 | 0.02 | 0.14 | 0.02 | 0.1 | 0.04 | 0.08 | 0.03 | 0.06 | 0.07 | 0.07 | 0.17 | 0.04 | 0.07 | 0.08 | 0.08 |
| HC1 | 0.02 | 0.02 | 0.03 | 0.02 | 0.02 | 0.04 | 0.02 | 0.03 | 0.03 | 0.03 | 0.07 | 0.04 | 0.04 | 0.07 | 0.02 | 0.03 |
| HC2 | 0.02 | 0.02 | 0.03 | 0.02 | 0.02 | 0.04 | 0.02 | 0.03 | 0.03 | 0.01 | 0.03 | 0.04 | 0.04 | 0.07 | 0.08 | 0.03 |
| HC3 | 0.09 | 0.1 | 0.14 | 0.16 | 0.1 | 0.08 | 0.08 | 0.14 | 0.03 | 0.07 | 0.07 | 0.08 | 0.04 | 0.07 | 0.08 | 0.09 |
| SC1 | 0.09 | 0.02 | 0.22 | 0.02 | 0.1 | 0.04 | 0.08 | 0.14 | 0.13 | 0.07 | 0.07 | 0.17 | 0.08 | 0.07 | 0.08 | 0.09 |
| SC2 | 0.02 | 0.02 | 0.03 | 0.02 | 0.1 | 0.04 | 0.02 | 0.03 | 0.03 | 0.01 | 0.01 | 0.04 | 0.04 | 0.03 | 0.02 | 0.03 |
| SC3 | 0.02 | 0.02 | 0.03 | 0.11 | 0.02 | 0.04 | 0.02 | 0.03 | 0.03 | 0.07 | 0.01 | 0.04 | 0.04 | 0.07 | 0.04 | 0.04 |
| λax = 16.46 C.I= 0.1 RI= 1.59 CR= CI/RI= 0.05. |
NC1-NC3, natural capital; PC1-PC3, physical capital; FC1-FC3, financial capital; HC1-HC3, human capital; SC1-SC3, social capital of Tehsil Shakargarh.
Table 6. Farmer's response regarding adaptive capacity components.
| Major components | Sub-components | Categories | Percentage of respondents | Rank | ||
|---|---|---|---|---|---|---|
| Narowal | Shakargarh | Narowal | Shakargarh | |||
| Natural capital | Land ownership | Land with ownership rights | 55.7 | 54.4 | 3 | 3 |
| Land without ownership rights | 44.3 | 45.6 | 1 | 1 | ||
| Soil quality | Fertile | 45.7 | 48.1 | 3 | 3 | |
| Partially infertile | 30 | 30.6 | 2 | 2 | ||
| Fully infertile | 24.3 | 21.3 | 1 | 1 | ||
| Land size | 1–3 hectare | 33.6 | 40.6 | 2 | 3 | |
| 4–5 hectare | 53.6 | 49.4 | 3 | 2 | ||
| More than 5 hectare | 12.9 | 10 | 1 | 1 | ||
| Physical capital | Type of irrigation | No irrigation | 36.4 | 20.6 | 1 | 2 |
| Canal water | 45 | 39.4 | 3 | 1 | ||
| Ground water | 18.6 | 40 | 2 | 3 | ||
| Livestock ownership | No livestock | 35.7 | 31.9 | 3 | 2 | |
| 1–5 livestock | 37.9 | 44.4 | 1 | 1 | ||
| 6–10 livestock | 26.4 | 23.8 | 2 | 3 | ||
| Main road accessibility | Not accessible | 17.9 | 15 | 3 | 2 | |
| partially accessible | 46.4 | 36.3 | 2 | 3 | ||
| fully accessible | 35.7 | 48.8 | 1 | 1 | ||
| Financial capital | Household income sources | Single source | 35 | 33.1 | 3 | 3 |
| Two sources | 44.3 | 43.8 | 2 | 2 | ||
| more than 2 sources | 20.7 | 23.1 | 1 | 1 | ||
| Farm insurance | Yes | 12.1 | 8.8 | 1 | 3 | |
| No | 87.9 | 91.3 | 3 | 1 | ||
| Access to loans | Yes | 52.1 | 51.9 | 1 | 1 | |
| No | 47.9 | 48.1 | 3 | 3 | ||
| Human capital | Farming experience | 1–10 years | 24.3 | 18.8 | 2 | 3 |
| 11–15 years | 34.3 | 43.8 | 3 | 2 | ||
| more than 15 years | 41.4 | 37.5 | 1 | 1 | ||
| Farmer education | Illiterate | 29.3 | 29.4 | 2 | 1 | |
| Distance to basic health facility | 6km | 37.1 | 43.8 | 3 | 2 | |
| 7–10km | 33.6 | 26.9 | 1 | 3 | ||
| More than 10km | 51.4 | 40.6 | 3 | 3 | ||
| Social capital | Support during hazards from local institutions | Yes | 25 | 28.1 | 2 | 2 |
| No | 23.6 | 31.3 | 1 | 1 | ||
| Farmers joined community organization | Yes | 47.9 | 46.9 | 1 | 3 | |
| No | 52.1 | 53.1 | 3 | 1 | ||
| Adaptive capacity indexa | 1= low | 0.48 | 0.36 | |||
| 2 = medium | ||||||
| 3 =high |
aACI is calculated from equation (6).
Overall, 28.34% of households exhibited high adaptive capacity, 39.1% moderate, and 32.56% low. Specifically, 55.02% of households reported high adaptive capacity linked to land ownership, soil quality, and irrigation access, while 25 and 19.8% fell into moderate and low categories, respectively.
The considerable variation in the levels of adaptive capacity was found between both tehsils, the farmers of both tehsils had high moderate and low adaptive capacity levels. Households are categorized by their mean adaptive capacity in Table 7.
Table 7. Calculation of IPCC-CVI for Tehsil Narowal and Shakargarh.
| Index dimensions of IPCC-CVI status of vulnerability | Tehsil Narowal | Tehsil Shakargarh |
|---|---|---|
| Adaptive Capacity > exposure= least vulnerable | (0.83-0.48) x 0.76= 0.26 | (0.62-0.36) x 0.71= 0.19 |
| (Adaptive Capacity = Exposure) = Moderately Vulnerable | ||
| (Exposure > adaptive capacity) = extremely vulnerable |
*IPCC-CVI = −1(least vulnerable) to +1 (extreme vulnerable).
After calculating the major components of climate vulnerability for agricultural households in both tehsils, a comparative analysis was conducted (Figure 4). The IPCC-CVI ranges from −1 (least vulnerable) to +1 (most vulnerable). Based on the IPCC-CVI formula, Tehsil Narowal recorded a vulnerability score of 0.26, while Shakargarh scored 0.19, indicating that Narowal is more vulnerable to climatic hazards due to lower adaptive capacity (Table 7). The vulnerability triangle (Figure 5) illustrates that Narowal had higher exposure (0.83) compared to Shakargarh (0.62). Both tehsils showed moderate sensitivity to climatic and natural hazards, with values of 0.76 and 0.72, respectively. Adaptive capacity was relatively low in both areas, with Narowal at 0.48 and Shakargarh at 0.36, reinforcing the overall vulnerability pattern.


3.4 Differential adaptation strategies preferences by the farmers
Farmers adopted several strategies to cope with climatic hazards (Figure 6). The survey revealed that 40.0% of farmers in Narowal and 48.8% in Shakargarh moderately changed their sowing calendars, while 52.8 and 41.9% had high capacity to adjust crop calendars. Fertilizer and pesticide adjustments were common, with 43.6 and 50.6% of respondents reporting moderate changes to maintain soil fertility. However, most farmers did not adopt climate-tolerant crop varieties, showing low adaptive capacity (48.6 and 50.8%). Switching from traditional cereals to cash crops was moderately practiced by 50.0 and 42.5% of farmers to manage financial risks. Tree planting was limited in Narowal (40.0% low capacity) but moderately practiced in Shakargarh (43.8%). A significant proportion of farmers (50.7 and 43.8%) left land fallow for one or more cycles to restore soil organic matter and moisture.

4 Discussion
The study assessed the vulnerability of agricultural households in District Narowal, Punjab, to multi-hazard climate variability using an extended analytical framework. The findings indicate high exposure to climatic hazards, while adaptive capacity remains limited due to inadequate knowledge, financial constraints, and weak institutional linkages. Socio-economic and demographic factors—particularly age, farming experience, landholding size, and ownership rights—emerged as key determinants of vulnerability and adaptation, consistent with previous studies (Abid et al., 2016; Ahmed et al., 2021; Bryan et al., 2013; Mazhar et al., 2021; Yohannes, 2016).
Age-related vulnerability was evident, as older farmers tend to rely on traditional practices and resource-intensive irrigation, increasing environmental stress. Gender dynamics further shape vulnerability: while farming is predominantly male-led, women contribute through off-farm activities such as handicrafts and home gardening (Zamasiya et al., 2017). However, mechanization and sociocultural barriers limit women's participation in modern agricultural practices, corroborating findings by Anjum and Fraser (2021). Education and financial capacity strongly influence adaptive responses. Low education levels constrain awareness of climate risks and preparedness for extreme events, as also observed by Tripathi and Mishra (2017).
Farmers' perceptions of climate change align with meteorological trends, particularly increasing temperature variability and erratic rainfall, which disrupt rice–wheat cropping cycles and reduce productivity (Gbetibouo, 2009). These findings mirror evidence from other regions, including Nigeria, Ethiopia, and Ghana, where similar climatic stresses adversely affect agricultural output (Ayanlade et al., 2017; Araro et al., 2020; Derbile et al., 2022; Ketema and Negeso, 2020).
Flooding emerged as the most significant hazard, causing recurrent crop and livestock losses, especially during post-monsoon periods (Ullah et al., 2018). In addition to direct climatic impacts, biological hazards—including pest outbreaks, invasive weeds, and livestock diseases—significantly exacerbate vulnerability. Farmers reported increased pest attacks, locust swarms, and persistent weeds despite pesticide use, consistent with findings across South Asia and Africa (Khan et al., 2020; Shukla et al., 2021; Dhungana et al., 2020). These pressures have intensified chemical use, degrading soil and water quality while increasing production costs and health risks.
Livestock—which plays a central role in rural livelihoods—has been adversely affected by disease outbreaks such as lumpy skin disease, reducing income and food security (Khan et al., 2021). Financial vulnerability is compounded by the absence of crop insurance, limited access to formal credit, and reliance on exploitative informal lending systems. Land tenure insecurity further constrains adaptation, as smallholders without ownership rights have limited access to government support and reduced incentives to invest in long-term resilience measures.
Social and institutional factors also influence adaptive capacity. Weak extension services, limited infrastructure, and insufficient early warning systems increase sensitivity to climate shocks. Conversely, community networks and social capital enhance resilience by facilitating information sharing and collective coping strategies, as demonstrated in other contexts (Bryan et al., 2013; Tran et al., 2022). Emerging local adaptation practices—such as adjusting sowing calendars, greenhouse cultivation, and short-duration crops—indicate adaptive potential but remain insufficiently scaled.
Comparative evidence from other regions reinforces the patterns observed in this study. For instance, (Mramba and Mapunda 2025) report that smallholder farmers in Tanzania adopt strategies such as crop diversification, early sowing, and the use of drought-tolerant varieties—demonstrating adaptive behaviors that closely parallel those observed among rice–wheat farmers in Punjab. These similarities highlight the broader relevance of localized adaptation strategies across diverse agro-ecological contexts.
Similarly, studies from Ethiopia identify key determinants of adaptation, including crop diversification, adjustments in planting dates, and livestock management practices. These findings provide important comparative insights into how smallholder farmers make adaptation decisions under climate stress, further supporting the evidence from Narowal that adaptive capacity is shaped by both environmental pressures and socio-economic conditions (Shukla et al., 2016; Paul et al., 2019; Parker et al., 2019; Thorlakson and Neufeldt, 2020; Jamshidi et al., 2020).
5 Conclusion and policy implications
This study demonstrates that rice–wheat farming systems in Narowal are subject to significant vulnerability arising from the combined effects of climatic variability, biological hazards, and socio-economic constraints. Climatic stressors, including floods, erratic rainfall, and temperature extremes, disrupt cropping patterns, reduce yields, and increase production costs. At the same time, limited financial capacity, inadequate institutional support, and uneven access to resources constrain farmers' ability to adapt effectively, with small and marginal farmers facing disproportionate risks.
Addressing these vulnerabilities requires a comprehensive and integrated policy response (Arbaiah et al., 2021). Expanding access to finance through crop insurance schemes and affordable credit can enhance farmers' adaptive capacity and reduce reliance on informal lending. At the same time, linking the economic trade-offs inherent in rice–wheat systems to actionable responses such as crop diversification, improved irrigation efficiency, climate-smart agriculture, and land restoration approaches, including Land Degradation Neutrality (LDN)—is critical to enhancing both sustainability and resilience. Such measures can help reduce pressure on water resources while maintaining agricultural productivity.
Investments in infrastructure, including efficient irrigation systems, improved market access, and reliable early warning services, are essential to strengthening system-wide resilience. In parallel, targeted capacity-building initiatives can equip farmers with the knowledge and skills required for climate adaptation and sustainable resource management. Ensuring gender-inclusive approaches is also crucial, as empowering women can enhance household resilience and broaden participation in climate-smart practices.
Integrating these interventions into provincial planning frameworks, while aligning them with national climate and agricultural policies, will be essential for positioning Narowal's farming systems within a broader resilience framework that balances productivity with long-term environmental sustainability under increasing water stress and climate uncertainty.
Generative AI statement
The author(s) declared that Generative AI was not used in the creation of this manuscript.
Any alternative text (alt text) provided alongside figures in this article has been generated by Frontiers with the support of artificial intelligence and reasonable efforts have been made to ensure accuracy, including review by the authors wherever possible. If you identify any issues, please contact us.
Publisher's note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Footnotes
Edited by: Naghmana Ghafoor, Lahore College for Women University, Pakistan
Reviewed by: Sadhon Chandra Swarnokar, Khulna University, Bangladesh
Zainab Nuhu, Bayero University Kano, Nigeria
Statements
- Conflict of interest
The author(s) declared that this work was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
- Funding
The author(s) declared that financial support was not received for this work and/or its publication.
- Data availability
The original contributions presented in the study are included in the article/supplementary material, further inquiries can be directed to the corresponding author.
- Ethics approval
The studies involving humans were approved by Government College University Ethics Board review Committee. The studies were conducted in accordance with the local legislation and institutional requirements. The ethics committee/institutional review board waived the requirement of written informed consent for participation from the participants or the participants' legal guardians/next of kin because Written informed consent was waived because no personally identifiable data were collected from participants, and they were informed that all data would be anonymized and handled in accordance with established ethical guidelines.
- Author contributions
LS: Writing – original draft, Writing – review & editing, Conceptualization, Supervision. MS: Formal analysis, Writing – original draft, Methodology. FS: Validation, Writing – review & editing, Investigation, Resources, Visualization. AS: Methodology, Software, Writing – review & editing. MH: Writing – review & editing, Formal analysis, Visualization, Validation. MF: Investigation, Software, Writing – review & editing. MD: Validation, Methodology, Writing – original draft, Software.
References
- Abid M., Schilling J., Scheffran J., Zulfiqar F. (2016). Climate change vulnerability, adaptation and risk perceptions at farm level in Punjab. Pakistan. Sci. Total Environ. 547, 447 – 460. doi: doi:10.1016/j.scitotenv.1125
- Adams R. A. (1999). On the search for the correct economic assessment method. Clim. Change 41, 363 – 370. doi: doi:10.1023/A:1005434215112
- Ahmad A., Ashfaq M., Rasul G., Wajid S. A., Khaliq T., Rasul F.. (2015). “Impact of climate change on the rice–wheat cropping system of Pakistan,” in Handbook of climate change and agroecosystems: The agricultural model intercomparison and improvement project integrated crop and economic assessments, Part 2, ed. Rosenzweig C. (Singapore: World Scientific), 219 – 258. doi: doi:10.1142/9781783265640_0019
- Ahmad M. I., Ma H. (2020). Climate change and livelihood vulnerability in mixed crop–livestock areas: the case of Province Punjab, Pakistan. Sustainability 12: 586. doi: doi:10.3390/su12020586
- Ahmed Z., Guha G. S., Shew A. M., Alam G. M. (2021). Climate change risk perceptions and agricultural adaptation strategies in vulnerable riverine char islands of Bangladesh. Land Use Policy 103: 105295. doi: doi:10.1016/j.landusepol.2021105295
- Anjum G., Fraser A. (2021). Vulnerabilities associated with slow-onset events (SoEs) of climate change: multi-level analysis in the context of Pakistan. Curr. Opinion Environ. Sustain. 50, 54 – 63. doi: doi:10.1016/j.cosust.02004
- Anwar A., Younis M., Ullah I. (2020). Impact of urbanization and economic growth on CO2 emission: a case of far east Asian countries. Int. J. Environ. Res. Public Health 17: 2531. doi: doi:10.3390/ijerph17072531
- Araro K., Legesse S. A., Meshesha D. T. (2020). Climate change and variability impacts on rural livelihoods and adaptation strategies in Southern Ethiopia. Earth Sys. Environ. 4, 15 – 26. doi: doi:10.1007/s41748-019-00134-9
- Arbaiah O., Daud A. R., Surinah A., Noorhaida U., Shaharom N. A. M. C. D., Rahim A.. (2021). Public health preparedness and response to flood disaster in Johore, Malaysia: challenges and lessons learned. Malaysian J. Commun. Health 15: 104870. doi: doi:10.1108/DPM-07-2013-0112
- Aryal J. P., Sapkota T. B., Rahut D. B., Krupnik T. J., Shahrin S., Jat M. L.. (2020). Major climate risks and adaptation strategies of smallholder farmers in coastal Bangladesh. Environ. Manag. 66, 105 – 120. doi: doi:10.1007/s00267-020-01291-8
- Ayanlade A., Radeny M., Morton J. F. (2017). Comparing smallholder farmers' perception of climate change with meteorological data: a case study from southwestern Nigeria. Weather Clim. Extre. 15, 24 – 33. doi: doi:10.1016/j.wace.12001
- Behera S. K. (2024). Understanding the impact of climate change on extreme events. Front. Sci. 2: 1433766. doi: doi:10.3389/fsci.2024.14337661433766
- Bennett N. J., Blythe J., Tyler S., Ban N. C. (2016). Communities and change in the anthropocene: understanding social-ecological vulnerability and planning adaptations to multiple interacting exposures. Region. Environ. Change 16, 907 – 926. doi: doi:10.1007/s10113-015-0839-5
- Birkmann J. (2006). “Indicators and criteria for measuring vulnerability: theoretical bases and requirements,” in Measuring Vulnerability to Natural Hazards: Towards Disaster Resilient Societies, ed. Birkmann J. (Tokyo: United Nations University Press), 55 – 77.
- Bryan E., Ringler C., Okoba B., Roncoli C., Silvestri S., Herrero M.. (2013). Adapting agriculture to climate change in Kenya: household strategies and determinants. J. Environ. Manag. 114, 26 – 35. doi: doi:10.1016/j.jenvman.10036
- Datta P., Behera B. (2022). Assessment of adaptive capacity and adaptation to climate change in the farming households of Eastern Himalayan foothills of West Bengal, India. Environ. Challenges 7: 100462. doi: doi:10.1016/j.envc.2022100462
- Derbile E. K., Bonye S. Z., Yiridomoh G. Y. (2022). Mapping vulnerability of smallholder agriculture in Africa: vulnerability assessment of food crop farming and climate change adaptation in Ghana. Environ. Challenges 8: 100537. doi: doi:10.1016/j.envc.2022.100537.
- Dhungana N., Silwal N., Upadhaya S., Khadka C., Regmi S. K., Joshi D.. (2020). Rural coping and adaptation strategies for climate change by Himalayan com- munities in Nepal. J. Mountain Sci. 17, 1462 – 1474. doi: doi:10.1007/s11629-019-5616-3
- Dogar M., Shahzad L., Khan S. (2025). “ Human-animal interactions in disaster management: insights from Pakistan and the global south,” in Palgrave Handbook of Human-Animal Interactions in the Global Context of Climate Change, Disasters, and Other Crises, eds. Wu H., Breen K., DeYoung S. E. (Cham: The Palgrave Macmillan). doi: doi:10.1007/978-3-032-00313-3_7
- Dorji T., Morrison-Saunders A. Blake, D. (2026). Navigating climate change in carbon negative Bhutan: insights from policy influencers and comparison to the wider Himalayan region. J. Mt. Sci. 23, 775 – 797. doi: doi:10.1007/s11629-025-9752-7
- Duffy C., Pede V., Toth G., Kilcline K., O'Donoghue C., Ryan M.. (2021). Drivers of household and agricultural adaptation to climate change in Vietnam. Clim. Develop. 13, 242 – 255. doi: doi:10.1080/17565529.20201757397
- Fahad S., Inayat T., Wang J., Dong L., Hu G., Khan S.. (2020). Farmers' awareness level and their perceptions of climate change: a case of Khyber Pakhtunkhwa province, Pakistan. Land Use Policy 96: 104669. doi: doi:10.1016/j.landusepol.2020104669
- Fahad S., Wang J., Hu G., Wang H., Yang X., Shah A. A.. (2018). Empirical analysis of factors influencing farmers crop insurance decisions in Pakistan: evidence from Khyber Pakhtunkhwa province. Land Use Policy, 75, 459 – 467. doi: doi:10.1016/j.landusepol.04016
- Faisal M., Abbas A., Cai Y., Ali A., Shahzad M. A., Akhtar S.. (2021). Perceptions, vulnerability and adaptation strategies for mitigating climate change effects among small livestock herders in Punjab, Pakistan. Int. J. Environ. Res. Public Health 18: 10771. doi: doi:10.3390/ijerph182010771
- FAO (2023). “ The state of food security and nutrition in the World 2023,” food and agriculture organization of the united nations. Available online at: (Accessed March 28, 2024).
- FAO (2023a). “ The impact of disasters on agriculture and food security”, food and agriculture organization of the united nations. Available online at: (Accessed March 28, 2024).
- FAO (2023b). “ Estimating global and country-level employment in agrifood systems. Available online at: (Accessed March 28, 2024).
- Field C. B., Barros V. R. (eds.) (2014). Climate change 2014–Impacts, Adaptation and Vulnerability: Regional Aspects. Cambridge: Cambridge University Press. doi: doi:10.1017/CBO9781107415379
- Fuglie K. O., Morgan S., Jelliffe J. (2024). World Agricultural Production, Resource Use, and Productivity, 1961–2020 (Report No. EIB-268). U.S. Department of Agriculture, Economic Research Service. Available online at:. doi: doi:10.32747/2024.8327789.ers
- Füssel H.-M., Klein R. J. T. (2006). Climate change vulnerability assessments: an evolution of conceptual thinking. Climat. Change 75, 301 – 329. doi: doi:10.1007/s10584-006-0329-3
- Gbadegesin T. K., Andrée B. P. J., Braimoh A. (2026). Effects of climate shocks on food security, food prices, and agricultural wages in Afghanistan. Environ. Dev. Sustain. doi: doi:10.1007/s10668-026-07834-z
- Gbetibouo G. A. (2009). Understanding Farmers' Perceptions and Adaptations to Climate Change and Variability: The Case of the Limpopo Basin, South Africa. Washington DC: International Food Policy Research Institute.
- Gbetibouo G. A., Ringler C., Hassan R. (2010). Vulnerability of the South African farming sector to climate change and variability: an indicator approach. Nat. Res. 34, 175 – 187. doi: doi:10.1111/j.1477-8947.2010.01302.x
- GOP (2019). Economic survey of Pakistan finance division, economic advisor s wing, Government of Pakistan, Islamabad. Available online at: (Accessed November 11, 2019).
- GOP (Government of Pakistan) (2020). Economic Survey of Pakistan 2019-20. Chapter: Agriculture. Available online at: (Accessed October 19, 2020).
- Hahn M. B., Riederer A. M., Foster S. O. (2009). The Livelihood Vulnerability Index: a pragmatic approach to assessing risks from climate variability and change-a case study in Mozambique. Glob. Environ. Change 19, 74 – 88. doi: doi:10.1016/j.gloenvcha.2008.11.002
- Hinkel J. (2011). “Indicators of vulnerability and adaptive capacity”: towards a clarification of the science–policy interface. Glob Environ Change-Human Policy Dimens. 21: 198 – 208. doi: doi:10.1016/j.gloenvcha.08002
- Hultgren A., Carleton T., Delgado M., Gergel D. R., Greenstone M., Houser T.. (2025). Impacts of climate change on global agriculture accounting for adaptation. Nature 642, 644 – 652 doi: doi:10.1038/s41586-025-09085-w
- Huong N. T. L., Bo Y. S., Fahad S. (2019). Economic impact of climate change on agriculture using Ricardian approach: a case of northwest Vietnam. J. Saudi Soc. Agricult. Sci. 18, 449 – 457. doi: doi:10.1016/j.jssas.02006
- IPCC (2007). “Climate change 2007: Impacts, adaptation and vulnerability,” in Contribution of Working Group II to the Fourth Assessment. Report of the Intergovernmental Panel on Climate Change, eds. Parry M. L., Canziani O. F., Palutikof J. P., van der Linden P. J., Hanson C. E. (Cambridge: Cambridge University Press), 976.
- IPCC (2023). “Summary for Policymakers,” in Climate Change 2023: Synthesis Report. Contribution of Working Groups I, II and III to the Sixth Assessment Report of the Intergovernmental Panel on Climate Change, eds. Core Writing Team H. Lee, Romero J. (IPCC: Geneva), 1 – 34.
- IPCC Climate Change (2014). Synthesis report. Contribution of working groups I, II and III to the fifth assessment report of the intergovernmental panel on climate change. IPCC, Geneva, Switzerland, 151.
- IPCC. (2018). Emissions scenarios: summary for policymakers. a special report of IPCC working group III. Published for the intergovernmental panel on climate change. Available online at: (Accessed April 10, 2025).
- Jamshidi O., Asadi A., Kalantari K., Movahhed Moghaddam S., DadrassJavan F., Azadi H.. (2020). Adaptive capacity of smallholder farm- ers toward climate change: evidence from Hamadan province in Iran. Clim. Develop. 12, 923 – 933. doi: doi:10.1080/17565529.20191710097
- Ketema A. M., Negeso K. D. (2020). Effect of climate change on agricultural output in Ethiopia. J. Persp. Pembiayaan Dan Pembangunan Daerah 8, 195 – 208. doi: doi:10.22437/ppd.v8i39076
- Khan M. A., Dogar M., Khalid M. (2024). Infrastructural development under CPEC and its impact on Pakistan's agricultural landscape. In 21st Century Landscape Sustainability, Development and Transformations: Geographical Perceptions, 1st Edn., Chapter 12. Cosenza: Il Sileno Edizioni.
- Khan N. A., Gao Q., Abid M., Shah A. A. (2021). Mapping farmers' vulnerability to climate change and its induced hazards: evidence from the rice-growing zones of Punjab, Pakistan. Environ. Sci. Pollut. Res. 28, 4229 – 4244. doi: doi:10.1007/s11356-020-10758-4
- Khan N. A., Gao Q., Iqbal M. A., Abid M. (2020). Modeling food growers' perceptions and behavior towards environmental changes and its induced risks: evidence from Pakistan. Environ. Sci. Pollut. Res. Int. 27, 20292 – 20308. doi: doi:10.1007/s11356-020-08341-y
- Khan N. A., Gong Z., Shah A. A. (2022). Synergy between climate risk perception, adaptation responses, and agricultural productivity: the case of rice farming communities in Pakistan. Environ. Sci. Pollut. Res. 29, 23750 – 23766. doi: doi:10.1007/s11356-021-17615-y
- Maddison D. (2007). The Perception of and Adaptation to Climate Change in Africa, Vol. 4308. Washington, DC: World Bank Publications. doi: doi:10.1596/1813-9450-4308
- Mazhar N., Shirazi S. A., Stringer L. C., Carrie R. H., Dallimer M. (2021). Spatial patterns in the adaptive capacity of dryland agricultural households in South Punjab, Pakistan. J. Arid Environ. 194: 104610. doi: doi:10.1016/j.jaridenv.2021104610
- Mehmood A. (2022). Narowal-Sialkot rail track washed away by flood from (Dawn Newspaper). Available online at: (Accessed February 21, 2024).
- Mehmood K., Rehman A., Khan A. (2021). Farmers' perceptions, awareness and adoption of improved groundnut varieties in Potwar Plateau of Pakistan. Sarhad J. Agric. 37, 1364 – 1376. doi: doi:10.17582/journal.sja/2021/37.4.1364.1376
- Mekonen A. A., Berlie A. B. (2021). Rural households' livelihood vulnerability to climate variability and extremes: a livelihood zone-based approach in the Northeastern Highlands of Ethiopia. Ecol. Proc. 10, 1 – 23. doi: doi:10.1186/s13717-021-00313-5
- Mramba R. P., Mapunda P. E. (2025). Assessing farmers' awareness of climate change impacts and adaptation practices in a semi-arid region in Tanzania. Discov. Sustain. 6: 877. doi: doi:10.1007/s43621-025-01831-4
- Naqvi S. A. A., Hasis ul Hassan R., Wu W., Shah A. A., Makhdum M. S. A., Shah S. A. R. (2020). Synergy between adaptations and resilience of livelihood from climate change vulnerability: a group-wise comparison of adapters and non-adapters. PloS ONE, 15:e0236794 doi: doi:10.1371/journal.pone.0236794
- National Disaster Management Authority (NDMA) (2018). National Disaster Management Plan. Government of Pakistan.
- Nida H., Kashif M., Janjua A. A., Aslam M., Cheema K. S., Ullah S. (2025). Impacts of climate change on Pakistan's weather patterns: a comprehensive study of temperature and precipitation trends. Environ. Monit. Assess. 197: 509. doi: doi:10.1007/s10661-025-13931-9
- OCHA (2019). Pakistan Drought Response Plan (Jan - Dec 2019). UNOCHA. Available online at: (Accessed August 17, 2024).
- Pakistan Bureau of Statistics (2025). Agricultural statistics of Pakistan; Government of Pakistan. Islamabad: Statistics Division.
- Pakistan Economic Survey (2024). Pakistan Economic Survey 2023-24, Finance Division, Government of Pakistan. Available online at: (Accessed August 12, 2024).
- Parker L., Bourgoin C., Martinez-Valle A., Läderach P. (2019). Vulnerability of the agri-cultural sector to climate change: the development of a pan-tropical climate risk vulnerability assessment to inform sub-national decision making. PLoS ONE 14: e0213641. doi: doi:10.1371/journal.pone.0213641
- Paul A., Deka J., Gujre N., Rangan L., Mitra S. (2019). Does nature of livelihood regulate the urban community's vulnerability to climate change? Guwahati city, a case study from North East India. J. Environ. Manag. 251: 109591. doi: doi:10.1016/j.jenvman2019109591
- Rehman A., Batool Z., Ma H., Alvarado R., Oláh J. (2024). Climate change and food security in South Asia: the importance of renewable energy and agricultural credit. Humanit Soc Sci Commun 11: 342. doi: doi:10.1057/s41599-024-02847-3
- Rehman A., Jingdong L., Shahzad B., Chandio A. A., Hussain I., Nabi G.. (2015). Economic perspectives of major field crops of Pakistan: an empirical study. Pacific Sci. Rev. B: Human. Soc. Sci. 1: 145 –158. doi: doi:10.1016/j.psrb.09002
- Riaz U., Shahzad L., Tonelli F. M. P. (2026). Concept of Resilience Related to Water Governance. In Water governance, 1st Edn. (Oakville, ON: Apple Academic Press), 9. doi: doi:10.1201/9781998511631-1
- Saaty T. L. (2008). Decision making with the analytic hierarchy process. Int. J. Serv. Sci. 1, 83 – 98. doi: doi:10.1504/IJSSCI.2008.017590
- Shahzad L. (2025). “Outbreak of epidemics during disaster response,” in The Routledge Handbook of Disaster Response and Recovery, eds. Paul B. K., Juran L. (Abingdon: Routledge), 202 – 216. doi: doi:10.4324/9781003414834-22
- Shahzad L., Ali M., Sharif F., Shedayi A. A. (2024). “ Managing Disasters in Mountains: Challenges in the Era of Global Warming,” in Warming Mountains, eds. Tripathi S., Bhadouria R., Garkoti S. C. (Cham: Springer). doi: doi:10.1007/978-3-031-62197-0_11
- Shahzad L., Shah M., Saleem M., Mansoor A., Sharif F., Tahir A.. (2021). Livelihood vulnerability index: a pragmatic assessment of climatic changes in flood affected community of Jhok Reserve Forest, Punjab, Pakistan. Environ. Earth Sci. 80, 1 – 16. doi: doi:10.1007/s12665-021-09562-1
- Shahzad L., Tahir A., Sharif F., Haq I. U., Mukhtar H. (2019c). Assessing the impacts of changing climate on forest ecosystem services and livelihood of Balakot mountainous communities. Pak. J. Bot. 51, 1405 – 1414. doi: doi:10.30848/PJB2019-4(1)
- Shahzad L., Tahir A., Sharif F., Hayyat M. U., Ghani N., Farhan M.. (2019a). Does livelihood vulnerability index justify the socio-economic status of mountainous community? A case study of post-earthquake ecological adaptation of Balakot population. Appl. Ecol. Environ. Res. 17, 6605 – 6624. doi: doi:10.4172/2157-7617.1000335
- Shahzad L., Tahir A., Sharif F., Khan W. U. D., Farooq M. A., Abbas A.. (2019b). Vulnerability, wellbeing and livelihood adaptation under changing environmental conditions: a case from mountainous region of Pakistan. Environ. Sci. Pollut. Res. 26, 26748 – 26764. doi: doi:10.1007/s11356-019-05880-x
- Shahzad L., Waheed A., Sharif F., Ghafoor G. Z., Rafique A. (2024) Understanding role of climatic parameters adaptation strategies in agriculture productivity of South Asian countries. Sustain. Environ. 10: 2345453. doi:10.1080/27658511.20242345453
- Shukla R., Gleixner S., Yalew A. W., Schauberger B., Sietz D., Gornott C.. (2021). Dynamic vulnerability of smallholder agricultural systems in the face of climate change for Ethiopia. Environ. Res. Lett. 16: 044007. doi: doi:10.1088/1748-9326/abdb5c
- Shukla R., Sachdeva K., Joshi P. K. (2016). An indicator-based approach to assess village-level social and biophysical vulnerability of agriculture communities in Uttarakhand, India. J. Mount. Sci. 13, 2260 – 2271. doi: doi:10.1007/s11629-016-4058-4
- Tate E. (2012). Social vulnerability indices: a comparative assessment using uncertainty and sensitivity analysis. Natural Hazards, 63, 325-347. doi: doi:10.1007/s11069-012-0152-2
- Teshome M. (2016). Rural households' agricultural land vulnerability to climate change in Dembia woreda. Northwest Ethiopia. Environ Syst Res. 5, 1 – 18. doi: doi:10.1186/s40068-016-0064-3
- Thakur S. B., Bajagain A. (2019). Impacts of climate change on livelihood and its adaptation needs. J. Agricult. Environ. 20, 173 – 185. doi: doi:10.3126/aej.v20i025067
- Thorlakson T., Neufeldt H. (2020). Reducing subsistence farmers' vulnerability to climate change: evaluating the potential contributions of agroforestry in western Kenya. Agricult. Food Secu. 1, 1 – 13. doi: doi:10.1186/2048-7010-1-15
- Tran P. T., Vu B. T., Ngo S. T., Tran V. D., Ho T. D. (2022). Climate change and livelihood vulnerability of the rice farmers in the North Central Region of Vietnam: a case study in Nghe An province, Vietnam. Environ. Challenges 7: 100460. doi: doi:10.1016/j.envc.2022100460
- Tripathi A., Mishra A. K. (2017). Knowledge and passive adaptation to climate change: an example from Indian farmers. Climate Risk Manag. 16, 195 – 207. doi: doi:10.1016/j.crm.11002
- Ullah F., Shah S. A. A., Saqib S. E., Yaseen M., Haider M. S. (2021). Households' flood vulnerability and adaptation: empirical evidence from mountainous regions of Pakistan. Int. J. Dis. Risk Reduct. 52: 101967. doi: doi:10.1016/j.ijdrr.2020101967
- Ullah H., Rashid A., Liu G., Hussain M. (2018). Perceptions of mountainous people on climate change, livelihood practices and climatic shocks: a case study of Swat District, Pakistan. Urban Clim. 26, 244 – 257. doi: doi:10.1016/j.uclim.10003
- UNICEF (2018) National Nutrition Survey 2018: Key Findings Report. United Nations Children's Fund. Available online at: (Accessed October 7 2025).
- Vincent K. (2007). Uncertainty in adaptive capacity and the importance of scale. Global Environ. Change 17, 12 – 24. doi: doi:10.1016/j.gloenvcha.11009
- Williams P. A., Crespo O., Abu M., Simpson N. P. (2018). A systematic review of how vulnerability of smallholder agricultural systems to changing climate is assessed in Africa. Environ. Res. Lett. 13: 103004. doi: doi:10.1088/1748-9326/aae026
- World Bank (2022). Pakistan: Flood Damages and Economic Losses Over USD 30 billion and Reconstruction Needs Over USD 16 billion - New Assessment. World Bank Group. Available online at:
- Yohannes H. (2016). A review on relationship between climate change and agriculture. J. Earth Sci. Clim. Change 7, 1 – 8. doi: doi:10.15666/aeer/1703_66056624
- Young O. R., Berkhout F., Gallopin G. C., Janssen M. A., Ostrom E., van der Leeuw S. (2009). The globalization of socio-ecological systems: an agenda for scientific research. Glob. Environ. Change 16, 304 – 316. doi: doi:10.1016/j.gloenvcha.2006.03.004
- Zamasiya B., Nyikahadzoi K., Mukamuri B. B. (2017). Factors influencing smallholder farmers' behavioural intention towards adaptation to climate change in transitional climatic zones: a case study of Hwedza District in Zimbabwe. J. Environ. Manag. 198, 233 – 239. doi: doi:10.1016/j.jenvman.04073
Citation
Shahzad, L., et al. (2026). Understanding climate vulnerability and economic impacts on rural livelihoods: lessons from Narowal, Punjab. Frontiers in Environmental Economics, Volume 5. https://doi.org/
- Received
- Accepted
- Published
- Volume
- 5 - 2026
Copyright
Copyright © 2026 Shahzad, Saleem, Sharif, Shedayi, Hayyat, Farhan and Dogar.
License: https://creativecommons.org/licenses/by/4.0/
*Correspondence: Laila Shahzad
Article metrics
- Views
- 0
- Downloads
- 0
- Citations
- 0
- Shares
- 0
