Weeds Enhance Multifunctionality in Arable Lands in South-West of France

The current challenge in agriculture is to move from intensively managed to multifunctional agricultural landscapes that can simultaneously provide multiple ecological functions (multifunctionality), thus ensuring the delivery of ecosystem services important for human well-being. There is evidence that biodiversity is the main driver of multiple ecosystem functions. However, how biodiversity, and which components of biodiversity are the sources of multifunctionality, remain elusive. In the present study, we explore the role of weed richness and weed abundance as possible sources of ecosystem multifunctionality of an intensive agricultural landscape. Weeds are a key component of the arable field ecosystem trophic network by supporting various ecological functions while being a possible threat for production. We combine empirical data on ten ecosystem functions related to pollination, pest control and soil fertility, and measured across 184 fields cultivated with winter cereal, oilseed rape or hays in the Long Term Socio-Ecological Research site Zone Atelier Plaine & Val de Sèvre. We found that weed diversity was a strong contributor to multifunctionality in all crop types, especially when using the threshold-based approach. The effects of weed diversity were less pronounced for individual ecological functions except for weed seed predation and urease activity. As weeds may have dual effects on yields, we also explored the relationship between ecosystem multifunctionality and yield considering weed abundance. We however found a neutral relationship between yield and ecosystem multifunctionality. These results suggest that field management that maintains high levels of weed diversity can enhance multifunctionality and most ecological functions. Understanding how to maintain weed diversity in agricultural landscapes can therefore help to design sustainable management favoring the delivery of multiple services while maintaining food production. The next challenge will therefore be to assess the relative contribution of management practices, landscape features and weed diversity on ecosystem multifunctionality and yield.

The current challenge in agriculture is to move from intensively managed to multifunctional agricultural landscapes that can simultaneously provide multiple ecological functions (multifunctionality), thus ensuring the delivery of ecosystem services important for human well-being. There is evidence that biodiversity is the main driver of multiple ecosystem functions. However, how biodiversity, and which components of biodiversity are the sources of multifunctionality, remain elusive. In the present study, we explore the role of weed richness and weed abundance as possible sources of ecosystem multifunctionality of an intensive agricultural landscape. Weeds are a key component of the arable field ecosystem trophic network by supporting various ecological functions while being a possible threat for production. We combine empirical data on ten ecosystem functions related to pollination, pest control and soil fertility, and measured across 184 fields cultivated with winter cereal, oilseed rape or hays in the Long Term Socio-Ecological Research site Zone Atelier Plaine & Val de Sèvre. We found that weed diversity was a strong contributor to multifunctionality in all crop types, especially when using the threshold-based approach. The effects of weed diversity were less pronounced for individual ecological functions except for weed seed predation and urease activity. As weeds may have dual effects on yields, we also explored the relationship between ecosystem multifunctionality and yield considering weed abundance. We however found a neutral relationship between yield and ecosystem multifunctionality. These results suggest that field management that maintains high levels of weed diversity can enhance multifunctionality and most ecological functions. Understanding how to maintain weed diversity in agricultural landscapes can therefore help to design sustainable management favoring the delivery of multiple services while maintaining food production. The next challenge will therefore be to assess the relative contribution of management practices, landscape features and weed diversity on ecosystem multifunctionality and yield.
Keywords: bee richness, carabids, haylands, pest control, pollination, soil organic carbon, soil enzyme activity, yield INTRODUCTION Knowledge of how biodiversity contributes to ecosystem functioning at multiple scales is critical to conserving, managing and restoring multifunctional landscapes, especially because the capacity of ecosystem to maintain multiple processes has been related to stability of ecosystems (Huang et al., 2019). Biodiversity is one of the main drivers of ecosystem functioning, along with land use or soil conditions (Lavorel et al., 2011;Balvanera et al., 2014). In general, a positive asymptotic relationship between biodiversity and single ecosystem function has been found (Isbell et al., 2011;Liang et al., 2016), suggesting that a tiny fraction of the species pool is necessary to support an individual function (Slade et al., 2019). However, to maintain the entire set of functions and services simultaneously, far more species are required (Isbell et al., 2011;Sasaki et al., 2019). There is evidence that multidiversity, i.e., the diversity of a variety of taxa, enhances multifunctionality (Wang et al., 2019). Indeed, different species enhance ecosystem functioning during different years, at different places, for different functions and under different environmental contexts (Isbell et al., 2011). However, the source of multifunctionality, i.e., how biodiversity, and which components of biodiversity are the sources of multifunctionality, remain elusive.
Weeds are a key component of the trophic network in arable fields, supporting various functions such as pollination, pest control or soil fertility (Marshall et al., 2003;Nicholls and Altieri, 2013;Bretagnolle and Gaba, 2015). Multifunctionality is here defined as ecosystem function multifunctionality (Manning et al., 2018), which has been widely studied in grasslands (Sasaki et al., 2019;Wang et al., 2019), and more recently in other ecosystems, such as forests (Zavaleta et al., 2010;Xie et al., 2018;Huang et al., 2019), soils (Valencia et al., 2018), and farmland habitats (Rallings et al., 2019;see Hölting et al., 2019 for a review). Agricultural land represents ∼40% of the Earth's land surface and therefore a better understanding of the drivers of multifunctionality as well as the relationship between multifunctionality and crop production may provide fill knowledge gaps to develop sustainable management strategies.
Species richness has been the main metric analyzed in regard to the role of biodiversity in multifunctionality (Huang et al., 2019). But in addition to species richness within communities, the functional properties of the most locally abundant species may drive ecosystem functioning following the mass-ratio hypothesis (Grime, 1998), or at least may provide a limited number of functions (Gamfeldt et al., 2008;Isbell et al., 2011). Locally rare species, alternatively, have been shown to be important diversity component for preserving high levels of ecosystem multifunctionality . Therefore, whether locally abundant species, locally rare species or both are required to maintain multifunctionality is still unclear. In addition, multifunctionality may be quantified by various methods [review in (Hölting et al., 2019)]: the most common consists in aggregating the multiple assessed functions into a single metric [>80% of studies reviewed in (Hölting et al., 2019)], most often by using the threshold approach (Gamfeldt et al., 2008) though the averaging approach accounts for 30% of studies.
In the present study, for the first time to our knowledge, we explore the role of weed diversity as a possible source of multifunctionality in an intensive agricultural landscape. To do so, we examine the relationship between weed diversity and ecosystem multifunctionality in 184 production fields cultivated with an arable crop or a hay in the Long Term Socio-Ecological Research site (LTSER) Zone Atelier Plaine & Val de Sèvre in South-West France. To understand how weed diversity may affect ecosystem multifunctionality, we explore the relationships between weed diversity and ecosystem multifunctionality using two different components of weed diversity (weed richness and weed total abundance) and two approaches to measure ecosystem multifunctionality (i.e., the averaging approach and the threshold approach). We also examine to which extent ecosystem multifunctionality is related to either locally abundant weed species or locally rare ones. In each field, biodiversity and various ecological functions [pollination, pest control (weed seed and aphid predation) and soil fertility (total organic carbon and enzyme activities involved in carbon, nitrogen, sulfur and phosphorus cycles)] were estimated to assess the effects of weed species diversity on ecosystem functioning. Finally, as crop production is the main service in agricultural landscapes, we investigate the relationship between crop production and ecosystem multifunctionality. Because a higher ecosystem multifunctionality suggests the provision of more ecological functions related to pest control, pollination and soil activities, we expect higher productivity in fields with higher ecosystem multifunctionality.

Study Site
The study site area, LTSER "Zone Atelier Plaine & Val de Sèvre" [hereafter ZA PVS; (Bretagnolle et al., 2018b)], covers 435 km² of agricultural land in the south of the city of Niort, in the Deux-Sèvres department in the Nouvelle-Aquitaine Region, France. The most common crops are winter cereal (33.8%), maize (9.6%), sunflower (10.4%), winter oilseed rape (8.3%), and haylands (13.5%, including mainly forage temporary crops such as alfalfa and to a lesser extent permanent grasslands). The field size is 4.5 ha in average.

Biodiversity Surveys
Biodiversity (weeds, ground beetles and bees) was surveyed in 184 fields (i.e., the scale of management decisions in farmland) including 78 winter cereals (winter wheat, winter barley) and 45 winter oilseed rape fields, and 61 haylands (i.e., temporary grasslands used as forage crops). For convenience and because most of haylands are temporary crops (generally 2-5 years, mainly alfalfa), we refer as crop types when referring to these three types of cultivated plants. All fields sown with winter cereals and winter oilseed rape were tilled before sowing. Weed species identities and occurrences were recorded in 80, 0.25 m² quadrat positioned in the core field along two transects of 10 quadrats after spring herbicide application (mid-April-end of May). Weed abundance was computed as the sum of each species occurrence in the 80 0.25 m² quadrats in each field. Ground beetles (carabid) diversity and abundance were assessed using two pitfall traps placed in the margin (first crop row) and two others at 10 and 25 m from the edge. Bee diversity was monitored in the same fields using colored six pan traps per field left for 4 days in the fields. Figure 1 depicts the survey design. Details on all protocols are available in Bretagnolle et al. (2018a).

Estimation of Ecological Functions
Pollination by insects was measured using two oilseed rape phytometer plants (Brassica napus sp.), one placed in the field edge and one in the field core. Phytometers were grown in an insect-proof greenhouse in the lab. At flowering, they were left for 4 days in fields and then placed again in the insect-proof greenhouse in the lab. For each phytometer, we assessed the fruiting success of the flowers that were opened during the 4 days in the fields. Here, we used the average fruiting success per field to estimate pollination (Perrot et al., 2018).
Biological pest control was quantified in each field using two types of predation cards: cards with weed seeds (Viola arvensis) and cards with aphids (Acyrthosiphon pisum Harris.). Fourteen cards were placed during 4 days in the first 10 m in each field along two transects. Predation rates were estimated by dividing the number of seeds that were predated over the number of seeds on the cards.
Soil cores were sampled to quantify soil carbon and soil enzyme activities, between the end of March and early April 2016. Five soil cores were collected in the topsoil layers (0-15 or 0-30 cm depending on the soil depth, i.e., the total soil depth never exceeds 30 cm in the study area). Position of cores was random in the field. The core samples were pooled per field and stored at ambient temperature until the analyses were performed. Soil samples were air-dried for estimation of total organic carbon (C) content that was measured by dry combustion. Fresh soils were used to measure oil enzyme activities involved in carbon (C), nitrogen (N), sulfur (S), and phosphorus (P) cycles. Arylamidase, β-glucosidase, urease, arylsulfatase and phosphatase were quantified by the platform Biochem-Env (Cheviron et al., 2018) using colorimetric methods according to the ISO (ISO, 2018) standard, with a slight modification for urease. All measurements were performed at the soil pH, in an unbuffered soil water solution (Lessard et al., 2013), except arylamidase, performed in Tris 100 mM pH 7.5 (Acosta- Martínez and Tabatabai, 2000).
Yields and the amount and type of fertilizers applied during the cropping season were collected by means of farmers' interviews. Data were collected for 100 fields, because some farmers did not accept face-to-face interviews. From these surveys, we derived the amount of N input. The quantity of inorganic nitrogen used was directly calculated from the fertilizer composition and the quantity applied, and the quantity of nitrogen mineralized in organic fertilizers was calculated using the method described by Jeuffroy and Recous (1999).
The general statistics of biodiversity, ecological functions and agricultural practices are given in Supplementary Material, Table S1.

Statistical Analysis
We first compared, among crop types, the magnitude of each individual ecological function, the number of functions achieved above each threshold, and ecosystem multifunctionality (EF-multifunctionality; Manning et al., 2018) using ANOVA. We used EF multifunctionality (also referred to averaging approach) which summarizes ecosystem multifunctionality and reflects the change in the average level of a bundle of ecosystem functions (Barnes et al., 2014). EF-multifunctionality per field is the sum of the standardized values of all the ecosystem functions. To avoid an overweight of the bee richness, soil carbon content and enzyme activities, values of these functions were standardized to range between 0 and 1. Pollination and pest control estimates, being rates, already vary between 0 and 1, and were thus not standardized. We also evaluated whether multiple functions are simultaneously performing at high levels using the threshold approach (Byrnes et al., 2014). This approach consists in computing the number of functions achieved above a given threshold. We selected six thresholds which corresponds to the percentages (i.e., 20, 30, 40, 50, 60, and 70%) of the maximum observed value of each function (Byrnes et al., 2014). If the number of functions greater than a threshold is always lower than the total number of functions, this highlights a tradeoff between functions i.e., one function is always maximized when another is minimized. For all tests, we checked for the prerequisites of homogeneity of variances using Bartlett test, normality using Shapiro test, and applied Kruskall-Wallis rank test when necessary. We used Tukey's post hoc tests to assess differences among crop types.
We then tested whether weeds (richness or abundance) would promote each individual ecological function, the number of functions achieved above each threshold, and EFmultifunctionality. We used ordinary least squares regression to determine how the various weed metrics influenced the dependent variables (each individual function, the number of functions achieved above each threshold, and EF-multifunctionality). We also included crop type and its interaction with weed richness or abundance. For pest control (aphid and seed predation rates), we also included ground beetle diversity and abundance as explanatory variables with, respectively, weed diversity and weed abundance. Two models per explanatory variable were built and compared to identify whether weed richness or abundance better explained each function, the number of functions achieved above each threshold and EF-multifunctionality. A model selection procedure based on Akaïke criterion (AIC) was performed to select the model that best explained the data (Burnham and Anderson, 2002).
To examine how EF-multifunctionality and each ecological functions were related to abundant vs. rare weed species, we built ten competing models (five with the number of abundant species and five with the number of rare species) and quantified the percentage of variance explained by the number of locally abundant species or of locally rare ones. The number of abundant species was arbitrary defined as the number of species over quantile 50, 60, 70, 80, and 90% of the distribution of the plant abundance per field whereas the number of rare species was in symmetry the number of species lower than quantile 50, 40, 30, 20, and 10% of that distribution. Models include crop type and its interaction with either the number of abundant species or the number of rare species.
Finally, we examined the relationship between EFmultifunctionality and productivity, i.e., crop yield in cereal and oilseed rape fields and aboveground biomass production in grasslands. We built a complete linear model with EFmultifunctionality, crop type, field area, weed abundance and the amount of N input (kg/ha) as covariate. We considered the amount of N input because it is usually positively related to productivity, and use field area to account for the discrepancy in the spatial coverage between productivity (measured at field scale) and ecosystem function (part of the field). To account for potential interactive effect between EF-multifunctionality on one hand and weed abundance, crop type and field area on the other hand, we included these three two-way interactions. Finally, because weed assemblies and the amount of N input vary among crops, we included the interaction between weed abundance and crop type, and between the mount of N input and crop type. Among the fields surveyed, 21 fields (8 winter cereals and 13 haylands) were managed as organic farming. Yields and weed abundances were standardized per crop type and farming system using z-scores. This transformation does not constrain the variability found in the raw data and allows focusing on each effect independently of the crop and farming system effects.
Analyses and plots were performed using packages MASS (Venables and Ripley, 2002) and ggplot2 (Wickham, 2009) on R v. 3.5.1 (R Core Team, 2018). R code are available in Supplementary Text.

Multifunctionality in Haylands and Annual Crops
EF-multifunctionality was higher in haylands [average 4.79 (sd ± 1.13)] than in annual crops (winter cereal: 4.42 ± 0.83 and oilseed rape: 4.08 ± 1.03; Figure 2A), although not significantly. No significant differences in bee diversity and oilseed rape fruiting success (a proxy of insect pollination) were found among crop types (Figures 2D,E). In contrast, pest control was significantly higher in haylands (weed seed predation rate: 0.82 ± 0.13 and aphid predation rate: 0.73 ± 0.22), and to a lesser extent in winter cereal (0.61 ± 0.21 and 0.64 ± 0.2) than in oilseed rape (0.56 ± 0.22 and 0.45 ± 0.18) (Figures 2B,C). A similar pattern was observed when considering belowground ecological functions: soil carbon and soil enzyme activities related to C, N, P, and S cycles were significantly higher in haylands compared to annual crop fields and no differences were observed among annual crops ( Figures 2F-K).

Weed Abundance Rather Than Weed Diversity Sustains Ecological Functioning
Weed richness and abundance significantly varied among crop types with highest richness and abundance in hay (21.5 ± 9.8 species and 308 ± 128 plants/m²) than in winter cereal (12.3 ± 9.75 species and 153 ± 176 plants/m²) and in oilseed rape (16.6 ± 4.26 species and 153 ± 83 plants/m²) fields, a difference that could not be accounted for by sown species in hay (e.g., alfalfa, ray-grass, Festuca spp.) that were withdrawn from analyses. The most common species were Epilobium tetragonum, Cirsium arvense and Poa pratensis in winter cereal, Epilobium tetragonum, Lolium multiflorum and Trifolium pratense in oilseed rape and Dactylis glomerata, Poa trivialis and Rubus fruticosus in haylands. The thresholdbased approach to multifunctionality revealed strong positive effects of weed abundance whatever the crop type (Figure 3A), especially for thresholds higher than quantile 30 ( Figure 3B; Table 1B). Weed abundance was also a strong contributor Frontiers in Sustainable Food Systems | www.frontiersin.org FIGURE 2 | Among crop type variation of EF-multifunctionality (A) and each ecological function (B-K). Soil enzyme activities are expressed in mU.g −1 dry soil, representing nanomoles of product released per minute and per g of equivalent dry soil. Significant differences are indicated by *(P < 0.05), **(P < 0.01), and ***(P < 0.001). NS indicates a non-significant effect.
to aboveground EF-multifunctionality (i.e., considering pest control, bee diversity, and pollination success; Table 1B). EFmultifunctionality tended to increase with both weed abundance and richness in all crop types ( Figure 3A; Table 1). The effects of weed richness on multifunctionality thresholds and aboveground EF-multifunctionality were less pronounced than those of weed abundance (Supplementary Table S1) and varied among crop types being positive in winter cereal and oilseed rape and slightly negative in hay (Supplementary Figure 1).
When considering each individual function, only slight differences in the goodness-of-fits were found between weed abundance or weed richness except for weed predation rate and urease activity (Supplementary Table S2). Weed abundance strongly and positively influenced aphid predation rate and soil enzyme activities (arylsulfatase, β-glucosidase, and urease), while weed richness affected weed seed predation rate and arylamidase activity ( Table 2). Weed diversity had no significant effects on bee diversity, oilseed rape fruiting success, and phosphatase activities. However, relationships between weed abundance or richness and each individual function varied strongly among crop types and between functions (Figure 4). In haylands, weed diversity (either abundance or richness) had a negative effect on soil carbon, and most enzyme activities, and a positive one on aphid and weed seed predation rates. In winter cereal fields, weed richness had a significant negative effect on weed seed predation and soil carbon. Finally, in oilseed rape, weed diversity generally had a positive effect on weed seed predation rate, soil carbon and functions related to soil activity, and a negative one on aphid predation (Figure 4).

The Contribution of Rare Weed Species to Multifunctionality
The relative contributions of rare vs. abundant species to multifunctionality or the single functions, varied with the quantile considered to classify the species. The highest explanatory power of EF-multifunctionality variance (i.e., % of variance explained) was observed when the number of rare and abundant species was defined based on respective quantiles 40 and 60% (Supplementary Figure 2). While we observed a higher contribution of rare weed species FIGURE 3 | Weed abundance effects on EF-multifunctionality (A) and across functional thresholds (B) measured as the number of functions achieved in each field above thresholds (T), where T is the quantile-based ranking of each function across all fields. EF-multifunctionality was computed using the averaging approach. Lines represent the predicted relationship from the statistical models and shade areas the 95% confidence interval. Colors indicate the different crop type (cereals in red, oilseed rape in green and haylands in blue). Weed abundance is log-transformed.
to EF-multifunctionality (35.2% of variance explained) than of abundant ones (30.1%; Figure 5). This suggests the importance of the threshold used to define species groups on the results. Among the 4.15 (±3.78) rare weed species, the most frequent species were Medicago lupulina, Galium aparine and Vicia sativa in winter cereal, Galium aparine, Helminthotheca echioides, Geranium rotundifolium in oilseed rape and Anthriscus caucalis, Achillea millefolium, and Geranium pusillum in haylands. The provision of aboveground ecological functions was also related to the rare weed species, although the difference was relatively low. An exception was observed for bee diversity which was strongly explained by abundant weed species. Surprisingly, whatever the threshold for defining rare and abundant species, neither abundant nor rare weed species had high explanatory power for belowground ecological functions (soil carbon content and the enzyme activities); the differences was mostly explained by crop types.

Neutral Relationship Between Productivity and Multifunctionality
Contrary to our expectation, productivity (yield in annual crops and biomass production in haylands) was not higher in fields with high EF-multifunctionality. Rather, we observed a neutral relationship between EF-multifunctionality and productivity ( Figure 6A; Table 3). The amount of N, crop type and field area did not affect productivity either. But most importantly, we did not find any significant effect of weed abundance on productivity ( Figure 6B; Table 3). Overall, the variables selected for the productivity model (Table 3) resulted in a model with a low explanatory power (R² = 22.07%, P = 0.098).

DISCUSSION
In this study, we quantified the role of weeds as a potential source of multifunctionality in real farming conditions in an intensive agricultural landscape. The relationship between weed abundance or diversity, and ecosystem multifunctionality was evaluated in three crop types, cereal, oilseed rape, and hay. We considered above-and below-ground functions related to regulation services. The magnitude of ecosystem multifunctionality and of most individual ecosystem functions were higher in haylands than in annual crop fields. We found that weed diversity was a strong contributor of ecosystem multifunctionality especially when considering the simultaneous TABLE 1 | F-Statistics for the effects of weed richness (A) or weed abundance (B), crop type and their interaction on ecosystem multifunctionality computed as the sum of all standardized ecological functions (EF-MF), aboveground ecological functions (EF-MF aboveground), belowground ecological functions (EF-MF belowground) and using a threshold approach (with threshold from 20 to 70%). Significant effects are indicated as follows: *P < 0.01, **P < 0.001 and ***P < 0.0001. R 2 are also given for each model. Significant effects are indicated as follows: *P < 0.01, **P < 0.001, and ***P < 0.0001. Abbreviations for soil enzyme activities correspond to arylamidase (ARM), β-glucosidase (GLU), urease (URE), arylsulfatase (ARS) and phosphatase (PHOS).

EF-MF EF-MF aboveground
delivery of functions (threshold approach) and when considering only aboveground ecological functions. Weed diversity had also an effect in the provision of most the individual ecosystem functions under study here, with the strongest effect observed on weed seed predation and urease activity. Interestingly, across crop types, rare species had a higher contribution to ecosystem multifunctionality than abundant weed species at a given threshold (i.e., % quantile) used to defined the species groups. Finally, we found that contrary to our expectation, fields with higher ecosystem multifunctionality did not show higher productivity in terms of yield or biomass. Rather, we found a neutral relationship between productivity and ecosystem multifunctionality.
Our finding that biodiversity correlates with ecosystem multifunctionality extends past work showing a positive relationship between plant diversity and ecosystem multifunctionality, especially in experimental studies where such positive relationship has been repeatedly found (Wardle, 2016). However, in naturally assembled communities, as the ones studied here, there is much less evidence and the picture is less clear (Allan et al., 2015;Soliveres et al., 2016;Pennekamp et al., 2018). Our results show that the positive relationship between weed diversity and ecosystem multifunctionality was even stronger when considering aboveground ecological functions. High weed diversity contributed to the regulation of pests by increasing weed seed and aphid predation rates. Fields with FIGURE 4 | Weed abundance or richness effect on each ecological function in winter cereals (red), oilseed rape (green) and haylands (blue). Ecological functions are: aphid predation rate (A), weed seed predation rate (B), bee richness (C), OSR fruiting success (D), soil carbon (E), arylsulfatase (F), phosphatase (G), arylamidase (H), β-glucosidase (I) and urease (J). Lines represent the predicted relationship from the statistical models. The black line shows the relationship across crop types and shade area the 95% confidence interval. Dotted lines indicate non-significant relationships. Weed abundance is log-transformed. Soil enzyme activities are expressed in mU.g −1 dry soil, representing nanomoles of product released per minute and per g of equivalent dry soil.
high weed diversity may shelter more pest natural enemies. Indeed, predator and parasitoïd species richness are generally positively related to plant species diversity (Scherber et al., 2010;Leles et al., 2017;Schuldt et al., 2019). More pest natural enemies may enhance pest regulation by improving predator (parasitoïd)-prey interactions (Letourneau et al., 2009). Such increase may also result from a higher pest mortality, or through a synergistic interaction between predators in their ability to suppress pest populations, with a mortality greater than the summed mortality caused by each natural enemy species on its own (Barbosa and Castellanos, 2005). Evidence of synergistic interaction has indeed been found between foliar-foraging and ground-foraging pest natural enemies in the suppression of pea aphid populations in alfalfa (Losey and Denno, 1998). In our study, we only sampled ground beetles and found that their richness and abundance increased with weed diversity (Supplementary Figure 3). Ground beetles have been suggested to be natural enemies of both aphids (Firlej et al., 2013) and weed seeds (Bohan et al., 2011), the consumption rate of the latter being variable with both ground beetle and weed species . Determining whether the increase of weed diversity promotes a higher diversity of pest natural enemies and whether this higher diversity results in synergistic or additive effects among these natural enemies should be investigated in future research through a trophic network approach allowing the quantification of the intensity of predation rate per guild of natural enemies.
Surprisingly, although a high percentage of variance in bee richness was explained by rare or abundance weeds, neither a higher weed richness nor a higher weed abundance promoted bee richness or improved fruiting success of OSR phytometers. A positive relationship was expected between weed diversity and bee richness, since weeds are important resource for insect pollinators (Rollin et al., 2013;Requier et al., 2015). Higher plant species richness should enhance pollinator richness because of plant species specific pollinator preferences and a greater temporal and spatial availability of pollen and nectar resources (Potts et al., 2003;Fontaine et al., 2006;Ebeling et al., 2008). Higher weed abundance should also increase bee richness by providing a higher resource abundance and increasing bee attractiveness (Papanikolaou et al., 2017). Rather we found neutral relationships between weed richness or abundance and the number of bee species and OSR fruiting success whatever the crop type. In oilseed rape fields, bee monitoring was performed during oilseed rape flowering period (from April to May 2016). Oilseed rape is a mass-flowering crop offering a resource bloom during its flowering period (Requier et al., 2015), and weed floral resource may be less attractive for bees at this period. A different mechanism may result in neutral relationship in cereals and hay fields: the increase in the number of species of plants may not FIGURE 5 | Percentage of variance on EF-multifunctionality and each ecological functions explained by the number of locally abundant (Blue) and rare (Green) weed species which were defined using the number of species with an plant abundance <60% quantile and lower than 40% quantile of the weed abundance distribution.
be sufficient to attract bees because such crops do not provide floral resource either in quantity, quality, or both. In our study, we however found that the number of abundant weed species explained a high proportion of bee diversity (>50, Figure 5). Weeds differ in their melliferous potential (Baude et al., 2016;Ion et al., 2019). Consequently, a functional approach using traits related to pollination (e.g., flower resource), pest control (e.g., seed energetic content) and soil activities functions [root traits; but see Gaba et al. (2017) for other examples] may be useful to go deeper in the relationships. This more comprehensive and quantitative assessment would require to compile existing databases on traits because there is currently a lack of data for many functional traits especially for rare species (but see Bourgeois et al., 2019).
Similarly, to bee richness, we failed to find a positive relationship between weed richness or abundance, and fruiting success of OSR phytometers. Such relation would have suggested an indirect positive effect of weed richness or abundance through an increase of bee richness. Actually very few studies investigated the relationship between plant diversity and pollination success [only 8 according to van der Plas (2019)], but all found a neutral relationship, as here. (van der Plas, 2019) suggests that this absence of relationship may reflect a more important effect of pollinator diversity for pollination success than the one of plant diversity. We found a negative relationship between bee richness and OSR fruiting success (Supplementary Figure 4). In our study area, a previous study revealed that honeybee and Lasioglossum sp. abundances increased OSR fruiting success in oilseed rape fields (Perrot et al., 2018). The use of phytometers here to estimate OSR fruiting success may explained this discrepancy. Honeybees are assumed to be strong competitor, and exploitative competition between honeybees and wild bees was recently evidenced (Henry and Rodet, 2018), though relationship between honeybees and wild bees were found to vary from being positive, negative or neutral, depending on the plant species (Nielsen et al., 2012). Such relationships further affect pollinator community structure, resulting in a distortion in plant-pollinator interactions (Diekötter et al., 2010). In our study, we could not explore the link between weed diversity, OSR fruiting success, bee richness and honeybee abundance because bee richness was estimated using pan traps which poorly estimate honeybee abundance (Westphal et al., 2008). Further research therefore needs the use of complementary methods, for instance pan traps and sweep nets (Westphal et al., 2008), to assess bee diversity and hence disentangle the complex interactions between weeds, bees, and pollination in agricultural landscapes [but see (Bretagnolle and Gaba, 2015) for a conceptual framework].
The storage of organic carbon in soils represents a key function of soils that is critical in mitigating climate change effects (Paustian et al., 2016). Positive relationships between plant diversity and organic soil carbon stocks have been shown in experimental studies (Cardinale et al., 2012) and in most of natural communities (Gamfeldt et al., 2013;Maestre et al., 2016) although some studies found negative or neutral relationships (van der Plas, 2019). In our study, aboveground ecological functions appeared to be affected by the composition and structure of weed assemblies, but not by individual species as locally abundance and rare species poorly explained their variations. Weed richness and abundance equally explained the variation in soil carbon contents in the three crop types. Higher soil carbon content in more rich and abundant weed communities may influence soil microbial activity. This is supported by the positive relationships between weed abundance and the activity of β-glucosidase and urease, those enzymes being involved in C (β-glucosidase) and N (urease) cycles. In contrast, the activity of phosphatase (an enzyme involved in P cycle) was not affected by weed diversity. These results suggest that different mechanisms affect soil enzyme activities. Higher activity of β-glucosidase in abundant weed communities suggests greater availability of cellulose for soil microorganisms in high abundance weed communities. These communities may produce more root biomass resulting in a higher root exudation (Paterson FIGURE 6 | Relationship between productivity (standardized) and (A) EF-multifunctionality or (B) weed abundance. Colors indicate the different crop type (cereals in red, oilseed rape in green and haylands in blue) and the shade area the 95% confidence interval.
TABLE 3 | Statistics for the effects of EF-multifunctionality, crop type, amount of N input, weed abundance, field area and the interaction between EF-multifunctionality and crop type EF-multifunctionality and weed abundance, EF-multifunctionality and field area, as well as the interaction between weed abundance and field area and crop type and N input on productivity (yields in cereal and oilseed rape fields and aboveground biomass production in haylands). Productivity and weed abundance are centered and reduced using z-score per crop type and farming systems. et al., 2007) which is known to affect the activity and composition of soil microbial community (Shahzad et al., 2015) resulting in an increase of soil carbon uptake (You et al., 2014). Weed abundance may also increase enzyme activities by providing a greater amount of enzyme substrate in soil (Geisseler and Horwath, 2009). A high number of studies has investigated the relationship between multifunctionality and productivity (i.e., biomass production and yield), mostly in experimental conditions (Wardle, 2016). While positive relationships are often expected, we were not able to detect such positive effect of ecosystem multifunctionality and productivity whatever the crop type. In contrast, our analysis revealed a neutral relationship between ecosystem multifunctionality and productivity. Various parameters affect cereal and OSR yield and biomass production in haylands. Although we included nitrogen input, we may have missed other important variables such as soil properties (i.e., organic matter) or pest control management. Indeed, in this study site, a recent study showed that agrochemical applications overall accounted for about 24% of the variance of the OSR yield (Catarino et al., 2019). The neutral relationship may also arise from a compensation between benefits to yield provided by ecological functions sustained by weeds such as pest control, and competition between crop and weeds resulting in yield loss. Indeed, although we detected a negative relationship between weed abundance on productivity, this was not significant. Ecosystem multifunctionality was assessed using the averaging approach. Similar results were obtained with the threshold approach (data not shown). In both approaches, ecological functions had similar weights on productivity. However, it is likely that the different ecological functions benefit differentially to productivity (both in strength and type of relation e.g., linear, saturated). An assessment of ecosystem multifunctionality in which ecological functions are weighted according to their benefit to productivity may give a better picture of the relation [but see Allan et al. (2015) for an example in hay]. However, this requires a deeper understanding of the relationship between each ecological functions and productivity, and of the interactive effects of a suite of ecological functions on productivity. Exploring the effects of specific ecological functions while considering the others (i.e., taking account the complex interaction), for instance using structural equation modeling, could be relevant way to assess the effect of ecosystem multifunctionality on productivity. This would also allow taking account environmental factors (e.g., soil properties) and management practices.

Sum
Finally, our study emphasizes the critical role of haylands, and presumably grasslands in general, for maintaining ecosystem multifunctionality in agricultural landscape. Ecosystem multifunctionality was higher in haylands than in winter cereal and oilseed rape fields. Except for bee richness and fruiting success of OSR phytometers, the amount of ecological functions was greater in haylands. Haylands showed higher plant diversity than arable crop fields, and are less intensively managed fields, with no tillage and a low agrochemical inputs. There is strong evidence that intensive management has substantially altered biodiversity and affected ecosystem. The next challenge will be to assess the relative contribution of management practices and plant diversity (both being closely related) on ecosystem multifunctionality and productivity, to understand how diversifying agroecosystems by introducing meadows in landscapes and in crop rotation sequences could efficiently improve their functioning and sustain both the provisioning of ecological functions and yields.

DATA AVAILABILITY STATEMENT
The datasets generated for this study are available on request to the corresponding author.

AUTHOR CONTRIBUTIONS
SG and VB designed the study. NC and SP conducted the soil and enzyme activities analyses. TP and J-LG performed the phytometers pollination experiments. SG prepared the data and conducted the statistical analysis and wrote the first draft of the manuscript. All authors substantially contributed to revisions.

FUNDING
This paper was produced with the support of CESAB-FRB as part of the activities of the DISCO-WEED Working Group. The project was also supported by the French Ministry of Ecology (project 2017-2020 Pollinisateurs), the INRAE MP ECOSERV (BIOSERV project) and the ANR IMAGHO project (ANR-18-CE32-0002). Biochem-Env is a service of the Investment d'Avenir infrastructure AnaEE-France, overseen by the French National Research Agency (ANR) (ANR-11-INBS-0001).

ACKNOWLEDGMENTS
We would like to express our thanks to all field workers for their help with weed surveys, insect identification, pollination experiment and soil sampling. We sincerely thank the farmers of the LTSER Zone Atelier Plaine & Val de Sèvre for their involvement in our research programs. The authors are also grateful to C. Mougin and F. Martin for comments on the final version of the manuscript. We also wish to thank the editor and the two reviewers for valuable comments throughout the reviewing process.