The Importance of Species Selection in Cover Crop Mixture Design

Abstract Cover crops are increasingly being included in crop rotations as a mechanism to promote diversity and provide agroecosystem services, including weed suppression. Recently, cover crop mixtures have increased in popularity in an attempt to provide a greater diversity in ecological services as compared with monocultures. Several recent studies, however, have failed to detect a positive effect of cover crop diversity on biomass production or weed suppression. Here we assessed biomass productivity and weed suppression in 19 cover crops seeded as monocultures and 19 mixtures of varying species composition and functional richness (two- and three-species mixtures) of full-season cover crops in Atlantic Canada. Cover crop biomass production and weed suppression varied by species identity, functional diversity, and species richness. As cover crop biomass increased regardless of diversity, weed biomass declined. Highly productive forbs and grasses provided the greatest weed suppression in monoculture. In line with previous observations, mixtures were not more productive or weed suppressive on average than the most productive monocultures. We observed that the inclusion of the highly productive species buckwheat (Fagopyrum esculentum Moench) and sorghum–sudangrass [Sorghum × drummondi (Nees ex Steud.) Millsp. & Chase] in a mixture increased stand evenness, productivity, weed suppression, and spatiotemporal stability. Taken together, our results suggest that effects of diversity on mixture productivity and weed suppression are species specific. This further demonstrates the importance of species selection in cover crop mixture design.


Introduction
Cover crops are increasingly being adopted to combat declines in agroecosystem biodiversity and health and as a tool for sustainable agricultural intensification through provision of agroecosystem services such as nutrient retention and scavenging, improving soil organic matter, and reducing erosion (Blanco-Canqui et al. 2015). In addition, the use of cover crops can be an important component of integrated weed management programs (Teasdale 1996) and as a herbicide-resistance management tool (Norsworthy et al. 2012). Cover crops fulfilling specific ecological goals are traditionally sown in monoculture, as species vary in functional traits such as rooting density, rate of leaf appearance, canopy development, and weed suppression, influencing agroecosystem processes (Wagg et al. 2021). Recently, diverse species mixtures have gained popularity due to the suggestion of enhanced productivity, synergism in agroecosystem services and improved year-to-year stability as compared with monocultures (Finney and Kaye 2017). Therefore, it is important to understand how increasing cover crop diversity from a range of species can improve productivity across diverse environments.
Cover crops can suppress weeds through a variety of mechanisms. Living cover crops directly compete with weeds for light, water, nutrients, and space, reducing weed germination, growth, and reproductive potential (Christina et al. 2021). The capacity of cover crops to suppress weeds, however, is species and taxon specific through differences in phenotypic traits such as canopy development and relative growth rate. Plant functional groups such as grasses, including cereal rye (Secale cereale L.) and oat (Avena sativa L.), as well as forbs such as buckwheat (Fagopyrum esculentum Moench), have rapid growth and canopy development to suppress weeds and prevent their establishment (Baraibar et al. 2018). Similarly, living Brassica spp. have been demonstrated to suppress weeds through canopy light interception (Lawley et al. 2012). In contrast, some legume monocultures are often cited as poor annual cover crops for weed suppression due to their slow establishment and growth (Elsalahy et al. 2019). Taken together, cover crop biomass production appears to be a key metric and mechanism of weed suppression.
Cover crop mixtures are purported to increase productivity and agroecosystem services through greater complementarity in resource partitioning (Finney and Kaye 2017). Therefore, it stands to reason that cover crop mixtures will provide greater weed suppression than monocultures. Increased diversity driving enhanced productivity and agroecosystem services has been demonstrated in a variety of experimental grassland systems (Loreau and Hector 2001). Despite this, numerous agricultural studies have failed to observe enhanced productivity or weed suppression by diverse cover crop mixtures (Baraibar et al. 2018;Florence et al. 2019;Florence and McGuire 2020;Smith et al. 2014Smith et al. , 2020. It has been suggested that the inclusion of a highly productive grass species in a mixture is primarily responsible for improved weed suppression. Other highly productive species such as buckwheat or crimson clover (Trifolium incarnatum L.), however, may provide similar weed suppression as grass species when included in mixture and provide additional agroecosystem services. Diverse cover crop mixtures that include multiple highly productive species may provide additional benefits, such as improved stability of cover crop productivity and weed suppression, providing insurance against widely fluctuating environments (Loreau et al. 2021).
Within Atlantic Canada, stagnant yields and declines in soil health have been a driving factor for the adoption of cover crops (Mann et al. 2021;Nyiraneza et al. 2017). Due to the region's long cool summer growing season and late harvest, there is limited opportunity for fall cover crop establishment, and soil is often left bare over the winter. As such, cover crops are typically grown in the year before or following a main cash crop, often potato (Solanum tuberosum L.), soybean [Glycine max (L.) Merr.], or small grains. In this system, cover crops are seeded in late June or early July and mowed in the fall, with residues remaining on the soil surface to overwinter until spring. In this study, we evaluated cover crop biomass productivity and weed suppression of a range of cover crops grown for a full growing season from various functional groups seeded in monoculture and in various two-and three-way species mixtures in Atlantic Canada. This study is part of a larger project aimed at evaluating cover crop suitability for adoption in Atlantic Canadian cropping systems. We developed three hypotheses to evaluate productivity and weed suppression by annual cover crops grown for a full growing season in monoculture and mixture. We hypothesized that (1) greater cover crop biomass productivity would provide greater weed suppression; (2) diverse cover crop mixtures would provide equivalent biomass productivity and weed suppression as compared with monocultures; and (3) diverse cover crop mixtures would be more stable in terms of biomass productivity and weed suppression than monocultures.

Experimental Design
Field experiments were conducted at the Charlottetown Research and Development Centre's experimental farm (Harrington, PE, Canada; 46.35°N, 63.15°W) for 4 yr from 2018 to 2021. Soil at Harrington is classified as a Orthic Humo-Ferric Podzol (Canadian System of Soil Classification; Humic Cryorthods, U.S. Soil Taxonomy; pH 6.5 and organic matter 3%). In Atlantic Canada, cover crops are sown in late June or early July and mowed in September. Monthly temperatures at Harrington, PE, during this period average a low of 13.1 C and high of 21.0 C; monthly average precipitation was 76.6 mm over the previous 10 yr (2011 to 2021; Supplemental Table 1). Fields varied across years, with the preceding crops being soybean in 2018 and 2019 and buckwheat double cropped as a cover crop in 2020 and 2021.
A total of 38 cover crop treatments were selected; 19 cover crops sown as monocultures and 19 cover crop mixtures (Table 1). Cover crop treatments ranged in species and functional richness from one to three species. In addition, a bare soil control was included for a total of 39 treatments. Individual cover crops were chosen from those available in the region, and which provide a range of ecological services (Clark 2012) representing four taxonomic groups (brassica, forb, grass, and legume). Two-and three-species mixtures were composed of one, two, or three functional groups and were constructed using species of contrasting growth habits and phenologies. In addition, mixtures were designed to include warm-and cool-season species, which are similarly productive in the Atlantic Canadian climate (Wagg et al. 2021). This resulted in two-and three-species mixtures containing only legumes, grasses, brassicas, and forbs as well as one-, two-, and three-species mixture for all possible combinations of functional groups (Table 1). Besides single functional group mixtures, no mixture had more than one species of a functional group. One species from each functional group (brassica: brown mustard [Brassica juncea (L.) Czern.]; forb: buckwheat; grass: sorghum-sudangrass [Sorghum × drummondi (Nees ex. Steud.) Millsp. & Chase]; legume: crimson clover) was included in multiple mixtures to specifically evaluate these species inclusion in mixture. These species were chosen because they provide agroecosystem services that are of high priority to producers in the Atlantic Canadian region. Brown mustard and buckwheat have been shown to reduce wireworm [Coleoptera: Elateridae (Leach)] populations in subsequent potato crops (Vernon and van Herk 2017), sorghum-sudangrass has been shown to reduce Verticillium wilt, caused primarily by Verticillium dahliae (Kleb.), in subsequent potato crops (Larkin et al. 2011), and crimson clover is known for its soil-building properties and capacity to reduce erosion (Adetunji et al. 2020).
In all years, fields were chisel plowed and rolled before planting with a seed drill (Plotseed XL, Wintersteiger, Saskatoon, SK, Canada) with 15.24-cm spacing and a depth of 2.5 cm. Seeding rates of monocultures (Table 1) were determined from Extension recommendations, converted to a plant population density (plants m −2 ), and adjusted each year for germination percentage of seed. Seeding rates were halved in two-and three-way mixtures, except for buckwheat, which was divided by three according to Extension recommendations (Clark 2012; Table 1). All legumes were treated with Rhizobium leguminosarum (Frank 1889) before seeding. Plots measured 1.82 m by 10 m and were arranged in a randomized complete block design with 3 blocks (n = 468). Cover crops were planted on July 5, 2018, June 25, 2019, June 14, 2020, and June 14, 2021. Due to poor stand establishment in 2021, cover crops were replanted on July 13, 2021. Fertility, where required, was applied before planting and incorporated according to local recommendations. As is common in the region, cover crops were mowed with a flail mower in late September or early October depending on the year.

Data Collection
Each year and within each plot, we quantified cover crop and weed diversity and biomass production when >75% of cover crop species were flowering and had achieved canopy closure in mid-August. This stage was chosen to reflect outcomes of early-season competition and to measure biomass productivity before reproductive growth and plant senescence. A 0.5 m by 0.5 m quadrat (0.25 m 2 ) was randomly placed within each plot, and density of emerged cover crops was determined. All cover crops and weeds were then clipped at the soil surface, sorted by species, and dried in a crop dryer to constant weight at 80 C.

Statistical Analysis
All statistical analyses were conducted in the R environment v. 4.0.3 (R Core Team 2020). Diversity was measured using a variety of metrics, including species richness (S), which is the number of species in a given area; Shannon diversity (H), which extends species richness to account for relative abundances of species; and evenness (E), which defines how numerically equally a community is distributed among its constituent species. Species richness (S), Shannon diversity (H), and evenness (E) were calculated for cover crops with the R package VEGAN (Oksanen et al. 2020) for each plot using observed species densities. Shannon diversity was calculated as: where p i is the proportion of the ith species. Evenness was calculated as: where H is the Shannon diversity index, and S is species richness. Weed suppression for each treatment was calculated using weed biomass in the control of the corresponding block as: Annual ryegrass þ pearl millet þ sorghum-sudangrass LOLMU þ PESGL þ SORSU Grass mix (1) 3 10 12.5 20 Galaga þ alfalfa þ crimson clover GAGOF þ MEDSA þ TRFIN Legume mix (1) 3 5 10 7 Brown mustard þ oilseed radish þ tillage radish BRSJU þ RAPSO þ RAPSR Brassica mix (1) 3 3 5 4 Brown mustard þ sorghum-sudangrass þ crimson clover BRSJU þ SORSU þ TRFIN Three-way mix (3) 3 3 20 7 Buckwheat þ sorghum-sudangrass þ crimson clover FAGES þ SORSU þ TRFIN Three-way mix (3) 3 3 15 7 a Note, seeding rates were adjusted for germination percentage, those listed represent unadjusted seeding rates. Species 1, 2, and 3 refer to order listed in treatment. Due to low weed biomass in control plots in 2018, data for that year were removed from analysis of weed suppression. Data points lying 1.5 × outside the interquartile range were treated as outliers and removed before analysis.
To evaluate our first hypothesis, we took three different approaches. First, to determine cover crop and weed biomass productivity by all cover crop treatments, a linear mixed-effects model was constructed with cover crop treatment (n = 38) as the fixed effect and replication and year treated as random effects (Baraibar et al. 2018). The effect of year was evaluated with a Wald's test, and data were combined. Second, to evaluate effects of increasing species richness (1 to 3) on cover crop and weed biomass productivity, a linear mixed-effects model was constructed with species richness as the fixed effect and replication and year treated as random effects (Baraibar et al. 2018). Cover crop treatment was included as a random effect for the analysis of species richness to isolate effects of richness from community composition (Schmid et al. 2018). Further, to account for the unbalanced design, type III sums of squares were used. Linear mixed-effects models were constructed with the R package NLME (Pinheiro et al. 2021). To compare cover crop sown richness, least-squares means were generated with the R package EMMEANS (Lenth et al. 2021), and means were separated with a Tukey's honest significant difference (HSD) test. Effects were considered significant with a type I error rate of α = 0.05. Finally, we estimated the effect size of increasing functional richness on cover crop biomass productivity and weed biomass suppression according to Florence et al. (2019). In this analysis, the following linear model is constructed, which relates the standard deviation of cover crop biomass production to its mean with and without the interaction of functional richness and cover crop biomass as: A significant effect of the interaction between functional richness and cover crop biomass indicates functional richness alters this relationship. The effect of increasing functional richness on weed suppression was estimated by fitting the following exponential model with and without functional richness: The control was removed before analysis in all cases. To evaluate our second hypothesis, correlation and regression were used to investigate relationships between cover crop and weed biomass, weed suppression, species richness, Shannon diversity, and evenness. Linear and quadratic regression was conducted on the entire data set as well as split by functional group (Baraibar et al. 2018), and models with the lowest Akaike information criterion values were selected. Pearson's correlation coefficients were generated in all cases. Monoculture cover crops were removed from diversity regressions. To assess variability in cover crop evenness, the coefficient of variation (CV) was generated for each treatment and compared.
To assess stability of mixtures versus monocultures and the third hypothesis, the CV was calculated for cover crop biomass (n = 12) and weed suppression (n = 9) for each treatment calculated from all replicates and all years. The CV represents the spatial (location) and temporal (year) variability of biomass production and weed suppression of each treatment (Smith et al. 2020). As CVs cannot be compared statistically, treatments were ranked from lowest CV (most stable) to highest CV (least stable).

Cover Crop Productivity Varied by Species, Functional Group, and Diversity
Cover crop biomass production varied by treatment and sown richness (Figure 1; Table 2). Within monoculture cover crops, brassicas and legumes tended to produce less biomass than grasses and forbs; however, there was large variability within each group ( Figure 1A). On average, legumes produced the least biomass of all cover crops and were highly variable, ranging from a low of 6 g m −2 in galega (Galega officinalis L.) to a high of 252 g m −2 in field pea (Pisum sativum L.) ( Figure 1A). The grasses sorghum-sudangrass, pearl millet [Pennisetum glaucum (L.) R. Br.], and oat and the forb buckwheat were the most productive of all monocultures, producing between 388 and 476 g m −2 of biomass.
In general, cover crop mixtures were more productive than monocultures but not more than the most productive monocultures (Figure 1). Increasing sown richness, however, did increase cover crop biomass (Figure 2A). Within functional groups, grass mixes, forb mixes, and two-and three-way mixtures produced more biomass than brassica or legume monocultures and legume mixtures (Figure 1). Of mixtures, the two-way mixtures of buckwheat þ sorghum-sudangrass and sorghum-sudangrass þ red clover (Trifolium pratense L.) and the three-way mixtures of annual ryegrass [Lolium perenne L. ssp. multiflorum (Lam.) Husnot],þ pearl millet þ sorghum-sudangrass and brown mustard þ buckwheat þ sorghumsudangrass were the most productive and ranged from 426 to 468 g m −2 of cover crop biomass ( Figure 1B). Cover crop biomass increased with increasing density (Table 3; Figure 2B). Overall, while cover crop biomass production was positively correlated with cover crop density, density explained minimal variation across all cover crop treatments (Table 3). Within specific functional groups, however, a stronger relationship between density and biomass production was observed. Cover crop biomass of forbs and the two-and three-species brassica, legume, and forb mixes, as well as the three-way mixtures, was positively correlated with density (Table 3). Together, these results suggest that increased biomass production can be partly attributed to increased density; however, this relationship is species and functional group specific.

Evenness of Cover Crop Mixtures
Within mixtures, cover crop species evenness varied by mixture composition (Table 2). Significant variability in evenness was observed across cover crop mixtures ranging from high to low evenness ( Figure 2C). Grass and legume mixtures tended to be the least even of all functional groups, with mixtures containing sorghum-sudangrass consistently among those with the largest variability in evenness observed ( Figure 2C). In general, mixtures containing buckwheat and brown mustard had the least variability in evenness ( Figure 2C). For example, the addition of a highly Weed Science productive species such as buckwheat and/or brown mustard to a mixture containing sorghum-sudangrass reduced variability in evenness. In addition to differences across cover crop mixtures, as sown cover crop species richness increased from 2 to 3, evenness increased (Table 2; Figure 2D). Despite variability in evenness across cover crop mixtures, no relationship was observed between cover crop biomass production and evenness ( Figure 2E).  Table 1 for species (EPPO) codes.   Table 1 for species (EPPO) codes. The weed community across all years and treatments was dominated by smooth crabgrass [Digitaria ischaemum (Schreb) Schreb. ex Muhl.), common lambsquarters (Chenopodium album L.), redroot pigweed (Amaranthus retroflexus L.), quackgrass [Elymus repens (L.) Gould], and corn spurrey (Spergula arvensis L.), which comprised 92% of the total weed biomass. Within these, D. ischaemum and C. album were the most common species and constituted 77% of total weed biomass. Weed biomass was affected by cover crop treatment and sown cover crop richness (Figure 3). Within monocultures, weed biomass was greater in legume plots compared with other cover crop species, with a few exceptions ( Figure 3A). Weed biomass in the non-legume monocultures of annual ryegrass, timothy (Phleum pratense L.), brown mustard, and sunflower (Helianthus annuus L.) ranged from 100 g m −2 to 147 g m −2 and was similar to that observed for the majority of legume cover crops. In comparison, weed biomass in the most suppressive cover crops of oat, buckwheat, phacelia (Phacelia tanacetifolia Benth.), pearl millet, oilseed radish (Raphanus sativus L. var. oleiformis Pers.), and sorghum-sudangrass ranged from 12 g m −2 to 44 g m −2 . In general, weed biomass was reduced in cover crop mixtures as compared with monocultures ( Figure 3). Significant variability, however, was observed across all other functional groups ( Figure 3B), and increasing functional richness did not improve cover crop biomass productivity or weed suppression (Tables 4  and 5). The most weed-suppressive mixtures were both forb mixes (buckwheat þ phacelia [7 g m −2 ] and buckwheat þ phacelia þ sunflower [15 g m −2 ]), whereas legume mixtures were the least weed-suppressive functional group and comparable to legume monocultures. For example, the least suppressive legume mixtures were galega (244 g m −2 ) and galega þ alfalfa (Medicago sativa L.) þ crimson clover (188 g m −2 ), respectively. No relationship was found between cover crop Shannon diversity or evenness and weed biomass (data not shown).
Weed biomass declined with increasing cover crop biomass ( Figure 4A), density ( Figure 4B), and cover crop sown richness ( Figure 4C). The relationship between cover crop biomass and weed biomass, however, varied by functional groups (Table 3). Biomass of all functional groups except brassica and legume mixtures was correlated to weed biomass. Grass and forb mixtures had the highest correlation between cover crop biomass and weed biomass. Overall, cover crop biomass explained little variability in weed biomass. In comparison, density explained less variability in weed biomass than cover crop biomass (Table 3). Only forb monocultures and the three-way mixtures displayed a significant negative correlation between density and weed biomass (Table 3). Similar to the effects of cover crop biomass on weed biomass, weed suppression rapidly increased with increasing cover crop productivity ( Figure 4D). The regression of cover crop biomass and weed suppression demonstrated that cover crop biomass of 120 g m −2 , 335 g m −2 , and 450 g m −2 was required to achieve 50%, 80%, and 90% suppression of weed biomass, respectively ( Figure 4D). Maximum weed suppression (97%) was achieved when cover crop biomass was 655 g m −2 .

Importance of Cover Crop Biomass and Diversity for Weed Suppression
Across cover crop species, functional groups, and diversity, we observed a strong relationship between cover crop biomass production and weed biomass suppression; the most productive cover crops tended to be the most weed suppressive. As such, our results support the notion that biomass production is a key mechanism of weed suppression in annual cover-cropping systems and support our initial hypothesis. This is in agreement with numerous previously published reports demonstrating cover crop biomass as a key predictor of weed suppression for summer annual cover crops (MacLaren et al. 2019;Smith et al. 2015Smith et al. , 2020. Our observation that species such as buckwheat, sorghum-sudangrass, and oats are highly productive and provide a high degree of weed suppression is similar to results in other studies from the Northeast-Atlantic region (Smith et al. 2014(Smith et al. , 2015Wagg et al. 2021). Buckwheat has consistently been demonstrated to provide superior weed suppression through rapid growth, canopy development, and potential allelopathy (Falquet et al. 2015). Similarly, several studies have highlighted the productivity (Smith et al. 2014(Smith et al. , 2015 and allelopathic potential (Weston et al. 1989) of sorghum-sudangrass and its use as a weed-suppressive cover crop. Indeed, as observed in our study and elsewhere (Baraibar et al. 2018;Elsalahy et al. 2019), slow-growing legumes tend to produce little biomass and are poor annual cover crops for weed suppression. In addition, several mixtures, in particular ones containing buckwheat, brown mustard, and sorghum-sudangrass, were among the most weed-suppressive cover crops evaluated. This is in line with previous results demonstrating mixtures containing these species can provide weed suppression equivalent to the most superior monoculture (Baraibar et al. 2018;Nelson et al. 2012).
A degree of caution, however, is required with using biomass as the sole metric of plant competition (see Trinder et al. 2013Trinder et al. , 2020. Biomass alone provides an incomplete picture of plant-plant interactions and resource capture, as other factors such as relative growth rate and canopy light interception can be equally as important (Christina et al. 2021;Kruidhof et al. 2008;Lawley et al. 2012). Interactions between species functional traits, including rate of leaf appearance, resource use efficiency and capture, and emergence pattern in mixture may enhance or hinder weed suppression, further limiting the usefulness of biomass production as an indicator Table 4. Parameter estimates and standard error for linear models relating standard deviation (SD) of cover crop biomass (CCBio) production to mean cover crop biomass production with and without the effect of functional richness.

Model
Parameter Estimate ± SE a SD~CCBio CCBio 0.75 ± 0.02* SD~CCBio þ CCBio*FR CCBio 0.75 ± 0.04* CCBio*FR 0.002 ± 0.02 a Asterisked values (*) indicate functional richness altered the relationship between mean and SD of cover crop biomass production (α = 0.05). Table 5. Parameter estimates and standard errors for exponential models relating cover crop biomass production to weed biomass suppression with and without the effect of functional richness (FR).
of weed suppression in cover crop mixtures. Future studies evaluating these interactions between species are required to fully disentangle the weed-suppressive contribution of species to diverse cover crop mixtures.
In support of the majority of studies (Baraibar et al. 2018;Florence et al. 2019;Florence and McGuire 2020;Nelson et al. 2012;Smith et al. 2020), we observed cover crop mixtures were not necessarily more weed suppressive than the most productive monocultures. We did, however, detect a positive effect of increasing cover crop species richness on biomass production and weed suppression, which supported our second hypothesis. Further, the inclusion of highly productive species such as buckwheat or sorghum-sudangrass were integral to driving this relationship in mixture. These results are in line with the large body of evidence within the ecological literature wherein ecosystem functioning increases with increasing plant diversity (Cardinale et al. 2006;Isbell et al. 2011). We provide several possible explanations for the incongruences between our results and previous research demonstrating a poor relationship between diversity and weed suppression. First, our cover crop mixtures were constructed in a partial substitutive design; seeding rates were reduced by half in both two-and threeway mixtures. While greater cover crop density in three-way mixtures providing greater weed suppression may explain our response, we failed to detect a strong relationship between density and weed suppression in all other mixtures. We observed that several plots with the highest densities had similar levels of weed biomass as low-density cover crop plots. Similarly, we observed no differences in cover crop density between two-and three-way mixtures. Increasing monoculture seeding rate has been associated with superior weed suppression (Akemo et al. 2000;Weiner et al. 2001), but not necessarily cover crop biomass production (Ryan et al. 2011), due to greater intraspecific competition. In addition, several studies concerning cover crop mixtures using a substitutive versus additive design have similarly failed to detect strong evidence for increased weed suppression by increasing cover crop density in mixture (Mohler and Liebman 1987;Poggio 2005;Smith et al. 2015), despite larger effect sizes reported in the ecological literature with additive designs (Balvanera et al. 2006). Our analysis attempted to control for variable species communities and densities across treatments to isolate effects of richness on productivity and weed suppression. While we cannot fully disentangle the effects of increased cover crop density and diversity on biomass production and weed suppression, our data support the notion that increased cover crop diversity can promote agroecosystem services.
Results from our study demonstrate that increased biomass productivity and weed suppression are not necessarily associated with the dominance of a single highly productive grass species.  Baraibar et al. (2018) found that inclusion of cereal rye or oats in a mixture provided weed suppression equivalent to the most suppressive monocultures. Similarly, MacLaren et al. (2019) found the proportion of cereal biomass within a mixture to have a greater effect on weed suppression than cover crop biomass alone. Mixtures in past studies focused on few species and primarily consisted of grass-legume polycultures. We found that, regardless of functional identity, if a species was highly productive in monoculture, it was similarly productive in mixture. For example, evenness decreased in functionally rich mixtures containing sorghumsudangrass, while maintaining a high level of weed suppression similar to sorghum-sudangrass monocultures. This effect on evenness could be buffered against if the mixture contained a similarly productive species such as buckwheat, pearl millet, or annual ryegrass, resulting in reduced variability in mixture evenness. Other mixtures with variable evenness resulted from poor establishment by a cover crop species that also had poor establishment in monoculture, suggesting cover crop evenness may be unpredictable from year to year. For instance, mixtures containing faba bean (Vicia faba L.) tended to be the least even, and monocultures of faba bean similarly had poor productivity. Taken together, our results demonstrate that the inclusion of multiple highly productive species in Figure 5. Spatiotemporal stability of (A) cover crop productivity (n = 12) and (B) weed suppression (n = 9). Values are coefficients of variation and ranked from most stable to least stable across sites and years. See Table 1 for species (EPPO) codes. cover crop mixtures is essential to maintain stand evenness as diversity increases and ensure maximum productivity and weed suppression.

Spatiotemporal Stability of Cover Crops
To assess spatiotemporal stability of cover crop biomass production and weed suppression across replications and years, the CV was calculated for each treatment and ranked ( Figure 5). In terms of biomass production of the 10 most stable cover crops, 6 were mixtures, and all contained either buckwheat and/or brown mustard ( Figure 5A). The only monocultures within this group were the forbs buckwheat and phacelia, and grasses oat and annual ryegrass. In contrast, half of the least stable biomass-producing cover crops were legumes in monoculture, and only two were mixtures, brown mustard þ faba bean and the legume mixture of galega þ alfalfa þ crimson clover ( Figure 5A). Similar to biomass production, the greatest stability in weed suppression was observed in mixtures containing buckwheat and/or brown mustard and monocultures of oat, buckwheat, and pearl millet ( Figure 5B). In general, legumes and mixtures containing legumes provided the least stability in weed suppression and represented 7 of the 10 least stable cover crops.
The diversity-stability theory purports that species-rich ecosystems will be more spatiotemporally stable against environmental fluctuations than less species-rich ecosystems (Loreau et al. 2021). Several authors, however, have failed to find evidence in support of this hypothesis in annual cover-cropping systems (Florence et al. 2019;Smith et al. 2014Smith et al. , 2020. In contrast, we found that the majority of highly stable cover crops were mixtures and, therefore, partial support for our third hypothesis that diverse cover crop mixtures would provide greater stability in biomass production and weed suppression than monocultures. Though in congruence with Smith et al. (2020), monocultures of buckwheat and grass species such as pearl millet and oats displayed the greatest stability in weed suppression across sites and years. Therefore, while mixtures may not necessarily provide the most consistent weed suppression, many were among the most stable crops. Of mixtures, buckwheat and brown mustard were overrepresented among the most stable cover crops, suggesting they may have a greater influence on stabilizing productivity and weed suppression than other species under the variable Atlantic Canadian climate. The increase in cover crop stability through inclusion of these species in mixture may provide a means of bet-hedging against biotic and abiotic fluctuations producers face in the region and a mechanism to maximize weed suppression while providing additional agroecosystem services.
Our results demonstrate a positive effect of cover crop diversity on weed suppression; however, in line with previous studies, not all mixtures were more weed suppressive than the most productive monocultures. Monocultures of buckwheat, oat, pearl millet, or sorghum-sudangrass were often more productive and weed suppressive than the average mixture. This consistent result across regions, seasons, mixture composition, and functional diversity demonstrates that mixtures are not necessarily more productive or weed suppressive than monocultures. Proponents of cover crop mixtures suggest they provide an opportunity for increased diversity within annual cropping systems and provide increased agroecosystem services as compared with monocultures (Finney and Kaye 2017). If weed suppression is the primary goal, however, consistent evidence suggests that cover crop mixtures may not necessarily be the most effective option. Our results highlight the importance of species selection in cover crop mixture composition. In addition, complementarity can limit the dominance of highly productive species and maximize cover crop biomass and weed suppression. Within Atlantic Canada, cover crops are typically selected for their soil-building properties and to offset declines in soil organic matter (Mann et al. 2021;Nyiraneza et al. 2017). An important consideration of future research will be to assess weed suppression in addition to other agroecosystem services across a range of species and mixture compositions to fully evaluate the contribution of cover crop diversity to Atlantic Canadian cropping systems.