Does adaptation to climate change and variability provide household food security? Evidence from Muger sub-basin of the upper Blue-Nile, Ethiopia

The purpose of this study is to evaluate the associated impact of adoption of adaptation options to climate change and variability on household food security in the Muger sub-basin of the upper Blue-Nile of Ethiopia using a comprehensive data of 442 sampled households from four representative districts’ in the sub-basin. The study used a propensity score matching approach to evaluate the impact of adaptation options on household food security. Results show that the decision to adopt adaptation options is found to be positively influenced by male household heads, family size, access to extension service, the size of landholding, and frequency of drought and floods over the past many years. The results further reveal that farmers adopting any of the adaptation options had higher food calorie intake per day per adult equivalent than those who did not. A policy that promotes the adoption of soil and water conservation measures, small-scale irrigation, agronomic practices, and livelihood diversification strategies should be central to food security strategy in the study area.


Introduction
The impact of climate change on food security in Sub-Saharan Africa have largely been explored by a plethora of authors using either agronomic model or Ricardian analysis (Thornton and Herrero 2014;Beddington et al. 2012; Thornton et al. 2012;Conway 2011;Deressa and Hassan 2010). Empirical literature shows that Ethiopia is the most vulnerable country owing to its least adaptive capacity and low diversified economies (Stige et al. 2006). This becomes even more complicated where Ethiopia's agricultural systems have largely relied on rain-fed that has been closely associated with climate (World Bank 2006). It has been noted that extreme climate events such as drought and floods reduced one third of Ethiopia's economic growth (World Bank 2006).
A large body of literature has recognized adaptation as one of the policy options in response to climate change impact (Smit et al. 1999;Smith and Lenhart 1996;UNFCCC 1992). It is seen as critical to a great extent to reduce the ultimate effect of climate change on agriculture so as to improve livelihoods and food security of rural households in the continent (van de Giesen et al. 2010;Vermeulen et al. 2012). "Adaptation to climate change refers to adjustment in natural or human systems in response to actual or expected climatic stimuli or their effects, which moderates harm or exploits beneficial opportunities" (FAO 2011;IPCC 2011).
Given serious problem posed by climate change, many potential adaptation options have been suggested for developing countries. For instance, soil and water conservation (SWC) practices have been suggested in response to soil erosion problem posed by climate Simane 2017a, 2017b;McCarthy et al. 2011;Amsalu and De Graaff 2007). Based on the data from a comprehensive survey of agricultural households across three agro-ecologies in Muger River sub-basin of the Blue-Nile Basin, Amare and Simane (2017b) identified soil and water conservation practices are most widely used adaptation option in response to climate change. The use of these adaptation options was found to reduce soil erosion associated with short but heavy rains. Farmers are adapting SWC practices to retain soil-water content and maintain humidity during dry spells through an improved soil structure (McCarthy et al. 2011). Similarly, the use of different agronomic practices is considered as the potential adaptation option to the adverse effects of climate change on agriculture. The analysis by Amare and Simane (2017b) in the Muger River sub-basin showed that using agronomic practices such as droughttolerant crop varieties, crop diversification, and improved crop varieties is an another dominant strategy that is found to be used by smallholder farmers in adapting to the negative effects of climate variability and change as well as resultant changes in crop pest and disease pressures. Improved varieties (drought-tolerant varieties and/ or short cycle) allow for increased productivity even during dry seasons (Lobell et al. 2008). Furthermore, Ellis and Freeman (2004) found that crop diversity is used as a strategy for risk avoidance due to sharp fluctuations in crop yield or prices.
It is also disclosed that because of the unreliable and erratic pattern of rainfall and repeated drought, farmers in the study area started using small-scale irrigation schemes on their farm. It is also observed that diverted streams, pond construction, and use of water pump are found to be the major means for small-scale irrigation in the study area (Amare and Simane 2017b). In response to the adverse effects caused by climate variability and change, smallholder farmers have been diversifying their sources of livelihood with an understanding of more diversified livelihood strategies lead them both to enhance incomes and spread the risk for smallholder farmers. This includes the use of on-farm, off-farm, and non-farm activities (Amare and Simane 2017b;Morton 2007). However, despite large investment in adaptation options in response to the adverse effect of recurrent drought and floods, households in the study area continue to suffer from food insecurity.
Studies have been undertaken to measure the impact of climate change on Ethiopia agriculture and explored possible adaptation options in response to its adverse effect (Deressa 2007;Kidane et al. 2006;NMSA 2001). Insights from these studies are crucial in appreciating the extent of the problem and designing appropriate adaptation options. Notwithstanding the upsurge in the promotion of such adaptation options in recent years, there have been limited empirical studies that attempted to analyze their impacts on household food security (Ali and Erenstein 2017;Asfaw et al. 2015;Gebrehiwot and Anne Van Der 2015). Farmers adopting more adaptation practices had higher food security levels than those who did not (Ali and Erenstein 2017). Similarly, the analysis by Gebrehiwot and Anne Van Der (2015) using food security package program that has been implemented in Tigray regional state as an adaptation option to changing climate showed that the program has had a significant effect on improving household food calorie intake. Using a survey conducted in Ethiopia, Asfaw et al. (2015) showed that the adoption of climate-smart agricultural practice have positive and significant impacts on the objective measure of food security (net crop income) but no impact is observed for the subjective food security indicator.
However, there is inadequate evidence to what extent that adoption of those adaptation options impacted household food security in Ethiopia in general and in Muger sub-basin in particular. The results of these previous studies are highly fragmented and are of little help in addressing local conditions in relation to adaptations to climate change. Moreover, the studies overlooked the likelihood household food security impact of adaptation options at the household level. Using household-level data collected from a random cross-sectional sample of 442 farmers in the study area, the aim of the current study is to provide a comprehensive analysis on the impact of adaptation options of climate change and variability in the Muger sub-basin of the upper Blue-Nile of Ethiopia. The specific objective of the paper is, therefore, to estimate the effect of adoption of soil and water conservation measures, small-scale irrigation, agronomic practices, and livelihood diversification strategies as climate change adaptation options on household food security status measured by household calorie intake/day per adult equivalent using propensity score matching techniques. Addressing this question provides empirical evidence on the importance of adaptation options in improving household food security. Furthermore, this study provides important insights and lessons on the importance of access to resources on the ability of the farm households to invest in climate change adaptation options that latter improve household's capacity to adopt adaptation options to changing the climate. In addition, the finding of the study can be used in designing better adaptation interventions that can accommodate the existing resource potentials.

Bio-physical setting
Muger sub-basin is part of the upper Blue-Nile basin and cover a total area of 8188 km 2 . Muger River flows from the southeast of the basin into Abbay River. The altitude in Muger sub-basin ranges between 953 and 3550 masl.
The highlands in the eastern and southern part of the sub-basin are higher in altitude, greater than 2600 m up to 3550 m. The lowlands along the Muger River have lower altitude less than 1700 masl (Denekew and Seleshi 2009).
Following Deressa et al. (2009) classification of agroecological zones of Ethiopia, Kolla agroecology (lowland) is characterized by relatively hotter and drier climate, whereas Weyina Dega (middle land) and Dega agroecology (Highland) are wetter and cooler. Using this classification, the sub-basin of the study area contains Kolla, Woyina Dega, and Dega agro-ecological zones (Fig. 1). Evidence revealed that farmers living in different agroecological settings have their own choice of adaptation methods (Deressa et al. 2009;Legesse et al. 2013;Tessema et al. 2013). Moreover, farmer's vulnerability to climate change varies among different agro-ecological zones because there exists spatial heterogeneity among agro-ecological zones in terms of a varying level of socioeconomic and infrastructure development, households access to resources, level of food insecurity, and the ability to cope (Amare and Simane 2017c). Diversity in topography, socio-economic, and environmental issues determine the types of adaptation options farmers used in response to climate change.
The sub-basin has an annual rainfall that varies between 833 and 1326 mm. Lower annual rainfall ranging from 833 mm up to 1000 mm is observed along the river and lowlands. Relatively high rainfall is found in the highlands of the sub-basin. The annual maximum and minimum temperature of the sub-basin varies between 16-31.5°C and 3-16.5°C, respectively. Temperature is higher along the river with a maximum of 28-31.5°C and minimum of 13-16.5°C. The sub-basin is characterized by tepid to cool moist highlands. The northwestern part of the lowlands is hot to warm moist lowlands (Denekew and Seleshi 2009).
Besides the biophysical characteristics of the sub-basin, Figs. 2 and 3 present the trend in the annual maximum and minimum temperature data from the year 1990 to for the three metrology stations found in the sampled districts namely Gebereguracha station for the Kuyu district, Sululta station for the Sululta district, and Kachisi station for Abuna Gindeberet district, respectively. Similarly, Fig. 4 illustrates the trend in the annual precipitation in the selected stations. The mean annual temperature and rainfall records over the period under consideration were computed from the daily temperature and rainfall records obtained from each metrology station. The rainfall trend suggests a fluctuating and general decline in rainfall values over the study area. The results show a noticeable decline is observed in Sululta district in rainfall values from 1991 to 1992 and from 2006 to 2007. The result further shows that a noticeable decline in rainfall is observed from 1997 to 1998 in Kuyu district. It is this variability that created a decline in crop production and productivity and/or a total crop failure due to insufficient rainfall during the production seasons in the study area. This problem is more prevalent in Abuna Gindeberet district where most of the area of the district fall under kolla agro-ecological zone. Although the mean of maximum temperature at Abuna Gindebert and Sululta districts is more or less constant over years, an increasing trend is observed at Kuyu district between the years 1996 and 2000.
The major soils of the sub-basin are Leptosols, Luvisols, Vertisols, Fluvisols, and Alisols. Leptosols represents the most widely occurring soils within the sub-basin. The second dominant soil is Luvisols. Small patches of Cambisols, Nitosols, and Rigosols are also found in some parts of the sub-basin.
The Muger sub-basin is predominantly rural in character, and the farmers are engaged in small-scale and subsistence mixed agriculture. The dominant sources of livelihoods in the sub-basin are land and livestock. Muger sub-basin falls under three major livelihood zones including (1) Abay-Muger-Jemma Sorghum and Teff Belt Livelihood Zone where its category of food economy is mixed farming and its basic economic activities lie on crop and livestock production; (2) Ambo Selale Ginde-Beret Teff and Wheat Livelihood Zone which have mixed farming system and the main economic activities are a crop, livestock, and livestock product production; (3) Selale-Ambo Highland Barley, Wheat and Horse bean Belt Livelihood Zone with an economy characterized by the cultivation of rain-fed cereals complemented by livestock rearing and sales of trees and fodder (Denekew and Seleshi 2009).

Sampling procedure
The research design was based on multi-stage sampling procedure. In the first stage, the whole sub-basin constituting 15 Districts were grouped into three strata (Kolla, Woyina Dega, and Dega agro-ecological zones) based on their agro-ecological characteristics including the rainfall, soil, and topography. Then, Abuna Gindeberet district and Sululta district were randomly selected from Kola and Dega agro-ecological zones respectively. Similarly, two districts, namely Meta Robi and Kuyu Districts were also selected from Woyina Dega agro-ecology using simple random sampling technique. In the second stage, only PAs found in the sub-basin in each sampled District were listed in consultation with agricultural experts in the area. This is mainly to exclude PAs which are not part of the sub-basin in that particular District. Then, four PAs were randomly selected from each selected district. The third stage was a random selection of 450 respondents from the sampling frame of all the selected PAs using random sampling technique on the basis of probability proportional to size (PPS). The sampling frame was the list of households which was obtained from the PAs administration. Households for Focused Group Discussions (FGD) were also drawn from each identified district, and a member of the group was identified with the help of development agents working in the area.

Types of data and methods of data collection
Focused Group Discussions (FGD) and household survey were used to collect qualitative and quantitative data, respectively. Data from the FGDs' were used to complement the information obtained through a household survey in order to have a better understanding of causes of food insecurity and the causal influence of different adaptation options on food security. There were eight focus-group discussions held in three agroecology of which four were in Woyina Dega and two in each of the other two agroecology. The focus groups were composed of six to eight elders. One focus group in each of agroecology was composed of women who were perceived to have a deep knowledge of food security challenges and impacts of climate change on food security.
Quantitative data were collected using household survey comprised the same sets of questions. The dataset consisted of food security variables (total grain produced, total grain purchased, total grain obtained through FFW, total relief food received, total crop used for seed, total marketed output); adaptation options (SWC, small-scale irrigation, agronomic practices, livelihood diversification strategies); asset ownership (landholding, livestock ownership); social capital (local institutions/organizations households are membership); and human capital and access to financial capital. The survey also covered data on several factors including households' demographic and socioeconomic characteristics (age, education, gender, and family size) as well as environmental variables such as access to the early warning system, the frequency of occurrence of drought, and their experience's in crop failure due to climate change and variability.

Treatment variable
The treatment variable is based on the question "what adaptation options are households used to ameliorate The results from our analysis showed that soil and water conservation measures, small-scale irrigation, agronomic practices, and livelihood diversification are found to be the dominant adaptation options that farmers used in response to climate change and variability.

Outcome variable
The outcome variable that we used in this analysis is household food security status measured in household food calorie intake/day per adult equivalent. We used a household food balance model to compute the status of household food security using total grain produced by household; total grain purchased by household; total grain obtained through Food for Work (FFW) by household; total relief food received by household; postharvest losses to household; total crop utilized for seed from the home by the household; and total marketed output by household in a year of the survey period. The method is originally developed by FAO (1996) and later modified and used by many authors to measure food security status of households (Mesay 2010;Eshetu 2000;Degefa 1996).

Methods of data analysis Impact estimation strategies
If the adaptation options were randomly assigned to farmers, we could assess the impact of their adoption on households' food security by comparing the average consumption of adapters and non-adapters. In such a case, the average treatment effect (ATE) can be computed as follows: This is based on the assumption that the output levels of the adapters before their adoption (E (Y0|D = 1) can reasonably be approximated by the output level of nonadapters during data collection (E (Y0|D = 0). Otherwise, estimation of ATE using the above equation is not possible since we do not observe E (Y0|D = 1) though we do observe E (Y1|D = 1) and (E (Y0|D = 0). However, adaptation options are rarely randomly assigned. Instead, adoption usually occurs through self-selection of farmers or, sometimes, through program placement. In the presence of self-selection or program placement, the above procedure may result in a biased estimation of the impacts of adaptation options since the treated group (i.e., the adapters) are less likely to be statistically equivalent to the comparison group (i.e., the non-adapters) in a nonrandomized setting.
The literature shows various methods to address this selection bias. Some studies have adopted the Heckman two-step method that assumes a normal distribution of unobserved variables. Another method employs instrumental variables (IV). This approach usually requires at least one variable in the treatment equation to serve as an instrument for the specification of the outcome equation. Moreover, both ordinary least square (OLS) and IV procedures restrict the model to take a linear functional form, implying that the coefficients on the control variables are similar for treatment and control groups (Ali and Abdulai 2010).
In the absence of experimental data, the PSM is a widely used method to account for this sample selection bias, which is also employed in this study. The PSM technique pairs the treatment (adapters) and control (non-adapters) groups based on the similarity of observable characteristics (Ali and Abdulai 2010). Unlike the OLS and IV techniques, the PSM technique relaxes the assumptions of functional form, normal distribution of unobserved covariates and finding instrumental variables for the specification of the outcome equation. It only requires a set of observable covariates for matching and to determine causal effects of treatment on the outcome variable (Heckman and Vytlacil 2007). One limitation of PSM is that it does not account for the unobservable variables directly; rather, it assumes that selection is based on observable variables. PSM can be a better choice when instruments are weak or not available (Ali and Abdulai 2010). PSM is defined as the conditional probability that a farmer adapts to climate change, given the pre-adaptation characteristics (Rosenbaum and Rubin 1983).
In order to measure the impact of adaptation options (soil and water conservation, agronomic practices, livelihood diversification, and small-scale irrigation) in response to climate change and variability, this research used Propensity score matching technique. The first step in estimating the treatment effect is to estimate the propensity score. To get this propensity scores, any standard probability model can be used (for example, logit, probit or multinomial logit) (Rajeev et al. 2007). Since the propensity to adoption is unknown, the first task in matching is to estimate this propensity. Any resulting estimates of program effect rest on the quality of the adoption estimate. This can be routinely carried out using a choice model. Which choice model is appropriate depends on the nature of the program being evaluated. If the program offers a single treatment, the propensity score can be estimated in a standard way using, for example, a probit or logit model, where the dependent variable is "adaptation" and the independent variables are the factors thought to influence adaptation and outcome (Getachew et al. 2011b).
After obtaining the predicted probability values conditional on the observable covariates (the propensity scores) from the binary estimation, matching was done using a matching algorithm that was selected based on the data at hand. Then, the effect of household's adaptation in response to climate change on a given outcome (outcome in this study is food security measured in food calorie intake) (Y) is specified as: where τ i is treatment effect (effect due to adaptation), Y i is the outcome on household i, and D i is whether household i has got the treatment or not (i.e., whether a household adapted to climate change or not). However, one should note that Y i (D i = 1) and Y i (D i = 0) cannot be observed for the same household at the same time. Depending on the position of the household in the treatment (adoption of adaptation options), either Y i (D i = 1) or Y i (D i = 0) is the unobserved outcome (called counterfactual outcome). Due to this fact, estimating individual treatment effect τ I is not possible and one has to shift to estimate the average treatment effects of the population than the individual one. Most commonly used average treatment effect estimation is the average treatment effect on the treated (τ ATT ) and specified as: As the counterfactual mean for those being treated, E[Y(0)/D = 1] is not observed, one has to choose a proper substitute for it in order to estimate the average treatment effect (ATT). One may think to use the mean outcome of the untreated individuals, E[Y(0)/D = 0] as a substitute to the counterfactual mean for those being treated, E[Y(0)/D = 1]. However, this is not a good idea, especially in nonexperimental studies, because it is most likely that components which determine the treatment decision also determine the outcome variable of interest.
In this particular case, variables that determine household's decision to adopt could also affect household's food security. Therefore, the outcomes of individuals from treatment and comparison group would differ even in the absence of treatment leading to a self-selection (Dehinenet et al. 2014;Getachew et al. 2011a). Subtracting E[Y(0)/D = 0] from Eq. 3 yields the following specification for ATT: when there is no self-selection bias). In non-experimental studies, one has to introduce some identifying assumptions to solve the selection problem.
The validity of the outputs of the PSM method depends on the satisfaction of two basic assumptions, namely, the Conditional Independence Assumption (CIA) and the Common Support Condition (CSC) (Becker and Ichino 2002). CIA (also known as Confoundedness Assumption) states given a set of observable covariates (X) which are not affected by treatment (in our case, adaptation), potential outcomes (food security measured by food calorie intake) are independent of treatment assignment (independent of how adaptation decision is made by the household). This assumption implies that the selection is solely based on observable characteristics, and variables that influence treatment assignment (adaptation decision is made by the household) and potential outcomes (food security) are simultaneously observed (Getachew et al. 2011b). The common support condition entails the existence of sufficient overlap in the characteristics of the treated and untreated units to find adequate matches (or a common support).
Three commonly used matching algorithms, namely nearest neighbor matching, radius matching, and kernelbased matching, were employed to assess the impact of adaptation options (soil and water conservation, water harvesting, and small-scale irrigation) on households' food security. The nearest neighbor matching (NNM) method matches each farmer from the adopter group with the farmer from the non-adopter group having the closest propensity score. The matching can be done with or without replacement of observations. NNM faces the risk of bad matches if the closest neighbor is far away. This risk can be reduced by using a radius matching (RM) method, which imposes a maximum tolerance of the difference in propensity scores. However, some treated units may not be matched if the dimension of the neighborhood (i.e., the radius) is too small to contain control units. The kernel-based matching (KM) method uses a weighted average of all farmers in the adopter group to construct a counterfactual. The major advantage of the KM method is that it produces ATT estimates with lower variance since it utilizes greater information; its limitation is that some of the observations used may be poor matches.

Adaptation options
The study explored four dominant adaptation options used by smallholder farmers in response to the adverse effects of climate change and variability in the study area. These include soil and water conservation (28.7%); small-scale irrigation (27.4%); livelihood diversification (14%); and agronomic practices (11.5%). The rest 18.8% of households are found to be non-adopter of the adaptation options. To evaluate the total impact of these adaptation options on household food security calorie intake, we created a new data set by driving a binary variable equal to 1 if households adopted either of the adaptation options and 0 if households did not adopt any of the adaptation options. Based on this categorization, we found that 81. 2% of the households are falling under adopter group, and the 18.8% of the households are non-adopter.

Comparison of adopters and non-adopters of adaptation options
Tables 1 and 2 present the t test and chi-square comparison of means of selected variables between adopters' categories for adaptation options, respectively. Results show that positive and significant difference is observed in the age of household head, family size, social capital, land holding, distance to the nearest market, livestock holding, and a number of drought events for adopters of adaptation options. Moreover, adopters of adaptation options are also distinct in terms of agroecology, the gender of the household head, access to training, and access to credit.
The results also reveal that the adopter group is also significantly distinguishable in terms of food security status, measured in terms of household food calorie intake/day per adult equivalent. As far as calorie intake is concerned, adopters of adaptation options tend to have higher calorie intake as compared to the non-adopters. The mean calorie intake/day per adult equivalent is 4508 and 1549 kcal for adopters and non-adopters, respectively (Table 1). The t test result also shows that there exists a significant and positive difference between the two adopter groups in their calorie intake/day per adult equivalent.

Determinants of adoption of adaptation options
Given that variety of differences exist between adopters and non-adopters of adaptation options, it is important to control for these potential underlying effects in order to ensure reliable impact estimates. In order to provide information on household's probability of adopting adaptation options on household food security, a probit model is used to match adopter and non-adopter households. The results of the probit model presented in Table 3 depict estimated parameters on the adoption of climate change adaptation options.
The results from the probit model estimates indicate that gender is positively and significantly related to the adoption of adaptation options at 1% significance level. The result is in conformity with the previous argument by showing that male-headed households had better opportunity to take an adaptation measure than female household mainly due to cultural and social barriers that limit women's access to land and information climate change (Deressa, Yehualashet, Rajan, 2014;Asfaw and Admassie, 2004). The model result of this study reveals that size of cultivated land significantly increases the likelihood of using adaptation options at 5% significance level. On average, the probability of using adaptation options increases by 10.7% as the proportion of farm size increases by 5%. Furthermore, access to extension services has a significant influence on adoption of adaptation options. As expected, the influence of access to extension service on farmer's decision to invest in adaptation is significantly positive (p < 0.000). Having access to extension service increases the probability of adoption of adaptation options by 22.8%.
Better access to early warning information about drought and flood before it happened has a significant and positive impact on the likelihood of using different adaptation options at 10% significance level. The result reveals that getting access to climate warning about drought and/or flood increases the likelihood of using adaptation option by 5.04%. This finding is in line with results of previous studies that showed a positive relationship between early warning information and farmer's decision to adopt adaptation option (Deressa, Yehualashet, Rajan, 2014;Phillipo, Bushesha, Mvena, 2015). On the contrary, livestock holding in TLU negatively influences household's decision to adopt adaptation options at 5% probability level. This result reveals that a unit increase in a number of livestock in TLU would result in a 3.35% decrease in the probability of adopting adaptation options. Finally, the results from the probit model show that farmers decision to adopt adaptation options is not significantly influenced by the variables agroecology, access to credit, distance to the nearest market, social capital, education status, and age of the household head.

Estimation results of propensity scores
After calculating the propensity scores, the nearest neighbor matching (NNM) method, Kernel Matching, Stratification Matching, and Radius Matching were employed to match the control group of individuals (non-adapters) to the treated group (adapters) based on similar propensity scores. During the matching process, all the matching methods employed discard the unmatched non-adapters, and hence, they lead to the reduction in sample size for the post-matching impact analysis. The region of common support is [.01697335, .99999403] and the balancing property also satisfied. The results confirmed there exist a considerable overlap in common support. Figure 5 shows the distribution of propensity scores of matched and unmatched individuals in both groups. The result guarantees a sufficient overlap in the distribution of the propensity score between adopters and non-adopters. The bottom half of the graph stands for the propensity score distribution for the non-adopters and the upper half refers to the adopters. The densities of the scores are on the y-axis. In the next step, we calculated the average adaptation effects on the status of household food security (Table 4).

Estimation of the treatment effect
The correlation between adoption of adaptation options and household food security outcome such as household calorie consumption/day per adult equivalent is theoretically complex and there are further empirical pitfalls regarding the impact evaluation problem. We estimated the calorie consumption effect of climate change adaptation options based on cross-sectional data available. A propensity score matching technique was used to address the research questions. Table 4 presents the estimated effects of the adaptation options on the household food security by nearest neighbor matching (NNM), Kernel-based matching (KBM), and stratification matching methods. The post-matching result from nearest neighbor matching reveals that adaptation tends to positively and significantly affect household calorie intake. The estimates of nearest neighbor matching show that adoption of adaptation options improves household food calorie intake by 2635 kcal/day per adult equivalent. Similarly, the Kernel matching result reveals adopting either of the adaptation options also guarantee favorable effect on food security. This  means that compared to similar/matched households that did not adopt any of the adaptation options, households that did adopt experienced a 2660 kcal/day per adult equivalent higher number of calorie intake (significant at the 1% level). This implies that adoption of adaptation options has a causal influence on household food security status. This indicates that the return to investment in soil and water conservation measures; small-scale irrigation; agronomic practices; and diversifying livelihood options each does generate reliable results, especially in areas where climate change and variability adversely affects agriculture which is considered as the main livelihood option of the households. This reaffirms the narrative from the results of the focus-group discussions in which not adopting soil and water conservation measures such as stone bund, soil bund, check dam, and terraces, coupled with the lack of small-scale irrigation practices, undermines the prospect of food security. This is consistent with the secondary literature that shows a positive effect of technology adoption on food security (Amare et al. 2012;Magrini and Vigani 2014).

Conclusions
We have analyzed the effects of adoption of adaptation options on food security among smallholder farmers in Muger sub-basin of the upper Blue-Nile basin of Ethiopia. We used the non-parametric method involving the application of a binary propensity score matching estimator. Use of soil and water conservation measures, small-scale irrigation, agronomic practices, and livelihood diversification strategies is considered as an adaptation option and found to be the dominant climate change adaptation options used by the farmers in the study area. Adaptation decisions are significantly affected by various socio-economic, institutional, environmental, and demographic factors. In particular, gender, family size, land ownership, access to extension service, livestock ownership, and frequency of droughts and floods that have happened in the past 25 years are the important factors influencing farmer's decision to adopt adaptation options.
The results show that adoption of these climate change adaptation options has generated a significant positive impact on household food security. Farm households that did adopt would benefit the most from adoption. These results generate a strong evidence for the positive impact of adoption of climate change adaptations that limit the adverse effect of climate change on households' livelihoods on alleviating food insecurity in the study area. However, farmers are yet unable to enjoy the full benefit of adaptation options due to various constraints. These results open a scope for policy to further promote the adoption of climate-change adaptation options that seem to be particularly important for households who have the least capacity to adapt.
Firstly, the results highlight the importance of awareness and knowledge about adaptation options and its riskreducing potential. Hence, policy should target increasing the awareness of climate change adaptation options and its benefits using agricultural extension services. Secondly,  the result points to the importance of access to resources on the ability of the farm households to invest in climate change adaptation options. The policy option should give attention to resource-constrained farmers. Thirdly, the result indicates the importance of gender of the household heads to adopt adaptation options. Hence, policy should target gender-based adaptation to climate change and variability. All these policy implications may lead to better adoption of adaptation options and may be able to support farmers to ensure food security. On the basis of the results, one might deduce that investment in climate change adaptation options does generate substantial benefit in terms of achieving household food security because adaptation options increase food availability because they boost crop productivity and promote food stability by reducing the risk of crop failure, making the household less vulnerable to negative shocks, and improving the resilience capacity. However, such an aggregation does not provide adequate evidence to make firm conclusions. This paper posed a concern to see the potential effect of each adaptation options on household food security and to identify the most successful one for further research. Lastly, in addition to the policy relevance of these findings to the study area, it could be replicated or applied in the area where similar problems might prevail, particularly in other developing countries.