Introduction
Cover cropping is a management strategy which adds potentially beneficial biodiversity to agroecosystems. Depending on specific cover crop management practices (Bergtold et al., Reference Bergtold, Ramsey, Maddy and Williams2019), farmers may leverage planned biodiversity to enhance regulating ecosystem services (Beillouin et al., Reference Beillouin, Ben-Ari, Malézieux, Seufert and Makowski2021; McClelland, Paustian, and Schipanski, Reference McClelland, Paustian and Schipanski2021; Tamburini et al., Reference Tamburini, Bommarco, Wanger, Kremen, van der Heijden, Liebman and Hallin2020), increase cropping system resilience (Reiss and Drinkwater, Reference Reiss and Drinkwater2018; Renwick et al., Reference Renwick, Deen, Silva, Gilbert, Maxwell, Bowles and Gaudin2021), reduce agricultural externalities, and support sustainable intensification (Wittwer et al., Reference Wittwer, Dorn, Jossi and van der Heijden2017). In particular, cover crops can increase beneficial insect populations (English-Loeb et al., Reference English-Loeb, Rhainds, Martinson and Ugine2003), reduce pest insect populations (Bugg, Reference Bugg1992; Bugg and Waddington, Reference Bugg and Waddington1994), support several aspects of soil health (Romdhane et al., Reference Romdhane, Spor, Busset, Falchetto, Martin, Bizouard, Bru, Breuil, Philippot and Cordeau2019; Unger and Vigil, Reference Unger and Vigil1998), reduce soil erosion (Novara et al., Reference Novara, Gristina, Saladino, Santoro and Cerdà2011), reduce pollutants in agricultural runoff (Dabney, Delgado, and Reeves, Reference Dabney, Delgado and Reeves2001), increase crop yield stability (Gaudin et al., Reference Gaudin, Tolhurst, Ker, Janovicek, Tortora, Martin and Deen2015), increase farm profitability (Correia et al., Reference Correia, Brito, Sampaio, Dias, Bacelar, Gonçalves, Ferreira, Moutinho-Pereira and Rodrigues2015), sequester atmospheric carbon (Novara et al., Reference Novara, Minacapilli, Santoro, Rodrigo-Comino, Carrubba, Sarno, Venezia and Gristina2019), and otherwise increase the ecological value of farmland.
Whereas a large body of research has rigorously addressed the ecological impacts of annual cover crops grown in the fallow period between two annual cash crops (e.g., Teasdale, Beste, and Potts Reference Teasdale, Beste and Potts1991), there is less information about cover crop impacts in perennial systems. In contrast to cover crops in annual cropping systems, orchard cover crops are grown under the tree canopy and rely on spatial separation to avoid competition with the main cash crop. Cover crops have demonstrated impacts on abiotic factors in perennial systems, including improving soil structure (Ramos et al., Reference Ramos, Benítez, García and Robles2010; Walsh et al., Reference Walsh, MacKenzie, Salmins and Buszard1996), increasing soil nutrition (Sánchez et al., Reference Sánchez, Giayetto, Cichón, Fernández, Aruani and Curetti2007), increasing water use (Monteiro and Lopes, Reference Monteiro and Lopes2007), and reducing orchard temperature (O'Connell and Snyder, Reference O'Connell and Snyder1999). More research is needed to fully understand how the biotic functioning of orchard cover crops affects horticultural management.
Such research would support ecological systems-based approaches to agricultural sustainability, such as integrated pest (or weed) management, which rely on biodiversity and regulating ecosystem services to support multiple aspects of the agroecosystem (Haring, Reference Haring2021). Integrated pest management highlights the practical importance of basic ecological knowledge. For example, weeds can indicate the absence of unfilled ecological niches within the orchard system (Smith, Mortensen, and Ryan, Reference Smith, Mortensen and Ryan2010). Therefore, reducing available ecological niches (such as by planting predictable, domesticated plants that reduce resource availability) could displace weeds (Kruidhof, Bastiaans, and Kropff, Reference Kruidhof, Bastiaans and Kropff2008; Mirsky et al., Reference Mirsky, Curran, Mortenseny, Ryany and Shumway2011), suppress herbicide-resistant weeds (Bunchek et al., Reference Bunchek, Wallace, Curran, Mortensen, VanGessel and Scott2020; Pittman, Barney, and Flessner, Reference Pittman, Barney and Flessner2019), and provide additional sustainability benefits (Liebman and Davis, Reference Liebman and Davis2000). Whereas conventional orchards have significant unused resource pools that lead to the need for intensive vegetation control, cropping systems with diverse ground covers limit weed proliferation by regulating resource availability and safe sites for seed germination (Adeux et al., Reference Adeux, Vieren, Carlesi, Bàrberi, Munier-Jolain and Cordeau2019).
Winter annual cover crops have a life cycle that is coincidental with winter rains in Mediterranean climates as well as the dormant period of deciduous orchard trees (Baumgartner, Steenwerth, and Veilleux, Reference Baumgartner, Steenwerth and Veilleux2008; Bugg et al., Reference Bugg, McGourty, Sarrantonio, Lanini and Bartolucci1996; DeVincentis et al., Reference DeVincentis, Solis, Rice, Zaccaria, Snyder, Maskey, Gomes, Gaudin and Mitchell2022). This phenology allows winter annual weeds to have significant temporal niche differentiation compared to the orchard crop while allowing exploitative competition between cover crops and orchard weeds. Other forms of interference, such as allelopathy or suppression of summer weed germination with cover crop residues (Creamer et al., Reference Creamer, Bennett, Stinner, Cardina and Regnier1996; Putnam, Defrank, and Barnes, Reference Putnam, Defrank and Barnes1983) may also be regulating weed suppression potential of cover crop communities. The relative contribution of these mechanisms varies seasonally and over the life cycle of the orchard as changes in resource availability often alter the phenology of competition (Pearson, Ortega, and Maron, Reference Pearson, Ortega and Maron2017).
Life cycle is just one example of many functional traits that affect how the functional ecology of a cover crop species can support weed suppression, and other functions such as nutrient acquisition or germination ecology could be critical (Smith and Gross, Reference Smith and Gross2007). Many researchers have investigated the performance of cover crop mixes based on the functional traits of component species (Haramoto and Pearce, Reference Haramoto and Pearce2019; Kaye et al., Reference Kaye, Finney, White, Bradley, Schipanski, Alonso-Ayuso, Hunter, Burgess and Mejia2019; McKenzie-Gopsill et al., Reference McKenzie-Gopsill, Mills, MacDonald and Wyand2022; Ramírez-García et al., Reference Ramírez-García, Carrillo, Ruiz, Alonso-Ayuso and Quemada2015; Smith, Atwood, and Warren, Reference Smith, Atwood and Warren2014; Trichard et al., Reference Trichard, Alignier, Chauvel and Petit2013). By combining multiple species with complementary functional traits that match specific management goals, multispecies cover crops can enhance the resilience of an agroecosystem while providing several ecosystem services. Studies in unmanaged ecosystems likewise highlight the importance of diverse, multifunctional plant communities for reducing invasibility through mechanisms like niche differentiation (e.g., through resource partitioning and phenological differences) and variation in developmental biology and phenotypic plasticity (Levine and D'Antonio, Reference Levine and D'Antonio1999; Naeem et al., Reference Naeem, Knops, Tilman, Howe, Kennedy and Gale2000; Tilman, Wedin, and Knops, Reference Tilman, Wedin and Knops1996).
Despite a growing body of research on the functional ecology of multispecies cover crops, the relationships between cover crop species in a mix are frequently distinctive and unpredictable (Stefan, Engbersen, and Schöb, Reference Stefan, Engbersen and Schöb2021). There remains uncertainty in how these patterns emerge at scales relevant to intensive, commercial agricultural given increased disturbance and decreased species richness in such systems (Bybee-Finley, Mirsky, and Ryan, Reference Bybee-Finley, Mirsky and Ryan2017; Smith, Warren, and Cordeau, Reference Smith, Warren and Cordeau2020). More practical and hands-on information about the functional diversity of multispecies cover crop mixes at large scales could improve how orchard managers employ complementarity to design species mixes that balance multiple management goals. This study aims to fill a critical knowledge gap in understanding the ecological effects of multispecies cover crops when they are integrated into existing large-scale, intensified orchard cropping systems.
Our goal for this study was to determine how cover crop mixes with different levels of diversity perform within existing orchard management systems, and our experiment and analyses were selected to focus on plant communities across a range of variable, real-world conditions found in commercial almond orchards. We hypothesize that both functionally uniform and diverse cover crop mixes can provide living ground cover that displaces weeds, but that functionally diverse cover crops will emerge and compete for resources more consistently across growing seasons and locations and are therefore better able to impact weed community composition. To evaluate these hypotheses, we examined four indicators of cover crop function: (1) area of bare ground underneath cover crops and weedy resident vegetation, (2) relative ground cover of cover crops and weeds, (3) stability of ground cover provided by cover crops over space and time, and (4) impacts on weed communities.
Materials and methods
We evaluated plant communities under two multi-species cover crop mixes, one that is functionally diverse and one that is functionally uniform, over two seasons in three commercial almond (Prunis dulcis (Mill.) D.A. Webb) orchards in the Central Valley of California. We used two different cover crop mixtures that were designed to fulfill different agroecological goals within almond orchards. Namely, we used a mix with high functional diversity that is used commercially for improving soil structure and limiting soil erosion because its constituent species come from plant families with contrasting root architecture among other characteristics, as well as a mix with lower functional diversity that is used for providing floral resources to foraging pollinators with species from one plant family that have overlapping flowering phenology. We evaluated their impact on weed population density and species communities across a wide geographical area in central California.
Experimental design and management
We created replicated, large-plot experiments in commercial almond orchards in Tehama, Merced, and Kern Counties in California (Table 1). These locations span nearly 600 km in the Central Valley of California. This region produces over 99% of the almonds in the United States (Anonymous, 2020). The experiment used a randomized block design with four replicates of three or four different ground cover treatments at each site. Ground cover treatments were implemented for two years on the same plots, beginning in the fall of 2017 and ending in the late summer of 2019 (Table 2). These treatments represented commercial-standard management practices and two different five-species cover crop mixes with different levels of functional diversity. Plots were 25.5 m wide at each site, encompassing four orchard alleys and three tree rows with another tree row between plots. In the perpendicular direction, plots extended the entire length of the orchard (195 m at the Tehama site, 385 m at the Merced site, and 320 m at the Kern site).
The study was designed to use commercially relevant spatial and temporal scales, and orchard management was determined by grower cooperators for agronomic relevance. All orchards were equipped with microsprinkler irrigation, and irrigation schedules were determined based on almond evapotranspiration models in accordance with local weather conditions and recommendations (Fereres and Goldhamer, Reference Fereres and Goldhamer2003). Irrigation, insecticide, fungicide, and fertilizer treatments and rates were determined by each grower and applied to the tree rows only. Tree rows were maintained with conventional herbicide programs to create vegetation-free zones at the base of trees. These herbicide applications were performed with shielded sprayers at each site, and we did not observe any herbicide injury that would indicate herbicide drift to the cover crops throughout the experiment. Each of the sites was subjected to regular traffic from machinery and farmworkers to complete normal orchard management operations throughout the cover crop growing season.
Ground cover treatments
Two different winter cover crop mixes were planted in orchard alleyways and compared against two different control treatments that reflect mainstream orchard vegetation management programs (Table 3). These cover crop mixes were selected because of their current commercial availability in California and their contrasting functional diversity. Both mixes have the same number of species but represent different levels of functional diversity, with the uniform mix having five quite similar species and the diverse mix having species from several plant families and contrasting biological characteristics. Each mix was established at rates and with methods based on recommendations from the seed supplier.
The ‘uniform’ mix consisted of five functionally similar species, and this mix is used commercially in California and distributed as ‘PAm Mustard Mix’ by the Project Apis m. (Salt Lake City, UT, USA) Seeds For Bees program because its component species provide abundant floral resources for pollinators but slightly different flowering phenologies. The ‘diverse’ mix consisted of five species from the grass, brassica, and legume groups that are commonly used together in functionally diverse cover crop mixes (Altieri et al., Reference Altieri, Lana, Bittencourt, Kieling, Comin and Lovato2011) to support soil health by providing abundant biomass, multiple kinds of root architectures, and complementary resource needs.
The ‘resident’ vegetation treatment involved winter vegetation management primarily with mowing. The ‘bare’ treatment involved more intensive vegetation management, including one to two broadcast herbicide applications in the winter. For the bare treatment, herbicide applications timing and product choice were determined by grower cooperators but included broad spectrum herbicides without residual activity, which are common in California almond production, including glyphosate and glufosinate. The Tehama site included only the resident treatment due to the preference of our cooperator and to better reflect standard practices in this region of California which has more abundant winter rainfall. The Merced and Kern sites featured both the resident and bare treatments to better reflect high-intensity production systems in these regions.
Data collection
Orchard alley plant communities were evaluated with point-intercept transects. Transect surveys coincided with cover crop flowering for most species as well as winter weed flowering for many endemic species in the study area (dates are in Table 2). Each plot was surveyed with a single 50 m long transect with points observed evenly at each meter along the transect, and observations from all 50 points were combined to estimate plant cover across each plot. Transects were placed in the same location in each plot, with the starting point located 75 m from the end of the plot and inside the second tree row from the side of the plot. Transects extended diagonally across a single orchard alley, starting and ending on opposite edges of the planted swath.
At each point along the transect, plant cover was observed for the top layer of vegetation as observed from above, which varied depending on the height of target species. The observed vegetation type (plant species or bare ground) was recorded. Plants were identified to species visually in the field, except in the case of the white and yellow mustards in the uniform mix which were identified as one operational taxonomic unit due to morphological similarities at the observed growth stage, and the observed plant species or bare ground was recorded. The total number of occurrences of each vegetation type along each transect were summed and converted to a percentage to estimate ground cover in each plot. Transects described plant community composition for each plot, as well as bare ground (the portion of each transect that was not covered by any kind of vegetation) and ground cover from cover crops or weeds (found by counting the number of occurrences of all of the species associated with each vegetation type).
Statistical analysis
Analyses were performed in R 4.2.3 (R Core Team, 2023). To evaluate the first objective (evaluating ground cover and reduction of bare ground), comparisons of bare ground among treatments were made with ANOVA. ANOVA assumptions were inspected visually with qqPlot from the car package (Fox and Weisberg, Reference Fox and Weisberg2019). We detected a heavy tailed distribution and subsequently repeated the analyses using an arcsine square root transformed response variable as well as with a generalized linear model with a Poisson distribution. However, results were similar among the different analyses, and we display results from the untransformed ANOVA below. One outlier was identified with the Bonferroni outlier test using outlierTest. This outlier value was excluded from further analyses because it was collected in the same plot at the Merced site that had been previously excluded because it had not been planted in 2017 (i.e., no data from this plot from either study year were included). The model we used included these fixed effects as predictors: treatment, year, site, the interaction between year and site, and replicate nested within site. ANOVA was performed with Anova from the car package using type II sums of squares. Multiple comparisons were made with least-squares means using the emmeans package (Lenth, Reference Lenth2021).
For the second objective (evaluating tradeoffs in relative ground cover between cover crops and weeds), we evaluated the general relationship between cover crop and weed cover within cover crop treatments (i.e., not including the bare or resident treatments) across sites and years in this study. The relative ground cover from cover crops and weeds were modeled with the lm function in base R, and we used Anova for hypothesis testing. We used weed cover as the response variable and cover crop cover and cover crop treatment as predictors.
To evaluate the third objective (evaluating stability of cover crop cover over space and time), we assessed cover crop stability by comparing coefficients of variation for cover of each cover crop mix as pooled across sites and years in this study. Pooled coefficients of variation and their 95% confidence intervals were calculated with the ci.cv function in the MBESS package (Kelley, Reference Kelley2022) before they were compared with the modified signed-likelihood ratio test as implemented in the cvequality package (Marwick and Krishnamoorthy, Reference Marwick and Krishnamoorthy2019). These comparisons are based on recommendations by Reiss and Drinkwater (Reference Reiss and Drinkwater2018). For the final objective (evaluating impacts on weed communities), weed communities in the different cover crop treatments were analyzed with nonmetric multidimensional scaling (NMDS). We evaluated weed community groupings based on treatment, year, and site, as well as treatment and year within each site. NMDS was based on Bray–Curtis dissimilarity and was calculated using the metaMDS function in the vegan package (Oksanen et al., Reference Oksanen, Blanchet, Friendly, Kindt, Legendre, McGlinn, Minchin, O'Hara, Simpson, Solymos, Stevens, Szoecs and Wagner2020). We evaluated grouping variables using anosim, also from vegan, with 9999 permutations.
Results
Bare ground
Cover crop treatment (F 3,65 = 93.23, P < 0.001), site (F 2,65 = 30.21, P < 0.001), and their interaction (F 5,65 = 10.56, P < 0.001) had significant effects on the amount of bare ground observed in orchard alleys, while year (F 1,65 = 1.34, P = 0.251) and block (F 9,65 = 0.84, P = 0.582) did not (Fig. 1). Overall, the uniform and diverse mixes resulted in similar levels of bare ground (P = 0.289), at 15.9 ± 3.09 and 11.4 ± 2.89%, respectively, when averaged across sites and years. These values are less than the 39.3 ± 2.89% bare ground in the resident vegetation treatment (P < 0.001 for both comparisons).
Within the Kern site, each of the cover crop treatments resulted in significantly different (P = 0.003) levels of bare ground from one another, with the diverse mix (6.0 ± 4.64%) resulting in less bare ground than the uniform mix (28.0 ± 4.64%). Within the Merced site, the diverse mix (25.0 ± 6.55%) resulted in a similar level of bare ground to both the uniform mix (14.8 ± 7.88%, P = 0.332) and the resident vegetation treatment (36.5 ± 6.55%, P = 0.228). The Tehama site had low levels of bare ground across all three treatments, which were similar to one another (P > 0.25 for all comparisons). Both cover crop mixes had somewhat differing performance within sites, but both cover crop mixes on average reduced bare ground compared to standard commercial management practices.
Weed and cover crop cover
Across this study, cover crop cover was negatively associated with weed cover (Fig. 2; slope = −0.74, R 2 = 0.83, P < 0.001). When including cover crop treatment as a predictor, we found that cover crop cover was significant (F 1,43 = 176.72, P < 0.001) while cover crop treatment was not (F 1,43 = 0.35, P = 0.56). Regardless of the cover crop mix, we observed that increased ground cover from cover crops resulted in reduced ground cover from weeds. This relationship was described by a line with a slope less steep than negative one, indicating that every increase in the cover crop canopy covered some amount of bare ground in addition to the displaced weed vegetation.
Cover crop stability
The coefficient of variation for cover crop cover from the diverse mix was 48.6%, significantly less variation than the 91.5% variation observed in the uniform mix (Fig. 3; P = 0.035). Across the experiment, the diverse mix resulted in more consistent levels of ground cover than the uniform mix, reliably creating ground cover across a range of geographical, environmental, and management-related variation.
Weed communities
The Tehama site had three to five times greater species richness than either of the other sites (Table 4). All the sites were dominated by small-statured, winter annual weed species such as annual bluegrass (Poa annua L.) and common chickweed (Stellaria media (L.) Vill.). The remainder of the species at all sites included both grass and broadleaf (dicotyledonous) species, with some additional summer annual and perennial species present at the Tehama site. Cover crops influenced weed communities but to different extents depending on the site and year. When using ANOSIM to test whether weed communities had similar constituent species (stress = 0.166), we found that weed communities were strongly associated with their site (R = 0.569, P < 0.001) and cover crop treatments (R = 0.091, P = 0.005). However, effect sizes were generally small and clusters can be difficult to identify visually when weed communities from each plot were plotted in nonmetric scaled space (Fig. 4). As described above, the differences between weed community diversity at each site were relatively clear, while there were few noticeable differences between weed communities in different cover crop treatments.
Species are listed in order of prevalence, based on the cumulative total number of observations at each site.
Because of the inherent differences between weed communities at each site, we also analyzed similarity of weed communities within each site individually. Weed communities were similar across years and treatments at Merced, and weed communities remained relatively sparse and homogenous throughout the experiment at that site (stress = 0.119). At the Tehama site, weed communities differed across years (R = 0.981, P < 0.001; stress = 0.150), which is logical given that cover crop establishment was very strong in 2018 to the point that we observed no weeds in the cover crop treatments that year. Weed communities differed with cover crop treatment at the Kern site (R = 0.316, P < 0.001; stress = 0.036), though we observed few qualitative differences in weed communities at the Kern site.
Discussion
Orchard cover crop mixes, as implemented in this study, were effective at establishing, reducing bare ground, and suppressing weeds. However, these effects were variable, and there is little evidence that the cover crop mixes we used had fixed impacts on the composition of orchard weed communities. Differences among sites, which could include climate and management factors, contributed to some of this variability. Our goal was to observe cover crops as they would be implemented by commercial almond growers in California, which entailed a variety of unique management decisions at each site. Growers will always implement some level of site-specific management that could affect cover crop performance, but it remains encouraging that cover crops that are sold and used commercially for different purposes could result in weed suppression across the three locations of the experiment. The diverse mix resulted in more consistent ground cover in this study, and ground cover led to greater weed suppression. However, cover crop cover, not the specific cover crop treatment, was the primary driver for weed suppression and other effects on weeds in this study.
Cover crop abundance
This study demonstrates the critical importance of cover crop abundance for weed competition. As hypothesized, the diverse cover crop mix more reliably created an abundant cover crop compared to the uniform cover crop mix. However, the diverse cover crop mix was not inherently more competitive, and both cover crop mixes suppressed weeds when abundant. The present study is currently the largest-scale study focusing on weed suppression and cover crop mixes in orchard systems and demonstrates that weed-suppressing cover crops can be integrated into existing commercial production systems.
The importance of cover crop abundance is consistent with numerous existing studies of weed suppressing cover crops in annual cropping systems (Bybee-Finley, Mirsky, and Ryan, Reference Bybee-Finley, Mirsky and Ryan2017; Creamer et al., Reference Creamer, Bennett, Stinner, Cardina and Regnier1996; Florence et al., Reference Florence, Higley, Drijber, Francis and Lindquist2019; MacLaren et al., Reference MacLaren, Swanepoel, Bennett, Wright and Dehnen-Schmutz2019; Smith, Warren, and Cordeau, Reference Smith, Warren and Cordeau2020). It is important to note that these studies primarily observed plant abundance by measuring cover crop and weed biomass, while the present study came to a similar conclusion by measuring abundance through ground cover. While cover and biomass are distinct from one another, both are useful measures of plant abundance (Wright, Reference Wright1991). The present study confirms through field experiments that ground cover is an important measure of cover crop abundance, which can support higher-throughput, non-destructive cover crop research compared to previous studies that rely on biomass collection.
Previous research has established that an abundant and competitive cover crop is likely to have features that contribute to weed suppression as well as other ecosystem services. For example, large and abundant cover crops are effective in exploitative competition due to asymmetric resource acquisition, such as root competition for soil nutrients (Weiner, Reference Weiner1990). Large and abundant root systems can also improve soil structure and increase soil organic matter (Pierret et al., Reference Pierret, Maeght, Clément, Montoroi, Hartmann and Gonkhamdee2016) or create traps for pest insects while creating safe sites for beneficial insects (Hassanali et al., Reference Hassanali, Herren, Khan, Pickett and Woodcock2008). This abundance could also cause other challenges for orchard managers. For example, cereal rye was included in the multifunctional mix in this study and is known to be an important component species for weed suppression (Akemo, Regnier, and Bennett, Reference Akemo, Regnier and Bennett2000; Barnes and Putnam, Reference Barnes and Putnam1983). However, large amounts of cereal rye residue are frequently reported by nut growers as a concern because of potential interference with on-ground nut harvest. While this experiment focused on winter management, and thus we did not evaluate almond yield, future development of orchard cover crops should aim to identify cover crop termination strategies that create acceptable conditions for nut harvest in the summer.
Cover crops and plant diversity
This study highlights that cover crop mixes with different levels of functional diversity can successfully support weed management across different orchard contexts. The addition of cover crop mixes in this study contributed significantly to orchard plant diversity, essentially doubling species richness in the mature orchards (Merced and Kern sites). Weed species richness was highest in the young orchard (Tehama site), where the orchard floor was relatively unshaded and still populated with many weed species carried over from the previous pasture system. This study implemented cover crops on a time scale relevant for adoption in contemporary orchards, but more research is needed to understand the cumulative effects of cover crop competition on weed community assembly over the decades-long lifespans of commercial orchards. Maintaining biodiversity is a major challenge for agroecosystems. While there is ongoing conflict between promoting functional diversity and achieving some vegetation management goals, this study reinforces the idea that cover crops are flexible tools that can support multiple management goals if planned appropriately (Crézé and Horwath, Reference Crézé and Horwath2021; De Leijster et al., Reference De Leijster, Santos, Wassen, Ramos-Font, Robles, Díaz, Staal and Verweij2019; Mia et al., Reference Mia, Massetani, Murri and Neri2020).
Cover crop mixes in this study were primarily selected for their existing commercial uses, which are intended to address almond management goals other than weed suppression. In this study, the diverse mix provided functional diversity that led to more stability. We intentionally focused on the end performance of existing multifunctional cover crops, but optimization of cover crop mixes through agronomic management programs and additional species selection could improve weed suppression as well as other ecosystem services (Haring and Hanson, Reference Haring and Hanson2022). Multifunctional cover crop mixes could be improved by additional recognition of the specific ecological relationships between constituent species and their resulting agroecosystem services (Baraibar et al., Reference Baraibar, Hunter, Schipanski, Hamilton and Mortensen2018; Finney and Kaye, Reference Finney and Kaye2017; Ingels et al., Reference Ingels, Horn, Bugg and Miller1994; Schipanski et al., Reference Schipanski, Barbercheck, Douglas, Finney, Haider, Kaye, Kemanian, Mortensen, Ryan, Tooker and White2014). Pest management is generally used to reduce biodiversity in cropping systems, but cover crops provide an opportunity to support pest management goals while simultaneously promoting biodiversity through functional vegetation management. Agricultural systems are designed to support ample plant growth, and this resource-rich environment could be more practical and efficient if it supported competitive, complementary, and useful biodiversity instead of unwanted weedy plants.
Data availability statement
Data are available from https://doi.org/10.25338/B8S63T
Acknowledgments
The authors gratefully acknowledge the contributions of our collaborators throughout this field project: Cynthia Crézé, Mae Culumber, Kent Daane, Amanda K Hodson, Danielle M Lightle, Jeffrey Mitchell, Andreas Westphal, Houston Wilson, Mohammad Yaghmour, and Cameron At Zuber. We also acknowledge the grower-cooperators that supported this project: Steve Gruenwald, Castle Farms, and Wegis & Young. This research took place on lands that are the ancestral homes of the Nomlaki, Patwin, and Yokuts peoples, and these peoples remain committed to the stewardship of these lands today.
Author contributions
All authors developed hypotheses, contributed to experimental design, and critically revised the manuscript. S. H. led field data collection, data analysis, and manuscript drafting. Stakeholders, including orchard growers and cooperative extension professionals, were included in project design and management, and specific parties are listed in the acknowledgements.
Funding statement
The Almond Board of California funded this research under project numbers 16-STEWCROP7 and 18-HORT12.
Competing interests
We have no conflicts interest to declare.