Meteorological impacts on evapotranspiration in different climatic zones of Pakistan
Shahzada ADNAN1,2,*, Kalim ULLAH1, Azmat H KHAN2, GAO Shouting3
1 Department of Meteorology, COMSATS Institute of Information Technology, Islamabad 44000, Pakistan
2 Pakistan Meteorological Department, Islamabad 44000, Pakistan
3 State key Laboratory of Cloud-Precipitation Physics and Severe Storms, Institute of Atmospheric Physics, Chinese Academy of Sciences, Beijing 100029, China
*Corresponding author: Shahzada ADNAN (E-mail: shaz_adnan@pmd.gov.pk)
Abstract

Arid regions are highly vulnerable and sensitive to drought. The crops cultivated in arid zones are at high risk due to the high evapotranspiration and water demands. This study analyzed the changes in seasonal and annual evapotranspiration (ET) during 1951-2016 at 50 meteorological stations located in the extremely arid, arid, and semi-arid zones of Pakistan using the Penman Monteith (PM) method. The results show that ET is highly sensitive and positively correlated to temperature, solar radiation, and wind speed whereas vapor pressure is negatively correlated to ET. The study also identifies the relationship of ET with the meteorological parameters in different climatic zones of Pakistan. The significant trend analysis of precipitation and temperature (maximum and minimum) are conducted at 95% confidence level to determine the behaviors of these parameters in the extremely arid, arid, and semi-arid zones. The mean annual precipitation and annual mean maximum temperature significantly increased by 0.828 mm/a and 0.014°C/a in the arid and extremely arid zones, respectively. The annual mean minimum temperature increased by 0.017°C/a in the extremely arid zone and 0.019°C/a in the arid zone, whereas a significant decrease of 0.007°C/a was observed in the semi-arid zone. This study provides probabilistic future scenarios that would be helpful for policy-makers, agriculturists to plan effective irrigation measures towards the sustainable development in Pakistan.

Key words: evapotranspiration; meteorological parameters; climatic zone; Penman Monteith method; Pakistan
1 Introduction

The evapotranspiration (ET) is one of the most important climatic parameter determining the crop water requirements, water balance, and climatic features of a region. It is difficult to directly measure ET and models have been therefore widely used to calculate it (Gad and El-Gayar, 2010). The ET is generally affected by four major weather parameters including temperature, solar radiation, vapor pressure/humidity, and wind speed (Isikwue et al., 2014). The chaotic climatic conditions have significant impacts on hydrological factors like evapotranspiration, soil moisture, ground water, and runoff (Bultot et al., 1988). The global annual mean temperature is rising and significant warming trends have been observed spatially and in terms of intensity (Hansen et al., 2012). The increase of greenhouse gas concentration is considered as one of the major reasons of global warming (IPCC, 2014). Studies identified the spatial and temporal effects of global warming on the water resources and hydrological cycle around the world (Iglesias et al., 1994; Feddema, 1999; Held and Soden, 2006; IPCC, 2007, 2014). The precipitation and ET are the main components of water cycle to quantify the available water (Sabziparvar and Tabari, 2010).

Increase in temperature can enhance the water vapor holding capacity of air and therefore affects evapotranspiration. Moreover, the saturation pressure increases exponentially with the increase in temperature as well. By keeping all the other meteorological factors constant, the increase of temperature can cause dryness and enhance ET. Furthermore, the ET also increases with the increase in solar radiation and decreases in cloudiness and vapor pressure (Gong et al., 2006). Dry and warm wind increases ET, while moist and cold wind decreases ET. Hot and dry wind blows in summer, while cold and dry wind is observed in winter in most of the arid regions of Central South Asia, unless the monsoon and western disturbances (winter) brings rainfall to the region (Adnan et al., 2016). The ET is less sensitive to temperature in winter than in summer. Increase in wind speed enhances the concentration of dust particles in the atmosphere which lowers the incoming solar radiation. The sensitivity of ET depends upon the climatic factors, including temperature, sunshine hours, wind speed, especially humidity and the characteristics of plants and resistance of their stomata in different climatic zones (Mamassis et al., 2014). Importantly, both high temperature and low humidity increase the ET, which exacerbate the drought condition in both arid and humid regions (Le Houerou, 1993; Smakhtin and Schipper, 2008). Moreover, the increase in evapotranspiration rate generates high irrigation water demand. In any country having both arid and humid regions, increasing evapotranspiration rate may lead to a severe shortage of irrigation water. For Pakistan, a developing country with economy mainly depending on irrigated agriculture, the problem is even more important. The climate of Pakistan is diversified as that two-thirds of total area is in arid and semi-arid areas, and 19% of the area is dominated by humid climate (Chaudhary and Rasul, 2003; Adnan, 2009; Haider and Adnan, 2014).

A proper water management policy is required in Pakistan to avoid severe water shortage in the near future. The agricultural activities are carried out by proper irrigation in the extremely arid and arid regions whereas the supplementary irrigation is required in the semi-arid and humid regions. The ET remains very high in most time of a year in Pakistan. However, monsoon brings rainfall to ease the moisture stress in these regions and helpful to fulfill the crop water demand to some extent. Therefore, comprehensive information is required about the future climatic conditions in view of changing climatic scenarios. Kruss et al. (1992) identified a cooling trend due to increase in cloudiness and rainfall over northern and southeastern Pakistan during monsoon season. Singh and Sontakke (1996) observed large variability in temperature in the country by employing a long-term reconstructed time series temperature data from 1876 to 1993. Moreover, a warming trend was also detected with a total change of 0.2° C at the beginning of the 20th century. Afzaal et al. (2009) observed 0.06° C/10a increase in temperature with a total change of 0.64° C in last 127 years, showing the warming trend in Pakistan. This study is motivated by the above mentioned considerations to assess sensitivity of the evapotranspiration due to ± 20% change in several meteorological parameters.

The variation in ET directly influences the crop water requirements and irrigation schedule. The crops in arid regions always remain at risk due to high evapotranspiration and crop water demands. It is necessary to evaluate the impacts of meteorological parameters on the sensitivity of ET against the global warming in arid regions of Pakistan. This study provides a range of changes in each of the parameters, i.e., temperature, solar radiation, vapor pressure, and wind speed that likely to affect both in seasonal and annual ET. Furthermore, the results of annual future trends analysis (temperature and precipitation) will help to understand the drought conditions, ET demands, and water balance in extremely arid, arid, and semi-arid regions of Pakistan. This study will be helpful for key stakeholders i.e., agriculturist, agronomist disaster agencies to develop future contingency plan for sustainable development towards agriculture in different climatic zones in Pakistan.

2 Materials and methods
2.1 Study area

Pakistan is located in South Asia (23° 39′ -37° 01′ N, 60° 49′ -77° 40′ E) with a total area of 0.796× 106km2 (Fig. 1). The mean annual precipitation ranges from 20 to 800 mm. The annual mean maximum temperature is 27.5° C-34.3° C and the annual mean minimum temperature is 12.3° C-17.3° C. Two main weather systems, southwestern monsoon during summer and western disturbance system during winter, approach Pakistan and contribute about 45% and 20% of annual precipitation during July to September (monsoon) and December to February (winter), respectively. July is the month with the maximum rainfall and June is the hottest month of a year. There are considerable precipitation falling in the crop seasons of Kharif from May to September and Rabi from October to the succeeding April (Adnan and Khan, 2009).

Pakistan can be divided into extremely arid, arid, semi-arid, and humid climatic zones, which accounted for 25%, 25%, 31%, and 19% of the total area of the country, respectively (Adnan et al., 2017). The two-thirds of the country is dominated by the arid climate (Haider and Adnan, 2014), most of which lies in southern Pakistan where the topography is characterized by a slow-lying plains. The elevation of western Pakistan is below 1500 m with barren mountains. This research selected the extremely arid, arid, and semi-arid zones to analyze the sensitivity of ET against the meteorological parameters.

According to the characteristics of precipitation and temperature, Shamshad (1988), Khan et al. (2010), and Zahid and Rasool (2011) divided a year into five seasons including spring (March and April), dry summer (May and June), monsoon season (wet summer; July to September), autumn (October and November), and winter (December to the following February) in Pakistan. The summer starts from May and ends in September consisting of dry and wet months. The temperature is very high in the dry months and drops down during wet months. The monsoon season is very wet and contributes 44% of the annual precipitation.

Fig. 1 Location of Pakistan and distributions of climatic zones and meteorological stations in Pakistan

2.2 Data

The data used in this research were obtained from Pakistan Meteorological Department. Observational monthly temperature, solar radiation, humidity, and wind speed from 50 stations (Fig. 1) for the period 1981-2010 were used for climatic zone classification. A long-term (1951-2016) data series of precipitation and temperature (maximum, minimum) was used to determine the time series trends (95% confidence level) on annual basis.

2.3 Methods

Martin et al. (1989) demonstrated the impacts of temperature, solar radiation, humidity, and wind speed on evapotranspiration in this study area. Liu et al. (2010) also pointed out that ET changes are driven by the combination of these climatic variables. The behavior of these parameters is quite different in each month and season and they influence the evapotranspiration in different manners.

The solar radiation and humidity data are used to calculate the net solar radiation and vapor pressure. Penman Monteith (PM) method, defined by Food and Agriculture Organization of the United Nations (FAO), is used to calculate the annual and seasonal ET. The sensitivity of ET is identified by meteorological parameters within a possible range of ± 20%, which is an approach applied in a dry region of India by Goyal (2004). The IPCC (2007) report suggests that in case of carbon dioxide emission doubled in the 21st century, the average global temperature would rise from 2.0° C to 4.5° C with the best assessment of 3.0° C, but it would not be uniform throughout the world. Future scenarios are generated by these parameters to identify the meteorological impacts on ET.

FAO-PM method is considered to be the best model for estimating ET in any climatic zones (Monteith, 1965; Allen et al., 1998) although several models have been used for calculating reference evapotranspiration (ET0) by Doorenbos and Wilson (1974), Wilson (1974), Pruitt (1977), Viessman et al. (1977), Eagleson (1978), Subramanya (1984), Mavi (1986), Michael (1986); and PM method was used to determine the impact of meteorological parameters on ET in sub-tropical semi-humid region of Bangladesh (Ali et al., 2009) and in a desert transition zone of Heihe River Basin in northwestern China (Zhao and Zhao, 2014). The modified FAO-PM method (Allen et al., 1998) is as Equation 1.

$ET0=\frac{0.408\Delta ({{R}_{n}}-G)+\gamma \left[ {}^{900}/{}_{\left( T +273 \right)} \right]{{u}_{2}}\left( {{e}_{s}}-{{e}_{a}} \right)}{\Delta +\gamma \left( 1.0.34{{u}_{2}} \right)}$, (1)

where ET0represents reference evapotranspiration (mm/day); ∆ is the slope of the saturation vapor pressure at air temperature T (kPa/° C); G is the soil heat flux density (MJ/(d• m2)); Rn is the net solar radiation (MJ/(d• m2)); T is the mean daily air temperature (° C); u2is the average daily wind speed (m/s); es and ea are vapor pressure deficits (kPa); and γ is the psychrometric constant (kPa/° C).

The water vapor pressure is directly associated with water quantity in the air and can be measured by Equation 2.

${{e}_{s}}=\left[ \frac{{{e}^{o}}\left( {{T }_{\max }} \right)+{{e}^{o}}\left( {{T }_{\min }} \right)}{2} \right]$, (2)

${{e}^{o}}\left( T \right)=0.6108\left[ \frac{17.27T }{T +237.3} \right]$, (3)

$\Delta =4098\left[ \frac{{{0.6108}_{\exp }}{\left( 17.27T \right)}/{\left( T +237.3 \right)}\; }{{{\left( T +237.3 \right)}^{2}}} \right]$, (4)

where es is the mean saturation vapor pressure (kPa); eo is the saturation vapor pressure (kPa); Tmax is the maximum temperature (° C); and Tmin is the minimum temperature (° C).

The mean air temperature is used to measure the slope curve of saturation vapor pressure. The pressure exerted by the water in air is called actual vapor pressure and is denoted by ea (Eq. 5).

${{e}_{a}}=\left[ \frac{{{e}^{o}}\left( {{T }_{\min }} \right)\left( {RH_{\max }}/{100}\; \right)+{{e}^{o}}\left( {{T }_{\max }} \right)\left( {R{H }_{\min }}/{100}\; \right)}{2} \right]$, (5)

where, ea is the actual vapor pressure; eo(Tmin) and eo(Tmax) are the saturation vapor pressures (kPa) at daily minimum and maximum temperatures (° C); and RHmax and RHmin are the maximum and minimum relative humidities.

The net solar radiation Rn is the difference between the incoming solar radiation (Rns) and outgoing solar radiation (Rnl; Eqs. 6 and 7).

${{R}_{n}}={{R}_{ns}}-{{R}_{nl}}$, (6)

${{R}_{ns}}=(1-\alpha ){{R}_{s}}$, (7)

where α is the surface albedo and Rs is the incoming shortwave solar radiation (MJ/(d• m2)) and can be calculated by the Angstrom’ s formula (Eq. 8; Doorenbos and Pruitt, 1977).

${{R}_{s}}=(0.25+0.50\frac{n}{N}){{R}_{a}}$, (8)

where n/N is the ratio of observed sunshine hours and the maximum possible sunshine hours; Ra is the extraterrestrial radiation ((MJ/(d• m2)); Allen et al., 1998). The net long wave radiation can be expressed by Equation 9.

$Rnl=\sigma \left[ \frac{Tmax{{K}^{4}}+T\min {{K}^{4}}}{2} \right](0.34-0.14\sqrt{ea})\left\{ 1.35\frac{Rs}{Rso}-0.35 \right\}$, (9)

where Rnl is the net long-wave radiation (MJ/(d• m2)); Tmax and Tmin are the maximum and minimum absolute temperatures (° C); K is the absolute mean temperature ((° C); K=° C+273.16 (Allen et. al., 1998)); σ is the Stefan Boltzmann constant; Rso is the clear sky shortwave solar radiation (MJ/(d• m2)).

The sensitivity of ET is calculated by changing one meteorological parameter between the ranges of ± 20% with keeping the other three parameters constant. This method will determine the percentage change in ET (increase/decrease) with respect to the changes in temperature, vapor pressure, net radiations, and wind speed during each season in each climatic zone. Similarly, a hypothetical future scenario (± 10% change) is generated to estimate the total demand of ET by taking consideration of various combinations of meteorological parameters. The percentage change in ET is calculated over the normal period 1981-2010 for each combination of the meteorological parameters. The annual time series (1951-2016) analysis of temperature (maximum and minimum) and precipitation are conducted to identify the future trends. The Mann-Kendall and Sen’ s Slope methods are used to identify the trend and trend magnitude at 95% confidence level, respectively (Mann, 1945; Sen, 1968; Kendal, 1975). The same methods were used to calculate the significant trend and related magnitude for South Central Asia (Adnan et al., 2016).

3 Results and discussion
3.1 Historical data analysis

The impact of meteorological parameters on ET is not easy to predict. The probability of scenario for different seasons and the respective sensitivity of ET demands expected under such assumed futuristic meteorological conditions are estimated in different climatic zones of Pakistan. The seasonal impacts on total ET with respect to the changes in meteorological parameters for the extremely arid region are identified (Table 1). The results show that ET has seasonal variations of 9.75-111.05, 7.64-76.38, 8.31-77.00, 4.47-34.74 mm and annual variations of 329.59, 236.11, 212.99, 125.7 mm due to the changes (± 20%) in temperature, net radiation, vapors pressure, and wind speed, respectively. The maximum changes in ET are observed in dry summer, monsoon, autumn, and winter; and the minimum changes in ET are observed in spring.

Figures 2a-d illustrate seasonal and annual relative percentage changes in ET due to the relative changes in temperature, vapor pressure, net radiations and wind speed in the extremely arid region. The percentage variation of ET ranged from -15.78% to 16.88% annually in responses to ± 20% changes in temperature. The minimum effect in temperature on ET is observed during spring (-14.42% to 15.05%) and the maximum during monsoon season (-16.80% to 18.17%).

In the extremely arid zone, the percentage changes (± 20%) of net solar radiation caused annual variations of -12.09% to 12.03% in ET. The maximum effects of net solar radiation on ET, with percentage variations of -13.21% to 13.25% were observed in autumn and the minimum appeared in summer (-11.22% to 11.24%). The annual percentage change in ET was ± 10.8% due to the changes in vapor pressure, with the maximum percentage change in winter (-13.0% to 13.2%) and the minimum appeared in spring (-9.1% to 9.3%). The ± 20% changes in wind speed caused annual percentage changes of -6.44% to 6.09% in ET, with the maximum percentage changes in summer (-7.04% to 6.66%) and the minimum in monsoon season (-5.69% to 5.36%).

Table 1 Impacts of seasonal percentage changes in meteorological parameters on evaportranspiration (ET) in the extremely arid region of Pakistan

In the arid zone, the ± 20% changes in temperature, net radiation, vapor pressure, and wind speed caused seasonal changes in ET of 5.90-90.83, 6.63-71.32, 5.02-53.52, and 2.79-100.62 mm and annual changes of 245.21, 208.79, 142.75, and 100.62 mm (Table 2), respectively. The maximum changes in ET appeared in dry summer, monsoon, autumn, and winter; and the minimum appeared in spring.

Fig. 2 Impacts of percentage changes in temperature (a), net radiations (b), vapor pressure (c), and wind speed (d) on evapotranspiration (ET) in the extremely arid region of Pakistan. W, winter; S, spring; Su, dry summer; M, monsoon; At, autumn; An, annual.

Table 2 Impacts of seasonal percentage changes in meteorological parameters on total ET in the arid zone of Pakistan

The impacts of the percentage changes in meteorological parameters on ET are presented in Figure 3. The ± 20% change temperature caused seasonal percentage changes of -13.98% to 14.65% in ET (Fig. 3a), with the smallest influence in winter (-11.49% to 12.34%) due to the cold climate. The maximum percentage changes of ET appeared in the monsoon season (-15.24% to 16.17%) as the soil moisture was lost under the high temperature.

The percentage changes of net solar radiation (± 20%) caused annual percentage changes in ET of -12.46% to 12.48%, with the maximum effects on ET (-14.15% to 14.62%) observed in autumn, and the minimum (-10.97% to 11.00%) appeared in summer. The percentage changes in vapor pressure result in an annual percentage change of 8.5% in ET, with the maximum of 10.2% appeared in winter and the minimum (-6.7% to 7.1%) in autumn. The percentage changes (± 20%) in wind speed caused an annual percentage changes in ET of -6.01% to 5.73%, with the maximum percentage changes in ET being -7.03% to 6.67% appeared in summer and the minimum (-5.10% to 4.94%) in autumn. The maximum increase in ET is caused by the hot and dry climatology where dry wind causes more ET.

Fig. 3 Impacts of percentage changes in temperature (a), net radiations (b), vapor pressure (c), and wind speed (d) on ET in the arid zone of Pakistan

In the semi-arid zone of Pakistan, the ± 20% changes in temperature, net radiation, vapor pressure and wind speed caused seasonal changes in ET of 4.30-72.23, 5.08-70.21, 1.87-29.16, 2.11-22.62 mm and annual variations of 189.60, 189.66, 77.68, 72.39 mm (Table 3), respectively. The maximum changes caused by temperature and net radiation appeared in monsoon season while the maximum changes in ET caused by the vapor pressure and wind speed appeared in winter and dry summer, respectively. The minimum changes in ET caused by temperature and net radiation occurred in winter while the minimum changes caused by vapor pressure and wind speed occurred in dry summer and monsoon season, respectively.

The seasonal and annual percentage changes in ET caused by the percentage changes of the meteorological parameters in the semi-arid zone of Pakistan are shown in Figures 4. The annual ET had percentage changes of -12.92% to 13.56% under the ± 20% changes in temperature. The temperature has less effects on ET (-10.67% to 11.66%) during winter. The maximum percentage changes in ET (-14.24% to 15.17%) is observed in the monsoon season as abundant moisture of the season and the high ET in this area. The percentage change (± 20%) in net solar radiation caused an annual percentage changes in ET of -13.50% to 13.53%, with the maximum percentage changes (-14.75% to 14.66%) in ET observed in monsoon season and the minimum (-12.16% to 12.26%) in summer. The result also shows that the net solar radiation has significant affects in the semi-arid zone as the vegetative cover, moisture convergence and cloud cover enhancing or reducing the ET, which agree with Claussen (1997).

The ± 20% changes of vapor pressure caused annual percentage changes of -5.6% to 5.3 % in ET, with the maximum changes in winter (-10.3% to 10.1%) and the minimum in summer (-2.8% to 2.6%). The ± 20% changes of wind speed caused annual percentage changes of -5.18% to 4.96% on ET, with the maximum percentage changes in summer (-6.10% to 5.84%) and the minimum in monsoon (-4.06% to 3.91%).

Table 3 Impact of seasonal percentage changes in meteorological parameters on total ET in the semi-arid zone of Pakistan

The results show that ET is highly sensitive to temperature and vapor pressure in the study area in monsoon season and winter, respectively. The solar radiation influences ET more in the extremely arid and arid zones in autumn and the semi-arid zone during monsoon season. Wind speed affects ET more in summer. ET is less sensitive to temperature decrease in comparison to the increase temperature that makes temperature a key triggering factor among these parameters irrespective of season (Xu et al., 2006). The high wind speed is associated with low humidity whereas the low wind speed is associated with high humidity. ET has fewer changes in dry summer when the wind is hot and dry (Allen et al., 1998). The wind speed has the minimum effects on annual ET and the maximum effects on ET appeared in winter because dry and cool winds cause more ET from the ground surface (Subramanya, 2008).

Based on seasonal and annual analysis, we measured the sensitivity of ET to temperature, solar radiation, vapor pressure and wind speed. The results show that ET is most sensitive to temperature in monsoon season and least sensitive to temperature in winter in the study area. The temperature is the main triggering factor enhancing the ET regardless the season. The changes of ET in monsoon season depict that abundant of moisture goes into atmosphere during the process of ET. The temperature and solar radiation cause the maximum change in ET in all the extremely arid, arid, and semi-arid zones during monsoon season. The maximum changes in ET caused by the winter vapor pressure and summer wind speed are observed in all the zones as well.

Fig. 4 Impacts of percentage change in temperature (a), net radiations (b), vapor pressure (c), and wind speed (d) on ET in the semi-arid zone of Pakistan

3.2 Scenario analysis

A scenario was developed to measure how ET will respond to the changes of the meteorological parameters in the extremely arid, arid, and semi-arid zones based on observational data (1981-2010). The scenario set the temperature at a 10% increase and the other three parameters, vapor pressure, wind speed, and net solar radiation, with ± 10% fluctuation (Table 4). Fourteen cases generated in the scenario are presented in Table 4. In Case 9, the ET in extremely arid zone varied from a marginal increase of 0.27% (-10% Rn, 10% ea, and -10% u2) to the maximum of 29.93% (10% Rn, -10% ea, and 10% u2). ET slightly decreased by 1.51% and 2.71%, and increased by 25.97% and 24.68% in ET was observed in the arid and semi-arid zones in Case 9 and 12, respectively. The reduction in ET is due to the 10% decrease of net solar radiation and wind speed, and 10% increase of vapor pressure. The results also depict that increasing in net solar radiation and wind speed by 10% over normal period with 10% decrease actual vapor pressure could increase ET in the entire study area. In different climatic zones, ET could increase by 4.19% to 9.23% with increasing in temperature by 10% and decreasing in the other parameters by 10%. These results are based on the assumption of percentage changes in meteorological parameters and their sensitivities to ET. The changes in ET can have profound effects on agriculture, water, and forest resources in Pakistan.

The results depict that the maximum increases in annual ET are observed under 10% increase of solar radiation and wind speed, and 10% decrease of vapor pressure (Case 9) in the study area. The maximum decreases in annual ET are observed under 10% decrease of solar radiation and wind speed, and 10% increase of vapor pressure but the changes are marginal. So, we may infer that increase in solar radiation and wind speed and decrease in vapor pressure are triggering factors to accelerate the ET process in the study area irrespective of the climatology. Furthermore, the results show that ET is directly proportional to the temperature, net solar radiation, and wind speed; and inversely proportional to the vapor pressure.

3.3 Trend analysis

The mean annual precipitation data calculated for the extremely arid, arid and semi-arid zones are 150, 250, and 550 mm, respectively. The annual precipitation shows a decreasing trend in the extremely arid zone and increasing trend in both the arid and semi-arid zones (Fig. 5). However, the significant increasing trend (0.828 mm/a) of annual precipitation is observed only in the arid zone at 95% confidence level. The annual precipitation variability in the arid zone is higher than those in the extremely arid and semi-arid zones. The increase in annual precipitation may help to reduce the moisture and water stresses in the arid zone and would benefit the agricultural activities.

Table 4 Annual changes in total estimated ET (ET0) with respect to change in meteorological parameters in different climatic zones of Pakistan

Fig. 5 Time series trend of annual precipitation in different climatic zones of Pakistan during 1951-2016. Zmk, M-K z-value; Q, Sen’ s slope estimator.

The annual mean maximum temperatures are 33.25° C, 30.11° C, and 26.41° C in the extremely arid, arid, and semi-arid zones, respectively. The mean annual maximum temperature shows an increasing trend over the entire study area (Fig. 6). The significant increasing trend (0.014° C/a) in annual mean maximum temperature (Tmax) is observed only in the extremely arid zone at 95% confidence level, which will enhance ET and aridity in the extremely arid zone and put pressures on water resources and agricultural activities.

Fig. 6 Time series trend of annual mean maximum temperature in different climatic zones of Pakistan during 1951-2016

The annual mean minimum temperatures are 18.37° C, 15.53° C, and 12.95° C in the extremely arid, arid and semi-arid zones, respectively. Likewise the maximum temperature, the long-term mean minimum temperature shows an increasing trend over the entire study area (Fig. 7). The significant increase in annual mean minimum temperature (Tmin) observed in the extremely arid zone is 0.017° C/a and 0.019° C/a in the arid zone, whereas it significantly decreases in the semi-arid zone by 0.007° C/a at 95% confidence level. The decrease in annual Tmin will lessen the moisture stress, crop water demands, and water loss from the ground surface and helpful for sustainable agriculture in the semi-arid zone.

Fig. 7 Time series trend of annual mean minimum temperature in different climatic zones of Pakistan during 1951-2016

The significant increase in temperature is an indication of global warming, which influences the climate in Pakistan, and will increase the ET and precipitation in the study area. The rise in temperature will enhance the convective activity that will increase the frequency and intensity of extreme rainfall events over these regions. The recent study shows that there is a strong relationship between global warming and increase in rainfall and arid areas will receive more rainfall around the world (Donat et al., 2017).

4 Conclusions

The study is conducted in the light of IPCC report (2014) for the 21st century in view of that global warming offers overview of future water demands in different climatic zones. This study describes and quantifies the impacts of meteorological parameters on seasonal and annual ET based on the climate changes in the extremely arid, arid, and semi-arid zones of Pakistan. The results show that temperature and vapor pressure are the most important meteorological factors affecting ET. The maximum changes in ET under the influences of temperature and vapor pressure are observed in the entire study area in monsoon season and winter, respectively. The maximum change in ET caused by solar radiation is observed in autumn in the extremely arid and arid zones and in monsoon season in the semi-arid zone whereas the maximum change in ET caused by wind speed is observed in summer in the extremely arid and arid zones and in winter in the semi-arid zone.

A significant increase in temperature, solar radiation, wind speed, and decrease in vapor pressure can enhance ET and put pressures on agricultural water demands, soil moisture and water resources of the study area. The increase in ET may enhance the risk of drought which ultimately would convert the agriculture land to barren land. The mean annual precipitation and annual mean maximum temperature significantly increase (at 95% confidence level) by 0.014° C/a in the extremely arid zone and 0.828 mm/a in the arid zone. The mean annual minimum temperature increases 0.017° C/a in the extremely arid zone and 0.019° C/a in the arid zone, and decreases by 0.007° C/a in the semi-arid zone at 95% confidence level. The study concludes that the increase in temperature can trigger the increase of ET over the extremely arid, arid, and semi-arid zones and will enhance the water demands for agriculture and create water scarcity problem for the region. These estimated results are important and would be helpful for future management of regional water resources and the irrigation system for sustainable development in Pakistan.

Acknowledgments

We would like to thank two anonymous reviewers for their valuable comments and suggestions.

The authors have declared that no competing interests exist.

Reference
[1] Adnan S. 2009. Agro-climatic classification of Pakistan. MSc Thesis. Islamabad, Pakistan: COMSATS Institute of Information Technology (CIIT). [2014-10-29]. https://www.researchgate.net/publication/261700354_Agroclimatic_Classification_of_Pakistan. [Cited Within:1]
[2] Adnan S, Khan A H. 2009. Effective rainfall for irrigated agriculture plains of Pakistan. Pakistan Journal of Meteorology, 6(11): 61-72. [Cited Within:1]
[3] Adnan S, Ullah K, Gao S. 2016. Investigations into precipitation and drought climatology in South Central Asia with special focus on Pakistan over the period 1951-2010 Journal of Climate, 29(16): 6019-6035. [Cited Within:2]
[4] Adnan S, Ullah K, Gao S, et al. 2017. Shifting of agro-climatic zones, their drought vulnerability, and precipitation and temperature trends in Pakistan. International Journal of Climatology, 37(Suppl. ): 529-543. [Cited Within:1]
[5] Afzaal M, Haroon M A, Zaman Q U. 2009. Interdecadal oscillations and the warming trend in the area-weighted annual mean temperature of Pakistan. Pakistan Journal of Meteorology, 6(11): 13-19. [Cited Within:1]
[6] Ali M H, Adham A K M, Rahman M M, et al. 2009. Sensitivity of Penman-Monteith estimates of reference evapotranspiration to errors in input climatic data. Journal of Agrometeorology, 11(1): 1-8. [Cited Within:1]
[7] Allen R G, Pereira L S, Raes D, et al. 1998. Crop evapotranspiration-guidelines for computing crop water requirements. In: FAO Irrigation and Drainage Paper 56. Rome, Italy: FAO, 300(9): D05109. [Cited Within:4]
[8] Bultot F, Dupriez G L, Gellens D. 1988. Estimated annual regime of energy-balance components, evapotranspiration and soil moisture for a drainage basin in the case of a CO2 doubling. Climatic Change, 12(1): 39-56. [Cited Within:1]
[9] Chaudhary Q Z, Rasul G. 2003. Agro-climatic classification of Pakistan. Science Vision, 9(1-2): 59-66. [Cited Within:1]
[10] Claussen M. 1997. Modeling bio-geophysical feedback in the African and Indian monsoon region. Climate Dynamics, 13(4): 247-257. [Cited Within:1]
[11] Donat M G, Lowry A L, Alexand er L V, et al. 2017. Addendum: More extreme precipitation in the world’s dry and wet regions. Nature Climate Change, 7(2): 154-158. [Cited Within:1]
[12] Doorenbos J, Pruitt W O. 1977. Crop water requirements. In: FAO Irrigation and Drainage Paper 24. Rome, Italy: FAO, 144. [Cited Within:1]
[13] Eagleson P S. 1978. Climate, soil, and vegetation: 1. Introduction to water balance dynamics. Water Resources Research, 14(5): 705-712. [Cited Within:1]
[14] Feddema J J. 1999. Future African water resources: interactions between soil degradation and global warming. Climatic Change, 42(3): 561-596. [Cited Within:1]
[15] Gad H E, El-Gayar S M. 2010. Effect of solar radiation on the crops evapotranspiration in Egypt. In: The 14th International Water Technology Conference. Cairo, Egypt: IWTC, 769-783. [Cited Within:1]
[16] Gong L B, Xu C Y, Chen D L, et al. 2006. Sensitivity of the Penman-Monteith reference evapotranspiration to key climatic variables in the Changjiang (Yangtze River) basin. Journal of Hydrology, 329(3-4): 620-629. [Cited Within:1]
[17] Goyal R K. 2004. Sensitivity of evapotranspiration to global warming: a case study of arid zone of Rajasthan (India). Agricultural Water Management, 69(1): 1-11. [Cited Within:1]
[18] Haider S, Adnan S. 2014. Classification and assessment of aridity over Pakistan provinces (1960-2009). International Journal of Environment, 3(4): 24-35. [Cited Within:2]
[19] Hansen J, Sato M, Ruedy R. 2012. Perception of climate change. Proceedings of the National Academy of Sciences of the United States of America, 109(37): E2415-E2423. [Cited Within:1]
[20] Held I M, Soden B J. 2006. Robust responses of the hydrological cycle to global warming. Journal of Climate, 19(21): 5686-5699. [Cited Within:1]
[21] Iglesias A, Lopez Corcoles H, Canadas W, et al. 1994. Current and future strategies for water use optimization in corn hybrids of different growth durations under climate change scenarios. In: Proceedings of the 3rd Congress of the European Society for Agronomy. Abano-Padova: Padova University, 374-375. [Cited Within:1]
[22] Intergovernmental Panel on Climate Change (IPCC). 2007. Summary for policymakers. In: Solomon S, Qin D, Manning M, et al. Climate Change 2007: the Physical Science Basis-Contribution of Working Group I to the Fourth Assessment Report of the Intergovernmental Panel on Climate Change. Cambridge: Cambridge University Press, 996. [Cited Within:2]
[23] Intergovernmental Panel on Climate Change (IPCC). 2014. 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. Geneva, Switzerland : IPCC, 151. [Cited Within:3]
[24] Isikwue C B, Audu O M, Isikwue O M. 2014. Evaluation of evapotranspiration using FAO penman-monteith method in Kano Nigeria. International Journal of Science and Technology, 3(11): 698-703. [Cited Within:1]
[25] Kendall M G. 1975. Rank Correlation Measures. London: Charles Griffin, 202. [Cited Within:1]
[26] Khan J A. 1993. The Climate of Pakistan. Karachi: Rehber Publishers, 79. [Cited Within:1]
[27] Khan S U, Hassan M, Khan F K, et al. 2010. Climate Classification of Pakistan. Ohrid: Balwois, 1-47. [Cited Within:1]
[28] Kruss P O, Khan K A Y, Malik F M Q, et al. 1992. Cooling over monsoonal Pakistan. In: Proceedings of the 5th International Meeting on Statistical Climatology. Toronto: Environment Canada, 27. [Cited Within:1]
[29] Le Houerou H N. 1993. Climate change and drought desertifisation. Revue Sécheresse, 4(2): 95-111. [Cited Within:1]
[30] Liu Q, Yang Z F, Cui B S, et al. 2010. The temporal trends of reference evapotranspiration and its sensitivity to key meteorological variables in the Yellow River Basin, China. Hydrological Processes, 24(15): 2171-2181. [Cited Within:1]
[31] Mamassis N, Panagoulia D, Novcovic A. 2014. Sensitivity analysis of Penman evaporation method. Global Network for Environmental Science and Technology, 16(4): 628-639. [Cited Within:1]
[32] Mann H B. 1945. Non-parametric tests against trend. Econometrica, 13: 245-259. [Cited Within:1]
[33] Martin P, Rosenberg N J, McKenney M S. 1989. Sensitivity of evapotranspiration in a wheat field, a forest, and a grassland to changes in climate and direct effects of carbon dioxide. Climatic Change, 14(2): 117-151. [Cited Within:1]
[34] Mavi H S. 1986. Introduction to Agrometeorology. New Delhi: Oxford and IBH Publishing Co. , 82-97. [Cited Within:1]
[35] Michael A M. 1986. Irrigation Theory and Practice. Sahibabad: Vani Educational Books, 512-537. [Cited Within:1]
[36] Monteith J L. 1965. Evaporation and environment. In: Fogg G E. The State and Movement of Water in Living Organisms. New York: Academic Press, 205-234. [Cited Within:1]
[37] Sabziparvar A A, Tabari H. 2010. Regional estimation of reference evapotranspiration in arid and semiarid regions. Journal of Irrigation and Drainage Engineering, 136(10): 724-731. [Cited Within:1]
[38] Sen P K. 1968. Estimates of the regression coefficient based on Kendall’s tau. Journal of the American Statistical Association, 63(324): 1379-1389. [Cited Within:1]
[39] Shamshad K M. 1988. The Meteorology of Pakistan: Climate and Weather of Pakistan. Karachi: Royal Book Company Publishers, 313. [Cited Within:1]
[40] Singh N, Sontakke N A. 1996. Climate variability over Pakistan and its relationship to variations over the Indian region. In: Abrol Y P, Gadgil S, Pant G B. Climate Variability and Agriculture. New Delhi: Narosa Publishing House, 67-95. [Cited Within:1]
[41] Smakhtin V U, Schipper E L F. 2008. Droughts: The impact of semantics and perceptions. Water Policy, 10(2): 131-143. [Cited Within:1]
[42] Subramanya K. 1984. Engineering Hydrology. New Delhi: Tata McGraw-Hill Publishing Company, 48-82. [Cited Within:1]
[43] Subramanya K. 2008. Engineering Hydrology (3rd ed. ). New Delhi: Tata McGraw-Hill Publishing Company, 60-80. [Cited Within:1]
[44] Viessman W Jr, Knapp W J, Lewis L G, et al. 1977. Introduction to Hydrology (2nd ed. ). New York: IEP, Dun-Donnelley Publishing Corporation, 704. [Cited Within:1]
[45] Wilson E M. 1974. Engineering Hydrology. Hampshir: MacMillan Education Ltd, 34-52. [Cited Within:1]
[46] Xu C Y, Gong L B, Jiang T, et al. 2006. Decreasing reference evapotranspiration in a warming climate-a case of Changjiang (Yangtze) River catchment during 1970-2000. Advances in Atmospheric Sciences, 23(4): 513-520. [Cited Within:1]
[47] Zahid M, Rasul G. 2011. Thermal classification of Pakistan. Atmospheric and Climate Sciences, 1(4): 206-213. [Cited Within:1]
[48] Zhao L W, Zhao W Z. 2014. Evapotranspiration of an oasis-desert transition zone in the middle stream of Heihe River, Northwest China. Journal of Arid Land , 6(5): 529-539. [Cited Within:1]