Differential Assessment of Strategies to Increase Milk Yield in Small-Scale Dairy Farming Systems Using Multi-Agent Modelling and Simulation
Abstract
:1. Introduction
Agent-Based Modelling of Farming Systems
2. Materials and Methods
2.1. Average Milk Yield Model
- a.
- Design concepts
- Agent cow: The agent required an amount of energy daily to achieve its potential production. The only source of feeds for the agent was assumed to be grass. Amount of grass eaten by the agent is converted to energy. Daily energy requirements and rates of conversion of grass into energy were obtained from literature [33,34]. The agent’s energy use was allocated for movement, maintenance, resting, milk production, and gain in weight [35]. The model design considered the effects of infectious diseases that can affect milk yield. In order to minimize complexity, a random distribution of clinical mastitis was adopted [36]. By calculating the total probability of catching the infection, the agent cow eating and drinking rates were determined;
- Agent observer: This agent controls changes in the layers, which in turn affects the behaviour of the cow and farmer agents. The observer agent controls grass growth rate and the distribution of mastitis base probability in rainy season and reduces in dry season;
- Agent farmer: This is the farm manager, and is responsible to fill feeds and water troughs for the cows. This agent inputs the dataset that defines how the farm is managed. A change in the input data for the farmer has an impact on the farm productivity, including feeding pattern and annual milk yield from the cow agents. It is the farmer’s set of choices and strategies that is altered to observe changes in the milk yield. The alteration is conducted by adjusting values of the evolvement strategies given in Table 1;
- Adaptation: The cow agent decision room had to decide when to eat and drink based on energy, hunger, and thirsty values. The cow could move to the feeds and water troughs every time either its energy value was below the minimum requirement, hungry or thirsty, observing its food value, current energy level, and hydration in respect to amount of food eaten. Since these values were stored in the cow’s memory, their alteration alerted an action by the cow.
- Objectives: The farmers’ main objective was to increase the amount of milk, which is measured in kilograms per year. The changes in daily milk yield were affected by the choices of management strategies indicated in Table 1; where different values were tested. The cows’ main objective was to increase total energy value measured in MJ to achieve daily requirements and milk yield;
- Interactions: Direct interactions were between individual cow and grass biomass layer, individual cow and feeds and water trough layer, observer and grass biomass layer, farmer and feed and water trough layer, and observer and mastitis layer. Indirect interactions were between individual cows competing for grass and water as they strived to reach their daily objective;
- Stochasticity: Increase in grass biomass, feed and water trough refilling, and distribution of probability for infectious mastitis in the environment;
- Observations: At each simulation time, observations were made on the cow agent (milk level, body weight, and probability of infection), grass layer (increase in biomass), feeds and water trough layer (increase in feeds and water level), and the mastitis layer (random growth and decline of infection probabilities).
- b.
- Initialisation
- Cow agent: for each cow agent, the following were the initial values as gathered from existing research [33,34]:
- ○
- Location: randomly distributed in 45 × 45 grid cells;
- ○
- Instances: 5;
- ○
- Energy (MJ) = 0;
- ○
- Milk level (kg) = 0.5;
- ○
- Probability of infection = ((0.16 × Mhistory)/5) + ((0.16 × parity)/5) + (0.16 × (1/breed));
- ○
- Body weight (kg) = initial value given for each breed in an input data file;
- ○
- Breed = fixed breed type given for each agent in an input data file;
- ○
- Mastitis history = initial value given for each agent in an input data file;
- ○
- Parity = initial value given for each agent in an input data file;
- ○
- Farmer agent: the farmer agent was randomly distributed in a 45 × 45 grid cells, and a set of values for each evolvement determinant (as indicated in Table 1) was given as an input file;
- ○
- Observer agent: for the observer agent, initial values were a random distribution in a 45 × 45 grid cells, and one instance.
- c.
- Model layersThe implementation of layers was used to have highly flexible model interactions and improve focus on the individual agent [31]. Figure 1 shows the organisation of three independent layers, in which the agents interact with each other.
- Grass layer: Each cell in a 100 × 100 grid contained 4 kg of grass at initialisation and increases/reduces randomly during simulation;
- Feeds and water trough: 4 × 9 cells for feeds and 9 × 9 cells for water out of 100 × 100 grid. The farmer agent randomly increased values of the cells from 0 (at initialisation) depending on management practices given in the farmer input file;
- Mastitis base probability: Each cell in a 100 × 100 grid had 0 base probability at initialisation and increased or reduced randomly during simulation.
- d.
- Input dataInput data were generated for the model layers, cow, and farmer agents. The raster files for the layers were in .asc format and each cell carried a value of 4 for the grass raster, and 0 for the feeds and water trough, and the mastitis raster files.
- Grass raster file (100 × 100);
- Feeds and water trough raster file (100 × 100);
- Mastitis probability raster file (100 × 100);
- Cow agent input csv file;
- Farmer agent input csv file.
2.2. Milk Yield
- Energy requirements (equation adapted from [37])
- Derived probability of infection
- Derived energy use
- Derived excess energy
- Derived weight gain
2.3. Simulation Scenarios
- To understand the effect of the evolvement strategies (Table 1) on average annual milk yield;
- To identify among the categories of evolvement strategies (farm, farmer, income, and infrastructure characteristics), which category has a higher likelihood of improving milk yield.
3. Results
3.1. Assessment of Farm Characteristics
3.2. Assessment of Farmer Characteristics
3.3. Assessment of Income Characteristics
3.4. Assessment of Infrastructure Characteristics
4. Discussion
5. Conclusions
Author Contributions
Funding
Institutional Review Board Statement
Data Availability Statement
Acknowledgments
Conflicts of Interest
References
- Chawala, A.; Banos, G.; Peters, A.; Chagunda, M. Farmer-preferred traits in smallholder dairy farming systems in Tanzania. Trop. Anim. Health Prod. 2019, 51, 1337–1344. [Google Scholar] [CrossRef] [PubMed] [Green Version]
- Mavura, F.; Pandhare, S.M.; Mkoba, E.; Nyambo, D.G. Rule-based engine for automatic allocation of smallholder dairy producers in preidentified production clusters. Sci. World J. 2022, 2022, 6944151. [Google Scholar] [CrossRef] [PubMed]
- Cook, R. World Cattle Inventory: Ranking of Countries (FAO). 2015. Available online: https://www.drovers.com/article/world-cattle-inventory-ranking-countries-fao.Otte (accessed on 30 May 2019).
- Smith, J.; Sones, K.; Grace, D.; MacMillan, S.; Tarawali, S.; Herrero, M. Beyond milk, meat, and eggs: Role of livestock in food and nutrition security. Anim. Front. 2013, 3, 6–13. [Google Scholar] [CrossRef] [Green Version]
- Tanzania. Tanzania Livestock Modernization Initiative; Ministry of Livestock and Fisheries Development: Dar es Salaam, Tanzania, 2016. [Google Scholar]
- Hemme, T.; Otte, J. Status and Prospects for Smallholder Milk Production: A Global Perspective; FAO: Rome, Italy, 2010. [Google Scholar]
- Nyambo, D. Use of Agent-Based Models in Characterizing Farm Types and Evolvement in Smallholder Dairy Systems [NM-AIST]. 2020. Available online: https://dspace.nm-aist.ac.tz/handle/20.500.12479/895 (accessed on 20 February 2023).
- Kingu, D.; Ndiege, B.O. Empowering small scale dairy farmers through the cooperatives model. J. Coop. Bus. Stud. 2018, 2, 1–12. [Google Scholar]
- Gitau, K.J. Factors Influencing Milk Production Among Small Scale Dairy Farmers in Mirangine in Nyandarua County and Mauche in Nakuru County. Ph.D. Thesis, University of Nairobi, Nairobi, Kenya, 2013. [Google Scholar]
- Nyambo, D.G.; Luhanga, E.T.; Yonah, Z.O. Characteristics of smallholder dairy farms by association rules mining based on apriori algorithm. Int. J. Soc. Syst. Sci. 2019, 11, 99–118. [Google Scholar] [CrossRef]
- Taneja, M.; Jalodia, N.; Malone, P.; Byabazaire, J.; Davy, A.; Olariu, C. Connected cows: Utilizing fog and cloud analytics toward data-driven decisions for smart dairy farming. IEEE Internet Things Mag. 2019, 2, 32–37. [Google Scholar] [CrossRef] [Green Version]
- Gondro, C.; Okeyo, A.M.; Ojango, J.; Rao, J.; Dessie, T.; Rege, J.E.O.; Gibson, J.P. Genetic diversity of the indigenous cattle of Kenya, Uganda, Ethiopia and Tanzania using high-density SNP data results and discussion. In Proceedings of the World Congress on Genetics Applied to Livestock Production, Volume Electronic Poster Session—Genetic Gain—Strategies for Local Breeds, Auckland, New Zealand, 11–16 February 2018; Volume 2, p. 777. [Google Scholar]
- Nyambo, D.; Luhanga, E.; Yonah, Z.; Mujibi, F.N. Application of multiple unsupervised models to validate clusters robustness in characterizing smallholder dairy farmers. Sci. World J. 2019, 12, 1020521. [Google Scholar] [CrossRef]
- Maleko, D.; Msalya, G.; Mwilawa, A.; Pasape, L.; Mtei, K. Smallholder dairy cattle feeding technologies and practices in Tanzania: Failures, successes, challenges and prospects for sustainability. Int. J. Agric. Sustain. 2018, 16, 201–213. [Google Scholar] [CrossRef]
- Mujibi, F.D.; Rao, J.; Agaba, M.; Nyambo, D.; Cheruiyot, E.K.; Kihara, A.; Zhang, Y.; Mrode, R. Performance evaluation of highly admixed Tanzanian smallholder dairy cattle using SNP derived kinship matrix. Front. Genet. 2019, 10, 375. [Google Scholar] [CrossRef]
- Urassa, J.; Raphael, E. The Contribution of Small Scale Dairy Farming to Community Welfare: A Case Study of Morogoro Municipality; Sokoine University of Agriculture: Morogoro, Tanzania, 2017. [Google Scholar]
- Manzo, G.; Matthews, T. Potentialities and limitations of agent-based simulations. Rev. Fr. Sociol. 2014, 55, 653–688. [Google Scholar] [CrossRef]
- Mack, G.; Huber, R. On-farm compliance costs and N surplus reduction of mixed dairy farms under grassland-based feeding systems. Agric. Syst. 2017, 154, 34–44. [Google Scholar] [CrossRef]
- Clemen, T.; Ahmady-Moghaddam, N.; Glake, D.; Lenfers, U.A.; Ocker, F.; Osterholz, D.; Ströbele, J. Toward a movement paradigm for artificial human agents. In Proceedings of the 2022 Annual Modeling and Simulation Conference (ANNSIM), San Diego, CA, USA, 18–20 July 2022; pp. 176–187. [Google Scholar]
- Seid, M.; Bridgeland, D.; Bridgeland, A.; Hartley, D.M. A collaborative learning health system agent-based model: Computational and face validity. Learn. Health Syst. 2021, 5, e10261. [Google Scholar] [CrossRef] [PubMed]
- Li, L.M.W.; Wang, S.; Lin, Y. The casual effect of relational mobility on integration of social networks: An agent-based modeling approach. Curr. Psychol. 2022, 1–17. [Google Scholar] [CrossRef] [PubMed]
- Macmillan, W.; Huang, H.Q. An agent-based simulation model of a primitive agricultural society. Geoforum 2008, 39, 643–658. [Google Scholar] [CrossRef]
- Nolan, J.; Parker, D.; Kooten, G.C.V.; Berger, T. An overview of computational modelling in agricultural and resource economics. Can. J. Agric. Econ. 2009, 57, 417–429. [Google Scholar] [CrossRef]
- Faliszewski, P.; Rothe, J. Control and bribery in voting. Chapter 7. In Handbook of Computational Social Choice; Cambridge University Press: New York, NY, USA, 2016; pp. 146–168. ISBN 978-1-107-06043-2. [Google Scholar]
- Hahn, B.; Richardson, F.; Hoffman, M.; Roberts, R.; Todd, S.; Carrick, P. A simulation model of long-term climate, livestock and vegetation interactions on communal rangelands in the semi-arid Succulent Karoo, Namaqualand, South Africa. Ecol. Model. 2005, 183, 211–230. [Google Scholar] [CrossRef]
- Fust, P.; Schlecht, E. Integrating spatio-temporal variation in resource availability and herbivore movements into rangeland management: RaMDry—An agent-based model on livestock feeding ecology in a dynamic, heterogeneous, semi-arid environment. Ecol. Model. 2018, 369, 13–41. [Google Scholar] [CrossRef]
- Mulindwa, H.; Galukande, E.; Wurzinger, M.; Ojango, J.; Okeyo, A.M.; Sölkner, J. Stochastic simulation model of Ankole pastoral production system: Model development and evaluation. Ecol. Model. 2011, 222, 3692–3700. [Google Scholar] [CrossRef]
- Schilling, C.; Kaye-Blake, W.; Post, E.; Rains, S. The importance of farmer behaviour: An application of Desktop MAS, a multi-agent system model for rural New Zealand communities. In Proceedings of the 2012 Conference New Zealand Agricultural and Resource Economics Society, Nelson, New Zealand, 31 August 2021. [Google Scholar]
- Oudendag, D.; Hoogendoorn, M.; Jongeneel, R. Agent-based modeling of farming behavior: A case study for milk quota abolishment. In Modern Advances in Applied Intelligence, Proceedings of the 27th International Conference on Industrial Engineering and Other Applications of Applied Intelligent Systems, IEA/AIE 2014, Kaohsiung, Taiwan, 3–6 June 2014; Lecture Notes in Computer Science; Springer: Berlin/Heidelberg, Germany, 2014; Volume 8481, pp. 11–20. [Google Scholar] [CrossRef]
- Shiflet, A.B.; Shiflet, G.W. An introduction to agent-based modeling for undergraduates. Procedia Comput. Sci. 2014, 29, 1392–1402. [Google Scholar] [CrossRef] [Green Version]
- Hüning, C.; Adebahr, M.; Thiel-Clemen, T.; Dalski, J.; Lenfers, U.; Grundmann, L. Modeling & simulation as a service with the massive multi-agent system MARS. In Proceedings of the Agent-Directed Simulation Symposium, Pasadena, CA, USA, 3–6 April 2016. [Google Scholar]
- Grimm, V.; Berger, U.; DeAngelis, D.L.; Polhill, J.G.; Giske, J.; Railsback, S.F. The ODD protocol: A review and first update. Ecol. Model. 2010, 221, 2760–2768. [Google Scholar] [CrossRef] [Green Version]
- Sulabo, R.C. Energy Evaluation. Lipids, 2011; 1–12. [Google Scholar]
- AHDB Dairy. The Eternal Question: “How Much Milk Can You Get from Grass?”. 2015. Available online: https://dairy.ahdb.org.uk/news/technical-articles/june-2015/how-much-milk-can-you-get-from-grass/#.W2QcI9gzbOQ (accessed on 3 August 2018).
- Moran, J. Nutrient requirements of dairy cows. In Tropical Dairy Farming: Feeding Management for Small Holder Dairy Farmers in the Humid Tropics; Landlinks Press: Collingwood, Australia, 2015; p. 312. [Google Scholar]
- Østergaard, S.; Sørensen, J.T.; Kristensen, A.R. A stochastic model simulating the feeding-health-production complex in a dairy herd. J. Dairy Sci. 2000, 83, 721–733. [Google Scholar] [CrossRef] [PubMed]
- Bruinenberg, M.; Zom, R.; Valk, H. Energy evaluation of fresh grass in the diets of lactating dairy cows. NJAS—Wagening. J. Life Sci. 2002, 50, 67–81. [Google Scholar] [CrossRef] [Green Version]
- Axegärd, C. Individual Drinking Water Intake of Dairy Cows in an AMS Barn; Swedish University of Agricultural Sciences: Uppsala, Sweden, 2017. [Google Scholar]
- Nielsen, B.L. On the interpretation of feeding behaviour measures and the use of feeding rate as an indicator of social constraint. Appl. Anim. Behav. Sci. 1999, 63, 79–91. [Google Scholar] [CrossRef]
- Pelegrini, D.F.; Lopes, M.A.; Demeu, F.A.; Rocha, Á.G.F.; Bruhn, F.R.P.; Casas, P.S. Effect of socioeconomic factors on the yields of family-operated milk production systems. Semin. Ciênc. Agrár. 2019, 40, 1199–1214. [Google Scholar] [CrossRef] [Green Version]
- Ahmad, D.; Afzal, M.; Abro, A.A. Impact of formal credit on subsistence farmers dairy production in Southern Punjab, Pakistan. Sarhad J. Agric. 2022, 38, 287–294. [Google Scholar] [CrossRef]
- Oliveira, M.; Cocozza, A.; Zucaro, A.; Santagata, R.; Ulgiati, S. Circular economy in the agro-industry: Integrated environmental assessment of dairy products. Renew. Sustain. Energy Rev. 2021, 148, 111314. [Google Scholar] [CrossRef]
- King’ori, S. Factors affecting dairy production among dairy cooperative societies in Kenya. A critical literature review. Am. J. Livest. Policy 2022, 1, 8–19. [Google Scholar]
- Dankar, I.; Hassan, H.; Serhan, M. Knowledge, attitudes, and perceptions of dairy farmers regarding antibiotic use: Lessons from a developing country. J. Dairy Sci. 2022, 105, 1519–1532. [Google Scholar] [CrossRef]
- Britt, J.; Cushman, R.; Dechow, C.; Dobson, H.; Humblot, P.; Hutjens, M.; Jones, G.; Ruegg, P.; Sheldon, I.; Stevenson, J. Invited review: Learning from the future—A vision for dairy farms and cows in 2067. J. Dairy Sci. 2018, 101, 3722–3741. [Google Scholar] [CrossRef] [Green Version]
Farm | Farmer | Income | Infrastructure | |
---|---|---|---|---|
Tanzania | Frequency for vaccination | Years of experience | Retail milk sale | Distance to buyers |
Frequency for watering | Years of schooling | Bulk milk sale | Distance to water sources | |
Total land | Formal training | Total crop sale | Extension visits per year | |
Land for fodder | Membership in groups | Distance to breeding | ||
Stall feeding in rainy | Distance to market | |||
Stall feeding in dry | ||||
Grazing in rainy | ||||
Grazing in dry | ||||
Milking cows | ||||
Herd size | ||||
Employed labour | ||||
Months of purchased fodder | ||||
Breeding method |
S/No | Feature Name | Type | Range |
---|---|---|---|
1 | Grazing in dry season | Boolean | 0 (no) or 1 (yes) |
2 | Grazing in rainy season | Boolean | 0 (no) or 1 (yes) |
3 | Stall feeding in dry season | Boolean | 0 (no) or 1 (yes) |
4 | Stall feeding in rainy season | Boolean | 0 (no) or 1 (yes) |
5 | Frequency for watering/day | Discrete | 0–4 |
6 | Distance to water source (km) | Continuous | 0–15 |
7 | Total land | Continuous | 0–100 |
8 | Land for fodder | Continuous | 0–80 |
9 | Number of employees | Discrete | 1–10 |
10 | Frequency for vaccination/year | Discrete | 0–6 |
11 | Membership in groups | Discrete | 0–5 |
12 | Years of experience | Discrete | 1–50 |
13 | Years of schooling | Discrete | 0–21 |
14 | Breeding method | Boolean | 0 (bull) or 1 (AI) |
15 | Distance to breeding (km) | Continuous | 0–100 |
16 | Frequency of extension visits/year | Discrete | 1–54 |
17 | Herd size | Discrete | 1–50 |
18 | Milking cows | Discrete | 1–20 |
19 | Peak milk production/day | Continuous | 1–40 |
20 | Bulk milk sale (lt) | Continuous | 1–100 |
21 | Retail milk sale (lt) | Continuous | 1–100 |
22 | Distance to buyers (km) | Continuous | 1–37 |
23 | Total crop sale | Continuous | 0–950,000 (Tsh) |
24 | Distance to market (km) | Continuous | 0–8 |
25 | Months of purchased fodder | Discrete | 1–12 |
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content. |
© 2023 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (https://creativecommons.org/licenses/by/4.0/).
Share and Cite
Nyambo, D.G.; Clemen, T. Differential Assessment of Strategies to Increase Milk Yield in Small-Scale Dairy Farming Systems Using Multi-Agent Modelling and Simulation. Agriculture 2023, 13, 590. https://doi.org/10.3390/agriculture13030590
Nyambo DG, Clemen T. Differential Assessment of Strategies to Increase Milk Yield in Small-Scale Dairy Farming Systems Using Multi-Agent Modelling and Simulation. Agriculture. 2023; 13(3):590. https://doi.org/10.3390/agriculture13030590
Chicago/Turabian StyleNyambo, Devotha G., and Thomas Clemen. 2023. "Differential Assessment of Strategies to Increase Milk Yield in Small-Scale Dairy Farming Systems Using Multi-Agent Modelling and Simulation" Agriculture 13, no. 3: 590. https://doi.org/10.3390/agriculture13030590
APA StyleNyambo, D. G., & Clemen, T. (2023). Differential Assessment of Strategies to Increase Milk Yield in Small-Scale Dairy Farming Systems Using Multi-Agent Modelling and Simulation. Agriculture, 13(3), 590. https://doi.org/10.3390/agriculture13030590