Next Article in Journal
The Effect of a Club in Making Differences in Knowledge, Attitude, and Practices on Family Planning Among Married Adolescent Girls in Urban Slums in Bangladesh
Next Article in Special Issue
Climate-Driven Adaptation, Household Capital, and Nutritional Outcomes among Farmers in Eswatini
Previous Article in Journal
Community-Based Healthcare Programs Sustainability Impact on the Sustainability of Host Organizations: A Structural Equation Modeling Analysis
Previous Article in Special Issue
Understanding the Operational Concept of a Flood-Resilient Urban Community in Jakarta, Indonesia, from the Perspectives of Disaster Risk Reduction, Climate Change Adaptation and Development Agencies
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Dairy Production under Climatic Risks: Perception, Perceived Impacts and Adaptations in Punjab, Pakistan

1
College of Economics and Management, Nanjing Agricultural University, Nanjing 210095, China
2
Institute of Agricultural & Resource Economics, University of Agriculture, Faisalabad 38000, Pakistan
3
Institute of Business Management Sciences, University of Agriculture, Faisalabad 38000, Pakistan
*
Author to whom correspondence should be addressed.
Int. J. Environ. Res. Public Health 2019, 16(20), 4036; https://0-doi-org.brum.beds.ac.uk/10.3390/ijerph16204036
Submission received: 17 September 2019 / Revised: 17 October 2019 / Accepted: 17 October 2019 / Published: 21 October 2019
(This article belongs to the Special Issue Climate Change Adaptation and Risk Reduction)

Abstract

:
The changing climatic conditions coupled with fodder availability have posed severe challenges and threats for the dairy sector in Pakistan. The current paper determines the influence of climate change on the dairy sector in Pakistan. Comprehensive data set was collected from 450 farmers. The majority of farmers experienced the climate change and its variability and explained that severity and frequency of climatic extreme events such as droughts, heat waves, floods, pests and diseases and humidity is increasing. The study found that farmers considered drought as one of the major climatic risks which severely affects all aspects of dairy production. Specifically, to estimate the perceived impacts of climatic extreme event on milk production, an ordered probit model was applied and identified that climate change had high adverse impact on milk quantity in the study area. Different adaptation practices, such as changing cropping pattern for fodder production, off-farm income activities, diversifying the farm and regular vaccination are mostly used by dairy farmers. The study recommends policy initiatives to be taken by government for long term developments in the dairy farming.

1. Introduction

Climate change refers to the long-term changes in temperature, precipitation, wind and other climatic features [1]. The Intergovernmental Panel on Climate Change (IPCC) in its 5th assessment report specifies that by the end of this century, average temperature on the surface of the planet may rise between 0.3 and 4.8 °C [2]. The increasing influence of climate change on crops and livestock rearing systems is more negative in nature [3]. Likewise crops, livestock are also almost sensitive that climate extreme events have extremely affected them in the past [4]. It is obvious that climate change and its severe impacts on productivity of livestock impose additional devastated risks to the worldwide food security. Latest annual report of IPCC (AR5) illustrates that the most tropical regions of world are not adapting the mitigated measures for climate change that exaggerate more negative effects on food production systems [3].
Climate change has imposed negative impacts on livestock productivity in addition to adulterated water resources, poor quality of feed and prone diseases of livestock [5,6,7,8,9]. Heat stress in livestock due to rising temperatures leads to an adverse impact on milk production [10,11], reproduction [12], health [13] and mortality rate of animals [14]. According to Houghton (2001) and Herbut (2018) [15,16] and Rust (2013) [17] air temperature, humidity and wind speed have significant direct effect on milk production and reproduction rate. Due to the persistent drought conditions, the lactation period of dairy cattle’s always shrinks. Likewise, milk production quantity and quality declines [18]. In Pakistan, livestock production is likely to reduce by 20 to 30% in the future because of increasing temperature spells. Ultimately, a crisis in dairy and meat industry will increase the prices of products that will not be in the reach of middle class consumers [19].
Pakistan is the state that have also been adversely influenced because of climate change, its vulnerability, low adaptive capacity of mitigating strategies and poor infrastructure at farm level [2]. There is future prediction that by the end of 2050, the temperature in Pakistan will increase by 2–3 °C, and the precipitation rate will also be significantly varied [20]. The Global Climate Risk Index (GCRI) ranks Pakistan at 8th position among the countries the most affected by climate change and extreme weather events from 1995 to 2014 [21].
As the livestock sector is the central pillar of Pakistan’s rural economy, it contributes around 60.5% in value addition of agriculture as well as 11.4% in overall Gross Domestic Product (GDP). The most marginal rural people consider livestock sector as the best source of income. More than eight million rural families are directly involved in raising livestock and derive up to 35% of their livelihoods from livestock [22]. The livestock sector is playing a vital role in the mitigation of poverty and in the form of gross value addition, enhancing 1333 billion rupees (8.495 billion USD) as foreign exchange reserves. Climate change and poor endowments of natural resources are the most significant development challenges confronting the dairy farming in Pakistan. Small and marginal farmers’ share in milk production is more than 70% and they are highly susceptible to unpredictable events of extreme weather [23].
Climate change is the main cause of weather-related hazards that adversely affect livestock production system especially in developing countries. Policy makers and researchers agree at this point that climate change significantly impact the livestock sector [5,24]. Unfortunately, in South Asia, limited research on potential impacts of climate related risks on livestock sector is available compared to the cropping sector [25,26].
Farmers’ activities and plans to take adaptation measures are more dependent on their perceptions about climate variability and risks as compared to actual climate scenarios [27,28,29]. Perception insights could be economical and technically feasible to design the comprehensive policies related to climate change and environmental protection [30,31,32]. Moreover, climate-related risk perceptions enable farmers to decide for effective adaptation of risk coping strategies [33,34,35,36].
Figure 1 illustrates the conceptual framework. The basic theme behind this conceptual framework is that climate change impact can be quantified by farmers’ perceptions and their capability to adopt the most effective risk coping strategies [37,38].
Considering the importance of livestock for small and marginal farmers, determination of perceived impacts of climate change on different aspects in dairy farming will be helpful to formulate adaptation policies for forthcoming challenges. Research studies in Pakistan on this topic are limited, so there is dire need to identify the perceptions of farmers about climatic risks and its influence on dairy production generally and on milk production specifically. Mostly, dairy farmers depend upon the pasture and crop residues for animal feeding in the South Asian region. Therefore, observational insights are deemed essential to identify the adaptation or coping strategies practiced by farmers, to reduce the climate related risks for their dairy farming system.
The main aim of the current study is to investigate the implications of the climate change on livestock productivity in Punjab, Pakistan including key sub-objectives as (1) to assess dairy farmers’ perceptions about climate related risks and its variability, and (2) to measure the perceived impact of climate related risks on dairy production systems, especially milk production.
This research article is organized in the following sections. Research data and methods are presented in Section 2. Furthermore, empirical results of the study are given in Section 3 while Section 4 has a detailed discussion of results. The last section concludes the study with some policy implications.

2. Material and Research Methods

2.1. Study Area

The research was carried out in the Punjab, which is the second largest and the most populated province in Pakistan and located in the semi-arid low land region [39]. The geographical area of Punjab has 20.63 million hectares, in which 59% area is under cultivation. Reason being selecting this province includes that Punjab contributes its share of 54% and 62% in the national GDP and national agricultural sector [40]. Punjab has the largest buffalo and cow population, at 16.019 million and 13.204 million in total livestock population, respectively. Punjab is also the major milk producing province with annual production of 36.23 million liters, which is about 67% of total milk production in Pakistan [22]. Mixed cropping along with livestock rearing is the common practice of farmers in Punjab [41,42]. In Punjab, the average minimum and maximum temperature per annum has been in range of 16.52 °C to 21.50 °C and 30.09 °C to 32.75 °C between the periods 1980–2018. In Punjab, 50% to 75% rainfall is linked with the monsoon season that lasts from June to August each year. The Pakistan Agriculture Research Council (PARC) map bifurcates Punjab into different agro-ecological zones, as depicted in Figure 2, and each zone has different rainfall pattern. Rain-fed and irrigated zones of Punjab receives the highest and lowest rainfall respectively [43]. Keeping this in mind, three districts from three different agro-ecological zones were selected for purposes of this study.
Most of the farmers in selected districts grow crops as well as rear livestock as the second best source of income. The rain-fed area of Punjab has the lowest number of dairy animals as we can see from Table 1.

2.2. Sampling Procedure and Data Collection Methods

A field survey was carried out to collect cross sectional data from September to December 2018. A well designed and structured topic guide/questionnaire was used to draw the valuable data by interviewing dairy farmers/respondents. A multi-stage stratified random sampling procedure was applied as shown in Figure 3. In the first stage, study region was selected and for this purpose, map of PARC indicating agro-ecological zones of Punjab province was used. The Punjab province is divided into four agro-ecological zones which are irrigated plains (low intensity rain-fed area), mixed cropping zone (moderate rain-fed area), Barani (high intensity rain-fed) region and Thal region (desert area). Researchers selected the first three zones. In the second stage, total three districts consisting of one district from each zone by considering the homogeneity and heterogeneity of climatic variability, cropping pattern and irrigation systems and cattle breed were selected. From each district, two random cities (tehsils/sub-district) were selected at third stage and in next stage, from each city; three union councils (UCs) were chosen for study purpose by considering distance between UCs and distance from the UC to the main city. According to the Government of Pakistan administrative system, each UC has some villages and in the last fifth stage, from each UC, two villages were chosen. Twelve to thirteen dairy farmers were approached for interview without considering the size of their land holdings and farm location. On the whole, 450 farmers and out of which 150 were interviewed from each selected district.

2.3. Data Analysis

2.3.1. Climatic Risk Perception

During the interview, respondents were probed about the awareness of climate change and the severity and intensity of different climate related risks. Questions regarding climate related extreme events and fluctuations in climate patterns were asked by answering categorically in the form of increasing, no change and decreasing during the last decade (2008–2018). The following relevant categories: (1) Animal health (2) Feeds and feeding (3) quality and quantity of milk product (4) Animal breeding (5) Herd Size and (6) Production cost were developed to determine the farmers’ perceived impact of climate change on cattle production. Likewise, a Likert scale from 1 (no impact at all) to 5 (very high impact) was used to determine the total influence on milk production. Descriptive statistics were applied to assess their perceptions about impacts of climate change. The graphs were developed by summarizing the frequencies of responses in%ages.

2.3.2. Empirical Model

On the basis of utility maximization theory, Random Utility Model (RUM) states that a choice decision may be made in which an individual i can choose among a set of alternative j [45]. In this regard, the j ordered categories were developed by five point Likert scale in the form of 1 with no impact on milk production, 2 with low, 3 with medium, 4 with high and 5 with very high impact. An individual farmer chooses the best alternative from a set of choices that maximizes utility. Moreover, the level of farmer’s utility is a function of socio-economic characteristics. An ordered probit model (discrete choice model) was applied to evaluate the factors influencing farmers’ perceptions of climate related risks on milk production in comparison with linear model. Because the values of dependent variable are not continuous in nature but discrete and ordered [46] and it also avoided to detect a non-existing effect as well as loss of power problem like fail to detect an effect [47]. The specified model is given in Equation (1).
M P × i = i = 1 14 X i β   +   ϵ i                   i = 1 ,   2 , . . ,   n
MPi* is the latent unobserved variable that corresponds perceived milk production affected by climate related risks, Xi is the vector of observed explanatory variables (climatic parameters, farmers’ and farm characteristics and institutional links) of the ith farmer, the unknown parameter is estimated as β while the random term of the latent utility function is denoted as ϵ i . The details of the model are given in the Appendix A.

3. Estimated Results

3.1. Dairy Farmer’s Socio-Economic and Farm Characteristics

Table 2 presents the socio-economic and farm characteristics of the sampled farmers. The farmer’s average age was around 45 years, while the mean dairy farming experience was of 21 years. It is clear that ranchers were elderly with extensive dairy experience. The growing age of Pakistani farmers is becoming a concern for the agricultural sector [48]. Some researchers use age as proxy variable of agricultural experience [49] and some believe that farmers observe climate changes definitely as their farming knowledge increases [50].
The family size of the respondents was large as main rural families of Pakistan lived in joint family structures [9]. The average size of family members was 8 with 3.6 family members’ variation. Dairy farmers earn about 42% of their income via dairy farming which is major one and the daily-based source of income, followed by crops, labor and public and private sector employment. These figures, confirmed by the Government of Pakistan in 2018, estimates that farmers derive 35% of their annual income from the livestock sector. Around 69% of the respondents are literate at school level, while only 11% attended the college or university and rest (20%) has never been in school.
The average farm size was 2.9 hectare and holds 4 to 5 milking animals with 14 L milk production per day which varied about 6.5 L due to different reasons. About 52% and 47% of the farms were located in irrigated region and rain-fed area of Punjab province respectively. Most of the dairy farmers (51%) have mixed breed of animals and 38% of farmers are still relying on indigenous/dessi breed, while only 11% are rearing crossbreed animals. Farmers believe that indigenous breeds are more resilient to climate related risks [51].

3.2. Dairy Farmers’ Perceptions about Climatic Risks and Variability

During the interview, dairy farmers were also queried about their perceptions that how they associate different sources of climate related risks with their regular farming practices as presented in Figure 4. Dairy farmers perceived that heat waves were the major threat for their dairy animals followed by droughts, pest and diseases, flood and humidity. About 105 (23%) farmers consider heat waves as one of the major climate related risk sources. Because temperature range (maximum and minimum) in Pakistan is increasing day by day, that severely affects the farming system.
The second major risk was drought, which also adversely impacts the livestock and fodder production system. Around 73 (16%) respondents reported that animals’ diseases’ severity and sensitivity has increased over the time period and animal health and their productivity and reproduction cycle has altered due to number of biological diseases emerged by fluctuating climatic conditions. Fourteen% of dairy farmers consider flood a significant risk due to a series of floods from 2010 to 2014. Nearly 12% and 8% of farmers identified humidity and heavy rainfall as major climate related risks that affect the dairy production system in negative manner.
The farmer’s perception about climate variability is observed as overall understanding of the following climate parameters. During the survey, the question was raised to the dairy farmers as to whether they have perceived some fluctuations in precipitation and temperature during last decade and either the length of summer season has increased, decreased or there has been no change during the last decade. Simultaneously, the same approach was used for getting perceptions about the temperature in the winter season, while for the rainy season, we inquired that either monsoon has been started earlier, later or continuing on the same pattern during the said time period. All the responses reported by the surveyed farmers are summarized in Figure 5.
It was a common perception of almost all farmers that the time interval in summer season had been extended. More than 85% farmers reached on the consensus that winter season temperature is increasing that means the frequency and intensity of winter has been decreased in the Punjab province. Around 25% of farmers observed that raining season has been starting earlier every year while 70% argued that it starts later and only 5% observes that there is no change in precipitation pattern.
On the other hand, 66% of farmers observed that monsoons period terminate earlier but 27% had conflicts in their views about ending of rainy season. Less than 20% of the dairy farmers were of the view that length of rainy season has increased and almost 75% considered that span of rainy season has decreased and 8% respondents observe no change in length of season. Most of the farmers (85%) observed that the amount/quantity of rainfall has been decreased over the period of the last 10 years.
Figure 6 verifies the perception of the farmers, as we can see that since 2000 onward temperature has an increasing trend with minor variation while the precipitation has too much variation which very disturbing for the dairy farmers as heavy rainfall cause of floods and very low precipitation increased the drought time span. On the whole, we can conclude that farmers correctly perceived the changes in climate.

3.3. Perceived Impacts of Climatic Risks on Dairy Production System

We measured the perceived impact of pre-defined climate related risks on some aspects (feeding of animals, health of animals, breeding of animals, herd size, milk quantity and quality and production cost) of dairy production system and results are presented in Figure 7 Almost all respondents (98%) reported that climate change has affected their overall dairy production system during last decade. One of the main climate related risks is drought, which affected every aspect of the dairy production significantly. Dairy farmers perceive that heat waves and pests and diseases also posed severe adverse impact on all components of dairy farming. More than 20% of farmers reported that whenever they experienced with floods, ultimately herd size, health and production cost was hurtled. Dairy farmers revealed that during floods, liver fluke infection is the main disease observed in dairy animals.
Most of the respondents believe that during heat stress, milking animals collapse with mastitis disease. During a humid period, milk production decreased noticeably and breeding of animals is also affected. Many farmers observe that during humidity, conception rate in animals is declined and also have the opinion that during heat stress and dry period female cattle are not ready for mating process, which disturbed pregnancy and reproduction rate. About 7% of farmers also recalled that heavy rainfall in unexpected time disturb the heated animals for mating purpose. Heavy rainfalls increased the cost of production at the farm because they bear some unexpected expenditures during continuous rainfall and not able to perform the routine activities for animals. Fog and smog affect the health of animals significantly that is one of the main aspects of dairy farming system. Farmers were on the view that almost all climatic risks increase their production costs because they invested massively to maintain milk production and buy some extra equipments at farm level.

3.4. Estimates of Ordered Probit Model

We used an ordered probit model to evaluate perceived impacts of climate related risks on milk production at the farm level. This robust model was tested with specific tests to check the possibility of multicollinearity, heteroscedasticity and goodness of fit [52,53]. Variance Inflation Factor (VIF) was applied for multicollinearity and showed that multicollinearity did not exist in the model. Value of Chi square was highly significant for log likelihood function and coefficients of all explanatory variables were significantly differed from zero. The link test showed that the model was correctly specified.
Wollni (2010) [54] suggested that ordered probit model coefficients are very complex to explain. Estimated parameters of explanatory variables in model only show the direction (positive or negative) of impact on dependent variable and do not measure the actual change in magnitude or probabilities. Thus, marginal effects of the model are however easy to interpret the results. From the five perception levels, ‘no impact’, ‘low’, ‘medium’, high’ and ‘very high’ impact on milk production; we estimated marginal effects for only first four categories because only five farmers out of 450 respondents chose the last perception ‘high impact’ and estimated results of marginal effect for this category were not significant for all explanatory variables so we did not include it in Table 3. Here, we interpreted the marginal effect estimates as effect of changes in explanatory variables on the expected likelihoods underneath the category (high impact) of dependent variable.
Overall, 14 predictors were examined in the model and out of which 10 predictors showed significance for perceived impact of climate change. All variables had the positive signs and all climatic risk variables increased the degree of perceived impacts. The probability of perceiving ‘medium’ and ‘high’ impact on milk production increased by 2% and 3% respectively due to drought, while the probability of ‘no impact’ and ‘low’ impact on milk production decreased about 1.4% and 4% respectively as drought risk increased. All the remaining climatic risks indicating the similar probable impact on milk production with minor variation except heavy rainfall where the results are insignificant.
Education significantly influenced the perception of dairy farmers for climate change effects. Likelihood of perceiving ‘medium’ and ‘high’ impact on milk production increased about 4% and 6% respectively among educated farmers, whereas increase in number of schooling years decreased the probability of ‘no impact’ and ‘low impact’ by 3% and 8%, respectively.
The farmers with good dairy farming experience likely perceive more impact on milk production due to climate change by 1%. While, Addition of 1 year in dairy farming experience reduce farmer’s probability for perceiving ‘low’ or ‘no impact’ on milk production about 0.5% and 0.2% respectively. Although farmer’s age showed negative signal for ‘no impact’ and ‘low impact’, but it was a positive sign for ‘medium’ and ‘high’ impact (indicating that aged farmers’ perceived more about climate change impacts on milk production). However, variable results were statistically insignificant.
The number of milking animals increased the production at farm significantly and enhanced the income of the farmers. If the respondent had more milking animals at farm, probably s/he perceived the more impact (‘medium’ and ‘high’) of climate change on milk production about 5% and 6% respectively. In contrast, the probability of perceiving ‘low’ or ‘no impact’ decreased by 0.7% and 0.3% accordingly. Breed of animals had a significant impact on milk production. As the cross breed of cattle had more productivity of milk as compared to local breed. In selected districts, farmers mostly had indigenous/local breed and they perceived
The minor influence of climate change on milk production is because the local breed was considered to be resilient against some climatic risks. However, the results of breed of animals’ variable were not statistically significant in the model.
Institutional links had a strong influence on the perceived impact of climatic risks as they think that their role is to disseminate the relevant information to the society. Farmers’ participation in social activities increased the probability of perceiving higher impact on milk production about 1% and its perceiving probability for ‘no’ or ‘low’ impact also decreased with 0.5% and 1.5%, respectively. However, it was a non-significant result. Sources of information had greater role about the farmer’s perception for climate related risks. As the farmers had quick and accurate information, then, they would be able to perceive the severity of impact timely and may try to diminish the negative impact of risks by suitable adaptation practices.
The probability of farmers about climatic risks impact ‘high’ on milk production increased by 2.3%, at the same time the perceiving probability of ‘no impact’ or ‘low impact’ decreased about 1% and 3% respectively. Dairy farmers in regular interaction with agricultural extension officers were well informed in advance about climatic risks and they perceive more accurately as compared to others. Contact/access to veterinary officer or extension officer significantly increased the probability of perceiving ‘higher’ impact on milk production by 5%. On the other hand, they perceived ‘low impact’ about 6.4% if they regularly approaching the extension officers.

3.5. Risk Coping Adaptation Strategies

Perception of climate related risks as a major challenge enable dairy farmers to adapt certain coping strategies [9]. Researchers observe that dairy farmers always try to adapt good measures against the negative impacts of climate change but not in a consistent way [55]. Most of the farmers capitalize their dairy farming as compared to cropping sector heavily affected by climatic extreme events [26,48,56]. Farmers solemnly believe that cropping sector is more sensitive than dairy production. That is why they consider this opportunity as good income source in extreme weather conditions and minimize crop failures. Farmers’ are keen to adapt ex-ante risk coping strategies to reduce the negative impacts of climate related risks because such responses enhance their earnings and stabilize their livelihoods [57].
A summary of adaptation strategies of dairy farmers is given in Figure 8. Generally, farmers adopted the following risk coping strategies: changing the cropping pattern for fodder production, migration, off-farm income activities, regular vaccination and diversification of farm and selling of weak/diseased animals. Most of the farmers applied these adaptation practices against drought situation because droughts severely affected their production at farm level.
More than 50% of respondents change the cropping pattern for fodder production during drought period and while more than 15% of farmers also practice this strategy in peak summer and rainy season [39]. About 62% of farmers migrated from one place to other with their herd in the search of fodder and water. During floods some were forced to migrate for their own and animals’ survival [9,58].
Most of the farmers, about 75%, earn an income from selling their labor to non-agricultural activities in cities, in this way they spread the risk and also enables them to invest in their farm in the near future. They diversified their income by working as daily laborer, petty traders, and construction worker or in some industry as wage earner, so income from these activities enables farmers to modernize their production system during unexpected hazards [59,60,61]. Heat waves, pests and diseases, droughts and floods are the major risks which cause the increase of infectious disease and parasite occurrence in cattle. Therefore, farmers use regular vaccination during these climatic risks. Famous livestock breeds and crops cultivars are more resistant against heat waves and droughts. Such strategy of diversification enhances animal productivity during high spells of rainfall and temperature. Furthermore, such diversification in variety of crops and animal species (for meat purpose) is an effective way to combat climate-related diseases and pest outbreaks [7,62,63]. About 45% of farmers sell their animals when they become weak or diseased, particularly during the drought situation when fodder and water shortage increased and they have to buy it from market. Marginal farmers adopt this strategy frequently as they were not able to feed and vaccinate their diseased animals due to lack of financial resources [9].

4. Discussion

4.1. Perceived Climatic Risks and Variability

In Pakistan, dairy farmers are very much concerned about climatic risks and its variability. Their past experience expose that they are well aware of occurrence and severity of climate. It is difficult to comprehend the nature of climate that poses a major threat for dairy production system as well as livelihoods of farmers [7,64]. According to [65], changes in the dairy production system that are associated with the climate change are important factors, for example, economically. Therefore, researchers believe that it is essential to know the perceptions of dairy farmers about climate related risks, as this may provide a profound familiarity with susceptible nature of farmers and their adaptation behavior.
Dairy farmers have brilliant memories of the events which cause severe damage to their farming system, as mentioned in some earlier studies [26,66,67]. Older farmers can easily relate the climate change to production changes. Some national studies [39,68,69,70] explained that farmers’ perceptions are based on past disaster experiences and information of predicted climate change scenarios. Farmers’ understanding of climate change could be seen in this context. Almost all dairy farmers in Pakistan exceptionally observe climate change and its variability. This reality advocates that farmers’ knowledge of climate change and their experience of extreme weather events can improve their perception [71,72].
In general, findings of this study confirm that farmers have good knowledge about climate change and its variability at farm level and their perceptions are similar according to the severity and frequency of floods [73] and sometimes are more accurate with the trends of temperature and precipitation [74].

4.2. Perception of Climatic Risks on Dairy Production System

Keeping in view the awareness level of farmers about climate variability, researchers further tried to find out that how farmers perceive the risks against dairy production system. In South Asian countries, the huge population (25% of Global Population) of dairy animals is exposed to climate change, but it has not been recorded properly [75]. Furthermore, farmers in these countries believe that indigenous breed of cattle and buffalo have the natural ability to adjust themselves with hot and humid environment. Although several researchers have focused their intentions on the biological functions of dairy animals and found that climatic factors like heat waves have more vulnerable effects [18,51,76,77]. In this study, farmers understand the nature of climate related risks by observing productivity and changes in physical conditions of animals. Dairy farmers are very experienced to identify the impact of climate related risks on dairy animals through their observation, as some of them quoted that “animals do not tell a lie”. They can easily observe some special signs that appeared in animals to identify either they are sick or trapped by heat stress because the “body language of animals always speaks to them”. Farmers are more certain to report that which particular aspect of dairy production system is affected by climate related risks that ultimately leads to enormous damage their farm productivity as well as their revenues. Feeding and animal’s health are the two major aspects of the dairy production system that are highly affected by climate-related risks. On the other hand, milk quantity and quality is severely diminished due to pests and infectious diseases, droughts and heat waves.
Climate-related risks enhance the probability of vulnerable feeding process of animals and affect the feed quality and quantity. Farmers perceive that feeding is the main and most expensive element of dairy farming, therefore many researchers have tried to study it [78,79,80]. Fodder production is severely decreased due to increasing temperature, high CO2 intensity and pests and diseases. In many tropical regions, length of fodder production period is also reduced accompanied by more significant changes in precipitation patterns and more frequent droughts [81]. Due to the limited land available for small farmers in Pakistan, the provision of pasture and fodder availability is becoming problematic day by day [9]. Some farmers purchase additional feed supplies from market, but most of them are not able to afford it. Local research focusing on the feed quality, availability and improvement in feeding methods for dairy animals is either very limited or under-reported.
Climate change has an adverse effect on livestock’s health conditions in direct and indirect ways. Vector-borne diseases and factors associated with soil, floods and humidity has directly impacted animal health. Furthermore, [82] and [83] found that there was indirect effect of climate change on animal’s health, because climate change encourage the adaptation of microbes that spread the vector-borne diseases and reduced the animal’s immunity against infectious virus due to water and fodder shortages. Floods are considered severe threat, because when animals are exposed to wet and submerged pastures, probability of occurrence of parasitic liver disease is increased [8,37]. About 3000 animals perished and 2 million animals were treated by vaccination in the 2014 drought in Thar-Parkar, in a remote area of Sindh province in Pakistan [84].
Heat waves are ranked one of the important climate related risks due to severe climatic conditions in region. Whenever asked about heat waves, farmers mistakenly mixed their perception because heat waves are mainly observed during droughts. They also recalled that milk quantity and quality is always affected by heat waves during the last decade [8,16,25,79,85,86]. More than 120 respondents confirm that breeding is highly affected by heat waves, humidity and drought period. Most of them opined that in the dry season, the fertility of cattle is reduced, while semen quality in bulls increased in silent heat. They also considered the poor conception rate as one of the major reasons for the slow breeding process [7,37].
Although climate related risks such as heat waves, humidity, pests and diseases and drought affect the dairy farming system severely, but indigenous breeds are more adaptive to climate change. In short, farmers perceive that climate related risks have adverse impact on production system of animals.

4.3. Impact on Milk Production

To identify which climatic risk severely affect the milk production at farm level according to farmers’ perception, we included six climatic parameters: drought, flood, heat waves, humidity, pests and diseases and heavy rainfalls. The estimated results showed that dairy farmers perceive drought as the major climatic risk, followed by heat waves, pests and diseases, flood and humidity. The frequency of droughts and floods vary substantially with change in the frequency and intensity of rainfall. [18] concluded that due to the persistent drought conditions, the lactation period of dairy cattle’s shrinks and milk production quantity and quality also declined. In a study from Tanzania about 95% cattle owners were worried about unpredicted rainfalls as well recurrent drought periods which adversely affect the milk production and growth in cattle [87].
Heat stress in livestock due to rising temperature lead to an adverse impact on milk production [10,11]. Similar results were determined by other authors that milk production decreased due to increasing temperature and humidity [15,16,37,88]. An Indian study [89] estimated that due to different diseases milk production and weight of cattle reduced by 46% and 9%, respectively. Some other studies also found that farmers’ perceiving higher degree of severity about diseases and parasite occurrence which reduced milk production at farm level [90,91]. According to [15,16] and [17], air temperature, humidity and wind speed have significant direct effect on milk production and reproduction rate. Whenever floods occur it severely damaged the agricultural productivity (crops) and livestock sector especially. In 2010, flood in Pakistan spread over 132 thousand square kilometers damaged the livestock sector of worth US$ 441 Million [92]. The overall economic losses in livestock sector in 2014 flood were estimated about 0.23 billion rupees.
Two out of three variables of farmer characteristics showed significant perceived impact on milk production. Both education and farming experience had positive significant effects, whereas age had insignificant impact. The estimated results suggested that highly educated farmers show better understanding of climatic risk perception and its impact on production. Several studies described the role of education to combat the adverse impact of climatic risks [68,93]. We also determined in the study that dairy farming experience and age had contrasting impact on perception of farmers. The older and less experienced farmers perceived more impact of climatic risks on their dairy production system. Our results contradicted with [94] who considered the age variable as a proxy for farming experience. Our results are consistent with [8,50,95] as they concluded that experienced farmers perceive impacts of climate change in true spirits. Besides, another study [96] pointed out that age was not the best predictor of awareness for climatic variability and experience was weakly correlated with age.
For perceiving the climatic risks impacts, some farm topographies had better linkage, such as the number of milking animals and breed of animals. Farmers’ with more milking heads perceived higher impacts of climatic risks. A relevant conclusion was reported in Philippine study [8]; that herd size affected by different climatic risks, particularly drought and floods affected the larger farms. Other study from East Africa had contradicting results, because they estimated that farmers with larger herd size had better capacity and adaptation strategies against climatic hazards impacts [97]. The crossbreed cattle were found to be more sensitive than local and mixed breed; an Indian study verified that temperature and humidity affected the lactation length and lactation yield of milk about 46.50% and 48.70% in ‘Holdeo’ crossbreed cattle [86]. While our study results for breed of animals were insignificant as in Pakistan farmers prefer local breed than crossbreed.
Farmers’ linkage with institutions had greater impact on perception level. Two (source of information and contact with extension services) out of three variables in this category significantly perceive the climatic risks impacts on milk production. Respondents who had quick and reliable sources of information, perceive more impact than those who have less information or outdated information about climate changes. Many studies estimated the similar results that weather information had a bigger role in shaping the farmers’ perception about climate change [8,98,99]. Frequently contact with extension services also enhance the perception of farmers for climatic impacts on their farms as well their livelihoods. Extension services disseminate the accurate and timely information about weather and also suggests appropriate measure (according to farmers’ financial status) to mitigate the adverse impact of climatic changes. Several researchers around the world reached at the same conclusion that provision of extension services improve resilient nature of farmers against risks [49,100,101,102,103]. To conclude, the farmers who perceived higher impacts of climatic changes, also had higher intent to adapt.

5. Conclusions

This study analyzed the climatic risks, farmers’ perceptions, and perceived impacts on dairy farming systems, specifically in terms of milk production and risk coping adaptations practiced by the dairy farmers. Clearly evident is that dairy farmers have good awareness about climate change and its relation to extreme weather events. Most of the farmers observed that the frequency and severity of climate related risks such as droughts, floods, heat waves, pest and diseases and humidity has increased in the study area and reflects their perceptions according to the actual trends of temperature and precipitation. Farmers perceived that some climatic events like heat waves, drought, pests and diseases, floods and humidity affected dairy production system severely. However, drought was considered the most dangerous climatic risk as it affected the feeding, herd size, milk production and production cost severely. An ordered probit model showed that most of the farmers perceived that climatic risks reduced the milk production by 20% to 30%. Similarly, their income level is also affected due to reduction in milk production. To sustain their livelihoods and production, farmers adopted some risk coping adaptation strategies. More than 50% of sampled farmers changed their cropping pattern for fodder production, practice migration and tried to maintain the farming activities as well as their livelihoods needs, they do off-farm work and sell their diseased/weak animals.
At present, Pakistan’s agricultural policy focuses on improving the resilient crop production systems against climate change. This research suggests that farmers’ reliance on livestock sector is increasing, so, policies should be designed in such way that could help the livestock farmers, especially the dairy farming system, and these policies must be long term climatic resilient and locally adaptable, particularly for small scale farmers. However, the provincial government must focus on the education of farmers with the collaboration of agricultural educational institutes as well research organizations, which facilitate the farmers with short courses relevant to livestock development. Finally, this study also suggests that government should enhance the role of extension services (training programs, exhibitions of model farms, farmer to farmer contact, etc.) and media (social media, TV programs and mobile phone applications) as they motivate farmers to adopt measures against climatic risks.

Author Contributions

Data curation, R.U.; Writing—original draft, Q.A.; Writing—review & editing, J.H. and A.A.

Funding

National Natural Science Foundation of China: 71973066 and Nanjing Agricultural University, China: KYGB201802

Acknowledgments

The authors acknowledge the sponsorship of “A Project Funded by the Priority Academic Program Development of Jiangsu Higher Education Institution (PAPD)”, KYGB201802 funded by Nanjing Agricultural University and Project “Study on sampling inspection efficiency, information disclosure mode and policy optimization of food safety regulation” (No. 71973066) supported by NSFC

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

The specified model is given in Equation (A1).
M P ˙ i = i = 1 14 X i β   +   ϵ i                                     i = 1 ,   2 , . . ,   n  
MPi* is the latent unobserved variable that corresponds perceived milk production affected by climate related risks, Xi is the vector of observed explanatory variables (climatic parameters, farmers’ and farm characteristics and institutional links) of the ith farmer, the unknown parameter is estimated as β while the random term of the latent utility function is denoted as ϵ i.
The details of the model are given in supplementary section.
MPi* is unobservable and we therefore observe:
MP i   = { 1   i f                           M P i *   μ 1 2   i f   μ   1   <   M P i *     μ 2   3   i f   μ   2   <   M P i *     μ 3 4   i f   μ   3   <   M P i *     μ 4 5   i f                                 M P i *     μ 4 .  
where MPi = when farmers perceive ‘no impact’ ‘low’ ‘medium’ ‘hh’ or ‘very high’ by observing the statement that milk production is affected by climate change at farm level.
The unknown parameters μs are jointly assessed with β-coefficients. It is assumed that the random term of the ordered probit model is trailed as standard normal distribution. The maximum likelihood estimation methods is used to estimate the model with the probability specified as follows;
P ( M P = 1 | X ) = F ( β X i ) P ( M P = 2 | X ) = F ( μ 1 β X i ) F ( β X i ) P ( M P = 3 | X ) = F ( μ 2 β X i ) F ( μ 1 β X i ) .   P ( M P = 4 | X ) = F ( μ 3 β X i ) F ( μ 2 β X i ) P ( M P = 5 | X ) = F ( μ 4 β X i ) F ( μ 3 β X i ) .  
where F is the cumulative standard function with normal distribution. P is the likelihood of farmer for choosing either ‘no impact’, ‘low’, ‘medium’, ‘high’ or ‘very high’ with the help of given X variables.
X is vector of independent variables that affect the perception levels of farmers.
β is vector of unknown parameters to be evaluated.

Determination of Marginal Effects

Ordered probit model narrates that as independent variables change, the coefficients from the regression are not the marginal changes in dependent variables. Therefore, it is a non-linear regression model. Marginal effects are calculated to evaluate the marginal changes. In this case, the marginal effects can be figured out as follows;
P ( M P = 1 | X ) X j = f ( μ 1 β X i ) β j P ( M P = 2 | X ) X j = f ( μ 1 β X i ) β j f ( μ 2 β X i ) β j P ( M P = 3 | X ) X j = f ( μ 2 β X i ) β j f ( μ 3 β X i ) β j P ( M P = 4 | X ) X j = f ( μ 3 β X i ) β j f ( μ 4 β X i ) β j .   P ( M P = 5 | X ) X j = f ( μ 4 β X i ) β j f ( μ 5 β X i )
where f (.) is a density function of standard normal variable.
P ( . ) X j = . The change in dependent variable given X as independent variable of the jth probability
MP i * ( impact )   =   β 1 op X 1 op   +   β 2 op X 2 op   +   β 3 op X 3 op   +   β 4 op X 4 op   +   β 5 op X 5 op   +   β 6 op X 6 op   +   β 7 op X 7 op   +   β 8 op X 8 op   +   β 9 op X 9 op   +   β 10   op X 10 op   +   β 11   op X 11 op   +   β 12   op X 12 op   +   β 13 op X 13 op   +   β 14 op X 14 op   +   μ i
where MPi*(impact) = Perceived impact on milk production: 1 with no impact, 2 with low, 3 with medium, 4 with high or 5 with very high impact. X1 =Severity of Drought, 1 less/not severe, 2 not sure/normal severe and 3 extreme severe; X2 = Severity of Flood, 1 less/not severe, 2 not sure/normal severe and 3 extreme severe; X3 = Severity of heat waves, 1 less/not severe, 2 not sure/normal severe and 3 extreme severe; X4 = Severity of Humidity, 1 less/not severe, 2 not sure/normal severe and 3 extreme severe; X5 = Severity of pest and disease, 1 less/not severe, 2 not sure/normal severe and 3 extreme severe; X6 = Severity of Heavy rainfall, 1 less/not severe, 2 not sure/normal severe and 3 extreme severe; X7 = Education (years); X8 = Age (years); X9 = Dairy Farming Experience (years); X10 = Milking Animals (No. of Heads); X11 = Breed of Animals (1 Indigenous/local breed, 2 Mixed breed and 3 Cross breed); X12 = Social Participation: Member of village committee or member of any farmers’ organization/milk cooperative society/ political/NGO (1 never participation, 2 sometimes and 3 frequently; X13 = Sources of information about weather related (1 from extension office, 2 friends & relatives, 3 Newspaper, 4 TV and 5 Internet/Mobile); X14 = Contact with extension services: contact with veterinary officer/extension officer (1 Never, 2 sometimes and 3 frequently) βs = Coefficient Parameters to be estimated; μ = Error term.

References

  1. US Environmental Protection Agency. Climate Change: Basic Information. In Clim. Chang. Div. Available online: https://www3.epa.gov/climatechange/basics/ (accessed on 23 June 2019).
  2. Stocker, T.F.; Qin, D.; Plattner, G.-K.; Tignor, M.; Allen, S.K.; Boschung, J.; Nauels, A.; Xia, Y.; Bex, V.; Midgley, P.M. Climate Change 2013: The Physical Science Basis; Cambridge University Press: Cambridge, UK, 2013. [Google Scholar]
  3. Porter, J.R.; Xie, L.; Challinor, A.J.; Cochrane, K.; Howden, S.M.; Iqbal, M.M.; Lobell, D.B.; Travasso, M.I. Food security and food production systems. In Climate Change 2014: Impacts, Adaptation, and Vulnerability; Cambridge University Press: Cambridge, UK, 2014. [Google Scholar]
  4. Food and Agriculture Organization of the United Nations. The Impacts of Disasters on Agriculture and Food Security. Available online: http://www.fao.org/resilience/ (accessed on 20 May 2019).
  5. Thornton, P.K. Livestock production: Recent trends, future prospects. Philos. Trans. R. Soc. B Biol. Sci. 2010, 365, 2853–2867. [Google Scholar] [CrossRef] [PubMed]
  6. Polley, H.W.; Briske, D.D.; Morgan, J.A.; Wolter, K.; Bailey, D.W.; Brown, J.R. Climate change and North American rangelands: Trends, projections, and implications. Rangel. Ecol. Manag. 2013, 66, 493–511. [Google Scholar] [CrossRef]
  7. Rojas-Downing, M.M.; Nejadhashemi, A.P.; Harrigan, T.; Woznicki, S.A. Climate change and livestock: Impacts, adaptation, and mitigation. Clim. Risk Manag. 2017, 16, 145–163. [Google Scholar] [CrossRef]
  8. Escarcha, J.F.; Lassa, J.A.; Palacpac, E.P.; Zander, K.K. Understanding climate change impacts on water buffalo production through farmers’ perceptions. Clim. Risk Manag. 2018, 20, 50–63. [Google Scholar] [CrossRef]
  9. Rahut, D.B.; Ali, A. Impact of climate-change risk-coping strategies on livestock productivity and household welfare: Empirical evidence from Pakistan. Heliyon 2018, 4, e00797. [Google Scholar] [CrossRef]
  10. Bohmanova, J.; Misztal, I.; Cole, J. Temperature-humidity indices as indicators of milk production losses due to heat stress. J. Dairy Sci. 2007, 90, 1947–1956. [Google Scholar] [CrossRef]
  11. Hammami, H.; Bormann, J.; M’hamdi, N.; Montaldo, H.H.; Gengler, N. Evaluation of heat stress effects on production traits and somatic cell score of Holsteins in a temperate environment. J. Dairy Sci. 2013, 96, 1844–1855. [Google Scholar] [CrossRef] [Green Version]
  12. Hansen, P.J. Effects of heat stress on mammalian reproduction. Philos. Trans. R. Soc. B Biol. Sci. 2009, 364, 3341–3350. [Google Scholar] [CrossRef]
  13. Sanker, C.; Lambertz, C.; Gauly, M. Climatic effects in Central Europe on the frequency of medical treatments of dairy cows. Animal 2013, 7, 316–321. [Google Scholar] [CrossRef] [Green Version]
  14. Vitali, A.; Segnalini, M.; Bertocchi, L.; Bernabucci, U.; Nardone, A.; Lacetera, N. Seasonal pattern of mortality and relationships between mortality and temperature-humidity index in dairy cows. J. Dairy Sci. 2009, 92, 3781–3790. [Google Scholar] [CrossRef] [Green Version]
  15. Houghton, J.T.; Ding, Y.; Griggs, D.J.; Noguer, M.; van der Linden, P.J.; Dai, X.; Maskell, K.; Johnson, C. Climate Change 2001: The Scientific Basis; The Press Syndicate of the University of Cambridge: Cambridge, UK, 2001. [Google Scholar]
  16. Herbut, P.; Angrecka, S.; Walczak, J. Environmental parameters to assessing of heat stress in dairy cattle—A review. Int. J. Biometeorol. 2018, 62, 2089–2097. [Google Scholar] [CrossRef] [PubMed]
  17. Rust, J.; Rust, T. Climate change and livestock production: A review with emphasis on Africa. South Afr. J. Anim. Sci. 2013, 43, 255–267. [Google Scholar] [CrossRef]
  18. Maurya, R. Alternate Dairy Management Practices in Draught Prone Areas of Bundelkhand Region of Uttar Pradesh. Master’s Thesis, IVRI, Izatnagar, India, 2010. [Google Scholar]
  19. Ministry of Finance Revenue and Economic Affairs. Islamabad. Pakistan Economic Survey. Available online: www.finance.gov.pk›survey_1617 (accessed on 13 July 2019).
  20. Gorst, A.; Dehlavi, A.; Groom, B. Crop productivity and adaptation to climate change in Pakistan. Environ. Dev. Econ. 2018, 23, 679–701. [Google Scholar] [CrossRef] [Green Version]
  21. Kreft, S.; Eckstein, D.; Melchior, I. Global Climate risk Index 2017: Who Suffers Most from Extreme Weather Events? Weather-Related Loss Events in 2015 and 1996 to 2015; Germanwatch Nord-Süd Initiative eV: Bonn, Germany, 2016. [Google Scholar]
  22. Government of Pakistan. Economic Survey of Pakistan 2018–19; Economic Advisor’s Wing, Finance Division, Government of Pakistan: Islamabad, Pakistan, 2018. Available online: http://www.finance.gov.pk/survey_1819 (accessed on 17 June 2019).
  23. Board, L.A.D.D. Dairy Production. Available online: www.lddb.org.pk (accessed on 15 July 2019).
  24. Naqvi, S.; Sejian, V. Global climate change: Role of livestock. Asian J. Agric. Sci. 2011, 3, 19–25. [Google Scholar]
  25. York, L.; Heffernan, C.; Rymer, C.; Panda, N. A deterministic evaluation of heat stress mitigation and feed cost under climate change within the smallholder dairy sector. Animal 2017, 11, 900–909. [Google Scholar] [CrossRef]
  26. Fahad, S.; Wang, J. Farmers’ risk perception, vulnerability, and adaptation to climate change in rural Pakistan. Land Use Policy 2018, 79, 301–309. [Google Scholar] [CrossRef]
  27. Weber, E.U. What shapes perceptions of climate change? Wiley Interdiscip. Rev. Clim. Chang. 2010, 1, 332–342. [Google Scholar] [CrossRef]
  28. Bohensky, E.; Kirono, D.; Butler, J.; Rochester, W.; Habibi, P.; Handayani, T.; Yanuartati, Y. Climate knowledge cultures: Stakeholder perspectives on change and adaptation in Nusa Tenggara Barat, Indonesia. Clim. Risk Manag. 2016, 12, 17–31. [Google Scholar] [CrossRef] [Green Version]
  29. Nguyen, T.P.L.; Seddaiu, G.; Virdis, S.G.P.; Tidore, C.; Pasqui, M.; Roggero, P.P. Perceiving to learn or learning to perceive? Understanding farmers’ perceptions and adaptation to climate uncertainties. Agric. Syst. 2016, 143, 205–216. [Google Scholar] [CrossRef]
  30. Arunrat, N.; Wang, C.; Pumijumnong, N.; Sereenonchai, S.; Cai, W. Farmers’ intention and decision to adapt to climate change: A case study in the Yom and Nan basins, Phichit province of Thailand. J. Clean. Prod. 2017, 143, 672–685. [Google Scholar] [CrossRef]
  31. Burnham, M.; Ma, Z. Climate change adaptation: Factors influencing Chinese smallholder farmers’ perceived self-efficacy and adaptation intent. Reg. Environ. Chang. 2017, 17, 171–186. [Google Scholar] [CrossRef]
  32. Sun, Y.; Han, Z. Climate change risk perception in Taiwan: Correlation with individual and societal factors. Int. J. Environ. Res. Public Health 2018, 15, 91. [Google Scholar] [CrossRef] [PubMed]
  33. Mercado, R.M. People’s risk perceptions and responses to climate change and natural disasters in BASECO compound, Manila, Philippines. Procedia Environ. Sci. 2016, 34, 490–505. [Google Scholar] [CrossRef]
  34. Mase, A.S.; Gramig, B.M.; Prokopy, L.S. Climate change beliefs, risk perceptions, and adaptation behavior among Midwestern US crop farmers. Clim. Risk Manag. 2017, 15, 8–17. [Google Scholar] [CrossRef]
  35. Rauf, S.; Bakhsh, K.; Abbas, A.; Hassan, S.; Ali, A.; Kächele, H. How hard they hit? Perception, adaptation and public health implications of heat waves in urban and peri-urban Pakistan. Environ. Sci. Pollut. Res. 2017, 24, 10630–10639. [Google Scholar] [CrossRef] [PubMed]
  36. Bakhsh, K.; Rauf, S.; Zulfiqar, F. Adaptation strategies for minimizing heat wave induced morbidity and its determinants. Sustain. Cities Soc. 2018, 41, 95–103. [Google Scholar] [CrossRef]
  37. Nardone, A.; Ronchi, B.; Lacetera, N.; Ranieri, M.S.; Bernabucci, U. Effects of climate changes on animal production and sustainability of livestock systems. Livest. Sci. 2010, 130, 57–69. [Google Scholar] [CrossRef]
  38. Sautier, M.; Piquet, M.; Duru, M.; Martin-Clouaire, R. Exploring adaptations to climate change with stakeholders: A participatory method to design grassland-based farming systems. J. Environ. Manag. 2017, 193, 541–550. [Google Scholar] [CrossRef]
  39. Abid, M.; Scheffran, J.; Schneider, U.; Ashfaq, M. Farmers’ perceptions of and adaptation strategies to climate change and their determinants: The case of Punjab province, Pakistan. Earth Syst. Dyn. 2015, 6, 225–243. [Google Scholar] [CrossRef]
  40. Pasha, H.A.; Growth of the Provincial Economies. Institute for Policy Reforms (IPR). 2015. Available online: http://ipr. org. pk/wp-content/uploads/2016/04/GROWTH-OF-PROVINCIAL-ECONOMICS-.Pdf (accessed on 12 March 2019).
  41. Ashfaq, M.; Razzaq, A.; Hassan, S. Factors affecting the economic losses due to livestock diseases: A case study of district Faisalabad. Pak. J. Agric. Sci. 2015, 52, 515–520. [Google Scholar]
  42. Razzaq, A.; Qing, P.; Abid, M.; Anwar, M.; Javed, I. Can the informal groundwater markets improve water use efficiency and equity? Evidence from a semi-arid region of Pakistan. Sci. Total Environ. 2019, 666, 849–857. [Google Scholar] [CrossRef] [PubMed]
  43. Mohammad, S. Supply Response of Major Crops in Different Agro-Ecologic, Zones in Punjab. Ph.D. Thesis, University of Agriculture: Faisalabad, Pakistan, 2005. [Google Scholar]
  44. Punjab Development Statistics (P.D.S). Bureau of Statistic, Government of Punjab. 2017. Available online: http://www.bos.gop.pk/developmentstat (accessed on 28 September 2019).
  45. McFadden, D. Modeling the Choice of Residential Location; Transportation Research Record: Washington, DC, USA, 1978. [Google Scholar]
  46. Greene, W.H. Econometric Analysis, 4th ed.; Prentice-Hall: Upper Saddle River, NJ, USA, 2003. [Google Scholar]
  47. Liddell, T.M.; Kruschke, J.K. Analyzing ordinal data with metric models: What could possibly go wrong? J. Exp. Soc. Psychol. 2018, 79, 328–348. [Google Scholar] [CrossRef] [Green Version]
  48. Ali, A.; Erenstein, O. Assessing farmer use of climate change adaptation practices and impacts on food security and poverty in Pakistan. Clim. Risk Manag. 2017, 16, 183–194. [Google Scholar] [CrossRef]
  49. Deressa, T.T.; Hassan, R.M.; Ringler, C.; Alemu, T.; Yesuf, M. Determinants of farmers’ choice of adaptation methods to climate change in the Nile Basin of Ethiopia. Glob. Environ. Chang. 2009, 19, 248–255. [Google Scholar] [CrossRef]
  50. Bryan, E.; Ringler, C.; Okoba, B.; Roncoli, C.; Silvestri, S.; Herrero, M. Adapting agriculture to climate change in Kenya: Household strategies and determinants. J. Environ. Manag. 2013, 114, 26–35. [Google Scholar] [CrossRef]
  51. Upadhyay, R.; Ashutosh, R.V.; Singh, S.; Aggarwal, P. Impact of Climate Change on reproductive functions of cattle and buffaloes. In Global Climate Change and Indian Agriculture; Aggarwal, P.K., Ed.; ICAR: New Delhi, India, 2009; pp. 107–110. [Google Scholar]
  52. StataCorp LP. Stata Multilevel Mixed-Effects Reference Manual; StataCorp LP: College Station, TX, USA, 2013. [Google Scholar]
  53. Hassan, R.M.; Nhemachena, C. Determinants of African farmers’ strategies for adapting to climate change: Multinomial choice analysis. Afr. J. Agric. Resour. Econ. 2008, 2, 83–104. [Google Scholar]
  54. Wollni, M.; Lee, D.R.; Thies, J.E. Conservation agriculture, organic marketing, and collective action in the Honduran hillsides. Agric. Econ. 2010, 41, 373–384. [Google Scholar] [CrossRef]
  55. Ullah, R.; Shivakoti, G.P.; Rehman, M.; Kamran, M.A. Catastrophic risks management at farm: The use of diversification, precautionary savings and agricultural credit. Pak. J. Agric. Sci. 2015, 52, 1139–1147. [Google Scholar]
  56. Ahmad, D.; Afzal, M.; Rauf, A. Analysis of wheat farmers’ risk perceptions and attitudes: Evidence from Punjab, Pakistan. Nat. Hazards 2019, 95, 845–861. [Google Scholar] [CrossRef]
  57. Cooper, P.; Dimes, J.; Rao, K.; Shapiro, B.; Shiferaw, B.; Twomlow, S. Coping better with current climatic variability in the rain-fed farming systems of sub-Saharan Africa: An essential first step in adapting to future climate change? Agric. Ecosyst. Environ. 2008, 126, 24–35. [Google Scholar] [CrossRef] [Green Version]
  58. Ma, L.; Chen, M.; Che, X.; Fang, F. Farmers’ Rural-To-Urban Migration, Influencing Factors and Development Framework: A Case Study of Sihe Village of Gansu, China. Int. J. Environ. Res. Public Health 2019, 16, 877. [Google Scholar] [CrossRef] [PubMed]
  59. Twongyirwe, R.; Mfitumukiza, D.; Barasa, B.; Naggayi, B.R.; Odongo, H.; Nyakato, V.; Mutoni, G. Perceived effects of drought on household food security in South-western Uganda: Coping responses and determinants. Weather Clim. Extrem. 2019, 24, 100201. [Google Scholar] [CrossRef]
  60. Udmale, P.D.; Ichikawa, Y.; Kiem, A.S.; Panda, S.N. Drought impacts and adaptation strategies for agriculture and rural livelihood in the Maharashtra State of India. Open Agric. J. 2014, 8, 41–47. [Google Scholar] [CrossRef]
  61. Akhtar, S.; LI, G.-C.; Nazir, A.; Razzaq, A.; Ullah, R.; Faisal, M.; Naseer, M.A.U.R.; Raza, M.H. Maize production under risk: The simultaneous adoption of off-farm income diversification and agricultural credit to manage risk. J. Integr. Agric. 2019, 18, 460–470. [Google Scholar] [CrossRef] [Green Version]
  62. Batima, P.; Natsagdorj, L.; Gombluudev, P.; Erdenetsetseg, B. Observed climate change in Mongolia. Assess Imp Adapt Clim Chang. Work Pap. 2005, 12, 1–26. [Google Scholar]
  63. Rota, A.; Sidahmed, A. IFAD’s Livestock Position Paper: Livestock Planning, Challenges and Strategies for Livestock Development in IFAD. IFAD: Rome, Italy, 2010. [Google Scholar]
  64. McCubbin, S.; Smit, B.; Pearce, T. Where does climate fit? Vulnerability to climate change in the context of multiple stressors in Funafuti, Tuvalu. Glob. Environ. Chang. 2015, 30, 43–55. [Google Scholar] [CrossRef]
  65. Mertz, O.; Mbow, C.; Reenberg, A.; Diouf, A. Farmers’ perceptions of climate change and agricultural adaptation strategies in rural Sahel. Environ. Manag. 2009, 43, 804–816. [Google Scholar] [CrossRef]
  66. Shrum, T.R.; Travis, W.R.; Williams, T.M.; Lih, E. Managing climate risks on the ranch with limited drought information. Clim. Risk Manag. 2018, 20, 11–26. [Google Scholar] [CrossRef]
  67. Ado, A.M.; Leshan, J.; Savadogo, P.; Bo, L.; Shah, A.A. Farmers’ awareness and perception of climate change impacts: Case study of Aguie district in Niger. Environ. Dev. Sustain. 2018, 1–15. [Google Scholar] [CrossRef]
  68. Ullah, W.; Nihei, T.; Nafees, M.; Zaman, R.; Ali, M. Understanding climate change vulnerability, adaptation and risk perceptions at household level in Khyber Pakhtunkhwa, Pakistan. Int. J. Clim. Chang. Strateg. Manag. 2018, 10, 359–378. [Google Scholar] [CrossRef]
  69. Akhtar, S.; LI, G.-C.; Ullah, R.; Nazir, A.; Iqbal, M.A.; Raza, M.H.; Iqbal, N.; Faisal, M. Factors influencing hybrid maize farmers’ risk attitudes and their perceptions in Punjab Province, Pakistan. J. Integr. Agric. 2018, 17, 1454–1462. [Google Scholar] [CrossRef]
  70. Bakhsh, K.; Sana, F.; Ahmad, N. Dengue fever in Punjab, Pakistan: Knowledge, perception and adaptation among urban adults. Sci. Total Environ. 2018, 644, 1304–1311. [Google Scholar] [CrossRef] [PubMed]
  71. Harrison, M.T.; Cullen, B.R.; Armstrong, D. Management options for dairy farms under climate change: Effects of intensification, adaptation and simplification on pastures, milk production and profitability. Agric. Syst. 2017, 155, 19–32. [Google Scholar] [CrossRef]
  72. Kais, S.M.; Islam, M.S. Perception of Climate Change in Shrimp-Farming Communities in Bangladesh: A Critical Assessment. Int. J. Environ. Res. Public Health 2019, 16, 672. [Google Scholar] [CrossRef] [PubMed]
  73. Saqib, S.E.; Ahmad, M.M.; Panezai, S.; Rana, I.A. An empirical assessment of farmers’ risk attitudes in flood-prone areas of Pakistan. Int. J. Disaster Risk Reduct. 2016, 18, 107–114. [Google Scholar] [CrossRef]
  74. Imran, M.; Shrestha, R.P.; Datta, A. Comparing farmers’ perceptions of climate change with meteorological data in three irrigated cropping zones of Punjab, Pakistan. Environ. Dev. Sustain. 2018, 1–20. [Google Scholar] [CrossRef]
  75. Siddiky, M.N.A. SAARC Dairy Outlook; SAARC Agriculture Centre: Dhaka, Bangladesh, 2015; pp. 1–160. [Google Scholar]
  76. Mandal, D.; SUBBA RAO, A.; SINGH, K.; Singh, S. Effects of macroclimatic factors on milk production in a frieswal herd@. Indian J. Dairy Sci. 2002, 55, 166–170. [Google Scholar]
  77. Younas, M.; Ishaq, K.; Ali, I. Effect of climate change on livestock production in Pakistan. In Proceedings of the International Seminar on Animal Industry, Jakarta, Indonesia, 5–6 July 2012. [Google Scholar]
  78. Morton, J.F. The impact of climate change on smallholder and subsistence agriculture. Proc. Natl. Acad. Sci. USA 2007, 104, 19680–19685. [Google Scholar] [CrossRef] [Green Version]
  79. Kumar, S.; Raju, B.; Ramarao, C.; Ramilan, T. Sensitivity of livestock production to climatic variability under indian drylands and future perspective. Curr. Agric. Res. J. 2015, 3, 142–149. [Google Scholar] [CrossRef]
  80. Hidosa, D.; Guyo, M. Climate Change Effects on Livestock Feed Resources: A Review. J. Fish. Livest Prod. 2017, 5, 259. [Google Scholar]
  81. Change, I.C. The Physical Science Basis, Contribution of Working Group I to the Fourth Assessment Report of the Intergovernmental Panel on Climate Change; IPCC: Geneva, Switzerland, 2007; p. 996. [Google Scholar]
  82. Lacetera, N. Impact of climate change on animal health and welfare. Anim. Front. 2018, 9, 26–31. [Google Scholar] [CrossRef] [Green Version]
  83. Grace, D.; Bett, B.K.; Lindahl, J.F.; Robinson, T.P. Climate and Livestock Disease: Assessing the Vulnerability of Agricultural Systems to Livestock Pests under Climate Change Scenarios; CCAFS Working Paper: Copenhagen, Denmark, 2015. [Google Scholar]
  84. National Disaster Management Authority Pakistan. Annual Report; National Disaster Management Authority: Islamabad, Pakistan, 2014.
  85. Hill, D.L.; Wall, E. Weather influences feed intake and feed efficiency in a temperate climate. J. Dairy Sci. 2017, 100, 2240–2257. [Google Scholar] [CrossRef] [PubMed]
  86. Mote, S.; Chauhan, D.; Ghosh, N. Effect of environment factors on milk production and lactation length under different seasons in crossbred cattle. Indian J. Anim. Res. 2016, 50, 175–180. [Google Scholar] [CrossRef]
  87. Kimaro, E.G.; Mor, S.M.; Toribio, J.-A.L. Climate change perception and impacts on cattle production in pastoral communities of northern Tanzania. Pastoralism 2018, 8, 19. [Google Scholar] [CrossRef] [Green Version]
  88. Dhakal, C.; Regmi, P.; Dhakal, I.; Khanal, B.; Bhatta, U.; Barsila, S.; Acharya, B. Perception, impact and adaptation to climate change: An analysis of livestock system in Nepal. J. Anim. Sci. Adv. 2013, 3, 462–471. [Google Scholar]
  89. Sharma, V.B.; Verma, M.R.; Qureshi, S.; Bharti, P. Effects of diseases on milk production and body weight of cattle in Uttar Pradesh. Int. J. Agric. Environ. Biotechnol. 2016, 9, 463–465. [Google Scholar] [CrossRef]
  90. Kasulo, V.; Chikagwa-Malunga, S.; Chagunda, M.; Roberts, D. The perceived impact of climate change and variability on smallholder dairy production in northern Malawi. Afr. J. Agric. Res. 2012, 7, 4830–4837. [Google Scholar] [CrossRef]
  91. Bett, B.; Kiunga, P.; Gachohi, J.; Sindato, C.; Mbotha, D.; Robinson, T.; Lindahl, J.; Grace, D. Effects of climate change on the occurrence and distribution of livestock diseases. Prev. Vet. Med. 2017, 137, 119–129. [Google Scholar] [CrossRef]
  92. Dorosh, P.; Malik, S.J.; Krausova, M. Rehabilitating agriculture and promoting food security after the 2010 Pakistan floods: Insights from the south Asian experience. Pak. Dev. Rev. 2010, 49, 167. [Google Scholar] [CrossRef]
  93. Barnes, A.; Islam, M.M.; Toma, L. Heterogeneity in climate change risk perception amongst dairy farmers: A latent class clustering analysis. Appl. Geogr. 2013, 41, 105–115. [Google Scholar] [CrossRef]
  94. Deressa, T.T.; Hassan, R.M.; Ringler, C.; Alemu, T.; Yesuf, M. Analysis of the Determinants of Farmers’ Choice of Adaptation Methods and Perceptions of Climate Change in the Nile Basin of Ethiopia [in Amharic]; International Food Policy Research Institute (IFPRI): Washington, DC, USA, 2008. [Google Scholar]
  95. Lasco, R.D.; Espaldon, M.L.O.; Habito, C.M.D. Smallholder farmers’ perceptions of climate change and the roles of trees and agroforestry in climate risk adaptation: Evidence from Bohol, Philippines. Agrofor. Syst. 2016, 90, 521–540. [Google Scholar] [CrossRef]
  96. Lee, T.M.; Markowitz, E.M.; Howe, P.D.; Ko, C.-Y.; Leiserowitz, A.A. Predictors of public climate change awareness and risk perception around the world. Nat. Clim. Chang. 2015, 5, 1014. [Google Scholar] [CrossRef]
  97. Shikuku, K.M.; Winowiecki, L.; Twyman, J.; Eitzinger, A.; Perez, J.G.; Mwongera, C.; Läderach, P. Smallholder farmers’ attitudes and determinants of adaptation to climate risks in East Africa. Clim. Risk Manag. 2017, 16, 234–245. [Google Scholar] [CrossRef]
  98. Nhemachena, C.; Mano, R. Assessment of the Economic Impacts of Climate Change on Agriculture in Zimbabwe: A Ricardian Approach; The World Bank: Washington, DC, USA, 2007. [Google Scholar]
  99. Ochenje, I.; Ritho, C.; Guthiga, P.; Mbatia, O. Assessment of farmers’ perception to the effects of climate change on water resources at farm level: The case of Kakamega county, Kenya. In Proceedings of the 2016 Fifth International Conference, Addis Ababa, Ethiopia, 23–26 September 2016. [Google Scholar]
  100. Nakano, Y.; Tsusaka, T.W.; Aida, T.; Pede, V.O. Is farmer-to-farmer extension effective? The impact of training on technology adoption and rice farming productivity in Tanzania. World Dev. 2018, 105, 336–351. [Google Scholar] [CrossRef] [Green Version]
  101. Mustafa, G.; Latif, I.A.; Bashir, M.K.; Shamsudin, M.N.; Daud, W.M.N.W. Determinants of farmers’ awareness of climate change. Appl. Environ. Educ. Commun. 2019, 18, 219–233. [Google Scholar] [CrossRef]
  102. Rondhi, M.; Fatikhul Khasan, A.; Mori, Y.; Kondo, T. Assessing the Role of the Perceived Impact of Climate Change on National Adaptation Policy: The Case of Rice Farming in Indonesia. Land 2019, 8, 81. [Google Scholar] [CrossRef]
  103. Pakistan Meteorological Department. Climate Change Scenerios Data for Pakistan. Available online: http://www.pmd.gov.pk (accessed on 30 September 2019).
Figure 1. Dairy farming context under climatic risks.
Figure 1. Dairy farming context under climatic risks.
Ijerph 16 04036 g001
Figure 2. Map of study area.
Figure 2. Map of study area.
Ijerph 16 04036 g002
Figure 3. Sampling framework.
Figure 3. Sampling framework.
Ijerph 16 04036 g003
Figure 4. Climate-Related risks perceived by the dairy farmers.
Figure 4. Climate-Related risks perceived by the dairy farmers.
Ijerph 16 04036 g004
Figure 5. Climatic variability perceived by the dairy farmers.
Figure 5. Climatic variability perceived by the dairy farmers.
Ijerph 16 04036 g005
Figure 6. Historical trend of climate change in Punjab, Pakistan Source: Pakistan Meteorological Department (PMD) Dataset.
Figure 6. Historical trend of climate change in Punjab, Pakistan Source: Pakistan Meteorological Department (PMD) Dataset.
Ijerph 16 04036 g006
Figure 7. Perceived impacts of climatic risks on dairy farming system.
Figure 7. Perceived impacts of climatic risks on dairy farming system.
Ijerph 16 04036 g007
Figure 8. Risk coping adaptation strategies practiced by dairy farmers.
Figure 8. Risk coping adaptation strategies practiced by dairy farmers.
Ijerph 16 04036 g008
Table 1. Number of Dairy Animals in the Study Area (millions).
Table 1. Number of Dairy Animals in the Study Area (millions).
TerritoryCattleBuffaloesSheepGoatsMajor Crops
Pakistan47.840.030.976.1
Punjab13.20416.0194.94217.392
Muzaffar-Garh
(from irrigated zone of Punjab)
1.10.90.51.3Cotton, wheat, maize, sugarcane
Jhang
(mixed cropping zone of Punjab)
0.91.20.41.0Sugarcane, cotton, wheat, maize, tobacco
Jhelum
(rain-fed zone of Punjab)
0.10.090.0130.13Wheat, rice, fruit, vegetables, fodder
Source: Pakistan Economic Survey 2018–19 [22]; Punjab Development Statistics 2017 [44].
Table 2. Summary statistics of surveyed dairy farmers (n = 450).
Table 2. Summary statistics of surveyed dairy farmers (n = 450).
VariablesMean ± SD or n (%)
Socio-economic characteristics
Age (years)45.1 ± 12.6
Dairy Farming Experience (years)21.2 ± 11.7
Size of Family (no. of heads)8.5 ± 3.6
Share of Dairy Income (%)42.8 ± 15.5
Educational level (years)
No Education88 (19.6)
Primary School137 (30.4)
High School174 (38.7)
College/University51 (11.3)
Farm characteristics
Farm Size (ha)2.9 ± 1.2
Land Allocated to Dairy Animals (ha)0.2 ± 0.01
No. of Milking Animals (no. of heads)4.5 ± 2.6
Milk Production/per day (liter)14.1 ± 6.5
Farm type (in number)
Irrigated234 (52)
Rain-fed210 (46.7)
Mixed6 (1.3)
Breed of dairy animals
Indigenous172 (38.2)
Cross49 (10.9)
Mixed229 (50.9)
Table 3. Model results with marginal effects.
Table 3. Model results with marginal effects.
VariablesCoefficientsMarginal Effects
Prob (Y = 1|X)Prob (Y = 2|X)Prob (Y = 3|X)Prob (Y = 4|X)
dY/dXdY/dXdY/dXdY/dX
Climatic Parameters
Drought0.1404 (0.064) **−0.014(0.007) **−0.040(0.018) **0.021(0.010) **0.035(0.014) **
Flood0.1461 (0.063) **−0.014(0.006) **−0.041(0.017) **0.022(0.001) **0.031(0.013) **
Heat Waves0.1515 (0.065) *−0.015(0.006) **−0.015(0.006) **0.023(0.011) **0.023(0.011) **
Humidity0.1433 (0.066) **−0.014(0.007) **−0.039(0.018) **0.021(0.010) **0.031(0.014) **
Pest & Diseases0.1530 (0.063) **−0.015(0.006) **−0.043(0.016) **0.023(0.011) **0.032(0.013) **
Heavy Rainfall0.0994 (0.064) ns−0.001(0.006) ns−0.028(0.017) ns0.015(0.001) ns0.021(0.014) ns
Farmers’ Characteristics
Education0.0275 (0.013) **−0.030(0.010) **−0.080(0.040) **0.040(0.020) **0.060(0.030) **
Age0.0113 (0.007) ns0.0113 (0.007) ns−0.003(0.002) ns0.002(0.001) ns0.002(0.002) ns
Dairy Farming Experience0.0201 (0.008) **−0.002(0.000) **−0.005(0.002) **0.008(0.004) **0.009(0.005) **
Farm Characteristics
Milking Animals0.026 (0.014) ***−0.003 (0.001) ***−0.007 (0.004) ***0.048 (0.024) ***0.060 (0.030) ***
Breed of Animals0.122 (0.080) ns−0.012 (0.008) ns−0.034 (0.022) ns0.020 (0.012) ns0.030 (0.020) ns
Institutional Link
Social Participation0.055 (0.072) ns−0.005 (0.007) ns−0.015 (0.020) ns0.008 (0.010) ns0.012 (0.012) ns
Source of Information0.1061(0.046) **−0.010 (0.004) **−0.030 (0.013) **0.016 (0.007) **0.023 (0.001) **
Contact with Extension0.232 (0.080) *−0.022 (0.008) ***−0.064 (0.023) *0.034 (0.013) *0.050 (0.017) *
Services
μ1 2.4518 (0.494) ***
μ2 3.7753 (0.496) ***
μ3 5.2023 (0.518) **
μ4 6.8030 (0.576) **
Observations 450
LR chi2(9) 94.90
Prob > chi2 0.0000
Log likelihood −512.1638
Standard Errors are given in parentheses; *, **, *** are 1%, 5% and 10% level of significance, ‘ns’ indicates not significant.

Share and Cite

MDPI and ACS Style

Abbas, Q.; Han, J.; Adeel, A.; Ullah, R. Dairy Production under Climatic Risks: Perception, Perceived Impacts and Adaptations in Punjab, Pakistan. Int. J. Environ. Res. Public Health 2019, 16, 4036. https://0-doi-org.brum.beds.ac.uk/10.3390/ijerph16204036

AMA Style

Abbas Q, Han J, Adeel A, Ullah R. Dairy Production under Climatic Risks: Perception, Perceived Impacts and Adaptations in Punjab, Pakistan. International Journal of Environmental Research and Public Health. 2019; 16(20):4036. https://0-doi-org.brum.beds.ac.uk/10.3390/ijerph16204036

Chicago/Turabian Style

Abbas, Qasir, Jiqin Han, Adnan Adeel, and Raza Ullah. 2019. "Dairy Production under Climatic Risks: Perception, Perceived Impacts and Adaptations in Punjab, Pakistan" International Journal of Environmental Research and Public Health 16, no. 20: 4036. https://0-doi-org.brum.beds.ac.uk/10.3390/ijerph16204036

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop