1. Introduction
The Crop legumes are an essential component of sustainable agriculture, given the ecological services they provide, especially their ability to fix atmospheric nitrogen [
1]. However, legume cultivation is hampered by biotic and abiotic stresses, which can seriously affect their yields [
2]. One of these is the competition of weeds, which is a serious challenge for managing these crops. Not only may weeds reduce crop yields, but they also may affect the quality of the grain [
3], serve as alternative hosts for diseases and pests [
4,
5], or even hinder harvesting. In general, weed control methods must be implemented for legume crops since these are inefficient competitors [
6]. Herbicides are the most widely used weed management method, but these present several drawbacks, such as high costs for the farmer, environmental damage, and the eventual appearance of resistance to the active ingredients in weeds. These resistances may be especially problematic, as no herbicides with new modes of action have appeared in recent years [
7]. Other strategies are mechanical control and management of the spatial arrangement of crops [
8,
9].
Faba bean (
Vicia faba L.) stands out as one of the most cultivated legumes worldwide for its use for both animal feed and human consumption [
10]. Weeds are one of the major biological constraints for faba bean, reducing yields by up to 50%, although they appear more competitive with weeds than other legumes [
11].
Intercropping is the practice of growing two or more crop species simultaneously in one field. It may help in the control of weeds in many crops [
12], complementing herbicides if necessary [
13]. In the case of legumes, weed reduction has been reported when combined with diverse crops, such as corn, wheat, oat, barley or rape [
14,
15,
16,
17,
18,
19]. There are different types of intercropping, depending on the arrangement and sowing densities of the combined crops; mixed intercropping, when plants are completely mixed; and alternate intercropping, when rows of each crop alternate regularly [
20]; in addition intercropping all crops are sown at their normal densities, while in replacement intercropping densities are adjusted for the presence of the other(s) crop(s) [
21].
Different mechanisms explain the smothering of weeds in intercropping systems. The most common is the competition of crops and weeds for natural resources such as light, water or nitrogen [
18]. Another one is allelopathy, which is a characteristic of plants that may be exploited for weed control [
22]. Allelopathy is defined as the capability of plants to exert positive or negative influence in the surrounding area by releasing chemicals [
23]. It has been known for a long time that several crops present allelopathic activity, such as oat, wheat, barley, and sorghum [
24].
The main objective of this work was to establish if it is possible to control weeds in faba bean by intercropping, determining the best crop combination if so. To this effect, field trials were performed across different years. Additionally, experiments under controlled conditions were carried out to investigate the possible role of allelopathy in weed control in intercrops of faba bean.
2. Materials and Methods
2.1. Field trials
Five field trials were carried out from 2014 to 2018 in Córdoba (South of Spain) to assess the effect of different crop combinations on the control of weeds in faba bean (
Table 1). The combinations tested were faba bean/wheat, faba bean/barley, and faba bean/pea as well as the monocrops of the four crops, and a monocrop of faba bean at 50% of sowing density (doubling the distance between rows) to check the effect of density on weed infestation. The cultivars used were 'Muchamiel' (faba bean), 'Califa' (wheat), 'Henley' (barley), and 'Messire' (pea). Sowing densities were 80 seeds/m
2 for legumes and 200 seeds/m
2 for cereals. Two different intercropping systems were evaluated: alternate with replacement at 50%, where rows of each crop are alternated, giving a final sowing rate of 50/50, i.e., half of the normal density of each crop; and alternate with addition, where rows are also alternated, but to a final rate of 100/100, i.e., the normal density of each crop. Replacement intercropping was evaluated in trials Córdoba-15-r and Córdoba-16-r, and addition intercropping was evaluated in trials Córdoba-16-a, Córdoba-17-a, and Córdoba-18-a. The experimental plots had a length of 3 m and comprised eight rows at 35 cm distance between them in the case of monocrops and replacement intercropping and 16 rows at 17.5 cm in the case of addition intercropping; the monocrop of faba bean at 50% had four rows at 70 cm per plot (
Figure 1). For each trial the experiment was designed as a randomized complete block with four replications.
Two areas of 0.7 m2 were randomly evaluated in the central part of each plot at the faba bean maturity stage to assess weed infestation. The weeds were counted and classified by botanic families. Two diversification indices were calculated in trials Córdoba-16-a, Córdoba-17-a, and Córdoba-18-a:
Richness (R) and
Relative density (D):
R: Richness; nF
i: number of weed families present in a given plot (i); nF
T: total number of weed families found in the whole trial
Dij: Relative density of weed family j in plot i; Nji: number of plants of weed family j present in plot i; NTi: total number of weed plants in plot i.
Additionally, two evaluators visually estimated the weed coverage as a percentage of the area covered by weeds. At the end of the experiment, in trials Córdoba-16-a and Córdoba-17-a, the weeds in a central area of 2 m2 in each plot were harvested, and their biomass determined by drying them in an oven at 60 ºC for three days and subsequently weighing them.
2.2. Controlled-Condition Experiments
Two experiments were carried out under controlled conditions to investigate a possible allelopathic effect of barley. Thus, four weeds that had grown in the field trials were selected for this: Polypogon monspeliensis (L.) Desf., Matricaria camomilla L., Sinapis arvensis L., and Medicago truncatula Gaertm, belonging to the botanical families Poaceae, Asteraceae, Cruciferae and Fabaceae, respectively. In the first one, 10 seeds of barley were sown on plastic pots (5× 5 × 10 cm) filled with a mixture of sand and peat at 1:3 rate (v:v) and grown in a growth chamber with a photoperiod of 12 h of visible light (150 μmol m−2 s−1 photon flux density) at 25 ºC, and 12 h of darkness at 20 ºC. Two weeks later, the barley plants were removed, and 20 seeds of the same weed species were sown in the same pots maintaining the substrate. Weed plants were counted every two weeks for 42 days. Then, weed plants were removed, dried in an oven at 60ºC for three days, and weighed to obtain their biomass. A treatment control was used, in which weeds were sown in pots where barley has not been previously grown. All this was done in separate pots for each one of the four weeds. The second experiment was carried out in the same way as in the first one, but barley plants remained in the pots when weeds were sown, and they grew together till the end of the experiment. The design in the two cases was randomized with seven replications. Each of these experiments was performed twice.
2.3. Statistical Analysis
Considering the type of data, generalized linear models (GLMs) were performed to study the influence of the treatment on weed coverage and biomass in field trials, following the next model:
Where
Y is the dependent variable,
βo,
β1 and
β3 denote the estimated parameters (weights), and
Ɛ is the residual. Model fits were evaluated through the residual plots. Subsequently, we assessed the contribution of the treatment, year, and treatment×year according to Cohen’s f factor [
25].
Given that the interactions treatment×year were significant (Cohen’s f > 0.4) for both coverage and biomass, and since the study's main objective was to compare the intercropping treatments independently of the year×crop interaction, a Fridman nonparametric two-way analysis of variance were used. Subsequently, Dunn’s test were performed to compare the treatment mean ranks at P ≤ 0.05. Additional ANOVA analyses were carried out to study the effect of the treatment on the weed richness and weed relative density in each experiment.
In the experiments under controlled conditions, the Poisson regression procedure was applied using the maximum likelihood estimation method to study how the plant counts depended on the presence or absence of barley plants. Thus, data were fitted to the following equation:
Where Y is the number of plants, βo is the independent term, β1 and β2 are the estimated associated with the treatment (presence or not of barley). Difference between treatments was evaluated according to the p-value of the Poisson regression. Additionally, a t-test was performed to calculate the barley presence's effect on the weed species' biomass and compare the proportion of reduction of weed plants and weed biomass between the two types of experiments. Finally, the proportion of weed plants and biomass reduction between weed species was assessed by ANOVA analysis.
Data were analyzed using Statistix software (Version 10; Statistix, Tallahassee, FL) and the R statistical software package.
4. Discussion
In this work, we have evaluated intercropping as a tool to control weeds in faba bean in the South of Spain. The level of control attained was very high, reaching reductions of 92.7% in weed coverage when combined with barley and 76.6% and 46.1% in weed biomass when mixed with barley and wheat, respectively. The mixture of faba bean and barley, then, might be effective enough to dispense with the application of herbicides in this agroecosystem.
Previously only one work had studied the effect of the combination of faba bean and barley on weeds: Dhima et al. [
26] reported that alternate replacement intercropping achieved levels of reduction of corn poppy (
Papaver rhoeas L.) of around 90%. As for other cereals, the faba bean and wheat mixture has proven to reduce weed biomass by 60% in alternate and mixed replacement intercropping [
27]. Equally, Boutagayout et al. [
28] found that an alternate intercrop of faba bean with wheat and oat decreased the weed biomass by around 50%. As far as we know, our work is the first to achieve such a level of weed supression in faba bean, validated across three field experiments. Equally, barley has proved to be successful in reducing broomrape (
Orobanche crenata Forsk.) infestation in faba bean [
29], and also to be effective in controlling rust disease in faba bean [
30]. In these cases, the intercropping systems were different from those assessed in this work, so it would be necessary to integrate them in the best way to maximize the benefits of faba bean / barley intercropping. Furthermore, it is remarkable that the combination of faba bean and wheat reduced 64% weed biomass compared to the wheat monocrop, making this mixture an exciting option for farmers that can be combined with some herbicide applications. Pea, however, has proved to be a poor competitor against weeds, both as monocrop and in mixture with faba bean, which is in line with the fact that the combination of grain legumes alone is not a good strategy to control weeds. Actually, these combinations are not found in the literature dealing with weed supression by intercropping [
31].
We did not find differences in weed diversity among all the tested treatments: no weeds belonging to a particular botanic family were more affected by any crop combination or monocrops than others. This is in contrast with what has been found in other intercropping systems [
32,
33], although there are also situations where weed diversity has not been influenced by intercropping [
33]. In our case, the weed community has proved to be relatively stable regardless of the crop or crops present.
Of the two types of intercropping we tested, alternate with replacement and alternate with addition, only the latter was adequate for controlling weeds. The relationship between crops and weeds is based on the competition for available resources such as water, nutrients or light [
34,
35]. Plant density is one key factor in the improvement of crop competitiveness against weeds [
36], and that is precisely the difference between both intercropping systems. In addition, intercropping plant density is doubled in comparison with alternate intercropping. The high weed pressure levels that we observed in the plots with faba bean sown at half density confirm the importance of plant density.
Plant density, however, is not the only mechanism that explains weed suppression in intercropping. If that were the case, we would have observed similar results with barley, wheat, or pea. Weed pressure levels for them as sole crops illustrate that not all crops have the same competitive ability against weeds: barley presents very low weed infestation compared to the other two, with pea ranking the highest. Other factors influence the performance in the presence of weeds, such as plant architecture, vigor, or allelopathy [
37,
38].
Barley has been described as one of the most competitive crops against weeds by different authors [
39,
40]. The rapid biomass accumulation and high growth rates that barley shows at the beginning of its cycle are some of the reasons for this [
40]. Another reason is the efficiency of barley in taking nitrogen: it has been reported to be more competitive for nitrogen than pea in intercrops, so depriving weeds of this nutrient. Beyond that, barley is considered a crop with high levels of allelopathy [
41], and as many as 44 potential allelochemicals have been identified so far [
42]. The two more important are the alkaloids Gramine and Hordeine, which appear in barley plants' leaves, roots, and roots exudates [
23,
42]. All this has made barley a common partner in crop diversification for weed control [
31].
The experiments under controlled conditions aimed at evaluating the role that barley allelopathy might have on our results. Different types of bioassays under controlled conditions may be used to assess the allelopathic ability of a plant species, such as testing extracts from the allelopathic plant [
43,
44], agar bioassays [
39] or pots screenings [
45], where plants are grown together in Petri dishes or pots with soil, respectively. We opted for pot screening because it may better reflect the conditions under which allelopathy operates. Besides, our design, in which removing barley plants is compared with the effect of not removing them, allows for discriminating allelopathy from competition effects.
The results from the pots where barley was removed before sowing the weeds show a substantial allelopathic effect against them. All four weed species presented a decrease in plant emergence and biomass that can only be explained by the presence in the soil of chemical compounds previously released by the barley plants. These results also suggest that the main allelopathic effect is related to the first stages of seed germination and seedling development. Remarkably, there were no significant differences between the final number of weed plants in the pots where barley had been removed and in those where it remained till the end. As for weed biomass, however, in the case of two species (
Sinapis arvensis and
Matricaria camomilla) the decrease was higher when barley remained than when it was removed. This is probably due to additional allelopathic effects, although it is more difficult to separate them from competition effects in this case. Barley allelopathy had been previously tested on
Sinapis arvensis [
46,
47]. Still, as far as we know, this is the first time it has been assessed on the other three weeds.
The fact that there has been no difference in the reduction of emerged plants and biomass between four weeds belonging to such different botanical families points to a global and non-discriminatory effect of barley in our case. These weeds are a sample of the ten families found in our area, and this global effect could explain the lack of differences in the composition of weed communities between the barley intercrops and the monocrops in our field experiments.
In conclusion, this research has established that the combination of faba bean and barley has great potential for the control of weeds in the agroclimatic region of the South of Spain. Further work should focus on developing the optimal intercropping strategy to obtain the best advantages from combining these two crops. It would also be of great interest to identify the allelochemicals behind weed suppression in our experiments and determine their mechanisms of action.
Figure 1.
Different cropping systems evaluated in the field trials. A: Faba bean monocrop. B: Alternate replacement intercropping. C: Alternate addition intercropping. D: Faba bean monocrop at 50% of sowing density. Row distances in A and B were 35 cm., in C were 17.5 cm. and in D 70 cm.
Figure 1.
Different cropping systems evaluated in the field trials. A: Faba bean monocrop. B: Alternate replacement intercropping. C: Alternate addition intercropping. D: Faba bean monocrop at 50% of sowing density. Row distances in A and B were 35 cm., in C were 17.5 cm. and in D 70 cm.
Figure 2.
Weed coverage in the different treatments evaluated across the three intercropping field trials where addition alternate intercropping was tested. Different letters mean significant differences according to Dunn’s test at p ≤ 0.05.
Figure 2.
Weed coverage in the different treatments evaluated across the three intercropping field trials where addition alternate intercropping was tested. Different letters mean significant differences according to Dunn’s test at p ≤ 0.05.
Figure 3.
Weed biomass (g/m2) in the different treatments evaluated across the two intercropping field trials where addition alternate intercropping was tested. Different letters mean significant differences according to Least Significant Difference (LSD) test at p ≤ 0.05.
Figure 3.
Weed biomass (g/m2) in the different treatments evaluated across the two intercropping field trials where addition alternate intercropping was tested. Different letters mean significant differences according to Least Significant Difference (LSD) test at p ≤ 0.05.
Table 1.
Field trials carried out for the study of control of weeds in faba bean by intercropping.
Table 1.
Field trials carried out for the study of control of weeds in faba bean by intercropping.
Trial |
Córdoba-15-r |
Córdoba-16-r |
Córdoba-16-a |
Córdoba-17-a |
Córdoba-18-a |
Season |
2014/15 |
2015/16 |
2015/16 |
2016/17 |
2017/18 |
Intercropping system |
Replacement |
Replacement |
Addition |
Addition |
Addition |
Maximum T (ºC) |
35.3 |
31.7 |
31.7 |
32.2 |
29.7 |
Minimum T (ºC) |
−3.3 |
-2.3 |
-2.3 |
−3.4 |
-3.4 |
Mean T (ºC) |
12.3 |
12.6 |
12.6 |
14.3 |
12.0 |
Precipitation (ml) |
150 |
336 |
336 |
143.8 |
444.6 |
Table 2.
Presence of the different botanical families to which the weeds identified belonged in the field trials Córdoba-16, Córdoba-17 and Córdoba-18.
Table 2.
Presence of the different botanical families to which the weeds identified belonged in the field trials Córdoba-16, Córdoba-17 and Córdoba-18.
|
Córdoba-16 |
Córdoba-17 |
Córdoba-18 |
Amarantaceae |
|
x |
x |
Asteraceae |
x |
|
|
Convolvulaceae |
x |
x |
x |
Cruciferae |
x |
|
x |
Cyperaceae |
|
x |
|
Fabaceae |
x |
x |
x |
Papaveraceae |
|
x |
x |
Poaceae |
x |
|
x |
Polygonaceae |
x |
x |
x |
Quenopodiaceae |
x |
|
|
Table 3.
Weed richness (R) for the different treatments tested in three field trials (Córdoba-16, Córdoba-17, and Córdoba-18). Different letters mean significant differences for each trial (Tukey test, p <0.05). Richness was not measured for barley, wheat and pea monocrops in Cordoba-16.
Table 3.
Weed richness (R) for the different treatments tested in three field trials (Córdoba-16, Córdoba-17, and Córdoba-18). Different letters mean significant differences for each trial (Tukey test, p <0.05). Richness was not measured for barley, wheat and pea monocrops in Cordoba-16.
|
Richness (R) |
|
Córdoba-16 |
Córdoba-17 |
Córdoba-18 |
Faba bean 100% |
35.7 a |
78.6 ab |
76.2 ab |
Faba bean / barley |
25.0 a |
50.0 ab |
42.3 a |
Faba bean /wheat |
25.0 a |
46.4 a |
61.9 ab |
Faba bean / pea |
50.0 a |
64.3 ab |
85.7 b |
Faba bean 50% |
46.4 a |
75.0 ab |
85.7 b |
Barley 100% |
|
83.3 b |
76.2 ab |
Wheat 100% |
|
75.0 ab |
95.0 b |
Pea 100% |
|
83.3 b |
80.9 ab |
Table 4.
Number of weed plants per pot for each treatment in the two experiments under controlled conditions: where barley was removed before sowing the weeds and where barley was not removed after sowing the weeds. Different letters mean significant differences for each weed in each type of experiment, i.e., with barley removed or not removed. Differences were significant in all cases according to Poisson’s p-value at p ≤ 0.05.
Table 4.
Number of weed plants per pot for each treatment in the two experiments under controlled conditions: where barley was removed before sowing the weeds and where barley was not removed after sowing the weeds. Different letters mean significant differences for each weed in each type of experiment, i.e., with barley removed or not removed. Differences were significant in all cases according to Poisson’s p-value at p ≤ 0.05.
|
Barley removed |
Barley not removed |
|
Control |
Barley sown |
Control |
Barley sown |
Polypogon monspeliensis |
19.9 a |
5.0 b |
19.7 a |
4.6 b |
Matricaria chamomilla |
18.9 a |
5.9 b |
18.6 a |
5.1 b |
Sinapis arvensis |
11.1 a |
3.8 b |
11.9 a |
2.7 b |
Medicago truncatula |
5.2 a |
1.4 b |
7.8 a |
1.4 b |
Table 5.
Biomass of weed plants per pot (mg) for each treatment in the two experiments under controlled conditions: where barley was removed before sowing the weeds and where barley was not removed after sowing the weeds. Different letters mean significant differences for each weed in each type of experiment, i.e., with barley removed or not removed, according to the t-student test at p ≤ 0.05.
Table 5.
Biomass of weed plants per pot (mg) for each treatment in the two experiments under controlled conditions: where barley was removed before sowing the weeds and where barley was not removed after sowing the weeds. Different letters mean significant differences for each weed in each type of experiment, i.e., with barley removed or not removed, according to the t-student test at p ≤ 0.05.
|
Barley removed |
Barley not removed |
|
Control |
Barley sown |
Control |
Barley sown |
Polypogon monspeliensis |
39.9 a |
20.2 b |
55.5 a |
24.9 b |
Matricaria chamomilla |
115.4 a |
49.0 b |
131.6 a |
25.6 b |
Sinapis arvensis |
70.6 a |
41.7 b |
104.6 a |
10.1 b |
Medicago truncatula |
16.3 a |
10.6 b |
35.0 a |
4.8 b |
Table 6.
Proportion of number of weed plants and biomass per pot as compared to the control for both types of experiments under controlled conditions (removing barley and not removing it just prior to sowing the weeds). Different letters mean significant differences for each weed and parameter: number and biomass of weeds (according to the t-student test at p ≤ 0.05).
Table 6.
Proportion of number of weed plants and biomass per pot as compared to the control for both types of experiments under controlled conditions (removing barley and not removing it just prior to sowing the weeds). Different letters mean significant differences for each weed and parameter: number and biomass of weeds (according to the t-student test at p ≤ 0.05).
|
No of weeds (%) |
Biomass of weeds (%) |
|
Barley removed |
Barley not removed |
Barley removed |
Barley not removed |
Polypogon monspeliensis |
25.2 a |
23.6 a |
70.1 a |
61.9 a |
Matricaria chamomilla |
31.1 a |
27.3 a |
65.5 a |
43.7 b |
Sinapis arvensis |
32.1 a |
23.0 a |
77.0 a |
29.3 b |
Medicago truncatula |
24.9 a |
18.4 a |
58.9 a |
31.9 a |