1. Introduction
Weeds have been a major cause of crop yield loss since the beginning of agriculture. Today, herbicide-based control methods play a key role in maximizing agrosystem productivity in the short term. However, the intensification of agriculture has led to undesirable negative consequences to both the environment and society. In this context, the combined implementation of preventive (legal, cultural) and curative (chemical, mechanical, physical, and biological) methods has been proposed as a way to mitigate externalities (soil and water contamination, biodiversity loss, ecotoxicity, etc.). Therefore, from a strategic viewpoint, an integrated weed-management (IWM) program should be based on a combination of preventive and curative methods applying knowledge-based principles. The use of cultural methods for weed management has proven to increase the competitive ability of crops, reducing their dependence on herbicides [
1]. However, integrated management approaches are still incipient in Argentina [
2].
The cost/benefit quantification of different IWM strategies is not a straightforward process due to the necessity of a large amount of information that requires further systematization to be implemented within a decision-making framework. In this context, simulation models provide an ideal approach for systematizing this type of analysis [
3,
4,
5].
A weed–crop simulation model was proposed by [
5] to support the IWM decision-making process in winter cereal crops of the semi-arid temperate region of Argentina. The model possesses a higher level of detail than similar models and, although it requires a relatively large amount of data, it could be easily adapted to represent diverse agrosystems. Therefore, the proposed model could be considered a flexible and adaptable tool.
This model uses bioecological and agronomic information as inputs, such as daily weather records, weed population dynamics data, weed-management tactics (chemical, mechanical, and cultural methods), and crops’ ecophysiological requirements. Typical results are the daily values of weed population dynamics, crop-growth/development dynamics, and the resulting weed–crop competitive interactions. At the end of each crop season, both bioecological and agronomic outputs are obtained (i.e., seed production, economic gross margin, environmental impact, etc.).
In this work, the model from Molinari et al. (2020) [
5] was extended to improve the economic and environmental evaluation of weed-management strategies. Specifically, the calculation of the present value of money was included to improve economic comparisons in multi-year simulations. Additionally, the quantification of the environmental impact was extended with the T index, which represents the soil-erosion risk associated with mechanical weed control [
6]. The P index [
6] was also added to quantify the environmental impact of pesticides, complementing the EIQ index calculations [
7].
In this study, the described model is applied to the agricultural system
Euphorbia davidii Subils in competition with soybean in the center of the Buenos Aires province (Argentina).
Euphorbia davidii belongs to the Euphorbiaceae Juss. family, represented by species of economic value and others considered to be weeds [
8,
9,
10]. Four species have been found in Argentina that behave as important weeds in summer crops (
Euphorbia serpens,
Euphorbia heterophylla,
Euphorbia dentata, and
Euphorbia davidii), sharing many common characteristics, which complicates their easy identification, and, therefore, the design of effective management strategies for each one [
11].
Euphorbia davidii has been reported as an invasive species in multiple regions of Europe [
12]. In the agrosystems of the central part of the Buenos Aires province,
E. davidii is considered a highly competitive weed that is difficult to control. In general, there is a close relationship between phenological stage, dose, and control efficacy [
13,
14]. According to [
14], under semi-controlled conditions, yield losses of 35–45% are observed in soybean crops at weed densities higher than 100 individuals.m
−2. Likewise, in the study area, field experiments indicate yield losses of 30% at 100 individuals.m
−2, with significant losses observed from 8–10 individuals.m
−2 onwards [
15].
It is hypothesized that a previously developed model [
5] can be adapted to the agricultural system composed of
Euphorbia davidii Subils in competition with soybean in the center of the Buenos Aires province, in order to support decision making for integrated weed management.
The objectives of this article are: (i) to extend the model proposed in [
5] with additional detail in the economic and environmental impact modules; (ii) to evaluate the model when applied to the soybean/
E. davidii agricultural system in the central-southern region of the Buenos Aires province; (iii) to generate annual and multiannual scenarios comparing different management strategies; and (iv) to evaluate the model’s advantages/weaknesses for its future adaptation to other agrosystems.
3. Calibration and Validation
To properly estimate the expected crop yield (Yld), parameters a and k of Equation (4) were tuned for the system under study [
3].
where Yld is the expected crop yield (as a proportion of the weed-free yield), Cs is the standard crop density, a is a crop-dependent constant, Ca is the actual crop sowing density, k is a constant reflecting the weed competitiveness of the crop, WC is the sum of the weed competitive effects on the crop at the end of the season, and Myl is the maximum yield loss proportion at high interspecific competition.
In this contribution, parameters a and k were calculated by solving a parameter estimation problem using experimental data reported in [
15,
28,
29,
30]. Field trials were conducted in the Azul district (36°47′00″ S; 59°51′00″ W), Buenos Aires province, Argentina. Different cultural management strategies (i.e., soybean crop varieties, sowing dates, row spacing, and sowing densities), as well as herbicides and mechanical control, were included. Field trials reported in [
15,
29,
30] were repeated over two crop seasons, while those reported in [
29] were carried out for a single crop season.
The available experimental data (N = 37) were divided 70/30% for calibration and validation, respectively (randomly selected). Parameters, a and k, that minimize the root-mean-square Error (RMSE) between the observed and simulated Yld were obtained using the solver add-on in the MS Excel spreadsheet version 2013 (Microsoft Corporation, Redmond, WA, USA) (a = 0 and k = 0.1, RMSE = 0.08) (
Figure 2a).
Next, we simulated the validation dataset and compared it with the observed data, obtaining an RMSE = 0.085 as shown in
Figure 2b.
5. Conclusions
In this contribution, a very detailed population-based model [
5] was extended by improving the multi-year economic calculations, and by adding new indexes to estimate pesticides and soil-erosion impact, to better compare and evaluate alternative weed-management strategies.
Following this, the model was adapted to a typical soybean/Euphorbia davidii agrosystem for the center of the Buenos Aires province, Argentina. Annual and multiannual case studies were simulated to analyze crop–weed interactions under different cultural measures and control actions (chemical and mechanical).
In general, the simulated results showed that, under high infestation conditions, it was necessary to combine: (i) an estimation of weed-emergence flowrates; (ii) the adoption of cultural management methods, such as delayed sowing times, higher sowing densities, and narrower distances between rows; (iii) chemical control methods, especially the use of a mixture of non-selective and residual herbicides, in combination with mechanical methods due to their high control rate of E. davidii at advanced development stages. By making such combinations, satisfactory agronomic outcomes could be obtained without having a high impact on gross margin and externalities due to chemical and/or mechanical actions.
While the proposed approach seems to provide a balance in terms of biological, agronomic, economic, and environmental details of the complex agrosystem under study, many improvements for future adaptations can be outlined. For example, it is known that E. davidii can coexist with several other weeds. The modelling of a multispecies agrosystem requires a great deal of specific information. Another extension that should be incorporated in future versions of the model is weed-resistance quantification, which should be considered in long-period studies for strategic and integrated weed management.