1. Introduction
Globally, weeds interfere with agricultural production activities and have a very serious impact on agriculture and animal husbandry [
1,
2]. Weeds and crops form a complex system in the ecosystem, competing directly for resources, including water, nutrients, light and space, which can affect the growth, yield and quality of crops, ultimately leading to substantial economic losses[
3]. Weeds caused 45% of yield loss; however, this percentage may rise to 94-96% in rice, 50% in pulses, 72% in sugarcane and nearly 90% in almost all vegetables with an increase in weed interference[
4]. In traditional agriculture, farmers had a proactive approach to weed removal, clearing them as soon as they were noticed[
5]. Nonetheless, hand weeding is time-consuming, labor-intensive and tedious, leading to higher production costs [
6]. Although chemical weed control using herbicides has been the most popular and effective method [
7,
8], it affects the quality and safety of crops and pollutes the environment. In addition, a single herbicide application cannot control the entire weeds over the growing season[
9]. Moreover, the evolution of herbicide resistance has also become a major problem, further complicating the process of weed control [
10]. All these factors can significantly impact the sustainable development of farmland ecosystems. Therefore, identifying a safe and reliable weed control strategy that can not only increase yield and production net income but also improve crop quality is of great significance.
Previous studies have shown that some weed control strategies can improve the yield and quality of cultivated crops. Weed diversity mitigates crop yield loss. For instance, yield loss decreased from 60% to 30% in wheat paddocks with an increase in the weed species richness from 7 to 20 [
11,
12]. As one of the common weeds of grain fields, corn cockle (
Agrostemma githago L.) promotes crop growth and improves the yield and quality in the wheat field by secreting gibberellin, agrosfemin and allantoin into the soil [
13,
14].
Proactive management of weed diversity may increase crop productivity[
15]. Most studies have shown that weeds impact crop yield, the yield is sub-correlated with the intensity of weed disturbance and constantly clearing all weeds during the entire growth period is the best approach to increase output. A previous study found that corn fields without weed interference had the highest grain yield[
16]. Continuous weeding significantly reduced weed density and weed dry matter and increased weed control efficacy, ultimately improving plant growth, yield attributes and the overall yield of groundnut [
17]. Similar results were reported in wheat[
7], rice [
18] and other crops[
19]. Some studies have also suggested that weeds not only interfere with crop yield but also affect quality. The yield of the paddy field without weed interference was about 50% higher than that of the weed field and rice processing, and appearance quality and seed nitrogen accumulation were also the best[
20]. Other studies have shown that although weed interference affects crop growth and yield, they have no significant effect on quality. Notably, an increasing duration of weed interference negatively affected crop height, head diameter and 1000-kernel weight but not seed oil content[
21]. Three times of artificial weeding in the early stage of the
Radix bupleuri field significantly increased the root weight but had little effect on the total amount of saponin A and saponin D[
22]. It can be seen that weed control has different effects on the yield and quality of cultivated crops. Photosynthesis [
23,
24], anti-stress enzyme activity [
25,
26] and nitrogen metabolism during growing season significantly affect crop yield and quality. However, the impact of weed interference frequency on these plant indicators has not been fully tested.
Licorice (
Glycyrrhiza uralensis Fisch) is a perennial herb of the genus
Glycyrrhiza in the family Leguminosae. It is one of the most popular traditional Chinese herbal medicines and has been used for over 2000 years in China. Its dried roots and rhizomes are widely used in Eastern and Western countries as medicine. This species contains numerous active ingredients [
27,
28], including liquiritin (LQ) and glycyrrhizic acid (GA), which are the main medicinal components specified in the Chinese Pharmacopoeia. The artificial cultivation of wild
Licorice began in the 1980s in China and has become one of the leading industries in recent years. Weed control is a main technical obstacle restricting the large-scale cultivation of
Licorice.
In the past 40 years,
Licorice has been comprehensively studied in terms of plant systematic classification[
29], chemical composition[
30], pharmacological efficacy [
31,
32], separation and extraction methods [
33,
34,
35], metabonomics [
36,
37], clinical application [
38,
39,
40], etc. Most studies on artificial
Licorice cultivation have largely focused on seed treatment [
41,
42], the best sowing and harvesting date [
43,
44], water stress [
37,
45], fertilization [
46] and pest control[
47], with only a handful of studies on weed control and herbicides. Herbicides pose significant challenges to the safe medical use of licorice.
We hypothesized that (i) different weed interference frequencies might affect the growth, photosynthetic physiology and enzyme activity of Licorice, and then change its yield and quality and (ii) the effects of different weed interference frequencies on the yield and quality might have a critical value, which is conducive to improving the economic benefits of its ecological planting. The present study aimed to (i) identify the differential efficacy of weed interference frequencies on yield and quality and (ii) determine the optimal interference frequency and the cost of ecological planting.
2. Materials and Methods
2.1. Study Site
The field experiment was conducted in 2021 and 2022 at Tianjizhang village of Licorice cultivation base, Yanchi County, Ningxia, in northwest China (Long. 107º16′48"E, Lat. 37°48′0"N, ca. 1427 m above sea level). The region is characterized by a typical continental monsoon climate with a mean annual temperature of 7.7°C, total annual precipitation of 290 mm and annual evaporation of 2132 mm. The soil type is mainly aeolian sandy soil, sierozem, which is barren and vulnerable to wind erosion. The region is a typical desert steppe area, dominated by xerophytes and mesoxerophytes.
2.2. Experimental Set-Up
One-year-old transplanted seedlings of Licorice were planted in early May at a depth of approximately 15 cm with a 30 × 12 cm spacing row and harvested in mid-October in each year (2021-2022). The planting density was about 27 plants m-1. Each plot had an area of 3 × 5 m.
The weed control test was started after the emergence of Licorice. The experiment consisted of seven treatments and was carried out as a randomized complete block (RCB) design with three replications. The treatments were (1) no weeding, marked as WF0; (2)-(7) artificial weeding by hoe once every 1, 2, 4, 6, 8, 10 weeks after emergence, marked as WF1, WF2, WF4, WF6, WF8 and WF10, respectively. All plots, whether weeding or not weeding, received the same agronomic management practice.
2.3. Weeds Investigation
Three quadrats (1 × 1 m) were randomly selected in the three plots without weeding. Weeds were removed from the ground and kept in ovens at 80℃ to a constant weight. Then, aboveground biomass accumulation (g/m2) was calculated. Weed species and density were also recorded.
2.4. Sampling and Measurements
2.4.1. Growth index
In mid-September of each year, 30 healthy plants with uniform growth were chosen for plant height (vertical height through the center of the crown) and crown size (length and width through the center of the crown) measurement using steel tape (accuracy of 0.1 cm) in different treatment plots. The ground diameter was measured using a vernier caliper (accuracy of 0.01 mm).
2.4.2. Photosynthetic parameters
Healthy leaves were selected at 2/3 of the plant from the ground on a sunny day at the end of August of each year for the measurement of the net photosynthetic rate (Pn, μmol/m2·s), transpiration rate (Tr, mmol/m2·s), stomatal conductance (G, mmol/m2·s) and intercellular carbon dioxide (CO2) concentration (Ci, μmol/mol) using a LI-COR 6800 portable photosynthesis system (LI-COR Biosciences, Lincoln, NE, USA). The reference (CO2) was set to 400 μmol/mol, light intensity was 1000 μmol/m2·s, temperature was 24°C and water vapor pressure difference (VPD) was 1.0 kPa. Three leaves were selected for each treatment as three replicates. The measurement time was between 9:00 and 11:00 A.M., and leaves were tiled to cover the whole leaf chamber (6 cm2). The chlorophyll soil plant analysis development (SPAD) value was measured using a portable chlorophyll SPAD meter (SPAD-502 Plus, Konica Minolta, Japan).
2.4.3. Antioxidant and nitrogen metabolism enzymes
Healthy leaves were collected at the end of August every year and stored in a dry ice sampling box. The leaves were taken back to the laboratory and placed in a refrigerator at -80℃ for determination. Antioxidant enzymes, including superoxide dismutase (SOD), peroxidase (POD) and catalase (CAT), and nitrogen metabolism enzymes, including glutamine synthetase (Gs), nitrate reductase (NR), nitrite reductase (NiR) and glutamate synthase (NG), were determined by spectrophotometry (BioPhotometer, Eppendorf, Germany) according to the instructions of the corresponding kits (provided by Suzhou Michy Biomedical Technology Co., Ltd., China).
2.4.4. Yield and effective components
All Licorice roots in each plot (3 × 5 m) were dug out in mid-October each year and dried naturally. Its dry weight was measured and the yield per hectare was calculated. The content of LQ and GA were determined by high-performance liquid chromatography (HPLC) according to the 2020 edition of Chinese Pharmacopoeia.
2.5. Data Analysis
Analysis of variance (ANOVA) was estimated using Statistical Analysis System software (SAS 8.1) and means were compared following a protected least significant difference (LSD) procedure at a 5% level of probability using the Duncan Multiple Range Test (DMRT). Data from the two seasons were examined independently. The box line chart, correlation analysis graph and principal component analysis (PCA) graph were performed using Metware Cloud — a free online platform for data analysis (
https://cloud.metware.cn ).
3. Results
3.1. Weeds Investigation
There were 9 species of weeds in the non-weeding plot, belonging to 4 families and 7 genera (
Table 1). Weeds with a higher density and aboveground biomass were
Setaria viridis and
Chenopodium album, with a density of 27 m
-2 and 7 m
-2 and biomass of 62.12 g·m
-2 and 236.71 g·m
-2, respectively. The weed densities in the plot without weeding was 46 m
-2, about 465000 ha
-1.
3.2. Growth of Licorice
Different weed interference frequencies had different effects on growth in both growing seasons (
Table 2). Generally, the growth index value showed a downward trend with a decrease in the weeding frequency. The plant height, ground diameter and crown width of
Licorice treated with WF1 were the highest, and the values in 2021 were 51.17 cm, 5.51 mm and 963.23 cm
2 respectively, which were 31.54, 31.19% and 94.55% higher than those of WF0 treatment (lowest values), and the difference was significant (P < 0.05). The values in 2022 were 42.67 cm, 5.13 mm and 747.52 cm
2 respectively, which were also significantly increased by 81.96, 95.05 and 137.39% compared with WF0 treatment (lowest values).
Meanwhile, no significant difference in plant height index was found between WF1 and WF2 treatments in 2021 (P > 0.05). Similarly, no significant difference was detected between WF1 and WF2 or WF4 treatments in 2022 (P > 0.05); however, the difference between WF1 and other treatments was significant in both growing seasons (P < 0.05). Overall, no significant differences in plant height, ground diameter and crown width were observed between WF2 and WF4, WF6 and WF8, WF10 and WF0 treatments (P > 0.05).
3.3. Photosynthetic Parameters
The effects of different weed interference frequencies on photosynthetic parameters were also different in both growing seasons (
Table 3). The effects of different treatments on the
Pn were WF1 > WF4 > WF2 > WF6 > WF8 > WF10 > WF0 in 2021 and WF1 > WF4 > WF2 > WF8 > WF6 > WF10 > WF0 in 2022. The effects on the
Tr were WF1 > WF4 > WF2 > WF8 > WF6 > WF10 > WF0 in 2021 and WF1 > WF4 > WF2 > WF6 > WF8 > WF10 > WF0 in 2022.
Pn, Tr, G and SPAD values in the WF1 treatment within 2 years were the highest. Compared with the WF0 treatment (lowest values), these parameters increased by 0.85, 0.47, 1.96 and 0.33 times and 0.96, 0.73, 1.50 and 0.39 times in 2021 and 2022 respectively, and the differences were significant (P < 0.05). However, no significant differences in these parameters were found between WF2 and WF4 treatments and among WF6, WF8 and WF10 treatments (P > 0.05).
3.4. Antioxidant Enzyme Activities
Generally, the values of antioxidant enzyme activity in different treatments showed an increasing and then a decreasing trend in both growing seasons (
Table 4). Notably, the WF1 treatment exhibited the lowest SOD, CAT and POD values.
In 2021, the SOD activity in WF10 and WF0 treatments were significantly higher than that in other treatments (P < 0.05), which was 4.01 and 3.74 times higher than that in WF1 treatment, respectively. The CAT activity in the WF8 treatment and POD activity in the WF10 treatment were significantly higher than those in other treatments (P < 0.05). In 2022, the SOD activity in WF8 treatment was the highest, and the change in CAT and POD activities was consistent with that in 2021.
3.5. Nitrogen Metabolism Enzyme Activities
Similarly, the activities of nitrogen metabolism enzymes in different treatments also exhibited an increasing and then a decreasing trend in both growing seasons (
Table 5). The value of higher weeding frequencies was significantly higher than that of lower weeding frequencies. NR, NiR, NG and Gs activities in WF1, WF2, WF4 and WF6 treatments were higher, while those in WF8, WF10 and WF0 treatments were lower. NR, NiR and NG activities in the WF4 treatment were the highest. Compared with the WF0 treatment, these parameters increased by 1.76, 0.67, and 2.19 times and 6.93, 0.56, and 2.33 times in 2021 and 2022 respectively, and the differences were significant (P < 0.05). Gs activities were the highest in WF6 and WF4 treatment in 2021 and 2022 respectively.
3.6. Yield and Effective Components
The yield showed a downward trend with a decrease in the weeding frequency in both growing seasons (
Figure 1). The highest yield was obtained in the WF1 treatment, and the dry weight reached 2536.46 and 2122.20 kg/ha in 2021 and 2022 respectively, which was 71.41% and 55.65% higher than that without weeding respectively, and the difference was significant (P < 0.05). However, no significant differences were observed among WF1, WF2 and WF3 treatments and between WF8 and WF10 treatments (P > 0.05).
The contents of LQ and GA in different treatments were much higher than those specified in the 2020 edition of Chinese Pharmacopoeia (LQ: 0.50%, GA: 2.00%) (
Figure 2 and
Figure 3). The highest LQ and GA were found in the WF4 treatment in both growing seasons, and the lowest value was found in the WF0 treatment. In 2021, LQ content increased by 35.67% and 53.24% in the WF4 treatment compared with the WF1 and WF0 treatments. The GA increased by 11.74% and 36.57% in the WF4 treatment compared with the WF1 and WF0 treatments. In 2022, LQ and GA increased by 19.23% and 21.30% in the WF4 treatment compared with the WF1 treatment respectively and by 24.71% and 35.48% compared with the WF0 treatment respectively.
3.7. Economic Benefit
The change trend in the gross income and dry weight of Licorice under different treatments was similar, while the net income changed differently in both growing seasons. WF4 treatment exhibited the largest increases in the net income among all treatments in both growing seasons, with respective increases up to 71.39 % and 78.81% compared with the WF0 treatment. The net income in the WF6 treatment was the second and that in the WF0 treatment was the lowest. The WF4 treatment did not differ significantly from the WF6 treatment (P > 0.05) but significantly differed from other treatments (P < 0.05).
Table 6.
Effect of weeds interference frequencies on Economic Benefit .
Table 6.
Effect of weeds interference frequencies on Economic Benefit .
Year |
Weeds interference frequency |
Weeding costs (CNY/ha) |
Gross income (CNY/ha) |
Net income (CNY/ha) |
2021 |
WF1 |
12000 |
40583.31±1276.38a |
16583.31±1276.38bc |
WF2 |
6000 |
33973.09±1986.02b |
15973.09±1986.03cd |
WF4 |
3000 |
35011.40±353.84 b |
20011.40±353.84a |
WF6 |
1950 |
32720.22±1834.66b |
18770.22±1834.66ab |
WF8 |
1500 |
28103.03±946.79c |
14603.03±946.79cd |
WF10 |
1200 |
27122.79±1186.75c |
13922.79±1186.75de |
WF0 |
0 |
23675.86±1413.25d |
11675.85±1413.25e |
2022 |
WF1 |
12000 |
33955.20±1312.96 a |
9955.20±1312.96d |
WF2 |
6000 |
33225.60±1743.35a |
15225.60±1743.35b |
WF4 |
3000 |
32548.80±1177.42 a |
17548.80±1177.42a |
WF6 |
1950 |
29336.00±1452.55 b |
15386.00±1452.55ab |
WF8 |
1500 |
26294.40±1086.57c |
12794.40±1086.57c |
WF10 |
1200 |
26325.60±1358.29c |
13125.60±1358.29bc |
WF0 |
0 |
21814.40±815.04d |
9814.40±815.04d |
3.8. Combined Correlation Analysis of Growth and Physiological Indexes with Yield and Effective Components
The correlation between growth and physiological indexes (photosynthetic indicators, antioxidant enzyme activity, nitrogen metabolism enzyme activity) and yield and effective components of Licorice was analyzed, and the Mantel test results were as follows.
In 2021, the dry weight of
G.uralensis was significantly positively correlated with GD,
Pn, SPAD, CS,
Tr and G (Sorted by the size of the correlation coefficient) and negatively correlated with SOD, POD and CAT. LQ was positively correlated with NG, NiR, Gs and NR. GA was significantly positively correlated with Gs and NR (
Figure 4 and
Figure 5).
In 2022, the dry weight of
G.uralensis was significantly positively correlated with Pn, Tr, SPAD, CS, PH and GD and negatively correlated with POD, CAT and SOD. Meanwhile, LI was significantly positively correlated with NR, NG and Gs. Additionally, GA was significantly positively correlated with Gs, NR, NiR and NG (
Figure 6 and
Figure 7).
3.9. PCA of Weed Interference Frequencies
To assess the difference between the seven weed interference frequencies in the field, PCA was performed on 18 indicators, including growth, photosynthetic indicators, antioxidant enzyme activity, nitrogen metabolism enzyme activity, yield and effective components of
G.uralensis (
Figure 8).
PCA showed that the total variance of the model described by PC1 and PC2 axes was 73.29% in 2021 and 67.44% in 2022, which suggested that there were differences among the seven treatments.