JDS
HOME HELP FEEDBACK SUBSCRIPTIONS ARCHIVE SEARCH TABLE OF CONTENTS
 QUICK SEARCH:   [advanced]


     


This Article
Right arrow Abstract Freely available
Right arrow Full Text (PDF)
Right arrow Interpretive Summary
Right arrow Alert me when this article is cited
Right arrow Alert me if a correction is posted
Services
Right arrow Similar articles in this journal
Right arrow Similar articles in PubMed
Right arrow Alert me to new issues of the journal
Right arrow Download to citation manager
Right arrow reprints & permissions
Citing Articles
Right arrow Citing Articles via HighWire
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Val-Arreola, D.
Right arrow Articles by France, J.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by Val-Arreola, D.
Right arrow Articles by France, J.
J. Dairy Sci. 89:1662-1672
© American Dairy Science Association, 2006.

Modeling Small-Scale Dairy Farms in Central Mexico Using Multi-Criteria Programming

D. Val-Arreola*, E. Kebreab{dagger},1 and J. France{dagger}

* School of Agriculture, Policy and Development, The University of Reading, Earley Gate, Reading, RG6 6AR, UK
{dagger} Centre for Nutrition Modelling, Department of Animal and Poultry Science, University of Guelph, Guelph, Ontario, N1G 2W1, Canada

1 Corresponding author: ekebreab{at}uoguelph.ca


    ABSTRACT
 TOP
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS
 DISCUSSION
 APPENDIX
 ACKNOWLEDGEMENTS
 REFERENCES
 
Milk supply from Mexican dairy farms does not meet demand and small-scale farms can contribute toward closing the gap. Two multi-criteria programming techniques, goal programming and compromise programming, were used in a study of small-scale dairy farms in central Mexico. To build the goal and compromise programming models, 4 ordinary linear programming models were also developed, which had objective functions to maximize metabolizable energy for milk production, to maximize margin of income over feed costs, to maximize metabolizable protein for milk production, and to minimize purchased feedstuffs. Neither multi-criteria approach was significantly better than the other; however, by applying both models it was possible to perform a more comprehensive analysis of these small-scale dairy systems. The multi-criteria programming models affirm findings from previous work and suggest that a forage strategy based on alfalfa, ryegrass, and corn silage would meet nutrient requirements of the herd. Both models suggested that there is an economic advantage in rescheduling the calving season to the second and third calendar quarters to better synchronize higher demand for nutrients with the period of high forage availability.

Key Words: small-scale dairy farm • multi-criteria programming • Mexico


    INTRODUCTION
 TOP
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS
 DISCUSSION
 APPENDIX
 ACKNOWLEDGEMENTS
 REFERENCES
 
In Mexico, national milk production has increased in the last 5 yr but demand has increased at a higher rate, resulting in about 35% of milk consumed being imported (Gallardo-Nieto et al., 2004). Such a situation within the North America Free Trade Agreement is causing the Mexican dairy sector to face major challenges to stay competitive. Small-scale dairy enterprises contribute between 35 and 40% of the national milk production. They represent an important source of employment and are enhancers of Mexican rural families’ livelihoods. They have several advantages that contrast with intensive systems: i) small-scale dairy enterprises have relatively low production costs; ii) they are highly integrated with crop production systems; and iii) they absorb most of the family labor and require a low investment in facilities (Arriaga-Jordan and Pearson, 2002; Wiggins, 2002). Several studies have pointed to small-scale dairy systems as an important means of increasing national milk production and at the same time being competitive with international standards (e.g., Espinoza-Ortega et al., 1997; Wiggins, 2002). New approaches are required to provide insights into improving the efficiency of this type of system. Previous work has suggested considering changes related to: i) selecting the most appropriate feeding strategy according to resource availability and herd needs (Val-Arreola et al., 2005); ii) deciding the number of cows calving in each period of the year (Val-Arreola et al., 2004a); and iii) selecting the best forage production strategy (Val-Arreola et al., 2004b). The previous work, however, considered different aspects of the system in isolation, unlike the current study, which considers the system holistically, leading to more specific recommendations.

Few studies related to modeling agricultural systems have been carried out in Mexico, particularly concerning dairy farms. Candler (1983) attempted to optimize a whole intensive dairy farm in northern Mexico using a linear programming (LP) model designed over a 12-mo horizon. The model was aimed at providing insights for government policy makers, but the detailed information required limited its applicability. Castelán-Ortega et al. (2003a) developed a simulation model to support the decision processes of small-scale dairy producers in central Mexico. However, the complexity of the model (3 submodels) and its size (more than 15,000 variables and 600 restrictions) makes the model matrix, besides being difficult to use and maintain, highly demanding of detailed information. Such information is difficult to obtain or is simply unavailable for this type of enterprise.

Ordinary LP has been used extensively for farm planning, but it is limited to dealing with a single objective (Piech and Rehman, 1993). Therefore, a multiple criteria decision-making (MCDM) framework should be more appropriate in dealing with the multidimensional nature of agricultural systems. Although application of such techniques is not common in agriculture, the techniques do allow more realistic models than ordinary LP (Rehman and Romero, 1993).

The primary aim of this study was to use the MCDM approach (considering objectives: maximize ME for milk production, maximize margin of income over feed costs, maximize MP for milk production, and minimize purchased feedstuffs) to assist the decision process in small-scale dairy farms in central Mexico. The decisions considered related to selecting the best forage production strategy capable of providing most of the nutrients required by the herd and determining the number of cows calving in each period of the year to match their requirements with nutrient availability through the most appropriate feeding strategy.


    MATERIALS AND METHODS
 TOP
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS
 DISCUSSION
 APPENDIX
 ACKNOWLEDGEMENTS
 REFERENCES
 
Two MCDM techniques were used, goal programming (GP) and compromise programming (CoP), to examine whether a single optimal solution or a set of nearly optimal solutions could better achieve the objectives proposed. Goal programming deals with different goals that are included in the GP model as approximate equalities; CoP seeks to determine a good compromise between different objectives rather than optimizing one of them.

Technical Coefficients
Technical coefficients correspond to a typical small-scale dairy farm located in Michoacan State, central Mexico. These enterprises comprise two basic components, one livestock production and the other crop production, with both being strongly interrelated (Val-Arreola et al., 2004a). Regarding the livestock component, herd size averaged 13 mature cows with a live weight of 550 kg and a milk yield of 14 L/d over a lactation period of 305 d (Méndez-y-Cazarín et al., 2000). Using this basic information, DM intake, ME, and MP requirements were calculated according to AFRC (1993). Requirements were calculated to meet maintenance, milk production, weight change, and pregnancy needs in each month of pregnancy and lactation. Val-Arreola et al. (2004c) reported that the lactation equation of Dijkstra et al. (1997) was the most appropriate for small-scale farms in Mexico and was adopted to determine milk production requirements. Maintenance requirements were calculated using live weight changes. Regarding the crop component, average land surface was assumed to be 5.9 ha, and crops were grown in 2 cropping cycles, spring-summer and autumn-winter. Table 1Go shows the cultivated crops and their technical coefficients considered in this study, which were taken from Val-Arreola et al. (2004b).


View this table:
[in this window]
[in a new window]
 
Table 1. Technical coefficients for crops grown on small-scale dairy enterprises in central Mexico
 
Constraints
The process of constructing the MCDM models began by defining a set of constraints using a standard LP model framework. The year was divided into 4 periods or quarters, each 91 d. Each cropping cycle covered 2 quarters.

Nutrient requirements were established for the cow in each quarter by summing monthly values. Stage of lactation was established by considering 4 periods: calved plus dry-off (comprising a dry-off period of 59 d plus the first 32 d of lactation), first one-third, second one-third, and last one-third (each one of 91 d).

Constraints are shown in Appendix 1Go. Both inequalities (1) and (2) are similar in that the diet must cover maintenance, pregnancy, lactation, and weight change requirements for ME and MP. The DM intake constraint (3) ensures that DM intake requirement is met by the feeds selected by the model.

Results of studies carried out on farm have shown that N excreted by these small-scale systems is more than 3 times the Mexican standard due to confined space for housing cows and deficiencies in manure management (Méndez-y-Cazarín et al., 2000). To incorporate this element into the model, the equation suggested by Kebreab et al. (2001) was used to estimate N excretion from N intake. Appendix 1Go shows the constraint formulated for N excretion. The left hand side of the inequality (4) was taken from Kebreab et al. (2001) and the constraint based on Méndez-y-Cazarín et al. (2000), who calculated that a typical small-scale dairy herd with cows averaging 585 kg in liveweight produces about 428 kg of N per quarter. Constraint (5) ensures that the amount of homegrown feed required does not exceed the capacity of crops grown on the land available.

Constraints related to crop production were formulated to ensure coherent land distribution from the solution provided by the model; that is, to avoid conflicts between different cropping areas suggested by the solution. Therefore, restrictions represented by equations (6) and (7) have the purpose of allocating the same land surface to a given crop in both quarters of the cropping cycle. Constraint (8) was formulated to distribute the area for perennial and annual crops realistically. If the solution selects alfalfa (hay or fresh cut), the area used for this crop need to be the same in all quarters, which is the purpose of constraint (9). The crop cost component in the model is regulated by (10).

Constraint (11) ensures the number of cows calving in each quarter does not exceed the herd size. Restriction (12) and (13) were formulated to ensure the number of cows at a given lactation stage in a given calendar quarter was not greater than in the previous stage in the previous quarter. Constraint (14) aims to maintain at least a minimum production of milk in each quarter.

Pay-Off Matrix and Initial Solutions
Definition of the pay-off matrix (Table 2Go) commenced by selecting 4 objectives to be optimized: i) maximization of ME for milk production, ii) maximization of margin of income over feed costs, iii) maximization of MP for milk production, and iv) minimization of the purchase of feedstuffs. These objectives were considered the most appropriate for modeling a representative small-scale dairy farm. Use of GP and CoP involves definition of targets and ideal points, respectively. The approach of Piech and Rehman (1993) was adopted, in which such values were obtained by solving conventional LP problems for each objective.


View this table:
[in this window]
[in a new window]
 
Table 2. Pay-off matrix for the 4 objective functions optimized1
 
Appendix 2Go shows the functions formulated for each of the objectives mentioned above. Objectives (1) and (2) were formulated assuming a sufficiency of nutrients provided by the crops to cover ME and MP requirements for maintenance, pregnancy, and weight change. The remaining ME will be utilized for milk production. Although all the objective functions are shown in Appendix 2Go, only one objective was optimized at a time when the single LP models were used to create the pay-off matrix and determine the best value achieved for each objective.

Goal Programming
Goal programming is one of the most widely used MCDM techniques. Its main characteristic is its capacity to deal with different goals, which are included in the model as approximate equalities. The right hand side values are the targets that may or may not be satisfied. Thus, the purpose of GP is to minimize the deviations between achievement of the goals and their targets. Therefore, GP includes positive (pj) and negative (nj) deviation variables that indicate over-achievement or under-achievement accordingly (Romero and Rehman, 1989). Goal programming generally involves a composite objective function, in which goal satisfaction may be traded-off using relative weights on deviations from the targets. Mathematically, the GP problem can be expressed as follows:


Formula

where wj is the weight of goal j, nj is the negative deviation variable, pj is the positive deviation variable, F is the feasible set, x is the vector of decision variables, gj is the objective target for the goal j, and Zj(x) is the objective function for the goal j as defined in Appendix 2Go. If all the goals have the same importance, then all the weights are the same; that is, w1 = w2 = ... = wn = 1. The general aim in GP is to minimize the deviations between achievement of the goals and their targets; thus, variable nj is included in the optimal basis if the gj target is underachieved and variable pj if the gj target is overachieved. Thus, the objective function in the model minimizes the sum of appropriate deviation variables for each goal.

Compromise Programming
Because it is most unlikely that all the objectives can be optimized simultaneously, an alternative is to seek a good compromise amongst the different objectives rather than optimizing only one of them. Compromise programming seeks to satisfy simultaneously a number of objectives subject to a given set of constraints. The objective is to find a subset of efficient solutions known as pareto-optimal or nondominated solutions from a set of feasible solutions (Romero and Rehman, 1989).

First, ideal points need to be defined using the optimal values underlined in the pay-off matrix (Table 2Go), to which the solution obtained by the CoP model must be as close as possible. To measure closeness, a distance function is introduced into the analysis, which minimizes the distance between the solution and the ideal points. Depending on the measure of distance used, a set of compromise solutions can be obtained (Romero et al., 1987). Two metrics were considered: L1 (≥0) representing the longest distance geometrically, which was minimized by solving the following LP problem (Cohon, 1978):


Formula

where F is the feasible set and x is the vector of decision variables, zj* and z*j are the ideal and antiideal values for objective j (the ideal value was assumed to be the best value from the pay-off matrix and the antiideal the worst), zj(x) is the jth objective function and wj (>0) is the weight attached to objective j.

The second metric used to obtain a compromise solution is L{infty}, the largest deviation from among the individual deviations. This is minimized, and the LP problem to be solved is the following (Cohon, 1978):


Formula

where d{infty} is the largest deviation. It was assumed in this study that w1 = w2 = ... = wn = 1 for both metrics.


    RESULTS
 TOP
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS
 DISCUSSION
 APPENDIX
 ACKNOWLEDGEMENTS
 REFERENCES
 
Ordinary Linear Programming
Table 2Go shows the pay-off matrix in which the optimal solutions for each objective function specified in the LP model are underlined, together with the corresponding values of the other objectives for each optimal. For example, the elements of the first row of Table 2Go show that the maximum ME available for milk production is 1,628,637 MJ, which corresponds to a margin of M$81,553, MP for milk production of 12,809 kg, and a cost of purchased feedstuffs of M$39,127. Table 2Go illustrates the conflict between maximizing ME and MP and maximizing margin. Also, there is conflict between maximizing ME and MP and minimizing the cost of purchased feeds.

Table 3Go shows feed cost, milk yield, feed cost per liter of milk, DM supply, and DM required when each objective is optimized individually. The lowest costs were obtained by maximizing margin over feed costs and minimizing purchase of external feeds, and the most expensive ration was obtained when ME and MP were maximized. In all LP models, the most expensive quarter was the second. The models that gave a much higher DM supply than DM required by the herd were those that maximized ME and MP. The table shows the conflict between maximizing nutrient production and maximizing farm profitability.


View this table:
[in this window]
[in a new window]
 
Table 3. Feed cost, milk yield, DM supply, and DM required for each single objective linear programming model
 
Optimum land distribution for each objective and forage production strategies suggested by the models are shown in Figure 1Go. In the case of the model that maximizes ME, the 3 crops alfalfa hay, fresh ryegrass, and corn silage were chosen. With the model that maximizes margin, the same crops were selected but less land was allocated. In the case of maximized MP, only fresh alfalfa was selected, and similarly, only one crop was selected when minimum purchased feeds was the objective. The models that maximized ME and MP provided a higher DM supply because they allocated a higher proportion of land to crops than the other 2 models. However, all the models indicate that it is possible to maintain current milk production utilizing only forage in the ration.


Figure 1
View larger version (27K):
[in this window]
[in a new window]
 
Figure 1. Land allocation from solutions of each single objective linear programming model: a) maximize ME for milk production, b) maximize margin over feed cost, c) maximize MP for milk production, and d) minimize purchased feeds.

 
Goal Programming
Table 4Go shows the deviational values of each goal specified in the GP model. The overachievement values obtained were margin over feed costs and cost of purchased feeds. The overachievement of purchased feed cost was less desired. However, overachievement of margin over feed cost represents 52% more than the desired goal. Considering the deviational values, the GP model showed a better performance than the single objective LP models. In terms of land distribution (Table 5Go), the GP model selects the crops fresh alfalfa, fresh ryegrass, and corn silage, like the LP model that maximizes ME for milk production but with fresh alfalfa replacing alfalfa hay (Figure 1Go). Obviously, this solution supplied higher quantities of DM and therefore higher quantities of ME, MP, and milk (Table 5Go). Excess of available DM may, in theory, be sold or offers the possibility of increasing the herd size of these farms.


View this table:
[in this window]
[in a new window]
 
Table 4. Optimal deviation variables from solving the goal-programming model
 

View this table:
[in this window]
[in a new window]
 
Table 5. Land allocation, feed cost, milk yield, DM supply, and DM required from the goal-programming model solution
 
Herd composition and milk production suggested by the GP model are shown in Table 6Go. The lowest milk production level was achieved during the second quarter and the largest in the fourth quarter. In comparing herd structure with milk yield in each quarter, it can be seen that the lowest milk production level was when most of the cows were in the calved or dry-off stage. To match forage availability with ME needs by manipulating the herd structure, the GP model suggests most of the cows should calve in the second and third calendar quarters.


View this table:
[in this window]
[in a new window]
 
Table 6. Herd structure and milk production suggested by the goal and compromise programming models1
 
Compromise Programming
Compromise solutions obtained from the CoP model are shown in Table 7Go. Using the L1 and L{infty} metrics, nearest values to the ideal point for each of the objectives specified in the pay-off matrix can be determined. In general there was underachievement with regard to the ideal points established in the pay-off matrix (Table 2Go). The only objective that was met by CoP was minimizing purchased feeds. The CoP model shows (similar to the GP model) that most of the cows should calve in the second and third quarters. However, the herd structure suggested by the CoP model gives higher milk production than the GP model (Table 6Go).


View this table:
[in this window]
[in a new window]
 
Table 7. Land allocation and most efficient compromise set for small-scale dairy farms according to the compromise programming model1
 
Table 7Go also shows the best set of compromise solutions for L1 and L{infty}. Like the GP model, the CoP model suggested a forage production strategy based on alfalfa, fresh ryegrass, and corn silage, and a similar pattern of land allocation. The land allocated differs from some of the LP solutions. The CoP model assigned the same land as the LP model that maximizes ME and more land than the one that maximizes margin. The feed costs, milk yield, DM supply, and DM required given by the compromise solutions are shown in Table 8Go. The milk cost generally showed more uniform values between quarters than those values obtained by the single objective LP models.


View this table:
[in this window]
[in a new window]
 
Table 8. Feed cost, milk yield, DM supply, and DM required according to the compromise programming model1
 

    DISCUSSION
 TOP
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS
 DISCUSSION
 APPENDIX
 ACKNOWLEDGEMENTS
 REFERENCES
 
Based on the results shown in the previous section, 2 main issues require examination. First, whether MCDM is a better approach than ordinary LP for modeling small-scale dairy enterprises, and second, whether the models provide a solid basis for decision-making aimed at improving the performance of such enterprises.

The pros and cons of GP and CoP have been pointed out in different works, with computational difficulties and lack of relevant information provided by the solutions noted as the main drawbacks experienced (Romero and Rehman, 1989; Piech and Rehman, 1993). When GP is used, there is a possibility of obtaining equivalent solutions to those from an LP version of the same problem. Barnett et al. (1982) applied GP to model a subsistence farm in Senegal. They also tried to develop a method of determining goal weights for the model. Their model did not improve on the performance of a profit-maximization LP model. The choice of a single objective function, however, is often equivocal, and in our case the GP model provided a coherent solution that differed from the single objective LP model solutions, which were quite diverse. Although all goals were assumed to have the same importance in this study, the GP model has the facility to incorporate weights that reflect preferences of the farmer or decision maker, giving an advantage over conventional LP.

The CoP approach represents a tool that offers flexibility in evaluating small-scale dairy enterprises through the L1 and L{infty} metrics that enable the decision maker to select the best compromise. Piech and Rehman (1993) considered CoP to be the best option amongst MCDM techniques, in that L1 and L{infty} can be solved using 2 standard LP models. However, in the literature reviewed, applications of CoP were scarce, making it difficult to identify possible limitations in technique applicability. An interesting application of the technique was the study by Zekri and Romero (1992) to assess irrigation agriculture in a village in Spain to measure the opportunity cost of maintaining the current situation from the point of view of a public and private perspective. Their approach comprised 2 phases. First, a GP model was used to approximate the current situation, and second, a CoP model was used to assess the current situation in terms of criteria relevant to the perspectives mentioned above. In our case, the CoP model was not significantly better than the GP one, but applying both improved upon using a single objective LP model.

Regarding the second issue of this discussion, the system can be improved by effecting 2 main changes. The first is related to selection of the most appropriate forage production strategy by considering land availability and cropping cycles. The second is by rescheduling the calving season to synchronize higher demand for nutrients with periods when forage availability is highest. According to the solutions provided by the GP and CoP models, an increase in nutrient (ME) availability by utilizing a larger land surface has a positive economic effect.

The forage production strategy selected by the models (basically comprising alfalfa, ryegrass, and corn silage) is similar to one (of several) identified in a previous work (Val-Arreola et al., 2004b). Dry matter supply by the forage production strategy suggested by our models was also in agreement with this identified strategy. The area proportion of forage to corn found in that work was in agreement with that suggested by the MCDM models (44:56; Table 7Go). These results were also consistent with the proportion given by Castelán-Ortega et al. (2003b).

In the work described herein, all the models selected similar ingredients to feed the animal, with most showing differences in the quantities recommended. Such quantities were dependent on the different objectives established and the type of model utilized to meet those objectives. For most of the cases, the feeding strategy was based on alfalfa, fresh ryegrass, and corn silage.

None of the models advocated supplementation with corn grain or commercial concentrate. Castelán-Ortega et al. (2003b) developed a simulation model for small-scale dairy farms in the highlands of Mexico that used corn straw, corn grain, and improved pasture as a feeding strategy. Previous work on small-scale dairy farms in Michoacan State have described feeding strategies based on fresh alfalfa, corn grain plus straw or commercial concentrate, according to availability and quality of forage, and in some cases a second forage such as oat hay (Val-Arreola, 1998; Méndez-y-Cazarín et al., 2000). Because our models considered factors like seasonality of forage availability, herd structure, herd requirements, and milk production, the solutions agree with studies that have pointed out that use of commercial concentrate in Mexico is determined by quality and availability of forage (Arriaga-Jordan et al., 2002; Castelán-Ortega et al., 2003b). To avoid such shortages, strategies that allow forage production over the whole year are needed. Arriaga-Jordan et al. (2002) suggest incorporating corn silage with the baseline forage. Schiere et al. (1999) suggest a way to adapt a farm is by matching quality and availability of resources with different types of animals to explore different adjustments between animals and feeds to achieve maximum system output. The way our models realized such adjustment was by allocating animals according to their lactation stage in each period of the year to match their nutritional requirements with available nutrients. Castelán-Ortega et al. (2003b) found using their model that the number of lactating cows was higher during the period with higher availability of feed of better quality (half of year). Our GP and CoP models allocated a higher number of lactating cows to the first and last quarters (Table 6Go).

Results obtained from the models provide information about how the forage strategy selected and the land utilized affect economic performance. Higher margins were obtained when most of the surface available was cultivated. Our models also provide information as to how forage shortage may be alleviated by modifying calving period to match nutrient availability with requirements at each lactation stage, with a minimum reduction in milk production. Excess of DM supply in relation to DM requirement in the case of the CoP and GP models agrees with previous work that suggested adopting a forage strategy based on alfalfa, ryegrass, and corn silage ensures that the nutrient requirements of the herd will be met (Val-Arreola et al., 2004b).


    APPENDIX
 TOP
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS
 DISCUSSION
 APPENDIX
 ACKNOWLEDGEMENTS
 REFERENCES
 

Appendix 1. Summary of constraints
Constraint Description of variables

1. Formula LEi = ME requirement for milk production in stage i of lactation (MJ per cow),
xij = number of cows at stage i in quarter j,
NLEi = ME requirement for maintenance, pregnancy and weight change in stage i of lactation (MJ per cow),
MEk = ME content of feed k (MJ/kg of DM),
yjk = quantity of feed k required in calendar quarter j; j=1,2,3,4 (kg of DM),
2. Formula LPi = MP requirement for milk production in stage i of lactation (g per cow),
NLPi = MP requirement for maintenance, pregnancy and weight change in stage i of lactation (g per cow),
MPk = MP supplied by feed k (g/kg of DM).
3. Formula DMIi = estimated DM intake for cows in stage i of lactation (kg per cow).
4. Formula Formula, quantity of nitrogen potentially excreted by consuming feed k (g/kg DM),
CP = crude protein content of feed k (g/kg of DM),
Nexcretion = maximum quantity of nitrogen that can be excreted by the cows (428 kg per quarter).
5. Formula DMPjk = DM production on farm of feed k in quarter j (kg/ha),
LND = total land available, 5.9 ha.
k summed over homegrown crops.
6. y2ky1k = 0 y2k = quantity required of homegrown crop k for second quarter (kg of DM),
y1k = quantity required of homegrown crop k for first quarter (kg of DM).
7. y4ky3k = 0 y4k = quantity required of homegrown crop k for fourth quarter (kg of DM),
y3k = quantity required of homegrown crop k for third quarter (kg of DM).
8. Formula k summed over homegrown perennial crops,
l summed over homegrown annual crops.
9. Formula Subscript 1 denotes fresh alfalfa,
Subscript 2 denotes alfalfa hay.
10. Formula cjk = cost of feed k in quarter j (M$/kg of DM),
UCap = upper limit on capital available in a given quarter (M$22,420).
11. Formula X = number of mature cows in herd (13 cows).
12. Formula I = 1,2,3; j = 1,2,3.
13. Formula j = 1,2,3.
14. Formula MYi = milk yield of cow in stage i of lactation (L),
LMilk = lower limit on milk production in a given quarter (L).


Appendix 2. Summary of objective functions
Objective Function Description of variables

1. Maximize ME for milk production. Formula
2. Maximize margin over feed cost Formula pj = milk price in quarter j (M$/kg)
3. Maximize MP for milk production Formula
4. Minimize purchased feeds Formula k summed over purchased feeds


    ACKNOWLEDGEMENTS
 TOP
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS
 DISCUSSION
 APPENDIX
 ACKNOWLEDGEMENTS
 REFERENCES
 
Financial support for D. Val-Arreola was provided by the National Council of Science and Technology (CONACyT) of Mexico.

Received for publication June 29, 2005. Accepted for publication November 14, 2005.


    REFERENCES
 TOP
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS
 DISCUSSION
 APPENDIX
 ACKNOWLEDGEMENTS
 REFERENCES
 


Agricultural and Food Research Council (AFRC). 1993. Energy and protein requirements of ruminants. CAB International, Wallingford, UK.

Arriaga-Jordan, C. M., B. Albarran-Portillo, A. Espinosa-Ortega, A. Garcia-Martinez, and O. A. Castelan-Ortega. 2002. On-farm comparison of feeding strategies based on forages for small-scale dairy production systems in the highlands of central Mexico. Exp. Agric. 38:375–388.

Arriaga-Jordan, C. M., and R. A. Pearson. 2002. The contribution of livestock to small-holder livelihoods. Pages 21–22 in Responding to Increasing Global Demand for Animal Products. Proc. Br. Soc. Anim. Sci., Am. Soc. Anim. Sci., Mexican Soc. Anim. Prod., Univ. Yucatan, Merida, Mexico. Br. Soc. Anim. Sci., Edinburgh, UK.

Barnett, D., B. Blake, and B. A. McCarl. 1982. Goal programming via multidimensional scaling applied to Senegalese subsistence farms. Am. J. Agric. Econ. 64:720–727.

Candler, W. V. 1983. A demonstration model for the dairy industry of La Laguna. Pages 511–552 in The Book of CHAC: Programming Studies for Mexican Agriculture. R. D. Norton and L. Solis, ed. The John Hopkins University Press, Baltimore, MD.

Castelán-Ortega, O. A., R. H. Fawcett, C. Arriaga-Jordán, and M. Herrero. 2003a. A decision support system for smallholder campesino corn-cattle production systems of the Toluca Valley in Central Mexico. Part I—Integrating biological and socio-economic models into a holistic system. Agric. Syst. 75:1–21.

Castelán-Ortega, O. A., R. H. Fawcett, C. Arriaga-Jordán, and M. Herrero. 2003b. A decision support system for smallholder campesino corn-cattle production systems of the Toluca Valley in Central Mexico. Part II—Emulating the farming system. Agric. Syst. 75:23–46.

Cohon, J. 1978. Multiobjective Programming and Planning. Academic Press, New York, NY.

Dijkstra, J., J. France, M. S. Dhanoa, J. A. Maas, M. D. Hanigan, A. S. Rook, and D. E. Beever. 1997. A model to describe growth patterns of the mammary gland during pregnancy and lactation. J. Dairy Sci. 80:2340–2354.[Abstract]

Espinoza-Ortega, A., C. M. Arriaga-Jordan, and O. A. Castelan-Ortega. 1997. Economic analysis of the peasant production system at the Toluca Valley. Pages 94–102 in Proc. Small-Scale Dairy Syst., Toluca, Mexico. Autonomous University of Mexico State.

Gallardo-Nieto, J. L., L. Villamar-Angulo, H. Perez-Frias, and E. Olivera-Cazares. 2004. Current situation and perspectives of the Mexican dairy industry 2004. Ministry of Agriculture Bulletin (SAGARPA). http://www.sagarpa.gob.mx/Dgg/estudio/sitlech04.pdf Accessed June 27, 2005.

Kebreab, E., J. France, D. E. Beever, and A. R. Castillo. 2001. Nitrogen pollution by dairy cows and its mitigation by dietary manipulation. Nutr. Cycl. Agroecosyst. 60:275–285.

Méndez-y-Cazarín, M. D., R. Tzintzun-Rascón, and D. Val-Arreola. 2000. Production evaluation, environmental effects and problems in small-scale dairy farms. Livest. Res. Rural Devel. 12(1). Online. Available: http://www.cipav.org.co/lrrd/lrrd12/1/manu121.htm Accessed Apr. 5, 2005. [Spanish with English abstract]

Piech, B., and T. Rehman. 1993. Application of multiple criteria decision-making methods to farm planning: A case study. Agric. Syst. 41:305–319.

Rehman, T., and C. Romero. 1993. The application of the MCDM paradigm to the management of the agricultural systems: some basic considerations. Agric. Syst. 41:239–255.

Romero, C., F. Amador, and A. Barco. 1987. Multiple objectives in agricultural planning: A compromise programming application. Am. J. Agric. Econ. 69:78–86.

Romero, C., and T. Rehman. 1989. Multiple Criteria Analysis for Agricultural Decisions. Elsevier, Amsterdam, The Netherlands.

Schiere, J. B., J. D. Wit, F. A. Steenstra, and H. V. Keulen. 1999. Design of farming systems for low input conditions: principles and implications based on scenarios studies with feed allocation in livestock production systems. Neth. J. Agric. Sci. 47:169–183.

Val-Arreola, D. 1998. Maximisation of the margin of incomes over feed cost in small-scale dairy farms. M.Sc. Diss., Univ. Michoacan, Morelia, Mexico.

Val-Arreola, D., E. Kebreab, J. A. N. Mills, and J. France. 2005. Analysis of feeding strategies for small-scale dairy systems in central Mexico using linear programming. J. Anim. Feed Sci. 14:607–624.

Val-Arreola, D., E. Kebreab, J. A. N. Mills, and J. France. 2004a. Optimising the calving pattern of herds in small-scale dairy systems in central Mexico. J. Anim. Feed Sci. 13:557–573.

Val-Arreola, D., E. Kebreab, J. A. N. Mills, S. Wiggins, and J. France. 2004b. Forage production and nutrient availability in small-scale dairy systems in central Mexico using linear programming and partial budgeting. Nutr. Cycl. Agroecosyst. 69:191–201.

Val-Arreola, D., E. Kebreab, J. Dijkstra, and J. France. 2004c. Study of the lactation curve in dairy cattle on farms in central Mexico. J. Dairy Sci. 87:3789–3799.[Abstract/Free Full Text]

Wiggins, S. 2002. Small-scale dairying as the key piece in the rural livelihoods. Pages 197–210 in Proc. New Trends in the Socioeconomic Analysis of the Dairy Industry in the Globalisation Context, Toluca, Mexico. Autonomous University of Mexico State.

Zekri, S., and C. Romero. 1992. A methodology to assess the current situation in irrigated agriculture: An application to the village of Tauste (Spain). Oxf. Agrar. Stud. 20:75–88.


This article has been cited by other articles:


Home page
J ANIM SCIHome page
E. Kebreab, J. Dijkstra, A. Bannink, and J. France
Recent advances in modeling nutrient utilization in ruminants
J Anim Sci, April 1, 2009; 87(14_suppl): E111 - E122.
[Abstract] [Full Text] [PDF]


This Article
Right arrow Abstract Freely available
Right arrow Full Text (PDF)
Right arrow Interpretive Summary
Right arrow Alert me when this article is cited
Right arrow Alert me if a correction is posted
Services
Right arrow Similar articles in this journal
Right arrow Similar articles in PubMed
Right arrow Alert me to new issues of the journal
Right arrow Download to citation manager
Right arrow reprints & permissions
Citing Articles
Right arrow Citing Articles via HighWire
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Val-Arreola, D.
Right arrow Articles by France, J.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by Val-Arreola, D.
Right arrow Articles by France, J.


HOME HELP FEEDBACK SUBSCRIPTIONS ARCHIVE SEARCH TABLE OF CONTENTS