Impact of Fertilizing Pattern on the Biodiversity of a Weed Community and Wheat Growth

Weeding and fertilization are important farming practices. Integrated weed management should protect or improve the biodiversity of farmland weed communities for a better ecological environment with not only increased crop yield, but also reduced use of herbicides. This study hypothesized that appropriate fertilization would benefit both crop growth and the biodiversity of farmland weed communities. To study the effects of different fertilizing patterns on the biodiversity of a farmland weed community and their adaptive mechanisms, indices of species diversity and responses of weed species and wheat were investigated in a 17-year field trial with a winter wheat-soybean rotation. This long term field trial includes six fertilizing treatments with different N, P and K application rates. The results indicated that wheat and the four prevalent weed species (Galium aparine, Vicia sativa, Veronica persica and Geranium carolinianum) showed different responses to fertilizer treatment in terms of density, plant height, shoot biomass, and nutrient accumulations. Each individual weed population exhibited its own adaptive mechanisms, such as increased internode length for growth advantages and increased light interception. The PK treatment had higher density, shoot biomass, Shannon-Wiener and Pielou Indices of weed community than N plus P fertilizer treatments. The N1/2PK treatment showed the same weed species number as the PK treatment. It also showed higher Shannon-Wiener and Pielou Indices of the weed community, although it had a lower wheat yield than the NPK treatment. The negative effects of the N1/2PK treatment on wheat yield could be balanced by the simultaneous positive effects on weed communities, which are intermediate in terms of the effects on wheat and weeds.


Introduction
Weeds are one of the major constraints to crop yields and quality [1,2]. However, as one of the primary producers within farming systems, weeds are of central importance to the arable system's food web. The weed community provides a range of resources for higher trophic groups, supports a high diversity of insect species and birds [3,4], and therefore plays an important role in the biological diversity of agroecosystems. Weeds also have other ecosystem functions, including nutrient cycling and soil preservation [5]. In addition, there are correlations between an impoverished landscape and the appearance of pests and diseases [6]. Moreover, an overreliance on herbicides has imposed a cost upon society and the environment [7]. The biodiversity of weed communities in a cropland can therefore be an important element for the reliable and sustainable provision of agroecosystem services. However, encouraging in-field biodiversity is unpopular among farmers because of the risk of decreased crop production as a result of weed competition. It is important to match crop production with conservation of biological resources to develop more sustainable systems [8].
Crop nutrient management practices may serve as an important component of more robust weed management programs [9][10][11]. Fertilization alters soil fertility, thus affecting weed density, nutrient uptake, and biomass yield, which in turn affects species composition and biodiversity [9,[12][13][14]. For example, Mahn [15] observed a general decline in the number of weed individuals and an increase in weed biomass with increasing rates of nitrogen (N) fertilizer. In another study, Digitaria ischaemum Shreb (smooth crabgrass) was found to be the dominant species under N + potassium (K) and non-fertilized treatments, Cyperus rotundus L (purple nutsedge) dominated under phosphorus (P)+K treatment, and more weed species and higher Shannon' s diversity (H9) values were detected in the balanced fertilization treatment [14]. These effects, however, indirectly depend on the light penetration caused by crop competition variation [16]. Balanced fertilization promotes the growth of crops, resulting in closed crop stands and light limitation for the weed communities growing underneath, thereby affecting weed species diversity [14,17].
Varying physiological responses of weed species to soil amendments are one of the explanations for weed community biodiversity [18]. Murphy and Lemerle [18] reported that the type and rate of fertilizers applied and the physiology of the species involved play an important role in weed population shifts. Haas and Streibig [19] also showed that high N levels will favor weed species that possess either physiological shade tolerance (esp. Stellaria media) or are able to climb into more favorable light conditions (esp. Galium aparine). Pyšek and Lepš [20] and Bengtsson et al. [21] also found that N fertilization favors shade-tolerant, climbing, and competitive weed species, but suppresses other kinds of species and results in a decrease of weed species diversity. With cessation of mineral fertilizer application during the period of conversion in a cropland, van Elsen [6] observed a gradual decline of weeds like G. aparine, which need a high nitrate level, and an increase in leguminous weeds (esp. Vicia spp.).
Documentation of the effects of a farming system on vegetation diversity is an important step toward understanding ecosystem function in agricultural landscapes [22], but only a handful of studies have examined the effect of fertilization on weed community biodiversity in agroecosystems. Biodiversity promotes ecosystem productivity, sustainability and stability in grassland [23][24]. Based on these reports, it is hypothesized that appropriate fertilization of agroecosystems would not only provide desirable crop productivity but also maintain the biodiversity of weed communities.
Long-term field experiments would ensure that proper data are accumulated and that confounding experimental effects are reduced, thereby providing better insights into the effects of prolonged fertilization over time [25]. It was also hypothesized that weed species could develop physiological adaptive mechanisms under long-term fertilization conditions. This study was conducted in a winter wheat field under continuous fertilization since 1994 with the objective of examining the effects of different N, P and K fertilization patterns on crop growth and yield. Effects on weeds were not examined at initiation, so baseline data are unavailable. Therefore, the objective of this study was to evaluate the cumulative effects of different fertilizing patterns on weed community biodiversity and their adaptive mechanisms with data from one-year sampling in this long-term experiment.

Site and Experimental Design
The long-term fertilization field experiment used in this study has been run by the Anhui Academy of Agricultural Sciences since 1994 and is located in Mengcheng county in the Anhui province of China (33u13938 N, 116u36958 E). The field studies did not involve endangered or protected species and no specific permits were required. This region has a warm temperate to sub-humid monsoon climate. The mean temperature, precipitation, evaporation and selected soil properties of the experimental site are given in Table 1. Since 1999, there has been a crop rotation of winter wheat (Triticum aestivum L.) and soybean (Glycine max [L.] Merr.).
Fertilizer treatments consisted of six combinations with different rates of N, P and K (Table 1) applied as urea, calcium superphosphate, and potassium chloride, respectively. All fertilizers were applied by soil surface broadcasting before the sowing of wheat. No fertilizer was applied during the soybean production seasons. The plot size was 20 m 2 (465 m) with three replications in a randomized complete block design. Soybean ('Zhonghuang 13') and winter wheat ('Yannong 19') were planted at a density of 249,800 seeds and 4,685,100 plants per hectare on 15 June and 20 October 2010, respectively.
Herbicides were used according to the weed spectrum and expert recommendations. All treatments had a broad application of herbicides, and manual weeding was implemented before 2008. To better examine the effects of fertilization on weed community biodiversity, only a broadcast application of herbicides was continued in each growing season, and without manual weeding after 2008. Tribenuron methyl (Anhui Research Institute of Chemical Industry, China) at 12 g ai.ha 21 for wheat and Acetochlor (Anhui Huilong Group Rmf Agrochemical Co., Ltd, China) at 1500 g ai.ha 21 for soybean were used once every growing season in all the experimental plots.

Sampling Procedures
Weed evaluations were conducted three times during wheat growing season on April 10 (177 days after sowing [d]), May 8 (205 d), and May 23 (221 d), 2010. At each sampling date (growing period), three 0.560.5 m quadrats were systematically positioned to avoid edge effects and re-sampling of the previously sampled areas in each plot. Five wheat plants and all weeds present in the quadrat were clipped, collected, sorted by species, counted, and oven dried at 70uC for 48 h before weighing. Wheat tillers were also counted, and plant heights of wheat and weed species were measured.
At 205 d and 221 d, the tallest plant of Galium aparine (a prevalent broadleaved weed species) in each quadrat was selected for the measurement of internode length. Light measurements were also made in the quadrats (using Minolta Illuminance Meter T-1H) above the wheat and on the soil surface. At 221 d, each individual weed species and wheat plant was sampled for total N, P and K analysis. The N was measured using the Kjeldahl method, and P and K were measured using an inductively coupled plasma optical emission spectrometer (ICP-OES) after the plant samples were digested [26].

Statistical Procedures
Weed density was calculated as the total plant number of a particular weed species per square meter, total weed density was calculated as the total of all weed species per square meter, and winter wheat density was calculated as the total number of tillers per square meter. The species diversity of weed communities was assessed by calculating different indices. For example, Species richness (S) was measured by the mean number of species in each treatment [27]. Shannon's diversity index (H9), Simpson index (D) and Pielou index (E) were calculated for the weed communities using the following equations: where N is the total number of individuals per plot and n is the number of individuals per species per plot. Shannon's diversity index (H9) takes not only the number of species (community richness) into account but also the relative distribution of each species in the community (community evenness). The Simpson index was applied to measure the degree of dominance of individual weed species.
Rank-abundance plots were used to display species relative abundance data. Abundance distributions provide a complete description of the community diversity and simultaneously show both components of species diversity, species number, and evenness. The relative abundance of a species indicates its degree of dominance or subordination in the weed community (i.e., the greater the relative abundance of a species in the weed community, the higher its dominance, [28]). Relative growth rate (RGR, g m 22 d 21 ), the increase of plant biomass per unit time, was used to compare crop growth with that of weed species.
Statistical analysis was done using SPSS version 16.0. Weed density and biomass data were log e (x+1) transformed before analysis to meet homogeneity of variance assumptions. For normally distributed parameters, the General Linear Model Univariate was used, and means were compared based on Tukey's multiple comparison tests P,0.05. Otherwise, nonparametric tests were used, and median values were presented.

Structural Changes and Biological Diversity at the Community Level
The results indicated that total weed density was influenced by fertilization (P,0.0001), but was not affected by growing periods (P = 0.59) (Fig. 1). A general decline in the number of weed individuals was observed in the different fertilization treatments except in the NK treatment where the number increased slightly. The PK treatment showed the highest total mean weed density (711 plants m 22 ) followed by the control and NK treatments. Weed densities in these three treatments were significantly higher than the 201 plants m 22 in the treatments with N and P (72% less than that in the PK treatment).
Fertilization increased weed biomass in the winter wheat field (Fig. 2). The PK treatment obtained the highest (P,0.05) total weed biomass of 108.03 and 117.94 g m 22 at 205 and 221 d, respectively. Treatments with N and P showed a lower weed shoot biomass.
The species diversity of the weed community was modified by fertilizing patterns. Species number and evenness were affected differently by the fertilization treatments, which were shown in (a) the rank-abundance plots and (b) the vertical bar chart (Fig. 3). Among all the recorded species (13 spp.), the species numbers in the PK (11 spp.) and N1/2PK (11 spp.) treatments reached the highest, followed by the control (9 spp.) and NK (9 spp.) treatments, while the lowest value (7 spp.) was observed in the NPK treatment.
The evenness of the weed community also varied with the fertilization treatments (Fig. 3a, b). Fertilizing patterns greatly shaped the equitability in the partitioning of total biomass among species in the community as shown by the slope of the rankabundance plots. A higher Shannon-Wiener Index (Fig. 3b) and weed community evenness (Pielou Index, Fig. 3a, b) was observed in the PK treatment, followed by the NK and control treatments. However, in these treatments, the Simpson Index was lower compared with that in the NP and NPK treatments (Fig. 3b). Compared with the NPK treatment, the half P rate treatment (N1/2PK) showed a higher equitability of total biomass partitioning ( Fig. 3a) with higher Pielou Index and Shannon-Wiener Index values (Fig. 3b).

Changes at the Population Level
Densities of Weed Population. Galium aparine L. (tender catchweed bedstraw), Vicia sativa L. (common vetch), Veronica persica Poir (iran speedwell), and Geranium carolinianum L. (carolina geranium) were the most dominant species in the weed community. These four dominant species comprised .90% of the total weed density. Average densities of the four prevalent broadleaved weeds during the three growing periods in the different fertilization treatments are shown in Fig. 4. Densities of these four weed species varied significantly among the treatments. Galium aparine densities were greatly affected by fertilization (P,0.0001) and growing period (P = 0.02). No significant difference was detected in the interactions of treatments with growing period (P = 0.43). The densities of G. aparine (138 plants m 22 ) were higher in the treatments with N and P fertilizers, while no G. aparine was found in the control treatment.
Vicia sativa density was significantly affected by fertilization (P,0.0001) and growing period (P = 0.008). However, in contrast to G. aparine, Vicia sativa grew better in treatments without N and/ or P with an average density of ,200 plants m 22 , while it had only 10-30 plants m 22 in the treatments with N and P. Geranium carolinianum density was also affected by fertilization (P,0.0001) and growing period (P,0.001), and the interactions of fertilization with growing period (P,0.001). Only three plants of G. carolinianum were found in the NP treatment at 205 d and no plants were found at 221 d. During the early growing period, G. carolinianum densities showed no significant differences among the treatments, with an average value of 72 plants m 22 . However, its densities responded differently to fertilization treatments in the latter growing period, with the same tendency as V. sativa in each treatment. Veronica persica grew well only in the PK treatment, and few or no plants were detected in the other treatments.
Plant Height of Weed Population. Among the four prevalent broadleaved weed species, V. persica grew underneath the crop canopy, and its maximum plant height was ,50 cm ( Fig. 5D). On the contrary, other species were erect for capturing adequate light. Fertilization significantly affected plant height of G. carolinianum (Fig. 5A), V. sativa (Fig. 5B) and G. aparine (Fig. 5C). Geranium carolinianum was taller in the PK and NK treatments, even taller than in the NPK treatment at 205 d. Plant height of V. sativa in the PK treatment was at a maximum, reaching 73.8 cm at maturity stage, followed by the NK treatment where it reached 64.1 cm, while lower values were observed in treatments with N and P. Conversely, G. aparine showed higher plant heights of 83.8 cm, 85.3 cm and 73.3 cm in the NP, NPK and N1/2PK treatments, respectively. The study confirmed that fertilization had little effect on the number of G. aparine shoot internodes. However, fertilization significantly modified the internodal length of G. aparine, especially from the 4 th to the 12 th internode ( Table 2).
Shoot Biomass of the Weed Population. Shoot biomass of the four weed species varied with different fertilization treatments (Table 3). Like the weed density, G. aparine shoot biomass in the NP treatment was highest at each growing period, followed by the  NPK and N1/2PK treatments. The relative growth rates of G. aparine in the treatments without N and/or P were lower, especially in the NK treatment (0.11 g m 22 d 21 ), than that in the NP treatment (0.83 g m 22 d 21 ). The highest proportion of G. aparine shoot biomass to total weed biomass was obtained in the NPK treatment (91%), followed by the NP and N1/2PK treatments at  89% and 77%, respectively, while treatments without N and/or P had the lowest value (8%). These results indicated that G. aparine had an absolute growth advantage under higher soil fertility conditions, and its growth could be significantly promoted by the application of N and P fertilizers.
Biomass accumulation of V. sativa followed an opposite trend to that of G. aparine for the different fertilization treatments. The average shoot biomass of V. sativa in the control, NK and PK treatments were 3.32, 14.81, and 25.17 g m 22 , respectively in contrast to values under the NP, NPK and N1/2PK treatments where the values were 0.63, 0.69, and 1.33 g m 22 at 177, 205 and 221 d, respectively. The relative growth rate showed that V. sativa also grew rapidly in the control, PK and NK treatments (0.50 g m 22 d 21 ). The highest proportion of V. sativa shoot biomass to the total weed biomass was observed in the control treatment with 68% and 41% at 205 and 221 d, respectively, followed by the NK and PK treatments, while treatments incorporating N and/or P fertilizers showed the lowest values (4%).
The average shoot biomass of G. carolinianum at 177 d reached 5.05 g m 22 , which was much higher than the corresponding values for G. aparine, V. sativa and V. persica (1.72, 1.97 and 0.55 g m 22 , respectively). The G. carolinianum shoot biomass under the NPK treatment (9.11 g m 22 ) was higher than that under control. However, in the middle and latter growing periods, G. carolinianum shoot biomass performed similarly to V. sativa with the different fertilization treatments. Specifically, biomass accumulation of V. persica was only higher in the PK treatment. Its biomass at 205 d accounted for 81% of the total biomass at 221 d, and it also accounted for .50% of the total shoot biomass of weeds present in the PK treatment.
Nutrient Uptake by Weeds. Nutrient uptake by the four weed species differed among treatments (Table 3). In general, nutrient uptake was higher in G. aparine and V. sativa than in G. carolinianum, suggesting that the two weed species were stronger nutrient competitors to wheat. Nutrient uptake of weed species in different treatments followed a pattern similar to shoot biomass in the corresponding treatments. Average N, P, and K uptakes by G. aparine in the NP, NPK, and N1/2PK treatments were 0.353, 0.090, and 0.447 g m 22 , respectively, which were 435, 375, and 425% higher than those under the PK treatment, and 45, 158, and 62% higher than those under the NK treatment. The N, P, and K uptakes under the control, PK, and NK treatments were 0.426, 0.076, and 0.361 g m 22 , respectively for V. sativa and 0.163, 0.052, and 0.187 g m 22 , respectively for G. carolinianum. For V. persica, the highest nutrient uptake was found under the PK treatment, and its P and K uptakes were 127 and 241% higher, respectively than those of V. sativa.

Winter Wheat Responses to Fertilization
Winter Wheat Density. Significant differences in wheat density were found among the fertilization treatments (P,0.001). Wheat density increased in the treatments that included N and P. Wheat density in the NP treatment was significantly (P,0.01) higher than that in the control, PK and NK treatments. However, there were no significant differences in winter wheat density within the NP, NPK and N1/2PK treatments. The average density was 524 tillers per m 2 in the latter three treatments, which was 82% higher than that in the control treatment (data not shown). Plant Height of Winter Wheat. Fertilization significantly increased wheat height at all the growing periods and the final heights ranged from 60.0 to 79.4 cm at maturity stage (Fig. 6). Wheat heights of the NPK (78.2 cm) and N1/2PK treatments (79.4 cm) were 31 and 33% higher, respectively than the control treatment.

Winter Wheat Biomass, Nutrient Uptakes and Grain
Yield. Fertilization significantly increased the shoot biomass accumulation of winter wheat (Table 4). Higher wheat shoot biomasses were observed in the treatments with N and P. The highest wheat shoot biomass was found in the N1/2PK treatment, and significant differences were also found in two of the three growing periods within the treatments containing N and P.
Wheat nutrient uptake is usually positively correlated with its shoot biomass accumulation. However, the wheat P uptake had no significant difference between the N and P treatments in this study ( Table 4). The highest N and K uptakes were in the N1/2PK treatment. They were 18 and 78% higher than those in the NP treatment and 35 and 31%, higher respectively than those in the NPK treatment. These results might be attributed to the lower shoot biomass and nutrient accumulations of the weed community in the N1/2PK treatment (Fig. 2, Table 4). It also suggested that reducing P application in the NPK treatment decreased the productivity of weeds, such as G. aparine (Table 3), and improved nutrient uptake and biomass accumulation of wheat.
The NPK treatment resulted in the highest wheat yield, followed by the N1/2PK treatment (Fig. 7). In contrast, the control treatment obtained the highest light penetration (55%). The average light transmittance rate in the treatments with N and P was 5% at 205 d, which was 92% lower than the light transmittance value of the control treatment (Fig. 7).

Discussion
Weed community structure changed sharply among the different fertilization treatments and the results concurred with previous research conducted under field conditions by Mahn [15]. The PK treatment showed the highest weed density and shoot biomass, followed by the NK and control treatments, while the treatments with N and P resulted in lower values. Yin et al. [14] and Nie et al. [17] also reported similar results. The structural variability of the weed community across treatments could be directly explained by the different soil fertility conditions. For example, soil available N concentration plays an important role in terminating dormancy by removing germination constraints and promoting the emergence of seedlings [29]. Soil available P was the primary nutrient regulating the species composition and floristic construction of the weed communities, followed by N and K [30]. Moreover, solar radiation reaching weeds was indirectly modulated by different fertilization treatments (Fig. 7). Wheat growth was promoted in the treatments with N and P, especially in the balanced fertilization treatments (Table 4, Fig. 7), increasing its ability to intercept solar radiation (Fig. 7), and thereby affecting the growth of weed species. The structural variability of the weed community may also be significantly affected by rainfall and other environmental factors, such as mean annual temperature. Ochoa-Hueso and Manrique [31] showed that apart from nutrients, plant germination and growth are also influenced by water availability. Pal et al. [32] and Š ilc et al. [33] also reported that environmental factors, such as mean annual precipitation and temperature showed significantly positive and negative correlation coefficients with weed species composition, respectively. The four prevalent weed populations and wheat differed substantially in their responses to different fertilization treatments. Galium aparine grew well in high fertility soil conditions, and was taller than wheat. Vicia sativa and G. carolinianum grew well under the control, NK, and PK treatments, and they were also taller than wheat in these treatments. Veronica persica grew well only in the PK treatment. These results show that the different weed species had formed their own adaptive mechanisms. Morphological plasticity enables plants to cope with a wide variety of ecological conditions, and nutrient availability has been widely shown to induce important anatomical responses [34]. For instance, G. aparine grew well under higher soil fertility conditions through climbing into more favorable light conditions by internode elongation (Table 2),  and responded to soil nutrient supply in the order of P.N.K. In contrast, V. sativa and G. carolinianum increased their heights for receiving more light under lower soil fertility and limited nutrient conditions. Li et al. [35] reported similar results for G. carolinianum, but Yin et al. [13] reported different results for G. aparine versus V. sativa and V. persica. These differences in the effects of fertilization on cropland weeds might be due to herbicide usage. A broadcast application of herbicides was used by Li et al. [35], while no herbicide was used by Yin et al. [13]. Furthermore, environmental conditions, such as plot location, temperature, crop type, precipitation, soil texture, neighboring habitat, and soil pH, which vary in time and space, may also cause variation in weed species responses [36][37][38].
The PK treatment showed the highest Shannon-Wiener and Pielou Indices, followed by the control and NK treatments. In contrast, these three treatments had lower Simpson Indices than the treatments with N and P. Solar penetration on the ground surface was also one of the reasons for the variation in weed community biodiversity, which is similar to the findings of other studies [16,39,28]. Biomass is often used as an indicator of the amount of resources captured by a crop. Higher biomass and grain yield of wheat in the NP, NPK and N1/2PK treatments indicated that wheat intercepted a larger proportion of light radiation (Fig. 7), thereby reducing the amount reaching the weeds. Wheat yield showed a significant negative correlation with light penetration (Pearson correlations .0.880, P,0.0001). The species richness, Shannon-Wiener and Pielou Indices and the Simpson Index of weed communities were proved to have significantly negative and positive linear function relationships with wheat yield [30]. Shading leads to thinning mortality in individuals of subordinate weed species, and consequently reduces evenness and biodiversity among weed species [39,40]. Only those weed species which possess either physiological shade tolerance or are able to climb into more favorable light conditions could grow well [18]. In contrast to other species, G. aparine had higher plant density, shoot biomass and nutrient accumulations in the NP, NPK and N1/2PK treatments with lower light penetration. Mahn [15] also reported that the mean biomass values for G. aparine in higher N fertilization treatments within a spring barley field were almost two times more than the values under the no N treatment. Under high soil fertility conditions in this study, G. aparine could enhance its plant height through internode elongation, especially from the 4 th to the 12 th internodes for capturing more sun light. Its plants were 14% and 9% taller than wheat plants in the NP and NPK treatments, respectively (Fig. 5). Therefore, weed vegetation under higher soil fertility conditions was predominantly composed of a few species, such as G. aparine, resulting in a higher Simpson Index and lower Pielou and Shannon-Wiener Indices.
More weed biomass and nutrient accumulations were detected in the PK treatment. Weed biomass under the NP, NPK and N1/ 2PK treatments was 26.4, 26.2, and 26.1%, respectively, of that under the PK treatment. Similarly, the N, P and K accumulations, as a proportion of that in the PK treatment, were 25, 14, and 13%, respectively in the NP treatment, 36, 19, and 24%, respectively in the NPK treatment and 25, 12, and 18%, respectively in the N1/ 2PK treatment (Table 4). These results indicated that fewer resources were captured by the weed community in the treatments with N and P, especially in the N1/2PK treatment. If a single weed species has a shoot biomass to total weed biomass ratio . 10%, it is regarded as a dominant weed species. More than two dominant weed species were found under the control, PK and NK treatments, with four species present in the PK treatment. Only one dominant weed species (G. aparine) was present in the balanced fertilization treatment. Thus, the results also indicated that balanced fertilization could not only maintain a stable crop yield, but also greatly reduce weed productivity.
Lower solar radiation under high soil fertility conditions resulted in fewer weed species compared with the control and PK treatments. However, the number of weed species was only 73% and 64% in the NP and NPK treatments, respectively and their number in the N1/2PK treatment was the same as in the PK treatment ( Table 4). These results demonstrated that G. aparine is sensitive to soil P, and reducing the P application rate indirectly improved the growth of other weed species. Blackshaw and Brandt [41] also reported that the competitiveness of the high Presponsive species was progressively improved as the P dose increased. Likewise, with a lower P dose in the N1/2PK treatment, its competitive ability would decrease and more opportunities would be available for other species. The N1/2PK treatment could be considered as a balanced fertilization treatment in adjusting the interactions of wheat and weeds in agroecosystems. Although wheat yield in the N1/2PK treatment was 31% less than that in the NPK treatment, its weed species number, Shannon-Wiener and Pielou Indices were significantly higher than those in the NPK treatment. It is noteworthy that a certain amount of G. aparine would have a high biodiversity value for invertebrates and is important for seed-eating birds [22]. Thus, applying half the P rate in the NPK treatment could be considered as a better strategy to provide desirable farming benefits and maintain the biodiversity of weed communities. The negative effects on wheat yield could be balanced by the simultaneous positive effects of reduced herbicide usage. Appropriate combinations and rates of fertilization would not only be helpful for keeping desirable crop productivity, but also maintaining the biodiversity of weed communities in an agroecosystem.

Author Contributions
Conceived and designed the experiments: FC DZW. Performed the experiments: CPC KYW RHL YT JFP LLT. Analyzed the data: LLT. Contributed reagents/materials/analysis tools: JX LLT. Wrote the paper: LLT.