Poultry farm distribution models developed along a gradient of intensification
Introduction
Population and income growth, urbanisation and technological advances have led to the intensification of livestock production systems over the past decades. Although this trend in intensification is observable globally, the level of intensification varies between countries and it mainly occurs today in low- and middle-income countries (LMICs). Intensification leads to large-scale and specialized production units, market-oriented and capital- and input-intensive. Such production units make heavy use of feed concentrates, antibiotics and vaccination. The increase in production meets the increasing demand for animal source food, but it has tremendous effects on livelihoods, animal and human health. Concentration of monogastric species have led to pollution of land and water due to the overload of improperly managed manure (Gerber et al., 2005). In terms of public health, intensive livestock production fosters the emergence of infectious diseases (Jones et al., 2013) and often leads to increased use of antibiotics, which is likely to increase antimicrobial resistance (Aarestrup, 2005; Vieira et al., 2011; Chantziaras et al., 2014; Van Boeckel et al., 2015). Livestock also traditionally represents an important source of livelihood for rural populations in LMICs but the intensification of livestock production seldom benefits the local population (FAO, 2018).
As production intensifies, the geographical distribution of monogastric animal farming changes. At a pre-industrial stage, chickens and pigs are mainly associated with human settlements for their role of waste converters (Perry et al., 1999; Steinfeld et al., 2006), with farms homogeneously distributed among rural populations. This situation is encountered in traditional production systems of LMICs (Steinfeld et al., 2006). Production intensification, in parallel to urbanisation, leads larger farms to locate within the peri-urban belt of major consumption centres (Steinfeld et al., 2006). Livestock products being perishable, the production profits from being located close to demand centres. However, as cities expand, land value in the periphery increases and, with improving transport infrastructures, farms move further from peri-urban belts, with the advantage of being closer to feed production or transportation areas (Steinfeld et al., 2006), and benefits from agglomeration economies (Roe et al., 2002a; Herath et al., 2005; Larue et al., 2011). This evolution is particularly pronounced for monogastric species, which are raised in landless production systems (Naylor, 2005; Steinfeld et al., 2006).
As the distribution of farms changes as it intensifies, so do the effects of livestock raising activities on their surroundings. Highly-detailed maps of livestock production allow the assessment of the local effects of intensive production (Robinson et al., 2011). In the case of disease outbreaks for instance, adequate control is planned based on mathematical modelling of infectious disease which require farm distribution and sizes (number of animals) (Hill et al., 2018; Tildesley and Ryan, 2012). Accurate data on the distribution and stock of farms rarely exist as agricultural censuses are time and resource consuming. In high-income countries, censuses are usually conducted but data access may be restricted due to confidentiality and privacy reasons. In the United States for instance, livestock data are only provided as a number of animal per county (Tildesley et al., 2010). In LMICs, censuses are not systematically carried out and data from occasional censuses are available as total numbers of animals per administrative units. However, the resolution at which these are provided vary greatly across countries (Robinson et al., 2011).
Downscaling procedures that increase the spatial resolution of existing livestock data, such as the Gridded Livestock of the World (GLW) are available (Wint et al., 2007; Robinson et al., 2014; Gilbert et al., 2018). However, these models predict livestock as a continuous, gradually varying, density of animals per pixel in a raster. Such data do not provide information at farm level, and densities of animal per pixel do not specify any information on the production systems. As production intensifies, livestock are fed on imported feed concentrates rather than on crops produced locally. Livestock production is hence less determined by environmental characteristics of the surrounding lands (Naylor, 2005; Steinfeld et al., 2006). Moreover, the spatial clustering of intensive farms is poorly explained by spatially continuous surfaces (Van Boeckel et al., 2012; Robinson et al., 2014). As suggested by Gilbert et al. (2015), a method accounting for the spatial clustering of intensive farms would be required to predict their distribution. In a previous study, we developed farm location models in Thailand, based on point pattern analysis methods (Chaiban et al., 2019). Other authors have predicted cattle population on individual farms in New Zealand with a zero-inflated Poisson regression, but only carried out an internal validation (van Andel et al., 2017). Burdett et al. (2015) modelled farm locations based on detailed geographical information as well as total animal and number of farms at county level. In this study, we present a methodology to predict both the location and the size of farms even where spatial information is scarce and which would account for the spatial clustering of farms.
We selected four countries to develop our farm distribution models, i.e. Nigeria, Thailand, Argentina and Belgium (hereafter referred to as the “sample countries”). These sample countries were assumed to have different levels of intensification as they have increasing levels of Gross Domestic Product (GDP) in the following order: Nigeria, Thailand, Argentina and Belgium. The GDP per capita (in purchasing power parity) is strongly associated with the level of intensification of livestock production of a country (Robinson et al., 2011; Gilbert et al., 2015). As farms tend to be more clustered as production intensifies (Abdalla et al., 1995; Roe et al., 2002a), we hypothesized that intensive farms would be spatially more clustered with increasing GDP. In addition to being at different stages of the intensification process, these countries present different ecological, demographic and socio-economic conditions. Developing farm distribution models (FDMs) along a wide range of conditions allows to assess how these model compare to each other in very different settings.
We first aimed to investigate the spatial distribution of commercial farms along the spectrum of intensification. Second, we built FDMs for intensive farms while accounting for the spatial clustering of farms. We then assessed the capacity of our modelling procedure to reproduce the farm distribution observed in Bangladesh, where farm locations and sizes were available.
Section snippets
Response variable and predictors
The modelling procedure used data on size (number of chickens per farm) and location of all intensive chicken farms in each country. As these data arose from different agricultural censuses, the criteria used to consider a farm as "intensive" may differ. Besides, the number of chickens recorded for each farm differed from one country to another. It can represent the number of chickens present at the time of the census or the maximum capacity. Details on the different datasets are presented in
Characterisation of spatial distribution of intensive farms
Intensive chicken farms were clustered in all countries. However, the level of clustering differed across countries, as highlighted by the L-function values (Fig. 2). Chicken farms in Argentina had the most clustered distribution. In Nigeria and in Thailand, chicken farms were also clustered but less so than in Argentina. While still clustered, chicken farms in Belgium were closer to a random distribution.
Internal validation - Farm location modelling
The combinations of predictors showing the lowest AIC in the iLGCP model (SM Fig. S2) by
Discussion
The paper aimed to (i) study the distribution of intensive farms along a gradient of intensification and (ii) investigate farm distribution models that would predict both location and size of chicken farms, while accounting for spatial clustering of farms in data-poor countries. Our model could reproduce only partly the location and size of farms in the observed data and poorly in external validation on Bangladesh. Despite some disappointing results, we believe our results identify potential
Conclusion
We attempted to develop a FDM approach predicting both the locations of individual farms and the flock sizes in the shape of a set of potential farm distributions in a given country. This methodology still has major issues, especially when extrapolated to a new area, which was the main purpose of the procedure. Further improvements should include further socio-economic and demographic covariates, and methodological modifications. Extrapolation should likely only be considered between countries
Acknowledgement
Authors acknowledge the “Fonds pour la formation à la Recherche dans l'Industrie et dans l'Agriculture” (FRIA) from the FNRS, the FNRS PDR “Mapping people and livestock” (PDR 477 T.0073.13), and the Walloon Institute for Sustainable Development (FRFS-WISD) PDR “Mapping livestock’s transition” (PDR-WISDX302317F), for supporting this project. This research also benefited from computational resources provided by the supercomputing facilities of the Université catholique de Louvain (CISM/UCL) and
References (54)
- et al.
Determining global population distribution: methods, applications and data
- et al.
Point pattern simulation modelling of extensive and intensive chicken farming in Thailand: accounting for clustering and landscape characteristics
Agric. Syst.
(2019) - et al.
Geographical determinants and environmental implications of livestock production intensification in Asia
(2005) - et al.
Mapping invasive plants using hyperspectral imagery and Breiman Cutler classifications (randomForest)
Remote Sens. Environ.
(2006) - et al.
Predicting farm-level animal populations using environmental and socioeconomic variables
Prev. Vet. Med.
(2017) - et al.
Predicting the distribution of intensive poultry farming in Thailand
Agric. Ecosyst. Environ.
(2012) Veterinary drug usage and antimicrobial resistance in bacteria of animal origin
Basic Clin. Pharmacol. Toxicol.
(2005)- et al.
What we know about historical trends in firm location decisions and regional shifts: policy issues for an industrializing animal sector
Am. J. Agric. Econ.
(1995) - et al.
Spatial Point Patterns: Methodology and Applications With R
(2015) Discussion on Dr Ripley’s paper
J. R. Stat. Soc. Ser. B
(1977)
Improved prediction accuracy for disease risk mapping using Gaussian process stacked generalization
J. R. Soc. Interface
Simulating the distribution of individual livestock farms and their populations in the United States: an example using domestic swine (Sus scrofa domesticus) farms
PLoS One
Correlation between veterinary antimicrobial use and antimicrobial resistance in food-producing animals: a report on seven countries
J. Antimicrob. Chemother.
GHS Built-Up Grid, Derived From Landsat, Multitemporal (1975-1990-2000-2014), R2018A
Downscaling livestock census data using multivariate predictive models: sensitivity to modifiable areal unit problem
PLoS One
World Livestock: Transforming the Livestock Sector Through the Sustainable Development Goals
Mapping global cropland and field size
Glob. Change Biol
High resolution population distribution maps for Southeast Asia in 2010 and 2015
PLoS One
Income disparities and the global distribution of intensively farmed chicken and pigs
PLoS One
Global distribution data for cattle, buffaloes, horses, sheep, goats, pigs, chickens and ducks in 2010
Sci. Data
High-resolution global maps of 21st-Century forest cover change
Science
Spatial dynamics of the livestock sector in the United States: do environmental regulations matter?
J. Agric. Resour. Econ.
The impact of surveillance and control on highly pathogenic avian influenza outbreaks in poultry in Dhaka division, Bangladesh
PLOS Comput. Biol.
Zoonosis emergence linked to agricultural intensification and environmental change
Proc. Natl. Acad. Sci.
Gridded global datasets for gross domestic product and Human Development Index over 1990–2015
Sci. Data
Positive and negative agglomeration externalities: arbitration in the pig sector
Spat. Econ. Anal.
Rural livelihood diversification in Sub-Saharan Africa: A literature review
J. Dev. Stud.
Cited by (3)
EPINEST, an agent-based model to simulate epidemic dynamics in large-scale poultry production and distribution networks
2024, PLoS Computational Biology