Introduction
Grasses such as giant foxtail (Setaria faberi Herrm.), yellow foxtail [Setaria pumila (Poir.) Roem. & Schult.], and large crabgrass [Digitaria sanguinalis (L.) Scop.], each resistant to several herbicide sites of action (Heap Reference Heap2019), are among the most common and problematic grass weeds in soybean [Glycine max (L.) Merr.] crop production systems in the United States (Van Wychen Reference Van Wychen2015, Reference Van Wychen2016). Barnyardgrass [Echinochloa crus-galli (L.) P. Beauv.] and jungle rice [Echinochloa colona (L.) Link] are two other troublesome monocot weeds in the midsouthern United States that have evolved resistance to seven and three herbicide mechanisms of action, respectively (Heap Reference Heap2019; Rouse et al. Reference Rouse, Roma-Burgos, Norsworthy, Tseng, Starkey and Scott2018; Schwartz-Lazaro et al. Reference Schwartz-Lazaro, Green and Norsworthy2017). Because herbicide options to control these weeds are limited, new management practices are urgently needed as weeds throughout the United States continue to develop herbicide resistance (Heap Reference Heap2019; Norsworthy et al. Reference Norsworthy, Walsh, Bagavathiannan, Bradley, Steckel, Kruger, Loux, Eubank, Davis and Johnson2014; Walsh et al. Reference Walsh, Broster, Schwartz-Lazaro, Norsworthy, Davis, Tidemann, Beckie, Lyon, Soni, Neve and Bagavathiannan2018).
Harvest weed seed control (HWSC), a nonchemical weed control approach that targets the collection and destruction of weed seeds during grain harvest, has helped Australian growers manage herbicide-resistant weed populations (Walsh et al. Reference Walsh, Newman and Powles2013). Potential effectiveness of HWSC systems depends upon seed retention of the target weed species at crop maturity, enabling its collection and processing at crop harvest and the effectiveness of the specific HWSC tactics employed (Walsh et al. Reference Walsh, Broster, Schwartz-Lazaro, Norsworthy, Davis, Tidemann, Beckie, Lyon, Soni, Neve and Bagavathiannan2018). Plants of many annual weed species shatter seeds at crop maturity in the United States (Davis Reference Davis2008; Norsworthy et al. Reference Norsworthy, Walsh, Bagavathiannan, Bradley, Steckel, Kruger, Loux, Eubank, Davis and Johnson2014; Schwartz-Lazaro et al. Reference Schwartz-Lazaro, Green and Norsworthy2017; Walsh et al. Reference Walsh, Broster, Schwartz-Lazaro, Norsworthy, Davis, Tidemann, Beckie, Lyon, Soni, Neve and Bagavathiannan2018). The efficacy of seed destruction necessary to reduce the soil seedbank using HWSC varies from 40% to 80% (Liebman and Davis Reference Liebman and Davis2009; Tidemann et al. Reference Tidemann, Hall, Harker and Alexander2016). Davis (Reference Davis2008) reported that S. faberi shattered 35% of seed in corn (Zea mays L.) and 45% of seed in soybean fields by harvest in east-central Illinois. Grass weeds such as jointed goatgrass (Aegilops cylindrica Host) and downy brome (Bromus tectorum L.) have been found to shatter a low proportion (<25%) of seed at crop maturity in eastern Colorado (Walsh et al. Reference Walsh, Broster, Schwartz-Lazaro, Norsworthy, Davis, Tidemann, Beckie, Lyon, Soni, Neve and Bagavathiannan2018). Preliminary field surveys of winter wheat (Triticum aestivum L.) fields near Pullman, WA, in 2013 found that 42% of Italian ryegrass [Lolium perenne L. ssp. multiflorum (Lam.) Husnot] seeds were shattered 15 cm (header height) above the soil surface at harvest (Walsh et al. Reference Walsh, Broster, Schwartz-Lazaro, Norsworthy, Davis, Tidemann, Beckie, Lyon, Soni, Neve and Bagavathiannan2018). However, lower seed retention (41% at harvest, and 32% at 1 mo later) of E. crus-galli in soybean was reported from Arkansas (Schwartz-Lazaro et al. Reference Schwartz-Lazaro, Green and Norsworthy2017).
In studies conducted in Alberta and Saskatchewan, Canada, green foxtail [Setaria viridis (L.) P. Beauv.], a common weed in the northern Great Plains, had high seed retention rates (≥80%) making it a suitable candidate for HWSC (Beckie et al. Reference Beckie, Blackshaw, Harker and Tidemann2017; Burton et al. Reference Burton, Beckie, Willenborg, Shirtliffe, Schoenau and Johnson2016). However, the lower seed retention (20%) for S. viridis observed in Minnesota cornfields at harvest (Forcella et al. Reference Forcella, Peterson and Barbour1996) limits the benefit of using HWSC for S. viridis in the region. The amount of weed seed retention at crop harvest varies among weed species and is influenced by agronomic factors and environmental conditions (Shirtliffe et al. Reference Shirtliffe, Entz and Van Acker2000; Taghizadeh et al. Reference Taghizadeh, Nicolas and Cousens2012). However, little research has been conducted to quantify the seed retention rates of various economically important monocot weeds in the United States, leaving the potential for HWSC systems to manage problematic grass weeds in U.S. cropping systems largely unknown. Here we present studies conducted to determine the seed retention of eight economically important grass weeds across the three major U.S. grain-producing regions.
Materials and Methods
Study Sites
We outlined a research protocol that included 11 states that were divided into three geographical regions: South-Central, Mid-Atlantic, and the North-Central regions. Field experiments were conducted in 2016 and 2017. Each state collected data both years, except for Tennessee, which only participated in 2016. Each location planted soybean using local standard practices described in local extension bulletins, including variety, seeding rate, row spacing, fertility, and other practices, and collected information on planting date, physiological maturity progression, and harvest date (Table 1).
a Regions include South-Central (SC): Arkansas (AR), Mississippi (MS), Tennessee (TN), and Texas (TX); North-Central (NC): Illinois (IL), Michigan (MI), and Missouri (MO); and Mid-Atlantic (MA): Delaware (DE), Maryland (MD), North Carolina (NC), and Virginia (VA).
b NA, data unavailable.
Data Collection
Sampling protocols were the same as the broadleaf species data collection in Schwartz-Lazaro et al. (Reference Schwartz-Lazaro, Shergill, Evans, Bagavathiannan, Beam, Bish, Bond, Bradley, Curran, Davis, Everman, Flessner, Harring, Jordan, Korres, Lindquist, Norsworthy, Sanders, Steckel, VanGessel, Young and Mirsky2021). Locally (within-state) problematic weeds were chosen for study for each state. Weeds that did not emerge from the soil seedbank were either seeded or transplanted into the crop. Transplanted weeds were of the same growth stage as those in the study field to mimic having germinated with the soybean crop. Weeds were transplanted in-row if the soil seedbank was not high enough to support a specific weed. A total of eight grass species were examined. Other than the individual weeds used in the studies, the soybean crop was kept weed-free throughout the growing season. Once the weeds began to flower, four seed-collection trays (F1721 Tray, T.O. Plastics, Clearwater, MN) were placed around the bottom of at least 10 randomly chosen plants to collect any seed shed from the plant. Trays were placed so that there was not a gap between the trays or the tray and the base of the plant. To help ensure trays captured shattered seed, if a plant spread over the outer edges of the trays during the course of the study, it was trained using twine and stakes to keep the entire plant over the trays. The greenhouse trays were emptied weekly using a portable vacuum and placed into envelopes for counting (see Schwartz-Lazaro et al. Reference Schwartz-Lazaro, Shergill, Evans, Bagavathiannan, Beam, Bish, Bond, Bradley, Curran, Davis, Everman, Flessner, Harring, Jordan, Korres, Lindquist, Norsworthy, Sanders, Steckel, VanGessel, Young and Mirsky2021). At the conclusion of the experiment, the plants were harvested to obtain a final seed count and determine the percentage of seed retention.
Data Processing and Statistical Analysis
Our analysis of grass species was conducted using the same methods as our broadleaf analysis. For details, readers should refer to the statistical methods in “Part 1: Broadleaf Species” (Schwartz-Lazaro et al. Reference Schwartz-Lazaro, Shergill, Evans, Bagavathiannan, Beam, Bish, Bond, Bradley, Curran, Davis, Everman, Flessner, Harring, Jordan, Korres, Lindquist, Norsworthy, Sanders, Steckel, VanGessel, Young and Mirsky2021). The emphasis of the analysis was to quantitatively and qualitatively describe the phenology of seed shatter at the site, species, and individual plant level in relation to soybean maturity. We will very briefly summarize the analysis here. All analyses were based on calculations of percent cumulative seed shatter over time, either at the site, species, or individual plant level. Seed shatter was calculated as the number of seeds that had shattered at a particular time point divided by the total seasonal seed production, including unshattered seed that was retained.
We plotted spatial heat maps to visualize regional to continental patterns in the rates of combined grass weed seed shatter during the weeks following soybean maturity. These were created using calculations of total cumulative seed shatter of all grass species studied within each state during the week of soybean maturity and at 2, 3, and 4 wk following maturity. States were only plotted on the map if they sampled during a given time interval. For example, if a state sampled within ±3 d of maturity (a 7-d window centered on the maturity date), we plotted it on the “week of maturity” map. To visualize how the distribution of seed shatter progressed at the species level, looking across states, the cumulative seed-shatter percentage values were converted to categorical groups and binned by increments of 10% (i.e., 0%<= shatter <10%, 10%<= shatter <20%, etc.), and the number of site-years in each bin was tallied for each species. These were then plotted as heat maps showing the percent of site-years for each species in each bin for each time interval. Finally, we calculated mean per capita daily seed rain rates (i.e., seeds plant−1 d−1) and mean per capita cumulative percent seed shatter for each species during the first 1 to 4 wk following maturity, accounting for site and year differences. The models used individual sample plants as the unit of replication with site, year, and their interaction as fixed effects. For analyses of seed rain rate, we used linear models with normally distributed errors. For analyses of percent seed shatter, we used generalized linear models with binomial errors for the proper fitting of proportion data. Because not all species were sampled in the same sites during both 2016 and 2017, the model structure had to be tailored to the data available for each species. The model structure and selection process are detailed in Schwartz-Lazaro et al. (Reference Schwartz-Lazaro, Shergill, Evans, Bagavathiannan, Beam, Bish, Bond, Bradley, Curran, Davis, Everman, Flessner, Harring, Jordan, Korres, Lindquist, Norsworthy, Sanders, Steckel, VanGessel, Young and Mirsky2021). We ran these tests with different model structures depending on data availability, because some species were not sampled in multiple sites, and others were only sampled in a single year. Although some species were studied in multiple sites and for multiple years, most were not studied in the same set of sites for both years. Thus, we were only able to fit site by year interactions in S. faberi (Schwartz-Lazaro et al. Reference Schwartz-Lazaro, Shergill, Evans, Bagavathiannan, Beam, Bish, Bond, Bradley, Curran, Davis, Everman, Flessner, Harring, Jordan, Korres, Lindquist, Norsworthy, Sanders, Steckel, VanGessel, Young and Mirsky2021). For one species, goosegrass [Eleusine indica (L.) Gaertn.], sampling was ended before the second week postmaturity and was only assessed at week 1. For two others, broadleaf signalgrass [Urochloa platyphylla (Munro ex C. Wright) R.D. Webster], and Texas millet [Urochloa texana (Buckley) R. Webster], sampling began after soybean maturity in 2016. All others were sampled at soybean maturity and at both 1 and 2 wk after maturity. All data processing and analyses were conducted in R (R Core Team 2018).
Results and Discussion
As in our study of broadleaf weed phenology (Schwartz-Lazaro et al. Reference Schwartz-Lazaro, Shergill, Evans, Bagavathiannan, Beam, Bish, Bond, Bradley, Curran, Davis, Everman, Flessner, Harring, Jordan, Korres, Lindquist, Norsworthy, Sanders, Steckel, VanGessel, Young and Mirsky2021), cumulative percent seed shatter was lowest in the southern U.S. regions and increased moving north through the states (Figure 1). This trend remained from soybean physiological maturity through maturity plus 4 wk. This result is consistent with previous studies that showed low seed retention (<40%) for S. faberi, E. crus-galli, A. fatua, spiny annual sow thistle [Sonchus asper (L.) Hill], and S. viridis in North American regions (Beckie et al. Reference Beckie, Blackshaw, Harker and Tidemann2017; Burton et al. Reference Burton, Beckie, Willenborg, Shirtliffe, Schoenau and Johnson2016; Forcella et al. Reference Forcella, Peterson and Barbour1996; Schwartz-Lazaro et al. Reference Schwartz-Lazaro, Green and Norsworthy2017; Shirtliffe et al. Reference Shirtliffe, Entz and Van Acker2000; Tidemann et al. Reference Tidemann, Hall, Harker, Beckie, Johnson and Stevenson2017). Further, the eight grass weed species’ percent shatter continually increased from soybean physiological maturity to harvest maturity (Tables 2 and 3; Figure 2). While the annual grass weeds have low seed retention, Johnsongrass [Sorghum halepense (L.) Pers.], a perennial grass weed, had high seed retention of >96% in Texas (Tables 2 and 3), which is similar to results reported by Walsh et al. Reference Walsh, Broster, Schwartz-Lazaro, Norsworthy, Davis, Tidemann, Beckie, Lyon, Soni, Neve and Bagavathiannan2018. This finding potentially indicates that the life cycle of the weed influences seed retention of a species.
a Seed rain rates were calculated from soybean physiological maturity to 1, 2, 3, or 4 wk after maturity. Cumulative seed shatter was calculated from the beginning of seed shatter through 1, 2, 3, or 4 wk after soybean physiological maturity. Cumulative shatter values are predicted from fitted logistic regressions for each species after accounting for differences due to states and years. χ2 values are from likelihood ratio tests comparing the fitted model with a null model. No test was performed for species with just a single site-year of data (indicated as “NA”), because we had already fit intercept-only null models to these.
b N is equivalent to the total number of plants for all sites and years.
c Model structures were dependent on the number of sites and years for each species. The model test used in seed rain rate analyses is determined by the model structure that was fit to each species: F-tests were used for seed rain rate models with state and/or year fixed effects; t-tests were used for intercept-only seed rain rate models; χ2 tests were used for likelihood ratio tests of binomial generalized linear models of seed shatter (%). No likelihood ratio tests were conducted for species with only 1 site-year of data.
a Regions include South-Central (SC): Arkansas (AR), Mississippi (MS), Tennessee (TN), and Texas (TX); North-Central (NC): Illinois (IL), Michigan (MI), and Missouri (MO); and the Mid-Atlantic (MA): Delaware (DE), Maryland (MD), North Carolina (NC), and Virginia (VA). A dashed line (—) indicates that there are no data for that time period.
The two most frequently examined grass species in the present study were E. crus-galli and S. faberi. Echinochloa crus-galli was examined by most of the South-Central region states and S. faberi by the North-Central and Mid-Atlantic regions (Figure 3). As for the broadleaf weeds (Schwartz-Lazaro et al. Reference Schwartz-Lazaro, Shergill, Evans, Bagavathiannan, Beam, Bish, Bond, Bradley, Curran, Davis, Everman, Flessner, Harring, Jordan, Korres, Lindquist, Norsworthy, Sanders, Steckel, VanGessel, Young and Mirsky2021), one of the more striking outcomes was the difference in variation across sites from year to year. In both species, there was little variation in seed-shatter progression between sites in 2016, while seed shatter was highly varied across sites in 2017. Both species were studied in more sites in 2017, but the pattern is noteworthy. Of the differences that were seen, more seed shatter occurred in 2017 overall than in 2016. The large range of percent seed shattered in these species could be due in part to annual differences in weather or regional differences. At soybean maturity the percent of seeds shattered ranged from 1% to 70% across species (Table 2). However, at 3 to 4 wk after soybean maturity, that range shifted to 5% to 100% (mean: 42%) seeds shattered. After accounting for site and year differences, E. crus-galli still retained over 80% of its seeds at 4 wk after soybean physiological maturity (Table 2) and considerably more during some sites and years (Table 3; Figure 3).
These results indicate that many summer annual grass species are likely not to be controlled consistently or adequately with HWSC, but S. halepense, a perennial, could be. While seed spread can be contained through HWSC, it will not manage rhizomes and other belowground perennial structures, so these must be managed by other means. Seed shatter in the annual grasses began before soybean maturity; thus, some additions to the soil seedbank had already been made by harvest. Soybean harvest can vary dramatically across regions, being earlier in the year in the southern United States and later (1 to 3 mo) in the northern United States. However, in some annual species (e.g., E. crus-galli), it may be possible to capture a significant amount (60% to 90%) of seed production with HWSC within 2 to 4 wk of soybean maturity during certain years, Variation within species across sites might indicate that the outcome will be variable between years and locations. While this would not result in eradication, it could lead to meaningful reductions in weed populations. Furthermore, header height, seed that is below the header, seed that is shattered at the header and not brought into the combine, and delayed harvest are all factors that could result in a limited amount of weed seed entering the combine to go through a HWSC tactic. More research is needed on what can be done to reduce inputs from grassy weeds to the soil seedbank as well as the amount of time that these weeds could begin to select for earlier shattering potential with the selection pressures of HWSC. Additionally, other economically important summer and cool-season annual grass weeds need to be evaluated for seed retention at harvest.
Acknowledgments
We would like to thank the staff and students at each university for helping conduct this research, specifically Kreshnik Bejleri, Sheri Heard, John Sanders, Barbara Scott, Annie Klodd, Zach Schaefer, Russ Garetson, Vitor Damiao, Matheus Martins, Camille Werner, Bruno Flaibam, and Camila Grassmann, and each institute’s research and experimental stations. The authors would also like to thank the U.S. Department of Agriculture, Agricultural Research Service Areawide program for funding and the HATCH Program of the National Institute of Food and Agriculture and the U.S. Department of Agriculture for providing partial funding for this work. No conflicts of interest have been declared.