Introduction
Demand for livestock products in Sub-Saharan Africa has been increasing and is projected to continue increasing due to population growth, rising incomes and urbanization (Thornton et al. 2007; FAO 2009; Thornton 2010; Robinson and Pozzi 2011; Ghimire et al. 2015). Supply of livestock products has not kept pace with demand, due primarily to low productivity and limited land area (Rakotoarisoa et al. 2011). Production of livestock products is further complicated by climate change (Thornton et al. 2007; Thornton 2010). One of the major factors behind the region's chronic low productivity is a lack of quality feed options with high nutrient content. Producers in mixed, rain-fed crop-livestock systems are particularly constrained by a shortage of feed resources during dry seasons, a situation that is systematically aggravated by increasing pressures from climate change and variability (Dzowela 1990; Thornton 2010; Rakotoarisoa et al. 2011).
Better use of the natural resource base offers tremendous potential to increase livestock production in the region (FAO 2009; Ghimire et al. 2015). Research programs such as 'Climate-smart Brachiaria' have begun developing strategies to tap into this potential (Djikeng et al. 2014). Such efforts are built around the development of drought-resistant Brachiaria grass varieties with climate change-mitigating properties (Ghimire et al. 2015; Maass et al. 2015). Building on our earlier work (González et al. 2016), in this study we present an ex-ante assessment of the potential welfare impacts of increasing milk production by introducing such technology to mixed rainfed crop-livestock systems (Table 1; Figure 3A) in Eastern Africa, using the economic surplus method previously described by Alston et al. (1995).
Brachiaria technology and milk production
The genus Brachiaria consists of roughly 100 species which grow in the tropics and subtropics. Most of these species are native to Africa, where they constitute important components of the natural savanna landscape (Ghimire et al. 2015). Outside of Africa, widespread commercial adaptation and adoption of Brachiaria species in non-native environments has enhanced livestock industries worldwide - notably in Latin America and the Caribbean, as well as in Asia and Australia - and has made Brachiaria the most extensively cultivated forage monoculture in the world (Jank et al. 2014; Ghimire et al. 2015).
The widespread appeal of Brachiaria lies in a diverse set of traits, depending on species and cultivar, including adaptability to infertile and acidic soils, resistance to drought, tolerance of shade and flooding and palatability. From an environmental perspective, it is also appealing because it transfers carbon from the atmosphere into the soil, makes efficient use of nitrogen, and helps to minimize groundwater pollution (Fisher et al. 1994; Fisher and Kerridge 1996; Rao et al. 1996; Subbarao et al. 2009; Rao 2014).
The success of Brachiaria in other parts of the world has motivated concerted efforts to introduce higher-performing, improved cultivars in Africa. The Brachiaria hybrids developed at CIAT over the course of the 1980s and 1990s for release in the Americas (Mulato and Mulato II) have been introduced to several African countries on an experimental basis since 2001. Limited uptake and diffusion of these hybrids has occurred through farmer-to-farmer transfer of planting material promoted by research programs (Maass et al. 2015). Much of this diffusion is associated with the spread of 'climate-adapted push-pull' farming systems (Midega et al. 2015). Based on estimates of seed sales, as of 2014, some 1,000 hectares of these hybrids were under cultivation in various African countries, primarily in East Africa (Maass et al. 2015).
While initial results of crop trials demonstrate a potential for positive return on investment (Kabirizi et al. 2013; Ghimire et al. 2015), these hybrids were developed specifically in response to biotic and abiotic stresses in Latin America. Their introduction in Africa has encountered biotic challenges which must be overcome before adoption and diffusion can be scaled up significantly (Maass et al. 2015).
A Swedish-funded program called 'Climate-Smart Brachiaria Grasses for Improving Livestock Production in East Africa' (referred to as CSB) is addressing these challenges (Djikeng et al. 2014; Ghimire et al. 2015). The program is led by the Biosciences Eastern and Central Africa-International Livestock Research Institute (ILRI) Hub, and is in partnership with the Kenyan Agricultural and Livestock Research Organization, the Rwanda Agricultural Board, CIAT and Grasslanz Technology Limited. The program is currently implemented in Kenya and Rwanda, with planned future expansion in Eastern Africa and beyond.
In advance of the CSB program, 10 Brachiaria cultivars - mostly from the brizantha species, but also including the hybrids Mulato and Mulato II - were tested in greenhouses at CIAT in Colombia against Eastern African baseline varieties such as Rhodes and Napier grass. Results were encouraging and, beginning in 2013, 8 of 10 cultivars were selected for field trials at multiple sites in Kenya and Rwanda. Of these 8 cultivars, B. brizantha cvv. Piatã, Marandu, La Libertad (also known as MG-4) and Toledo (also known as Xaraés), B. decumbens cv. Basilisk and the hybrid Mulato II emerged as the best performing varieties. Mulato II and Marandu were subsequently removed from trials after they proved susceptible to local pest infestation. On-farm evaluation of the remaining 4 cultivars began in 2014 and is ongoing at the time of this study (Ghimire et al. 2015; CSB 2016).
CSB trials also included a focus on the quality of the grass as animal feed. Preliminary data from recent trials indicate that adoption of these mostly B. brizantha cultivars has the potential to increase baseline milk production of 3-5 L/cow/d on participating farms in Kenya by 15-40%. A farm trial in Rwanda reported a 30% increase in milk production and a 20% increase in meat production (CSB 2016). No meat production data were available from the Kenya trials.
Brachiaria grasses tend to be drought-resistant and resilient in infertile soils, and produce well with relatively low levels of fertilizer inputs. They are also resistant to many diseases affecting baseline varieties in Eastern Africa, particularly Napier stunt and smut disease (Ghimire et al. 2015; Maass et al. 2015). Brachiaria production can be further enhanced by intercropping with deep-rooted, nitrogen-fixing legumes such as Centro and Clitoria (Kabirizi et al. 2013), which themselves are useful sources of nutrition for animals.
Though the dry matter yields of the Brachiaria forages under evaluation tend to be lower than those of baseline varieties found in Eastern Africa, their leaf areas are relatively larger, effectively increasing palatability and nutrition per unit of dry matter. The protein concentration of Brachiaria, ranging from 8 to 17% at harvest, remains stable for a relatively long time as compared with that of baseline varieties, where protein concentration diminishes after about 4 months (CSB 2016). Surplus Brachiaria not immediately consumed can thus be dried and conserved as hay for sale or future use.
The advantages and disadvantages of improved Brachiaria grasses relative to baseline varieties also tend to vary seasonally. While Brachiaria outperforms baseline varieties during dry seasons, the baseline varieties exhibit certain advantages during rainy seasons (Kabirizi et al. 2013). On many farms, it may make sense to introduce the improved Brachiaria grasses as a dry season complement to the baseline grasses. Kabirizi et al. (2013) point out that small farms, which introduce Brachiaria in such a complementary role, would probably have to displace a cash crop in order to make room for the new addition, and should thus consider the opportunity cost in terms of potential forgone revenue from the cash crop.
As of May 2016, at least 4,000 farmers in Kenya and Rwanda had reported planting one of the Brachiaria cultivars under CSB evaluation (CSB 2016). Experts at CIAT report that participating farmers appear to prefer B. brizantha cv. Piatã of the 4 cultivars currently under CSB evaluation (J. A. Cardoso pers. comm.). The already substantial numbers of farmers using the technology and the corresponding return on investment and increased resilience for the forage systems suggest that there is substantial potential for impact of these forages in Eastern Africa. Using data collected from a number of sources, we evaluate, ex-ante, the potential impact of improved forages in the region.
Modelling the plausible outcome space
In every ex-ante impact study, there is an implicit tradeoff between the accuracy of model parameterization and the time and budget within which this can reasonably be accomplished. In the vein of demand-driven modelling (Antle et al. 2017), we aim not to maximize accuracy, but rather to maximize accuracy subject to the constraints and needs of the stakeholders motivating the study. These stakeholders include a variety of public and private sector actors, all of whom are ultimately motivated by the needs of smallholder farmers who are the end users of the research product. Considering that order-of-magnitude accuracy is often a sufficient premise on which to base policy decisions, and that stakeholders need assessments of potential impact in a timely manner, we take a parsimonious approach based on existing data and consultations with regional experts. We present our modelled outputs not as a final conclusion, but as a map of plausible outcomes intended to aid readers in their navigation towards a conclusion based on their own understandings of forage systems. We further distill this map into a single envelope equation by which the reader can easily generate model outputs for any level of impact, adoption rate and production increase he/she wishes to consider. Finally, we conduct sensitivity analysis on several key parameters.
The Model
When assessing return on investment in research products, the whole process from research outcomes through release and uptake of the new agricultural technology must be considered. The economic benefit for each country in the study area is thus defined as the net present value (NPV) of the cost-benefit stream extending from year one of research up to the point where the adoption ceiling is reached. Program-level costs occur from the initial year of research until release of the new technology. Subsequent costs associated with production of planting materials, marketing and distribution are typically incurred by private sector actors and thus excluded from the calculation, although we do account for minimal country-level diffusion costs incurred by public sector actors from the year of release over an initial phase of adoption.
The tool we use to calculate the benefit stream is Alston et al.'s (1995) economic surplus model for closed economies. This model, summarized in Figure 1, measures benefits in a given year as the increase in total surplus resulting from a research-induced shift in the supply curve for a given commodity of interest (the shaded area). The total surplus can be divided further into benefits accruing to producers (producer surplus, the shaded area above line P 1 b ) and benefits accruing to consumers (consumer surplus, the shaded area below line 𝑃1𝑏).
The commodity of interest for this study is fresh cows' milk. We evaluate one such model for each country in the study zone, for each year from release of the technology to the year of maximum adoption. The markets are said to be 'closed' because we assume no cross-border trade of fresh milk. Note that these benefit streams understate the true benefit to some degree since they take no account of positive impacts on meat production, which as stated earlier were 20% increases in Rwanda.
Model parameters
In order to calculate the cost stream and the total surplus represented by the shaded area in Figure 1, we require as input the parameters in Table 2. As in most economic surplus studies, estimates of the supply and demand elasticities for the precise commodity and geographical area in question are difficult to acquire. We set the milk supply and demand elasticities to 0.7 and -0.5, respectively, in accordance with an estimate for all of Sub-Saharan Africa obtained by Elbasha et al. (1999). The research time horizon, annual research cost and depreciation factor were set based on consultation with a breeding expert (M. Peters pers. comm.). Based on the success of past CIAT forage research programs for release in other parts of the world, we feel justified in setting the probability of success at 80%. We set the interest rate at 10% to reflect the opportunity cost of not investing the research funds in a stock portfolio of comparable risk.
As discussed earlier, preliminary trial results suggest that adoption of the new technology can increase cow milk production by 15-40%. Another key advantage of the improved varieties is that they are robust on infertile soils, which implies a decrease in the variable costs associated with fertilizer applications. We assume that this potential cost decrease will be insignificant in Eastern Africa, where fertilizer use is already notoriously low.
On the other hand, as mentioned in the same section, many smallholder farmers who introduce the new technology in a complementary role may have to displace a cash crop, thus incurring an opportunity cost in the form of forgone revenue. However, the new technology is most likely to appeal to mixed rainfed smallholder systems within the study zone, where soils are marginal and where opportunity costs are, consequently, low.
For this study, we therefore assume that, on average, the potential variable cost reductions and opportunity costs associated with the new technology would either be negligible or offsetting, resulting in a percentage change in variable costs equal to zero.
Fixed capital costs associated with adoption of the new technology are not accounted for in this model.
After release of the new technology, it is typically acquired by a private sector actor which then accepts any subsequent costs associated with marketing and diffusion. We exclude these costs from our calculation of NPV since they are not incurred by the research institution nor governments. Nonetheless, as a conservative measure, we do include a minimal yearly diffusion cost to public sector actors for the period of initial release and uptake, modelled as a marginally diminishing function of the target industry size. The target industry size is measured as the number of cattle in the country's mixed, rainfed crop-livestock systems (N MR ). The parameter values 0.10 and 0.97 in Table 2 are chosen because they generate diffusion cost magnitudes commensurate with the types of promotional, training and outreach activities that are typical of country-level diffusion efforts in the study area. The diffusion cost magnitudes produced by this formula are presented for each country in Table 3. Though diffusion costs reflect an approximate cost based on industry size, they do not specifically take into account the nuances of the technology adoption environment in each country.
Regional expert survey: The technology adoption environment
Technology adoption varies depending on a number of factors and local conditions. Adoption of the new Brachiaria technology is modelled using a logistic curve (see Alston et al. 1995 for details). This 2 parameter curve reflects the typical slow rate of adoption initially, followed by a period of rapid diffusion, and then a tapering off of uptake as the adoption rate ceiling is reached (Figure 2). The curve parameters are calculated based on the duration of the uptake period.
In order to assess local conditions influencing technology adoption, we sent questionnaires to regional experts in the study zone. The responses we received, summarized in Tables 4-6, confirm that the Brachiaria cultivars under evaluation are most likely to appeal to mixed, rainfed systems. They convey moderate optimism about technology uptake in these systems, but also acknowledge considerable impediments, e.g. access to finance, quality inputs and extension services and infrastructure, which may hamper diffusion and uptake of the new technology. For these reasons, rather than present results for a single rate, we present outcomes for all adoption rate levels (at 5% intervals), giving the reader freedom to examine the outcomes that seem most likely to him/her based on his/her own experience and interpretation of the survey responses made available here.
1No response.
2Respondent gave a verbal response - 'modest' - which we have interpreted numerically as 2.
3Respondent gave an actual adoption rate - 25% - which we have assigned a scale rating of 2.
Most respondents indicated a moderate to long uptake period, where the terms 'moderate' and 'long' are subject to a great deal of interpretation. Our interpretation for this study is that the overall uptake period, including the diffusion period, would last 20 years in all countries.
Producer prices
In addition to the parameters summarized above, contemporary producer prices are required in order to calculate the total surplus stream. These were obtained both from regional experts in the study zone and from FAO. While neither of these sources on its own offered complete price data for all countries involved in this study, together they provide a more robust picture.
FAO reports recent producer milk prices for Kenya, Ethiopia and Rwanda. For these countries, we used the average over 2010-2012, which is the most recent consecutive period for which FAOSTAT reports price data for all 3 countries.
Field experts provided price data for Tanzania, Ethiopia, Uganda and Rwanda. In order to be consistent with the prices obtained from FAOSTAT, we again use the 2010-2012 average for these countries, except Uganda. The Uganda respondent reported prices for only years 2013-2015, so the Uganda producer milk price is averaged over this period. Respondents reported prices in local currency per kilogram, so we converted these prices to USD per metric tonne (MT) using historical exchange rates retrieved for 15 June in each respective year.
For Rwanda and Ethiopia, price data were available from both FAOSTAT and local experts. In these cases we used the lesser of the 2 prices. No price data were obtained for Burundi from any source, so we set Burundi's producer price equivalent to that found in Rwanda.
Quantity of production affected
The final piece of information required for calculation of the total surplus area depicted in Figure 1 is the quantity of production affected by the new technology. This is the baseline production already occurring in areas where the new technology is likely to appeal to producers. The Brachiaria varieties under evaluation in the CSB program are expected to appeal primarily to producers in mixed, rainfed crop-livestock systems, where baseline forage varieties currently fail to generate a sufficient feed supply during dry seasons (An Notenbaert pers. comm.). Under the Seré and Steinfeld classification map in Figure 3A, these production systems are designated as MRA, MRH and MRT (see Table 1 for definitions) (Robinson et al. 2011). These systems are characterized by their small size and marginal soils.
Baseline cow milk production data are available from FAO at the country level, but the production system levels defined by Seré and Steinfeld cut across national boundaries. In order to obtain a baseline production figure for each production system within each country, we first calculate the number of cattle within each system within each country by overlaying a production system map (Figure 3A) onto the latest available cattle density map (Figure 3B) to give the numbers presented in Table 8. We then generate modelled estimates of milk production for each system within each country as a function of total cattle based on the empirical relationships observed in Figure 4.
Figure 4A suggests that a log-linear relationship exists between total cattle and production, but that the y-intercept varies by region. This is drawn out more explicitly in Figure 4B, where regions are plotted separately.
Plots for other years in the FAO database exhibit the same log-linear relationship. In Figure 5, we see that the parameter values for this relationship are stable - albeit over the time periods 1961-2001 and 2006-2014, with a transition period in between . For a reasonable approximation, we conclude that, for a given region, the following scale invariant relationship exists between milk production (P) and total cattle (N)
where: for the Sub-Saharan Africa region, the mean values of a and 𝛽 over 2006-2014 are 1.23 and -6.563 (with standard deviations 0.01 and 0.152), respectively.
Since the relationship is scale invariant, we then apply this model (Equation 1) to the 2010 calculated numbers of cattle per production system within each country (Table 8) to determine milk production at the production system level. We fit parameters a and 𝛽 for each country such that they are close to their region-wide means of 1.23 and -6.563 above, and such that the total production in each country adds up to within 10% of the corresponding FAO 2010 country level totals. For most countries in the study zone, this results in values for a and 𝛽 that fall within 2 or 3 standard deviations of the region-wide means, although for Kenya and Rwanda the values are 4 standard deviations from the means (still reasonably close considering that the standard deviations are very small). These modelled approximations of baseline milk production at the production system level are presented for each country in Table 9. Finally, in each country we add up the modelled production in the mixed rainfed production systems. These figures (the 'MR Subtotal' in Table 9) represent the baseline production potentially affected by the new Brachiaria technology.
Results
Ex-ante approaches offer a forward-looking view of potential return on investment in an agricultural technology. The previous sections illustrate how the economic surplus model of Alston et al. (1995) can be parameterized, even in relatively sparse data environments. With the model parameterized, we can now populate the outcome map based on aforementioned adoption and benefit criteria.
Plausible outcomes map
Below we present NPV estimates based on a wide range of potential production increases resulting from adoption of the new Brachiaria technology in Eastern Africa (Figures 6-8). For each potential production increase, we also present results over the range of all possible adoption rates (0-100% at 5% intervals). Outcomes are calculated in terms of producer, consumer and total surplus. Each map cell is colored in accordance with the NPV it contains. Lower values are redder, higher values are greener; and the 50th percentile of NPVs is colored yellow.
NPV outcomes isoquant map and envelope formula
Results are also presented in an isoquant format in Figure 9. Analogous to isobars on a weather map or elevation contours on a terrain map, each isoquant represents an NPV outcome level, and each point on an isoquant indicates the production increase and adoption rate necessary to reach that NPV outcome. Equations fitted to these isoquants are of the form:
for the i th NPV isoquant, where A max is the adoption rate, E[y] is the expected increase in production resulting from adoption, and y¡ is a parameter to be fitted. This equation implies a one-to-one tradeoff between the adoption rate and the expected percentage increase in production. If the increase in production falls some percentage below expectations, the same level of NPV will still be achieved so long as the associated adoption rate is the same percentage above expectations.
Plotting the y¡ values against the log of their associated NPV values in Figure 9 reveals an interesting linear relationship (Figure 10) that permits us to reduce all possible NPV isoquants to a single formula (Equation 3).
This envelope formula encapsulates the model such that, for a given NPV outcome, the adoption rate (A max ) and expected change in production (E[y]) are allowed to vary, while the other parameters are held constant at their values in Table 2, encoded in the fitted parameters 0.744 and 6.775. Using this formula, the reader may determine the adoption rate necessary to achieve any given NPV outcome for any given percentage increase in production (or vice versa).
Sensitivity analyses
In any model, results may be sensitive to inaccuracies in input parameter values. It is therefore important to assess how sensitive the results presented above are to inaccuracies in key parameters, especially those parameters which are most uncertain. Sensitivity to fresh cow milk supply and demand elasticity values in particular warrant close scrutiny, as these were defined for all of Sub-Saharan Africa. In Figure 11, we present sensitivity analyses on these plus 2 other parameters.
In these sensitivity maps, an absolute value of 1 means that the NPV outcome for that scenario is as accurate as the parameter value. In other words, if the parameter value is off by 10%, then the NPV will also be off by 10%. Figures 11a and 11d indicate this kind of 1:1 model sensitivity to inaccuracy in the supply elasticity and producer price/quantity affected parameters for most scenarios, with sensitivity becoming extreme for a few of the low adoption scenarios on the fringe of the plausible outcomes space. Figure 11C indicates more moderate sensitivity to inaccuracy in the change in input cost parameter, and Figure 11B indicates very little sensitivity to inaccuracy in the demand elasticity.
Broadly speaking, the modelled NPV outcomes are about as accurate as the parameter values for supply elasticity, producer prices or quantity affected. The model is also moderately sensitive to inaccuracy in the change in input costs parameter. However, for a wide range of plausible scenarios, even a substantial inaccuracy in any single one of these would mean the difference between an 8th order result ($100s of millions) and a 7th order result ($10s of millions). Major inaccuracies would have to occur in several parameters simultaneously in order to critically skew the model output.
Discussion and Conclusions
The results of this economic surplus analysis suggest that investment in a research program involving the development of improved forage varieties for release in Eastern Africa would be a low-risk, high-reward endeavor. Preliminary data from ongoing multi-site trials in Kenya and Rwanda suggest that release and uptake of improved forages would increase milk production by 1540%. On a producer surplus basis alone, NPV outcomes are positive across this entire range so long as the adoption rate is at least 10%, and rise quickly into the tens of millions of dollars for a wide range of plausible adoption rates. When consumer side benefits are added in, the NPV outcomes are much greater still, reaching half a billion dollars for a wide range of plausible scenarios.
As far as the inner workings of the model are concerned, the overwhelmingly positive assessment is due in large part to the massive pool of potential beneficiaries in the study area (reflected in the baseline milk production), and because we assume there is no increase in input costs associated with adoption of the new technology. The relatively brief research period, compared with prior CIAT forage research programs, also contributes to this result.
When interpreting these results, it should be kept in mind that the economic surplus model employed in this study is a parsimonious, minimum data approach. This approach thus simplifies many important features of the underlying reality. In particular, we ignore any fixed capital improvements and other transition costs that might be associated with adoption of the new technology, e.g. terrain preparation, fencing, etc. The model employed in this study also makes no allowance for the often complex nature of land tenure in Eastern Africa, and the many ways this and other heterogeneous farm characteristics can vary across landscapes in the study zone. In other words, the model assumes that the percentage increase in production is the same for all adopting farms, regardless of variation in local conditions and factor endowments. Finally, we do not account for potential delays in diffusion due, for example, to production of planting materials by private sector actors subsequent to release of the research product. These simplifications in representation may bias our NPV outcomes upward, depending on the structure of the heterogeneity present in the region. We also assume that the supply and demand elasticities, adoption rate ceilings and uptake period durations are the same across all countries and across all production systems, although it is not clear in which direction these assumptions might drive the results.
On the other hand, our results are conservative in some respects. For example, we have taken no account of the additional benefits that might arise from increased meat production, enhanced production from associating the grass with a forage legume, the storage and/or sale of hay, the spread of climate-adapted push-pull systems, and potential multiplier effects on the broader economy.
The model results are presented in a heat map format that covers a broad range of potential outcomes, allowing the reader to compensate for the aforementioned potential biases by choosing an adoption rate consistent with his/her own level of optimism/pessimism regarding these sources of uncertainty, and with his/her interpretation of the regional expert opinions in Tables 5 and 6. The model envelope equation is also presented (Equation 1), whereby readers can calculate, for any given production increase that seems feasible to them, the modelled adoption rate required for a desired level of NPV (on a total surplus basis). This reporting format is intended to invite exploratory 'what-if questions and inter-comparison of scenarios which can be further refined with new data as they become available.