Research article Topical Sections

Vegetation structure of bio-belts as agro-environmentally-climatic measures to support biodiversity on arable land: A case study

  • Loss of biological diversity is one of the greatest challenges that our civilization must face nowadays. Reaction to the diminishing biodiversity of agricultural landscapes is various measures promoting free-living organisms. The study deals with the vegetation composition and structure of agro-environmental-climatic measures applied on arable land in operating conditions (intensively farmed regions of the Czech Republic). Additional study focus was applied to a popular measure of the feeding bio-belts. Bio-belts are not only hiding places for free-living animals but can provide them a rich food offer in the period from the harvest of main crops until winter. Thanks to the bio-belts, the landscape gains in biodiversity, and sloping sites can be protected from soil erosion. The vegetation of land parts used as bio-belts was assessed using phytocoenological relevés. Dominant plant species sown in the bio-belts were Avena sativa, Panicum miliaceum, Brassica oleracea var. acephala, Fagopyrum esculentum, Phacelia tanacetifolia, and Pisum arvense. Apart from the sown plants, there were also weeds occurring in the bio-belts, of which the most abundant were Chenopodium album, Amaranthus retroflexus, Setaria verticillata, Cirsium arvense, Equisetum arvense, etc. Risks connected with the realization of feeding bio-belts in respect of weeds occurring on arable land are negligible. Weeds from bio-belts have only a limited potential to spread to adjacent arable land. A potential spreading of weeds from the bio-belts to adjacent arable land was not demonstrated. On the contrary, thanks to its composition, the vegetation of bio-belts has the potential to extend the food offer for animals. Thus, bio-belts are useful for supporting biodiversity in regions intensively used for agriculture.

    Citation: Helena Hanusová, Karolína Juřenová, Erika Hurajová, Magdalena Daria Vaverková, Jan Winkler. Vegetation structure of bio-belts as agro-environmentally-climatic measures to support biodiversity on arable land: A case study[J]. AIMS Agriculture and Food, 2022, 7(4): 883-896. doi: 10.3934/agrfood.2022054

    Related Papers:

    [1] Tian Xue, Liu Li, Liu Shuang, Du Zhiping, Pang Ming . Path planning of mobile robot based on improved ant colony algorithm for logistics. Mathematical Biosciences and Engineering, 2021, 18(4): 3034-3045. doi: 10.3934/mbe.2021152
    [2] Xuewu Wang, Bin Tang, Xin Zhou, Xingsheng Gu . Double-robot obstacle avoidance path optimization for welding process. Mathematical Biosciences and Engineering, 2019, 16(5): 5697-5708. doi: 10.3934/mbe.2019284
    [3] Xue Li, Lei Wang . Application of improved ant colony optimization in mobile robot trajectory planning. Mathematical Biosciences and Engineering, 2020, 17(6): 6756-6774. doi: 10.3934/mbe.2020352
    [4] Baoye Song, Shumin Tang, Yao Li . A new path planning strategy integrating improved ACO and DWA algorithms for mobile robots in dynamic environments. Mathematical Biosciences and Engineering, 2024, 21(2): 2189-2211. doi: 10.3934/mbe.2024096
    [5] Liwei Yang, Lixia Fu, Ping Li, Jianlin Mao, Ning Guo, Linghao Du . LF-ACO: an effective formation path planning for multi-mobile robot. Mathematical Biosciences and Engineering, 2022, 19(1): 225-252. doi: 10.3934/mbe.2022012
    [6] Yuzhuo Shi, Huijie Zhang, Zhisheng Li, Kun Hao, Yonglei Liu, Lu Zhao . Path planning for mobile robots in complex environments based on improved ant colony algorithm. Mathematical Biosciences and Engineering, 2023, 20(9): 15568-15602. doi: 10.3934/mbe.2023695
    [7] Won Man Park, Young Joon Kim, Shaobai Wang, Yoon Hyuk Kim, Guoan Li . Investigation of lumbar spine biomechanics using global convergence optimization and constant loading path methods. Mathematical Biosciences and Engineering, 2020, 17(4): 2970-2983. doi: 10.3934/mbe.2020168
    [8] Zhen Yang, Junli Li, Liwei Yang, Qian Wang, Ping Li, Guofeng Xia . Path planning and collision avoidance methods for distributed multi-robot systems in complex dynamic environments. Mathematical Biosciences and Engineering, 2023, 20(1): 145-178. doi: 10.3934/mbe.2023008
    [9] Teng Fei, Xinxin Wu, Liyi Zhang, Yong Zhang, Lei Chen . Research on improved ant colony optimization for traveling salesman problem. Mathematical Biosciences and Engineering, 2022, 19(8): 8152-8186. doi: 10.3934/mbe.2022381
    [10] Jinzhuang Xiao, Xuele Yu, Keke Sun, Zhen Zhou, Gang Zhou . Multiobjective path optimization of an indoor AGV based on an improved ACO-DWA. Mathematical Biosciences and Engineering, 2022, 19(12): 12532-12557. doi: 10.3934/mbe.2022585
  • Loss of biological diversity is one of the greatest challenges that our civilization must face nowadays. Reaction to the diminishing biodiversity of agricultural landscapes is various measures promoting free-living organisms. The study deals with the vegetation composition and structure of agro-environmental-climatic measures applied on arable land in operating conditions (intensively farmed regions of the Czech Republic). Additional study focus was applied to a popular measure of the feeding bio-belts. Bio-belts are not only hiding places for free-living animals but can provide them a rich food offer in the period from the harvest of main crops until winter. Thanks to the bio-belts, the landscape gains in biodiversity, and sloping sites can be protected from soil erosion. The vegetation of land parts used as bio-belts was assessed using phytocoenological relevés. Dominant plant species sown in the bio-belts were Avena sativa, Panicum miliaceum, Brassica oleracea var. acephala, Fagopyrum esculentum, Phacelia tanacetifolia, and Pisum arvense. Apart from the sown plants, there were also weeds occurring in the bio-belts, of which the most abundant were Chenopodium album, Amaranthus retroflexus, Setaria verticillata, Cirsium arvense, Equisetum arvense, etc. Risks connected with the realization of feeding bio-belts in respect of weeds occurring on arable land are negligible. Weeds from bio-belts have only a limited potential to spread to adjacent arable land. A potential spreading of weeds from the bio-belts to adjacent arable land was not demonstrated. On the contrary, thanks to its composition, the vegetation of bio-belts has the potential to extend the food offer for animals. Thus, bio-belts are useful for supporting biodiversity in regions intensively used for agriculture.



    The logistics industry has strongly promoted the development of the national economy under the impetus of information technology. However, this has led to high energy consumption and carbon emissions. Only by insisting on economic development, clean development and safe development, that is, sustainable development, can we achieve sound and rapid economic development. Sustainable development is an economic model that focuses on long-term development, while energy conservation and emission reduction are necessary measures for sustainable development, and the two complement each other.

    Energy consumption of the logistics industry in China accounts for 10% of the total national energy consumption [1], which is not conducive to the development concept of green GDP [2]. Green GDP is the cost of economic loss caused by environmental pollution, degradation of natural resources, poor education, uncontrolled population, and poor management of GDP [3]. Therefore, the logistics industry is facing the formidable pressure of energy saving and emission reduction. The logistics system established with the help of modern information technology can effectively improve operational efficiency, and achieve cost reduction, energy saving, and consumption reduction [4]. Specifically, multi-center joint distribution can break the limitation of the distribution area, reduce the round trip, improve the efficiency of distribution, and achieve the goal of energy saving and emission reduction, which is an important choice for the low-carbon and green development of logistics enterprises [5].

    Malik and Kim [6] considered the production rate affects carbon emissions generation in production, i.e., generally, higher production rates result in more emissions. Sarkar et al. Introduced a three-echelon sustainable supply chain model with a single supplier, single manufacturer, and multiple retailers. In this supply chain, the main task of the manufacturer is to reduce defective products and control carbon emissions to maintain sustainability. Numerical experiments find that the model obtains the global optimum solution at the optimum values of the decision variables and the reduction of carbon emission has been proved [7].

    Ghanbarpour and Gustafsson [8] proposed that customers' positive perceptions of firm actions do not directly impact financial earnings; however, they do impact earnings through customer satisfaction. Tayyab et al. analyze a cleaner multi-stage production management system for carbon emissions reduction and active participation in corporate social responsibility activities, while also advancing the system economically. According to the model results, the lower batch size is the optimal strategy in response to increasing holding and CSR activity costs, but otherwise, for the setup cost in a serial production system [9].

    Dantzig and Ramser [10] proposed the vehicle routing problem (VRP) in 1959. Since its inception, VRP has sparked considerable interest. In particular, the research findings on VRP considering carbon emission, time window and customer satisfaction have been fruitful. In order to solve minimum total travel time and minimum fuel consumption in VRP, Cinar et al. used a variable travel distance and load weight to construct the cumulative function, and used K-means clustering and C & W algorithm to form a two-stage algorithm to improve the algorithm search capability [11]. Xiao et al. considered vehicle type, speed, load weight, and time window in the vehicle path optimization process, and the optimization can reduce carbon emissions by 8% [12]. Deng et al. considered both the soft time windows and hard time windows and established a multi-modal transport route optimization model with the objectives of minimizing total cost and carbon emissions [13]. Zhang et al. introduced a cold chain logistics route optimization model considering the low-carbon economy, an ant colony optimization was used to avoid the influence of unreasonable parameter selection of algorithm performance [14]. Tao et al. proposed a distribution route optimization method to minimize carbon emission cost as well as comprehensive cost [15]. Zhang et al. proposed a multi-objective route optimization model of the instant distribution system taking into account total cost and customer time satisfaction, testing results showed that a slight rise in delivery cost can achieve more on-time delivery demand [16]. Li and Zhou presented a logistics distribution center location model which considered the influence of customer satisfaction, dynamic and static carbon emission, construction and operation cost, and proposed penalty cost under a time window to measure customer satisfaction [17].

    With the development of green logistics, a single distribution center can no longer meet the distribution needs, and scholars have gradually conducted research on the joint distribution model of multiple distribution centers. Rabbani et al. presented a deterministic vehicle routing problem model taking into account multiple middle depots [18]. Adelzadeh et al. studied the vehicle routing problem of multi-distribution center joint distribution under time window constraints, providing a reference value for enterprises to reduce their distribution costs and improve their service level [19]. Li et al. established an optimization model of joint distribution terminal distribution path with carbon emission, taking into account various costs with a time window, with the total distribution cost optimized. The model was solved using a genetic algorithm [20]. Golestani et al. investigated a green hub location problem with the objectives of minimizing the system's total cost and maximizing the quality of the delivered product to the customer [21].

    At present, there exists a number of literature on the research of carbon emission optimization, analyzing the role of low carbon and other factors in the distribution systems, and using intelligent optimization algorithms such as heuristic algorithms and exact algorithms [22,23,24] to solve similar VRP problems [25,26], and some also study the customer satisfaction problems in logistics and distribution processes such as centralized distribution, joint distribution, and cold chain distribution. However, there are few logistics distribution models that consider both carbon emissions and customer satisfaction. Logistics enterprises, can reduce carbon emissions through vehicle path optimization and achieve cost reduction. At the same time, economic benefits play a critical role. Therefore, it is meaningful to achieve both economic and environmental benefits for enterprises by maximizing customer satisfaction and reducing carbon emissions in the distribution process. Based on the above analysis, this paper presents a joint distribution network with multiple distribution centers, constructing a multi-objective planning model taking into account distribution cost, customer dissatisfaction, and carbon emission. An improved ant colony optimization (ACO) is proposed to solve the model by using an elite strategy and saving strategy to obtain an optimal solution.

    The logistics industry has a long history of high energy consumption and carbon emissions [27]. Reducing energy consumption and carbon emission in the logistics process helps enterprises to realize both economic and social benefits. This paper solves the problem of joint distribution with multiple distribution centers and customer time window requirements, i.e., it is an MDVRPTW problem. In the multi-centers joint distribution network, the customer order is no longer handled by a fixed distribution route, but the appropriate distribution vehicle and distribution route are decided according to the overall planning of the logistics system. This distribution method can provide rapid response to customer time window demand and maximize delivery efficiency. The sharing of vehicles and inventory resources can effectively solve the problems of urban traffic congestion and vehicles running empty. The vehicle selects the original distribution center for replenishment, then continues to deliver until all orders are completed, and finally returns to the original distribution center. In this paper, carbon emissions and customer dissatisfaction are firstly measured. The carbon emission in operation is calculated by the vehicle's actual load rate and vehicle driving distance, and the customer dissatisfaction is calculated by the fuzzy mathematical affiliation function. Secondly, the joint distribution model of multiple distribution centers considering carbon emission and time window is constructed. The model takes into account actual vehicle loading rate, vehicle capacity, driving speed, loading and unloading time, number of vehicles, and number of distribution centers in the distribution network. As a result, this paper presents a multi-objective path optimization model with the objectives of minimum cost, minimum customer dissatisfaction, and minimum carbon emission.

    Assumptions in this work are listed as follows.

    1) The location of the distribution center is known, the reserve capacity can well serve each demand point, and the transport capacity is able to complete the distribution task.

    2) Each vehicle departs from its respective distribution center and returns to the same departure point after completing the distribution task.

    3) The vehicle type is consistent.

    4) Each node is served by one vehicle once only, and the node requirements can be satisfied in a one-time service.

    5) The time window requirements of the node must be satisfied.

    6) The loading quantity of a vehicle can satisfy the total demand in its corresponding path.

    7) Each vehicle is not reassigned in the middle of the distribution.

    The meaning of the parameters in the model are presented as follows.

    M: vehicle set in distribution centers.

    Vm: orders served by vehicle m.

    N: order set.|N|=n, there are n orders.

    i,j: indicates order number, also used to mark customer number of the order.

    dij: distance from customeri to j, dij=dji.

    xmij: decision variable, whether vehicle m serves customer j directly after serving customer i.

    ymi: decision variable, whether customer i is served by vehicle m or not.

    Q: vehicle loading capacity and capacity of all vehicles are the same.

    di: demand of customer i.d0=0, which means the demand of the distribution center is 0.

    TD: the maximum traveling distance of the vehicle.

    ti: the arrival time of vehicle m reaches customer i.

    μi(ti): dissatisfaction rate of customer i when served at time ti.

    pfuel: unit fuel price.

    cep: penalty cost per unit of time if the vehicle arrives early.

    clp: penalty cost per unit of time if the vehicle arrives late.

    C1: distribution cost.

    C2: customer dissatisfaction.

    C3: carbon emission.

    The costs generated by distribution are composed of vehicle driving fuel cost, and penalty cost if the vehicle arrives early or late.

    C11 presents vehicle driving fuel cost, which is calculated by following equation.

    C11=pfuelmMiVmjVmiε(Qij)dijxmij (1)

    The penalty cost consists of early and late arrival penalty cost. The early arrival penalty cost is Tem,i=max(0,aiti),mM,iV and the late one is Tlm,i=max(0,tibi),mM,iV. Therefore, the penalty cost can be calculated as following equation.

    C21=mMiVmymicepTem,i+mMiVmymiclpTlm,i (2)

    At last, the vehicle distribution cost can be obtained from formula (3).

    C1=C11+C21 (3)

    The VRPTW only considers how to complete distribution tasks with the least number of vehicles and the shortest driving distance to save cost. However, in the retail industry, with the diversification of cargo and the increase of real-time distribution requirements from customers, multiple small-batch distributions are becoming more and more popular. Customers make delivery appointments at any time, and they have a certain tolerance for early or late arrival. In the actual situation, untimely delivery, in the long run, will damage enterprises' reputations. Consequently, considering customer dissatisfaction with delivery efficiency is of great practical importance.

    The method of determining customer dissatisfaction is shown in Figure 1. ti is the vehicle arrival time at customer i. The agreed time window is [ETi,LTi]. If the vehicle arrives at the customer's during this period, customer dissatisfaction is 0. The tolerable early arrival time of the customer is ETiai, and the tolerable delayed arrival time isbiLTi. If the vehicle arrives at the customers within the tolerable time window, although the order still would be received, the customer would be dissatisfied depending on the degree of order advancement or delay. The customer's intolerable early arrival time isai, and the intolerable delay time is bi. If the vehicle arrives outside of the time period [ai,bi], customer dissatisfaction is 1.

    Figure 1.  Customer delivery time window.

    Therefore, based on above analysis, when customerireceives order at time ti, the dissatisfaction rate μi(ti) can be presented as follows.

    μi(ti)={1,ti<aiETitiETiai,aiti<ETi0,ETitiLTitiLTibiLTi,LTi<tibi1,ti>bi (4)

    Figure 2 shows customer dissatisfaction when receiving orders at different moments.

    Figure 2.  Customer dissatisfaction.

    Customer dissatisfaction consists of the dissatisfaction of all customers, i.e., total customer dissatisfaction rate is

    C2=iNμi(ti) (5)

    FCR during vehicle driving is closely related to the actual load factor of the vehicle, this factor is an important indicator to measure the effective utilization of the vehicle. The higher the actual load factor of the vehicle, the higher the cost of fuel consumption, but the lower the transportation cost, the higher the utilization rate of the vehicle. According to the FCR regression equation [28], the fuel consumption per unit distance can be calculated as below.

    ε=0.00556Q+0.254,Q[0t,22.4t] (6)

    Where, ε is the fuel consumption per unit distance, Q is vehicle load capacity, and Q is the maximum vehicle load. When the vehicle is empty, ε0=0.254, when it's fully loaded, ε=0.37944.

    Let actual load factor pi=QQ, FCR is displayed in Eq (7) [29].

    ε(Q)=ε0+(εε0)pi=ε0+(εε0)QQ (7)

    dij is the traveling distance between customer i to customer j, Qij is the vehicle load during this driving interval, ε(Qij) is FCR, and cfuelij is the fuel consumption, then cfuelij can be expressed as below [30].

    cfuelij=ε(Qij)dij=(ε0+(εε0)QQ)dij (8)

    Emissions of CO2 during vehicle operation and fuel consumption cfuelij on the combustion conversion factor F is linearly related [31]. Carbon emission is measured by Eq (9).

    C3=mMiVjV{i}xmijcfuelijF (9)

    Different fuels have different combustion conversion factors. This work assumes that minivans are used for distribution, according to literature [32], the combustion conversion factor F = 2.3, which means that the weight of CO2 released per liter of gasoline combustion is 2.3 kg.

    Based on the above sub-item analysis, the multi-objective optimization model for joint distribution considering carbon emission and customer dissatisfaction is constructed by using a dictionary order multi-objective planning method.

    Equation (10) is the multi-objectives.

    minL=(C1,C2,C3) (10)

    Equation (11) is the priority of different objectives, which can be adjusted according to the objective priority during decision-making.

    C1C2C3 (11)

    The loading capacity of the vehicle and the maximum driving distance are considered during distribution.

    s.t.

    iVm{j}xmij=ymj,jVm,mM (12)
    jVm{i}xmij=ymi,iVm,mM (13)
    mMymi=1,iN (14)
    iVmjVm{i}dijxmijTD,mM (15)
    iNdiymiQ,mM (16)
    iVmjVm{i}xmij|Vm|1,mM (17)
    tibi,iN (18)
    xmij{0,1},i,jN,mM (19)
    ymi{0,1},iN,mM (20)

    Constraint (12) indicates that vehicle m can only serve one customer i before serving customer j. Constraint (13) means that after serving customer i, vehicle m should serve customer j immediately. Constraint (14) indicates that each customer can only be served by one vehicle. Constraint (15) indicates that the total distance traveled by the vehicle does not exceed its maximum traveling limit. Constraint (16) presents that the sum of customer demands delivered by the vehicle does not exceed its maximum capacity. Constraint (17) is to avoid loops during delivery. Constraint (18) means the time constraint of the vehicle at each customer point.xmij and ymi are decision variables. When vehicle m serves customer j immediately after serving customer i, xmij=1, otherwise, xmij=0. When customer i is served by vehicle m, ymi=1, otherwise, ymi=0.

    This paper studies a multi-objective optimization problem, which is an NP-HARD problem, and heuristic algorithms are widely used in solving such problems [33]. Among them, the ant colony algorithm uses a distributed parallel computing mechanism, which has the advantages of positive feedback, good robustness, and obtaining a more satisfactory and feasible solution within an acceptable time frame. The basic idea of the ant colony algorithm: ants release pheromones along the way during crawling, and the pheromone will gradually weaken with time. When other ants pass the path, the pheromone is strengthened, and the ants tend to crawl in the direction of high pheromone intensity during the crawling process [34]. The more ants walk on a certain route, the higher the concentration of pheromone left on the route. A large number of ants form a positive information feedback mechanism, and finally, a route with the largest number of crawling ants is formed, which is also the optimal route [35].

    However, it has disadvantages of early maturity and low solution efficiency. Therefore, an improved ant colony algorithm is proposed in this paper.

    Considering both vehicle load and farthest driving distance in joint distribution is a new VRP. The taboo table design is different from the traditional solution. It needs to find the original warehouse for replenishment, and then continue. What's more, the vehicle has a limit of maximum traveling distance. In order to solve the problem described above, this paper designed a specific taboo table coding method of ACO, which is shown in Figure 3.

    Figure 3.  Coding scheme.

    For example, in Figure 3, numbers 23, 24 and 25 represent warehouses, and numbers 1 to 22 represent customers.

    The route of a certain ant contains two distribution vehicles, where the driving path of vehicle 1 is 23→2→11→15→5→7→23→13→17→1→23→18→21→9→10→16→23, while the driving path of vehicle 2 is 23→6→12→20→19→23→3→8→14→4→22→23. Distribution vehicle 1 departs from warehouse number 23 and delivers orders according to route 23→2→11→15→5→7. After completing the delivery task to customer 7, there are no more goods for the next customer, therefore, vehicle 1 goes to the original warehouse number 23 for replenishment, after reloading, it continues delivery along route 23→13→17→1. And goods are out again at customer 1, so vehicle 1 returns to the original warehouse number 23 for reloading, and then continues distribution until the completion of all orders assigned to vehicle 1. Vehicle 2 adopts a similar distribution method as vehicle 1. The driving routes of all vehicles constitute an ant's route, and the routes of all ants constitute a taboo table.

    The distribution center needs to select a vehicle for delivery orders. Since the demand of each customer is unequal, and the distance between two adjacent customers is different, in order to avoid uneven free and busy times of vehicles, the distribution center gives priority to the vehicle with lower busyness for delivery according to a certain probability.

    Because it is in accordance with the probability to choose the distribution vehicle, there will be some vehicles arranged for more distribution tasks, and some less. For example, vehicle 1 delivers more goods to customers far away, while vehicle 2 delivers fewer goods to customers nearer, then, vehicle 1 takes more time and vehicle 2 takes less. In order to avoid this problem, the busy level of distribution vehicles is determined in accordance with the current distribution time of each vehicle. Such as, if vehicle 1 spends more time, and vehicle 2 spends less, then vehicle 1 is busier than vehicle 2, and vehicle 2 will be given priority for new assignments. The vehicle selection rules are shown in Figure 4.

    Figure 4.  Vehicle selection rules.

    After completing vehicle selection, it is time to assign customers to be served by specific vehicles under constraints of maximum load capacity, the farthest driving distance, and other factors. When assigning customers, the transfer probability of the vehicle traveling from the current customer to the remaining customers is firstly calculated according to the indexes such as pheromone and visibility, then the customer with maximum transfer probability is selected. If the customer with the maximum transfer probability is not selected, then roulette is constructed according to the transfer probability to select the customer to be served. After completing customer selection, whether the customer's demand exceeds the maximum load capacity of the vehicle is judged, if it exceeds the maximum load capacity, then go back to the last step, if it still does not meet the maximum load capacity of the vehicle, then select the original warehouse to the current customer for replenishment, and continue to select the customer to visit according to the probability after replenishment. Customer selection rules are shown in Figure 5.

    Figure 5.  Customer selection rules.

    In the transfer rule of the basic ant colony algorithm, the transfer probability takes into account the pheromone concentration and visibility factors. After departure, whether a distribution vehicle visits point j from point i is determined by the probability. In order to accelerate the convergence speed of the ant colony algorithm, the maximum transfer probability customer point is selected with a certain probability, and if the probability is not satisfied, the next point to be visited is selected according to the probability of each point. The pseudo-random rule is: rand is a random number between 0 and 1, whenrand<=maxshift_p, the next point to be visited is selected according to the maximum probability, otherwise a probability roulette is constructed and the point to be visited is selected according to the roulette. The probability of transfer from point i to point j is:

    (21)

    In order to improve the searchability of the basic ant colony algorithm, the saving factor is added to the transfer probability by using the idea of the saving algorithm as a reference, and the transfer probability is presented in Eq (22).

    (22)

    allowedm indicates the set of points that m vehicle can be accessed from current point i. Still, α,β indicate the importance of pheromone, and the importance of visibility respectively, and τij indicates the pheromone concentration from point i to point j, ηij is the visibility from point i to point j. μij is the amount of savings of edge(i,j), i.e., μij=hib+hbjhij(bD). The savings of the edge reflects the importance of the path(i,j) [36].

    After the ant finishes constructing paths of all delivery vehicles, it needs to update the pheromones on the paths. Since the problem is a multi-objective optimization with cost as the first priority, and the cost is mainly affected by driving distance and full load rate, therefore, the pheromone update for the ant k should refer to travel distance Lk of ant k. The pheromone update equation is as the following equation.

    τkij=QLk (23)
    τij=(1ρ)τij+τkij (24)

    Where ρ(0<ρ<1) is the pheromone volatility factor, τkij is the pheromone increment of ant k from point i to point j, Q is a constant, and fk is the multi-objective weighted sum of ant k, and τij is the pheromone concentration from point i to point j.

    The pheromone of the ant gradually evaporates during walking, and leaves a new pheromones, this process guides the optimal solution, which will be gradually strengthened and kept stable through repeated operations.

    The ant colony algorithm relies on probability to guide the search when selecting the access point. It has been proved that the ant colony algorithm converges to the optimal solution with probability 1 after repeated iterations, and in order to speed up the convergence of the ant colony algorithm, the convergence of the algorithm can be accelerated by enhancing the pheromone of good ants [37]. In this paper, the elite individuals of each generation are preserved by the percentage of elitist_coef, which are merged with the next generation to enhance the number of excellent individuals and enhance the pheromone on the path of the excellent solution.

    The algorithm optimization process considers the busy degree of vehicles, gives priority to the less busy vehicles for distribution, and improves the local search ability of the algorithm through elite strategy. The vehicle selects the original warehouse for replenishment and uses the dictionary order method in multi-objective processing. The flow chart of the improved algorithm is shown in Figure 6.

    Figure 6.  Flow chart of improved ant colony optimization.

    Using the simulation data, the validity of the model and the improved ACO are verified. In the set of this simulation experiment, enterprise M is a professional logistics company that provides service to stores and supermarkets. In this paper, 41 stores served by enterprise M are selected for data collection. The location coordinates of 41 stores and 3 warehouses, the time window requirements of customers, and the demand of each customer point are shown in Table 1.

    Table 1.  Information table of each node.
    No. Coordinate (X) Coordinate (Y) Demand (tons) Time Window Tolerance time
    1 113.696 36.588 1.524 [1,9] [0, 14]
    2 114.484 36.59 1.257 [1,9] [0, 14]
    3 114.517 36.767 1.139 [1,9] [0, 14]
    4 114.805 36.563 1.348 [1,9] [0, 14]
    5 114.545 36.626 1.227 [1,9] [0, 14]
    6 115.165 36.29 0.489 [1,9] [0, 14]
    7 114.468 36.636 0.536 [1,9] [0, 14]
    8 114.378 36.388 0.078 [1,9] [0, 14]
    9 114.506 36.624 1.032 [1,9] [0, 14]
    10 114.499 36.602 0.125 [1,9] [0, 14]
    11 114.38 36.365 1.703 [1,9] [0, 14]
    12 114.699 36.452 1.061 [1,9] [0, 14]
    13 114.527 36.641 0.221 [1,9] [0, 14]
    14 114.241 36.424 1.171 [1,9] [0, 14]
    15 114.509 36.613 0.481 [1,9] [0, 14]
    16 114.209 36.694 1.107 [1,9] [0, 14]
    17 114.555 36.614 1.151 [1,9] [0, 14]
    18 114.524 36.631 1.261 [1,9] [0, 14]
    19 114.506 36.633 0.431 [1,9] [0, 14]
    20 114.478 36.614 0.511 [1,9] [0, 14]
    21 114.534 36.613 0.06 [1,9] [0, 14]
    22 114.489 36.598 0.055 [1,9] [0, 14]
    23 114.528 36.606 0.491 [1,9] [0, 14]
    24 114.552 36.615 0.611 [1,9] [0, 14]
    25 114.502 36.369 1.171 [1,9] [0, 14]
    26 114.932 36.362 1.101 [1,9] [0, 14]
    27 114.803 36.556 1.23 [1,9] [0, 14]
    28 114.22 36.479 1.57 [1,9] [0, 14]
    29 114.68 36.45 1.61 [1,9] [0, 14]
    30 114.625 36.334 0.4 [1,9] [0, 14]
    31 114.538 36.627 0.921 [1,9] [0, 14]
    32 114.506 36.594 0.09 [1,9] [0, 14]
    33 114.454 36.622 1.421 [1,9] [0, 14]
    34 114.938 36.367 1.33 [1,9] [0, 14]
    35 114.949 36.359 0.431 [1,9] [0, 14]
    36 114.542 36.563 1.794 [1,9] [0, 14]
    37 114.486 36.602 0.068 [1,9] [0, 14]
    38 114.542 36.612 0.037 [1,9] [0, 14]
    39 114.511 36.605 1.597 [1,9] [0, 14]
    40 114.731 36.541 1.344 [1,9] [0, 14]
    41 114.616 36.347 1.342 [1,9] [0, 14]
    42 114.549 36.599
    43 114.205 36.698
    44 114.512 36.582

     | Show Table
    DownLoad: CSV

    As the operating hours of convenience stores are 8:00–22:00, in order to facilitate the management of goods, the ideal time to receive goods is 9:00–17:00. During the rest of the time, goods can be received at 8:00–9:00 and 17:00–22:00 with dissatisfaction generated due to staff shortage. Dissatisfaction of receiving goods before 8:00 and after 22:00 is 1. The time when the vehicle starts to deliver is set as 0, and the delivery time and time window are accumulated in turn.

    The distribution process is set up in the distribution center with three distribution vehicles, each with a maximum load of 5 tons, and the proposed algorithm is implemented in MATLAB2017a. The parameters used in this paper are set according to the data of China Emission Accounts and Datasets and reference [38], as shown in the following Table 2.

    Table 2.  Description of relevant parameters.
    Parameter Value Parameter Value
    Unit fuel price pfuel 6.99 RMB/liter Weight of cost ω3 0.7
    Penalty factor per unit time if the vehicle arrives early ce 100 RMB/h Number of ants ant_n 98
    Penalty factor per unit time if the vehicle arrives late cl 25 RMB/h Probability of selecting the distribution vehicle with the lowest load capacity minloadcar_p 0.45
    Service hour per customer t0 0.25 h Probability of selecting a client according to the maximum transfer probability maxshift_p 0.003
    Time required for distribution vehicle to be loaded at a warehouse t1 0.5 h Importance degree of pheromone α 1.4
    Travel speed of distribution vehicles 55 Km/h Visibility importance degree β 1.7
    Combustion conversion factorF 2.3 Pheromone evaporation rate ρ 0.2
    Weight of carbon emissions ω1 0.2 Maximum number of iterations 200
    Weight of customer dissatisfaction ω2 0.1 Constant Q 0.5
    Percentage of the number of ants saved by the elite 2%

     | Show Table
    DownLoad: CSV

    Compare the improved ant colony algorithm and genetic algorithm with the basic ant colony algorithm. Both of them use joint distribution method, and there are 3 distribution vehicles. Take the best results of the three algorithms as the optimal solution for comparison. Table 3 displays the comparison of the number of replenishment, full load rate, cost, dissatisfaction and carbon emissions. As shown in Table 3, the improved ant colony algorithm has 95.49, 92.4 and 86.42% full load rates, which are better than the other two algorithms, respectively. It can also be seen that the cost and carbon emission of the improved ant colony algorithm is the smallest due to the highest vehicle loading efficiency of the improved ant colony algorithm, which also improves customer satisfaction.

    Table 3.  Comparison of the effectiveness of the improved ant colony algorithm.
    Comparison of algorithm results Number of replenishment Full load rate Cost Dissatisfaction Carbon Emissions
    Improved ant colony algorithm 6 95.49% 92.4% 86.42% 583.91 1.948 45.15
    Basic ant colony algorithm 6 89.81% 93.71% 89.98% 636.71 2.224 45.75
    Genetic algorithm 6 93.97% 92.4% 84.93% 692.60 2.08 47.29

     | Show Table
    DownLoad: CSV

    Besides, the comparison reveals that the improved ant colony algorithm outperforms the basic ant colony algorithm and genetic algorithm in all metrics.

    After the algorithm enters the stagnation phase, the first ranked ant route is selected as the satisfactory solution, and since the model deals with the actual problem, the satisfactory solution can be used as an effective solution for this problem. The solutions of this multi-objective optimization are shown in Table 4, while the carbon emission and customer satisfaction are shown in Table 5. According to the results shown in Table 4, the full load rate of all vehicles exceeds 85%, and the full load rate of vehicle 1 reaches 95.49%. The lowest cost obtained is 583.91 Yuan, which solves the cost loss problem to a greater extent; the lowest customer satisfaction is 1.948, and the vehicle full load rate improves the efficiency of cargo arrival, thus improving customer satisfaction; the lowest carbon emission is 45.15 Kg, and the increase of vehicle full load rate also leads to the reduction of the number of times the vehicle is used, thus suppressing the increase of carbon emission.

    Table 4.  Solutions for multi-objective optimization.
    No. of vehicle Distribution route Carbon emissions (kg) Fuel cost Number of replenishment Full load rate Cost
    (Yuan)
    Dissatisfaction Carbon Emissions
    (kg)
    1 42→1→36→34→
    26→42→29→40→4→27→2→42→23→10→37→20→7→22→42
    22.7031 59.23 2 95.49% 583.91 1.948 45.15
    2 43→6→35→41→
    12→43→30→11→8→9→39→16→43
    11.7817 30.73 1 92.4%
    3 44→28→14→25→17→44→3→33→
    32→15→13→44→18→19→31→5→
    24→44→21→38→
    44
    10.6698 27.83 3 86.42%

     | Show Table
    DownLoad: CSV
    Table 5.  Delivery time and customer dissatisfaction.
    No. Of customer Customer
    dissatisfaction
    Carbon Emissions No. Of customer Customer dissatisfaction Carbon Emissions
    1 0.14693 0.51499 22 0 1.8113
    2 0
    1.7562 23 0 1.7552
    3 0 0.62484 24 0 0.82553
    4 0 1.5623 25 0.36476 0.38062
    5 0 0.8287 26 0 1.2858
    6 0 0.6159 27 0 1.5627
    7 0 1.8047 28 0.69037 0.1871
    8 0 1.256 29 0 1.4612
    9 0 1.4425 30 0 1.1057
    10 0 1.7584 31 0 0.81908
    11 0 1.2866 32 0 0.73796
    12 0 1.0383 33 0 0.72384
    13 0 0.7709 34 0 1.2874
    14 0.63149 0.2208 35 0 0.74885
    15 0 0.75562 36 0 1.0318
    16 0 1.645 37 0 1.7638
    17 0.11409 0.53059 38 0 0.8306
    18 0 0.79463 39 0 1.4719
    19 0 0.79126 40 0 1.5157
    20 0 1.7894 41 0 0.96403
    21 0 0.8263

     | Show Table
    DownLoad: CSV

    The routes of three delivery vehicles are obtained according to calculation results, which are shown in Figures 79, and the overall routes of three vehicles is presented in Figure 10.

    Figure 7.  Distribution route of vehicle 1.
    Figure 8.  Distribution route of vehicle 2.
    Figure 9.  Distribution route of vehicle 3.
    Figure 10.  Distribution routes of vehicle 1, vehicle 2, and vehicle 3.

    From delivery routes presented in above figures, it is obvious to find that vehicle 1 is mainly replenished in distribution center 42 during the distribution process, with 2 replenishments. Vehicle 2 is mainly replenished in distribution center 43 with 1 replenishments. And vehicle 3 is mainly replenished in distribution center 44, with 3 replenishment. The cross-regional distribution characteristics of vehicles 1, 2 and 3 are clearly displayed in the distribution process.

    The full load rate of the first vehicle is 95.49%, its carbon emission is 22.7031 KG, and the fuel cost is 59.23 Yuan. It has been replenished twice. The full load rate of the second vehicle is 92.4%, its carbon emission is 11.7817 KG, and its fuel cost is 30.73 Yuan. It is replenished once. The full load rate of the third vehicle is 86.42%, its carbon emission is 10.6698 KG, and its fuel cost is 27.83 Yuan. It has been replenished three times. Five of the 41 stores failed to deliver at the specified time, and their total dissatisfaction was 1.948.

    This section discusses the industrial benefits for industry managers and are as follows:

    The implementation of a reliable logistics and transportation system is a necessary condition for environmental protection. At the same time, a reliable logistics transportation system can also respond more quickly to the transportation route, reduce the transportation cost, and thus improve the core competitiveness of enterprises.

    Based on above results analysis, it can be concluded that the solution of proposed algorithm achieves lower distribution cost with higher customer satisfaction when considering fuel consumption cost, penalty cost, and other factors in route planning. Considering the carbon emissions of trucks, it should pay attention to the full load ratio of trucks. The higher the carrying capacity, the more carbon emissions, the lower the carrying capacity, and the higher the cost. Consequently, VRP optimization with considerations of carbon emissions calculated by the actual weight of the vehicle and the driving mileage helps to realize green logistics.

    For enterprises, multi-center joint distribution considering cost, customer dissatisfaction and carbon emissions can satisfy diversified customer needs, meeting requirement of multi-species, small-lot, and individualized orders and improve customer experience. This implementation makes contributions of reducing enterprise operation cost, improving logistics and distribution efficiency, as well as realizing good social value.

    To reduce damage to the environment during operation, logistics enterprises pay more attention to green logistics, thus, energy consumption and carbon emission as well as customer dissatisfaction are core issue to solve in distribution system. In this paper, a multi-objective planning model considering distribution cost, customer dissatisfaction and carbon emission is constructed, and an intelligent optimization algorithm is designed to solve the problem. Based on the characteristics of slow convergence and easy premature of traditional ant colony algorithm, elite strategy and saving strategy are used to improve the ant colony algorithm. The case study results show that the improved ACO has better performance compared to the basic ACO, which can effectively solve this multi-objective optimization problem and can be applied to similar multi-objective optimization. At the same time, the joint distribution of multiple distribution centers can break the barrier that resources cannot be shared between regions, realize the purpose of reducing costs and increasing efficiency of enterprises. In the process of multi-center joint distribution, considering distribution cost, customer dissatisfaction and carbon emissions can realize both economic and social benefits, which makes it be more in line with the requirements of green GDP regulations.

    The future development direction of this paper is as follows:

    Multi objective issues should be considered more comprehensively. This paper considers distribution cost, customer dissatisfaction and carbon emissions, but only considers distance, weight, full load ratio and other factors when quantifying these factors, lacking comprehensive consideration of more factors.

    The solving efficiency of the improved ant colony algorithm should also be improved. The problem studied in this paper is NP hard problem, which is difficult to solve. When solving this multi-objective optimization problem, the ant colony algorithm is improved to improve the solution efficiency of the algorithm, but the search ability of the ant colony still needs to be enhanced.

    Therefore, the next development direction is as follows:

    The driving speed of the vehicle shall be considered when calculating carbon emissions.

    When calculating customer dissatisfaction, consider the problem of traffic congestion when driving.

    Vehicles are required to go to the nearest warehouse for replenishment.

    The pheromone updating function in ant colony algorithm needs to be improved.

    The proposed algorithm has no verification and exploration of multiple scenarios or multiple data sets.

    This paper was supported by the Social Science Grand Research of the Hebei Education Department (ZD202105). The authors also gratefully acknowledge the helpful comments and suggestions of the reviewers, which have improved the presentation.

    The authors declare that there are no conflicts of interest.

    The research data used to support the findings of this study are included within the article, the data used are shown in Tables 1 and 2.



    [1] Andreasen C, Stryhn H, Streibig JC (1996) Decline of the flora of Dutch arable fields. J Appl Ecol 33: 619–626. https://doi.org/10.2307/2404990 doi: 10.2307/2404990
    [2] Sutcliffe OL, Kay QON (2000) Changes in the arable flora of central southern England since the 1960s. Biol Conserv 93: 1–8. https://doi.org/10.1016/S0006-3207(99)00119-6 doi: 10.1016/S0006-3207(99)00119-6
    [3] Fried G, Petit S, Dessaint F, et al. (2009) Arable weed decline in Northern France: Crop edges as refugia for weed conservation? Biol Conserv 142: 238–243. https://doi.org/10.1016/j.biocon.2008.09.029 doi: 10.1016/j.biocon.2008.09.029
    [4] Wilson S, Mitchell GW, Pasher J, et al. (2017) Influence of crop type, heterogeneity and woody structure on avian biodiversity in agricultural landscapes. Ecol Indic 83: 218–226. http://dx.doi.org/10.1016/j.ecolind.2017.07.059 doi: 10.1016/j.ecolind.2017.07.059
    [5] Lomba A, da Costa JF, Ramil-Rego P, et al. (2022) Assessing the link between farming systems and biodiversity in agricultural landscapes: Insights from Galicia (Spain). J Environ Manage 317: 115335. https://doi.org/10.1016/j.jenvman.2022.115335 doi: 10.1016/j.jenvman.2022.115335
    [6] Thomine E, Mumford J, Rusch A, et al. (2022) Using crop diversity to lower pesticide use: Socio-ecological approaches. Sci Total Environ 804: 150156. https://doi.org/10.1016/j.scitotenv.2021.150156 doi: 10.1016/j.scitotenv.2021.150156
    [7] Singh A (2017) Managing the environmental problems of irrigated agriculture through the appraisal of groundwater recharge. Ecol Indic 92: 388–393. https://doi.org/10.1016/j.ecolind.2017.11.065 doi: 10.1016/j.ecolind.2017.11.065
    [8] Marshall EJP, Brown VK, Boatman ND, et al. (2003) The role of weeds in supporting biological diversity within crop fields. Weed Res 43: 7–89. https://doi.org/10.1046/j.1365-3180.2003.00326.x doi: 10.1046/j.1365-3180.2003.00326.x
    [9] Nagy GG, Ladányi M, Arany I, et al. (2017) Birds and plants: Comparing biodiversity indicators in eight lowland agricultural mosaic landscapes in Hungary. Ecol Indic 73: 566–573. https://doi.org/10.1016/j.ecolind.2016.09.053 doi: 10.1016/j.ecolind.2016.09.053
    [10] Yvoz S, Cordeau S, Ploteau A, et al. (2021a) A framework to estimate the contribution of weeds to the delivery of ecosystem (dis) services in agricultural landscapes. Ecol Indic 132: 108321. https://doi.org/10.1016/j.ecolind.2021.108321 doi: 10.1016/j.ecolind.2021.108321
    [11] Chancellor RJ (1985) Changes in the weed flora of an arable cultivated field for 20 years. J Appl Ecol 22: 491–501. https://doi.org/10.2307/2403180 doi: 10.2307/2403180
    [12] Hald AB (1999) The impact of changing the season in which cereals are sown on the diversity of the weed flora in rotational fields in Denmark. J Appl Ecol 36: 24–32. https://doi.org/10.1046/j.1365-2664.1999.00364.x doi: 10.1046/j.1365-2664.1999.00364.x
    [13] Krähmer H, Andreasen C, Economou-Antonaka G, et al. (2020) Weed surveys and weed mapping in Europe: State of the art and future tasks. Crop Prot 129: 105010. https://doi.org/10.1016/j.cropro.2019.105010 doi: 10.1016/j.cropro.2019.105010
    [14] Preston CD, Pearman DA, Dines T (2002) New Atlas of the British and Irish Flora. Oxford University Press, Oxford.
    [15] Hyvönen T, Ketoja E, Salonen J, et al. (2003) Weed species diversity and community composition in organic and conventional cropping of spring cereals. Agric, Ecosyst Environ 97: 131–149. https://doi.org/10.1016/S0167-8809(03)00117-8 doi: 10.1016/S0167-8809(03)00117-8
    [16] Wilson PJ, Aebischer NJ (1995) The distribution of dicotyledonous arable weeds in relation to distance from the field edge. J Appl Ecol 32: 295–310. https://doi.org/10.2307/2405097 doi: 10.2307/2405097
    [17] Kleijn D, van der Voort LAC (1997) Conservation headlands for rare arable weeds: The effects of fertiliser application and light penetration on plant growth. Biol Conserv 81: 57–67. https://doi.org/10.1016/S0006-3207(96)00153-X doi: 10.1016/S0006-3207(96)00153-X
    [18] Solé-Senan XO, Juárez-Escario A, Conesa JA, et al. (2014) Plant diversity in Mediterranean cereal fields: Unraveling the effect of landscape complexity on rare arable plants. Agric, Ecosyst Environ 185: 221–230. https://doi.org/10.1007/s10661-019-8042-7 doi: 10.1007/s10661-019-8042-7
    [19] Yvoz S, Cordeau S, Zuccolo C, et al. (2020) Crop type and within-field location as sources of intraspecific variations in the phenology and the production of floral and fruit resources by weeds. Agric, Ecosyst Environ 302: 107082. https://doi.org/10.1016/j.agee.2020.107082 doi: 10.1016/j.agee.2020.107082
    [20] Critchley CNR, Allen DS, Fowbert JA, et al. (2004) Habitat establishment on arable land: assessment of an agri-environment scheme in England, UK. Biol Conserv 119: 429–442. https://doi.org/10.1016/j.biocon.2004.01.004 doi: 10.1016/j.biocon.2004.01.004
    [21] Sotherton NW (1991) Conservation headlands: A practical combination of intensive cereal farming and conservation. In: Firbank LG, Carter N, Darbyshire JF, et al. (Eds.), The Ecology of Temperate Cereal Fields, Blackwell, London, 373–397.
    [22] Yvoz S, Petit S, Cadet E, et al. (2021b) Taxonomic and functional characteristics of field edge weed communities along a gradient of crop management intensity. Basic Appl Ecol 57: 14–27. https://doi.org/10.1016/j.baae.2021.10.001 doi: 10.1016/j.baae.2021.10.001
    [23] Firbank LG, Smart SM, Crabb J, et al. (2003) Agronomic and ecological costs and benefits of set-aside in England, Agric, Ecosyst Environ 95: 73–85. https://doi.org/10.1016/S0167-8809(02)00169-X doi: 10.1016/S0167-8809(02)00169-X
    [24] Kleijn D, Joenje W, Le Coeur D, et al. (1998) Similarities in vegetation development of newly established herbaceous strips along contrasting European field boundaries. Agric, Ecosyst Environ 68: 13–26. https://doi.org/10.1016/S0167-8809(97)00098-4 doi: 10.1016/S0167-8809(97)00098-4
    [25] Critchley CNR, Fowbert JA (2000) Development of vegetation on set-aside land for up to nine years from a national perspective. Agric, Ecosyst Environ 79: 159–173. https://doi.org/10.1016/S0167-8809(99)00155-3 doi: 10.1016/S0167-8809(99)00155-3
    [26] Ryspekov T, Jandák J, Balkozha M, et al. (2021) Vegetation of abandoned fields on soil types of kastanozems in Northern Kazakhstan (Kostanay Region). J Ecol Eng. 22: 176–184. https://doi.org/10.12911/22998993/142188 doi: 10.12911/22998993/142188
    [27] Weibull AC, Ostman O, Granqvist A (2003) Richness in agroecosystems: The effect of landscape, habitat and farm management. Biodivers Conserv 12: 1335–1355. https://www.jstor.org/stable/3838497
    [28] Burel F, Butet A, Delettre YR, et al. (2004) Differential response of selected taxa to landscape context and agricultural intensification. Landscape Urban Plann 67: 195–204. http://dx.doi.org/10.1016/S0169-2046(03)00039-2 doi: 10.1016/S0169-2046(03)00039-2
    [29] Bianchi FJJA, Booijand CJH, Tscharntke T (2006) Sustainable pest regulation in agricultural landscapes: a review on landscape composition, biodiversity and natural pest control. Proc R Soc B 273: 1715–1727. https://doi.org/10.1098/rspb.2006.3530 doi: 10.1098/rspb.2006.3530
    [30] Geiger F, Waeckers FL, Bianchi FJJA (2009) Hibernation of predatory arthropods in semi-natural habitats. Biocontrol 54: 529–535. https://doi.org/10.1007/s10526-008-9206-5 doi: 10.1007/s10526-008-9206-5
    [31] Rusch A, Valantin-Morison M, Sarthou JP, et al. (2010) Biological control of insect pests in agroecosystems: Effects of crop management, farming systems, and semi-natural habitats at the landscape scale: A review. Adv Agric 109: 19–259. https://doi.org/10.1016/B978-0-12-385040-9.00006-2 doi: 10.1016/B978-0-12-385040-9.00006-2
    [32] Storkey J (2006) Afunctional groupapproachto themanagement of UK arable weeds to support biological diversity. Weed Res 46: 513–522. https://doi.org/10.1111/j.1365-3180.2006.00528.x doi: 10.1111/j.1365-3180.2006.00528.x
    [33] Bastiaans L, Kropff MJ, Goudriaan J, et al. (2000) Design of weed management systems with a reduced reliance on herbicides poses new challenges and prerequisites for modeling crop—weed interactions. Field Crops Res 67: 161–179. https://doi.org/10.1016/S0378-4290(00)00091-5 doi: 10.1016/S0378-4290(00)00091-5
    [34] Marshall EJP, Brown VK, Boatman ND, et al. (2003) The role of weeds in supporting biological diversity within crop fields. Weed Res 43: 77–89. https://doi.org/10.1046/j.1365-3180.2003.00326.x doi: 10.1046/j.1365-3180.2003.00326.x
    [35] Petit S, Boursault A, Le Guilloux M, et al. (2011) Weeds in agricultural landscapes: A review. Agron Sustain Dev 31: 309–317. https://doi.org/10.1051/agro/2010020 doi: 10.1051/agro/2010020
    [36] Gaba S, Fried G, Kazakou E, et al. (2014) Agroecological weed control using a functional approach: A review of cropping systems diversity. Agron Sustainable Dev 34: 103–119. https://doi.org/10.1007/s13593-013-0166-5 doi: 10.1007/s13593-013-0166-5
    [37] Schumacher M, Dieterich M, Gerhards R (2020) Effects of weed biodiversity on the ecosystem service of weed seed predation along a farming intensity gradient. Global Ecol Conserv 24: e01316. https://doi.org/10.1016/j.gecco.2020.e01316 doi: 10.1016/j.gecco.2020.e01316
    [38] Dieleman JA, Mortensen DA (1999) Characterizing the spatial pattern of Abutilon theophrasti seedling patches. Weed Res 39: 455–467. https://doi.org/10.1046/j.1365-3180.1999.00160.x doi: 10.1046/j.1365-3180.1999.00160.x
    [39] Marshall EJP (1988) Field-scale estimates of grass weed populations in arable land. Weed Res 28: 191–198. https://doi.org/10.1111/j.1365-3180.1988.tb01606.x doi: 10.1111/j.1365-3180.1988.tb01606.x
    [40] Rew LJ, Cussans GW, Mugglestone MA, et al. (1996) A technique for mapping the spatial distribution of Elymus repens, with estimates of the potential reduction in herbicide usage from patch spraying. Weed Res 36: 283–292. https://doi.org/10.1111/j.1365-3180.1996.tb01658.x doi: 10.1111/j.1365-3180.1996.tb01658.x
    [41] Winkler J, Lukas V, Smutný V (2015) Analysis of the heterogeneity of weed infestation in cereal stands. Acta Univ Agric Silvic Mendelianae Brun 63: 153–161. http://dx.doi.org/10.11118/actaun201563010153 doi: 10.11118/actaun201563010153
    [42] Bourgeois B, Gaba S, Plumejeaud C, et al. (2020) Weed diversity is driven by complex interplay between multi-scale dispersal and local filtering. Proc R Soc B 287: 20201118. https://doi.org/10.1098/rspb.2020.1118 doi: 10.1098/rspb.2020.1118
    [43] Hussain RI, Brandl M, Maas B, et al. (2021) Re-established grasslands on farmland promote pollinators more than predators. Agric, Ecosyst Environ 319: 107543. https://doi.org/10.1016/j.agee.2021.107543 doi: 10.1016/j.agee.2021.107543
    [44] Šálek M, Zámečník V (2020) Historic overview and conservation perspectives of the Czech grey partridge (Perdix perdix) population. In: The Changing Status of Arable Habitats in Europe, Springer, Cham, 227–243. https://doi.org/10.1007/978-3-030-59875-4_15
    [45] Hejcmanová-Nežerková P, Hejcman M (2006) A canonical correspondence analysis (CCA) of the vegetation–environment relationships in Sudanese savannah, Senegal. S Afr J Bot 72: 256–262. https://doi.org/10.1016/j.sajb.2005.09.002 doi: 10.1016/j.sajb.2005.09.002
    [46] Ter Braak CJF, Šmilauer P (2012) Canoco reference manual and user's guide: Software for ordination (version 5.0). Ithaca: Microcomputer Power, 2012.
    [47] Chytrý M, Danihelka J, Kaplan Z, et al. (2021) Pladias database of the Czech flora and vegetation. Preslia 93: 1–87. https://doi.org/10.23855/preslia.2021.001 doi: 10.23855/preslia.2021.001
    [48] Durka W (2002) Blüten- und Reproduktionsbiologie. In: Klotz S, Kühn I, Durka W (Eds), BIOLFLOR—Eine Datenbank mit biologisch-ökologischen Merkmalen zur Flora von Deutschland, Schriftenreihe für Vegetationskunde, 38: 133–175.
    [49] Sádlo J, Chytrý M, Pergl J, et al. (2018) Plant dispersal strategies: A new classification based on themultiple dispersal modes of individual species. Preslia 90: 1–22. https://doi.org/10.23855/preslia.2018.001 doi: 10.23855/preslia.2018.001
    [50] Šálek M, Hula V, Kipson M, et al. (2018) Bringing diversity back to agriculture: Smaller fields and non-crop elements enhance biodiversity in intensively managed arable farmlands. Ecol Indic 90: 65–73. https://doi.org/10.1016/j.ecolind.2018.03.001 doi: 10.1016/j.ecolind.2018.03.001
    [51] Otieno NE, Jacobs SM, Pryke JS (2022) Small-scale traditional maize farming fosters greater arthropod diversity value than conventional maize farming. J Insect Conserv 26: 477–489. https://doi.org/10.1007/s10841-021-00330-x doi: 10.1007/s10841-021-00330-x
    [52] Tarjuelo R, Benítez-López A, Casas F, et al. (2020) Living in seasonally dynamic farmland: The role of natural and semi-natural habitats in the movements and habitat selection of a declining bird. Biol Conserv 251: 108794. https://doi.org/10.1016/j.biocon.2020.108794 doi: 10.1016/j.biocon.2020.108794
    [53] Jachuła J, Denisow B, Wrzesień M, et al. (2022) The need for weeds: Man-made, non-cropped habitats complement crops and natural habitats in providing honey bees and bumble bees with pollen resources. Sci Total Environ 2022: 156551. https://doi.org/10.1016/j.scitotenv.2022.156551 doi: 10.1016/j.scitotenv.2022.156551
    [54] Ekwealor KU, Echereme CB, Ofobeze TN, et al. (2019) Economic importance of weeds: A review. Asian J Plant Sci 3: 1–11. https://doi.org/10.9734/aprj/2019/v3i230063 doi: 10.9734/aprj/2019/v3i230063
    [55] Tissier ML, Marchandeau S, Habold C, et al. (2019) Weeds as a predominant food source: A review of the diet of common hamsters Cricetus cricetus in farmlands and urban habitats. Mamm Rev 49: 152–170. https://doi.org/10.1111/mam.12149 doi: 10.1111/mam.12149
    [56] Ouvrard P, Jacquemart AL (2018) Agri‐environment schemes targeting farmland bird populations also provide food for pollinating insects. Agric For Entomol 20: 558–574. https://doi.org/10.1111/afe.12289 doi: 10.1111/afe.12289
    [57] Crochard L, Julliard R, Gaba S, et al. (2022) Weeds from non-flowering crops as potential contributors to oilseed rape pollination. Agric, Ecosyst Environ 336: 108026. https://doi.org/10.1016/j.agee.2022.108026 doi: 10.1016/j.agee.2022.108026
    [58] Bretagnolle V, Gaba S (2015) Weeds for bees? A review. Agron Sustainable Dev 35: 891–909. https://doi.org/10.1007/s13593-015-0302-5 doi: 10.1007/s13593-015-0302-5
    [59] Gaba S, Cheviron N, Perrot T, et al. (2020) Weeds enhance multifunctionality in arable lands in south-west of France. Front Sustainable Food Syst 4: 71. https://doi.org/10.3389/fsufs.2020.0007 doi: 10.3389/fsufs.2020.0007
    [60] Kleiman B, Koptur S, Jayachandran K (2021) Beneficial interactions of weeds and pollinators to improve crop production. J Res Weed Sci 4: 151–164.
    [61] Marshall EJP (1989) Distribution patterns of plants associated with arable field edges. J Appl Ecol 26: 247–257. https://doi.org/10.2307/2403665 doi: 10.2307/2403665
    [62] Alignier A, Uroy L, Aviron S (2020) The role of hedgerows in supporting biodiversity and other ecosystem services in intensively managed agricultural landscapes. In Reconciling agricultural production with biodiversity conservation, Burleigh Dodds Science Publishing, 177–204.
    [63] Martín‐Vélez V, van Leeuwen CHM, Sánchez I, et al. (2021) Spatial patterns of weed dispersal by wintering gulls within and beyond an agricultural landscape. J Ecol 109: 1947–1958. https://doi.org/10.1111/1365-2745.13619 doi: 10.1111/1365-2745.13619
    [64] Marshall EJP, Smith BD (1987) Field margin flora and fauna: Interaction with agriculture. In: Way JM, Greig-Smith PW (Eds.), Field Margins, British Crop Protection Council, UK, 23–33.
    [65] Mante J, Gerowitt B (2009) Learning from farmers' needs: Identifying obstacles to the successful implementation of field margin measures in intensive arable regions. Landscape Urban Plann 93: 229–237. https://doi.org/10.1016/j.landurbplan.2009.07.010 doi: 10.1016/j.landurbplan.2009.07.010
    [66] Cordeau S, Reboud X, Chauvel B (2011) Farmers' fears and agro-economic evaluation of sown grass strips in France. Agron Sustain Dev 31: 463–473. https://doi.org/10.1007/s13593-011-0004-6 doi: 10.1007/s13593-011-0004-6
    [67] Pellegrini E, Buccheri M, Martini F, et al. (2021) Agricultural land use curbs exotic invasion but sustains native plant diversity at intermediate levels. Sci Rep 11: 1–10. https://doi.org/10.1038/s41598-021-87806-7 doi: 10.1038/s41598-021-87806-7
  • This article has been cited by:

    1. Xiaolin Wang, Liyi Zhan, Yong Zhang, Teng Fei, Ming-Lang Tseng, Environmental cold chain distribution center location model in the semiconductor supply chain: A hybrid arithmetic whale optimization algorithm, 2024, 187, 03608352, 109773, 10.1016/j.cie.2023.109773
    2. Binghui Qie, Xun Weng, Zhiwei Sun, Minyu Jin, Runfeng Yu, Location-routing and cost-sharing models under joint distribution, 2024, 27, 1386-7857, 5879, 10.1007/s10586-024-04282-0
    3. Huaixia Shi, Qinglei Zhang, Jiyun Qin, Cold Chain Logistics and Joint Distribution: A Review of Fresh Logistics Modes, 2024, 12, 2079-8954, 264, 10.3390/systems12070264
    4. Bin Yue, Junxu Ma, Jinfa Shi, Jie Yang, A Deep Reinforcement Learning-Based Adaptive Search for Solving Time-Dependent Green Vehicle Routing Problem, 2024, 12, 2169-3536, 33400, 10.1109/ACCESS.2024.3369474
    5. Changlu Zhang, Liqian Tang, Jian Zhang, Liming Gou, Optimizing Distribution Routes for Chain Supermarket Considering Carbon Emission Cost, 2023, 11, 2227-7390, 2734, 10.3390/math11122734
    6. Xianlong Ge, Yonghong Liang, Yuanzhi Jin, Chunbing Song, Proactive dynamic vehicle routing problems considering cooperation services for the store-depot-integrated retailer, 2023, 20, 1551-0018, 18030, 10.3934/mbe.2023801
  • Reader Comments
  • © 2022 the Author(s), licensee AIMS Press. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0)
通讯作者: 陈斌, bchen63@163.com
  • 1. 

    沈阳化工大学材料科学与工程学院 沈阳 110142

  1. 本站搜索
  2. 百度学术搜索
  3. 万方数据库搜索
  4. CNKI搜索

Metrics

Article views(2104) PDF downloads(153) Cited by(2)

Figures and Tables

Figures(5)  /  Tables(3)

/

DownLoad:  Full-Size Img  PowerPoint
Return
Return

Catalog