Supporting Beneﬁcial Insects for Agricultural Sustainability: The Role of Livestock-Integrated Organic and Cover Cropping to Enhance Ground Beetle (Carabidae) Communities

: Ground beetles (Carabidae) are beneﬁcial insects providing ecosystem services by regulating insect pests and weed seeds. Despite several studies conducted on ground beetles worldwide, there is a lack of knowledge on how these insects are a ﬀ ected by di ﬀ erently managed organic systems (e.g., tillage-based versus grazed-based) compared to that of chemical-based no-tillage conventional cropping systems. In a 5-year (2013–2017) study, we assessed the ground beetle communities in cover crops and winter wheat ( Triticum aestivium L.) in Montana, USA, with three contrasting cropping systems: a chemically managed no-tillage, a tillage-based organic, and a livestock-integrated organic with reduced tillage. The ﬁrst three years (i.e., 2013–2015) corresponded to the transition to organic period, while the last two (i.e., 2016–2017) were conducted in United States Department of Agriculture (USDA) organic-certiﬁed tillage-based and livestock-integrated organic systems. The experiment was designed with three management systems across three blocks as the whole plot variable and 5-year rotation of crop phases as the subplot variable. Using pitfall traps, we sampled ground beetles across all cover crop and winter wheat subplots for ﬁve years ( n = 450). The data were analyzed using mixed e ﬀ ects models and PERMANOVA and visualized with non-metric multidimensional scaling ordination. Our study indicated that organically managed farms, whether tilled or grazed, enhance activity density, species richness, diversity, and evenness of ground beetles in the dryland row crop productions. Also, irrespective of farming system, cover crops supported higher species richness, diversity, and evenness of ground beetles than winter wheat. The ground beetle communities were mostly similar during the transition to organic period. However, during the established organic phase, cropping systems acted as contrasting ecological ﬁlters and beetle communities became dissimilar. Cover cropping a ﬀ ected ground beetle communities positively not only in organically managed systems but also in chemical-based conventional systems. Our study provides evidence supporting the adoption of ecologically-based cropping systems such as crop-livestock integration, organic farming, and cover cropping to enhance beneﬁcial insects and their pest-regulation services.


Introduction
Industrialized crop production relies heavily on the use of off-farm chemical and mechanical inputs to control pest populations, maintain soil fertility, and prepare fields for planting. While these among three contrasting agricultural management systems: chemical-based conventional no-till, tilled organic, and integrated-crop-livestock organic production with reduced tillage. Targeted grazing was used to terminate cover crops and manage weeds during pre-seeding and post-harvest. This study encompassed five years with two temporal periods in the organic systems: the three years of transition phase and the first two years of the United States Department of Agriculture (USDA) organic certified phase. We hypothesized that because conventional no-till, tilled organic, and integrated crop-livestock organic cropping systems use different farm inputs and management, they act as distinct ecological filters affecting the ground beetle communities differently. Also, based on previous studies [29][30][31], we expected that both organic systems and cover cropping, compared to conventional and cereal cropping, would enhance ground beetle abundance and diversity.  Table S1). The site is underlain with Blackmore silt loam (fine-silty, mixed, superactive, frigid Typic Agruistoll) [32].

Site Description and Cropping History
Prior to 2004, the study site was planted with perennial grasses (Bromus inermis L., Thinopyrum intermedium (Host) Barkworth and D.R. Dewey, and Poa compressa L.). Between 2004 and 2009, the experimental site followed either continuous spring wheat (Triticum aestivum L.), or spring wheat-fallow, or winter wheat-fallow crop rotation. From 2009 to 2012, the study site followed either a continuous alfalfa (Medicago sativa L.) or a three-year crop rotation consisting of spring wheat in the first year followed by pea (Pisum sativum L.), and hay barley (Hordeum vulgare L.) in the second and third years, respectively [16,33]. In the spring of 2012, the entire experimental site was planted with glyphosate tolerant rapeseed (Brassica napus L.) and treated with herbicide. The rapeseed was tilled to a depth of 15 cm in July 2012 and planted in September 2012 following the experimental design described below.

Experimental Design
Our experiment followed a split-plot design, with cropping management systems as the whole plot variable and crop phases as the subplot variable. We assigned each whole plot randomly to one of three management systems: conventional no-tillage row crop production, tilled organic row crop production, and integrated crop-livestock organic row crop production. The conventional no-tillage used herbicides, while the tilled organic used tillage to terminate cover crops and control weeds (Supplementary Table S2). In grazed organic system, sheep (Ovis aries L.) were used to eliminate or reduce tillage intensity, terminate cover crops, and manage weeds. Both tilled and grazed organic systems received the USDA organic certification on July 11, 2015 making that year's crop harvest eligible for the organic label. Thus, this study comprises the three years of transition to organic period (July 2012 to July 2015) and the first two years (2016-2017) of the established organic systems. The agronomic management details are provided in Supplementary Table S2.
Each subplot (13 m × 90 m each), separated by a 1 m fallow strip, was assigned randomly to one of the five crops in the 2012-2013 growing season and then followed a 5-year rotation of crop phases sequentially as (Yr1) safflower (Carthamus tinctorius L.) under sown with yellow sweet clover [Melilotus officianalis (L.) Lam], (Yr2) yellow sweet clover cover crop, (Yr3) winter wheat (Triticum aestivum L.), (Yr4) lentils (Lens culinaris Medik), and (Yr5) winter wheat. Due to logistics, we compared the impacts of cropping systems beginning with the yellow sweet clover cover crop (Yr2) and winter wheat (Yr3) phases in all five years (2013-2017). Austrian winter pea (Pisum sativum subsp. arvense.) was planted in fall 2012 for the first year (2013) of the experiment because of the biennial nature of yellow sweet clover which would have required a previous year of seeding (Supplementary Table S2).

Ground Beetle Community Sampling
We estimated ground beetle activity density and community structure using pitfall traps, the most common method of sampling for ground dwelling insects [34]. In the center of one half of each subplot, we placed four pitfall traps arranged in a 2 m × 2 m rhombus pattern (Supplementary Figure S1). To construct the pitfall traps, we dug 10 cm × 12 cm holes using a post-hole auger, placed two stacked 0.5 L plastic cups (Solo Cup Co, Lake Forest, IL, USA) flush with the soil surface, and filled approximately 25% of the top cup with a propylene glycol-based antifreeze as a killing agent (Arctic Ban Antifreeze, Camco Manufacturing, Greensboro, NC, USA). The traps were open for a 48-h interval in every month from early spring (May) to late summer (September) for five consecutive years (2013-2017). We transferred all ground beetles caught in the pitfall traps from each subplot (n = 450) to plastic bags (Whirl-Pak ® , Nasco Inc., Fort Atkinson, WI, USA), cleaned in the lab, preserved in 70% ethyl alcohol, and identified to species following Lindroth [35]. The few ground beetles that we could not positively identify to species (<1% of the sample) were identified to genus and recorded as morpho-species. Nomenclature follows Bousquet [36].

Data Analysis
Metrics of ground beetle community structure included abundance (measured as activity density), species richness, α-diversity, and evenness. We calculated abundance as the number of individual sampled, species richness (S) as the total number of species present, and α-diversity as Simpson's Diversity Index (1-D) in each subplot: where p i is proportion of individuals of ith species in a community and S is the total number of species in that community [37,38]. Similarly, species evenness (J ) was calculated using Pielou's evenness [39]: where H is the Shannon-Wiener diversity index, Hmax = log(S), and S is the total number of species or species richness. The Shannon-Wiener diversity index is calculated as − s i p i log p i where p i is the proportion of each species [40].
We compared these metrics of ground beetle diversity among cropping systems, between crop phases, and through time using linear and generalized linear mixed effects models with a random intercept for year and block. Cropping system, crop phase, and sampling year were treated as fixed effects and block, year, and month were treated as random effects with month nested within year. Residuals were examined for homogeneity of variance and qq-plots were used to check for normality. The R packages 'lme4 and 'lmerTest' were used for analysis [41,42]. The standard lmerTest provides p-values for type III analysis of variance (ANOVA) for the models in lme4 via Satterthwaite's approximation of degrees of freedom method. Estimated marginal means and post-hoc tests were made using the R package "emmeans" [43].
Temporal differences in community structure of ground beetles across cropping systems and crop phases were assessed using non-metric multidimensional scaling (NMDS) ordination and permutational multivariate analysis of variance (PERMANOVA) on a Bray-Curtis dissimilarity matrix. Pairwise dissimilarity was computed using the Bray-Curtis index: where BC jk is the dissimilarity between communities j and k, a ij and a ik are the abundances of species i in sites j and k, respectively, and S is the combined total number of species in beetle communities [44,45].
To assess whether beetle communities were different among cropping systems and between crop phases across the five years of this study, we performed "adonis" test of PERMANOVA. Finally, we identified indicator species for cropping systems and crop phases using "IndVal" [46], a procedure to detect key species in a particular site based on their relative abundance (specificity) and relative frequency (fidelity) [47]. Indicator species analysis was performed to one-way (i.e., each system alone) and two-way combinations (i.e., systems paired) of farming systems. All statistical analyses were performed using R 3.2.4 [48]. Multivariate analyses and ordination graphics were conducted using vegan [49] packages of R. We produced all other graphics using the "ggplot2" and "sciplot" packages in R.

Results
We collected a total of 2929 beetle specimens over the five-year study period, representing 71 species; 1544 specimens and 60 species occurred in cover crop subplots, and 1385 specimens and 48 species were occurred in winter wheat subplots. Among the 60 beetle species that were sampled in cover crops, 20 (315 specimens) were collected in the conventional no-till system, 47 (674 specimens) in the tilled organic, and 42 (555 specimens) in the reduced-till grazed organic (Supplementary Table S3). Similarly, among the 48 beetle species that were sampled in winter wheat, 20 (168 specimens) were collected in the conventional no-till system, 34 (437 specimens) in the tilled organic, and 31 (780 specimens) in the reduced-till grazed organic (Supplementary Table S4).
In cover crops, the five most abundant beetle species sampled represented about 80% of the total capture in the conventional no-till system, 62% in tilled organic system, and 62% in reduced-till grazed organic system ( Figure 1; Supplementary Table S3). Similarly, in winter wheat, the five most abundant beetle species sampled represented about 70% of the total capture in the conventional no-till system, 63% in the tilled organic system, and 80% in the reduced-till grazed organic system ( Figure 1; Supplementary Table S4).
Agronomy 2020, 10, x FOR PEER REVIEW 5 of 16 relative frequency (fidelity) [47]. Indicator species analysis was performed to one-way (i.e., each system alone) and two-way combinations (i.e., systems paired) of farming systems. All statistical analyses were performed using R 3.2.4 [48]. Multivariate analyses and ordination graphics were conducted using vegan [49] packages of R. We produced all other graphics using the "ggplot2" and "sciplot" packages in R.

Results
We collected a total of 2929 beetle specimens over the five-year study period, representing 71 species; 1544 specimens and 60 species occurred in cover crop subplots, and 1385 specimens and 48 species were occurred in winter wheat subplots. Among the 60 beetle species that were sampled in cover crops, 20 (315 specimens) were collected in the conventional no-till system, 47 (674 specimens) in the tilled organic, and 42 (555 specimens) in the reduced-till grazed organic (Supplementary Table  S3). Similarly, among the 48 beetle species that were sampled in winter wheat, 20 (168 specimens) were collected in the conventional no-till system, 34 (437 specimens) in the tilled organic, and 31 (780 specimens) in the reduced-till grazed organic (Supplementary Table S4).
In cover crops, the five most abundant beetle species sampled represented about 80% of the total capture in the conventional no-till system, 62% in tilled organic system, and 62% in reduced-till grazed organic system ( Figure 1; Supplementary Table S3). Similarly, in winter wheat, the five most abundant beetle species sampled represented about 70% of the total capture in the conventional notill system, 63% in the tilled organic system, and 80% in the reduced-till grazed organic system ( and winter wheat (right panel) across three cropping systems (first row: conventional no-till, second row: reduced-till grazed organic, third row: tilled organic). Overall, cover crops and organic systems, compared to winter wheat and conventional no-till system, helped to increase evenness of beetle species (see Table 1 and text for details).
The effects of cropping system on activity density, diversity, and evenness, depended on year (Year × System; Table 1), but the overall trend was that these impacts were lower during the transition to organic period when compared to established organic systems (Table 1; Figure 2). The effects of crops on activity density, richness, diversity, and evenness of ground beetles also depended on year (Year × Crop) (Table 1). Overall, cover crop subplots had greater beetle richness, diversity, and evenness, but not activity density, than winter wheat subplots (Table 1; Figures 1 and 2; p < 0.01 in all cases). The effects of crops and cropping system on activity density, richness, diversity, and evenness did not depend on year (Year × Crop × System; Table 1).  In both cover crops and winter wheat subplots, the mean ground beetle activity density, diversity, and richness were greater in tilled and reduced-till grazed organic than in conventional notill, but generally the organic systems were not different from each other (Table 1; Figure 2). In cover crop subplots, the activity density was not different among cropping systems in the years 2013, 2016, Figure 2. Ground beetle activity density/sampling (A,B), species diversity (C,D), and species richness (E,F) in cover crops and winter wheat, across conventional no-till, reduced-till grazed organic, and tilled organic systems. Error bars represent the standard errors of the mean, while letters above bars indicate significant differences between systems at p < 0.05 within each year. Vertical black lines separate the transition to organic (to the left) and organic certification (to the right) phases.
In both cover crops and winter wheat subplots, the mean ground beetle activity density, diversity, and richness were greater in tilled and reduced-till grazed organic than in conventional no-till, but generally the organic systems were not different from each other (Table 1; Figure 2). In cover crop subplots, the activity density was not different among cropping systems in the years 2013, 2016, and 2017 (Supplementary Table S3 and S-I; Figure 2). Interestingly, beetle species richness and diversity were similar among cropping systems during transition to organic period (2013-2015), and these metrics were larger in established organic systems in 2016 and 2017 than in conventional no-till systems (Figure 2). Similarly, in winter wheat subplots, reduced-till grazed organic had greater beetle activity density in 2015, 2016, and 2017 and tilled organic had greater activity density in 2016 and 2017 than in conventional no-till ( Figure 2). Beetle richness and diversity were similar among systems during the transition to organic period (2013-2015), but greater in organic systems in 2016 and 2017 than in the conventional no-till after the organic certification was achieved (Figure 2).
The PERMANOVA based on activity density and species composition indicated that the community composition of ground beetle in the two studied crops did not differ among cropping systems during transition to organic period (Cover crops: Figure 3 . Indicator species analysis on beetle community matrix of five years' data showed that Apristus pugetanus was indicative of tilled organic system, but no species was indicative of either conventional no-till or reduced-till grazed organic systems alone (Supplementary Table S5). However, when the systems were paired, P. melanarius was indicative of conventional no-till and tilled organic whereas A. placidum, A. cupreolata, A. littoralis, Clivina fossor, Harpalus fuscipalpis, Poecilus corvus, and P. cursitor were indicative of both organic systems, but no species was indicative of conventional no-till and reduced-till grazed organic systems (Supplementary Table S5).
As opposed to the impact of cropping systems, the ground beetle community composition was different between crop types during the initial two years of our study, marginally different during the third and fourth years, and similar during the final year ( Indicator species analysis on beetle community matrix of five years' data showed that Amara apricaria, A. patruelis, Apristus pugetanus, Bradycellus congener, and P. scitulus were indicative of cover crops whereas only Loricera pilicornis was indicative of the winter wheat crop (Supplementary Table S5).
Agronomy 2020, 10, x FOR PEER REVIEW 8 of 16 pugetanus was indicative of tilled organic system, but no species was indicative of either conventional no-till or reduced-till grazed organic systems alone (Supplementary Table S5). However, when the systems were paired, P. melanarius was indicative of conventional no-till and tilled organic whereas A. placidum, A. cupreolata, A. littoralis, Clivina fossor, Harpalus fuscipalpis, Poecilus corvus, and P. cursitor were indicative of both organic systems, but no species was indicative of conventional no-till and reduced-till grazed organic systems (Supplementary Table S5). As opposed to the impact of cropping systems, the ground beetle community composition was different between crop types during the initial two years of our study, marginally different during the third and fourth years, and similar during the final year ( Figure 5

Discussion
To the best of our knowledge, this is the first study comparing the effects of tilled and reducedtill grazed transition to organic and established organic systems with a chemical-based no-till conventional system on ground beetle communities. With the goal of minimizing the negative effects of tillage, we used sheep to reduce tillage intensity, terminate cover crops, and manage weeds. Our study indicated that organically managed diversified systems, particularly with livestock integration, enhance ground beetle communities in the dryland crop production systems. Specifically, we observed the highest ground beetle activity density, species richness, species diversity, and species evenness in the reduced-till grazed organic, intermediate in tilled organic, and the lowest in no-till conventional systems. Also, we observed greater species richness and diversity of ground beetles in cover crop phases than in winter wheat. Hence, our study highlights the advantages of adopting ecologically-based agriculture such as organic cropping, integrated livestock, and cover cropping to enhance beneficial insects such as ground beetles that potentially provide ecosystem services of pest regulation.
Consistent with the previous studies conducted in the region [17,50,51], ground beetle communities in our study site were dominated by a small subset of species from the commonly occurring genera in the Northern Hemisphere such as Agonum, Amara, Carabus, Harpalus, and Pterostichus. In agreement with previous studies [51,52], but in contradiction to other studies [53,54], we found a greater activity density, richness, and diversity of ground beetles in established organic systems compared to those in no-till conventional, suggesting that in the semi-arid Northern Great Plains organic cropping systems provides diverse and suitable habitat conditions [55]. The

Discussion
To the best of our knowledge, this is the first study comparing the effects of tilled and reduced-till grazed transition to organic and established organic systems with a chemical-based no-till conventional system on ground beetle communities. With the goal of minimizing the negative effects of tillage, we used sheep to reduce tillage intensity, terminate cover crops, and manage weeds. Our study indicated that organically managed diversified systems, particularly with livestock integration, enhance ground beetle communities in the dryland crop production systems. Specifically, we observed the highest ground beetle activity density, species richness, species diversity, and species evenness in the reduced-till grazed organic, intermediate in tilled organic, and the lowest in no-till conventional systems. Also, we observed greater species richness and diversity of ground beetles in cover crop phases than in winter wheat. Hence, our study highlights the advantages of adopting ecologically-based agriculture such as organic cropping, integrated livestock, and cover cropping to enhance beneficial insects such as ground beetles that potentially provide ecosystem services of pest regulation.
Consistent with the previous studies conducted in the region [17,50,51], ground beetle communities in our study site were dominated by a small subset of species from the commonly occurring genera in the Northern Hemisphere such as Agonum, Amara, Carabus, Harpalus, and Pterostichus. In agreement with previous studies [51,52], but in contradiction to other studies [53,54], we found a greater activity density, richness, and diversity of ground beetles in established organic systems compared to those in no-till conventional, suggesting that in the semi-arid Northern Great Plains organic cropping systems provides diverse and suitable habitat conditions [55]. The divergence of the beetle communities after three years of initiating this study indicates that the three studied cropping systems acted as distinct ecological filters. A study done by Rivers et al. [56] also found an increase in ground beetle activity density and richness during their 3-year transition to organic experiment. However, previous studies were done either only for a short and transition to organic period [29,54,57] or in the established organic phase [30,58], which may not be representative of the temporal changes occurring during both transition to and established organic phases.
In accordance with previous studies [59], ground beetle species were more evenly distributed in organically managed crops than in no-till conventional ones, which were dominated by a few species including P. melanarius, A. cupreolata, and H. amputatus. In addition to these three and other omnivorous species (e.g., Agonum placidum, Amara cupreolata), organic systems also had predominantly granivorous species such as H. herbivagus, A. idahona, and A. quenseli and carnivorous species such as Bembidion quadrimaculatum, C. fossor, and, P. corvus [23] that were either rare or absent in the conventional system. Hence, the enhancement of species evenness in organic cropping has the potential for increasing biological control of pests [60], including mealworms, cricket, aphids, fly maggots, grubs, slugs, and lepidopteran eggs [21,59]. In addition to consuming insect pests [61], many ground beetle species are known to consume substantial amounts of weed seeds [62,63], having the potential to help regulate weed population [64,65]. As such, the increase of natural enemies observed in organic cropping systems [66,67] could help design ecologically-based insect pests and weeds management programs [68].
The ecology of many of the ground beetle species sampled in this study is largely unknown. Previous studies have shown that large bodied ground beetle species such as P. melanarius prefer undisturbed ground such as reduced-or no-till systems [54,69,70]. We observed P. melanarius to be associated with both conventional no-till and tilled organic systems, but less so with reduced-till grazed organic crops, suggesting that this species might prefer habitats that has no animal manures. Similarly, small-sized beetles such as Amara and Bembidion are known to be associated with tilled systems [54,70], which is not consistent with our results that small-bodied beetles were abundant in both till and reduced till organic systems. For example, while Apristus pugetanus has been reported in dry sand or gravel near banks [35], in our study it was strongly associated with reduced-till grazed organic, suggesting its affinity to animal manure and urine. Agonum placidum, Amara cupreolata, A. littoralis, C. fossor, P. corvus, and P. cursitor were sampled in the two studied organic systems and were rare in the conventional system. These species occur in cultivated soil with weedy vegetation [35], and the associations we observed could be due to the weedier characteristics of organic systems in Montana [51,71,72].
Previous studies showed that increased soil disturbance due to tillage reduces ground beetle habitat suitability [53,64] and increase the beetle mortality [73]. In accordance, the greater beetle activity density and richness observed in the reduced-till grazed organic system than in the tilled organic suggests that soil disturbance creates a less favorable environment for the overall beetle communities. In contrast with these observations, Hatten et al. [69] indicated that crop residues in no-till conventional fields may impede carabid movement, but tillage in organic fields reduces crop residues and facilitates beetle activity density. Additionally, Jowett et al. [74] found insignificant effects of tillage on ground beetle communities.
Consistent with previous studies [31,57,75], our study showed that cover crops support greater species richness and diversity of ground beetles than winter wheat, which could have positive effects on weed seed and insect pest regulations. Interestingly, as opposed to the effects shown by cropping systems, the ground beetle communities were different between crop phases during the initial two years of this study, but they were similar during the final three years. The exact reason of the similarity in established phases is unknown, but based on previous studies [29,56,76,77], cover crops, which were followed by winter wheat in our rotation scheme, could have legacy effects in affecting ground beetle communities in the subsequent crops. Cover crops provide ecosystem services such as soil erosion control, water-quality regulation, accumulation of soil organic matter, increase soil microbial biomass, increase in subsequent crop yield, colonization of arbuscular mycorrhizal fungi, and weed suppression [78,79], but the research on effects of cover crops on insect pests, natural enemies, and biological control were lacking [80]. Hence, while cover crops can help organic and conventional growers improve soil fertility and suppress weeds [79,81], our study shows that they can also play an important role in supporting ground beetles in semi-arid agroecosystems.
Our study found that A. apricaria, A. patruelis, Apristus pugetanus, Bradycellus congener, and Poecilus scitulus were indicative of cover crops, but ecological drivers of these associations are, to our understanding, unknown. Only one species, Loricera pilicornis, was indicative of wheat subplots. Previous studies have found that L. pilicornis prefers shade and soils rich in organic matter [36], and that it primarily consumes Collembola in cereal fields [22]. A previous study has shown a significantly lower Collembola in cover crops compared to control treatment plots [75] so, it is possible that the wheat subplots had more Collembola, than cover crops subplots but assessing underground arthropods was beyond the scope of our study.

Conclusions
Management practices such as tillage, cover cropping, crop rotation, and pesticides alter habitat suitability for ground beetles [82,83]. Our study suggests that, by enhancing generalist predators like ground beetles, ecologically-based farming such as diversified organic systems with integrated livestock and cover crops, can enhance conservation biological control and, ultimately, help to avoid chemical-based pest control tactics. Integrated crop-livestock systems can be a viable option not only to terminate cover crops, reduce soil erosion, and increase economic benefits to the producers [84,85], but, as shown in this study, to support beneficial insects. However, the benefits of integration of livestock into cropping systems may come with the cost of increased weed pressure [16], requiring more effective weed management practices. With the benefits of cover cropping in increasing soil health, suppressing weeds [79], and supporting beneficial insects, either a close collaboration of ranchers and growers for the grazing lease to terminate cover crops [84] or a local integration of cropping with livestock systems [18] could provide new approaches to increase the sustainability of the agricultural industry.
Supplementary Materials: The following are available online at http://www.mdpi.com/2073-4395/10/8/1210/s1: Figure S1: Experimental design showing three blocks (i.e., reps), plots (i.e., farming systems), and split plots (i.e., crops) at Fort Ellis Experimental center, Bozeman, MT. Table S1: Precipitation and temperature data for Fort Ellis Research and Extension Center in Bozeman, MT, USA. Table S2: Agronomic management details in cover crops and winter wheat across conventional no-till, reduced-till grazed organic, and tilled organic between 2013 to 2017 at Fort Ellis Research and Extension Center in Bozeman, MT, USA. Table S3: List of ground beetle species collected from Fort Ellis Experimental center, Bozeman, MT, 2013-2017 from cover crops across three farming management. Table S4: List of ground beetle species collected from Fort Ellis Experimental center, Bozeman, MT, 2013-2017 from winter wheat across three farming management. Table S5