Fifth IAERE Annual Conference 16 17 February 2017, Rome Climate change impact on Ethiopian small holders production efficiency Solomon Asfaw, FAO, Agricultural Development Economics Division (ESA) Sabrina Auci, University of Palermo, European Studies and International Relations Department (DEMS) Manuela Coromaldi, University of Rome Niccolò Cusano
Building blocks Motivation Country background and climate variability Data description Empirical strategy Estimation results Conclusions 2
Motivation (1) Many of the anticipated adverse effect of climate change (CC) such as increases in global temperature, sea level rise, enhanced monsoon precipitation and increase in drought intensity aggravate the development of agriculture-based economies such as Sub-Saharian countries the aim of this study is to provide an analysis of the impact of weather risk on small farmers technical efficiency climate change effects analysis of Ethiopia is interesting because of the presence of different microclimates 3
Motivation (2) Besides the important policy implications that can be derived from the investigation of these issues, we focus on weather-related risk for the following two reasons: the growing availability of high quality geo-referenced data on weather CC is one of the most important exogenous factor affecting income and consumption of small rural households 4
Country background Ethiopia, like many other sub Saharan countries, weightily relies on agricultural sector and consequently the economic impact of CC is crucial for small-scale farmers food security and welfare the agriculture sector contributes about 46% of total GDP agricultural production is completely dependent on rainfall 85% of population dependent on rain-fed agriculture for their livelihood (ACCRA, 2010) 5
Climate variability Since the 1970s the severity and the frequency of droughts have increased and the areas affected by drought and desertification are expanding (World Bank 2013, 2014) Major floods, which caused significant damage, with numerous livestock deaths and damage to planted crops and stored food, occurred in different parts of Ethiopia in 1988, 1993, 1994, 1995, 1996, 2006 and 2012 (NAPA, 2007; UN OCHA-Ethiopia, 2013) Regional projections of climate models not only predict a substantial rise in frequency of both extreme flooding and droughts, but also suggest an increase in mean temperatures over the 21 st century due to global warming (EACC, 2010) 6
Data description Two main sources of data: socio-economic data from World Bank Ethiopia Living Standards Measurement Study - Integrated Surveys on Agriculture (LSMS-ISA) 2011/2012 and 2013/2014; Historical analysis data on rainfall and temperature from the National Oceanic and Atmospheric Administration (NOAA) and the European Centre for Medium Range Weather Forecasts (ECMWF). 7
Table 1. Descriptive statistics Model Var Description mean p50 sd min max N Production function inputs Inefficiency model factors Y L Chemfert Seed H Rainfall_sd Tmax_sd Yields=Crop harvested (maize, sorghum, barley, teff, wheat) /Operated land (Kg/Ha) Days of work employed / Operated land Use of Chemical Fertilizer in kilos / Operated land (Kg/Ha) Seeds used in kilos / Operated land (Kg/Ha) Average years of education in HH Standard deviation of rainfall in the last 5 year Standard deviation of maximum temperature in the last 5 year 26,276.71 12,353.93 72,458.73 464.79 3,248,693 4366 751.17 405.47 2,033.61 14.29 93,232.91 4366 153.48 88.80 406.09 3.67 18,622.16 4366 167.05 85.85 495.61 2.61 18,699.85 4366 2.54 2.00 2.84 0.00 20.00 4366 31.20 28.51 8.10 13.07 55.35 4406 2.18 2.06 0.81 0.86 4.16 4406 8
Table 2. Descriptive statistics by Ethiopian regions (mean value) Region Y L Chemfert Seed H Rainfall_sd tmax_sd Tigray 34,124.57 916.19 229.52 213.34 2.86 37.25 2.93 Afar 38,378.25 1262.14 220.24 238.87 2.19 28.07 1.59 Amhara 19,212.16 498.12 142.53 120.69 2.52 40.20 2.56 Oromiya 27,476.06 746.98 153.18 146.91 2.58 29.72 2.20 Somalie 13,238.40 417.91 63.05 79.25 1.79 23.37 1.31 Benishangul Gumuz 37,381.07 1251.26 240.45 245.94 2.45 32.51 3.18 SNNP 26,921.05 735.13 170.63 146.36 2.62 24.51 1.82 Gambella 50,199.76 1471.72 282.42 301.41 3.96 24.34 2.43 Harari 24,223.16 1080.31 174.91 153.43 1.72 26.57 1.22 Dire Dawa 24,194.53 891.09 133.32 144.66 1.91 28.89 1.05 9
Empirical strategy (1) The applied SF models allow : to control for time-invariant inefficiency (Pitt and Lee, 1981) to disentangle time-invariant heterogeneity from time-varying inefficiency ( true fixed- and random- effects; Greene, 2005) to control for territorial heterogeneity clustering our estimation by regions to control for heteroskedasticity: specifying explanatory variables for the inefficiency variance function (Greene, 2007) 10
Empirical strategy (2) The use of exogenous variables such as climatic variables in the stochastic frontier approach (SFA) allows to capture the effects of the environment in which a farmer produces his agricultural output (Kumbhakar and Knox Lovell, 2000) We include the weather variables in the inefficiency model to estimate technical efficiency The incorporation of exogenous variables within the estimation of technical efficiency means the incorporation of features beyond the control of the farmer and the separation of the production frontier model from the inefficiency equation 11
Empirical strategy (3) Pitt and Lee (1981) assume that the inefficiency term is farmer-specific so that technical efficiency is constant over time (time-invariant model). The production function can be expressed as: (1) Y it = + x it +(v it u i ) i=1,..., N and t=1,, T where v it ~ N(0,σ v2 ) and u i is distributed as a half-normal random variable N + (0,σ u2 ). 12
Empirical strategy (4) The limitation of the time-invariant model can be overcome following Greene (2005): a time-varying SFA with a household-specific fixed or stochastic term i is applied. The simultaneous estimation of production input coefficients and inefficiency factor parameters is based on a two-stage maximum-likelihood procedure (Kumbhakar et al., 1991). The production function can be expressed as: (2) Y it =x it +(v it -u it )+ i i=1,..., N and t=1,, T where v it ~ N(0,σ v2 ) and u it is distributed as a half-normal random variable N + (0,σ u2 ). 13
Empirical strategy (5) To analyze the inefficiency determinants related to observed heterogeneity of Ethiopian households which directly affects the inefficiency u it, we include efficiency determinants in the variance of the inefficiency error as in Greene (2007) and Hadri et al. (2003): (3) 2 uit exp z it where z it is a px1 vector of variables that may have an indirect effect on the production function of an Ethiopian household; γ is a 1xp vector of parameters to be estimated 14
Model specification (1) The standard the Cobb-Douglas (C-D) production function is applied: lny it = 0 lnseed 4 + lnl it 1 + v it it lnh u it 2 i it lnchem_fer 3 t it where Y it is the yield per hectare in kilos L it are days of work employed H it is average years of education in HH Chem_fert it is chemical fertilizer in kilos Seed it is the amount of seeds used in kilos 15
Model specification (2) The inefficiency model can be specified as follows: 2 it exp 0 + 1Rainfall_sd it + = Tmax_sd 2 itit it Rainfall_sd: the standard deviation of September-March rainfall average for the last five years by enumeration area Tmax_sd: the standard deviation of September-March maximum temperature average for the last five years by enumeration area Following Dell, Jones, and Olken (2012), we take the lags of weather variability 16
Technical efficiency The technical efficiency of the i-th Ethiopian household using the Battese and Coelli (1988) estimation is given by : TE it e ( u it ) ( z ti ti e ) Technical in/efficiency values will oscillate between 0 and 1 If TE i <1 then the observable output is less than the maximum feasible output, meaning that the statistical unit is not efficient 17
Table 3. Stochastic frontier estimation results PITT and LEE (1981) TFE without inefficiency model (GREENE, 2005) TRE without inefficiency model (GREENE, 2005) TFE with inefficiency model (GREENE, 2007) TRE with inefficiency model (GREENE, 2007) Production Function Frontier Model lnl 0.567*** 0.560*** 0.369*** 0.432*** 0.351*** (6.697) (3.39) (4.44) (2.63) (4.362) lnseed 0.211** 0.0757 0.0619 0.208 0.0603 (2.038) (0.575) (0.675) (1.597) (0.677) lnqchem_fert 0.184*** 0.229*** 0.140*** 0.234*** 0.148*** (3.517) (3.098) (2.842) (4.729) (3.01) H 0.0235*** 0.0923*** 0.0258*** 0.0899*** 0.0266*** (2.901) (5.605) (3.04) (3.816) (3.117) Constant 4.212*** 5.850*** 5.996*** (31.2) (20.51) (20.09) Note: t statistics in parentheses. * significant at 10%; ** significant at 5%; *** significant at 1% 18
Table 3. Stochastic frontier estimation results (continued) Inefficiency model Rainfall_sd -0.0407** -0.125*** (-2.127) (-4.049) Tmax_sd 0.648** 1.554*** (2.268) (3.442) Constant -0.966-4.206*** (-1.38) (-2.982) Obs 4,366 4,366 4,366 4,366 4,366 N. of households 2,183 2,183 2,183 2,183 2,183 Wald chi2 11056*** 3092*** 783.3*** 1422*** 570.4*** Prob > chi2 0.00 0.00 0.00 0.00 0.00 N_clust 10 10 10 10 10 Log Likelihood -3961-1631 -3834-1360 -3817 sigma_u 0.000232 0.694 0.0016 sigma_v 0.599 0.0277 0.582 0.0111 0.576 avg_sigmau 0.677 0.106 Note: t statistics in parentheses. * significant at 10%; ** significant at 5%; *** significant at 1% 19
Table 4. ANOVA Analysis of technical inefficiency by Ethiopian regions True random effects model without inefficiency model True random effects model with inefficiency model Region Code Mean efficiency Region Code Mean efficiency Harari 0.9981986 Gambella 0.8587259 Oromiya 0.9981990 Benishangul Gumuz 0.8658033 SNNP 0.9981990 Tigray 0.909941 Dire Dawa 0.9981990 SNNP 0.9150778 Tigray 0.9981991 Oromiya 0.916757 Amhara 0.9981992 Somalie 0.9371281 Benishangul Gumuz 0.9981992 Amhara 0.9384624 Somalie 0.9981993 Afar 0.9412216 Afar 0.9981995 Harari 0.9535999 Gambella 0.9981995 Dire Dawa 0.9647118 ANOVA test F 1.40 261.16 Prob > F 0.1833 0.000 20
Conclusion (1) This paper contributes to the CC literature: by using a novel data set that combines information coming from two large-scale household surveys with geo-referenced historical rainfall and temperature data by estimating the technical efficiency of Ethiopian households in the period 2011-12 and 2013-14 by analyzing the CC determinants on farmers efficiency taking into account the heteroskedasticity by considering the heterogeneity of climate areas (desert, tropical forest and plateau) by ranking the Ethiopian regions from the less efficient to the most efficient 21
Conclusion (2) The results of the production functions are generally conform to our expectations and to literature In the inefficiency model, we find that short-term weather variables computed on precipitations and maximum temperatures in the cropping season affect negatively and positively Ethiopian farmers inefficiency respectively Regions in the west side of the country (such as Tigray, Benishangul Gumuz, Gambella) are less efficient than regions in the north-west side (such as Afar, Dire Dawa, Harari) 22
Further developments A series of intra- and inter-seasonal indicators to measure both the level and the variability of rainfall and maximum temperature should be calculated: o intra-seasonal indicators measure weather fluctuations such as Walsh and Lawler (1981) Seasonality Index for Precipitations o inter-seasonal indicators measure climate shocks such as Coefficient of Variation computed for the past 10, 5 and 3 years o shortfall indicators measure how much less rain fell in the year t during the cropping season compared to average of the last 25, 10, 5, 3 years 23
Thank you! 24
Table 3. Stochastic frontier estimation results PITT and LEE (1981) True Random Effects (TRE) without inefficiency model (GREENE, 2005) True Random Effects (TRE) with inefficiency model (GREENE, 2007) Production Function Frontier Model lnl 0.567*** 0.567*** 0.351*** (6.697) (6.696) (4.362) lnseed 0.211** 0.211** 0.0603 (2.038) (2.037) (0.677) lnqchem_fert 0.184*** 0.184*** 0.148*** (3.517) (3.518) (3.01) H 0.0235*** 0.0235*** 0.0266*** (2.901) 2.901 (3.117) Constant 4.212*** 4.214*** 5.996*** (31.2) 15.70 (20.09) Note: t statistics in parentheses. * significant at 10%; ** significant at 5%; *** significant at 1% 25
Table 3. Stochastic frontier estimation results (continued) Inefficiency model Rainfall_sd 0.125*** ( 4.049) Tmax_sd 1.554*** (3.442) Constant 4.206*** ( 2.982) Obs 4,366 4,366 4,366 N. of households 2,183 2,183 2,183 Wald chi2 11056*** 11050 *** 570.4*** Prob > chi2 0.00 0.00 0.00 N_clust 10 10 10 Log Likelihood 3961 3961 3817 sigma_u 0.000232 0.00226 sigma_v 0.599 0.599 0.576 avg_sigmau 0.106 Note: t statistics in parentheses. * significant at 10%; ** significant at 5%; *** significant at 1% 26