Climate change is one of the greatest challenges in the 21st century and the agriculture sector is very vulnerable to this phenomenon. Since wheat is the most important cereal crop in Iran, we aim to analyze the potential impact of climatic variables (temperature and precipitation) on rainfed wheat productivity in Hamedan Province, Iran. For this purpose, generalized additive models have been used to model yields of rainfed wheat based on climatic variables during 2004–2012. Then, based on sensitivity of rainfed wheat to temperature and precipitation in this period, we predict the potential effects of climate change on rainfed wheat yield under the IPCC SRES A1FI and B1 climate change scenarios. Results suggest that yields of rainfed wheat would decrease in all Hamedan's counties primarily because of decreasing October to June precipitation and higher temperature. As a result, it is predicted that the yield of rainfed wheat in Hamedan under the A1F1 and B1 scenarios will fall by 41.3% and 20.6%, respectively, in the 2080s. In other words, according to the A1F1 scenario, in the 2080s, Hamedan Province's rainfed wheat production will decline from 1090 kg/ha to 639 kg/ha and under the B1 scenario to 865 kg/ha.
INTRODUCTION
Climate change is one of the greatest challenges of the 21st century and agriculture may be particularly vulnerable to climate change (IPCC 2007; Brown & Funk 2008; Khan et al. 2009; Audsley et al. 2010; Hille et al. 2012) because of its dependence on natural weather patterns and climate cycles for productivity. The impacts of climate change on agricultural production are both direct and indirect. Important direct effects will be through changes in temperatures, precipitation and radiation, more frequent extreme weather (like droughts, floods, and wind storms), atmospheric carbon dioxide concentration (Bannayan et al. 2007; Smith et al. 2007; Lobell et al. 2008), length of growing season (Saarikko & Carter 1996) and modifying evaporation, runoff, and finally soil moisture (Rosenzweig & Hillel 1998). Indirect effects will include potentially deterministic changes in diseases, pests, weeds, and sea level rise (Keane et al. 2009), the effects of which have not yet been quantified in most studies. All these factors and variables can change yield and agricultural productivity (Battisti & Naylor 2009). Although there will be gains in some crops in some regions of the world, the overall impacts of climate change on agriculture are expected to have a negative effect, threatening global food security. Predictions show that global agricultural productivity would fall by 15.9% in the 2080s if global warming continues unabated (Cline 2007) and developing countries lying in the tropical and sub-tropical regions would face catastrophic results (Morton 2007; Morison & Morecroft 2008; Mann et al. 2009).
In Iran, agriculture is the second largest sector of the economy, after the service industry, accounting for about 26 of GDP and 26% of non-oil exports. It also provides nearly 23% of the total country's employment, and 80% of food production (Dabiri et al. 2013). Nonetheless, Iran has arid or semiarid climates with long dry summers and low winter rainfall, average annual precipitation around 250 mm and high potential evapotranspiration. Such climatic characteristics associated with phenomena such as desertification, drought, water table reduction and flooding increment and vulnerability of land resources (Momeni 2003) become very vulnerable to climate change. Findings of recent studies have confirmed and documented the occurrence of climate change (IRIMO 2006a, 2006b; Tabari & Talaee 2011) and its damaging potential future effects on the Iranian agricultural sector (see Bolle 2003; Kouchaki & Nasiri 2008; Dastorani 2012; Dehghanpour et al. 2014).
Wheat is the most important cereal crop in Iran, accounting for 35% of the food grain production of the country (12 million tonnes in 2004). Rainfed wheat accounts for about 60–65% of the land area under wheat production in Iran and contributes 30–35% of wheat production in the country. The share of wheat in the daily energy supply is 47% (Ministry of Agriculture Jahad of Iran 2000) and the demand for wheat is predicted to be over 20 million tonnes in 2020 (Sharifi 2001). Rainfed agriculture is much more vulnerable to climate change (Trethowan & Pfeiffer 1999; Ludwig & Asseng 2006). For these reasons, we aimed to analyze the impacts of climate change on rainfed wheat production in the Hamedan Province as one of the main wheat producer provinces in the country.
MATERIALS AND METHODS
Study area
Location of the study area within Iran (left) and its existing land use/cover (right).
Location of the study area within Iran (left) and its existing land use/cover (right).
Average temperature (°C) and total precipitation (mm) for the period 2004–2012 in Hamedan Province.
Average temperature (°C) and total precipitation (mm) for the period 2004–2012 in Hamedan Province.
Data set and climate model
Although radiation, soil moisture and carbon dioxide (CO2) concentration are all important variables to determine wheat productivity and assessing potential impacts of climate change, temperature and precipitation effects are more significant (Warrick 1988; Chiotti & Johnston 1995; Wang et al. 2009). For this reason, and also because of the lack of accurate information about other variables, the temperature and rainfall data has been used to analyze the effects of climate change on rainfed wheat production in Hamedan Province, for current and future years. For this purpose, climatic data series (minimum and maximum temperature and precipitation) of synoptic stations of Hamedan Province were obtained for the period 2004–2012 from the Iran Meteorological Organization. Wheat production data series used in this study cover the same period from 2004 to 2012 and include data on cultivated areas, production and yield (tonnes/ha) under rainfed conditions. The choice of these periods was justified by the availability of wheat productions statistics for the areas at those periods.
The SRES A1FI (highest future emission trajectory) and B1 (lowest future emission trajectory) climate scenarios for three time slices, namely 2020s (2010–2039), 2050s (2040– 2069) and 2080s (2070–2099), were chosen to evaluate climate change impacts (Table 1).
Projected changes in surface air temperature and precipitation for west Asia (IPCC 2007)
. | 2010–2039 . | 2040–2069 . | 2070–2099 . | |||||||||
---|---|---|---|---|---|---|---|---|---|---|---|---|
. | T (C°) . | P (%) . | T (C°) . | P (%) . | T (C°) . | P (%) . | ||||||
Season . | A1F1 . | B1 . | A1F1 . | B1 . | A1F1 . | B1 . | A1F1 . | B1 . | A1F1 . | B1 . | A1F1 . | B1 . |
DJF | 1.26 | 1.06 | −3 | −4 | 3.1 | 2 | −3 | −5 | 5.1 | 2.8 | −11 | −4 |
MAM | 1.29 | 1.24 | −2 | −8 | 3.2 | 2.2 | −8 | −9 | 5.6 | 3 | −25 | −11 |
JJA | 1.55 | 1.53 | 13 | 5 | 3.7 | 2.5 | 13 | 20 | 6.1 | 2.7 | 32 | 13 |
SON | 1.48 | 1.35 | 18 | 13 | 3.6 | 2.2 | 27 | 29 | 5.7 | 3.2 | 52 | 25 |
. | 2010–2039 . | 2040–2069 . | 2070–2099 . | |||||||||
---|---|---|---|---|---|---|---|---|---|---|---|---|
. | T (C°) . | P (%) . | T (C°) . | P (%) . | T (C°) . | P (%) . | ||||||
Season . | A1F1 . | B1 . | A1F1 . | B1 . | A1F1 . | B1 . | A1F1 . | B1 . | A1F1 . | B1 . | A1F1 . | B1 . |
DJF | 1.26 | 1.06 | −3 | −4 | 3.1 | 2 | −3 | −5 | 5.1 | 2.8 | −11 | −4 |
MAM | 1.29 | 1.24 | −2 | −8 | 3.2 | 2.2 | −8 | −9 | 5.6 | 3 | −25 | −11 |
JJA | 1.55 | 1.53 | 13 | 5 | 3.7 | 2.5 | 13 | 20 | 6.1 | 2.7 | 32 | 13 |
SON | 1.48 | 1.35 | 18 | 13 | 3.6 | 2.2 | 27 | 29 | 5.7 | 3.2 | 52 | 25 |
DJF = Dec., Jan., Feb.; MAM = Mar., Apr., May; JJA = June, Jul., Aug.; SON = Sep., Oct., Nov. T (C°) = Temperature (Celsius); P (%) = Precipitation (percentage); A1FI = highest future emission trajectory; B1 = lowest future emission trajectory.
Generalized additive models
Many researchers have evaluated the possible impact of climate change on crop yields using mainly indirect crop simulation models (see Rosenzweig & Parry 1994; Reilly et al. 2003; Izaurralde et al. 2003; Nassiri et al. 2006; Wang et al. 2013). There are relatively few direct assessments on the impact of observed climate change based on past crop yield (see Carter & Zhang 1998; Naylor et al. 2002; Peng et al. 2004; De la Casa & Ovando 2014; Nara et al. 2014) and then prediction of the future based on fitted trends. In this study we have used the generalized additive models (GAMs) to model past climatic trends and their effects on rainfed wheat yield. Then, based on the sensitivity of rainfed wheat to temperature and precipitation in the period 2004–2012, we predicted the potential effects of climate change on rainfed wheat in Hamedan Province.
Generalized additive models (Hastie & Tibshirani 1986) extend the well-known generalized linear models (GLMs) by allowing for nonlinear covariate effects through the incorporation of additive non-parametric smooth functions. In the univariate case, assume the distribution of the response variable belongs to the exponential family and that its mean
related to explanatory variables
,
, … through the use of a link function
. GAMs are very flexible and there is no need for any assumption about the parametric form for the dependence of the response variable
on variables
,
, … When the mean relationships are complex and cannot be easily modeled by specific linear or nonlinear functions, a GAM can be a natural choice. Due to their flexibility in model specification, GAMs have seen a variety of applications. GAMs have been widely used in environmental, biological, and ecological studies.



Box plot of wheat yield and October to June precipitation for the period 2004–2012 in Hamedan Province.
Box plot of wheat yield and October to June precipitation for the period 2004–2012 in Hamedan Province.
RESULTS AND DISCUSSION
Fitting the proposed model
Fitted partial smooth terms of generalized additive model showing the effects of the year, rain and temperature on yield in Hamadan Province from 2004 to 2012. The dotted curves are 95% confidence intervals of smooth functions.
Fitted partial smooth terms of generalized additive model showing the effects of the year, rain and temperature on yield in Hamadan Province from 2004 to 2012. The dotted curves are 95% confidence intervals of smooth functions.
Analyzing climate change impacts on wheat production
Impacts of climate change on rainfed wheat in Hamedan's counties under the SRES A1FI and B1 climate change scenarios, where Y.Kg/H is yield (kg/ha) and P/T is production (tonnes)
. | . | Asadabad . | Bahar . | Hamedan . | Kaboodarahang . | Malayer . | Nahavand . | Razan . | Tuyserkan . | Hamedan Province . | Total reduction (%) . | ||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|
Period . | Scenario . | Y.Kg/H . | P/T . | Y.Kg/H . | P/T . | Y.Kg/H . | P/T . | Y.Kg/H . | P/T . | Y.Kg/H . | P/T . | Y.Kg/H . | P/T . | Y.Kg/H . | P/T . | Y.Kg/H . | P/T . | Y.Kg/H . | P/T . | Y.Kg/H . | P/T . |
Current (2004–2012) | A1F1 | 1,393.9 | 131,144 | 1,051.3 | 282,020 | 829.4 | 564,039 | 981.6 | 833,972 | 948.2 | 281,774 | 1,179.3 | 95,526 | 1,123.2 | 583,799 | 1214.3 | 72,103 | 1090.2 | 2,844,375 | − | − |
2020s | 948.1 | 92,392 | 948.1 | 264,468 | 948.1 | 677,443 | 930.1 | 792,400 | 894.8 | 267,463 | 1,015.8 | 81,837 | 930.8 | 4,479,694 | 1119.3 | 67,555 | 980.4 | 2,723,250 | 10 | 4.2 | |
2050s | 837.1 | 78,493 | 837.1 | 224,914 | 837.1 | 576,235 | 783.21 | 667,995 | 744.5 | 223,215 | 872.1 | 70,051 | 783.9 | 404,422 | 981.3 | 58,977 | 834.6 | 2,304,301 | 23.4 | 19 | |
2080s | 664.5 | 60,316 | 664.5 | 173,225 | 664.5 | 443,904 | 593.7 | 507,652 | 541 | 163,335 | 670 | 53,544 | 953.9 | 307,620 | 779.7 | 46,373 | 639 | 1,755,968 | 41.3 | 38.2 | |
2020s | B1 | 1,003.4 | 94,206 | 1,003.4 | 269,560 | 1,003.4 | 690,606 | 947.61 | 807,332 | 914.6 | 273,360 | 1,038.8 | 83,694 | 948.6 | 488,680 | 1147.4 | 69,258 | 1001 | 2,776,721 | 8.1 | 2.3 |
2050s | 936.2 | 87,860 | 936.2 | 251,526 | 936.2 | 644,431 | 883 | 752,394 | 842.4 | 252,139 | 968.8 | 77,982 | 883.7 | 455,562 | 1078 | 64,953 | 933.1 | 2,586,840 | 14.4 | 9 | |
2080s | 868.9 | 81,507 | 868.1 | 233,507 | 868.1 | 598,199 | 816.15 | 695,860 | 775.8 | 232,429 | 901.5 | 72,474 | 816.8 | 421,326 | 1008.3 | 60,649 | 865.7 | 2,395,954 | 20.6 | 15.7 |
. | . | Asadabad . | Bahar . | Hamedan . | Kaboodarahang . | Malayer . | Nahavand . | Razan . | Tuyserkan . | Hamedan Province . | Total reduction (%) . | ||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|
Period . | Scenario . | Y.Kg/H . | P/T . | Y.Kg/H . | P/T . | Y.Kg/H . | P/T . | Y.Kg/H . | P/T . | Y.Kg/H . | P/T . | Y.Kg/H . | P/T . | Y.Kg/H . | P/T . | Y.Kg/H . | P/T . | Y.Kg/H . | P/T . | Y.Kg/H . | P/T . |
Current (2004–2012) | A1F1 | 1,393.9 | 131,144 | 1,051.3 | 282,020 | 829.4 | 564,039 | 981.6 | 833,972 | 948.2 | 281,774 | 1,179.3 | 95,526 | 1,123.2 | 583,799 | 1214.3 | 72,103 | 1090.2 | 2,844,375 | − | − |
2020s | 948.1 | 92,392 | 948.1 | 264,468 | 948.1 | 677,443 | 930.1 | 792,400 | 894.8 | 267,463 | 1,015.8 | 81,837 | 930.8 | 4,479,694 | 1119.3 | 67,555 | 980.4 | 2,723,250 | 10 | 4.2 | |
2050s | 837.1 | 78,493 | 837.1 | 224,914 | 837.1 | 576,235 | 783.21 | 667,995 | 744.5 | 223,215 | 872.1 | 70,051 | 783.9 | 404,422 | 981.3 | 58,977 | 834.6 | 2,304,301 | 23.4 | 19 | |
2080s | 664.5 | 60,316 | 664.5 | 173,225 | 664.5 | 443,904 | 593.7 | 507,652 | 541 | 163,335 | 670 | 53,544 | 953.9 | 307,620 | 779.7 | 46,373 | 639 | 1,755,968 | 41.3 | 38.2 | |
2020s | B1 | 1,003.4 | 94,206 | 1,003.4 | 269,560 | 1,003.4 | 690,606 | 947.61 | 807,332 | 914.6 | 273,360 | 1,038.8 | 83,694 | 948.6 | 488,680 | 1147.4 | 69,258 | 1001 | 2,776,721 | 8.1 | 2.3 |
2050s | 936.2 | 87,860 | 936.2 | 251,526 | 936.2 | 644,431 | 883 | 752,394 | 842.4 | 252,139 | 968.8 | 77,982 | 883.7 | 455,562 | 1078 | 64,953 | 933.1 | 2,586,840 | 14.4 | 9 | |
2080s | 868.9 | 81,507 | 868.1 | 233,507 | 868.1 | 598,199 | 816.15 | 695,860 | 775.8 | 232,429 | 901.5 | 72,474 | 816.8 | 421,326 | 1008.3 | 60,649 | 865.7 | 2,395,954 | 20.6 | 15.7 |
Rainfed wheat yield reduction in Hamedan's counties based on A1FI and B1 scenarios in the 2080s.
Rainfed wheat yield reduction in Hamedan's counties based on A1FI and B1 scenarios in the 2080s.
Prediction for yield in Hamadan based on the data in Table 1 in terms of mean temperature in JFM.
Prediction for yield in Hamadan based on the data in Table 1 in terms of mean temperature in JFM.
Prediction for yield in Hamadan based on the data in Table 1 in terms of mean temperature in AMJ.
Prediction for yield in Hamadan based on the data in Table 1 in terms of mean temperature in AMJ.
Prediction for yield in Hamadan based on the data in Table 1 in terms of mean temperature in JAS.
Prediction for yield in Hamadan based on the data in Table 1 in terms of mean temperature in JAS.
Prediction for yield in Hamadan based on the data in Table 1 in terms of mean temperature in OND.
Prediction for yield in Hamadan based on the data in Table 1 in terms of mean temperature in OND.
CONCLUSIONS
Climate change is now recognized as one of the most challenging and complex problems facing the globe. Literature review shows that climate change is global; likewise its impacts; but the most adverse effects will be felt mainly by developing countries. The agriculture sector, especially in terms of rainfed crop yields, may be particularly vulnerable to this phenomenon because of the great dependence on natural weather patterns and climate cycles for productivity.
Since wheat is the most important cereal crop in Iran in terms of nutrition and food security, the present study was designed to determine the impact of variation in temperature and rainfall on rainfed wheat in Hamedan Province. In order to quantify the effect of annual precipitation and temperature on rainfed wheat, precipitation and temperature are usually assumed to be linearly related to the indicators of wheat production in the literature, but the relation to the calendar time and weather variables is not assumed to be parameterized. For this purpose, the generalized additive model has been used to analyze past climatic trends and their effects on rainfed wheat yield.
Then, based on sensitivity of rainfed wheat to temperature and precipitation in the period 2004–2012, we predict the potential effects of climate change on rainfed wheat yield in Hamedan Province. To take into account the seasonal effect, we also included the year as a covariate in the proposed model. We also studied the effect of latitude, longitude and elevation of the eight counties of Hamadan Province but since they are not statistically significant we have removed them from the final model. For the proposed model fitting, we used the mgcv package in R statistical software which provides functions for generalized additive modeling and generalized additive mixed modeling.
The results derived from the analysis by R statistical software indicate that the proposed model is a reasonable model with normally distributed residuals. Approximately 67.4% of the variation in yield is well described by the fitted model. Finally, based on the SRES A1FI (highest future emission trajectory) and B1 (lowest future emission trajectory) climate change scenarios we analyzed the potential impacts of climate change on rainfed wheat in Hamedan Province. Results suggest that yields of rainfed wheat would decrease in all Hamedan's counties primarily because of decreasing October to June precipitation (rainfed wheat-growing seasons) and higher temperature under the climate change scenarios. Future research needs to take into account the all-important variables such as radiation, soil moisture and carbon dioxide (CO2) concentration to determine wheat productivity and assess potential impacts of climate change in Hamedan Province.
ACKNOWLEDGEMENTS
The authors would like to thank the two referees and the editor whose comments have been very useful in improving the manuscript.