Next Article in Journal
Applications of a Novel Method of Ecosystem Services Assessment into Local Policy Making in the River Blackwater Estuary, Ireland
Previous Article in Journal
A Systematic Review of Air Quality Sensors, Guidelines, and Measurement Studies for Indoor Air Quality Management
Previous Article in Special Issue
Glyphosate Resistance in Sonchus oleraceus and Alternative Herbicide Options for Its Control in Southeast Australia
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Economic Costs of Sharing the Harvester in the Control of an Invasive Weed

1
Unidad de Protección Vegetal, Centro de Investigación y Tecnología Agroalimentaria de Aragón (CITA), Avda. Montañana 930, 50059 Zaragoza, Spain
2
Instituto Agroalimentario de Aragón—IA2 (CITA-Universidad de Zaragoza), 50013 Zaragoza, Spain
3
Charles H. Dyson School of Applied Economics & Management, Cornell University, Ithaca, NY 14853, USA
4
Department of Economic Analysis, University of Zaragoza, Gran Vía 2–4, 50004 Zaragoza, Spain
*
Author to whom correspondence should be addressed.
Sustainability 2020, 12(21), 9046; https://doi.org/10.3390/su12219046
Submission received: 17 September 2020 / Revised: 23 October 2020 / Accepted: 27 October 2020 / Published: 30 October 2020
(This article belongs to the Special Issue Weed Management in Sustainable Agricultural Production)

Abstract

:
Spatial externalities, such as the sharing of harvesting equipment by many farmers, have an impact on the control of invasive species in the agricultural environment. In these cases, the regulator must design a set of measures to promote coordinated control by affected parties. We aim to analyze the determinants of private versus collective control efforts in the case of a particular invasive species (teosinte) occurring as a weed in corn fields throughout North-Eastern Spain. Using a simple discrete space-dynamic framework, we model the effect of the decisions made by the farmer of an infested plot on a noninfested plot, with the harvester being the only potential pathway for the invader to spread and assuming a one-way invasion. The results reveal that failure to adopt optimal cooperative strategies causes losses to other plots if they become infested amounting to an annual average of EUR 322/ha, when the infestation is low, and EUR 364/ha, when it is high. Results suggest that cleaning the harvester, a measure currently recommended by the regulatory agency in low-infestation cases but that does not guarantee that the machine is completely clean, is not socially optimal if monocropping practices are permitted in the region.

1. Introduction

The emergence of a new invasive species that behaves as a weed in crop fields is a major challenge for the agents involved in controlling it (farmers and regulators). Rapidly understanding how the new species adapts to local environmental conditions is vital in the early stages of detection to design the most appropriate prevention, containment and eradication strategies. Identifying possible dispersal pathways—which are normally affected by control externalities—among neighboring plots is also of great importance.
Evaluating these externalities is crucial in the practical framework so that regulators can adequately guide farmers in their control decisions. The primary reason why these externalities exist in the case of agricultural fields is that the individuals (or stakeholders) involved in control base their decisions on their own farms [1], in other words, on a subset of the total area at risk of invasion. This private perspective ignores the link between the control efforts made by neighboring fields, while spatial externalities are likely when control efforts affect the spread of an invasive species across the landscape [2,3].
Control decisions made by farmers are also often based on measures after the invader has already become established in the environment, thus neglecting prevention measures in what is known as “myopic” behavior [4]. In these cases, the role of the regulator, who approaches the problem from a collective standpoint, is to design a set of measures that promote coordinated control by the affected parties.
Ultimately, the decisions made by the parties (farmers and regulators) involved in controlling the invader are influenced by an attempt to minimize the damage the invader causes by trying to reduce the probability of infestation through prevention measures, or by investing in control measures once the invader has become established [5]. The fact that every species has its own particular dispersal method is an additional difficulty in adopting preventive measures that must be addressed using a specific procedure.
The problem of evaluating externalities in the control of an invasive species and the available solutions is not new in the literature. In this context, solutions to internalize these externalities, many based on the spirit of [6], are diverse and include the introduction of side payments between producers [7], landowners’ bilateral negotiation [2,8], and voluntary contributions to a cooperative control district [3]. In this vision, a farmer might offer to share eradication or containment costs with a neighbor invaded by a weed to avoid or postpone being invaded.
Advances in species knowledge, optimization techniques and computational tools have driven the design of simulation models that can serve to assist farmers in controlling weeds in their fields. Examples in the literature include weeds in natural environment [9,10,11,12], and in agroecosystems [13,14,15]. Usually, these models are designed as a decision tool for weed control in the fields, making recommendations based on weed densities and available alternative treatments. The inclusion of economic variables in these models makes it possible to calculate and compare expected benefits for the different control strategies and to identify the optimal ones according to the farmer’s objective (e.g., minimum cost, maximum benefit or multicriteria based) [15,16,17].
Recently, some studies consider dynamic and spatial aspects in bioeconomic models of weed control in various agroecosystems [16,18]. These studies address the trade-offs between current weed control and future profitability of the land, and the effect of control effort in a particular plot on the neighboring land. This paper contributes to this growing literature by combining new knowledge on the spatial dispersion of an invasive weed with its impacts on economic costs. To our knowledge, there is no previous work on invasive weed control that incorporates a public cost function associated with monitoring.

Teosinte as an Invasive Plant in the Ebro Valley

In 2014, teosinte (Zea mays subspp.) infestations in corn fields in the Ebro Valley (Spain) prompted researchers to learn about unknown aspects of its biological behavior, such as germination, seed survival capacity [19] and available chemical and manual control methods in affected fields [20].
Since teosinte was first detected in Aragon, the Centro de Sanidad y Certificación Vegetal (CSCV), which is the region’s Plant Protection Service Agency, has surveyed more than 7000 ha each year, and monitored and recorded the number of infested plots, the total affected area and the infestation incidence (low or high infestation levels). Table 1 contains the information from 2014 to 2018, which will be used in this analysis.
The data recorded in that period show an increase in both the number of plots and the total area affected, although the annual growth rate decreased from 2014 to 2017. The area with high infestation and its proportion of the total decreased significantly from the first year to 2017, from 93% in 2014 to 4.2% in 2017. The most recent survey data, from 2018, show that new infested plots have continued to appear, although they differ from previously identified plots [22]. Moreover, there is an increase in the area with high infestation levels, suggesting that spatial dispersal is not being fully controlled despite the CSCV’s dissemination, surveying and monitoring efforts.
This evidence has led to new research aimed to identify the invader’s possible spatial dispersal mechanisms in the area. Technicians considered two main probable dispersal pathways: shared harvesters and stubble sheep grazing in affected areas [23]. Although both have been shown to be potential sources of teosinte dispersal, studies have concluded that harvesters have played a determining role in dispersing teosinte into new plots [22], as it occurs with other weed species, such as Avena sterilis L. [24] and Lolium rigidum [25].
Based on the survey data, the CSCV has established some mandatory phytosanitary measures to control intraplot infestation level and interplot spread of teosinte [26]. These measures comprise a set of cultural controls such as false seedbed technique, manual control, harvester cleaning protocols and rotations without corn. The first two strategies are only recommended for plots with a low infestation level, while harvester cleaning protocols are mandatory when infestation is low and rotations are adopted only when infestation is high. In addition, stubble sheep grazing is prohibited in all infested plots until the infestation has been completely eradicated.
From an economic perspective, the shared use of a harvester by infested and noninfested fields involves the existence of an externality in the control of the invasive weed. In a previous study, [27] obtained the optimal control strategies in a model including weed and seed bank dynamics and the effect of available control methods on those dynamics. An issue that affects the dispersal of this invasive species and still needs to be studied is how a farmer’s control decisions influence neighboring fields, in other words, identifying and analyzing externalities.
Consequently, we specifically focus on the shared use of a harvester in the area identifying how private control decisions on one farm affect the level of infestation on neighboring farms when infestation likelihood is uncertain. With this aim, we construct a bio-economic model comprising dynamic and spatial dispersal dimensions to identify profit-maximizing strategies to tackle the problem of teosinte. In addition, we quantify externalities in controlling the invasive species and we evaluate the possibility of establishing economic compensation mechanisms among farmers through voluntary agreements (cooperation) or enforced by the regulator (no cooperation).

2. Materials and Methods

2.1. Theoretical Model

Our aim is to compare private versus social optimal farmer behavior when negative externalities of teosinte control are included in a dynamic profit-maximization problem and to identify available options to encourage private farmers to undertake a socially optimal level of control. Consequently, we first developed a general theoretical model and then illustrated it with a numerical application using real data obtained in the study area. This allowed us to quantify externalities and to internalize them.
We chose to use a profit-maximization framework because data on market prices and margins of the focal crops were available. One of the advantages of this approach is that it allowed us to carry out sensitivity analyses of the results to these model parameters. In addition, this approach allowed direct estimation of the loss of profit associated with the no-control strategy.
We defined a general discrete space and time model to include the main biological processes of the teosinte, its economic impact and interactions among farms in their efforts to control infestations. We supposed there were j = 1, …, J different farms in the area where the same harvester is used. Each farmer j can choose among a set of control strategies i = 1, …, I involving a variety of efforts ei for the plot, which are related to a given cost function ci(·). The control effort ei was measured by the number of hectares controlled under strategy i. We defined the functions of the model depending on two state variables: w (weeds, in plants per m2) and s (seed bank, in number of seeds per m2); thus, we omitted the other inputs and costs related to the production process to focus attention on weed-crop interaction. The private decision problem of farmer j was the following:
M a x e i j B t j = i = 1 I v i , t j ( w i , t j ) c i , t j ( e i , t j )
s . t .   w i , t + 1 j = W ( s i , t j , e i , t j )
s i , t + 1 j = S ( s i , t j , w i , t j , e i , t j , k j J s i , t k j ( e i , t k ) )
where Bjt is the benefit from the control effort ej made by farmer j in period t; vi is the profit margin obtained from the crop production under strategy i (profit margin depends on weed); cji is the cost of control on site j, which depends on the control effort made by j; wij is the weed function, which depends on seeds and control effort ei; sji is the seed dynamics function, which depends on the control effort made on farm j, and also on the effort made on the other farms ( e i , t k , with kj). The externality in teosinte control is expressed by the first derivative of function S(·) with respect to ek, which is considered negative, in other words, the greater the control effort in plot k, the smaller the seed bank in plot j. Thus, the increase in the control effort in site k in period t affects j by increasing its margin through a reduction in sj at period t + 1, which in turn affects the number of plants wj. The first-order conditions of this problem are:
v j w j × w j s j × s j e j + v j w j × w j e j c j e j = 0
implying that, in each period, farmer j will equate the private marginal benefits and costs of the control effort.
The social decision problem takes the form:
M a x e i j S B t = j = 1 J i = 1 I v i , t j ( w i , t j ) c i , t j ( j = 1 J e i , t j ) D ( j = 1 J e i , t j )
subject to Equations (2) and (3), where SBt denotes the social benefit in the total area of control and function D(·) includes the public costs accruing to the control program to manage teosinte infestations set by the regulatory agency. We included these public costs by formulating a linear D(·) function, which depends on the control effort made by all the farmers and will be explained in detail in the numerical application below. The necessary first-order conditions to solve this problem include:
v j w j × w j s j × s j e j + k j J v j w j × w j s j × s j s k × s k e k   + v j w j × w j e j c j e j j = 1 J D e j = 0 ,     j
where the term k j J v j w j × w j s j × s j s k × s k e k measures the impact that the control effort made in farm k has on the margin of all other neighboring farms, while the term j = 1 J D e j is the sum of the impacts of individual control efforts on the public cost incurred by the regulator (i.e., marginal public costs avoided). We then obtained a classical result of the economics of invasion: the spread of any invasive species depends on the control efforts made by all those affected. In the private decision problem, the decision makers have no incentive to consider the effects of their actions on others, while solving the social problem requires that they do. The implication is again a well-known result: it is socially desirable to encourage farmers to consider all the impacts of their control efforts and to adopt the socially optimal level. In our specific case, in which we consider the shared use of the harvester as a source of expansion of the invader to neighboring farms, we explored existing incentives to encourage farmers to adopt optimal control strategies with and without the need for intervention by a central regulator.
To solve the problem numerically, we needed to estimate the marginal external benefits of control, in other words, the value of the term v j w j × w j s j × s j s k × s k e k and the specification of the functions included in the model.

2.2. Numerical Illustration: Data and Study Area

Our study was based on data obtained from 2014 to 2018 in the infested areas of Aragon, where the presence of teosinte was first detected. This region is one of the main corn-monocropping zones in Spain, with 20% of the country’s total production.
The numerical solution of the bio-economic models required specifying the biological functions and parameters. For a better understanding of the dispersal mechanisms, we simplified our theoretical model to the case of J = 2 different farms, j and k. We supposed that field k is initially infested at a given infestation rate (low or high), while field j is initially noninfested. Thus, the initial values of weeds and seeds considered are w0k = 0.001 and s0k = 0 for a low infestation rate, and w0k = 0.1 and s0k = 0.007 for a high infestation rate, while w0j = 0, s0j = 0. We then focused on the effects that the control decisions made by farmer k have on farm j only considering the case when the harvester is used by farmer k and then by farmer j in the same period t. Figure 1 schematically illustrates the main biological processes in annual teosinte population dynamics.
The schematic diagram includes the three plant phenological stages of teosinte (S: seeds; Sl: seedlings; and W: adult plants) and the corresponding biological growth processes (e: emergence; d: development; and F: seed production). In addition, the life cycle contains the evolution of the soil seed bank, which is dominated by emergence (e) and seed survival parameters (ss). This dynamic is affected by the available control strategies which affect the development process (d) and the seed survival capacity (ss) in different ways [27].
Assuming that the invader’s dispersal source is the harvester used in the first infested plot k, we considered there is a probability piin that the harvester will move teosinte seeds from the infested field k to the clean field j, that is, the probability of infestation is piin ≠ 0 when corn is grown and piin = 0 when rotations or alfalfa are cropped.
In accordance with the data recorded in the study area, the probability of infestation was drastically reduced (but not completely eliminated) when the harvester is cleaned after use in an infested plot and before harvesting the next. These probabilities of infestation (with and without cleaning the harvester) were, however, unknown. Therefore, farmer j makes a decision by calculating the expected benefit (EBtj), which is the average benefit weighted by the probability of infestation. Consequently, the variable Btj of Equation (1) is substituted by EBtj, which is defined as follows:
E B t j = B t j ( i n f e s t e d ) × p i i n + B t j ( n o n i n f e s t e d ) × ( 1 p i i n )
Since the probability of infestation piin depends on the cropping plan of farmer k, when piin = 0, we obtained the case in which farmer j maximizes benefits without infestation, given that w 0 j = 0 and s 0 j = 0, in other words, the expected benefit of j is the same as without infestation E B t j = B t j ( n o n i n f e s t e d ) .
We incorporated the decision to clean the harvester linked to the control alternatives i that include the cultivation of corn (i.e., no control, false seedbed and manual control) and we supposed that farmer k covers the cleaning cost. Thus, the model comprises 9 alternative control strategies (6 with corn crop and 3 rotations without corn):
(1)
No control–no cleaning (corn crop).
(2)
No control–cleaning (corn crop).
(3)
False seedbed technique–no cleaning (corn crop).
(4)
False seedbed technique–cleaning (corn crop).
(5)
Manual control–no cleaning (corn crop).
(6)
Manual control–cleaning (corn crop).
(7)
Barley–sunflower rotation.
(8)
Pea–sunflower rotation.
(9)
Alfalfa.
We initially set a value piin = 0.5 for strategies in which the harvester is not cleaned (i = 1, 3, 5), piin = 0.1 for strategies in which the harvester is cleaned (i = 2, 4, 6) and piin = 0 for crop rotations (i = 7, 8, 9).
The possibility of plot j becoming infested was incorporated into the model through the specification of the seed bank function S(·) (Equation (3) in the model). In particular, we modified the formulation from [17] by incorporating the fact that the seed bank in site j is replenished by seeds from the infested field k. Thus, Equation (3) of the model takes the specific form:
s i , t + 1 j = S ( s i , t j , w i , t j , e i , t j , k j J s i , t k j ( e i , t k ) ) = g ( s i , t j , w i , t j ) × e i , t j + ε i × s i , t k × e i , t k
where function g(·) represents the dynamics of the seed bank under farmer j control, while the second part of the sum includes the effect of the control decisions made by farmer k on plot j (i.e., the spatial externality). Parameter εi denotes the portion of teosinte seeds that leave plot k in the harvester and may infest plot j. A study of the infestation patterns of the plots in the area showed that in 66% of the fields teosinte very probably arrived by means of the harvester [22]. Although the value of ε is unknown, [23] found that only 1% of teosinte seeds come out of the harvester hopper, which may indicate that a smaller percentage could be transported in the harvester slits to other plots. Therefore, we set a value of ε = 0.5%. When the harvester is cleaned after use in an infested plot, this value is considered significantly reduced, although it was difficult to quantify. We set ε = 0.25% in that case, in other words, cleaning the harvester results in a 50% reduction of the seeds that are left in the harvester and can potentially infest plot j.
Hence, we consider that cleaning the harvester affects 2 parameters in the model: the value of ε, and the probability of infestation pin. In subsequent sections, we subject our results from both parameters to a sensitivity analysis.
For our numerical illustration, we considered a planning horizon of T = 15 years—to capture the primary biological and economic aspects of controlling teosinte—and a discount rate of 3%. Moreover, we included additional restrictions concerning average farm size (8 ha) in the models and, in the case of social decision problems, a crop rotation restriction, which is a mandatory measure introduced by the CSCV only in areas with highly infested plots and it means that a crop cannot be planted in the same plot for more than 1 year in a row, with the exception of alfalfa, which remains for 5 seasons. Furthermore, harvester cleaning is mandatory for plots with low infestation levels in the social model to capture the current situation. According to the data obtained in the study area, harvester cleaning costs EUR 120 [28].
Table A1 in Appendix A shows the complete description of the model functions and the values of the biological and economic parameters related to each control strategy. Both private and social problems were solved with the CONOPT2/GAMS algorithm (General Algebraic Modeling System) [29].
The biological data used in the model were obtained from the experimental trials started in 2014 to investigate the biology of teosinte in the growing conditions found in the study area under a research project funded by the Spanish National Agriculture Research Institute (INIA). Some research from this effort has been already published ([19,20,21,22,23]) and further related work is ongoing. In addition, the economic model incorporates actual data obtained from 2014 to 2018 by the CSCV on infested areas, farmer behavior, actual evolution of the invasive species in the affected regions and actual costs of monitoring.
Thus, the model was calibrated to capture the main phenological stages of teosinte and the results were validated using actual data on farmers’ behavior and teosinte evolution in the infested areas. In addition, the robustness of the results has been strengthened by performing a sensitivity analysis on the parameters linked to uncertainty.

3. Results

Obtaining the numerical solution of the models allows us to: (i) identify and compare the pattern of private versus social optimal control strategies; (ii) quantify the value of the externality caused by farmer k to farmer j; (iii) explore the possibility of cooperation among farmers to internalize the external effect through side payments; and (iv) estimate the loss associated with the presence of teosinte in the area and the avoided costs when socially optimal control strategies are adopted.
Achieving the specific objectives (ii) to (iv) requires the quantification of the expected benefits of farmer j (i.e., Equation (7)) under each infestation scenario, which involves obtaining the farmer’s optimal strategies when infested and without teosinte. In the real context, monocropping practices are allowed if there is no evidence of infestation; therefore, in the absence of teosinte evidence, farmer j maintains corn monocropping. Consequently, we must focus on the response of farmer j to the infestation caused by seeds from neighboring plot k. For this purpose, we calculate the optimal strategies of farmer j by assuming that plot j always becomes infested when the corn is grown by farmer k, in other words, piin = 1 for i = 1, …, 6 and the values of parameter εi are constant. Subsequently, the expected benefit of j will be calculated by modifying the probability of infestation as indicated in the previous section.

3.1. Private vs. Social Optimal Control Strategies

The problem defined in Equations (1) and (3) is solved for farmers j and k to provide the noncooperative solution, that is, the optimal private decision rule for farmers under both initially low- and high-infestation scenarios assuming that plot j is initially noninfested in t = 1. Figure 2 shows the optimal strategy pathway when farmer k adopts the controls that maximize private benefits, in other words, when the external costs/benefits caused to farmer j—whose plot is not initially infested—are not considered and cleaning the harvester is mandatory with low infestation levels.
In accordance with Figure 2, with low infestation (top cells in Figure 2), farmer k would select the “no control–cleaning” strategy for the first 3 years and then adopt “manual control–cleaning” in year four. From year five to year nine, corn would be substituted by alfalfa and then farmer k would return to corn monocropping in year 10. As a consequence of this noncooperative behavior of k, at the end of year one farmer j would receive teosinte seeds when using the harvester, and in year three this farmer would identify the presence of teosinte plants in the field. Figure 2 shows that farmer j adopts rotations one year after farmer k.
In the case of high infestation (bottom cells in Figure 2), the model establishes that farmer k would select the “no control–cleaning” strategy for 2 years, “false seedbed technique–cleaning” in the third year and alfalfa for 5 years. Afterwards, corn monocropping would be restored in year nine, because rotations are not mandatory in this context. Consequently, farmer j would detect the presence of teosinte plants in year three and would select “no control–cleaning”, “false seedbed–cleaning” and rotations in year five.
These noncooperative strategies in plots (with an initial low infestation level in plot k) cause weed and seed density to increase until year five (data not shown), when all farmers have to introduce alfalfa due to the strong competition between teosinte and corn. The complete elimination of teosinte is achieved in year 10. Compared to low-infestation strategies, highly-infested plots adopt alfalfa one year earlier, thus allowing the eradication of teosinte in year nine on all farms.
These results are consistent with those obtained in [27], in the context of the optimal private decision, and demonstrate that farmers only pay the cleaning cost when it is mandatory.
However, our results show that when we include in the model the existence of an external effect on the control of teosinte, the timing when farmers adopt those control methods changes. Thus, the timing of the adoption of controls affects the period in which total eradication of the invader is achieved in a particular plot. This result may be a plausible explanation for the staggered detection of new infested fields in the study area over time. Of course, this issue influences the costs to be borne by the regulator and highlights the myopic behavior of agents when they do not consider the consequences of the spatial expansion in their decision-making process.
The social model is solved for two cases: (1) In the current situation, in which rotations are mandatory only when the plots have high infestation levels. (2) Considering the case in which crop rotations would be mandatory in all the plots for agronomic reasons, that is, with the aim of preventing future infestations, diseases and other phytosanitary problems. The second case serves to draw conclusions about the potential loss of profits for farmers due to rotations, and about the costs avoided by these measures in potential infestation events such as teosinte.
Figure 3 presents the results when the social problem defined in Equations (5) and (6) is solved with mandatory rotations only for highly infested plots and monocropping is permitted under no infestation. In this case, the model indicates that farmer k should adopt rotations in the first year, both with low and high infestation levels, since growing corn in plot k would imply a nonzero probability of infestation in plot j; therefore, if cleaning the harvester does not ensure that the risk of infestation is eliminated, then rotation is the only way to avoid teosinte dispersal and the associated public costs.
As a result, farmer j would not receive teosinte seeds and could, therefore, continue with corn monocropping. This means that cleaning the harvester is not socially optimal.
Figure 4, however, illustrates the optimal control strategies if crop rotations are always mandatory. In this case, plot k with low infestation selects the “no control–cleaning” strategy in year one and adopts rotations from years two to six with half the area devoted to alfalfa and the other half alternating with pea–sunflower and barley–sunflower. In year seven, the corn crop can be planted again because teosinte and its seed bank have been eradicated. From year eight, the area allocated to alfalfa is planted with pea–sunflower and corn in alternating years.
In addition, farmer j would replicate farmer k’s strategies, but may include corn in the rotations from the beginning, since cleaning the harvester would reduce the teosinte seeds k receives and rotations would prevent their proliferation.
When plot k is initially highly infested, the results suggest that plot k would adopt rotations starting in year one and could return to a corn crop by year six. In that case, farmer j would not receive teosinte seeds due to the immediate incorporation of rotations by k.
Comparing these results with those of the private noncooperative model (Figure 2) leads to two important conclusions: first, under the assumption that cleaning the harvester reduces the probability of infestation, the obligation to clean in infested plots only delays the detection of infestation in other initially noninfested plots and, therefore, from a social perspective, this measure is not appropriate in the current situation, in which rotations are not always mandatory; and second, the obligation to adopt rotations in highly infested plots ensures that the invader is totally eradicated after 5 years and also prevents the spread to other plots. In agreement with [27], the results also confirm that manual control and false seedbed techniques are not socially optimal. Whether mandatory rotations are adopted has relevance to the control costs covered by the regulator, which are quantified in the section below.
Obviously, the above results depend significantly on the value of parameter ε, in other words, the percentage of seeds considered to enter plot j from plot k, which, as already mentioned, is difficult to determine precisely. For this reason, conducting a sensitivity analysis of the results of this parameter seems appropriate.
First, we vary the value of εi from 0.001% to 1% for no cleaning strategies (i = 1, 3, 5) and then we modify it between 10% and 90% (multiplying by 0.9 to 0.1) of its initial value and assign it to the cleaning strategies (i = 2, 4, 6). We record the effects of these variations on private strategies.
The analysis shows that the value of εi only influences the period in which plot j detects the presence of teosinte plants but not the optimal private strategies. In particular, values between 0.2% and 0.8% increase the number of seeds coming from plot k but do not influence the timing of rotations. For the private model, we find that only the values of εi ≥ 0.9% anticipate the moment when the presence of teosinte plants is detected in plot j and, thus, rotations are advanced one period (until year five for low infestation and six for high infestation) due to the detection of more teosinte plants in the field. Our analysis shows that only for εi ≥ 2% does plot j present low infestation (w ≥ 0.001 plants m−2) two periods earlier (year three and four for low and high infestation, respectively) and the rotations are adopted two periods earlier than with the initial parameter value. However, as mentioned above, it seems unrealistic to expect εi ≥ 1%.
Conversely, values of εi ≤ 0.01% delay the adoption of rotations for one period (until year seven with low infestation or year eight with high infestation) and εi ≤ 0.001% delay rotations for two periods (until year eight or nine for low-and high-infestation scenarios, respectively).
Regarding the sensitivity of the results to the εi value for cleaning strategies, we consider variations from 10% to 90% (i = 2, 4, 6) with respect to the value for no cleaning strategies. The analysis shows that if cleaning the harvester reduces the value of εi by less than 60%, then rotations in plot j are adopted in the same year. In contrast, if cleaning reduces the value of εi by more than 60%, then the rotations are delayed by one year. Values of εi ≤ 0.001% would delay rotations for two periods. For example, considering the initial value for no cleaning strategies εi = 0.5% (i = 1, 3, 5), if cleaning the harvester reduces εi (i = 2, 4, 6), given that 0.001 ≤ εi ≤ 0.2, then the rotations are delayed for one period (to year six for the low-infestation scenario and to year seven for the high-infestation scenario), while if εi < 0.001%, rotations are delayed for two periods (to year seven or eight). Table 2 summarizes the previous results and shows the sensitivity of the benefits of j to changes in the values of parameter εi with and without cleaning the harvester.
In summary, the value of parameter εi influences the staggered detection of teosinte, which leads to the timing of control measures, total eradication of the invader and, consequently, public costs incurred in monitoring and controlling the invader as well.
The analysis shows that the use of shared harvesters, a common practice in many cereal-growing areas in the study area, makes it difficult to control infestations due to the risk of spatial dispersal. Although cleaning the machines reduces the risk of dispersal to neighboring plots, the results show that this practice does not ensure total eradication of teosinte and should be replaced by mandatory rotations also under low infestation if monocropping is not prohibited in the whole area. Alternatively, since crop rotations are the only measure capable of eradicating the invader and its spread, the regulator should insist on measures to promote the replacement of monoculture as a way of protecting against invaders and other potential phytosanitary hazards.

3.2. Quantifying and Solving Externalities

The estimate of the benefits obtained by farmers j and k with different infestation situations and its comparison with the case without infestation provides the quantification of the externalities generated in the control of teosinte.
As shown in the theoretical section above, the total external cost comprises two parts: a private cost covered by j and a public cost covered by the regulator. To quantify the former (the private external cost), we need to focus on the case in which farmer k selects the first strategy (no control–no cleaning) and how this behavior affects the benefit obtained by farmer j; in other words, we need to numerically calculate the losses farmer k causes to farmer j at both low and high infestation levels. This situation, in which farmer k selects the optimal private strategies, is shown in Figure 2. The private external cost is the difference between the benefits that j would obtain with no infestation minus those that j obtains after being infested by k. The public cost is calculated through the damage function defined in Equation (5).
Table 3 shows the average annual benefits per hectare obtained by farmers with low and high infestation levels and optimal private and social strategies, associated public costs and the quantification of external costs.
Results show that adopting optimal private strategies generates EUR 1041.40/ha for farmer k and EUR/ha for farmer j with low infestation and EUR 1006.30/ha and EUR 1009.80/ha with high infestation, respectively. This means farmer j would have an average annual loss of EUR 322/ha when farmer k initially has a low infestation, and EUR 364.20/ha when farmer k initially has a high infestation. These amounts correspond to the quantification of the private external cost when j becomes infested assuming that farmer k adopts the appropriate optimal private strategies (Figure 2). Furthermore, strategies adopted by k imply an annual public cost amounting to EUR 49.50/ha with low infestation and EUR 37.60/ha with high infestation, because rotations are adopted earlier with high infestation. In the case of farmer j, the public costs amount to EUR 37.60/ha and EUR 21.10/ha for low and high infestation levels, respectively. Public costs are lower for farmer j compared to k due to delayed infestation.
When adopting socially optimal strategies in the current context (column 2 in Table 3), farmer k would obtain an average annual benefit of EUR 999.50/ha with both low and high infestation levels, while the benefit of farmer j is not to have teosinte. Public costs due to teosinte would be eliminated in this situation.
Finally, when rotations are mandatory, average annual benefit without infestation amounts to EUR 920/ha since the model includes crop rotation restrictions even without infestation. In this case, if farmer k adopts optimal social strategies, benefit is reduced to EUR 693.50/ha for low infestation and to EUR 623.40/ha for high infestation, but the private external cost borne by j is eliminated. These strategies only generate public costs with low infestation levels (EUR 12.90/ha) in plot k because the immediate adoption of rotations with high infestation levels in plot k eliminates public costs and prevents spread to plot j.
The estimates shown in Table 3 make it possible to explore several policy responses to the externalities identified in teosinte control. The externalities resolution theory establishes that the agent causing the negative (positive) externality must internalize all the costs (benefits) of the activity.
As we indicated above, these external costs can be eliminated in several ways: by introducing mandatory rotations (even in the case of low infestation), by designing taxes for farmers who plant corn and suffer teosinte infestations, etc. However, here we focus on the possibility of establishing agreements among farmers to reduce/eliminate the risk of teosinte dispersal.
If we look at the case of noncooperation in Table 3, the optimal strategies of farmer k cause farmer j’s private losses of EUR 322/ha in the event of infestation, and no losses with no infestation. Since the expected losses of j amount to EUR 161/ha when the probability of infestation is 0.5, we look at the possibility of both farmers agreeing to reduce or eliminate these losses in the context of uncertainty.
Given that cleaning the harvester reduces the probability of infestation pin at a certain but unknown level, we estimate the average annual expected benefits for farmer j and the losses with respect to the noninfestation scenario considering different values of pin ranging from 0 (no probability of infestation) to 1 (total probability of infestation) with an initial interval of 0.1. Next, we expand the analysis for an interval of 0.01 after verifying that the results are sensitive to this range. Although the complete results of the analysis are shown in Appendix A (Table A2), Table 4 contains the results for low infestation and pin variations of 0.1 under the initial values for parameter εi and also for pin = 0.48 and pin = 0.49 to illustrate what happens with pin variations of 0.01.
By calculating the expected losses of j with respect to the situation of noninfestation, we can evaluate this farmer’s economic incentive to reduce/avoid the risk of infestation coming from plot k. For example, if cleaning the harvester reduces the probability of infestation to pin = 0.1, then the total discounted expected losses of j would be EUR 3834 compared to the current EUR 19,320 (pin = 0.5). Since the total discounted cost of cleaning the harvester amounts to EUR 459.42 (EUR 120 for the first four periods), j has an incentive to cover the total cost of cleaning.
This type of evaluation allows us to identify the pin values for which farmer j would have incentives to pay for the cleaning of the harvester. Table 4 shows that farmer j has incentives to promote an agreement whenever cleaning the harvester reduces the probability of an infestation in the range of 0.1 in probability of occurrence. When cleaning the harvester eliminates the probability of infestation (pin = 0), j receives a total discounted benefit of EUR 164,420.50 if the cleaning cost is covered, which is a higher benefit than in any situation if the harvester is not cleaned.
In the case of smaller variations of probability (0.01), we find that the incentive to establish agreements by j does not occur if the probability of infestation is reduced by 0.01 (one interval step). For example, if pin = 0.5 when the harvester is not cleaned and pin = 0.49 when the harvester is cleaned, the expected benefit from cleaning for farmer j amounts to EUR 386, while the cleaning costs come to EUR 459.42; therefore, j would not be willing to pay all these costs.
Since the possibility of an agreement depends on the quantification of externalities and losses caused by the presence of teosinte, which in turn depend on the value of parameters pin and εi, we subjected our results from both parameters to a sensitivity analysis, as shown in Appendix A (Table A2).
The data in Table A2 demonstrate that the incentive to cooperate is maintained even if the infestation in plot j is delayed as a result of the smaller value of εi. In particular, reductions of more than 0.01 in the probability maintain the incentive for j to cover the total cleaning cost, while reductions of 0.01 maintain it only partially. The analysis reveals that there are several chances of a cooperative agreement among farmers but these are linked to the effectiveness of cleaning in reducing the probability of infestation.
In contrast, if the regulator decides that the cost of cleaning should be covered by farmer k, then the results of the Table A2 allow the identification of additional agreement possibilities among farmers, since farmer j would be willing to pay at least part of the cleaning cost provided that the expected losses resulting from becoming infested are higher than the cleaning cost. For example, if cleaning the harvester implies that pin is reduced from 0.5 to 0.49, then farmer j would be willing to pay a positive amount with a maximum of EUR 386, which is the difference in the expected losses, while farmer k would be willing to accept any contribution to the cost of cleaning.

3.3. Estimating the Loss Associated with the Presence of Teosinte

Finally, the numerical data obtained in Table 3 for farmer k make it possible to estimate the losses associated with the presence of teosinte in the study area. For this purpose, we combine this information (Table 3) with the data provided by the CSCV regarding the number of infested hectares and the infestation level (Table 1).
Table 5 shows that private strategies lead to total losses over the whole period of EUR 18,271.40·103 relative to a noninfestation scenario; EUR 16,120.40·103 of this amount are private costs and EUR 2,151·103 are public costs covered by the regulator. Additionally, these strategies imply costs for the neighboring plots using the same harvester. These losses can be assessed by calculating the losses expected by farmer j when first infested (results shown in Table 4), to which the associated public costs must be added. Based on our model’s parameter values, these annual expected profit losses amount to EUR 161 and EUR 166.30/ha for low and high infestation, respectively; the expected public costs amount to EUR 18.80/ha and EUR 10.50/ha (half the total public cost shown in Table 3) for low and high infestation, respectively.
Adopting social strategies in the current context leads to losses with respect to the situation of noninfestation amounting to EUR 17,543.40·103, which corresponds to a loss of private profits, as the public costs would be eliminated. However, in this scenario, the risk of future weed infestations or other diseases is not avoided, since monocropping is still permitted. In contrast, the data show that permanent rotations lead to minor losses with respect to the noninfestation situation (EUR 12,095.10·103), while minimizing the risk associated with future diseases.
Obviously, this option implies private benefit losses of EUR 21,266.60·103 with respect to noninfestation, which explains why the regulator usually recommends permanent rotations without imposing them. In fact, [4] describe that managers’ increased risk aversion may lead to less prevention and more control in situations in which there is uncertainty regarding the probability of being affected by an invader or a disease, because prevention does not eliminate the invasion, it only reduces the likelihood of becoming invaded. The economic explanation is that the manager values euros spent on control (with a certain benefit now) more highly than euros spent on prevention (with an uncertain benefits) [4]. Our results confirm this behavior in the area of study.
The results show that monoculture is a practice that many farmers adopt in their fields due to profit-maximizing behavior and the industrial nature of crop production in the study area. However, considering a more global ecological and economic context, rotations protect against invasive infestations, pests, diseases and help maintain soil fertility. Thus, rotations reduce the use of agrochemicals, which increases farmers’ long-term benefits and results in greater sustainability of the agroecosystem [30,31]. In addition, as our results show, crop rotations substantially reduce the need for monitoring and control by the regulator, reducing public costs as a result.
To enhance the sustainability of the system, the regulator must incentivize farmers to adopt crop rotations. Our work estimates that the adoption of rotations would save EUR 5,448.3·103 in the control of teosinte from 2014 to 2018, without taking into account other advantages linked to soil fertility, pest and other weed controls.

4. Conclusions

The comparison of the results of the present study with those obtained in the previous work showed that a profit-maximizing farmer would never cover the cost of cleaning the harvester after using it in the infested plot when it is not mandatory.
Assuming that cleaning the harvester does not completely eliminate the probability of an infestation, we found that this measure is not socially optimal in the current situation, where corn-monocropping practices are not prohibited in noninfested plots. In this case, the model indicated that infestation and detection of the invader in initially noninfested fields is delayed over time, so that eradication is also delayed. Therefore, the pattern of temporal and spatial teosinte expansion observed in the affected areas can be explained, at least in part, by the use of the harvester, which causes seeds to enter new plots and increases the overall costs of controlling and eradicating teosinte.
The data suggested that if a monocrop is allowed, infested plots should adopt mandatory rotations, even when infestation is low, as this is the only way to avoid spatial dispersal of the invader. In addition, our results confirmed that manual control and seedbed techniques are not socially optimal so the regulator should reconsider these measures to control and eradicate teosinte.
The sensitivity analysis showed that a change in the value of parameters associated with the probability of infestation and the number of seeds before and after cleaning the harvester does not change the optimal strategies but can delay/advance their adoption.
Adopting socially optimal strategies would eliminate the externalities caused to the neighboring plots and also the public costs for the regulator, although monoculture maintains the risk of sustaining future phytosanitary problems. Quantifying the benefits and costs associated with several situations allowed us to affirm that adopting permanent rotations would avoid an annual average cost of EUR 1235.20·103.
Of course, the results depend heavily on the model’s ability to represent reality and on the data available for the study area. However, our methodology can be easily adapted to the study of other crops, weeds, pests or diseases. On the other hand, the data used in this work correspond to the only case of teosinte infestation reported in Europe for which there are available biological and economic data, so the results can be useful for other areas where teosinte could appear in the future. Another promising future research area is including in the model (through a function or specific parameters) the relationship between rotations and long-term indirect benefits discussed above (i.e., improved fertility, ease of controlling other diseases, reduced need for regulator control, etc.).

Author Contributions

Conceptualization, G.P. and Y.M.; methodology, G.P., M.I.G. A.C. and Y.M.; software, Y.M.; validation, G.P., A.C. and Y.M.; formal analysis, M.I.G. and Y.M.; investigation, G.P., A.C. and Y.M.; resources, G.P. and A.C.; data curation, G.P. and Y.M.; writing—original draft preparation, G.P., A.C. and Y.M.; writing—review and editing, G.P. and Y.M.; funding acquisition, G.P., M.I.G., A.C. and Y.M. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Spanish Instituto Nacional de Investigación Agraria (INIA), grant numbers ERTA-2014-00011-C02 and RTA2017-00082-00-00, and the Spanish Ministerio de Economía, Industria y Competitividad, grant number ECO2016-75927-R.

Acknowledgments

We would like to thank Fernando Arrieta, Jorge Pueyo, Joaquín Aibar and José A. Alins for their help in collecting experimental data in the area of study.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

Table A1. Model functions and parameters.
Table A1. Model functions and parameters.
FunctionSpecificationParameter ValuesSource
Weed dynamics w i , t + 1 j = w * [ 1 e x p ( α 0 ( α 1 + s i , t j ) ) ] · e i , t j w * = 22[27]
α0 = 0.0704
α1 = 0.1876
Seed bank dynamics g ( s i , t j , w i , t j ) = { β 1 s i , t j + β 2 w i , t j   if   s i , t < s * s *   if   s i , t j s * s * = 31.8[27]
β1 = 0.0738
β2 = 98.97
Profit margin v i , t j ( w i , t j ) = δ 0 , i + δ 1 , i w i , t j δ0 = 11.334 if i = 1, …, 6[23,27]
δ0 = 374.12 if i = 7
δ0 = 505.36 if i = 8
δ0 = 547.76 if i = 9
δ1 = −0.5456 if i = 1, …, 6
δ1 = 0 if i = 7,8,9
Control costs c i , t j ( e i , t j , e i , t k ) = γ 0 , i j + j = 1 J γ 1 , i j · e i j γ 1 , i j = 0 if i = 1,2Own from available data
γ 1 , i j = 546.7 if i = 3, 4
γ 1 , i j = 142.8 if i = 5,6
γ 1 , i j = 0 if i = 7, 8, 9
γ 0 , i k = 120 if i = 2, 4, 6
γ 0 j = 0, jk
Public costs D i ( e i ) = d 0 + d 1 , i · ( e i j + e i k ) d0 = 1900 if i = 1, …, 6Own from available data
d1 = 160 if i = 1, …, 6
d1 = 0 if i = 7, 8, 9
Total individual land restriction i = 1 9 e i j = E ¯ j     j E ¯ j = 8 Own from available data
Rotation restriction e i , t j z = 1 8 e z , t 1 j   w i t h   i z     i , z = 1 , , 9 [27]
Table A2. Results of sensitivity analysis for εi and pin parameters.
Table A2. Results of sensitivity analysis for εi and pin parameters.
εi Values When No-Cleaning (i = 1, 3, 5) and Multiplier ValuesProbability of Infestation (pin) When No-CleaningProbability of Infestation (pin) When CleaningAnnual Average Benefits under Infestationn (EUR/ha)Expected Annual Average Benefits(EUR/ha)Total Discounted Average Benefits (EUR)Total Expected Losses with Respect to Non-Infestation (EUR)
0.1 < εi < 1 (i = 1, 3, 5) and multiplier ≥50% or 0.1 < εi < 0.9 and multiplier <50%101054.71374164,8800
0.990.011054.71370.80164,496.84383.16
0.980.021054.71367.61164,113.68766.32
0.970.031054.71364.42163,730.521149.48
0.960.041054.71361.22163,347.361532.64
0.950.051054.71358.03162,964.201915.80
0.940.061054.71354.84162,581.042298.96
0.930.071054.71351.64162,197.882682.12
0.920.081054.71348.45161,814.723065.28
0.910.091054.71345.26161,431.563448.44
0.90.11054.71342.07161,048.403831.60
0.80.21054.71310.14157,216.807663.20
0.70.31054.71278.21153,385.2011,494.80
0.60.41054.71246.28149,553.6015,326.40
0.50.51054.71214.35145,72219,158
0.40.61054.71182.42141,890.4022,989.60
0.30.71054.71150.49138,058.8026,821.20
0.20.81054.71118.56134,227.2030,652.80
0.10.91054.71086.63130,395.6034,484.40
011054.71054.70126,56438,316
0.9 < εi < 1or 0.001 < εi ≤ 0.001 with any multiplier1010641374.00164,8800
0.990.0110641370.90164,508372
0.980.0210641367.80164,136744
0.970.0310641364.70163,7641116
0.960.0410641361.60163,3921488
0.950.0510641358.50163,0201860
0.940.0610641355.40162,6482232
0.930.0710641352.30162,2762604
0.920.0810641349.20161,9042976
0.910.0910641346.10161,5323348
0.90.110641343.00161,1603720
0.80.210641312.00157,4407440
0.70.310641281.00153,72011,160
0.60.410641250.00150,00014,880
0.50.510641219.00146,28018,600
0.40.610641188.00142,56022,320
0.30.710641157.00138,84026,040
0.20.810641126.00135,12029,760
0.10.910641095.00131,40033,480
0110641064.00127,68037,200
1010731374.00164,8800
0.0001 < εi < 0.001 with any multiplier0.990.0110731370.99164,518.8361.20
0.980.0210731367.98164,157.6722.40
0.970.0310731364.97163,796.41083.60
0.960.0410731361.96163,435.21444.80
0.950.0510731358.95163,0741806
0.940.0610731355.94162,712.82167.20
0.930.0710731352.93162,351.62528.40
0.920.0810731349.92161,990.42889.60
0.910.0910731346.91161,629.23250.80
0.90.110731343.90161,2683612
0.80.210731313.80157,6567224
0.70.310731283.70154,04410,836
0.60.410731253.60150,43214,448
0.50.510731223.50146,82018,060
0.40.610731193.40143,20821,672
0.30.710731163.30139,59625,284
0.20.810731133.20135,98428,896
0.10.910731103.10132,37232,508
0110731073128,76036,120
εi < 0.0001 with any multiplier101081.781374.00164,8800
0.990.011081.781371.07164,529.34350.66
0.980.021081.781368.15164,178.67701.32
0.970.031081.781365.23163,828.011051.99
0.960.041081.781362.31163,477.341402.65
0.950.051081.781359.38163,126.681753.32
0.940.061081.781356.46162,776.022103.98
0.930.071081.781353.54162,425.352454.64
0.920.081081.781350.62162,074.692805.31
0.910.091081.781347.70161,724.023155.97
0.90.11081.781344.77161,373.363506.64
0.80.21081.781315.55157,866.727013.28
0.70.31081.781286.33154,360.0810,519.92
0.60.41081.781257.11150,853.4414,026.56
0.50.51081.781227.89147,346.817,533.20
0.40.61081.781198.66143,840.1621,039.84
0.30.71081.781169.44140,333.5224,546.48
0.20.81081.781140.22136,826.8828,053.12
0.10.91081.781111.00133,320.2431,559.76
011081.781081.78129,813.6035,066.40

References

  1. González-Andújar, J.L. Co-operative versus non-co-operative farmers’ weed control decisions in an agricultural landscape. Weed Res. 2018, 58, 327–330. [Google Scholar] [CrossRef]
  2. Fenichel, E.P.; Richards, T.J.; Shanafelt, D.W. The Control of Invasive Species on Private Property with Neighbor-to-neighbor Spillovers. Environ. Resour. Econ. 2013, 59, 231–255. [Google Scholar] [CrossRef] [Green Version]
  3. Epanchin-Niell, R.S.; Wilen, J.E. Individual and cooperative management of invasive species in human-mediated landscapes. Am. J. Agric. Econ. 2014, 97, 180–198. [Google Scholar] [CrossRef] [Green Version]
  4. Finnoff, D.; Shogren, J.F.; Leung, B.; Lodge, D. Take a risk: Preferring prevention over control of biological invaders. Ecol. Econ. 2007, 62, 216–222. [Google Scholar] [CrossRef]
  5. Perrings, C.; Dalmazzone, S.; Williamson, M. The economics of biological invasions. In Invasive Alien Species: A New Synthesis; Mooney, H.A., Mack, R.N., McNeely, J.A., Neville, L.E., Schei, P.J., Waage, J.K., Eds.; Island Press: Covelo, CA, USA, 2005; Chapter 2. [Google Scholar]
  6. Coase, R. The Problem of Social Cost. In Classic Papers in Natural Resource Economics; Palgrave Macmillan: London, UK, 1960; pp. 87–137. [Google Scholar]
  7. Liu, Y.; Sims, C. Spatial-dynamic externalities and coordination in invasive species control. Resour. Energy Econ. 2016, 44, 23–38. [Google Scholar] [CrossRef]
  8. Wilen, J. Economics of Spatial dynamic processes. Am. J. Agric. Econ. 2007, 89, 1134–1144. [Google Scholar] [CrossRef]
  9. Visintin, C.; Briscoe, N.J.; Woolley, S.N.C.; Lentini, P.E.; Tingley, R.; Wintle, B.A.; Golding, N. steps: Software for spatially and temporally explicit population simulations. Methods Ecol. Evol. 2020, 11. [Google Scholar] [CrossRef] [Green Version]
  10. Shea, K. Models for Improving the Targeting and Implementation of Biological Control of Weeds 1. Weed Technol. 2004. [Google Scholar] [CrossRef]
  11. An, C.; Liu, C.; Bi, S. Stability in Distribution and Optimal Control in an Impulsive Toxin Input Bioeconomic System with Stochastic Fluctuations and Time Delays. IEEE Chin. Control Decis. Conf. 2019. [CrossRef]
  12. Dehnen-Schmutz, K.; Perrings, C.; Williamson, M. Controlling Rhododendron ponticum in British Isles: An Economic Analysis. J. Environ. Manag. 2004. [Google Scholar] [CrossRef]
  13. Buhler, D.D.; King, R.P.; Swinton, S.M.; Gunsolus, J.L.; Forcella, F. Field Evaluation of a Bioeconomic Model for Weed Management in Soybean (Glycine max). Weed Sci. 1997. [Google Scholar] [CrossRef]
  14. Renner, K.A.; Swinton, S.M.; Kells, J.J. Adaptation and Evaluation of the WEEDSIM Weed Management Model for Michigan. Weed Sci. 1999. [Google Scholar] [CrossRef]
  15. Colbach, N.; Biju-Duval, L.; Gardarin, A.; Granger, S.; Guyot, S.H.M.; Mézière, D.; Munier-Jolain, N.M.; Petit, S. The Role of Models for Multicriteria Evaluation and Multiobjective Design of Cropping Systems for Managing Weeds. Weed Res. 2014. [Google Scholar] [CrossRef] [Green Version]
  16. González-Díaz, L.; Blanco-Moreno, J.M.; González-Andújar, J.L. Spatially Explicit Bioeconomic Model for Weed Management in Cereals: Validation and Evaluation of Management Strategies. J. Appl. Ecol. 2015. [Google Scholar] [CrossRef] [Green Version]
  17. Gonzalez-Diaz, L.; Bastida, F.; Gonzalez-Andujar, J.L. A Bioeconomic Model for the Analysis of Control Strategies for Lolium rigidum and Avena sterilis ssp. Ludoviciana in Winter Wheat. Int. J. Plant Prod. 2020, 14. [Google Scholar] [CrossRef]
  18. Grimsrud, K.M.; Chermak, J.M.; Hansen, J.; Thacher, J.A.; Krause, K. A Two-Agent Dynamic Model with an Invasive Weed Diffusion Externality: An Application to Yellow Starthistle (Centaurea solstitialis L.) in New Mexico. J. Environ. Manag. 2008. [Google Scholar] [CrossRef]
  19. Cirujeda, A.; Pardo, G.; Marí, A.I.; Fuertes, S.; Aibar, J. Emergencia de teosinte en cultivos diferentes a maíz. In Proceedings of the XVI Congress of the Sociedad Española de Malherbología, Pamplona, Spain, 25–27 October 2017. [Google Scholar]
  20. Pardo, G.; Fuertes, S.; Marí, A.I.; Aibar, J.; Cirujeda, A. Evaluación de distintos herbicidas en el control de teosinte en cultivos diferentes al maíz. In Proceedings of the XVI Congress of the Sociedad Española de Malherbología, Pamplona, Spain, 25–27 October 2017. [Google Scholar]
  21. CSCV, Centro de Sanidad y Certificación Vegetal. Current Infestation Status of Teosinte in Aragón; Information Day for Farmers, Escuela Politécnica Superior de Huesca, University of Zaragoza: Zaragoza, Spain, 2018. [Google Scholar]
  22. Montull, J.M.; Pardo, G.; Aibar, J.; Llenes, J.M.; Marí, A.I.; Taberner, A.; Cirujeda, A. Aspectos de la dispersión y viabilidad de las semillas de teosinte (Zea mays ssp.) en el Valle del Ebro. ITEA Inf. Tecnica Econ. Agrar. 2020, 116, 227–240. [Google Scholar] [CrossRef]
  23. Cirujeda, A.; Pardo, G. Aparición de una Nueva mala Hierba en el Cultivo del Maíz en Aragón: El Teosinte. Caracterización Biológica y Estudio de Métodos para su Control; Working Paper; CITA: Zaragoza, Spain, 2019. [Google Scholar]
  24. Shirtliffe, S.J.; Entz, M.H. Chaff collection reduces seed dispersal of wild oat (Avena fatua) by a combine harvester. Weed Sci. 2005, 53, 465–470. [Google Scholar] [CrossRef]
  25. Blanco-Moreno, J.M.; Chamorro, L.; Masalles, R.M.; Recasens, J.; Sans, F.X. Spatial Distribution of Lolium rigidum seedlings following seed dispersal by combine harvesters. Weed Res. 2004, 44, 375–387. [Google Scholar] [CrossRef]
  26. Pardo, G.; Cirujeda, A.; Aibar, J.; Fernández-Cavada, S.; Rodríguez, E.; Fuertes, S.; Perdiguer, A. El Teosinte (Zea mays, ssp.); Informaciones Técnicas, 4/2014; Centro de Sanidad y Certificación Vegetal, Gobierno de Aragón: Zaragoza, Spain, 2014. [Google Scholar]
  27. Martínez, Y.; Cirujeda, A.; Gómez, M.I.; Marí, A.I.; Pardo, G. Bioeconomic model for optimal control of the invasive weed Zea mays subspp. (teosinte) in Spain. Agric. Syst. 2018, 165, 116–127. [Google Scholar] [CrossRef] [Green Version]
  28. Taberner, A.; University of Lérida, Lérida, Spain. Personal communication, 2020.
  29. Brooke, A.; Kendrick, D.; Meeraus, A.; Raman, R. GAMS Tutorial by R. Rosenthal; GAMS Development Corporation: Washington, DC, USA, 2018. [Google Scholar]
  30. Lengwati, D.M.; Mathews, C.; Dakora, F.D. Rotation Benefits from N2-Fixing Grain Legumes to Cereals: From Increases in Seed Yield and Quality to Greater Household Cash-Income by a Following Maize Crop. Front. Sustain. Food Syst. 2020, 4, 94. [Google Scholar] [CrossRef]
  31. Iocola, I.; Angevin, F.; Bockstaller, C.; Catarino, R.; Curran, M.; Messéan, A.; Schader, C.; Stilmant, D.; Van Stappen, F.; Vanhove, P.; et al. An Actor-Oriented Multi-Criteria Assessment Framework to Support a Transition towards Sustainable Agricultural Systems Based on Crop Diversification. Sustainability 2020, 12, 5434. [Google Scholar] [CrossRef]
Figure 1. Life cycle diagram of teosinte.
Figure 1. Life cycle diagram of teosinte.
Sustainability 12 09046 g001
Figure 2. Optimal private control strategies.
Figure 2. Optimal private control strategies.
Sustainability 12 09046 g002
Figure 3. Optimal social control strategies in the current situation.
Figure 3. Optimal social control strategies in the current situation.
Sustainability 12 09046 g003
Figure 4. Optimal social strategies when rotations are always mandatory.
Figure 4. Optimal social strategies when rotations are always mandatory.
Sustainability 12 09046 g004
Table 1. Number of infested plots, affected area and infestation incidence in the study area.
Table 1. Number of infested plots, affected area and infestation incidence in the study area.
YearNumber of Infested PlotsNumber of New Infested PlotsArea with Low Infestation (ha)Area with High Infestation * (ha)Total Infested Area (ha)
2014444427358 (93%)385
20156327441192 (30%)633
2016701462128 (4.3%)649
2017721363428 (4.2%)662
20187340419375 (47.2%)794
* In brackets, percentage of the total area. Source: [21].
Table 2. Sensitivity analysis of the expected annual average benefit of j to the parameter εi.
Table 2. Sensitivity analysis of the expected annual average benefit of j to the parameter εi.
Values of εi for No Cleaning Strategies
(i = 1, 3, 5)
Reduction in Values of εi for Cleaning Strategies (i = 2, 4, 6)Annual Average Benefit of j under Low Infestation (in EUR /ha)Annual Average Benefit of j under High Infestation
(in EUR /ha)
0.1 < ε < 1≤50%10521009.8
0.1 ≤ ε ≤ 0.9>50%1061.31029.1
0.9 < ε ≤ 1Any value10521009.8
0.001 < ε < 0.1Any value1061.31029.1
0.0001 < ε ≤ 0.001Any value1070.41039.7
ε ≤ 0.0001Any value1079.31051.4
Table 3. Average annual benefits per hectare with cooperative vs. noncooperative strategies EUR/ha).
Table 3. Average annual benefits per hectare with cooperative vs. noncooperative strategies EUR/ha).
Optimal Private StrategiesOptimal Social Strategies (Rotations Only with High Infestation)Optimal Social Strategies (with Mandatory Rotations)
Farmer jFarmer kFarmer jFarmer kFarmer jFarmer k
(1) Benefits, noninfestation 1374137413741374920920
(2) Benefits, low infestation (losses relative to noninfestation)1052 (322)1041.4 (332.6)1374 (0)999.5 (374.5)920 (0)693.5 (365.1)
(3) Public costs, Low infestation37.649.50 *0 *012.9
(4) Total benefits, Low infestation
(4) = (2) − (3)
1014.4991.91374999.5920680.6
(5) Benefits, high infestation (losses relative to noninfestation)1009.8 (364.2)1006.3 (367.7)1374 (0)999.5 (374.5)920 (0)623.4 (374.5)
(6) Public costs, high infestation21.137.60 *0 *00
(7) Total benefits, high infestation
(7) = (5) − (6)
988.7968.71374999.5920623.4
* In this case, costs associated with the risk of future infestations with other weeds or pests are not considered.
Table 4. Effect of the probability of infestation on expected benefits and losses of farmer j under low infestation.
Table 4. Effect of the probability of infestation on expected benefits and losses of farmer j under low infestation.
Probability of Infestation (pin)Expected Annual Average Benefits (EUR/ha)Total Expected Discounted Average Benefits (EUR)Total Expected Losses with Respect to Noninfestation (EUR)
01374164,8800
0.11341.8161,0163864
0.21309.6157,1527728
0.31277.4153,28811,592
0.41245.2149,42415,456
0.481219.4146,33318,547
0.491216.2145,94618,934
0.51213145,56019,320
0.61180.8141,69623,184
0.71148.6137,83227,048
0.81116.4133,96830,912
0.91084.2130,10434,776
11052126,24038,640
Table 5. Estimates of total economic impacts in the study area from 2014 to 2018.
Table 5. Estimates of total economic impacts in the study area from 2014 to 2018.
Total Discounted Benefit/Cost (in EUR 103)
Privately OptimalSocially Optimal (Current Situation)Socially Optimal (with Mandatory Rotations)
(1) Benefits, non-infestation64,36564,36543,098.4
(2) Benefits, low-infestation area32,772.831,454.221,824.445
(3) Public costs, low-infestation area1557.70405.9
(4) Total benefit, low-infestation area
(4) = (2) − (3)
31,215.131,454.221,418.4
(5) Benefits, high-infestation area15,471.815,367.39584.7
(6) Public costs, high-infestation area593.300
(7) Total benefit, high-infestation area
(7) = (5) − (6)
14,878.515,367.39584.7
(8) Losses relative to non-infestation
(8) = (1) − (4) − (7)
18,271.417,543.412,095.1
(9) Annual average losses relative to non-infestation (9) = (8)/53654.23508.62419
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Pardo, G.; Gómez, M.I.; Cirujeda, A.; Martínez, Y. Economic Costs of Sharing the Harvester in the Control of an Invasive Weed. Sustainability 2020, 12, 9046. https://doi.org/10.3390/su12219046

AMA Style

Pardo G, Gómez MI, Cirujeda A, Martínez Y. Economic Costs of Sharing the Harvester in the Control of an Invasive Weed. Sustainability. 2020; 12(21):9046. https://doi.org/10.3390/su12219046

Chicago/Turabian Style

Pardo, Gabriel, Miguel I. Gómez, Alicia Cirujeda, and Yolanda Martínez. 2020. "Economic Costs of Sharing the Harvester in the Control of an Invasive Weed" Sustainability 12, no. 21: 9046. https://doi.org/10.3390/su12219046

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop