Skip to main content

ORIGINAL RESEARCH article

Front. Plant Sci., 16 July 2021
Sec. Crop and Product Physiology

Modeling the Effects of Elevated Temperature and Weed Interference on Rice Grain Yield

  • 1Department of Agriculture, Forestry, and Bioresources, Research Institute of Agriculture and Life Sciences, College of Agriculture and Life Sciences, Seoul National University, Seoul, South Korea
  • 2Institute of Plasma Technology, Korea Institute of Fusion Energy, Gunsan, South Korea
  • 3Mushroom Research Division, Department of Herbal Crop Research, National Institute of Horticultural and Herbal Science, Rural Development Administration, Eumseong, South Korea
  • 4Crop Protection Division, Department of Agro-Food Safety and Crop Protection, National Institute of Agricultural Sciences, Rural Development Administration, Wanju, South Korea

A 3-year phytotron study was conducted in Suwon (37.27°N, 126.99°E), Korea, to evaluate and model the effects of elevated temperature on rice-weed competition. The dry weight and the number of panicles in rice were the most susceptible components to weed interference during the early growth of rice, regardless of weed species, while other yield components, including the number of grains, % ripened grain, and 1000-grain weight, were more susceptible to elevated temperature. A rectangular hyperbolic model well demonstrated that rice grain yield was affected by weed interference under elevated temperature, showing that the competitiveness of late watergrass (Echinochloa oryzicola) and water chestnut (Eleocharis kuroguwai) increased under elevated temperature conditions. Quadratic and linear models well described the effects of elevated temperature on the weed-free rice grain yield and weed competitiveness values of the rectangular hyperbolic model for the two weed species, respectively. Thus, a combined rectangular hyperbolic model incorporated with the quadratic and linear models well demonstrated the effects of elevated temperature and weed interference on rice grain yield across years. Using the combined model and estimated parameters, the rice grain yields were estimated to be 58.9, 48.5, 41.3, and 35.9% of the yields under weed-free conditions for 80 plants m−2 of late watergrass and 86.8, 64.3, 51.1, and 42.3% of the yields under weed-free conditions for 80 plants m−2 of water chestnut at 1,300, 1,500, 1,700, and 1,900°C·days of accumulated growing degree days (GDD; from transplanting to flowering, 89 days), respectively. The combined model developed in this study can provide an empirical description of both the elevated temperature and weed interference effects on rice yield and can be used for predicting rice grain yields due to weed interference under future elevated temperature conditions.

Introduction

According to the Representative Concentration Pathways (RCP) scenarios presented by the fifth assessment report (AR5) of the Intergovernmental Panel on Climate Change (IPCC), anthropogenic climate change challenges current and future global crop production due to the global warming exceeding 1.4°C (RCP6.0) and 4.8°C (RCP8.5) of mean air temperature in the 1980s by the 2080s, along with unprecedented extreme heat stress and high potential CO2 fertilization effects (Intergovernmental Panel on Climate Change [(IPCC)], 2013; Deryng et al., 2014). Rice is mainly cultivated in areas with tropical and temperate climates located between 53° northern latitude (Heilongjiang, China) and 35° southern latitude (New South Wales, Australia) (Koo et al., 2013; Kraehmer et al., 2016). The projected warming by the 2080s will have more serious consequences for rice production. For example, exposure to extreme events will reduce rice yield by damaging the development of grain spikelets during the anthesis and grain filling stages of rice (Welch et al., 2010; Nguyen et al., 2014; Sánchez et al., 2014; Jung et al., 2015; Lee et al., 2015). Even worse, the projected warming will negatively affect rice yield with increasing weed competitiveness due to more enhanced weed biomass compared with the biomass of rice (Alberto et al., 1996).

Late watergrass (Echinochloa oryzicola Vasing.) and water chestnut (Eleocharis kuroguwai Ohwi) are troublesome weed species that inhabit rice paddy fields, mainly in East Asian countries, including Korea, Japan, and Vietnam (Kraehmer et al., 2016; Song et al., 2016, 2017a). These weed species cause significant yield loss of rice due to the higher canopy height of late watergrass and stronger root competitiveness of water chestnut than those of rice during the entire growing season (Yamasue, 2001; Moon et al., 2010; Song et al., 2016). Even when herbicides are applied, late watergrass and water chestnut cannot always be controlled due to herbicide-resistant late watergrass (Fischer et al., 2000; Song et al., 2017a) and the tuberization characteristics of water chestnut in deep soil (Watanabe, 2011; Moon et al., 2014; Song et al., 2016). Late watergrass is a representative C4 weed species that can adapt to some environmental changes in rice cultivation, showing various adaptive characteristics, such as an enhanced photosynthetic capacity under high temperature (Bouhache and Bayer, 1993), enhanced tolerance to salinity (Nguyen et al., 2005), and germination and seedling growth under anaerobic conditions (Kennedy et al., 1980; Nah et al., 2015). Water chestnut can also adapt to such environmental conditions, although it is a C3 weed species (Ueno and Takeda, 1992) that shows tuber propagation under submerged soil conditions (Kobayashi and Ueki, 1983; Chun and Shin, 1994) and fast growth and high biomass under elevated temperatures (Kim et al., 2010). If these weed species are not properly managed under future climate change, then they will become more competitive against rice, resulting in greater yield losses than under current climate conditions.

Many studies have focused on the potential effects of climate change on rice production (Bachelet and Gay, 1993; Kim et al., 1996; Matsui et al., 1997; Matthews et al., 1997; Sheehy et al., 2006a; Mahajan et al., 2012; Nguyen et al., 2014) and paddy weed adaptation (Kim et al., 2010; Park et al., 2010; Rodenburg et al., 2011; Bir et al., 2014; Ramesh et al., 2017), but, relatively, few studies have examined the impacts of climate change on weed-rice competition. Previous studies reported the effects of elevated CO2 concentrations on the competitive interactions between weeds and rice (Alberto et al., 1996; Ziska et al., 2010; Zeng et al., 2011) and the effects of elevated CO2 or air temperature on the interactions between weeds and other crops, including soybeans (Flint and Patterson, 1983; Patterson et al., 1984; Ziska, 2000, 2001; Davis and Ainsworth, 2012), wheat (Thompson and Woodward, 1994), cotton (Flint et al., 1983), and sorghum (Ziska, 2003). A modeling study on the impacts of climate change and weed interference on crop yield under field conditions has not been reported to our knowledge (Ramesh et al., 2017). Late watergrass and water chestnut are dominant in rice cultivation in Korea and Japan but have not been investigated for their competitive effects on rice yield under future climate change. The evaluation of both the effects of climate change and weed interference on rice is essential to support decision-making for weed management and to establish effective weed management practices under the changing climate conditions. Although climate change can be induced by both elevated temperature and elevated CO2, it is difficult to evaluate and model both effects on rice-weed competition. Considering that plant growth is more sensitive to temperature than elevated CO2, it will be more effective for establishing weed management practices in rice under future climate conditions to evaluate and model the effects of elevated temperature on rice-weed competition.

Therefore, this study was conducted to evaluate the effects of elevated temperature on late watergrass-rice and water chestnut-rice competition and to develop a model to predict rice yield in the presence of weed interference under elevated temperature conditions.

Materials and Methods

Pot Experiments Under Phytotron

Pot experiments under phytotrons were conducted in 2014, 2015, and 2016 at the experimental farm station of Seoul National University, Suwon (37.27°N, 126.99°E), Korea, to evaluate the temperature effects on late watergrass-rice and water chestnut-rice competition. As shown in Table 1, four sunlit temperature-controlled phytotrons were adjusted to the target temperature conditions, ambient (A), A+1.5, A+3.0, and A+5.0°C, by using a heating and ventilation control system linked to a CR10x data logger (Camp-bell Sci., USA) (Nguyen et al., 2014; Jung et al., 2015; Kim et al., 2017; Kim and Lee, 2019). The ambient temperature was maintained with no sidewalls and only plastic film-covered roofs in a plastic house. The ambient temperature differed in different years (Supplementary Figure 1); thus, it was necessary to convert the daily mean air temperature (T) into accumulated growing degree days (GDD) based on the base temperature (Tbase) (McMaster and Wilhelm, 1997) as follows:

GDD=(T-Tbase)×days    (1)

The base temperature (Tbase) was taken as 10°C, which is commonly used in Korean rice-growing systems (e.g., Lee et al., 2015). Twenty-day-old rice seedlings (Oryza sativa cv. Chucheong) were transplanted into 30-cm-deep plastic pots (0.1 m2), containing paddy soil at a density of 30 hills m−2, which is equivalent to ~90 rice seedlings m−2, on June 1 over 3 years. The paddy soil in the pot had a loam texture with an organic matter content of 19.5 g kg−1, total nitrogen concentration of 1.10 g kg−1, NH4+-N concentration of 2.29 mg kg−1, and NO3--N concentration of 288.97 mg kg−1. Late water grass and water chestnut were manually sown on the paddy soil immediately after harrowing, randomly thinned to the target plant density for individual weed species immediately after rice transplanting and grown around rice plants until harvest at each of the four sunlit temperature-controlled phytotrons. The target densities for competition between rice and weed species were 0, 10, 20, 40, and 80 plants m−2 in 2014 and 0, 20, 40, 80, and 160 plants m−2 in 2015 and 2016 for both weed species. Other weed species were removed manually. Fertilization was performed with the basal dressing of N–P2O5-K2O = 55–45–45 kg ha−1 before harrowing and two sequential top-dressings of 27.5 kg ha−1 of N at both rice tillering and panicle initiation stages. Other cultivation practices, such as irrigation and pest management, were performed, following the standard practice of the Rural Development Administration of Korea (RDA) (Rural Development Administration [RDA], 2014). All the pots were arranged in a completely randomized block design with three replicates in each phytotron.

TABLE 1
www.frontiersin.org

Table 1. Mean air temperature during the period from transplanting (June 1) to tillering (August 10), during the period from transplanting (June 1) to flowering (August 28), and during the whole growth period (June 1 to October 15) at four sunlit temperature-controlled phytotrons adjusted to the target temperature conditions, such as ambient (A), A+1.5°C, A+3.0°C, and A+5.0°C, in 2014, 2015, and 2016.

In October of each year, rice was harvested at maturity to assess the yield components, including dry weight, number of panicles, number of grains, percentage of ripened grains, and thousand-grain weight. The grain weight and moisture content were also assessed, and the grain yield was calculated by adjusting the grain moisture content to 12%.

Model Development

A rectangular hyperbolic model (equation 2; Cousens, 1985; Kim et al., 2002; Song et al., 2017b) was used to describe the relationship between rice grain yield (Y) and initial single weed density (X) under different air temperature conditions.

Y=Yo1+βX    (2)

where Yo is the weed-free rice grain yield (t ha−1), and β is a measure of weed competitiveness. When rice and a single weed were grown in a mixture under elevated temperature, the parameters Yo and β for the model can change with the level of elevated temperature. Their relationships under elevated temperature are slightly unknown and thus must be parameterized separately at each level of elevated temperature (i) as follows:

Y=Yoi1+βiX    (3)

The parameter Yoi was further investigated based on the increasing air temperature level (i). A quadratic model was used to describe the relationship between the elevated temperature and the resultant rice grain yield (e.g., Sheehy et al., 2006a). Thus, equation 3 can be rewritten into equation 4 with the quadratic model.

Y=a+bGDD+cGDD21+βiX    (4)

where a, b, and c are unknown parameters for weed-free rice grain yield, and GDD is the accumulated growing degree days under elevated temperature. The parameter βi may also change with an increasing air temperature level (i). A linear model was used to describe the relationship between elevated temperature and resultant weed competitiveness (e.g., Kim et al., 2006). Equation 4 was further modified into equation 5 with the linear model.

Y=a+bGDD+cGDD21+(l+mGDD)X    (5)

where l and m are unknown parameters for weed competitiveness. Equations 4 and 5 were tested to determine whether the quadratic and linear models were appropriate for the change in weed-free grain yield and weed competitiveness, respectively, with GDD resulting from elevated temperature.

Statistical Analyses

All of the data were initially subjected to three-way analysis of variance (ANOVA). For weed biomass, rice yield components, and grain yield, data were analyzed with year, elevated temperature, and weed density as main factors and the replicate as a random factor. Correlation analyses were conducted for their relationships among yield components, elevated temperature, and weed density to evaluate the effects of elevated temperature and weed interference on the rice yield components and grain yield. Nonlinear regression analyses were also conducted to fit the rectangular hyperbolic model (equation 3) to the yield data corresponding to the plant density of each weed species for each air temperature condition. The quadratic and linear models were regressed on weed-free rice grain yield (Yo) and weed competitiveness (β) affected by elevated temperature in the pooled 3-year data, respectively. Then, by incorporating the models selected for the two parameters, Yo and β, into the rectangular hyperbolic model (equations 4 and 5), the combined models were evaluated, using the F-test (equation 6).

F=RSSi+1-RSSidfi+1-dfiRSSfdff    (6)

where RSS and df are the residual sum of the square and the degree of freedom, respectively; i + 1 is the model one-step reduced from its predecessor (i); and f is the full model. If the F-value was lower than the tabulated F-value (a 5% level) with (dfi+1-dfi, dff) degrees of freedom, then the reduced model was accepted. The performance of all the models was evaluated by the pseudo-R2 of the models and the root mean square error (RMS) of the prediction. All the statistical analyses were conducted, using Genstat (Genstat Committee, 2002).

Results

Effects of Elevated Temperature and Weed Interference on Rice

The three-way ANOVA revealed that both weed and rice were significantly affected by not only elevated temperature and weed density but also the year of experiment in a variety of ways (Supplementary Table 1). Late watergrass was affected solely by year and its plant density, while water chestnut was affected by each of the three main factors and the interaction of year and its plant density. In the case of rice, year, elevated temperature, and weed density interactively affected rice during the entire growing season. In particular, weed density more strongly affected rice vegetative growth and its yield component (e.g., number of panicles) determined at earlier growth stages, while the interactive effect of year and elevated temperature more strongly affected rice reproductive growth and its yield components (e.g., number of grains, % ripened grain, 1,000-grain weight) determined at later growth stages (Supplementary Table 1). As the ambient temperature differed in different years (Supplementary Figure 1), accumulated growing degree days (GDD) was expressed as a function of year and elevated temperature. Correlation analyses described the relationships between weed interference and rice yield components and between GDD and rice yield components across years (Figure 1). Dry weight and the number of panicles were negatively affected by weed density regardless of weed species, while the other yield components, including the number of grains, % ripened grains, and 1,000-grain weight, were not always affected by weed density (Figure 1). In contrast, the number of grains, % ripened grains, and 1,000-grain weight in the rice yield components were negatively affected by GDD, while the number of panicles was not always affected by GDD. Thus, rice grain yield was associated with weed interference and GDD, especially up to the flowering stage of rice.

FIGURE 1
www.frontiersin.org

Figure 1. Schematic representations of the correlations among growing degree days (transplanting to tillering, transplanting to flowering, and transplanting to maturity), plant densities of late watergrass (A) and water chestnut (B), rice yield components (dry weight, number of panicles, number of grains, % ripened grain, 1,000-grain weight), and rice grain yield over 3 years. Significance is indicated as follows: ns, not significant; *, P < 0.05.

Rice grain yield in 2016 was much lower than those in 2014 and 2015 (Figures 2, 3). This finding was attributed to the significant reduction in % ripened grains by high air temperatures over 30°C during the heading and flowering periods (Supplementary Figure 2). The % ripened grain was dependent on the ambient air temperature observed over 3 years, so it decreased with increasing ambient air temperature, giving 97.3, 96.0, and 88.8% ripened grains in 2014, 2015, and 2016, when the ambient (A) air temperature was 24.1, 26.4, and 29.0°C, respectively (Supplementary Figure 2). Elevated temperature further reduced the % ripened grain in combination with the ambient air temperature. In 2014, when the ambient air temperature was low, the % ripened grain remained high even at an elevated temperature of A+5.0°C, giving 93.4% ripened grains. In contrast, in 2016, when the ambient air temperature was high, the % ripened grain decreased significantly, resulting in 50.2, 4.9, and 1.7% ripened grains at elevated temperatures of A+1.5, A+3.0, and A+5.0°C, respectively. Due to the high ambient air temperature of 29.0°C in 2016, an elevated temperature of A+1.5°C resulted in over 30.5°C exposure during the heading and flowering periods (Table 1). This finding thus clearly demonstrates that the % ripened grain depends on the air temperature resulting from the combined ambient + elevated temperature that rice was exposed to during the heading and flowering periods, leading to a significant reduction of the % ripened grain if the combined ambient + elevated temperature is over 30°C.

FIGURE 2
www.frontiersin.org

Figure 2. Rice grain yield as a function of the weed density of late watergrass during accumulated growing degree days (GDD; from transplanting to flowering, 89 days) at four sunlit temperature-controlled phytotrons adjusted to ambient (A), A+1.5°C, A+3.0°C, and A+5.0°C in 2014 (•), 2015 (□), and 2016 (▾). The lines are fitted values calculated, using the rectangular hyperbolic model (equation 3) and parameter estimates (Table 2).

FIGURE 3
www.frontiersin.org

Figure 3. Rice grain yield as a function of the weed density of water chestnut during accumulated growing degree days (GDD; from transplanting to flowering, 89 days) at four sunlit temperature-controlled phytotrons adjusted to ambient (A), A+1.5°C, A+3.0°C, and A+5.0°C in 2014 (•), 2015 (□), and 2016 (▾). The lines are fitted values calculated, using the rectangular hyperbolic model (equation 3) and parameter estimates (Table 2).

Response of Weed-Free Rice Grain Yield and Weed Competitiveness to Elevated Temperature

Rice grain yield was regressed on single-species weed density at different GDDs, using a rectangular hyperbolic model (Table 2). The yield data for each species showed year-to-year variations due to the different GDDs (from transplanting to flowering, 89 days) at four sunlit temperature-controlled phytotrons adjusted to A, A+1.5°C, A+3.0°C, and A+5.0°C in 2014, 2015, and 2016 (Figures 2, 3). Quadratic and linear models were regressed on the weed-free rice grain yield and weed competitiveness affected by GDD (from transplanting to flowering, 89 days) over 3 years (Figures 4, 5). For the weed-free rice grain yield, the pooled value increased by up to 5.56 t and 5.72 t ha−1 for late watergrass and water chestnut, respectively, and, thereafter, decreased (Figure 4). Weed-free grain yield was associated with a lower percentage of ripened grains due to the elevated temperature during its heading and flowering periods (Supplementary Figure 2). For weed competitiveness, the pooled value exhibited a linear increase regardless of weed species with increasing GDD (Figure 5). Weed competitiveness was associated with a reduced number of panicles as a result of rice-weed competition during the period from the transplanting to tillering stages of rice. The weed competitiveness for the number of panicles showed a similar increase to the competitiveness for rice grain yield (Supplementary Figure 3).

TABLE 2
www.frontiersin.org

Table 2. Parameter estimates for the rectangular hyperbolic model of rice grain yield as a result of single-weed interference caused by late watergrass and water chestnut at accumulated growing degree days (from transplanting to flowering, 89 days) in 2014, 2015, and 2016.

FIGURE 4
www.frontiersin.org

Figure 4. Weed-free rice grain yield (Yo) as a function of the accumulated growing degree days (from transplanting to flowering, 89 days) at four sunlit temperature-controlled phytotrons adjusted to ambient, +1.5°C, +3.0°C, and +5.0°C in 2014, 2015, and 2016. The lines are fitted values calculated, using the quadratic model and parameter estimates (Supplementary Table 2) in the pooled 3-year data of late watergrass (A) and water chestnut (B).

FIGURE 5
www.frontiersin.org

Figure 5. Weed competitiveness (β) as a function of the accumulated growing degree days (from transplanting to flowering, 89 days) at four sunlit temperature-controlled phytotrons adjusted to ambient, +1.5°C, +3.0°C, and +5.0°C in 2014, 2015, and 2016. The lines are fitted values calculated, using the linear model and parameter estimates (Supplementary Table 3) in the pooled 3-year data of late watergrass (A) and water chestnut (B).

Combined Model for How Rice Grain Yield Is Affected by Elevated Temperature and Weed Interference

The quadratic and linear models were sequentially incorporated into the rectangular hyperbolic model to describe the relationships between weed-free rice grain yield and GDD and weed competitiveness and GDD, respectively. When the relationship between weed-free rice grain yield (Yoi) versus GDD was replaced by the quadratic model, there was a significant difference between equation 3 (the rectangular hyperbolic model) and equation 4 (the intermediate model) to describe the rice grain yield affected by late watergrass and water chestnut under elevated temperature (Table 3). When the relationship between weed competitiveness (βi) and GDD was replaced by the linear model, there was a significant difference between equation 4 (the intermediate model) and equation 5 (the combined rectangular hyperbolic model) in describing the same rice grain yield (Table 3). Although there was a significant difference between the tested models, the quadratic and linear models seemed to well describe the weed-free rice grain yield and weed competitiveness affected by GDD, respectively, with a small root mean square error and high pseudo-R2 values (Supplementary Tables 2, 3). Thus, these models can be incorporated in the rectangular hyperbolic model to describe the effects of elevated temperature on weed-free rice grain yield and weed competitiveness. Our results revealed that the combined rectangular hyperbolic model could describe both the elevated temperature and weed competition effects on rice grain yield regardless of weed species.

TABLE 3
www.frontiersin.org

Table 3. A summary of the nonlinear regression analysis and lack-of-fit test to compare the models for the regression of rice grain yield as a result of single-weed interference caused by late watergrass and water chestnut at accumulated growing degree days (from transplanting to flowering, 89 days) in the pooled 3-year data.

Using equation 5 (the combined rectangular hyperbolic model) and the estimated parameters (Table 4), the rice grain yield was predicted (Figure 6). The rice grain yields without interference were estimated to be 5.1 to 5.7 t, 5.6 to 5.8 t, 4.8 to 4.9 t, and 2.5 to 3.1 t ha−1 at 1,300, 1,500, 1,700, and 1,900°C days of GDD (from transplanting to flowering, 89 days), respectively. When late watergrass was interfered with 80 plants m−2, the model predicted that the rice grain yields would be reduced by 58.9, 48.5, 41.3, and 35.9% of the yields under weed-free conditions at 1,300, 1,500, 1,700, and 1,900°C days of GDD (from transplanting to flowering, 89 days), respectively. When water chestnut was applied at 80 plants m−2, the rice grain yields were reduced by 86.8, 64.3, 51.1, and 42.3% of the yields under weed-free conditions at 1,300, 1,500, 1,700, and 1,900°C days of GDD (from transplanting to flowering, 89 days), respectively. Our results demonstrated that rice grain yield was more reduced by weed interference under elevated temperature conditions under which weed competitiveness increases, but the resultant growth of rice decreases more than under ambient conditions. Therefore, proper weed management needs to be performed to prevent rice yield loss from both weed interference and elevated temperatures.

TABLE 4
www.frontiersin.org

Table 4. Parameter estimates for the combined rectangular hyperbolic model of rice grain yield as a result of single-weed interference caused by late watergrass and water chestnut at accumulated growing degree days (from transplanting to flowering, 89 days) in the pooled 3-year data.

FIGURE 6
www.frontiersin.org

Figure 6. Predicted rice grain yield as a function of the weed density and accumulated growing degree days (GDD; from transplanting to flowering, 89 days) under late watergrass (A) and water chestnut (B) interference and elevated temperature conditions. The mesh is fitted values calculated, using the combined rectangular hyperbolic model (equation 5) and parameter estimates (Table 4).

Discussion

Effects of Elevated Temperature on Rice-Weed Competition

Our results clearly demonstrated that year, elevated temperature, and weed interference interactively affected rice during the entire growing season, but each effect on rice yield components was different. Weed interference more strongly reduced rice vegetative growth and its yield component (e.g., number of panicles) determined at earlier growth stages, while the interactive effect of year and elevated temperature more strongly reduced rice reproductive growth and its yield components (e.g., number of grains, % ripened grain, 1,000-grain weight) determined at later growth stages (Supplementary Table 1; Figure 1). Previous studies reported that the dry weight and number of panicles in rice are determined earlier than other yield components because early rice is the most susceptible to weed interference (Moon et al., 2010, 2014). Other yield components related to grains, such as the number of grains, % ripened grain, and 1,000-grain weight, are the most susceptible to elevated temperature during the late growth stage of rice, particularly from the heading to flowering stages (Jung et al., 2015). Our study also revealed that rice grain yield was reduced more significantly as a result of decreased spikelet fertility and 1,000-grain weight due to the elevated temperature, resulting in a greater than 30°C mean air temperature during the heading and flowering periods. Plotting the % ripened grain against the mean air temperature during its heading and flowering periods revealed a significant decrease in the % ripened grain from the point when the mean air temperature increased over 30°C, resulting in a logistic decrease with increasing mean air temperature (Supplementary Figure 2). Previous studies also reported that high air temperatures over 30°C due to elevated temperature affect spikelet fertility and grain filling, resulting in a reduced % ripened grain and 1,000-grain weight, respectively (Kim et al., 2011; Nguyen et al., 2014). Therefore, plotting the observed weed-free rice grain yield against GDD resulting from elevated temperature showed continuous grain yield responses to GDD, and this relationship can be described, using the quadratic model, as shown in Figure 4.

Weed competitiveness was enhanced by elevated temperature, but its effects on rice grain yield were mainly due to a decreased number of panicles (Supplementary Figure 3). Late watergrass was more competitive than water chestnut at ambient temperature, but the response of weed competitiveness to elevated temperature was similar (Figure 5). Late watergrass has greater competitiveness against rice than water chestnut, as it intercepts more photosynthetically active radiation due to its high canopy height (Bouhache and Bayer, 1993; Lindquist and Kropff, 1996). In contrast, water chestnut has shown variable competitiveness under different environmental conditions, depending on the soil type and its relative seedling emergence time (Kobayashi and Ueki, 1983; Chun and Shin, 1994). With respect to elevated temperature, weed competitiveness was not significantly different between late watergrass and water chestnut, showing the same linear increase in their competitiveness against rice with GDD (Figure 5, Supplementary Table 3). This finding thus suggests that the increase in weed competitiveness with GDD can be described, using the linear model, and can be incorporated into the model describing rice-weed competition on rice grain yield.

Model Application for Future Weed Management

The combined model developed here can provide predictive information for rice grain yield caused by both elevated temperature and weed interference under future climate conditions. Our results indicated that the weed-free rice grain yield would increase by up to 5.8 t ha−1 at a GDD of 1,500 (°C days) and, thereafter, decrease with increasing GDD (Figure 6). The decrease in the weed-free rice grain yield was mainly attributed to high temperature-induced spikelet sterility during its heading and flowering periods (Table 1, Supplementary Figure 2). Similar results were previously reported for spikelet sterility of rice (Nguyen et al., 2014; Jung et al., 2015). The quadratic function of our combined model accurately described the impact of elevated temperature on the weed-free rice grain yield in the pooled 3-year data. The quadratic model represents the temperature effect on the weed-free yield of crops, including rice (Sheehy et al., 2006a,b). Our combined rectangular hyperbolic model based on GDD may propose a potential period for safe rice cultivation in consideration of weed interference under future elevated temperature conditions. However, the combined rectangular hyperbolic model needs to be further improved to describe the complex effects of elevated temperature on the spikelet sterility of rice during its heading and flowering periods. With respect to weed competitiveness, our results indicated that late watergrass and water chestnut can increase their competitiveness with an increase in the daily mean air temperature under future elevated temperature conditions. A linear function of our combined rectangular hyperbolic model provided an accurate description of the weed competitiveness response to elevated temperature. A study on weed competitiveness reported that the linear model accurately explained the response of weed competitiveness to other environmental factors, such as nitrogen fertilizer, in wheat (Kim et al., 2006). Weed competitiveness in our study increased less with increasing elevated temperature than the competitiveness affected by nitrogen fertilizer (e.g., Kim et al., 2006; Song et al., 2021). Although the elevated temperature had fewer effects on both crop growth and weed biomass than nitrogen fertilizer, the temperature effects on the interference relationships were quite considerable. Water chestnut was a less competitive weed species than late watergrass under ambient conditions but increased its competitiveness over rice to a value similar to the competitiveness for late watergrass under elevated temperature. Water chestnut can outcompete rice through root inhibition due to various factors, such as the dispersion of tuber propagation, earlier seedling emergence, and faster growth development under elevated temperature, although it is not a C4 weed species (Kim et al., 2010). In contrast, the same Echinochloa species as late watergrass showed a slight increase in aboveground biomass compared with the biomass of rice under elevated temperature (Alberto et al., 1996). These results indicated that the interference relationships under elevated temperature are dependent on the interactions among crop species, weed species, and air temperature during the growing seasons. Further validation of the combined rectangular hyperbolic model might be required to describe the complex effects of elevated temperature on rice-weed competition, using a more independent data set under various locations and years. Although the combined rectangular hyperbolic model and estimated parameters in this study need to be further validated, our combined model aids decision-making for timely weed control at the early rice growth stage under future elevated temperature conditions.

Conclusion

The combined rectangular hyperbolic model developed by incorporating the quadratic model for weed-free rice grain yield and the linear model for weed competitiveness under elevated temperature conditions accurately describes both the effects of elevated temperature and weed interference on rice grain yield over 3 years. Our combined model can support decision-making for late watergrass and water chestnut management in rice production under future elevated temperature conditions. Our results demonstrate that rice grain yield is significantly reduced by weed interference under elevated temperature. Therefore, proper weed management based on a prediction model is required to prevent rice yield loss from both weed interference and elevated temperature. For example, rice should not only be cultivated to avoid spikelet sterility under future elevated temperature conditions but should also be kept weed free by timely weed control at its early growth stages. Further studies should be conducted to validate the parameters of the combined rectangular hyperbolic model for spikelet sterility of rice during its heading and flowering periods and enhanced weed competitiveness during the vegetative stage of rice as a result of elevated temperature. Although further validation in a broad range of years and locations is required for the combined rectangular hyperbolic model and estimated parameters, our study can provide empirical information for describing both the effects of elevated temperature and weed interference on crop yield and give decision support for weed management under a future climate.

Data Availability Statement

The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.

Nomenclature

Late watergrass, Echinochloa oryzicola Vasing.; water chestnut, Eleocharis kuroguwai Ohwi; ambient temperature, A; and accumulated growing degree days, GDD

Author Contributions

J-SS, J-HI, and Y-HP were involved in conducting the experiments, gathering data, performing statistical analyses, interpretation of data, developing the figures and tables, and drafting the manuscript. D-SK was involved in designing the experiments and supervising the research as well as the overall revision of the manuscript. S-HL, M-JY, J-WK, and B-WL were involved in the revision of the manuscript. All the authors approved the final 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-SS was supported by the R&D Program of “Plasma Advanced Technology for Agriculture and Food (Plasma Farming, Project No. 1711124797)” through the Korea Institute of Fusion Energy funded by Government funds, Republic of Korea.

Conflict of Interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Acknowledgments

We thank Dr. Yean-Uk Kim and Dr. Beom-Seok Seo for their technical support.

Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fpls.2021.663779/full#supplementary-material

References

Alberto, A. M. P., Ziska, L. H., Cervancia, C. R., and Manalo, P. A. (1996). The influence of increasing carbon dioxide and temperature on competitive interactions between a C3 crop, rice (Oryza sativa) and a C4 weed (Echinochloa glabrescens). Funct. Plant Biol. 23, 795–802. doi: 10.1071/PP9960795

CrossRef Full Text | Google Scholar

Bachelet, D., and Gay, C. A. (1993). The impacts of climate change on rice yield: a comparison of four model performances. Ecol. Modell. 65, 71–93. doi: 10.1016/0304-3800(93)90127-E

CrossRef Full Text | Google Scholar

Bir, M. S. H., Eom, M. Y., Uddin, M. R., Park, T. S., Kang, H. W., Kim, D. S., et al. (2014). Weed population dynamics under climatic change. Weed Turf. Sci. 3, 174–182. doi: 10.5660/WTS.2014.3.3.174

CrossRef Full Text | Google Scholar

Bouhache, M., and Bayer, D. E. (1993). Photosynthetic response of flooded rice (Oryza sativa) and three Echinochloa species to changes in environmental factors. Weed Sci. 41, 611–614. doi: 10.1017/S0043174500076402

CrossRef Full Text | Google Scholar

Chun, J. C., and Shin, H. S. (1994). Effect of depth of tuber burial, soil temperature, and soil moisture on tuber sprouting of Eleocharis kuroguwai Ohwi. Kor. J. Weed Sci. 14, 49–55.

Google Scholar

Cousens, R. (1985). A simple model relating yield loss to weed density. Ann. Appl. Biol. 107, 239–252. doi: 10.1111/j.1744-7348.1985.tb01567.x

CrossRef Full Text | Google Scholar

Davis, A. S., and Ainsworth, E. A. (2012). Weed interference with field-grown soyabean decreases under elevated [CO2] in a FACE experiment. Weed Res. 52, 277–285. doi: 10.1111/j.1365-3180.2012.00915.x

CrossRef Full Text | Google Scholar

Deryng, D., Conway, D., Ramankutty, N., Price, J., and Warren, R. (2014). Global crop yield response to extreme heat stress under multiple climate change futures. Environ. Res. Lett. 9, 034011. doi: 10.1088/1748-9326/9/3/034011

CrossRef Full Text | Google Scholar

Fischer, A. J., Ateh, C. M., Bayer, D. E., and Hill, J. E. (2000). Herbicide-resistant Echinochloa oryzoides and E. phyllopogon in California Oryza sativa fields. Weed Sci. 48, 225–230. doi: 10.1614/0043-1745(2000)048(0225:HREOAE)2.0.CO;2

CrossRef Full Text | Google Scholar

Flint, E. P., and Patterson, D. T. (1983). Interference and temperature effects on growth in soybean (Glycine max) and associated C3 and C4 weeds. Weed Sci. 31, 193–199. doi: 10.1017/S004317450006882X

CrossRef Full Text | Google Scholar

Flint, E. P., Patterson, D. T., and Beyers, J. L. (1983). Interference and temperature effects on growth in cotton (Gossypium hirsutum), spurred anoda (Anoda cristata), and velvetleaf (Abutilon theophrasti). Weed Sci. 31, 892–898. doi: 10.1017/S0043174500070946

CrossRef Full Text | Google Scholar

Genstat Committee (2002) Reference Manual (Genstat 5 Release 6.1). Oxford, UK: VSN International.

Intergovernmental Panel on Climate Change (IPCC) (2013). “Climate change 2013: the physical science basis,” in Contribution of Working Group I to the Fifth Assessment Report of the Intergovernmental Panel on Climate Change, eds T. F. Stocker, D. Qin, G.-K. Plattner, M. Tignor, S. K. Allen, J. Boschung, et al. (Cambridge: Cambridge University Press), 1535.

Google Scholar

Jung, W. S., Lee, K. J., and Lee, B. W. (2015). Responses of spikelet fertility to air, spikelet, and panicle temperatures and vapor pressure deficit in rice. J. Crop Sci. Biotech. 18, 209–218. doi: 10.1007/s12892-015-0116-7

CrossRef Full Text | Google Scholar

Kennedy, R. A., Barrett, S. C. H., Zee, D. V., and Rumpho, M. E. (1980). Germination and seedling growth under anaerobic conditions in Echinochloa crus-galli(barnyard grass). Plant Cell Environ. 3, 243–248. doi: 10.1111/j.1365-3040.1980.tb00796.x

CrossRef Full Text | Google Scholar

Kim, D. S., Brain, P., Marshall, E. J. P., and Caseley, J. C. (2002). Modeling herbicide dose and weed density effects on crop-weed competition. Weed Res. 42, 1–13. doi: 10.1046/j.1365-3180.2002.00253.x

CrossRef Full Text | Google Scholar

Kim, D. S., Marshall, E. J. P., Brain, P., and Caseley, J. C. (2006). Modelling the effects of sub-lethal doses of herbicide and nitrogen fertilizer on crop-weed competition. Weed Res. 46, 492–502. doi: 10.1111/j.1365-3180.2006.00543.x

CrossRef Full Text | Google Scholar

Kim, H. Y., Horie, T., Nakagawa, H., and Wada, K. (1996). Effects of elevated CO2 concentration and high temperature on growth and yield of rice. Jpn. J. Crop Sci. 65, 644–651. doi: 10.1626/jcs.65.644

PubMed Abstract | CrossRef Full Text | Google Scholar

Kim, J., Shon, J., Lee, C. K., Yang, W., Yoon, Y., Yang, W. H., et al. (2011). Relationship between grain filling duration and leaf senescence of temperate rice under high temperature. Field Crops Res. 122, 207–213. doi: 10.1016/j.fcr.2011.03.014

CrossRef Full Text | Google Scholar

Kim, J. W., Moon, B. C., Lim, S. H., Chung, J. H., and Kim, D. S. (2010). Prediction of seedling emergence and early growth of Eleocharis kuroguwai Ohwi under elevated temperature. Kor. J. Weed Sci. 30, 94–102. doi: 10.5660/KJWS.2010.30.2.094

CrossRef Full Text | Google Scholar

Kim, Y. U., and Lee, B. W. (2019). Differential mechanisms of potato yield loss induced by high day and night temperatures during tuber initiation and bulking: photosynthesis and tuber growth. Front. Plant Sci. 10. doi: 10.3389/fpls.2019.00300

PubMed Abstract | CrossRef Full Text | Google Scholar

Kim, Y. U., Seo, B. S., Choi, D. H., Ban, H. Y., and Lee, B. W. (2017). Impact of high temperatures on the marketable tuber yield and related traits of potato. Eur. J. Agron. 89, 46–52. doi: 10.1016/j.eja.2017.06.005

CrossRef Full Text | Google Scholar

Kobayashi, H., and Ueki, K. (1983). Phenotypic variation and adaptation in Eleocharis kuroguwai Ohwi, a paddy perennial cyperaceous weed. J. Weed Sci. Technol. 28, 179–186. doi: 10.3719/weed.28.179

CrossRef Full Text | Google Scholar

Koo, B. H., Yoo, S. C., Park, J. W., Kwon, C. T., Lee, B. D., An, G., et al. (2013). Natural variation in OsPRR37 regulates heading date and contributes to rice cultivation at a wide range of latitudes. Mol. Plant 6, 1877–1888. doi: 10.1093/mp/sst088

PubMed Abstract | CrossRef Full Text | Google Scholar

Kraehmer, H., Jabran, K., Mennan, H., and Chauhan, B. S. (2016). Global distribution of rice weeds—a review. Crop Prot. 80, 73–86. doi: 10.1016/j.cropro.2015.10.027

CrossRef Full Text | Google Scholar

Lee, K. J., Kim, D. I., Choi, D. H., and Lee, B. W. (2015). Rice grain-filling characteristics under elevated air temperature in a temperate region. J. Crop Sci. Biotech. 18, 231–236. doi: 10.1007/s12892-015-0100-2

CrossRef Full Text | Google Scholar

Lindquist, J. L., and Kropff, M. J. (1996). Applications of an ecophysiological model for irrigated rice (Oryza sativa)-Echinochloa competition. Weed Sci. 44, 52–56. doi: 10.1017/S0043174500093541

CrossRef Full Text | Google Scholar

Mahajan, G., Singh, S., and Chauhan, B. S. (2012). Impact of climate change on weeds in the rice-wheat cropping system. Curr. Sci. 102, 1254–1255.

Google Scholar

Matsui, T., Namuco, O. S., Ziska, L. H., and Horie, T. (1997). Effects of high temperature and CO2 concentration on spikelet sterility in Indica rice. Field Crops Res. 51, 213–219. doi: 10.1016/S0378-4290(96)03451-X

CrossRef Full Text | Google Scholar

Matthews, R. B., Kropff, M. J., Horie, T., and Bachelet, D. (1997). Simulating the impact of climate change on rice production in Asia and evaluating options for adaptation. Agric. Syst. 54, 399–425. doi: 10.1016/S0308-521X(95)00060-I

CrossRef Full Text | Google Scholar

McMaster, G. S., and Wilhelm, W. W. (1997). Growing degree-days: one equation, two interpretations. Agric. For. Meteorol. 87, 291–300. doi: 10.1016/S0168-1923(97)00027-0

CrossRef Full Text | Google Scholar

Moon, B. C., Cho, S. H., Kwon, O. D., Lee, S. G., Lee, B. W., and Kim, D. S. (2010). Modelling rice competition with Echinochloa crus-galli and Eleocharis kuroguwai in transplanted rice cultivation. J. Crop Sci. Biotechnol. 13, 121–126. doi: 10.1007/s12892-010-0066-z

CrossRef Full Text | Google Scholar

Moon, B. C., Kim, J. W., Cho, S. H., Park, J. E., Song, J. S., and Kim, D. S. (2014). Modelling the effects of herbicide dose and weed density on rice-weed competition. Weed Res. 54, 484–491. doi: 10.1111/wre.12102

CrossRef Full Text | Google Scholar

Nah, G., Im, J. H., Kim, J. W., Park, H. R., Yook, M. J., Yang, T. J., et al. (2015). Uncovering the differential molecular basis of adaptive diversity in three Echinochloa leaf transcriptomes. PLoS ONE 10, e0134419. doi: 10.1371/journal.pone.0134419

PubMed Abstract | CrossRef Full Text | Google Scholar

Nguyen, D. N., Lee, K. J., Kim, D. I., Anh, N. T., and Lee, B. W. (2014). Modeling and validation of high-temperature induced spikelet sterility in rice. Field Crops Res. 156, 293–302. doi: 10.1016/j.fcr.2013.11.009

CrossRef Full Text | Google Scholar

Nguyen, H. T. T., Shim, I. S., Kobayashi, K., and Usui, K. (2005). Effects of salt stress on ion accumulation and antioxidative enzyme activities of Oryza sativa L. and Echinochloa oryzicola Vasing. Weed Biol. Manag. 5, 1–7. doi: 10.1111/j.1445-6664.2005.00150.x

CrossRef Full Text | Google Scholar

Park, M. W., Kim, J. W., Lim, S. H., Lee, I. Y., and Kim, D. S. (2010). Prediction of seedling emergence and early growth of Monochoria vaginalis and Scirpus juncoides under elevated temperature. Kor. J. Weed Sci. 30, 103–110. doi: 10.5660/KJWS.2010.30.2.103

CrossRef Full Text | Google Scholar

Patterson, D. T., Flint, E. P., and Beyers, J. L. (1984). Effects of CO2 enrichment on competition between a C4 weed and a C3 crop. Weed Sci. 32, 101–105. doi: 10.1017/S0043174500058586

CrossRef Full Text | Google Scholar

Ramesh, K., Matloob, A., Aslam, F., Florentine, S. K., and Chauhan, B. S. (2017). Weed in a changing climate: vulnerabilities, consequences, and implications for future weed management. Front. Plant Sci. 8, 95. doi: 10.3389/fpls.2017.00095

PubMed Abstract | CrossRef Full Text | Google Scholar

Rodenburg, J., Meinke, H., and Johnson, D. E. (2011). Challenges for weed management in African rice systems in a changing climate. J. Agric Sci. 149, 427–435. doi: 10.1017/S0021859611000207

CrossRef Full Text | Google Scholar

Rural Development Administration [RDA]. (2014). Nongsaro: Agricultural Technology Information. Available online at: http://www.nongsaro.go.kr. (accessed January 1, 2014).

Sánchez, B., Rasmussen, A., and Porter, J. R. (2014). Temperatures and the growth and development of maize and rice: a review. Glob. Change Biol. 20, 408–417. doi: 10.1111/gcb.12389

PubMed Abstract | CrossRef Full Text | Google Scholar

Sheehy, J. E., Mitchell, P. L., Allen, L. H., and Ferrer, A. B. (2006b). Mathematical consequences of using various empirical expressions of crop yield as a function of temperature. Field Crops Res. 98, 216–221. doi: 10.1016/j.fcr.2006.02.008

CrossRef Full Text | Google Scholar

Sheehy, J. E., Mitchell, P. L., and Ferrer, A. B. (2006a). Decline in rice grain yields with temperature: models and correlations can give different estimates. Field Crops Res. 98, 151–156. doi: 10.1016/j.fcr.2006.01.001

CrossRef Full Text | Google Scholar

Song, J. S., Im, J. H., Kim, J. W., Kim, D. G., Lim, Y., Yook, M. J., et al. (2021). Modeling the effects of nitrogen fertilizer and multiple weed interference on soybean yield. Agronomy11:515. doi: 10.3390/agronomy11030515

CrossRef Full Text | Google Scholar

Song, J. S., Kim, J. W., Im, J. H., Lee, K. J., Lee, B. W., and Kim, D. S. (2017b). The effects of single- and multiple-weed interference on soybean yield in the far-eastern region of Russia. Weed Sci. 65, 371–380. doi: 10.1017/wsc.2016.25

CrossRef Full Text | Google Scholar

Song, J. S., Lim, S. H., Yook, M. J., Kim, J. W., and Kim, D. S. (2017a). Cross-resistance of Echinochloa species to acetolactate synthase inhibitor herbicides. Weed Biol. Manag. 17, 91–102. doi: 10.1111/wbm.12123

CrossRef Full Text | Google Scholar

Song, J. S., Park, Y. S., Park, M. W., Lee, J. D., and Kim, D. S. (2016). Selectivity of tefuryltrione between rice and Eleocharis kuroguwai. Weed Turf. Sci. 5, 191–195. doi: 10.5660/WTS.2016.5.4.191

CrossRef Full Text | Google Scholar

Thompson, G. B., and Woodward, F. I. (1994). Some influences of CO2 enrichment, nitrogen nutrition and competition on grain yield and quality in spring wheat and barley. J. Exp. Bot. 45, 937–942. doi: 10.1093/jxb/45.7.937

CrossRef Full Text | Google Scholar

Ueno, O., and Takeda, T. (1992). Photosynthetic pathways, ecological characteristics, and the geographical distribution of the Cyperaceae in Japan. Oecologia 89, 195–203. doi: 10.1007/BF00317218

PubMed Abstract | CrossRef Full Text | Google Scholar

Watanabe, H. (2011). Development of lowland weed management and weed succession in Japan. Weed Biol. Manag. 11, 175–189. doi: 10.1111/j.1445-6664.2011.00419.x

CrossRef Full Text | Google Scholar

Welch, J. R., Vincent, J. R., Auffhammer, M., Moya, P. F., Dobermann, A., and Dawe, D. (2010). Rice yields in tropical/subtropical Asia exhibit large but opposing sensitivities to minimum and maximum temperatures. Proc. Natl. Acad. Sci. U.S.A. 107, 14562–14567. doi: 10.1073/pnas.1001222107

PubMed Abstract | CrossRef Full Text | Google Scholar

Yamasue, Y. (2001). Strategy of Echinochloa oryzicola Vasing. for survival in flooded rice. Weed Biol. Manag. 1, 28–36. doi: 10.1046/j.1445-6664.2001.00008.x

CrossRef Full Text | Google Scholar

Zeng, Q., Liu, B., Gilna, B., Zhang, Y., Zhu, C., Ma, H., et al. (2011). Elevated CO2 effects on nutrient competition between a C3 crop (Oryza sativa L.) and a C4 weed (Echinochloa crus-galli L.). Nutr. Cycl. Agroecosyst. 89, 93–104. doi: 10.1007/s10705-010-9379-z

CrossRef Full Text | Google Scholar

Ziska, L. H. (2000). The impact of elevated CO2 on yield loss from a C3 and C4 weed in field-grown soybean. Glob. Change Biol. 6, 899–905. doi: 10.1046/j.1365-2486.2000.00364.x

CrossRef Full Text | Google Scholar

Ziska, L. H. (2001). Changes in competitive ability between a C4 crop and a C3 weed with elevated carbon dioxide. Weed Sci. 49, 622–627. doi: 10.1614/0043-1745(2001)049(0622:CICABA)2.0.CO;2

CrossRef Full Text | Google Scholar

Ziska, L. H. (2003). Evaluation of yield loss in field sorghum from a C3 and C4 weed with increasing CO2. Weed Sci. 51, 914–918. doi: 10.1614/WS-03-002R

CrossRef Full Text | Google Scholar

Ziska, L. H., Tomecek, M. B., and Gealy, D. R. (2010). Competitive interactions between cultivated and red rice as a function of recent and projected increases in atmospheric carbon dioxide. Agron. J. 102, 118–123. doi: 10.2134/agronj2009.0205

CrossRef Full Text | Google Scholar

Keywords: crop-weed competition, elevated temperature, modeling, late watergrass, water chestnut, rice

Citation: Song J-S, Im J-H, Park Y-H, Lim S-H, Yook M-J, Lee B-W, Kim J-W and Kim D-S (2021) Modeling the Effects of Elevated Temperature and Weed Interference on Rice Grain Yield. Front. Plant Sci. 12:663779. doi: 10.3389/fpls.2021.663779

Received: 03 February 2021; Accepted: 16 June 2021;
Published: 16 July 2021.

Edited by:

Brian Farneti, Fondazione Edmund Mach, Italy

Reviewed by:

Doğan Işik, Erciyes University, Turkey
D.K. Lee, University of Illinois at Urbana-Champaign, United States

Copyright © 2021 Song, Im, Park, Lim, Yook, Lee, Kim and Kim. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Do-Soon Kim, dosoonkim@snu.ac.kr

ORCID: Jong-Seok Song orcid.org/0000-0002-1900-7806
Do-Soon Kim orcid.org/0000-0002-7388-4308

Disclaimer: All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.