Next Article in Journal
Advances in Barley Breeding for Improving Nitrogen Use Efficiency
Next Article in Special Issue
Effects of High Temperature on Quality of Japonica Rice at Early and Middle Heading Stage under Different Planting Modes
Previous Article in Journal
Spatial and Temporal Evolution of Sowing and the Onset of the Rainy Season in a Region of Large Agricultural Expansion in Brazil
Previous Article in Special Issue
Exogenous Application of Methyl Jasmonate at the Booting Stage Improves Rice’s Heat Tolerance by Enhancing Antioxidant and Photosynthetic Activities
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Effects of Precise K Fertilizer Application on the Yield and Quality of Rice under the Mode of Light, Simple, and High-Efficiency N Fertilizer Application during the Panicle Stage

Rice Research Institute, Agronomy College, Shenyang Agricultural University, Shenyang 110866, China
*
Author to whom correspondence should be addressed.
Submission received: 8 June 2022 / Revised: 11 July 2022 / Accepted: 14 July 2022 / Published: 15 July 2022
(This article belongs to the Special Issue Advances in Rice Physioecology and Sustainable Cultivation)

Abstract

:
Light, simple, and high-efficiency fertilization is currently an effective method used to improve quality and increase yield. Most research has only focused on the yield or quality of rice, and no effective or in-depth studies exist on the key stage of panicle, which is essential for these two aspects. This study aimed to determine the effect of applying Nitrogen (N) and Potassium (K) fertilizers on the yield and quality at different leaf ages. The rice cultivar, Shennong 265, was grown in the field-tube condition at the 5-leaf age with K fertilizer at different panicle stages. Five K fertilizer and four N fertilizer levels were applied simultaneously during two growing seasons in 2020–2021. The application of K fertilizer at different panicle stages significantly affected the yield and quality. The application of K and N fertilizers at the 12th-leaf-age increased the number of panicles and grains per panicle, thereby increasing the yield with an average increase of 23.20% over local farmer’s fertilization model (CK) in two years. Application at the 10th-leaf age reduced the starch and protein content of the rice and improved the nutritional quality and taste, with an average increase of 11.08 points compared to CK in two years. The processing quality was the best at the panicle K fertilizer application rate of 47.81–64.69 kg ha−1, and the starch and protein contents were the lowest at the panicle K fertilizer application rate of 56.25 kg ha−1. Under different panicle K fertilizer application rates, N fertilizers had great differences in processing quality but had slight differences in nutritional quality; that is, the starch and protein content (total protein and four protein components) significantly increased. The application of panicle N fertilizer mainly affected starch pasting viscosity (RVA profile characteristics). When panicle K fertilizer was applied after panicle N fertilizer, the yield and quality showed a trend of synergistic improvement. Although this improvement was a low-level synergy, it can still be a direction used to explore the synergistic improvement of output and quality.

1. Introduction

Rice is one of the three most important cereal crops in the world. Asian countries account for 90% of the global rice production [1]. China, located in East Asia, is the world’s largest producer and consumer of rice. Within China, Northeast China is an important commercial grain base and the largest high-quality japonica rice-producer [2]. Therefore, it is very important to ensure the yield and quality of rice in Northeast China [3].
It is well known that nitrogen (N) plays a crucial role in improving rice yield. Potassium (K) is one of the essential macronutrients used for plant growth and development. It is involved in plant photosynthesis, cell osmotic regulation, enzyme activation, protein synthesis, grain quality, and yield in rice. The application of N and K at a reasonable ratio can effectively increase plant height, effective panicle number, and leaf area index; promote photosynthetic product synthesis and accumulation; increase rice yield; and improve rice quality [4,5,6,7]. Over the past 70 years, global rice production has gradually increased owing to increased N inputs [8]. However, with the rapid development of agriculture, N fertilizers are increasingly used in rice production to increase efficiency at a low cost. To ensure a high yield of rice, most rice farmers apply more N fertilizers than necessary. This coupled with runoff, denitrification, and leaching has caused a large amount of N loss, resulting in low N recovery efficiency and serious environmental pollution [9]. Moreover, because of the problem of resource distribution in China, K resources are relatively in short supply and are imported in large quantities. Thus, to save cost, K fertilizer application is neglected in agricultural production, resulting in an unbalanced and unscientific application of N and K fertilizers [10]. In recent years, with the development of China’s economy and the improvement of the living standards of the people, demand for rice has gradually changed from “eat enough” to “eat well”. Therefore, it is necessary to improve the quality of rice while stabilizing the output [11].
The number of various fertilizers required for rice is as follows: 2.0–2.4 kg of N, 0.9–1.4 kg of P2O5, and 2.5–2.9 kg of K2O per 100 kg of rice. The N uptake by rice is relatively average in the early, middle, and late stages, but it is the most significant in the booting and grain filling stages [12,13]. However, the absorption of K in rice is low in the early stages, 35% from the tillering to booting stages and 65% from the booting stage to maturity [14]. Therefore, to establish a balance between the output and quality of rice, more attention should be paid to the key stage affecting the yield and quality of rice—the fertilization management during the panicle stage. A large number of studies have focused on the effects of the total amount of N fertilization, the ratio of N fertilization in each period, the frequency and time of N fertilization on rice population indicators, the N use efficiency (NUE), and yield [15,16,17,18,19]. The results showed that increasing the proportion of nitrogen fertilizer application at a late growth stage could significantly optimize the population quality, N utilization, and yield of rice. However, these studies only focused on a single fertilizer and not on a more detailed and accurate ratio of N and K fertilizers during the panicle stage and the precise application time for the certain key period of rice growth and development and ignored the related effects on rice quality. Ample research has been conducted on the simplified application of N fertilizer in the early stage [20,21]. The results showed that a ratio of 6:4 for base + tiller N fertilizer to panicle N fertilizer has the best results. When the ratio of panicle N fertilizer is increased, the tiller N application period is postponed (the leaf age index is 60%), and the panicle N fertilizer application period is preponed (the leaf age index is 80%), which can significantly improve N recovery efficiency (NRE) and N physiological efficiency (NPE) and lead to high yield and efficiency [22,23]. In previous studies, increasing rice yield inevitably reduced rice quality; yield and quality have never been balanced. In this study, the tube planting method, which better simulates the environmental conditions of the field than the pot planting method, was used for two consecutive years of traditional cultivation methods to precisely control the time and amount of nitrogen and potassium fertilizer application. We focused on how to improve the quality traits of rice as much as possible without reducing the yield and to improve the efficient utilization of fertilizers without increasing the input (fertilizer or other economical inputs), to balance yield and quality. The study will provide ideas for establishing a yield and quality evaluation system in high-quality japonica rice-producing areas in Northeast China.

2. Materials and Methods

2.1. Plant Materials and Site Description

Field experiments were conducted on a farm at Shenyang Agricultural University (41°49′ N, 123°34′ E; 41.6 m above sea level) during the rice-growing season (May to October) in 2020 and 2021. The site was located within the main high-yielding rice region of the Liaohe River basin and has a temperate, semi-humid, continental climate. The daily mean temperatures, precipitation, relative humidity, and total sunshine hours during the growing seasons in 2020 and 2021 are provided (Figure 1). The soil type of the field was brown loam. The total N, available P, available K, organic matter, and pH in 0–20 cm of the soil were 1.07 g kg−1, 35.3 mg kg−1, 143.2 mg kg−1, 25.1 g kg−1, and 7.05, respectively; soil fertility is moderate. The test variety was Shennong 265 (15 leaves on the main stem). The typical plant height is 100–105 cm, the plant type is compact, the leaves are wide and firm, the stems and leaves are dark green, the tillering ability is medium, and the growth period is 155 days. It was a typical erect panicle japonica rice cultivar. It has the characteristics of high yield, high NUE, and a reasonable canopy structure [24]. It is one of the most widely planted japonica rice varieties in the rice-growing region in central Liaoning Province, China [25,26].

2.2. Experimental Design

The experiment was designed in a completely randomized manner with three replications. K fertilizer was applied in five leaf-growth stages in the panicle stage: 10th-leaf, 11th-leaf, 12th-leaf, 13th-leaf, and 14th-leaf. Five K fertilizer levels and four N fertilizer levels were set. The specific fertilization amount and fertilization period are shown in Table 1. Urea with an N content of 46.4% was used as N fertilizer, superphosphate with 12.0% P2O5 as P fertilizer, and potassium chloride with 60.0% K2O as K fertilizer. There were 101 treatments in the experiment (5 different periods of K fertilization × 5 K fertilizer levels × 4 spikes of nitrogen fertilizer application amount, a total of 100 treatments, plus 1 CK), and each treatment was performed with one barrel, which was repeated three times, for a total of 303 barrels. Planting in the field was carried out using the tube planting method. After raking the ground evenly, a bottomless PVC cylinder with an inner diameter of 30 cm, a wall thickness of 0.5 cm, and a cylinder height of 50 cm was buried. The cylinder nails were pressed to the bottom of the plow (35 cm). The seeds were sown on 24 April 2020 and 21 April 2021. In the four-leaf heart stage, seedlings with similar growth were selected and transplanted on 23 May 2020 and 27 May 2021, and they were harvested on 8 October 2020 and 7 October 2021. There were two seedlings per barrel at a hill spacing of 15.5 cm × 31.5 cm. Except for the fertilization period, conventional management methods were adopted. Basal fertilizer was applied on 22 May 2020 and 26 May 2021; tillering fertilizer was applied on 23 June 2020 and 1 July 2021; panicle N fertilizer was applied on 22 July 2020 and 26 July 2021. The first K application in the panicle stage was on 12 July and then every 6 days or so (the application time was based on the specific leaf age), and the last application in the panicle stage was on 1 August 2020 and 4 August 2021.
As there were many factors in this experiment, the analysis was divided into two parts. In the first part, we only analyzed the difference between the application ratio of N and K fertilizers at different leaf stages during the panicle stage and CK (the same amount of fertilizer). The different stages were (10,12), (11,12), (12,12), (13,12), and (14,12), which indicate the application times of the panicle K and N fertilizers. That is, the panicle K fertilizer was applied on the 10th-leaf, 11th-leaf, 12th-leaf, 13th-leaf, or 14th-leaf, respectively, and only once, and the panicle N fertilizer was applied on the 12th-leaf. In the second part, we only analyzed the difference in the amount of N and K fertilizers during different panicle stages, regardless of the time of application of panicle fertilizer. Among them, there were three levels of N (N1, N2, and N3) and five levels of K (K1, K2, K3, K4, and K5). N0 was used only for the calculation of the regression equation and does not appear in the analysis as an N application level.

2.3. Sampling and Measurements

2.3.1. Determination of Basic Physical and Chemical Properties of Soil

Before the spring test in 2020, soil samples of 0–20 cm in the experimental field were collected using the “S”-shaped sampling method and transferred to the lab. After air drying the soil samples, they were ground and passed through 20- and 100-mesh sieves. The total N content was determined using the elemental analyzer (EA 3000). The content of available P was determined using the molybdenum antimony scandium colorimetric method [27]. The content of available K was determined using flame photometry [27]. Oxidation of the soil organic matter was assessed using the potassium dichromate-external heating method [27]. The soil pH was measured using an electronic pH meter (PHSJ-3F), and the soil-liquid ratio was 1:2.5 [27].

2.3.2. Grain Yield

All ears of rice were collected and planted in tubes, and each hole was individually marked and packaged. After air-drying, an indoor seed test was conducted, and the number of panicles, number of grains per panicle, seed setting rate, and 1000-grain weight were measured.

2.3.3. Determination of Milling Quality of Rice Grain

The seeds were threshed and naturally dried and then stored for 3 months for grain quality measurement. Milling quality, including brown rice rate, milled rice rate, and head milled rice rate were measured according to the GB/T 17891-1999 (1999) standard.

2.3.4. Determination of Amylose and Protein Contents of Rice Grain

The amylose and protein content of brown rice was measured using the Foss 1241 near-infrared grain analyzer (Infratec TM 1241, FOSS, Hillerød, Denmark). The protein component was determined according to Luthe and Liu [28,29], based on Bradford’s BCA protein quantification method.

2.3.5. Determination of Starch Pasting Viscosity of Rice Grain

Pasting properties of rice flour were measured using a Rapid Viscosity Analyzer (RVA-TecMaster, Perten, Sweden), according to the American Society for Cereal Chemistry Operating Procedures [30]. The RVA spectrum characteristic value was mainly represented by peak viscosity (PKV), hot paste viscosity (HPV), cool paste viscosity (CPV), breakdown viscosity (BDV = PKV − HPV), and setback viscosity (SBV = CPV − PKV) [31].

2.3.6. Determination of Taste Quality of Rice Grain

The taste quality of the rice grain was measured using the STA1A rice taste meter produced by SATAKE in Japan. An amount of 30 g of milled rice was weighed, rinsed until the water became clear, and then soaked in a stainless-steel tank for 30 min. Water was added in a ratio of 1:1.2 (rice:water); the rice was put in a steamer to cook for 40 min and cooled for 1 h 50 min, and a taste meter was used to measure the indicators of appearance, hardness, viscosity, balance degree, and taste value. The process was repeated thrice for each sample.

2.4. Statistical Analysis

For experimental variables, two-factor and three-way analysis of variance (ANOVA) was used to assess differences among treatments using the IBM SPSS Statistics 22.0 software (Softonic International, Barcelona, Spain). Significant differences (p < 0.05) between treatments were indicated by different letters according to Fisher’s LSD. Graphs were drawn with the Origin 2021 software (OriginLab, Northampton, MA, USA).

3. Results

3.1. Grain Yield and Yield Components

The grain yield first increased and then decreased with the increase of leaf age of K application at the panicle stage. Compared with those of CK, the K fertilizer applied to 12th-leaf and 13th-leaf at the panicle stage substantially increased the yield to 1302.76 g∙m−2 and 1135.91 g∙m−2 in 2020 and 1174.86 g∙m−2 and 1070.92 g∙m−2 in 2021, respectively, with increases of 8.08–23.92% in 2020 and 11.58–22.40% in 2021, respectively. The yield was increased by substantially increasing the number of panicles per plant. The seed-setting rate and 1000-grain weight increased with the increase of the leaf age of the K fertilizer application at the panicle stage and reached the maximum value at the 14th-leaf stage, which was the trend for both years (Figure 2). With the increase of the amount of N fertilizer applied at the panicle stage, the yield increased and then decreased, and it reached the highest yield at the application rate of K2 and K3, which was substantially higher than that of K1, K4, and K5, mainly because the spike number and grain number per spike were substantially increased. With the increase of the N fertilizer application rate at the panicle stage, the panicle number, grain number per panicle, and yield gradually increased under low K (K1 and K2) conditions as N3 > N2 > N1, and there was no substantial difference between the seed setting rate and 1000-grain weight. However, under the condition of medium and high K (K3, K4 and K5), a gradually decreasing trend was observed as N3 < N2 < N1. The interaction of year and N and K application at the panicle stage was not substantial (Table 2).

3.2. Processing Quality

With the increase in leaf age of the K fertilizer application rate at the panicle stage, the processing quality first decreased, then increased, and then decreased. The brown rice rate and milled rice rate reached the maximum at the 12th-leaf stage, while the head rice rate reached the maximum at the 13th-leaf stage. The trend of both years was consistent, and the interaction was not substantial (Figure 3). With the increase of the K fertilizer application rate at the panicle stage, the difference in the brown rice rate and milled rice rate of K1–K4 and K5 reached a substantial level, but the difference between K1 and K4 was not substantial. The head rice rates of K2, K3, and K4 were substantially higher than those of K1 and K5, indicating that too much or too little K fertilizer application would reduce the processing quality. With the increase of the N fertilizer application rate, at the K1 and K2 levels, the brown rice rate, the milled rice rate, and the head rice rate showed a gradually increasing trend, and at the K3, K4, and K5 levels, it showed a gradually decreasing trend. The interaction of N and K both reached a substantial level, and the interaction between N and K and year was not substantial. The interaction between the two factors of N and K reached a substantial level. The interaction of year and N and K application at the panicle stage was not substantial (Table 3).

3.3. Nutritional Quality

With the increase of the leaf age of the K fertilizer application at the panicle stage, the amylose content first decreased and then increased, and all the content, except the 10th-leaf stage, was substantially lower than that of CK; the protein content showed a gradual upward trend, and only 10th-leaf stage had a substantially lower protein content than that of CK. In terms of protein components, the effect of the K fertilizer application on different leaf ages on albumin and prolamin was not obvious, whereas the effect of globulin and glutenin increased first and then decreased, and they both reached the highest at the 13th-leaf stage. Compared with that of CK, albumin in the 11th–14th-leaf, globulin in the 13th-leaf, and glutenin in the 13th- and 14th-leaf were substantially higher than those in CK; albumin in the 10th- and 11th-leaf, prolamin in the 10th–14th-leaf, and the gluten of the 10th-leaf was substantially lower than that of CK. The two-year trend remained the same; the differences among treatments reached an extremely substantial level, and the inter-annual differences were not substantial, indicating that the application of spike K fertilizers at different leaf ages could change the amylose and protein content and could affect the differences between protein components (Figure 4). With the increase of the K fertilizer application rate at the panicle stage, amylose and protein contents first decreased and then increased, and they reached the lowest value at K3. The change in the K fertilizer application rate can cause a change in protein composition. With the increase in the K application rate, the contents of albumin, prolamin, and glutenin first decreased and then increased, while globulin first increased and then decreased. The four proteins reached the lowest under the K3 treatment and reached the highest under the K5 treatment. Amylose and protein content showed gradual upward trends with the increase of the N fertilizer application. The analysis of the protein components showed that, under the low K condition (K1), the effect of the N fertilizer on gluten did not reach a substantial level, but the effect on the other three proteins reached a substantial level. Under the condition of medium K (K3), except for albumin, the three proteins had substantial differences only under the condition of low N (N1) but not under the condition of medium and high N (N2, N3). Under medium-high K (K4) conditions, the N application rate had the least effect on protein components. Under high K (K5) conditions, high N (N3) substantially increased the contents of the four proteins and severely reduced nutritional quality. The differences of N and K levels reached extremely substantial levels, and the differences in protein components of N and K interactions reached substantial or extremely substantial levels. The two-year trend was consistent. The interaction of year and N and K application at the panicle stage was not substantial (Table 4).

3.4. RVA Profile Characteristics

With the increase of the leaf age of the K fertilizer application at the panicle stage, PKV, CPV, and SBV showed gradual upward trends; HPV and BDV showed a trend of increasing first, then decreasing, and then increasing. Compared with that of CK, the 10th-leaf application substantially reduced PKV, HPV, and CPV, the 11th-leaf and 12th-leaf applications substantially increased SBV, the 13th-leaf application substantially decreased SBV, and the 14th-leaf application substantially increased CPV and SBV. In 2020, all the differences except PKV reached a substantial level, and in 2021, except for PKV and BDV, the differences reached a substantial level, and the inter-annual differences were not substantial (Figure 5). The two-year trends of the RVA profile differed slightly. The effect of the N application and K application on the RVA profile was quite different, and the effect of the K application on the RVA profile was small. Only the difference between BDV and SBV in 2020 reached a substantial level, and the others were not substantial. The effect of the N application on the RVA profile value was greater, all reaching the substantial or extremely substantial level. With the increase of the N application rate at the panicle stage, PKV, HPV, CPV, and BDV gradually decreased, while SBV increased gradually. The three-factor interaction difference between N and K and annual BDV and SBV reached a substantial level (Table 5).

3.5. Cooking Quality

With the increase of the leaf age of K fertilizer application at the panicle stage, the taste value gradually decreased. Compared with that of CK, 10th-leaf substantially improved the taste value by increasing the appearance, viscosity, and balance degree and by decreasing hardness. The 14th-leaf substantially reduced the taste value by reducing viscosity and balance degree and increasing hardness. Other leaf ages were not substantially different. The differences between treatments of each index reached extremely substantial levels; the inter-annual differences were not substantial; the differences between treatments and inter-annual interactions were not substantial (Figure 6). With the increase of the K fertilizer application at the panicle stage, the measure of the appearance, viscosity, and balance degree first increased and then decreased, whereas the hardness first decreased and then increased, thus affecting the taste value. The trend of taste value first increased and then decreased, reaching the maximum value of 49.21 points in 2020 and 50.30 points in 2021 at K2. It affected the taste value by substantially increasing appearance, viscosity, and balance degree. The trend was consistent for two years, but the difference in the effect on hardness was inconsistent for the two years, which needs to be further explored. With the increase of the N fertilizer application at the panicle stage, the taste value showed a downward trend. Under medium and high K conditions (K3, K4, K5), appearance, viscosity, and balance degree were substantially reduced and hardness increased substantially, which in turn affected the taste value. Under the condition of low K (K1), the effect of N fertilizer was lower, the difference in taste value was not substantial, and the two-year rule was generally consistent. The interaction of N and K factors reached an extremely substantial level, and the interaction of year and N and K application at the panicle stage and year was not substantial (Table 6).

3.6. Correlations between the Traits

Figure 7 shows the correlation between grain yield and K application at different leaf ages at the panicle stage. The two-year average value is related to the output in a univariate nonlinear manner, and the regression equation is y = −0.26 x2 + 29.26 x + 211.60 (R2 = 0.62, p < 0.05). This shows that the yield increased first and then decreased with the increase of the K application at the panicle stage. When the K application at the panicle stage amount reached 56.27 kg ha−1, the yield reached the maximum, which was 1034.82 g∙m−2. The trend was different among leaf ages.
Figure 8 shows the correlation between grain yield and N application at different leaf ages at the panicle stage. The two-year average value is related to the output in a univariate nonlinear manner, and the regression equation is y = −0.0012 x2 + 0.22 x + 38.86 (R2 = 0.95, p < 0.05). This shows that the yield increased first and then decreased with the increase of N application at the panicle stage. When the N application at the panicle stage amount reached 93.96 kg ha−1, the yield reached the maximum, which was 1008.55 g∙m−2. The trend was the same among leaf ages.
Figure 9 shows the correlation between taste value and K application at different leaf ages at the panicle stage. The two-year average value is related to the output in a univariate nonlinear manner, and the regression equation is y = −0.0023 x2 + 0.18 x + 46.83 (R2 = 0.93, p < 0.05). This shows that the yield increased first and then decreased with the increase of K application at the panicle stage. When the K application at the panicle stage amount reached 39.13 kg ha−1, the yield reached the maximum, which was 52.97. The trend was the same among leaf ages.
The correlation between taste value and the N application rate at different leaf ages at the panicle stage is shown in Figure 10. The two-year average is linearly correlated with taste value, and the regression equation is y= −0.088 x + 56.84, (R2 = 0.96, p < 0.05). This indicated that there was a substantial negative correlation between the N application at the panicle stage and taste value. The trend was the same among leaf ages.
The correlation analysis of the yield and quality indicators of different leaf ages applying K fertilizer during the panicle stage (Figure 11) reveals that the premature application of K fertilizer during the panicle stage (10th, 11th-leaf) had a very substantial positive correlation between the number of panicles per plant and the number of grains per panicle and yield; the PKV, HPV, BDV, appearance, viscosity, and balance degree were positively correlated with taste value; 1000-grain weight, albumin content, globulin content, SBV, and hardness were substantially negatively correlated with taste value (a). When K fertilizer was applied at the 12th-leaf age, its yield-related effects were the same as those at the 10th- and 11th-leaf ages. However, its effect on the taste value was different from that of the 10th- and 11th-leaf. The influence of amylose, viscosity, and balance degree on the taste value was extremely substantially positive; the protein content, globulin content, and SBV were extremely substantially negatively correlated with taste value (b). The correlation between the indicators changed during the panicle stage (13th and 14th leaves) when the K fertilizer was applied too late. The number of panicles per plant, the number of grains per panicle, brown rice rate, milled rice rate, globulin content, and viscosity were substantially positively correlated with yield; total protein content and prolamin content were extremely substantially negatively correlated with yield; brown rice rate, milled rice rate, head rice rate, albumin content, globulin content, appearance, viscosity, and balance degree were substantially positively correlated with taste value; protein content and hardness were substantially negatively correlated with taste value (c).

4. Discussion

4.1. Effects of the Quantitative Application of N and K Fertilizers at Different Leaf Ages on Yield

A large number of studies have shown that adjusting the amount of N fertilizer, optimizing the proportion of N fertilizer in the fertilizer, and appropriately reducing the total amount of N fertilizer application can be more effective than conventional fertilization and improve rice yield [6,20]. It might increase yield by improving the N use efficiency and photosynthesis [32,33,34]. Previous studies have found that Chinese rice farmers used 56–85% of the total N in the season 10 days before transplanting, and most of the N was applied at the seedling and tillering stages, resulting in a waste of N [35]. Delaying fertilization time and reducing fertilization amount in the tillering stage can improve grain yield and NUE [36]. Appropriately delaying the application of tiller fertilizer and advancing the application of panicle fertilizer can increase the tillering panicle rate and the number of spikelets, thereby promoting the increase of tiller and maintaining the number of spikelets [21]. Since K cannot be covalently bound to organic molecules, it cannot be absorbed by plants such as N can [37]. However, K+ is the most abundant cellular cation and can regulate stomatal movement, osmotic regulation, and homeostatic enzyme activation and can participate in membrane protein transport [38]. Studies have shown that the response of K fertilizer to agronomic yield is the greatest when the application rate of K fertilizer reaches 130 mg kg−1 [39]. Atapattu A [40] reported that corn and soybeans showed positive returns at low K levels, whereas in moderate K-level soils, K returns were mostly negative. Singh B [41] pointed out that in only 12.5% of the paddy fields in the sample survey, the K application had a significant impact on rice yield. Another study showed that applying K fertilizer 50% (18.75 kg ha−1) more than the recommended amount at the heading stage and five days later than the recommended time could increase the yield of direct-seeded rice [42]. Compared with that on N, the research on the application rate and application period of K is slightly limited. Most of the studies are regarding the improvement of the disease resistance and stress resistance of rice to increase the yield, and the research on its direct effect on yield, especially the research on N and K combined application in the panicle fertilizer stage, is not in-depth or systematic. This study showed that K fertilizer application at the 12th-leaf (80% leaf age) during the panicle stage can significantly increase the number of effective tillers and increase the yield of rice by 23.20%. Application at the 13th- and 14th-leaf (90–95% leaf age) can significantly increase the seed setting rate and 1000-grain weight. With the increase of the N fertilizer and K fertilizer application amounts at the panicle stage, their effects on yield showed a trend of increasing first and then decreasing. According to the results of the equation simulation, the highest yield can be obtained when the application rates of N and K during the panicle stage reach 91.67 kg ha−1 and 55.00 kg ha−1, respectively.

4.2. Effects of the Quantitative Application of N and K Fertilizers at Different Leaf Ages on Various Quality Indicators

Previous studies show that the cooking and processing quality of rice can be improved by changing the main agronomic practices (fertilizing at the right time, rationally improving the amount of fertilizer, etc.) [42]. In addition, the application of comprehensive nutrients (N, P, K, organic fertilizers, etc.) to paddy soil has a significant impact on rice quality. K uptake by field crops is generally faster than N uptake, and K also plays an important role in ensuring efficient N use. K is not as easily absorbed into the organic matter as N and P, but K facilitates the transport of photosynthetic products and other plant metabolites, thereby contributing to improved grain quality, including increased rice aroma, smoothness, and whiter appearance [43]. In this study, it was found that the application of the appropriate panicle K fertilizer at the right time improved the processing quality of rice, that is, the application of K fertilizer at a dose of 47.81–64.69 kg ha−1 at the 12th- and 13th-leaf (80–90% leaf age) was when the processing quality was the best. When the application rate of the K fertilizer per panicle was lower than 47.81 kg ha−1 or higher than 64.69 kg ha−1, the processing quality of the rice was reduced; with the increase of panicle N fertilizer, the processing quality gradually improved, but at high K levels (higher than 56.25 kg ha), the processing quality gradually decreased with the increase of panicle N fertilizer.
Rice grains are made up of about 80–85% starch, 4–10% protein, etc. [44]. The content of amylose, the structure of amylopectin, and the interaction of starch with proteins and lipids are important factors that determine the quality of rice [45]. Among them, amylose content may be the most important factor in determining the cooking taste [46]. Rice with low amylose is usually soft and waxy, whereas white rice with high amylose is hard [47]. Although the protein content of rice is low, in terms of taste quality, some studies have shown that the total protein content is related to the texture of cooked rice [48], and some studies have shown that the difference in the RVA profile between varieties with similar amylose is because of the total protein content. This in turn affects the eating quality, and protein components may also affect eating quality [49]. At present, there are some studies on protein components, but the results are not uniform. Some studies show that globulin has no relationship with RVA eigenvalues, but its variation at very low levels plays a very important function; some studies show that the difference between the content of gluten and prolamin is the largest, and the ratio of these two proteins affects the taste quality [50]. In addition to the above two viewpoints, some believe that the interaction between different types of proteins is the root cause of their influence on the taste value of cooking [51]. A large number of studies have shown that the N application can significantly increase the protein content and starch content of rice [52] and significantly reduce the quality of rice [53], but the effect of K on rice quality is less. The 10th-leaf (70% leaf age) was applied with ear fertilizer, which significantly reduced the starch and protein content and improved the nutritional quality of rice. For the protein component, the early application of ear fertilizer (70% leaf age) decreased gluten content and the too-late application (90% leaf age) increased gluten content. With the increase of K application, the content of albumin and gluten was consistent; that is, it first decreased and then increased, reaching the minimum value at K3. While the content of globulin increased gradually, the content of gliadin changed irregularly. With the increase in the N application rate, albumin, globulin, prolamin, and glutenin all showed an increasing trend. However, there were some differences in the change of each protein component under the action of different K application rates, indicating that different N and K ratios had an impact on the protein components. From the RVA profile characteristics, PKV, CPV, and SBV gradually increased with the increase of leaf age, and the regularity of HPV and BDV was poor. When 10th-leaf (70% leaf age) was applied with ear fertilizer, appearance, viscosity, and balance degree significantly increased and hardness decreased, which in turn significantly improved the cooking and eating taste score, which was 11.08 points higher than that of CK. However, 14th-leaf application (90% leaf age) significantly reduced the cooking and eating taste score, which was 9.11 points lower than that of CK, and the difference between the 11th- and 13th-leaf application and CK was not significant. When the amount of K fertilizer applied at the panicle stage was lower than 56.25 kg ha−1, the content of starch and protein gradually decreased with the increase of application amount and increased gradually when it was higher than 56.25 kg ha−1. The higher the amount of N fertilizer applied at the panicle stage, the higher the starch content and the protein content and the lower the nutritional quality. In the two-year test results, the effect of K application at the panicle stage on the RVA profile characteristics was small and insignificant. However, with the increase of the N application at the panicle stage, it decreased PKV, HPV, CPV, and BDV; increased SBV; and decreased the cooking and eating taste score. When the amount of K applied at the panicle stage was less than 47.81 kg ha−1, the taste score also increased with the increase of the K application. When the K application at the panicle stage was higher than 47.81 kg ha−1, the taste value decreased with the increase of the K application, which affected the taste value by significantly changing appearance, viscosity, and balance degree. With the increase of the N application, the taste value showed a decreasing trend, which was consistent with previous studies, but under low K conditions, with the increase of the N application, the change of the taste value was small and insignificant.

4.3. Study on the Synergistic Improvement of Yield and Quality by the Quantitative Application of N and K Fertilizers at Different Leaf Ages during the Panicle Stage

At present, a large number of studies at this stage have separated rice yield and rice quality and have proposed many methods to increase yield or quality. Using advanced molecular techniques, it is indeed possible to find some genes that regulate yield and quality [54,55,56]. However, there are few studies at present regarding cultivation, especially whether it can have a synergistic effect on yield and quality through the use of fertilizers and precise fertilization. This study integrated the previous research on N management in this subject and drew some preliminary conclusions for the two-year-early N and K fertilizer management: N fertilizer was applied on the 12th-leaf, and K fertilizer was applied on 10th-, 11th-, and 12th-leaf. There was no synergistic correlation between yield and quality, and the indicators affecting the two were different. When the K fertilizer was applied on the 13th- and 14th-leaf, the relationship between yield and quality changed a little, and some indicators had the same effect on both and reached a significant level. This shows that the relationship between yield and quality can be coordinated under different fertilization methods at the early stage, so that the two can be improved synergistically. However, this improvement is still incomplete and low, and further research is needed. Although some preliminary conclusions were drawn in this study, it cannot be further explored due to too much processing. In the future, through reducing the treatment, more precise and in-depth research can be conducted on how to improve the yield and quality together from the perspective of cultivation and reveal its physiological mechanism.

Author Contributions

W.Z. and J.G. designed the study and provided experimental materials. X.W. and Y.L. (Yuzhuo Liu) performed the determination of protein components. Y.L. (Yuqi Liu), Y.F. and Y.Z. performed the determination of quality-related indicators. L.C. and W.X. analyzed the results and prepared the figures and tables. L.C. wrote the paper. All authors discussed the results and commented on the manuscript. All authors have read and agreed to the published version of the manuscript.

Funding

This study was supported by LiaoNing Revitalization Talents Program (No. XLYC2002073; No. XLYC2007169).

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Normile, D. Agricultural research. Reinventing rice to feed the world. Science 2008, 321, 5887. [Google Scholar] [CrossRef] [PubMed]
  2. Li, W.; Younger, P.; Cheng, Y.; Zhang, B.; Zhou, H.; Liu, Q. Addressing the CO2 emissions of the world’s largest coal producer and consumer: Lessons from the Haishiwan Coalfield, China. Energy 2015, 80, 400–413. [Google Scholar] [CrossRef] [Green Version]
  3. Xin, F.; Xiao, X.; Dong, J.; Zhang, G.; Li, B. Large increases of paddy rice area, gross primary production, and grain production in Northeast China during 2000-2017. Sci. Total Environ. 2019, 711, 135183. [Google Scholar] [CrossRef] [PubMed]
  4. Gilbert, P. Eutrophication, harmful algae and biodiversity—Challenging paradigms in a world of complex nutrient changes. Mar. Pollut. Bull. 2017, 124, 591–606. [Google Scholar] [CrossRef]
  5. Wang, Y.; Chen, Y.; Wu, W. Potassium and phosphorus transport and signaling in plants. J. Integr. Plant Biol. 2021, 63, 34–52. [Google Scholar] [CrossRef]
  6. Ye, T.; Li, Y.; Zhang, J.; Hou, W.; Zhou, W.; Lu, J.; Xing, Y.; Li, X. Nitrogen, phosphorus, and potassium fertilization affects the flowering time of rice (Oryza sativa L.). Glob. Ecol. Conserv. 2019, 20, e00753. [Google Scholar] [CrossRef]
  7. Armengaud, P.; Breitling, R.; Amtmann, A. The Potassium-Dependent Transcriptome of Arabidopsis Reveals a Prominent Role of Jasmonic Acid in Nutrient Signaling. Plant Physiol. 2004, 1, 2556–2576. [Google Scholar] [CrossRef] [Green Version]
  8. Kalaji, H.M.; Oukarroum, A.; Alexandrov, V.; Kouzmanova, M.; Brestic, M.; Zivcak, M. Identification of nutrient deficiency in maize and tomato plants by invivo chlorophyll a fluorescence measurements. Plant Physiol Biochem. 2014, 81, 16–25. [Google Scholar] [CrossRef]
  9. Wang, H.; Lu, Y.; Xu, J.; Liu, X.; Sheng, L. Effects of additives on nitrogen transformation and greenhouse gases emission of co-composting for deer manure and corn straw. Environ. Sci. Pollut. Res. 2021, 28, 13000–13020. [Google Scholar] [CrossRef]
  10. Christian, Z.; Mehmet, S.; Edgar, P. Potassium in agriculture—Status and perspectives. J. Plant Physiol. 2014, 171, 656–669. [Google Scholar] [CrossRef]
  11. Cheng, F.; Zhong, L.; Wang, F.; Zhang, G. Differences in cooking and eating properties between chalky and translucent parts in rice grains. Food Chem. 2005, 90, 39–46. [Google Scholar] [CrossRef]
  12. Zhang, Z.; Chu, G.; Liu, L.; Wang, Z.; Wang, X.; Zhang, H.; Yang, J.; Zhang, J. Mid-season nitrogen application strategies for rice varieties differing in panicle size. Field Crops Res. 2013, 150, 9–18. [Google Scholar] [CrossRef]
  13. Yoshida, H.; Horie, T.; Shiraiwa, T. A model explaining genotypic and environmental variation of rice spikelet number per unit area measured by cross-locational experiments in Asia. Field Crops Res. 2006, 97, 337–343. [Google Scholar] [CrossRef]
  14. Tan, P.; Anh, T.; Luat, N.; Puckridge, D. Yield trends of a long-term NPK experiment for intensive rice monoculture in the Mekong River Delta of Viet Nam. Field Crops Res. 1995, 42, 101–109. [Google Scholar] [CrossRef]
  15. Yao, X.; Zhang, Y.; Lu, Z.; Zhou, Y. A bibliometric review of nitrogen research in eutrophic lakes and reservoirs. J. Environ. Sci. 2017, 66, 274–285. [Google Scholar] [CrossRef]
  16. Zhao, L.; Li, K.; Wang, Q.; Song, X.; Su, H.; Xie, B. Nitrogen Starvation Impacts the Photosynthetic Performance of Porphyridium cruentum as Revealed by Chlorophyll a Fluorescence. Sci. Rep. 2017, 7, 8542. [Google Scholar] [CrossRef]
  17. Zhang, H.; Wu, G.; Huo, Z.; Xu, K. Precise Postponing Nitrogen Application and Its Mechanism in Rice. Acta Agron. Sin. 2011, 37, 1837–1851. [Google Scholar] [CrossRef]
  18. Ling, Q.; Zhang, H.; Dai, Q.; Ding, Y.; Wang, S. Study on Precise and Quantitative N Application in Rice. Sci. Agric. Sin. 2005, 38, 2457–2467. [Google Scholar]
  19. Yin, C.; Zhang, H.; Zhang, Q.; Wei, H.; Dai, Q.; Huo, Z. Preliminary Study on Parameters of Precise and Quantitative Nitrogen Application in Rice Varieties with Different Growth Period Durations. Acta Agron. Sin. 2010, 43, 39–50. (In Chinese) [Google Scholar] [CrossRef]
  20. Yi, J.; Gao, J.; Zhang, W.; Zhao, Y.Z.; Xin, W. Delayed timing of tillering fertilizer improved grain yield and nitrogen use efficiency in japonica rice. Crop Sci. 2020, 60, 1021–1033. [Google Scholar] [CrossRef]
  21. Wang, Y.; Yi, J.; Gao, J.; Zhang, L.; Yang, J.; Zhao, Y.; Xin, W.; Zhen, X.; Zhang, W. Effects of precision leaf age fertilization on yield and nitrogen utilization of japonica rice. Crop J. 2020, 46, 102–116. (In Chinese) [Google Scholar] [CrossRef]
  22. Xin, W.; Zhang, L.; Zhang, W.; Gao, J.; Yi, J.; Zhen, X.; Du, M.; Zhao, Y.; Chen, L. Morphological and Physiological Characteristics of Rice Cultivars with Higher Yield and Nitrogen Use Efficiency at Various Nitrogen Rates. Agronomy 2022, 12, 358. [Google Scholar] [CrossRef]
  23. Du, M.; Zhang, W.; Gao, J.; Liu, M.; Zhou, Y.; He, D.; Zhao, Y. Improvement of Root Characteristics Due to Nitrogen, Phosphorus, and Potassium Interactions Increases Rice (Oryza sativa L.) Yield and Nitrogen Use Efficiency. Agronomy 2022, 12, 23. [Google Scholar] [CrossRef]
  24. Hirooka, Y.; Homma, K.; Shiraiwa, T.; Makino, Y.; Liu, T.; Xu, Z. Yield and growth characteristics of erect panicle type rice (Oryza sativa L.) cultivar, Shennong265 under various crop management practices in Western Japan. Plant Prod. Sci. 2018, 21, 1–7. [Google Scholar] [CrossRef] [Green Version]
  25. Huang, X.; Qian, Q.; Liu, Z.; Sun, H.; He, S.; Luo, D. Natural variation at the DEP1 locus enhances grain yield in rice. Nat. Genet. 2009, 41, 494–497. [Google Scholar] [CrossRef]
  26. Sun, H.; Qian, Q.; Wu, K.; Luo, J.; Wang, S. Heterotrimeric G proteins regulate nitrogen-use efficiency in rice. Nat. Genet. 2014, 46, 652–656. [Google Scholar] [CrossRef]
  27. Bao, S. Soil and Agricultural Chemistry Analysis; China Agriculture Press: Beijing, China, 2000; pp. 33–35, 44–49, 56–58, 81, 105–107. (In Chinese) [Google Scholar]
  28. Luthe, D. Storage protein accumulation in developing rice (Oryza sativa L.) seeds. Plant Sci. Lett. 1983, 32, 147–158. [Google Scholar] [CrossRef]
  29. Liu, Z.; Cheng, F.; Cheng, W.; Zhang, G. Positional variations in phytic acid and protein content within a panicle of japonica rice. J. Cereal Sci. 2005, 41, 297–303. [Google Scholar] [CrossRef]
  30. Champagne, E.; Bett, K.; Vinyard, B.; Mcclung, A.; Barton, F.; Moldenhauer, K. Correlation Between Cooked Rice Texture and Rapid Visco Analyser Measurements. Cereal Chem. 1999, 76, 764–771. [Google Scholar] [CrossRef]
  31. Cao, X.; Sun, H.; Wang, C.; Ren, X.; Liu, H.; Zhang, Z. Effects of late-stage nitrogen fertilizer application on the starch structure and cooking quality of rice. J. Sci. Food Agric. 2017, 98, 2332–2340. [Google Scholar] [CrossRef]
  32. Hou, W.; Yan, J.; Jákli, B.; Lu, J.; Ren, T.; Cong, R. Synergistic effects of nitrogen and potassium on quantitative limitations to photosynthesis in rice (Oryza sativa L.). J. Agric. Food Chem. 2018, 66, 5125–5132. [Google Scholar] [CrossRef] [PubMed]
  33. Fageria, N. Nutrient management for improving low land rice productivity and sustainability. Adv. Agron. 2003, 80, 63–152. [Google Scholar] [CrossRef]
  34. Srivastava, P.; Srivastava, P.; Singh, U.; Shrivastava, M. Effect of integrated and balanced nutrient application on soil fertility, yield and quality of Basmati rice. Arch. Agron. Soil Sci. 2009, 55, 265–284. [Google Scholar] [CrossRef]
  35. Sui, B.; Feng, X.; Tian, G.; Hu, X.; Guo, S. Optimizing nitrogen supply increases rice yield and nitrogen use efficiency by regulating yield formation factors. Field Crops Res. 2013, 150, 99–107. [Google Scholar] [CrossRef]
  36. Yan, S.; Wm, D.; Ls, A.; Ys, A.; Lw, C. Controlled-release fertilizer enhances rice grain yield and N recovery efficiency in continuous non-flooding plastic film mulching cultivation system. Field Crops Res. 2019, 231, 122–129. [Google Scholar] [CrossRef]
  37. Popp, M.; Slaton, N.; Norsworthy, J.; Dixon, B. Rice yield response to potassium: An economic analysis. Agron. J. 2020, 113, 287–297. [Google Scholar] [CrossRef]
  38. Mallarino, A.; Webb, J.; Blackmer, A. Soil test values and grain yields during 14 years of potassium fertilization of corn and soybean. J. Prod. Agric. 1991, 4, 560–567. [Google Scholar] [CrossRef]
  39. Fryer, M.; Slaton, N.; Roberts, T.; Ross, W. Validation of Soil-Test-Based Phosphorus and Potassium Fertilizer Recommendations for Irrigated Soybean. Soil Sci. Soc. Am. J. 2019, 83, 825–837. [Google Scholar] [CrossRef] [Green Version]
  40. Atapattu, A.; Prasantha, B.; Amaratunga, K.; Marambe, B. Increased rate of potassium fertilizer at the time of heading enhance the quality of direct-seeded rice. Chem. Biol. Technol. Agric. 2018, 5, 1–9. [Google Scholar] [CrossRef]
  41. Singh, B.; Singh, Y.; Imas, P.; Xie, J. Potassium Nutrition of the Rice–Wheat Cropping System. Adv. Agron. 2003, 81, 203–259. [Google Scholar] [CrossRef]
  42. Armengaud, P.; Sulpice, R.; Miller, A.; Stitt, M.; Amtmann, A.; Gibon, Y. Multilevel Analysis of Primary Metabolism Provides New Insights into the Role of Potassium Nutrition for Glycolysis and Nitrogen Assimilation in Arabidopsis Roots. Plant Physiol. 2009, 2, 772–785. [Google Scholar] [CrossRef] [Green Version]
  43. Wang, Q.; Zhen, R.; Ding, Y.; Ji, Z.; Huang, P. Effects of potassium fertilizer application rates on plant potassium accumulation and grain quality of japonica rice. Sci. Agric. Sin. 2004, 37, 1444–1450. (In Chinese) [Google Scholar]
  44. Bahmaniar, M.; Ranjbar, G. Response of rice (Oryza sativa L.) cooking quality properties to nitrogen and potassium application. Pak. J. Biol. Sci. 2007, 10, 1880–1884. [Google Scholar] [CrossRef]
  45. Balindong, J.; Ward, R.; Liu, L.; Rose, T.; Pallas, L.; Ovenden, B.; Waters, D. Rice grain protein composition influences instrumental measures of rice cooking and eating quality. J. Cereal Sci. 2018, 79, 35–42. [Google Scholar] [CrossRef] [Green Version]
  46. Umemoto, T.; Aoki, N.; Lin, H.; Nakamura, Y.; Maruyama, S. Natural variation in rice starch synthase IIa affects enzyme and starch properties. Funct. Plant Biol. 2004, 31, 671–684. [Google Scholar] [CrossRef]
  47. Waters, D.; Fitzgerald, M.; Henry, R. Gelatinization temperature of rice explained by polymorphisms in starch synthase. Plant Biotechnol. J. 2006, 4, 115–122. [Google Scholar] [CrossRef]
  48. Noda, T.; Nishiba, Y.; Sato, T.; Suda, I. Properties of starches from several low-amylose rice cultivars. Cereal Chem. 2003, 80, 193–197. [Google Scholar] [CrossRef]
  49. Lyon, B.; Champagne, E.; Vinyard, B.; Windham, W.; Barton, F.; Webb, B. Effects of degree of milling, drying condition, and final moisture content on sensory texture of cooked rice. Cereal Chem. 1999, 76, 56–62. [Google Scholar] [CrossRef]
  50. Baxter, G.; Zhao, J.; Blanchard, C. Albumin Significantly Affects Pasting and Textural Characteristics of Rice Flour. Cereal Chem. 2010, 87, 250–255. [Google Scholar] [CrossRef]
  51. Zhou, Z.; Robards, K.; Helliwell, S.; Blanchard, C. Effect of storage temperature on rice thermal properties. Food Res. Int. 2010, 43, 709–715. [Google Scholar] [CrossRef]
  52. Bharali, A.; Baruah, K. Effects of integrated nutrient management on sucrose phosphate synthase enzyme activity and grain quality traits in rice. Physiol. Mol. Biol. Plants 2020, 28, 383–389. [Google Scholar] [CrossRef]
  53. Müller, A.; Coradi, P.; Nunes, M.; Grohs, M.; Bressiani, J.; Teodoro, P.; Flores, E. Effects of cultivars and fertilization levels on the quality of rice milling: A diagnosis using near-infrared spectroscopy, X-ray diffraction, and scanning electron microscopy. Food Res. Int. 2021, 147, 110524. [Google Scholar] [CrossRef]
  54. Ahmad, I.; Mian, A.; Maathuis, F. Overexpression of the rice AKT1 potassium channel affects potassium nutrition and rice drought tolerance. J. Exp. Bot. 2016, 67, 2689–2698. [Google Scholar] [CrossRef] [Green Version]
  55. Wu, Y.; Hu, Y.; Xu, G. Interactive effects of potassium and sodium on root growth and expression of K/Na transporter genes in rice. Plant Growth Regul. 2009, 57, 271–280. [Google Scholar] [CrossRef]
  56. Yang, T.; Zhang, S.; Hu, Y.; Wu, F.; Hu, Q.; Chen, G. The Role of a Potassium Transporter OsHAK5 in Potassium Acquisition and Transport from Roots to Shoots in Rice at Low Potassium Supply Levels. Plant Physiol. 2014, 166, 945–959. [Google Scholar] [CrossRef] [Green Version]
Figure 1. Mean temperature (°C), relative humidity (%), precipitation (mm), and sunshine hours (h) during the rice-growing season in 2020 and 2021.
Figure 1. Mean temperature (°C), relative humidity (%), precipitation (mm), and sunshine hours (h) during the rice-growing season in 2020 and 2021.
Agronomy 12 01681 g001
Figure 2. Effects of N and K application on rice yield and yield components at different leaf ages patterns in 2020 and 2021, with the results of the two-factor ANOVA. The (10,12), etc., represent the patterns of different applications of N and K at the panicle stage; CK is the fertilization pattern of local farmers. Error bars indicate standard errors of three replicates. Different lowercase letters above the columns indicate significant differences (p < 0.05) among different panicle fertilizer application patterns in the same year.
Figure 2. Effects of N and K application on rice yield and yield components at different leaf ages patterns in 2020 and 2021, with the results of the two-factor ANOVA. The (10,12), etc., represent the patterns of different applications of N and K at the panicle stage; CK is the fertilization pattern of local farmers. Error bars indicate standard errors of three replicates. Different lowercase letters above the columns indicate significant differences (p < 0.05) among different panicle fertilizer application patterns in the same year.
Agronomy 12 01681 g002
Figure 3. Effects of N and K application on rice processing quality at different leaf ages patterns in 2020 and 2021, with the results of the two-factor ANOVA. The (10,12), etc., represent the patterns of different applications of N and K at the panicle stage; CK is the fertilization pattern of local farmers. Error bars indicate standard errors of three replicates. Different lowercase letters above the columns indicate significant differences (p < 0.05) among different panicle fertilizer application patterns in the same year.
Figure 3. Effects of N and K application on rice processing quality at different leaf ages patterns in 2020 and 2021, with the results of the two-factor ANOVA. The (10,12), etc., represent the patterns of different applications of N and K at the panicle stage; CK is the fertilization pattern of local farmers. Error bars indicate standard errors of three replicates. Different lowercase letters above the columns indicate significant differences (p < 0.05) among different panicle fertilizer application patterns in the same year.
Agronomy 12 01681 g003
Figure 4. Effects of N and K application on rice amylose, protein component at different leaf ages patterns in 2020 and 2021, with the results of the two-factor ANOVA. The (10,12), etc., represent the patterns of different applications of N and K at the panicle stage; CK is the fertilization pattern of local farmers. Error bars indicate standard errors of three replicates. Different lowercase letters above the columns indicate significant differences (p < 0.05) among different panicle fertilizer application patterns in the same year.
Figure 4. Effects of N and K application on rice amylose, protein component at different leaf ages patterns in 2020 and 2021, with the results of the two-factor ANOVA. The (10,12), etc., represent the patterns of different applications of N and K at the panicle stage; CK is the fertilization pattern of local farmers. Error bars indicate standard errors of three replicates. Different lowercase letters above the columns indicate significant differences (p < 0.05) among different panicle fertilizer application patterns in the same year.
Agronomy 12 01681 g004
Figure 5. Effects of N and K application on RVA profile characteristics under different leaf age patterns in 2020 and 2021, with the results of the two-factor ANOVA. The (10,12), etc., represent the patterns of different applications of N and K at the panicle stage; CK is the fertilization pattern of local farmers. Error bars indicate standard errors of three replicates. Different lowercase letters above the columns indicate significant differences (p < 0.05) among different panicle fertilizer application patterns in the same year.
Figure 5. Effects of N and K application on RVA profile characteristics under different leaf age patterns in 2020 and 2021, with the results of the two-factor ANOVA. The (10,12), etc., represent the patterns of different applications of N and K at the panicle stage; CK is the fertilization pattern of local farmers. Error bars indicate standard errors of three replicates. Different lowercase letters above the columns indicate significant differences (p < 0.05) among different panicle fertilizer application patterns in the same year.
Agronomy 12 01681 g005
Figure 6. Effects of N and K application on cooking quality under different leaf age patterns in 2020 and 2021, with the results of the two-factor ANOVA. The (10,12), etc., represent the patterns of different applications of N and K at the panicle stage; CK is the fertilization pattern of local farmers. Error bars indicate standard errors of three replicates. Different lowercase letters above the columns indicate significant differences (p < 0.05) among different panicle fertilizer application patterns in the same year.
Figure 6. Effects of N and K application on cooking quality under different leaf age patterns in 2020 and 2021, with the results of the two-factor ANOVA. The (10,12), etc., represent the patterns of different applications of N and K at the panicle stage; CK is the fertilization pattern of local farmers. Error bars indicate standard errors of three replicates. Different lowercase letters above the columns indicate significant differences (p < 0.05) among different panicle fertilizer application patterns in the same year.
Agronomy 12 01681 g006
Figure 7. Correlation between yield and K fertilizer application at different leaf ages during the panicle stage. The (10,12), etc., represent the patterns of different applications of N and K at the panicle stage. * p < 0.05.
Figure 7. Correlation between yield and K fertilizer application at different leaf ages during the panicle stage. The (10,12), etc., represent the patterns of different applications of N and K at the panicle stage. * p < 0.05.
Agronomy 12 01681 g007
Figure 8. Correlation between yield and N fertilizer application at different leaf ages during the panicle stage. The (10,12), etc., represent the patterns of different applications of N and K at the panicle stage. * p < 0.05.
Figure 8. Correlation between yield and N fertilizer application at different leaf ages during the panicle stage. The (10,12), etc., represent the patterns of different applications of N and K at the panicle stage. * p < 0.05.
Agronomy 12 01681 g008
Figure 9. Correlation between taste value and K fertilizer application at different leaf ages during the panicle stage. The (10,12), etc., represent the patterns of different applications of N and K at the panicle stage. * p < 0.05.
Figure 9. Correlation between taste value and K fertilizer application at different leaf ages during the panicle stage. The (10,12), etc., represent the patterns of different applications of N and K at the panicle stage. * p < 0.05.
Agronomy 12 01681 g009
Figure 10. Correlation between taste value and N fertilizer application at different leaf ages during the panicle stage. The (10,12), etc., represent the patterns of different applications of N and K at the panicle stage. * p < 0.05.
Figure 10. Correlation between taste value and N fertilizer application at different leaf ages during the panicle stage. The (10,12), etc., represent the patterns of different applications of N and K at the panicle stage. * p < 0.05.
Agronomy 12 01681 g010
Figure 11. Correlation analysis of K fertilizers applied at different leaf ages during the panicle stage on yield and quality indexes. * p < 0.05. Note: X1: NO. of panicle per plant, X2: NO. of grains per panicle, X3: seed-setting rate, X4: 1000-grain weight, X5: brown rice rate, X6: milled rice rate, X7: head rice rate, X8: amylose contents, X9: protein contents, X10: albumin, X11: globulin, X12: prolamin, X13: glutelin, X14: peak viscosity, X15: hot paste viscosity, X16: cool paste viscosity, X17: breakdown viscosity, X18: setback viscosity, X19: appearance, X20: hardness, X21: viscosity, X22: balance degree, X23: taste value, X24: grain yield per plant.
Figure 11. Correlation analysis of K fertilizers applied at different leaf ages during the panicle stage on yield and quality indexes. * p < 0.05. Note: X1: NO. of panicle per plant, X2: NO. of grains per panicle, X3: seed-setting rate, X4: 1000-grain weight, X5: brown rice rate, X6: milled rice rate, X7: head rice rate, X8: amylose contents, X9: protein contents, X10: albumin, X11: globulin, X12: prolamin, X13: glutelin, X14: peak viscosity, X15: hot paste viscosity, X16: cool paste viscosity, X17: breakdown viscosity, X18: setback viscosity, X19: appearance, X20: hardness, X21: viscosity, X22: balance degree, X23: taste value, X24: grain yield per plant.
Agronomy 12 01681 g011
Table 1. The amount and timing of fertilization.
Table 1. The amount and timing of fertilization.
K
Treatment
N
Treatment
Basal Fertilizer
(kg ha−1)
Tillering Fertilizer
(kg ha−1)
Panicle Fertilizer
(kg ha−1)
Application Time
NPKNNK
K1N081.00112.5056.2554.00039.38Basel fertilizer was applied one day before transplanting. Tillering fertilizer was applied on 9th-leaf. Panicle N fertilizer was applied on 12th-leaf. Panicle K fertilizer was applied on 10th-leaf, 11th-leaf, 12th-leaf, 13th-leaf, or 14th-leaf, respectively, and only once.
N181.00112.5056.2554.0072.0039.38
N281.00112.5056.2554.0090.0039.38
N381.00112.5056.2554.00108.0039.38
K2N081.00112.5056.2554.00047.81
N181.00112.5056.2554.0072.0047.81
N281.00112.5056.2554.0090.0047.81
N381.00112.5056.2554.00108.0047.81
K3N081.00112.5056.2554.00056.25
N181.00112.5056.2554.0072.0056.25
N281.00112.5056.2554.0090.0056.25
N381.00112.5056.2554.00108.0056.25
K4N081.00112.5056.2554.00064.69
N181.00112.5056.2554.0072.0064.69
N281.00112.5056.2554.0090.0064.69
N381.00112.5056.2554.00108.0064.69
K5N081.00112.5056.2554.00073.13
N181.00112.5056.2554.0072.0073.13
N281.00112.5056.2554.0090.0073.13
N381.00112.5056.2554.00108.0073.13
CK81.00112.5056.2554.0090.0056.25Basal fertilizer was applied one day before transplanting. Tiller fertilizer was applied 7 days after transplanting. Panicle fertilizer was applied 60 days after transplanting.
Table 2. Effects of different N, K, and Y combined application during the panicle stage on rice yield and yield components in 2020 and 2021, with the results of the three-way ANOVA.
Table 2. Effects of different N, K, and Y combined application during the panicle stage on rice yield and yield components in 2020 and 2021, with the results of the three-way ANOVA.
YearTreatmentNO. of Panicle per PlantNO. of Grains per PanicleSeed-Setting Rate
(%)
1000-Grain Weight
(g)
Grain Yield (g∙m−2)
2020K1N112.53 b155.83 a93.84 a23.12 a840.71 b
K1N213.47 ab152.55 a94.34 a23.60 a929.47 b
K1N314.60 a157.34 a93.49 a23.65 a1087.32 a
K2N112.87 b149.71 b94.47 a24.23 a958.99 b
K2N213.9 ab156.67 ab93.07 a24.22 a1097.98 a
K2N314.47 a165.59 a93.36 a24.30 a1151.49 a
K3N115.73 a178.96 a91.91 a23.70 a1080.56 a
K3N214.13 b173.01 a93.54 a24.06 a1080.56 a
K3N313.93 b168.88 a95.53 a23.69 a1009.22 a
K4N116.07 a156.71 a93.58 b24.38 a1089.37 a
K4N214.60 b155.49 ab94.30 ab24.38 a973.14 b
K4N314.07 b152.74 b95.45 a24.09 a924.55 b
K5N113.13 a136.44 a95.47 a25.41 a1038.74 a
K5N213.00 a134.43 ab96.94 a25.98 a982.36 ab
K5N311.87 a125.47 b96.16 a25.58 a894.21 b
2021K1N111.73 b152.13 a94.12 a24.07 a914.51 a
K1N212.52 ab153.12 a94.03 a24.32 a939.31 a
K1N313.46 a155.27 a93.78 a24.62 a945.05 a
K2N111.92 b150.27 b94.75 a24.72 a965.76 b
K2N212.91 ab153.41 ab94.42 a24.75 a1009.01 ab
K2N313.88 a159.57 a94.07 a24.73 a1064.98 a
K3N115.01 ab160.23 a93.21 b24.31 a1040.17 a
K3N213.92 b160.52 a94.33 a24.26 a1022.54 a
K3N312.73 b163.22 a94.36 a24.25 a998.56 a
K4N115.32 a150.23 a94.37 c24.00 a925.99 a
K4N214.09 b154.16 a94.89 b24.78 a877.40 b
K4N313.93 b150.36 a95.41 a24.89 a867.77 b
K5N112.42 a140.17 a96.31 b25.78 a954.07 a
K5N212.27 a137.24 ab96.48 ab25.72 a911.23 a
K5N312.00 a124.25 b97.00 a25.99 a925.78 a
F-valueK**********
N**nsnsnsns
Y**nsns*
K × Nns****ns*
N × Ynsnsnsnsns
K × Yns*nsnsns
K × N × Ynsnsnsnsns
Note: Different letters indicate differences between different N levels at the same K level. * p < 0.05; ** p < 0.01; ns, not substantial. K, potassium rate; N, nitrogen rate; Y, year.
Table 3. Effects of different N and K combined application rates during the panicle stage on rice processing quality in 2020 and 2021, with the results of the three-way ANOVA.
Table 3. Effects of different N and K combined application rates during the panicle stage on rice processing quality in 2020 and 2021, with the results of the three-way ANOVA.
YearTreatmentBrown Rice Rate
(%)
Milled Rice Rate
(%)
Head Rice Rate
(%)
2020K1N167.67 b60.57 b56.82 b
K1N268.25 b61.26 ab57.17 b
K1N369.96 a62.00 a59.40 a
K2N168.65 b61.25 b57.50 b
K2N270.00 ab63.01 ab60.15 a
K2N371.02 a64.60 a62.00 a
K3N171.54 a64.23 a61.09 a
K3N271.32 a64.11 a60.32 a
K3N370.02 b63.92 a59.62 b
K4N170.23 a64.15 a61.01 a
K4N268.39 b63.08 a59.25 ab
K4N366.57 b62.52 a57.56 b
K5N168.32 a63.10 a59.95 a
K5N266.40 a61.00 ab54.52 b
K5N364.63 b60.03 b53.97 b
2021K1N168.66 b62.63 b57.37 b
K1N270.14 b64.56 ab59.17 b
K1N372.30 a65.94 a61.40 a
K2N169.31 b63.27 b58.45 b
K2N270.40 ab64.01 ab61.15 a
K2N371.92 a66.60 a62.40 a
K3N172.64 a66.19 a62.13 a
K3N272.31 a65.92 a61.09 ab
K3N369.41 b65.73 a59.43 b
K4N172.41 a65.45 a62.06 a
K4N270.39 b64.18 a60.35 ab
K4N369.33 b63.74 a59.56 b
K5N170.54 a64.10 a60.95 a
K5N270.18 a61.87 ab56.60 b
K5N367.53 b60.21 b54.77 b
F-valueK*****
Nnsnsns
Y******
K × N*****
N × Ynsnsns
K × Y**ns
K × N × Ynsnsns
Note: Different letters indicate differences between different N levels at the same K level. * p < 0.05; ** p < 0.01; ns, not substantial. K, potassium rate; N, nitrogen rate; Y, year.
Table 4. Effects of different N and K combined application rates during the panicle stage on rice amylose and protein components in 2020 and 2021, with the results of the three-way ANOVA.
Table 4. Effects of different N and K combined application rates during the panicle stage on rice amylose and protein components in 2020 and 2021, with the results of the three-way ANOVA.
YearTreatmentAmylose Contents (%)Protein Contents (%)Albumin
(mg/grain)
Globulin
(mg/grain)
Prolamin
(mg/grain)
Glutelin
(mg/grain)
2020K1N118.20 a7.10 b2.17 c1.50 b0.40 c6.60a
K1N218.35 a7.16 ab2.65 b1.72 b0.72 b7.11a
K1N318.48 a7.31 a2.95 a2.23 a1.05 a7.11a
K2N117.90 b7.23 a2.25 b1.94 c0.25 c5.95c
K2N218.13 a7.29 a2.30 b2.36 b0.31 b6.50b
K2N318.23 a7.32 a2.61 a2.71 a0.90 a7.00 a
K3N117.33 b6.75 b1.85 c2.12 b0.15 b5.58 b
K3N217.60 a6.85 ab2.17 b2.56 a0.32 a6.29 a
K3N317.70 a7.01 a2.50 a2.73 a0.52 a6.76 a
K4N118.10 b7.15 a2.06 c2.46 b0.28 b6.75 a
K4N218.52 a7.19 a2.42 b2.87 b0.47 a6.89 a
K4N318.78 a7.25 a3.10 a3.12 a0.60 a6.99 a
K5N118.57 b7.10 b2.30 c2.59 c0.43 c6.68 b
K5N218.78 a7.25 b2.71 b3.41 b0.79 b6.89 b
K5N318.90 a7.57 a3.25 a4.00 a1.12 a7.61 a
2021K1N117.73 a6.97 b2.11 c1.56 b0.42 c6.63 a
K1N217.90 a7.03 ab2.72 b1.75 b0.65 b7.08 a
K1N318.00 a7.20 a3.01 a2.36 a1.11 a7.10 a
K2N117.47 b7.03 a2.20 b1.98 c0.20 c5.93 c
K2N218.00 a7.13 a2.31 b2.35 b0.37 b6.52 b
K2N318.07 a7.17 a2.66 a2.75 a0.89 a7.13 a
K3N117.13 b6.70 b1.88 c2.08 b0.14 b5.62 b
K3N217.80 a6.83 b2.16 b2.57 a0.35 a6.29 a
K3N317.90 a7.03 a2.52 a2.69 a0.51 a6.78 a
K4N117.80 b7.07 a2.02 c2.44 b0.26 b6.75 a
K4N218.23 a7.13 a2.50 b2.66 b0.46 a6.45 a
K4N318.40 a7.17 a3.60 a3.04 a0.59 a6.86 a
K5N118.07 b7.13 b2.40 c2.57 c0.41 c6.71 b
K5N218.37 a7.27 b2.81 b3.37 b0.72 b6.89 b
K5N318.40 a7.60 a3.34 a3.96 a1.09 a7.62 a
F-valueK************
N************
Y****nsnsnsns
K × Nns********
N × Ynsnsns*nsns
K × Yns*nsnsnsns
K × N × Yns*nsnsnsns
Note: Different letters indicate differences between different N levels at the same K level. * p < 0.05; ** p < 0.01; ns, not substantial. K, potassium rate; N, nitrogen rate; Y, year.
Table 5. Effects of different N and K combined application rates during the panicle stage on RVA profile characteristics in 2020 and 2021, with the results of the three-way ANOVA.
Table 5. Effects of different N and K combined application rates during the panicle stage on RVA profile characteristics in 2020 and 2021, with the results of the three-way ANOVA.
YearTreatmentPeak Viscosity (cP)Hot Paste Viscosity (cP)Cool Paste Viscosity (cP)Breakdown Viscosity (cP)Setback Viscosity (cP)
2020K1N12741.70 a2232.70 a3351.10 a509.00 a599.40 b
K1N22704.50 ab2145.20 a3307.70 a502.30 a603.20 ab
K1N32600.90 b2106.60 a3249.30 b500.30 a648.40 a
K2N12842.50 a2213.20 a3502.80 a629.30 a660.30 a
K2N22800.50 ab2183.00 a3442.20 ab617.50 ab671.70 a
K2N32701.50 b2127.40 a3304.80 a574.10 b703.30 a
K3N12800.80 a2202.70 a3402.10 a648.10 a601.30 a
K3N22782.10 a2152.20 a3357.50 a639.90 a605.40 a
K3N32724.20 a2124.90 a3332.60 a599.30 b638.40 a
K4N12824.40 a2222.10 a3498.40 a602.30 a674.00 b
K4N22801.20 ab2200.70 a3461.30 a600.50 a680.10 b
K4N32713.30 b2109.30 a3439.70 a613.00 a726.40 a
K5N12871.40 a2249.70 a3509.60 a621.70 a638.20 b
K5N22853.30 a2212.50 a3462.20 ab621.80 a638.90 b
K5N32702.30 b2106.30 a3430.50 b596.00 b728.20 a
2021K1N12836.10 a2251.80 a3499.60 a594.80 a645.50 a
K1N22761.30 ab2166.50 ab3470.70 ab584.30 a663.50 a
K1N32622.90 b2057.80 b3323.10 b565.10 a709.40 a
K2N12902.50 a2137.00 a3586.80 a642.80 a660.70 a
K2N22828.50 ab2183.00 a3489.20 ab581.50 b684.30 a
K2N32718.50 b2256.30 a3364.00 b572.20 b700.20 a
K3N12834.80 a2244.90 a3481.40 a636.50 a646.60 a
K3N22805.40 a2168.90 a3479.40 a589.90 ab661.10 a
K3N32734.30 a2156.90 a3395.40 a577.40 b674.00 a
K4N12866.80 a2236.30 a3512.90 a613.60 a650.80 b
K4N22816.10 ab2211.70 a3541.80 a604.40 a712.20 ab
K4N32690.30 b2123.30 a3402.50 a567.00 a725.70 a
K5N12893.00 a2270.70 a3537.70 a622.30 a644.70 b
K5N22871.20 a2262.50 a3554.20 a608.70 a683.00 ab
K5N32752.00 a2166.30 a3470.50 a585.70 a718.50 a
F-valueKnsnsnsnsns
N********
Y***ns*
K × Nns*nsnsns
N × Ynsnsnsns*
K × Ynsnsns**
K × N × Ynsnsns**
Note: Different letters indicate differences between different N levels at the same K level. * p < 0.05; ** p < 0.01; ns, not substantial. K, potassium rate; N, nitrogen rate; Y, year.
Table 6. Effects of different N and K combined application rates during the panicle stage on cooking quality in 2020 and 2021, with the results of the three-way ANOVA.
Table 6. Effects of different N and K combined application rates during the panicle stage on cooking quality in 2020 and 2021, with the results of the three-way ANOVA.
YearTreatmentAppearanceHardnessViscosityBalance DegreeTaste Value
2020K1N13.50 a7.82 a3.92 a3.32 a49.13 a
K1N23.21 a7.93 a3.52 b3.22 ab47.62 ab
K1N33.21 a8.01 a3.32 b3.18 b46.22 b
K2N13.72 a7.72 b4.11 a3.69 a50.25 a
K2N23.41 ab7.91 a3.96 ab3.53 ab49.13 ab
K2N33.20 b8.03 a3.72 b3.24 b48.25 b
K3N13.42 a7.81 b4.02 a3.41 a49.23 a
K3N23.44 a8.02 ab3.72 b3.22 ab47.14 b
K3N32.96 b8.17 a3.53 b3.01 b47.25 b
K4N13.62 a7.61 c4.12 a3.32 a51.12 a
K4N23.07 b8.03 b3.45 b3.17 a49.31 b
K4N32.71 c8.32 a3.02 c2.71 b47.11 c
K5N13.52 a7.71 b3.86 a3.41 a50.92 a
K5N22.72 b8.03 a3.21 b2.82 b47.22 b
K5N32.42 b8.25 a2.66 c2.51 c44.22 c
2021K1N13.61 a8.01 a4.03 a3.52 a50.90 a
K1N23.30 a8.15 a3.62 b3.15 b50.65 a
K1N33.32 a8.19 a3.59 b3.18 ab48.70 a
K2N13.77 a8.00 b4.29 a3.73 a52.15 a
K2N23.55 ab8.02 b3.88 b3.23 b48.75 b
K2N33.24 b8.23 a3.79 b3.43 ab50.00 ab
K3N13.51 a8.08 b4.06 a3.51 a50.65 a
K3N23.50 a8.17 b3.50 b2.89 b46.95 b
K3N32.92 b8.38 a3.77 ab3.26 a49.25 a
K4N13.81 a7.90 c4.23 a3.77 a52.40 a
K4N23.31 b8.15 b3.54 b3.16 b48.65 b
K4N32.59 c8.55 a3.09 c2.54 c44.30 c
K5N13.65 a7.97 b3.96 a3.57 a51.20 a
K5N22.92 b8.33 a3.13 b2.75 b46.40 b
K5N32.58 b8.51 a2.79 c2.41 b43.85 c
F-valueK**********
N**********
Ynsns*nsns
K × N**********
N × Ynsnsnsnsns
K × Ynsnsnsnsns
K × N × Ynsns*nsns
Note: Different letters indicate differences between different N levels at the same K level. * p < 0.05; ** p < 0.01; ns, not substantial. K, potassium rate; N, nitrogen rate; Y, year.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Chen, L.; Zhang, W.; Gao, J.; Liu, Y.; Wang, X.; Liu, Y.; Feng, Y.; Zhao, Y.; Xin, W. Effects of Precise K Fertilizer Application on the Yield and Quality of Rice under the Mode of Light, Simple, and High-Efficiency N Fertilizer Application during the Panicle Stage. Agronomy 2022, 12, 1681. https://0-doi-org.brum.beds.ac.uk/10.3390/agronomy12071681

AMA Style

Chen L, Zhang W, Gao J, Liu Y, Wang X, Liu Y, Feng Y, Zhao Y, Xin W. Effects of Precise K Fertilizer Application on the Yield and Quality of Rice under the Mode of Light, Simple, and High-Efficiency N Fertilizer Application during the Panicle Stage. Agronomy. 2022; 12(7):1681. https://0-doi-org.brum.beds.ac.uk/10.3390/agronomy12071681

Chicago/Turabian Style

Chen, Liqiang, Wenzhong Zhang, Jiping Gao, Yuzhuo Liu, Xue Wang, Yuqi Liu, Yingying Feng, Yanze Zhao, and Wei Xin. 2022. "Effects of Precise K Fertilizer Application on the Yield and Quality of Rice under the Mode of Light, Simple, and High-Efficiency N Fertilizer Application during the Panicle Stage" Agronomy 12, no. 7: 1681. https://0-doi-org.brum.beds.ac.uk/10.3390/agronomy12071681

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