Bulletin of Mathematical Biology (2014) 76:997-1016. Preprint Manuscript, DOI: 10.1007/s11538-014-9943-9 Mathematical Programming Models for Determining the Optimal Location of Beehives Maica Krizna A. Gavina a , Jomar F. Rabajante a,c, 1 and Cleofas R. Cervancia b,c a Institute of Mathematical Sciences and Physics, University of the Philippines Los Ba˜ nos, Laguna, 4031 Philippines b Institute of Biological Sciences, University of the Philippines Los Ba˜ nos, Laguna, 4031 Philippines c University of the Philippines Los Ba˜ nos Bee Program, University of the Philippines Los Ba˜ nos, Laguna, 4031 Philippines Abstract. Farmers frequently decide where to locate the colonies of their domesticated eusocial bees, especially given the following mutually exclusive scenarios: (1) there are limited nectar and pollen sources within the vicinity of the apiary that cause competition among foragers; and (2) there are fewer pollinators compared to the number of inflorescence that may lead to suboptimal pollination of crops. We hypothesize that optimally distributing the beehives in the apiary can help address the two scenarios stated above. In this paper, we develop quantitative models (specifically using linear programming) for addressing the two given scenarios. We formulate models involving the following factors: (1) fuzzy preference of the beekeeper; (2) number of available colonies; (3) unknown-but-bounded strength of colonies; (4) probabilistic carrying capacity of the plant clusters; and (5) spatial orientation of the apiary. Keywords: beekeeping, pollination management, Apis mellifera, stingless bees, bumblebees, crop 1 Corresponding author, E-mail: [email protected]1
31
Embed
Mathematical Programming Models for Determining the Optimal Location of Beehives
This document is posted to help you gain knowledge. Please leave a comment to let me know what you think about it! Share it to your friends and learn new things together.
Transcript
Bulletin of Mathematical Biology (2014) 76:997-1016.
M. K. A. Gavina et al. Determining the Optimal Location of Beehives
pollination
AMS subject classification: 92B05, 90B80
1. Introduction
Pollinators, such as bees, contribute to biodiversity maintenance, especially of plants requiring
outcrossing [28, 38, 60]. Bees can serve as indicator of environmental stress [36], and are important
in food production [11, 29, 49, 54]. Hived bee colonies have been used to pollinate orchard plants,
such as mangoes and apples [26, 52], coffee [37, 46], almonds [13, 49], cucurbits [15, 22], Chinese
cabbage and radish [16], and mangroves [4, 56]. Bee colonies are also managed for the production
of honey, pollen and beeswax.
By looking at the current trend, populations of some bee species are inadequate to sustain
global demands [1, 2, 31, 70]. The decline in the number of bee pollinators is due to various
factors, such as prevalence of diseases [7, 33, 40] and pesticides [23, 45, 53], and starvation [69].
Minimizing the exposure to these factors can reduce bee stress and can avoid colony collapse
[20, 30, 50, 51, 72].
Best practices in beekeeping and pollination management help increase agricultural yield [3,
11, 44, 73]. Placement of hives in farmlands is crucial in optimizing the efficiency of bees in
pollination and production of hive products [6, 12, 14, 35, 43]. It is better to have a systematic way
of locating new beehives or relocating existing beehives than haphazardly placing them in any area.
The hives should be positioned in an apiary in such a way that competition for forage is minimized
and access of bees to target flowers are facilitated. Various factors need to be considered, such as
strength of available colonies [21, 59] and carrying capacity of the foraging sites (plant clusters or
food patches) [55, 59].
Operations Research (OR) techniques have been employed to solve combinatorial problems
2
M. K. A. Gavina et al. Determining the Optimal Location of Beehives
[5, 10, 27, 67], such as choice selection and facility location. Determining the optimal place-
ment of bee colonies in a large apiary can be regarded as a combinatorial problem solvable using
mathematical programming tools. In this study, we formulate mixed-integer linear programs that
can be used in determining the best location of beehives taking into account the preferred loca-
tion of the beekeeper, number and strength of available colonies, carrying capacity of the plant
clusters, maximum flight distance that the bees can travel, and spatial orientation of the apiary. We
present deterministic models and then modify them by including fuzzy preference of the beekeeper,
unknown-but-bounded colony strength, and probabilistic carrying capacity of the plant clusters.
In this paper, colony strength denotes the number of foragers, and the carrying capacity of the
plant cluster denotes the number of foragers it can support or sustain (in view of the abundance of
pollen and nectar sources). The beekeeper’s preferred location may depend on the proximity of the
hives to his house or workplace, water sources, incidence of predation, insecticides or diseases, and
target plants to be pollinated. Environmental factors, type of bee species and floral choice of bees
should be considered in the modeling process [24, 48, 57, 58, 61]. The models are applicable to
a variety of social insects (such as honeybees, stingless bees and bumblebees) and to any farming
designs (such as forest farming, monoculture and permaculture) that fit our model assumptions.
The mathematical models are grouped into two sections considering the following scenarios
(1) minimizing competition among foragers due to limited food sources, and (2) maximizing the
use of the species as pollinators given small number of foragers. Note that these two scenarios are
mutually exclusive, that is, the first scenario happens when food from plant patches are limited,
and the latter scenario happens when there are more plant patches to pollinate compared to the
available number of bee foragers. The latter scenario is intended for crop pollination management.
In addition, the model can be used in determining the minimum number of hives to support optimal
pollination of crops.
3
M. K. A. Gavina et al. Determining the Optimal Location of Beehives
2. Minimization of Competition for Limited Resources
In this section, we present mathematical models for identifying the optimal spatial distribution of
beehives in an apiary for minimizing competition among colonies. In determining the best location
site for a beehive, we need to consider (1) the distance of the hive from the food sources and (2) the
maximum flight distance that the bee species can travel [8, 39, 64, 68]. These two considerations
can be incorporated in the network characterizing the spatial orientation of the apiary (see Figure
1). We assume that a forager can only visit a plant cluster within the maximum flight distance. We
can make the latter assumption stricter if we want to minimize the flight stress of the bees.
Figure 1: A sample network characterizing the spatial orientation of the apiary: The rectangularnodes represent the possible location sites; the circular nodes represent the plant clusters; and anedge denotes that the connected plant cluster is contained in the area that the foragers from theconnected location site can visit (i.e., within the maximum flight distance of the bees).
Our goal is to minimize the “overpopulation” of bees in a locale. In this paper, we define
“overpopulation” as the number of foragers that cannot anymore be accommodated by any plant
cluster in the foraging area. If the total number of bees foraging in nearby plant clusters exceeds
the total carrying capacity of the plant clusters then overpopulation arises. This overpopulation is
caused by the mismatch in the capacity of the habitat and the population size. Likewise, overpopu-
4
M. K. A. Gavina et al. Determining the Optimal Location of Beehives
lation stimulates competition among foragers. We hypothesize that minimizing overpopulation by
optimally distributing the beehives will lessen this competition.
In this paper, we use the following notations for the decision variables and parameters. Strength
of a colony denotes the number of foragers in the colony, and the carrying capacity of a plant cluster
denotes the number of foragers that the plant cluster can accommodate.
Decision Variables:
xij : portion of foragers in location site i that can be accommodated by plant cluster j, xij ∈ R⊕
Zki =
1 if colony k is located to site i
0 otherwise
Ei : portion of foragers in location site i that cannot be accommodated by nearby plant clusters,
Table 2: Carrying capacity of plant clusters (×10, 000). Values are derived by determining thenumber of flowers and the number of foragers that can be accommodated by each flower. Supposea survey of carrying capacity is done for 12 weeks. “PC j” denotes Plant Cluster j and “Wtj”denotes Week t corresponding to PC j.
Table 5: Optimal solution to the mathematical program. Pi ×10, 000.
by increasing or decreasing the number of colonies in the model. We do this to determine the
reasonable number of hives to support optimal pollination. In this example, we have P1 = 0.05.
Suppose the beekeeper can avail an additional colony (Colony 7) with strength b7 −∆b7 = 1. The
solution to the corresponding mathematical program suggests that colonies 1, 3, 4 and 5 be placed
at site 1; and colonies 2, 6 and 7 be placed at site 3. The additional colony results in Pj = 0 ∀j.
3.3. An alternative model
This model is specific for cases where we assume that food from plant cluster j is allocated to
nearby location sites. For simplicity, we assume that the beekeeper prefers to evenly allocate the
carrying capacity of a plant cluster to all connected location sites (if the number of bee colonies
is enough, this will force each location site to have at least one bee colony). The major difference
between this alternative model and the previous model (with IFD) is the construction of Constraint
1.
Let Pi represent the approximate carrying capacity of nearby plant clusters inadequately uti-
lized by colonies located in site i, and γj be the number of location sites connected to plant cluster
j less the number of location sites connected to hazardous places.
Objective function : Maximize∑j∈S+
wjMgj −∑
i∈Dj , j∈S+
Pi (3.7)
20
M. K. A. Gavina et al. Determining the Optimal Location of Beehives
subject to
Constraint 1 :n∑k=1
bkZki + Pi ≥∑j∈Ci
yjγjgj, ∀i ∈ Dj, j ∈ S+ (3.8)
Constraint 2 :∑
i∈Dj , j∈S+
Zki = 1, ∀k = 1, 2, ..., n. (3.9)
Take note that Pj and Pi are not equivalent since we use different assumptions in their calcula-
tions. Figures 2 and 3 illustrate the assumptions used in the model with IFD and in this alternative
model, respectively. In the model with IFD, total strength of colonies located in a certain site is
distributed among the connected plant clusters; while in this alternative model, carrying capacity
(food) of a plant cluster is allocated to the connected possible location sites. Moreover, the auxil-
iary deterministic program that represents the stochastic version of this alternative model may lead
to a nonlinear program [42] (because of the term∑
j∈Ciyj/γj in Constraint 3.8).
Figure 2: A sample network showing the distribution of the total strength of colonies to the con-nected plant clusters following the Ideal Free Distribution theory.
21
M. K. A. Gavina et al. Determining the Optimal Location of Beehives
Figure 3: A sample network showing the allocation of the carrying capacity of a plant cluster tothe connected possible location sites. Site 3 is not a possible location site since it is connected to aplant cluster situated in a hazardous area (e.g., with insecticides).
4. Concluding Remarks
In this paper, we formulated two auxiliary mixed-integer linear programs for (1) minimizing “over-
population” of bee species and for (2) maximizing pollination of flowers, where some parameters
are either fuzzy or stochastic. Keen observation of the spatial orientation of the apiary is necessary
in writing the appropriate mathematical program, and a close collaboration between beekeepers
and mathematical modelers is indispensable. The derived mathematical program can be solved us-
ing any Operations Research (OR) software capable of dealing with mixed-integer programming
problems. After obtaining the results, it is advised to do sensitivity analysis by identifying the
changes in the optimal solution when parameter values are perturbed. Our model can be extended
by considering a weighted network (spatial orientation of the apiary), where a weight is a function
of distance between nodes.
The optimal solution derived from the model can help beekeepers in selecting the best beehive
location sites. The optimal distribution of bee colonies in an apiary can help minimize the foraging
22
M. K. A. Gavina et al. Determining the Optimal Location of Beehives
stress of the bees, which may result in increased production of honey, pollen and wax. Optimal
distribution of beehives can also aid in augmenting the yield of high-value crops. However, the
computed optimal solution is not permanent. When ecological season changes, a new spatial
distribution of beehives may be necessary, bearing in mind the cost of transferring beehives.