Problem description and model formulation
Assuming that crop demands is the total demands for the whole period; plots are homogeneous, thus the productivity of fertilizer is the same;
The plan must be sustainable for over the long term not just for one year.
For agricultural production, a vegetable rotation is a planting calendar presenting how to arrange the sequence of vegetables on a plot in the given time horizon. The problem can be described as follows:
A given set of C vegetables C 1,2,…, C is to be cultivated on a given set of H
plotsH1,2,…,Hwithin a predetermined durationTof T periods, T=1,2,…T.
The H plots are assumed to be homogeneous. Each vegetable i , iC , is planting repeatable and has a production duration ti It including the time for preparing the soil and
harvesting. In addition, all of the crops are assumed to have only one harvest occurring in the last period of its production duration. The fixed production cost of vegetable i , i C , on plot h without using chemical inputs at the period t is ci,t,h , and this cost is charged at the period
vegetable i planted. The vegetable i , iC , belongs to F botanic families, F = 1, 2, …, F , and two vegetables from the same botanic family cannot be cultivated on
one plot immediately one after the other. Besides, vegetables from the same botanic family also cannot be cultivated at the same time on adjacent plots.
For example: There are four vegetables C1, C2, C3 and C4 with different planting date and production duration. C1 and C3 are from brassicaceae botanic family, C2 and C4 are from solanaceae botanic family and cucurbitaceae botanic family, respectively. Green manure and fallow period are takes as two special vegetables with one month production duration and all year round planting date. Fig. 1 shows a sustainable vegetable rotation schedules for two adjacent plots. For plot 1, fallow period is planted in the first month, C1, green manure and C3 are orderly planted. For plot 2, C2 is cultivated in the first month and harvested in the fifth month, C4, fallow period and green manure are orderly planted after that. The vegetable rotation schedules for plot 1 and plot 2 satisfy the constraints of (a)-(c).
Plot 1 Plot 2
Due to the increasing concerns about sustainable development in agricultural production, green manures and fallow periods are suggested to be respected in each rotation period (Santos et al., 2010; Santos et al., 20111; Santos et al., 2015; Alfandari et al., 2014). In the model, U
unit green manures iG=C +1,…,Gand V unit fallow periods must be arranged in each rotation on plot h H . Green manures and fallow period are taken as special vegetables.
The production duration of green manures and fallow period is assumed to one period, indicatingthatti 1foriC1,…,CG1.ThesetofvegetablesC1,2,…,C,
green manures G=C +1,…, G , and fallow periods all belong to the set of
1 Santos Lana Mara R. dos, Philippe Michelon, Marcos Nereu Arenales, Ricardo Henrique Silva Santos. Crop rotation scheduling with adjacency constraints, Annals of Operation Research, 2011,190:165-180.
F
C1
G
C3
C2
C4
F
G
1 2 3 4 5 6 7 8 9 10 11 12 Fig.1. Sustainable vegetable rotation schedules for plot 1 and plot 2
1
I=1,2,…, C G +1. If vegetablei,iI , is planted on plot h at period tT , decision variable xi,t,h takes the value 1, otherwise, takes the value zero.
Inaddition,chemicalinput jJ=1,2,…, J(e.g.fertilizer/pesticide)playsavitalrole in vegetable growth. Under desirable level, unit chemical input j can boost additional output of vegetable i , iC , on unit area of plot h by ai, j,h units. However, if the quantity of
chemical inputs exceeds the maximum admissible level, it will result in many side-effects such as water pollution, damage to human health (Santos et al., 2010; Radulescu et al., 20142), soil contamination and emission of NO2 into the air (Wu, 20113). Thus, in our model, the quantity of chemical inputs j for vegetable i cannot exceeds the maximum admissible level of Li, j and
cost cfp j will be charged for buying unit chemical input j . Green manures and fallow period do not need any chemical inputs, while fixed cost cgi,h will be charged for
i C 1,…, C G 1 per unit area of plot h .
In practice, the total production costs in each rotation must be smaller than budget B , and the production quantity must be meet the demand Di .
Notation
i Indexes of vegetables, green manures and fallow period, where i=,,…,C denotes vegetables, C is the set of all vegetables; i= C 1,…, G denotes green manure, G is the set of green manures, where I = , ,…, G .
j Index of chemical inputs, where j =, ,… , j J , where J is the set of chemical
inputs.
t h
Index of periods, where t=,,…,tT , whereT is the set of time periods.
Index of plots, where h , , … , h H , where H is the set of plots. Parameters
F p
Si aith bijh
cith
Set of botanic family, p N , where N denotes the number of botanic family, and NC.
Set of time periods in which vegetable i can be cultivated.
Quantity of vegetable i produced on unit area of plot h at period t without chemical
inputs, where iC ,tT and hH .
Quantity of vegetable i produced on unit area of plot h with unit chemical input j ,
whereiC, jJand hH.
Production cost of vegetable i on plot h at the period t without chemical inputs, where iI, tT andhH.
2 Radulescu Marius, Constanta Zoie Radulescu, Gheorghita Zbaganu. A portfolio theory approach to crop planning under environmental constrains, Annals of Operation Research, 2014,219:243-264.
3 Wu Yanrui. Chemical fertilizer use efficiency and its determinants in China’s farming sector. China Agricultural Economic Review, 2011, 3: 117-130.
2
cfpj Cost of unit chemical input j , where j J .
pit Price of unit vegetable i at period t , where i C and t T .
pdi Production duration of vegetable i from preparing the land to being harvest, where
iC.
Di Demand of vegetable i , where i C .
B Available budget for each rotation.
Lij Maximum admissible level of chemical input j for vegetable i , where i C , and
jJ.
M A large number.
tg Minimum time delays for green manures and fallows.
eh,h’ If plot hand plot h’are adjacent, eh,h’ takes the value of one, otherwise taking the
value of zero. The value of eh,h’ depends on the distribution of plots. Decision variables
yijth Quantity of chemical input j applied to vegetable i on plot h at period t , where iC, jJ , tT and hH .
xith xith if vegetable, green manure or fallow period i is planted on plot h at period t ; xith otherwise.WhereiI,tTandhH.
zith Production quantity of vegetable i on unit area of plot h at period t , where i C , tTand hH.
Model formulation
Max pitzith cfpy
iI tT hH zi,tpdi,h
t ‘ pdi 1
tt’ jJ
ijh ijth i
jJ
j itjh iI tT hH
cx ith ith
zith Di,iC,tSi,tpdi TandhH (4) tT hH
jJ tt
iC tT hH
subject to: cx+
ith ith
(1) (a x b y ),iC,tSandhH (2)
iC jJ tT hH
cfpy B j ijth
ith ith
zi,tpdi,h Mxith,iC,tSi,tpdi TandhH (3)
t pdi 1
y Mx ,iC,tS,tpdT,andhH (5)
ijt h ith i i
yijth Lij,iC,jJ,tTandhH (6)
xith ,tT and hH iI
(7) x x M,tT,tpdTandhH (8)
iI tt
t pdi 1
it h ith i
3
xith U ,hH (9) iG tT
xith xith’ eh,h’,pN,tSi,h,h’H,hh’ iF(p)
(12)
xith ,,iI,tT and hH (13)
yijth and zith ,iC,tSi,jJ and hH (14)
The objective function is to maximum the profit. Constraint (1) ensures that total production costs are no more than the budgets. Constraint (2) defines the yield of vegetable i planted in period t for unit area of plot h . Constraint (3) enables that the yield of vegetable i in harvest period t+pdi is zero if it is not planted in period t on plot h . Constraint (4)
ensures that yield of vegetable i must satisfy its demand. Constraint (5) guarantees that the total quantity of fertilizer for vegetable i in production duration is zero if it is not planted in period t on plot h . Constraint (6) ensures that the quantity of chemical input j is under the maximum admissible level. Constraint (7) ensures that vegetable i , i I cannot occupy the land at the same time in each rotation. That is, at most one crop is planted at each period and each plot. Constraint (8) guarantees that if a crop is planted, no other crop (including green manure) can be planted during the time interval of production until the crop being harvested. Both constraint (9) and Constraint (10) enforce that one green manure must be planted in each rotation. Constraint (11) forbids that vegetables from the same botanic family are cultivated immediately one after the other, while vegetables from the different botanic family, green manures and fallow period are allowed. Constrain (12) ensures that vegetables from the same botanic family are not cultivated on the adjacent plots. Constrain (13) and (14) enforce the restrictions on
decision variables.
t’tg x
=tT, ttgTand hH (10) xith,pN,tSi ,tpdiTandhH (11)
C G ,t,h tt’
t’pdi iFp tt’
4