1. Introduction
Soybean (
Glycine max (L.) Merr.) is an important oilseed and protein crop on a global scale [
1,
2,
3]. The production of this crop has increased, particularly in the past several decades, when farming practices have changed dramatically, especially regarding the methods used to control weeds.
Before the introduction of herbicides, mechanical and cultural methods were the only option for weed control. Soil-incorporated and pre-emergence herbicides began to replace tillage and cultivation practices and, after the 1960s, became the dominant weed control method, followed by mechanical cultivation until soybean canopies closed and shaded the weeds [
4]. The development of post-emergence herbicides after the 1980s allowed farmers to control weeds in season and became a predominant weed control treatment that included single or multiple post-emergence herbicide applications with less reliance on soil-applied herbicides [
5].
Moreover, the availability of post-emergence herbicides gave farmers an alternative weed control tool allowing them to plant soybeans in narrow row spacing and eliminate the need for cultivation [
6]. Traditionally, soybeans have been planted in 70 cm rows, followed by late-season cultivation [
7]; however, with rows as narrow as 25 cm, cultivation becomes impossible. Soybeans planted in narrow rows mean that the canopy can close quicker, providing better competition against weeds. Weed emergence following herbicide application decreases since the solar radiation that stimulates weed seed germination and weed development is intercepted by the crop canopy [
8,
9]. In addition, the research results demonstrated that the increasing soybean seeding rate for narrow-row soybeans is 20–45% greater than for wide-row soybeans [
10,
11]. However, no difference in weed-free soybean yield at a low, moderate, and high population in 38 and 76 cm rows was observed [
12].
Soybean yield has increased over the past 35 years due to crop improvement through plant breeding [
13] and due to significantly changed production systems, such as diverse crop rotations, reduced tillage, precision seeding, optimal fertilization rate, etc. In addition, the introduction of glyphosate (Roundup Ready
®) and glufosinate (Liberty Link
®)-resistant soybeans in 1996 and 2009, respectively, resulted in a shift in practices towards postemergence herbicide application [
14,
15]. Transgenic, herbicide-tolerant varieties of soybean have had a very rapid adoption rate because of reduced weed control costs and increased yields compared to conventional varieties [
16]. Moreover, they represent a revolutionary solution in weed management since they allow farmers to manage a broad spectrum of weeds without crop injury or crop rotation restrictions.
Although herbicides have revolutionized weed management, replacing manual labor and mechanical weed control, they face some challenges, such as safety, environmental issues, and the development of herbicide-resistant weeds, particularly in countries growing genetically modified (GM) soybean varieties [
17].
Many countries in the European Union, including Croatia, restrict or prohibit the cultivation of GM crops and the active ingredients of some extremely and highly hazardous pesticides in their territory due to their cumulative long-term effects on human health and the environment [
18]. Farmers are, therefore, encouraged to grow non-GM soya in crop rotations, using good agricultural practices for soil and nutrient management and the best weed management options. To implement all the above-mentioned initiatives into a common practice in Croatia, farmers need to be informed of the benefits and risks of implementing all the required sustainability criteria [
1,
19].
Therefore, the main objectives of this research were: (i) to describe the response of weed communities to the above-mentioned agronomic practices in soybeans and analyze the differences in floristic composition that exist between the crops at different row spacings; (ii) to evaluate weed biomass production in competition with soybeans; (iii) to determine if soybean yield and its attributes will increase and the weed yield decrease as row spacing decreases; and (iv) to simulate a risk analysis by building models of possible results of economic return in soybeans growing in 25, 50, and 70 cm rows.
2. Materials and Methods
2.1. Site Description and Experimental Set Up
The field experiments were conducted in clay-loam soil near the city of Vukovar (45°21′ N 18°59′ E), at the family farm “Zeleno polje” in the Vukovar-Syrmia county situated in the northeastern part of the Republic of Croatia. This is an open and flat region with agriculture as the main economic sector where soybean production is important. Climatically, this region experiences a warm to moderate dry lowland climate with an average yearly temperature of 11.4oC and average yearly rainfall of 699 mm, with the highest spring regime in June.
A soybean cultivar IKA, maturity group 0-I (Agricultural Institute Osijek), was sown to achieve a population of 500,000 plants/ha
−1 to a soil depth of 4 cm on 7 May 2014 (after sunflower), on 26 April 2015 (after sugar beet), and on the 1 May 2016 (after barley). Cropping practice, typical for the local practice of soybean production in this region, consisted of primary tillage in autumn, followed by spring plowing and harrowing [
20]. Fertilizers were applied as follows: 300 kg ha
−1 (NPK 7:20:30) in fall during the primary tillage, and 350 kg/ha
−1 (NPK 15:15:15) at sowing to achieve 74 kg N ha
−1, 113 kg P
2O
5 ha
−1 and 143 kg K
2O ha
−1. No irrigation was applied, and rainfall ranged from 293.8, 357.1 to 595.4 in 2015, 2016, and 2014, respectively (
Table 1).
All treatments were arranged in a split-plot design with four replications. Two different studies were conducted during the experiment: the first one determined the weed interference, weed biomass accumulation, and yield as well as yield attributes in the soybean, and the second one simulated a risk analysis by building models of the possible results of achieving a profit as a result of weed control.
In the first study, the three main plots included different soybean row spacings: 25 cm, 50 cm, and 70 cm. Sub-plots consisted of eleven weed removal timings, where the weeds were allowed to grow until the crops reached V2 (second trifoliate), V4 (fourth trifoliate), R1 (beginning bloom, first flower), R2 (full flower), R3 (beginning pod), R4 (full pod), R5 (beginning seed), R6 (full seed), and R7 (beginning maturity). Weedy and weed-free control treatments were also included. The weeds were removed by hand pulling and hoeing. The plot size was 2.4 × 3.5 m. The plots were separated by 0.5 m and blocks by 1.5 m unplanted distances.
In the second study, as an input for the simulation model, yield losses due to weeds in soybeans were estimated from the plots where row spacing was the main plot, and different weed control treatments were subplots with four replications of each treatment. Six weed control options, including the pre-emergence herbicide metribuzin + flufenacet at 2 kg ha −1 in all of them, were evaluated for: (i) pre-emergence application in 25 cm of soybean; (ii) pre-emergence application in 50 cm of soybean; (iii) pre-emergence application in 50 cm of soybean and cultivation at V1; (iv) pre-emergence application in 70 cm of soybean; (v) pre-emergence application in 70 cm of soybean followed by cultivation at V1; (vi) pre-emergence application in 70 cm of soybean followed by cultivation at V1 and at R1.
2.2. Data Assessment and Statistical Analysis
To characterize the weed community in the soybeans, the density of each species was counted in all weedy control plots at the R2 stage (16 July, 4 July, and 6 July in 2014, 2015, and 2016, respectively) from sixteen 0.5 by 0.5 quadrats per each soybean row spacing. To overcome non-uniform weed distribution, a relative abundance was calculated [
21]. This synthetic importance value included density and frequency components and was calculated by the plot for each weed species as follows: (relative density + relative frequency)/2. The variation in the species composition was analyzed with canonical correspondence analysis (CCA) using CANOCO 5 [
22]. The statistical significance of fitting the CCA axes was tested using a global permutation test (Monte-Carlo test) of the species data at 1000 iterations. The forward selection of explanatory variables was tested with Monte-Carlo permutations and was also used in determining the statistical significance for each explanatory variable singly (simple effect) and in order for additionally explained variance (conditional effect).
Weed biomass accumulation was estimated at each soybean growth stage (described in the first study) by harvesting all the weeds within each of the four 1 m
2 randomly located treatment plots. The weeds were clipped at the soil surface and dried at 70 °C to constant moisture content. Dry weed weight was converted to a g m
2 basis. The relationship between the treatments and the weed biomass accumulation was described by using PROC REG in SAS, version 9.4 [
23]. To determine the type of relationship, a Schumacher’s model [
24] was fitted to the weed-infested treatment and weed biomass accumulation.
where Y is the weed dry weight (g m
−2), e is a constant, a is the maximum weed biomass, b is the asymptote of the curve, and x is the duration of weed infested period expressed in growing degree days (GDD). GDD was used as an explanatory variable in the regression analysis, and for that purpose, the temperature was converted to GDD by using the following equation [
25]:
where T
max and T
min are the daily maximum and minimum air temperatures (°C), and T
base is the base threshold temperature, which was set at 10 °C.
Crop yield and yield components (the number of pods per plant and 1000 kernel weight) were recorded from each plot by hand harvesting the two middle rows (on 1 October, 30 September, and 3 October in 2014, 2015, and 2016, respectively), shelled and adjusted to 11% moisture. A mixed model (PROC MIXED in SAS) was used to evaluate the interference of the weeds in soybeans planted at 25, 50, and 70 cm row spacing on the soybean yield and yield components. The analysis involved three factors (row spacing, weed interference, and year) and repeated measures of the ANOVA model, with year as a repeated measure. Significance was assumed at p < 0.05.
A Monte-Carlo simulation model was constructed to forecast the distribution of the difference in profit based on the data from study 2. For that purpose, data for the seeds, fertilizers, and herbicide prices were obtained from the local seed and agrichemical dealers. Prices for fuel, services, and soybean markets were obtained from the TISUP [
26]. Weed management costs were the sum of the herbicide and their application cost. The gross margin for the treatments was determined and represented the difference between the gross receipt (product of crop yield and assumed market price) and production costs (seeds, fertilizers, fuel, services, and weed management).
By using a Monte-Carlo simulation model [
27], the inputs are not simply mean values of the estimated parameters but show the variability of these estimates by using their simulated distributions. Calculations were performed by using a Risk Solver
® platform in Excel.
3. Results
3.1. Weed Community Characteristics
The weed community comprised 34 grass and dicot species of varied perennation (
Table 2). They belong to 30 genera and 19 families, with
Asteraceae (seven weed species) and
Poaceae (with four weed species) as the leading families. The weed community was composed of twenty-six species in 2014, eighteen species in 2015, and fifteen species in 2016. The dominant weeds during this study were
Amaranthus retroflexus L.,
Ambrosia artemisiifolia L.,
Chenopodium album L.,
Datura stramonium L., and
Setaria viridis (L.) PB., and
Sorghum halepense (L.) Pers.
The weed relative abundance data over all the study years are presented in
Table 2 to give an overview of the community structure and the overall effects of the various row spacing in soybeans. However, in the analysis of the comprehensive data using CCA, specific weed species’ responses to the explored external variables proved to be statistically significant (
Figure 1). A Monte-Carlo permutation test showed both the first and all the CCA axes together to be statistically significant (test of significance of first canonical axis: eigenvalue = 0.1667; F-ratio = 30.2908;
p < 0.001; test of significance of all canonical axes: eigenvalue = 0.3195; F-ratio = 12.6323;
p < 0.001).
The first CCA axis captured 52.9% of the variation in the species composition and was explained more by the seasonal aspect than management factors. The biplot score for 2015 and 2016 had a longer vector length and was in opposite orientation in the ordination than for 2014 (
Figure 1). During the growing seasons of 2015 and 2016 (April–September), the amount of rainfall (293.8 mm and 357.1 mm for 2015 and 2016, respectively) was significantly lower than in 2014 (595.4 mm). A weaker association were observed for the crop rows. The soybeans planted in narrow rows (25 cm) were located in the same orientation space in the driest season in 2015, compared to a wider row spacing with less seasonal influence.
The weed relative abundance indicated that species, such as
Abutilon theophrasti,
Setaria verticilata,
Euphorbia helioscopia,
A. retroflexus, and
D. stramonium, are associated with a drier environment in soybean crops growing in 25 cm rows (
Table 3).
The second axis explains a further 29.4% of the species composition variability, and
Rumex crispus,
Rorippa sylvestris, and
Hordeum murinum appeared only in 2016 in wider crop rows (
Table 3).
The results from the forward selection of explanatory variables indicate that most variance in the species data, when examined singly (simple effects), can be explained by the weather conditions in very humid 2014 and very arid 2015, and moderate 2016 (
Table 4). Next, the variation explained by the environmental variables are in order of their inclusion in the model, i.e., the conditional effect indicated that only 2014 and 2015 explain significant portions of the variation in the data since the additional variance was explained by each variable at the time it was included.
3.2. Weed Biomass Accumulation
Weed biomass accumulation increased in each row spacing as the duration of the weed-infested period increased (
Figure 2). The highest total dry weed biomass (3329.2 g m
−2) was recorded in 2016 in soybeans growing in 70 cm rows, following 2015 (3082.4 g m
−2) where soybeans grew in 25 cm rows and 2014 (2521.7 g m
−2) where soybeans grew in 50 cm rows.
There was no reduction in the weed biomass accumulation with a reduction in the row spacing from 70 to 50 and 25 cm. Dominant weed species
A. retroflexus,
A. artemisiifolia,
C. album,
D. stramonium,
S. viridis, and
S. halepense germinated and emerged with the soybeans, formed the main biomass, and were spread over all the row spacings (
Table 2).
Figure 2.
The effect of increasing duration of weed interference on weed dry weight accumulation in soybeans growing in 25, 50, and 70 cm rows. Dots indicate observed data. Parameter values for response curves calculated using Schumacher’s model are presented in
Table 5.
Figure 2.
The effect of increasing duration of weed interference on weed dry weight accumulation in soybeans growing in 25, 50, and 70 cm rows. Dots indicate observed data. Parameter values for response curves calculated using Schumacher’s model are presented in
Table 5.
3.3. Weed Community Interference on Soybean Yield and Yield Components
Weed competition with soybeans can cause significant reductions in yield potential (
Table 6). The soybean yield, number of pods per plant, and 1000 kernel weight were significantly influenced by row spacing (
p < 0.001). The duration of weed interference affected soybean yield and its components (
p < 0.001). Significant yearly differences (
p < 0.001) were also recorded and can be attributed to the rainfall amount periodicity. The first and third years of the experiment received more rainfall during the vegetative growth compared to the second experimental year (
Table 1).
The significantly highest soybean yield was observed in weed-free crops growing in 70 cm rows, while in 50 cm rows, the yield was reduced by 10%, and in 25 cm rows, the yield reduction was 20% (
Table 7). However, when weeds were present throughout the growing season, the yield reductions were 52, 53, and 60% in 50 cm, 25 cm, and 70 cm soybean rows, respectively.
The number of pods per plant in weed-free soybeans did not significantly differ but allowed weeds to interfere the whole season, causing a dramatic reduction (70%) in the number of pods per plant to be observed in the 70 cm rows. The 1000 kernel weight was also not significantly different among the crop rows in weed-free plots. The presence of weeds reduced the 1000 kernel weight to about 10%.
Weed infestation until the second trifoliate (V2) stage had no detrimental effect on the soybean yield, regardless of the row spacing. The number of pods per plant started to significantly decrease at the same stage in 25 and 50-cm rows, but in 70-cm wide soybean rows, this process started later, at four unfolded trifoliate leaves (V4 stage). A 1000 kernel weight was less sensitive to weed infestation and significantly decreased at full flowering (R2 stage) in the 25 and 70cm rows, while soybeans growing in 50 cm rows decreased at four unfolded trifoliate leaves (V4 stage).
3.4. Effect of Weed Management in Soybean on Economic Return
Following the experimental results (the second study), a Monte-Carlo simulation was constructed to forecast the distribution of the difference in obtaining the profit between the six weed control strategies in soybeans growing in 25, 50, and 70 cm rows.
Figure 3 shows a flow diagram representing the sequence of the calculations performed. The differences in profits between the weed control strategies were presented for each study year separately. However, in the second year (very dry growing season), all examined weed management strategies failed, and variable costs (dotted line in
Figure 3) were higher than the gross margin.
A pre-emergence herbicide application was not able to obtain positive financial results in any of the three years of study. In 50 cm rows, the probability of achieving a profit was 20% in 2014 and 40% in 2016, while in 70 cm rows, the positive financial results with a pre-emergence herbicide application had a 60% probability in 2014 and 30% in 2016.
The combination of the pre-emergence herbicide and one inter-row cultivation at V1 in 50 cm rows only received a positive financial result in 2014 (a 20% probability). Soybeans growing in 70 cm rows with one inter-row cultivation and with two inter-row cultivation were the best options for achieving profit. There was a 50% and 70% possibility for positive financial results with pre-emergence herbicides and one inter-row cultivation in 2014 and 2016, respectively. The best option in this study was the strategy with a pre-emergence herbicide application followed by two inter-row cultivations. In this treatment, the probability of not receiving a profit was only 10% in 2014.
4. Discussion
The spontaneous weed community that developed in the soybean crops during this study was typical of row crops flora in this region [
28,
29]. The dominant weeds that appear in the study are, likewise, a major limiting factor of optimum soybean production in the region [
30,
31] as well as worldwide [
32,
33,
34].
It is well documented that management practices and environmental conditions affect soybean and weed competition. Weeds impact soybeans by competing for limited light, water, nutrient resources, and space. The manipulation of row spacing, as one of the management tactics, can have a sizable impact on weeds in soybeans [
35]. However, our results suggest that environmental factors, followed by crop management practices, had the most significant effects on weed community composition and structure and on soybeans. The variation in weed community composition was significantly explained by seasonal conditions. The amount of precipitation during the growing season was the most important factor determining weed species composition. Similar conclusions can be found in other studies [
36,
37,
38,
39], where year (as a factor) was found to show greater influence on the weed community than soil management practices, indicating fluctuational rather than directional or consistent changes in the weed community. However, some other studies suggest that management factors play more important roles than environmental ones; [
40,
41] revealed that different crop types and their associated management have more influence on weed composition than the relative importance of climatic variables.
The weed biomass and accumulation rate increased with the increasing duration of weed infestation. The absence of herbicide application or other weed control measures in these weeded plots gave the weeds competitive advantages over the crops growing in all three-row spacings. Since dominant weeds
A. retroflexus,
A. artemisiifolia,
C. album,
D. stramonium,
S. viridis, and
S. halepense were spread over all the row spacings, they formed the main biomass because of their size. Moreover, they emerged with the crop resulting in vigorous competition for a weak competitor as soybeans are in their early growth stages [
42]. It is well documented that weeds that germinate and emerge with soybeans are usually better competitors for light, water, and nutrients [
43,
44]. Broadleaf weeds, in particular, such as
A. retroflexus,
A. artemisiifolia,
C. album, and
D. stramonium, are more competitive with soybeans than grass weeds and many late emerging weeds [
45]. Moreover, soybeans are very sensitive to moisture deficiencies in late summer, and even a few large weeds left in the field can severely reduce the yield potential [
46].
Weed infestation is one of the main causes of low soybean yield [
47]. Narrow row spacing is a cultural practice that has been reported as a management tactic which reduces the amount of light that reaches the soil surface and that reduces the amount of time needed for the soybeans to reach a full canopy closure [
48,
49]. A number of studies showed a yield increase when soybeans were planted in narrow rows [
35,
50]. Our results are in contrast to reports where narrow row spacing improved weed control efficacy, productivity, and the profitability of the soybean. The reason for this lay in the fact that we did not manipulate with seeding density (the soybeans were sown to achieve the recommended population of 500,000 plants per hectare).
In a given crop population, soybeans planted in 70 cm spaced rows gave the best economic results, according to the Monte Carlo simulation. The probability distribution of achieving a profit showed the best results for soybeans growing in 70 cm rows, with preemergence herbicide applications and two inter-row cultivations (at soybean stages V1 and R1). Individual weed control methods generally do not provide complete control of the weeds [
51], and therefore, integrated weed management is often considered to be the most effective approach, as was confirmed by this research as well. Pre-emergence herbicide application followed by mechanical cultivation can increase the soybean yield and receive a positive financial result. Harder et al. [
6] reported that gross margins were usually greater in 19 and 38-cm soybean rows, but in 76-cm rows, at a low soybean population, the gross margin was the greatest. The combination of pre-emergence herbicide application and mechanical weed control is not a fixed process and needs to be adjusted depending on the crop type, farming operation, and seasonal conditions [
52].