Next Article in Journal
Traditional Sweet Chestnut and Hybrid Varieties: Chemical Composition, Morphometric and Qualitative Nut Characteristics
Previous Article in Journal
Can Soil pH Correction Reduce the Animal Supplementation Needs in the Critical Autumn Period in Mediterranean Montado Ecosystem?
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Modeling the Effects of Nitrogen Fertilizer and Multiple Weed Interference on Soybean Yield

1
Department of Agriculture, Forestry, and Bioresources, Research Institute of Agriculture and Life Sciences, College of Agriculture and Life Sciences, Seoul National University, Seoul 08826, Korea
2
Institute of Plasma Technology, Korea Institute of Fusion Energy, Gunsan 54004, Korea
3
Mushroom Research Division, National Institute of Horticultural and Herbal Science, Rural Development Administration, Eumseong 27709, Korea
4
Crop Protection Division, National Institute of Agricultural Sciences, Rural Development Administration, Wanju 55365, Korea
*
Author to whom correspondence should be addressed.
Agronomy 2021, 11(3), 515; https://doi.org/10.3390/agronomy11030515
Submission received: 15 January 2021 / Revised: 5 March 2021 / Accepted: 6 March 2021 / Published: 10 March 2021
(This article belongs to the Section Weed Science and Weed Management)

Abstract

:
Understanding the effects of nitrogen (N) fertilizer on soybean-weed competition is essential for establishing a practical tool for N application and weed management. A two-year field experiment was conducted in a soybean field located in Bogatyrka (43.82° N, 131.6° E), Primorsky krai, Russia, to investigate the effects of N fertilizer and multiple-weed interference on soybean (Glycine max) yield and to model these effects. Soybean yield loss caused by the interference of multiple weeds including common ragweed (Ambrosia artemisiifolia), barnyard grass (Echinochloa crus-galli), and American slough grass (Beckmannia syzigachne) at different levels of N fertilizer was accurately described by a combined model incorporating inverse quadratic and exponential models into the rectangular hyperbolic model for two parameters Y0 and β, respectively. The combined model used in our study indicated that the application of N up to 36 kg N ha−1 can increase weed-free soybean yield by 2.2 Mg ha−1 but soybean yield under multiple-weed interference can sharply decrease with increasing total density equivalent, particularly at 36 kg N ha−1. These results, including the combined model, thus can support decision making for weed management under different N uses in soybean cultivation.

1. Introduction

Nitrogen (N) fertilizer plays an important role in soybean production although soybean requires much less N fertilizer than other crops. Soybean is less dependent on N fertilizer because symbiotic rhizobia provide soybean with ammonium (NH4) by biologically fixing N2 gas in the air. Nonetheless, soybean cultivation requires a significant amount of N fertilizer, particularly at its early growth stage when rhizobia have not established a symbiotic relationship with soybean. On average, 55% of soybean N demand is reported as being met by biological N2 fixation [1]. Therefore, a full basal dressing of approximately 16 to 30 kg N ha−1 is artificially applied for soybean cultivation to supply insufficient N by biological N2 fixation in order to promote early vegetative growth and thus to attain maximum soybean yield in practice [2,3,4].
In the region of Primorsky krai, Russia, soybean is intensively cultivated during four months from late May until late September due to the short frost-free period, resulting in a short growing season. Due to this, N application is essential to promote early soybean growth and thus to secure soybean yield in the region. However, soybean growth and yield have been severely interfered with regional weed species, although N application increases soybean yield. In the region, soybean yield is as low as 64.5% of the average yield (1.69 Mg ha−1) in neighboring regions, including Heilongjiang province of China [5,6,7]. Ambrosia artemisiifolia (common ragweed), Echinochloa crus-galli (barnyard grass), and Beckmannia syzigachne (American slough grass) are dominant in the soybean fields of Primorsky krai [8]. Their competition effects on soybean yield under N fertilization have not been fully investigated. Ambrosia artemisiifolia and E. crus-galli are known as nitrophilous weed species which increase their biomass along with increasing N level [9,10]. When weeds are not adequately managed, N application can increase weed growth; rather, crop yield can be offset by weeds outcompeting for N. Many studies have reported that crop yield can be reduced more by nitrophilous weed species at high N fertilization [11,12,13,14,15,16]. For soybean, a study investigating the effects of N fertilizer and duration of weed interference on yield reported that a competitive advantage for soybean was detected under N fertilization [17]. More information for soybean yield versus weed interference under N fertilization is needed to evaluate the effects of N fertilizer and weed interference on soybean yield. As multiple weed species are present in real fields, it is necessary to consider weed interference via a mixture of multiple weed species.
A rectangular hyperbolic model has provided a tool to evaluate how crop yield is influenced by weed interference under field conditions [8,18,19,20,21]. Kim et al. [14] developed a modified version of the rectangular hyperbolic model to describe crop yield loss caused by single-weed interference at a given N level. However, the model cannot account for the effect of multiple-weed interference on crop yield. Crop yield loss often occurs as a result of multiple-weed interference under commercial-scale production [8]. A model for the effects of N fertilizer and multiple-weed interference on soybean yield would be useful to support decision making for weed management in combination with a N fertilizer use scenario.
Therefore, this study was conducted to investigate the effects of N fertilization and multiple-weed interference on soybean yield and to model these effects based on the rectangular hyperbolic model [19].

2. Materials and Methods

2.1. Field Experiments

Field experiments were conducted in 2013 and 2014 to evaluate the effects of N fertilizer and multiple-weed interference on soybean yield in Bogatyrka, Primorsky krai, Russia (43.82° N, 131.6° E). The respective mean daily temperature and total rainfall during the growing season (May to October) were 16.7 °C and 557.1 mm in 2013 and 15.3 °C and 408.5 mm in 2014, respectively [8,22]. The soil in this region had a silty-loam texture with a cation exchange capacity of 22.6 cmol kg−1, organic matter content of 29.6 g kg−1, total N concentration of 1.58 g kg−1, available phosphorus (P) concentration of 18.2 mg kg−1, and a pH of 6.6 [8,22]. N fertilizer was applied in the form of urea (CO(NH2)2) at 0, 12, 24, and 36 kg N ha−1 on 27 May 2013 and 22 May 2014. The other basal fertilizers were applied yearly at a rate of 31–31 (kg ha−1) = PO4-K2O before N application.
Soybean (Glycine max cv. Heinong 48) was planted at a seeding rate of 80 kg ha−1 with a row width of 70 cm at different N fertilizer plots. Weed seedlings that emerged naturally within rows were manually thinned to a target plant density for individual weed species A. artemisiifolia, E. crus-galli, and B. syzigachne and different plant density combinations of these weeds up to the V6 stage of soybean in pre-marked plots. Other weed species were manually removed. The maximum plant densities and the number of pre-marked plots for competition between soybean and single weed species were as follows: A. artemisiifolia, 49, 24, 20, and 17 plants m−2 in 17, 11, 9, and 10 plots at 0, 12, 24, and 36 kg N ha−1 in 2013, respectively, and 49, 43, 43, and 45 plants m−2 in 17, 17, 17, and 14 plots in 2014; E. crus-galli, 91, 91, 16, and 23 plants m−2 in 9, 9, 11, and 13 plots in 2013, respectively, and 100, 139, 139, and 139 plants m−2 in 11, 13, 13, and 16 plots in 2014; B. syzigachne, 26, 71, 23, and 31 plants m−2 in 11, 11, 9, and 10 plots in 2013, respectively, and 110, 56, 70, and 100 plants m−2 in 18, 13, 15, and 16 plots in 2014 (Table S1). The maximum total plant densities and the number of pre-marked plots for competition between soybean and various combinations of three weed species were as follows: 91, 279, 304, and 31 plants m−2 in 39, 35, 33, and 37 plots at 0, 12, 24, and 36 kg N ha−1 in 2013, respectively, and 110, 139, 139, and 139 plants m−2 in 68, 67, 69, and 71 plots in 2014 (Table S1). In the main plot of each N fertilization, pre-marked plots were laid out in a completely randomized design with a single replicate. The size of each main plot was 100 m by 56 m, while the size of each pre-marked plot was 4.2 m by 2 m including a buffer area. The sampling area for harvest was 2.1 m by 1 m in each pre-marked plot, and three soybean rows were harvested by hand at maturity in October of each year. The soybean seed yield was adjusted to 14% moisture content after its seed weight and moisture content were measured at harvest.

2.2. Model Development

A rectangular hyperbolic model was used to describe the competition effects of three weed species in a single stand (single-weed interference) or in a mixture (multiple-weed interference) on soybean yield (Equation (1)) [8,19,20].
Y = Y 0 1 + β 1 ( X 1 + β 2 β 1 X 2 + β 3 β 1 X 3 )
where Y0 is weed-free soybean yield (Mg ha−1) and β is a measure of weed competitiveness (β1, β2, and β3 are the competitiveness of three weed species; 1 is a reference weed species). In Equation (1), β2/β1 and β3/β1 are the relative competitiveness of each species in order to convert its original density into a relative density, i.e., the so-called density equivalent [18]. In this study, A. artemisiifolia is the reference weed species. Thus, X1+(β2/β1)X2+( β3/β1)X3 is the total density equivalent, which is the sum of the relative densities of A. artemisiifolia, E. crus-galli, and B. syzigachne, calculated based on the density equivalent [18]. If N fertilizer is applied to a soybean field at different levels, two parameters Y0 and β for the rectangular hyperbolic model may change with applied N (i). Therefore, at each N level, Equation (1) can be rewritten as follows:
Y = Y 0 i 1 + β 1 i ( X 1 + β 2 i β 1 i X 2 + β 3 i β 1 i X 3 )
where Y0i and βi are the parameters for the i-th level of applied N. However, Equation (2) requires many parameters to predict soybean yield loss caused by weed interference for a wide range of N levels. Weed-free soybean yield (Y0i) may change with increasing N level as a function of the inverse quadratic curve since excess N may be toxic [14,23]. Weed competitiveness (βi) may change as a function of the exponential curve [14]. To describe the relationship between the parameters and applied N, inverse quadratic and exponential models were used as follows:
Y 0 i = a + b N 1 + c N + d N 2
β i = l m N
where a, b, c, and d are unknown parameters for weed-free soybean yield, l and m are unknown parameters for weed competitiveness, and N is the level of applied nitrogen. Therefore, the inverse quadratic and exponential models for two parameters Y0i and βi, respectively, were incorporated into Equation (2) to give a combined model (Equation (5)) as follows:
Y = a + b N 1 + c N + d N 2 1 + l 1 m 1 N ( X 1 + l 2 m 2 N l 1 m 1 N X 2 + l 3 m 3 N l 1 m 1 N X 3 )
where X1 + (l2m2N/l1m1N)X2 + (l3m3N/l1m1N)X3 is the total density equivalent (plants m−2).

2.3. Statistical Analyses

All measurements were initially subjected to analysis of variance (ANOVA). Non-linear regression analyses were then conducted to fit the rectangular hyperbolic model (Equation (2)) to the yield data, corresponding to the plant density combination of three weed species at different N levels. Inverse quadratic and exponential models were regressed on the parameters Y0i and βi affected by applied N, respectively. The performance of the models was evaluated by the pseudo-R2. Parameter estimates were compared between years using dummy variables [24]. Significant difference (5% level) between values was determined based on whether or not the confidence intervals of the dummy variables contain zero. Lack-of-fit of the rectangular hyperbolic model (Equation (2); Full model) was tested to check that the inverse quadratic and exponential models (Equations (3) and (4)) were appropriate for each parameter change with increasing N. The combined model (Equation (5); Reduced model) was compared with the Full model by calculating the F-value as follows:
F = R S S r R S S f d f r d f f R S S f d f f
where RSS and df are the residual sum of square and the degree of freedom, respectively, and r and f are the Reduced and Full models. If the F-value was lower than the tabulated F-value (5% level) with (dfrdff, dff) degrees of freedom, the Reduced model could be accepted. All the statistical analyses were conducted using Genstat [25].

3. Results

3.1. Effect of Nitrogen on Soybean-Weed Competition under Single Weed Interference

Soybean yield at different levels of N fertilizer was regressed on single species weed density using the rectangular hyperbolic model (Table 1). The Y0 and β parameter values at different levels of N fertilizer showed little or no difference between years within weed species (Table S2), so the two-year data for each weed species were pooled and the model was regressed on the pooled data (Table 1). Weed-free soybean yield (Y0) and weed competitiveness (β) values at no N application were 1.61 Mg ha−1 and 0.093, 1.55 Mg ha−1 and 0.086, and 1.56 Mg ha−1 and 0.046 for A. artemisiifolia, E. crus-galli, and B. syzigachne, respectively. The weed-free soybean yield value increased with increasing N level; approximately 1.77 Mg, 2.05 Mg, and 2.26 Mg ha−1 at 12, 24, and 36 kg N ha−1, respectively, equivalent to 12.7%, 30.6%, and 43.9% yield increases compared with the yield at no N application, 1.57 Mg ha−1 (Table 1). Thus, the weed-free soybean yield showed the inverse quadratic increase with N fertilizer (Figure S1). The weed competitiveness value also increased with increasing N level. For A. artemisiifolia, weed competitiveness value was 0.127, 0.203, and 0.529 at 12, 24, and 36 kg N ha−1, respectively, equivalent to approximately 1.4-fold, 2.2-fold, and 5.7-fold increases compared with the value at no N application; for E. crus-galli, 0.142, 0.384, and 0.369 at 12, 24, and 36 kg N ha−1, respectively, equivalent to approximately 1.7-fold, 4.5-fold, and 4.3-fold increases; for B. syzigachne, 0.063, 0.090, and 0.100 at 12, 24, and 36 kg N ha−1, respectively, equivalent to approximately 1.4-fold, 2.0-fold, and 2.2-fold increases (Table 1; Figure S2). Thus, the weed competitiveness of single species showed an exponential increase with N fertilizer (Table 2; Figure S2).

3.2. Effect of Nitrogen on Soybean-Weed Competition under Multiple Weed Interference

The competitiveness of individual weed species at different N levels was converted into a density equivalent by dividing with the competitiveness of A. artemisiifolia as the reference species (Table 1; Table S3). In the pooled data, the density equivalents of A. artemisiifolia, E. crus-galli, and B. syzigachne at no N application were 1.00, 0.92, and 0.49, respectively. These values for A. artemisiifolia, E. crus-galli, and B. syzigachne in N applied plots were 1.00, 1.12, and 0.50 at 12 kg N ha−1 and 1.00, 1.89, and 0.44 at 24 kg N ha−1 and 1.00, 0.70, and 0.19 at 36 kg N ha−1, respectively (Table S3). Various plant density combinations of the three weed species were converted into a total density equivalent by multiplying the actual density of each weed with its density equivalent. Soybean yield at different levels of N fertilizer was regressed on the total density equivalent of multiple weed species using the rectangular hyperbolic model (Figure 1). A pairwise comparison of parameter estimates showed that the Y0 and β parameter values at different levels of N fertilizer showed little or no difference between years (Table S2), so the two-year data were pooled and the model was regressed on the pooled data (Figure 1). Weed-free soybean yield and multiple-weed competitiveness values at no N application were 1.63 Mg ha−1 and 0.079, respectively (Table 3). These values increased as follows: weed-free soybean yield, 1.75 Mg, 2.07 Mg, and 2.25 Mg ha−1 at 12, 24, and 36 kg N ha−1, respectively, equivalent to approximately 7.4%, 27.0%, and 38.0% yield increases compared with the yield at no N application; multiple-weed competitiveness, 0.121, 0.295, and 0.463 at 12, 24, and 36 kg N ha−1, respectively, equivalent to approximately 1.5-fold, 3.7-fold, and 5.9-fold increases (Table 3).

3.3. Responses of Weed-Free Soybean Yield and Multiple-Weed Competitiveness to Nitrogen

Inverse quadratic and exponential models were regressed on the weed-free soybean yield and multiple-weed competitiveness affected by applied N in the pooled data, respectively (Figure 2). These models showed good performance with high pseudo-R2. Weed-free soybean yield slightly increased with increasing N level from 0 to 36 kg N ha−1, and multiple-weed competitiveness exponentially increased with an increase of N up to 36 kg N ha−1. Thus, inverse quadratic and exponential models can be sufficient to describe N effects on weed-free soybean yield and multiple-weed competitiveness, respectively.

3.4. Combined Model for Soybean Yield Affected by Multiple-Weed Competitiveness and Applied Nitrogen

The rectangular hyperbolic model was incorporated with inverse quadratic and exponential models to describe the relationships between weed-free soybean yield and applied N, and between multiple-weed competitiveness and applied N, respectively. In order to check that inverse quadratic and exponential models were appropriate, the combined model was compared with the rectangular hyperbolic model by lack-of-fit test. Statistical analyses revealed that there was no significant difference between the rectangular hyperbolic model (Full model) and the combined model (Reduced model) to describe soybean yield affected by both multiple-weed competitiveness and applied N (Table S4). Thus, the rectangular hyperbolic model can be modified to the combined model by replacing the parameters Y0 and β with the inverse quadratic and exponential models, respectively.
Using the combined model (Equation (5)) and estimated parameters (Table 4), soybean yield was predicted in Figure 3. A steep decline of soybean yield was predicted with increasing total density equivalent, particularly at 36 kg N ha−1 although weed-free soybean yield increased with increasing N up to 36 kg N ha−1. Weed-free soybean yield value was predicted to be 1.65 Mg, 1.84 Mg, 2.03 Mg, and 2.23 Mg ha−1 at 0, 12, 24, and 36 kg N ha−1, respectively. However, when interfered with multiple-weed species at 12 plants m−2, soybean yield was reduced to 52%, 33%, 20%, and 11% of weed-free plot at 0, 12, 24, and 36 kg N ha−1, respectively. Thus, in a N applied condition, proper weed management should be made to achieve maximum soybean yield in the Far-Eastern regions of Russia.

4. Discussion

N fertilizer was applied to promote early soybean growth and thus to secure soybean yield due to the short frost-free period in the region of Primorsky krai, Russia. N application at 36 kg N ha−1 was effective for high soybean production in the region. In our study, soybean yield increased to 2.2 Mg ha−1 as N application increased to 36 kg N ha−1 over years at the Bogatyrka field. In similar weather conditions (Northern Great Plains in USA), low N application of 16 kg N ha−1 at planting enhanced early vegetative soybean growth and resultant seed yield compared to untreated plot [2]. Similarly, for high soybean production, N application in a range from 22.5 to 30 kg N ha−1 was effective in Northeast China, closely linked to the Far-Eastern regions [3,4]. Other regional studies also reported a high increase of soybean yield with increasing N level up to 36 kg N ha−1 [26,27,28]. However, soybean yield can be reduced more by enhanced weed competitiveness under N-fertilized conditions where weed biomass increases but resultant growth of soybean decreases. Particularly, N application up to 36 kg N ha−1 increased multiple-weed competitiveness to soybean, resulting in a significant yield loss more than 89% at 12 plants m−2 of the total density equivalent (Figure 3). Previous studies also indicated that N fertilizer increases weed competitiveness against crops including rice [11], wheat [14,29] and oilseed rape [30]. The N effect on weed competitiveness with crops might be dependent on interference relationships, since N application changes both crop growth and weed biomass. Our study indicates that A. artemisiifolia and E. crus-galli more increased their competitiveness over soybean by five-fold and three-fold increases compared with the value of B. syzigachne, respectively, as N application increased up to 36 kg N ha−1 (Table 1; Figure S2). Thus, when A. artemisiifolia and E. crus-galli are not adequately managed, N application can increase their competitiveness; rather, soybean yield can be offset by nitrophilous weed species at high N fertilization.
The combined rectangular hyperbolic model used in our study to describe N effects on weed-free soybean yield and multiple-weed competitiveness is useful in supporting decision making for weed management in combination with a N fertilizer use scenario. The combined model suggests weed threshold values for soybean yield, in consideration of weed interference under N-fertilized conditions. If target soybean yield is decided at a given N level, economic thresholds (ETs) of multiple-weed species can be predicted by comparing the cost of controlling weed species with the benefit gained by herbicide application (e.g., Song et al. [8], Cousens [31], Zanin et al. [32]). In this region, the sequential application of herbicides, S-metolachlor followed by bentazon + imazamox, is recommended for respective pre- and post-emergence weed control (e.g., Song et al. [22]). Soybean price and herbicide cost are US $650 Mg−1 and US $87.2 ha−1, respectively, and N costs are US $16.2, US $32.3, and US $48.5 ha−1 for 12, 24, and 36 kg N ha−1, respectively (e.g., Song et al. [22]). Using these values, the proportion of yield loss caused by unit total density equivalent obtained in our study, and assuming a 90% herbicide efficacy, the calculated ETs are 1.46, 0.90, 0.57, and 0.38 plants m−2 at 0, 12, 24, and 36 kg N ha−1, respectively (Table 5). Therefore, thorough weed management should be made to achieve maximum soybean yield in the region in combination with N application. In a combination of 36 N kg ha−1, the sequential application of herbicides at ET greater than 0.38 plants m−2 may secure the maximum soybean yield of 2.2 Mg ha−1.

5. Conclusions

The combined model developed here can support decision making for weed management under different N uses in soybean cultivation across Primorsky krai and neighboring regions. It can be used to advise both herbicide and N applications for a target soybean yield. Nonetheless, the combined model and parameter estimates need to be further validated for a broad range of N rate, years, and locations.

Supplementary Materials

The following are available online at https://www.mdpi.com/2073-4395/11/3/515/s1, Table S1: Weed species and densities for soybean-weed competition at different nitrogen levels in 2013 and 2014, Table S2: A pairwise comparison of parameter estimates between years using dummy variables at different nitrogen levels, Table S3: Density equivalents of Echinochloa crus-galli and Beckmannia syzigachne in comparison with Ambrosia artemisiifolia as the reference weed species at different nitrogen levels in 2013 and 2014 and for the pooled two-year data, Table S4: Summary of non-linear regression analysis and F-test to compare models for soybean yield as a function of the total density equivalent and applied nitrogen in the pooled two-year data, Figure S1: Weed-free soybean yield (Y0) as a function of applied nitrogen in 2013 (●) and 2014 (□) and for the pooled two-year (—) data presented in Table 1. The lines are fitted values calculated using inverse quadratic model, Figure S2: Weed competitiveness of Ambrosia artemisiifolia, Echinochloa crus-galli, and Beckmannia syzigachne as a function of applied nitrogen in 2013 (···) and 2014 (---) and for the pooled two-year (—) data presented in Table 1. The lines are fitted values calculated using exponential model.

Author Contributions

Conceptualization, D.-S.K.; methodology, D.-S.K. and J.-S.S.; formal analysis, J.-S.S.; investigation, J.-S.S., J.-H.I., J.-W.K., D.-G.K., Y.L., M.-J.Y. and S.-H.L.; data curation, J.-S.S.; Writing—Original draft preparation, J.-S.S.; Writing—Review and editing, D.-S.K. and J.-S.S.; visualization, J.-S.S.; supervision, D.-S.K.; project administration, D.-S.K.; funding acquisition, D.-S.K. All authors have read and agreed to the published version of the manuscript.

Funding

This work was carried out with the support of the “Cooperative Research Program for Agriculture Science & Technology Development (Project No. PJ015723)”, Rural Development Administration, Republic of Korea. J.-S.S. is supported by the R&D Program of the “Plasma Advanced Technology for Agriculture and Food (Plasma Farming, Project No. 1711124797)” through the Korea Institute of Fusion Energy funded by the Government funds, Republic of Korea.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Salvagiotti, F.; Cassman, K.G.; Specht, J.E.; Walters, D.T.; Weiss, A.; Dobermann, A. Nitrogen uptake, fixation and response to fertilizer N in soybeans: A review. Field Crops Res. 2008, 108, 1–13. [Google Scholar] [CrossRef] [Green Version]
  2. Osborne, S.L.; Riedell, W.E. Starter nitrogen fertilizer impact on soybean yield and quality in the Northern Great Plains. Agron. J. 2006, 98, 1569–1574. [Google Scholar] [CrossRef] [Green Version]
  3. Li, J.R. The scientific knowledge of applying fertilizer technique for soybean. J. Jili Agric. Sci. Technol. Coll. 2005, 14, 32–35. (In Chinese) [Google Scholar]
  4. Liu, X.; Jin, J.; Wang, G.; Herbert, S.J. Soybean yield physiology and development of high-yielding practices in Northeast China. Field Crops Res. 2008, 105, 157–171. [Google Scholar] [CrossRef]
  5. Development of Agriculture and Market Regulation of Agricultural Products, Raw Materials and Food: Improving Living Standards of Rural Population of Primorsky Krai. Available online: http://www.agrodv.ru/ (accessed on 18 May 2016).
  6. China June Crop Area and Production Report. Available online: http://www.informaecon.com/samplereports/ChinaSampleReport.pdf (accessed on 23 May 2016).
  7. Sinegovskii, M.; Yuan, S.; Sinegovskaya, V.; Han, T. Current status of the soybean industry and research in the Russian Federation. Soybean Sci. 2018, 37, 1–7. [Google Scholar]
  8. Song, J.S.; Kim, J.W.; Im, J.H.; Lee, K.J.; Lee, B.W.; Kim, D.S. The effects of single- and multiple-weed interference on soybean yield in the Far-Eastern region of Russia. Weed Sci. 2017, 65, 371–380. [Google Scholar] [CrossRef]
  9. Chauhan, B.S.; Abugho, S.B. Effects of water regime, nitrogen fertilization, and rice plant density on growth and reproduction of lowland weed Echinochloa crus-galli. Crop. Prot. 2013, 54, 142–147. [Google Scholar] [CrossRef]
  10. Leskovsek, R.; Datta, A.; Simoncic, A.; Knezevic, S.Z. Influence of nitrogen and plant density on the growth and seed production of common ragweed (Ambrosia artemisiifolia L.). J. Pest Sci. 2012, 85, 527–539. [Google Scholar] [CrossRef]
  11. Ampong-Nyarko, K.; De Datta, S.K. Effects of nitrogen application on growth, nitrogen use efficiency and rice-weed interaction. Weed Res. 1993, 33, 269–276. [Google Scholar] [CrossRef]
  12. Buchanan, G.A.; McLaughlin, R.D. Influence of nitrogen on weed competition in cotton. Weed Sci. 1975, 23, 324–328. [Google Scholar] [CrossRef]
  13. Dhima, K.V.; Eleftherohorinos, I.G. Influence of nitrogen on competition between winter cereals and sterile oat. Weed Sci. 2001, 49, 77–82. [Google Scholar] [CrossRef]
  14. Kim, D.S.; Marshall, E.J.P.; Brain, P.; Caseley, J.C. Modelling the effects of sub-lethal doses of herbicide and nitrogen fertilizer on crop-weed competition. Weed Res. 2006, 46, 492–502. [Google Scholar] [CrossRef]
  15. Naderi, R.; Ghadiri, H. Competition of wild mustard (Sinapis arvense L.) densities with rapeseed (Brassica napus L.) under different levels of nitrogen fertilizer. J. Agric. Sci. Tech. 2011, 13, 45–51. [Google Scholar]
  16. Wells, G.J. Annual weed competition in wheat crops: The effect of weed density and applied nitrogen. Weed Res. 1979, 19, 185–191. [Google Scholar] [CrossRef]
  17. Shafagh-Kolvanagh, J.; Zehtab-Salmasi, S.; Javanshir, A.; Moghaddam, M.; Nasab, A.D.M. Effects of nitrogen and duration of weed interference on grain yield and SPAD (chlorophyll) value of soybean (Glycine max (L.) Merrill.). J. Food Agric. Environ. 2008, 6, 368–373. [Google Scholar]
  18. Berti, A.; Zanin, G. Density equivalent: A method for forecasting yield loss caused by mixed weed populations. Weed Res. 1994, 34, 327–332. [Google Scholar] [CrossRef]
  19. Cousens, R. A simple model relating yield loss to weed density. Ann. Appl. Biol. 1985, 107, 239–252. [Google Scholar] [CrossRef]
  20. Kim, D.S.; Marshall, E.J.P.; Caseley, J.C.; Brain, P. Modelling interactions between herbicide dose and multiple weed species interference in crop-weed competition. Weed Res. 2006, 46, 175–184. [Google Scholar] [CrossRef]
  21. Lindquist, J.L.; Dieleman, J.A.; Mortensen, D.A.; Johnson, G.A.; Wyse-Pester, D.Y. Economic importance of managing spatially heterogeneous weed population. Weed Technol. 1998, 12, 7–13. [Google Scholar] [CrossRef]
  22. Song, J.S.; Chung, J.H.; Lee, K.J.; Kwon, J.; Kim, J.W.; Im, J.H.; Kim, D.S. Herbicide-based weed management for soybean production in the Far Eastern region of Russia. Agronomy 2020, 10, 1823. [Google Scholar] [CrossRef]
  23. Goyal, S.S.; Huffaker, R.C. Nitrogen toxicity in plants. In Nitrogen in Crop Production; Hauck, R.D., Ed.; ASA-CSSA-SSSA Inc.: Madison, WI, USA, 1984; pp. 97–118. [Google Scholar]
  24. Chism, W.J.; Birch, J.B.; Bingham, S.W. Nonlinear regressions for analyzing growth stage and quinclorac interactions. Weed Technol. 1992, 6, 898–903. [Google Scholar] [CrossRef]
  25. Genstat Committee. Reference Manual (Genstat Release 6.1); VSN International: Oxford, UK, 2002. [Google Scholar]
  26. Afza, R.; Hardason, G.; Zapata, F.; Danso, S.K.A. Effects of delayed soil and foliar N fertilization on yield and N2 fixation of soybean. Plant Soil 1987, 97, 361–368. [Google Scholar] [CrossRef]
  27. Caliskan, S.; Ozkaya, I.; Caliskan, M.E.; Arslan, M. The effects of nitrogen and iron fertilization on growth, yield and fertilizer use efficiency of soybean in a Mediterranean-type soil. Field Crops Res. 2008, 108, 126–132. [Google Scholar] [CrossRef]
  28. Taylor, R.S.; Weaver, D.B.; Wood, C.W.; van Santen, E. Nitrogen application increases yield and early dry matter accumulation in late-planted soybean. Crop Sci. 2005, 45, 854–858. [Google Scholar] [CrossRef]
  29. Blackshaw, R.E.; Brandt, R.N. Nitrogen fertilizer rate effects on weed competitiveness is species dependent. Weed Sci. 2008, 56, 743–747. [Google Scholar] [CrossRef]
  30. Kazemeini, S.A.; Naderi, R.; Aliabadi, H.K. Effects of different densities of wild oat (Avena fatua L.) and nitrogen rates on oilseed rape (Brassica napus L.) yield. J. Ecol. Environ. 2013, 36, 167–172. [Google Scholar] [CrossRef] [Green Version]
  31. Cousens, R. Theory and reality of weed control thresholds. Plant Prot. Q. 1987, 2, 13–20. [Google Scholar]
  32. Zanin, G.; Berti, A.; Toniolo, L. Estimation of economic thresholds for weed control in winter wheat. Weed Res. 1993, 33, 459–467. [Google Scholar] [CrossRef]
Figure 1. Soybean yield as a function of the total density equivalent at different nitrogen levels in 2013 (···) and 2014 (---) and for the pooled two-year (—) data. The lines are fitted values calculated using the rectangular hyperbolic model and parameter estimates (Table 3).
Figure 1. Soybean yield as a function of the total density equivalent at different nitrogen levels in 2013 (···) and 2014 (---) and for the pooled two-year (—) data. The lines are fitted values calculated using the rectangular hyperbolic model and parameter estimates (Table 3).
Agronomy 11 00515 g001
Figure 2. Weed-free soybean yield (a) and multiple-weed competitiveness (b) as a function of applied nitrogen in the pooled two-year data.
Figure 2. Weed-free soybean yield (a) and multiple-weed competitiveness (b) as a function of applied nitrogen in the pooled two-year data.
Agronomy 11 00515 g002
Figure 3. Soybean yield as a function of the total density equivalent and applied nitrogen in the pooled two-year data. The mesh is fitted values calculated using the combined model and parameter estimates (Table 4).
Figure 3. Soybean yield as a function of the total density equivalent and applied nitrogen in the pooled two-year data. The mesh is fitted values calculated using the combined model and parameter estimates (Table 4).
Agronomy 11 00515 g003
Table 1. Parameter estimates for the rectangular hyperbolic model for the regression of soybean yield as a result of single-weed interference caused by Ambrosia artemisiifolia (A), Echinochloa crus-galli (B), and Beckmannia syzigachne (C) at different nitrogen levels in 2013 and 2014 and for the pooled two-year data. The numbers in parentheses are standard errors.
Table 1. Parameter estimates for the rectangular hyperbolic model for the regression of soybean yield as a result of single-weed interference caused by Ambrosia artemisiifolia (A), Echinochloa crus-galli (B), and Beckmannia syzigachne (C) at different nitrogen levels in 2013 and 2014 and for the pooled two-year data. The numbers in parentheses are standard errors.
Nitrogen (kg N ha−1)Weed SpeciesParameter EstimatesPseudo R2
Y0β
20132014Pooled20132014Pooled20132014Pooled
0A1.621.441.610.0520.0850.0930.860.910.81
(0.095)(0.086)(0.096)(0.0104)(0.0132)(0.0159)
B1.521.391.550.0660.0860.0860.840.970.89
(0.169)(0.058)(0.085)(0.0259)(0.0114)(0.0177)
C1.641.401.560.0470.0630.0460.800.920.65
(0.101)(0.079)(0.131)(0.0122)(0.0089)(0.0130)
12A1.731.811.770.1300.1180.1270.860.840.84
(0.122)(0.138)(0.090)(0.0324)(0.0249)(0.0206)
B1.741.821.780.1350.1480.1420.960.860.90
(0.085)(0.163)(0.093)(0.0353)(0.0582)(0.0354)
C1.721.821.750.0400.0890.0630.890.810.80
(0.107)(0.157)(0.102)(0.0079)(0.0251)(0.0125)
24A2.101.962.040.2320.2200.2030.950.980.94
(0.104)(0.048)(0.070)(0.0396)(0.0175)(0.0236)
B2.071.962.020.4130.3070.3840.860.990.94
(0.170)(0.041)(0.080)(0.1170)(0.0477)(0.0664)
C2.181.962.080.0820.0930.0900.880.950.92
(0.141)(0.073)(0.077)(0.0205)(0.0106)(0.0112)
36A2.392.162.260.5450.5940.5290.980.980.98
(0.097)(0.058)(0.056)(0.0839)(0.1040)(0.0603)
B2.402.162.250.4380.2910.3690.910.940.92
(0.184)(0.105)(0.096)(0.1010)(0.0568)(0.0557)
C2.412.172.270.1130.0930.1000.930.960.94
(0.157)(0.074)(0.075)(0.0218)(0.0095)(0.0100)
Table 2. Parameter estimates for the exponential model for the regression of weed competitiveness of Ambrosia artemisiifolia (β1i), Echinochloa crus-galli (β2i), and Beckmannia syzigachne (β3i) as a function of applied nitrogen in 2013 and 2014 and for the pooled two-year data presented in Table 1. The numbers in parentheses are standard errors.
Table 2. Parameter estimates for the exponential model for the regression of weed competitiveness of Ambrosia artemisiifolia (β1i), Echinochloa crus-galli (β2i), and Beckmannia syzigachne (β3i) as a function of applied nitrogen in 2013 and 2014 and for the pooled two-year data presented in Table 1. The numbers in parentheses are standard errors.
YearParameter β
Ambrosia artemisiifolia (β1i)Echinochloa crus-galli (β2i)Beckmannia syzigachne (β3i)
l1m1R2l2m2R2l3m3R2
20130.0521.0670.990.1141.0410.720.0371.0320.87
(0.0071)(0.0044)(0.0606)(0.0178)(0.0082)(0.0077)
20140.0461.0730.980.1191.0280.650.0721.0090.49
(0.0144)(0.0099)(0.0438)(0.0130)(0.0082)(0.0043)
Pooled0.0531.0650.970.1241.0340.760.0501.0210.94
(0.0185)(0.0111)(0.0598)(0.0165)(0.0053)(0.0039)
Table 3. Parameter estimates for the rectangular hyperbolic model for the regression of soybean yield as a result of multiple weed interference expressed as the total density equivalent at different nitrogen levels in 2013 and 2014 and for the pooled two-year data. The numbers in parentheses are the standard errors.
Table 3. Parameter estimates for the rectangular hyperbolic model for the regression of soybean yield as a result of multiple weed interference expressed as the total density equivalent at different nitrogen levels in 2013 and 2014 and for the pooled two-year data. The numbers in parentheses are the standard errors.
Nitrogen (kg N ha−1)ParametersPseudo R2
Y0β
20132014Pooled20132014Pooled20132014Pooled
01.651.381.630.0660.0760.0790.750.730.81
(0.105)(0.093)(0.071)(0.0129)(0.0097)(0.0084)
121.711.791.750.0940.1410.1210.870.800.82
(0.092)(0.106)(0.073)(0.0129)(0.0182)(0.0117)
242.111.992.070.2540.3130.2950.880.890.89
(0.110)(0.080)(0.066)(0.0348)(0.0256)(0.0213)
362.402.142.250.4720.4540.4630.850.870.86
(0.153)(0.093)(0.082)(0.0630)(0.0391)(0.0341)
Table 4. Parameter estimates for the combined model for the regression of soybean yield as a result of multiple-weed interference expressed as the total density equivalent and applied nitrogen in the pooled two-year data. The numbers in parentheses are standard errors.
Table 4. Parameter estimates for the combined model for the regression of soybean yield as a result of multiple-weed interference expressed as the total density equivalent and applied nitrogen in the pooled two-year data. The numbers in parentheses are standard errors.
ParametersPseudo R2
Y0iβi
abcdlm
1.650.0162.96 × 10−181.05 × 10−200.0811.050.86
(0.064)(0.357)(0.207)(0.0013)(0.0062)(0.003)
Table 5. Economic thresholds (ETs) of multiple-weed species including Ambrosia artemisiifolia, Echinochloa crus-galli, and Beckmannia syzigachne at different nitrogen levels in the pooled data.
Table 5. Economic thresholds (ETs) of multiple-weed species including Ambrosia artemisiifolia, Echinochloa crus-galli, and Beckmannia syzigachne at different nitrogen levels in the pooled data.
Nitrogen (kg N ha−1)Parameters and ETs a
Cn
($ ha−1)
Ch
($ ha−1)
Ca
($ ha−1)
Yo
(Mg ha−1)
P
($ Mg−1)
LHET b
0087.217.71.656500.0750.901.46
1216.287.217.71.846500.1250.900.90
2432.387.217.72.036500.2020.900.57
3648.587.217.72.226500.3110.900.38
a ET = C n + C h + C a Y o P β H (Cn is urea cost (US $ ha−1); Ch is herbicide cost (US $ ha−1); Ca is application cost (US $ ha−1); Yo is weed-free soybean yield (Mg ha−1); P is value per unit soybean (US $ Mg−1); L is a proportional loss per unit total density equivalent; H is herbicide efficacy (% weed control/100). Parameters Ch, Ca, and P were obtained from Song et al. [22]. At a given nitrogen level, parameters Yo and L were calculated using the combined rectangular hyperbolic model and parameter estimates described in Table 4. b ET of multiple weed species at different nitrogen levels is expressed as the total density equivalent m−2.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Song, J.-S.; Im, J.-H.; Kim, J.-W.; Kim, D.-G.; Lim, Y.; Yook, M.-J.; Lim, S.-H.; Kim, D.-S. Modeling the Effects of Nitrogen Fertilizer and Multiple Weed Interference on Soybean Yield. Agronomy 2021, 11, 515. https://doi.org/10.3390/agronomy11030515

AMA Style

Song J-S, Im J-H, Kim J-W, Kim D-G, Lim Y, Yook M-J, Lim S-H, Kim D-S. Modeling the Effects of Nitrogen Fertilizer and Multiple Weed Interference on Soybean Yield. Agronomy. 2021; 11(3):515. https://doi.org/10.3390/agronomy11030515

Chicago/Turabian Style

Song, Jong-Seok, Ji-Hoon Im, Jin-Won Kim, Dong-Gil Kim, Yeonhwa Lim, Min-Jung Yook, Soo-Hyun Lim, and Do-Soon Kim. 2021. "Modeling the Effects of Nitrogen Fertilizer and Multiple Weed Interference on Soybean Yield" Agronomy 11, no. 3: 515. https://doi.org/10.3390/agronomy11030515

APA Style

Song, J. -S., Im, J. -H., Kim, J. -W., Kim, D. -G., Lim, Y., Yook, M. -J., Lim, S. -H., & Kim, D. -S. (2021). Modeling the Effects of Nitrogen Fertilizer and Multiple Weed Interference on Soybean Yield. Agronomy, 11(3), 515. https://doi.org/10.3390/agronomy11030515

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop