Skip to main content
Advertisement
Browse Subject Areas
?

Click through the PLOS taxonomy to find articles in your field.

For more information about PLOS Subject Areas, click here.

  • Loading metrics

Stability and applicability of the leaf value model for variable nitrogen application based on SPAD value in rice

  • Jie Li,

    Roles Data curation, Formal analysis, Investigation, Methodology, Project administration, Visualization, Writing – original draft

    Affiliations College of Agronomy, Guizhou University, Guiyang, China, Qiandongnan Vocational and Technical College for Nationalities, Kaili, China

  • Yuehua Feng ,

    Roles Conceptualization, Funding acquisition, Methodology, Project administration, Resources, Supervision, Visualization, Writing – review & editing

    fengyuehua2006@126.com

    Affiliations College of Agronomy, Guizhou University, Guiyang, China, Laboratory of Plant Resource Conservation and Germplasm Innovation in Mountainous Region (Ministry of Education), Guizhou University, Guiyang, China

  • Xiaoke Wang,

    Roles Investigation, Project administration

    Affiliation College of Agronomy, Guizhou University, Guiyang, China

  • Jinfeng Peng,

    Roles Investigation

    Affiliation College of Agronomy, Guizhou University, Guiyang, China

  • Dinghua Yang,

    Roles Investigation

    Affiliation Qiandongnan Vocational and Technical College for Nationalities, Kaili, China

  • Guiling Xu,

    Roles Investigation

    Affiliation College of Agronomy, Guizhou University, Guiyang, China

  • Qiangxin Luo,

    Roles Project administration

    Affiliation College of Agronomy, Guizhou University, Guiyang, China

  • Lingli Wang,

    Roles Investigation

    Affiliation College of Agronomy, Guizhou University, Guiyang, China

  • Da Ou,

    Roles Investigation

    Affiliation College of Agronomy, Guizhou University, Guiyang, China

  • Wei Su

    Roles Investigation

    Affiliation College of Agronomy, Guizhou University, Guiyang, China

Abstract

Many fertilization models have been created to scientifically determine the amount of fertilization. With the same purpose, we constructed a nitrogen (N) application model, the leaf value model, which can make N fertilizer decisions in a timely, fast and nondestructive manner during rice planting. However, only one area (A1, Jiuzhou Town, Xixiu District, Guizhou Province) and one cultivar (Qyou6) were involved in the construction of the leaf value model. Its stability and applicability could not be well evaluated. Thus, we chose another area (A2, Jiuzhou Town, Huangping County, Guizhou Province) in Guizhou Province and carried out the experiment by using four cultivars (Nie5you5399, Qyou6, Yixiangyou2115 and Zhongzheyou8) for the leaf value model construction. Compared with the average value of apparent total N uptake (Nz) obtained in 2 years in the A1 area, that in the Qyou6 leaf value model in the A2 area increased by 12%, reaching 635.72 kg ha−1, whereas the corresponding target yield changed slightly, reaching 10,999.90 kg ha−1. Simultaneously, the linear relationship between several good SPAD value-derived indexes (Ys) and apparent N supply of the field (Nx) was still significant or extremely significant in the Qyou6 leaf value model. Compared with the A1 area, it slightly differed, and the R2 of SPADL1 was higher than that of SPADL3×L4/mean. In the leaf value model of the other three cultivars, the relationship between yield and Nx and that between Ys and Nx were significant or extremely significant. The Nz of Yixiangyou2115 and Zhongzheyou8 (618.33 and 617.76 kg ha-1) were close to that of Qyou6 and the corresponding target yields were 10313.36 and 10301.99 kg ha-1, respectively. The Nz and target yield of Nie5you5399 were lowest at 546.63 and 10680.24 kg ha-1, respectively. In general, this study showed that relationships used in the construction of leaf value model had certain stability and applicability to difference areas and cultivars. The leaf value model can be considered in N fertilizer decision-making of rice planting management.

Introduction

Nitrogen (N) application is one of the most important artificial measures in crop planting management and is closely related to yield and environment. If the N application rate is too low to meet the needs of rice plant growth, then the ideal yield cannot be obtained; too much N application rate will lead to waste of resources and environmental pollution [12]. For this reason, researchers have constructed dozens of fertilization models. The construction of fertilization models is based on the diagnosis of nutrition and demand for yield. Nutrition diagnoses include the determination of soil nutrients and the analysis of plant nutrition indexes [3]. Previous research focused on the analysis of soil nutrients, with development of instruments, and later period focused on the analysis of plant nutrition indexes. In plant nutrition diagnosis, especially N nutrition diagnosis, physical diagnosis methods gradually replaced the traditional chemical diagnosis methods, such as the application of hyperspectral equipment, SPAD meter, Greenseeker and digital camera [48]. Physical diagnosis methods, such as the SPAD meter application, could evaluate the N nutrition in a timely, quick and non-destructive manner [9]. The SPAD value represents the relative chlorophyll content, which can be related to photosynthesis and N nutrition status. Thus, it is widely used in plant research.

Many reports on N diagnosis based on the SPAD value are available. These reports focus on the determination site, determination period and relationship with N nutrition [10]. In determining the N application rate based on the SPAD value, the threshold model is used. If the actual measured SPAD value exceeds the set threshold SPAD value, fertilization or a small amount of fertilizer will be applied; otherwise, fertilization will not be applied [1113]. The threshold model cannot well quantify the amount of N fertilizer, which is not easy to obtain. However, some fertilization models, such as the target yield model, involve some parameters, e.g. soil N supply and fertilizer recovery efficiency, which are variable and difficult to calibrate [14]. Therefore, in 2015–2016, based on a large number of experiments, we constructed an N application model, leaf value model, which had a linear relationship between the SPAD value and N application rate,and could guide N topdressing in rice [15]. The leaf value model can be used to estimate the N content of soil quickly and simply by the SPAD meter, and it also has the intuitiveness of the fertilizer effect function model.

At that time of first construction, however, only one area (A1, Jiuzhou Town, Xixiu District, Guizhou Province) and one cultivar (Qyou6) were involved in the construction of the leaf value model. We could not determine whether the relationship involved was tenable in another area and other cultivars. We could not further evaluate the model’s stability and applicability. Therefore, in this study, we chose another experiment site (A2, Jiuzhou Town, Huangping County, Guizhou Province) and constructed the leaf value model for four cultivars to test whether the relationship in model are still significant or extremely significant. Meanwhile, we carried out the variable N application experiment based on the leaf value model at the former test site to reinvestigate the variable N application effect of the leaf value model in that place.

Materials and methods

Experiment site description

Field experiment was conducted in Jiuzhou town (A1), Xixiu District, Guizhou Province, China (26°10′58.8″–26°17′56.4″N, 106°2′42″–106°11′42″E) and Jiuzhou town (A2), Huangping County, Guizhou Province, China (26°54′54″–27°5′49.2″N, 107°38′45.6″–107°51′39.6″E). The climate of A1 and A2 in the field test area of this study belongs to the north subtropical monsoon humid climate. The altitude, annual average temperature, a frost-free period and annual rainfall in A1 and A2 are as follows: 1275 m, 14.3°C, 269.9 days and 1360 mm; and 698 m, 15.7°C, 268 days and 1200 mm, respectively. In 2015–2016, we conducted experiments in A1 and constructed the leaf value model for the first time. One paddy field (F1) was used in in A1, and five paddy fields (F2, F3, F4, F5 and F6) were used in A2. The environmental differences between the two areas were obvious, and both areas have large paddy fields in Guizhou Province.

Experiment design

In F1 of A1, random block design was used with three treatments, blank control, average N fertilization and variable N fertilization. The blank control did not apply N fertilizer, and three plots were available. N fertilizer was split-applied at a rate of 150 kg N ha−1 as 35% basal and 20% at early-tillering in variable N fertilization. When fertilization was performed at the panicle fertilizer stage, the leaf value model was used to determine the top-dressing fertilizer rate. Five plots existed in variable N fertilization. Different from variable N fertilization, the N application rate of panicle fertilizer in average N fertilization was the average rate of the N fertilizer of five plots in variable N fertilization. Four plots were observed in average N fertilization. All plots were 12 m2 in F6.

In each paddy field of A2, the experiments were in a split-plot design with N application rates as the main plots and cultivars as the subplots. The N application rate was N0 (0 kg ha−1), N1 (75 kg ha−1), N2 (150 kg ha−1) and N3 (225 kg ha−1), and the rice cultivars included Nie5you5399, Qyou6, Yixiangyou2115 and Zhongzheyou8. The N fertilizer was split-applied in treatments with 35% as basal, 20% at early tillering (7 days after transplanting), 30% at flower-promoting fertilizer and 15% at flower fertilizer. The sizes of the main plot were 37.38, 34.20, 43.80, 33.60 and 34.20 m2 in F2, F3, F4, F5 and F6, respectively. The experiment was replicated twice in F2 and thrice in other fields. However, in the course of the experiment, all field plots were only repeated twice due to the relatively large workload.

In addition, Phosphorus (P) and potassium (K) fertilizers were applied in the same manner in all treatments. P fertilizer was applied as basal at a rate of 96 kg P2O5 ha−1, and K fertilizer was split-applied at a rate of 180 kg K2O ha−1 as 50% basal and 50% booting. Urea, superphosphate and potassium chloride were used as N, P and K fertilizers, respectively.

Soil and crop management

Before ploughing, soil samples were taken from the upper 20 cm of the soil to determine the chemical properties of six paddy fields. The soil chemical properties are listed in Table 1. Seeds were sown on late April and transplanting were conducted on late May. Before transplanting, the basal fertilizer was applied according to the fertilization design. When transplanting, the rice seedlings in A1 and A2 were transplanted at a hill spacing of 16.7 cm × 30 cm and 20 cm × 30 cm, with one seedling per hill. 7 days after transplanting, tillering fertilizer was applied. And panicle fertilizer was applied according to the fertilization design. It was fine management of water and fertilizer during rice growth, and insects were intensively controlled by chemicals to avoid yield loss.

thumbnail
Table 1. Soil chemical properties of the 0–20 cm layer of the experimental plots.

https://doi.org/10.1371/journal.pone.0233735.t001

Plant sampling

At maturity, according to the standard of the average tiller number (10 points), the plants were sampled from six hills in each subplot. The plant samples were separated into straw (including rachis) and filled and unfilled spikelets and oven-dried at 80°C to constant weight. The total biomass production was the summation of straw, filled and unfilled spikelet dry matter. Then, the N content of straw (including rachis), filled and unfilled spikelets was determined.

Introduction of leaf value model

According to the previous studies (Li et al., 2017), a general expression for leaf value model is (1)

Where:

Nw, total N application rate (kg ha−1);

Nz, apparent total N uptake (kg ha−1);

Ys, SPAD value derivative index;

k and b, the slope and origin of the linear relationship between the apparent N supply of field (Nx) and Ys;

Ng, N application rate before topdressing (kg ha−1).

Nx is the sum of the Nt (kg ha−1) and Ng. Nt is the product of alkali-hydrol N and soil density (2.25×106 kg ha−1). Nz is calculated by the curve fitting relationship between the yield and Nx and is the Nx that reaches the target yield at heading. Ys was obtained by processing the SPAD value of the rice canopy leaves through certain mathematical relations. In addition to the SPAD value of a single leaf, three Yss were used in this study, that is, SPADaverage (average SPAD value of the rice canopy leaves), SPADL2×L3/average (product of SPAD value of the second leaf and SPAD value of the third leaf divided by the average SPAD value of the rice canopy leaves), SPADL3×L4/average (product of the SPAD value of the third leaf and the SPAD value of the fourth leaf divided by the average SPAD value of the rice canopy leaves).

Measurements and calculations

SPAD measurements.

At jointing, five main stems of rice were randomly selected from each plot, and the SPAD value of the three fully expanded leaves at the top was determined by the SPAD meter (SPAD-502, Minolta Camera Co., Osaka, Japan), whereas the SPAD value of the four leaves at heading was determined. From the top down, the heart leaf of the main stem with more 80% extension was divided as the top 1 leaf (L1). When the heart leaf did not reach 80% extension, the leaf as its next leaf is divided as L1. During measurement, the narrow side on both sides of the main vein was selected for the measurement, and it was 1/2 of the width of narrow side of the leaf. Three measurement sites were available on the narrow side of a leaf, i.e. 1/2 of the length of the leaf and 3 cm up and down.

Yield measurements.

In A1, The grain yield was determined from 90 hills in each plot and adjusted to the standard moisture content of 0.135 g H2O g−1 at maturity. In A2, as the yield data, the grain yield was determined from 40 hills in each subplot and adjusted to the standard moisture content of 0.135 g H2O g−1.

N application rate of variable N fertilization and average N fertilization in A1.

At jointing, SPADL3 was measured as Ys. Then, the total N application rate of each plot was calculated through the leaf value model determined by the test data of 2015–2016: (2)

The results showed that the final N application rate of the five plots in F1 was 153.12, 175.09, 227.11, 238.48 and 182.99 kg ha−1. And the N application rate of each plot in average N fertilization was 195.36 kg N ha−1.

N use efficiency.

N agronomic efficiency (kg grain kg−1 N) = (grain yield in N fertilization condition—grain yield in N-unfertilized condition)/N application rate, N recovery efficiency (%) = 100 × (total aboveground plant N accumulation in the plot received N fertilizer—total aboveground plant N accumulation in the zero-N control)/N application rate, N partial productivity (kg grain kg−1 N) = grain yield in N fertilization condition/N application rate, N contribution rate (%) = 100 × (grain yield in N fertilization condition—grain yield in N-unfertilized condition)/grain yield in N fertilization condition.

Harvest index.

100 × filled spikelet weight/total biomass.

Statistical analysis

Data were organized and analyzed by Microsoft Excel 2010. The SPAD value data in all plots was the average value of 15 data (3 sites on one leaf, 5 leaves in total). The coefficient of determination (R2) was used to evaluate the relationship between yield and Nx, and between Ys and Nx by regression analysis (P = 0.01, P = 0.05). The maximum value, minimum value, average value, standard deviation (SD) and coefficient of variation (CV) were cited to compare and evaluate the effect of variable N fertilization and average N fertilization.

Results

Relationship between yield and Nx

In A2 area, an extremely significant curve fitting relationship was observed between the yield of four rice cultivars and Nx (Fig 1). According to this relationship, the target yields of Nie5you5399, Qyou6, Yixiangyou2115 and Zhongzeyou8 were 10680.24, 10999.90, 10313.36 and 10301.99 kg ha−1, respectively. The corresponding values of Nx, that was their Nz values, were 546.63, 635.72, 618.33 and 617.76 kg ha−1, respectively. The target yield and Nz of Yixiangyou2115 were close to those of Zhongzeyou8. Nz of Nie5you5399 was 14.01%, 11.60% and 11.51% lower than that of other cultivars. Whereas its target yield was only lower than those of Qyou6 about 2.91%, and was 3.56% and 3.67% higher than those of Yixiangyou2115 and Zhongzeyou8, respectively. In general, the target yields of the four cultivars were almost the same.

thumbnail
Fig 1. Curve fitting between the apparent N supply (Nx) and yield at heading.

** indicates extremely significant correlation at 0.01 probability level. Nx is the sum of the Nt (kg ha−1) and Ng. Nt is the product of alkali-hydrol N and soil density (2.25×106 kg ha−1).

https://doi.org/10.1371/journal.pone.0233735.g001

Relationship between Ys and Nx

Similarly, in A2 area, whether at jointing or heading, significant or extremely significant relationships were observed between several Yss of the four rice cultivars and Nx (Table 2). Yixiangyou8 all reached an extremely significant relationship at jointing and heading. In addition, R2 of SPADL1 and SPADL2 was higher at jointing than other Yss, whereas it was SPADL3 and SPADaverage at heading. As can be seen from the slope and origin of the linear fitting equation, slope at heading was generally higher than at jointing, and the origin was smaller at heading than at jointing. On the whole, R2 at heading was higher at heading than at jointing. These results were consistent with the experiment results of A1 in 2015 and 2016.

thumbnail
Table 2. Linear fitting relationship between the apparent N supply (Nx) and SPAD value derivative index (Ys) at the jointing and heading in rice.

https://doi.org/10.1371/journal.pone.0233735.t002

Mathematical expression of the leaf value model for the four rice cultivars

According to the relationship between Nx and yield, Nx and SPADL3, the leaf value models of the four rice cultivars were constructed (Table 3). When SPADL3 increased by 1, the top-dressing rates of Nie5you5399, Qyou6, Yixiangyou2115 and Zhongzeyou8 would decrease by 71.94, 67.57, 60.61 and 48.78 kg ha−1 at jointing and by 34.72, 48.78, 40.49 and 31.15 kg ha−1 at heading, respectively. Moreover, when SPADL3 was 46.20, 46.36, 46.07 and 48.43, respectively, N fertilizer could not be applied at jointing in the four rice cultivars.

thumbnail
Table 3. General mathematical expression of the leaf value model at the jointing and heading in the four rice cultivars.

https://doi.org/10.1371/journal.pone.0233735.t003

Effect of variable N application based on the leaf value model

As illustrated in Table 4, The yield of variable N fertilization was 4.39% higher than that of average N fertilization, and the CV of yield of variable N fertilization was significantly lower than that of average N fertilization. In harvest index, variable N fertilization was slightly higher than average N fertilization. Among N use efficiency, the N agronomic efficiency, N partial productivity and N contribution rate of variable N fertilization were 26.98%, 7.52% and 21.05% higher than that of average N fertilization, respectively. In general, when the leaf value model was used to guide the implementation of variable N fertilization, the differences in yield, harvest index, N agronomic efficiency, N recovery efficiency and N contribution rate between area units could be reduced. Meanwhile, which could improve the yield, N agronomic efficiency, N partial productivity and N contribution rate.

thumbnail
Table 4. Yield, harvest index and N use efficiency under the leaf value model based on the variable rate fertilization in Qyou6.

https://doi.org/10.1371/journal.pone.0233735.t004

Discussions

Determination of Nz

The relationship between the Nx and yield is an important basis in the construction of the leaf value model. Combined with the previous research, we can know that the relationship between the Nx and yield is extremely significant in the leaf value model of different cultivars in different areas and years.

In addition, we paid considerable attention to the change of Nz of different cultivars in different areas and years. Previous studies in A1 area showed that Nz of Qyou6 in 2015 and 2016 was 575.27 and 546.71 kg ha−1, respectively, with a difference of only 28.56 kg ha−1. The corresponding target yield was 9264.93 and 11167.97 kg ha−1, with a difference of 1903.04 kg ha−1 [15]. In this study, Nz and yield of Qyou6 were 635.72 and 10999.90 kg ha−1, respectively, which had approximately 13.32% and 7.67% higher than the average of 2 years in the A1 area. Yield is affected by cultivar and environment, and soil N and artificial N supply are only one of the factors. Therefore, in previous studies, we tend to assume that Nz for cultivars to obtain its target yield was fixed, whereas the yield would be different due to different environmental factors, such as light and temperature. Previous experiment results also supported this hypothesis. To date, the results in A1 and A2 areas showed a 13.32% change in Nz, whereas this did not affect our hypothesis. To make the application of leaf value model more convenient, we still think that this hypothesis is tenable and that Nz can be determined in 500–550 kg ha-1. According to the 3-year experiment data of Qyou6 and experiment data of the three other cultivars, the average value of Nz was 590.07 kg ha−1, and the corresponding yield was 10,454.73 kg ha−1. When Nz was determined to be 550 and 500 kg ha−1, the corresponding yield was 10390.40 and 10233.67 kg ha−1. Nz decreased by 6.79% and 15.26%, but its corresponding yield only decreased by 0.62% and 2.11%, with minimal change.

Selection of Ys

This study also showed that the relationship between Nx and Ys also had stability and applicability. In the four rice cultivars, significant or extremely significant relationship between Nx and many Yss at jointing and heading was observed, and the correlation at heading was still higher than at jointing. In combination with previous studies [15], under the same Ys, the slope and origin of the linear relationship between Nx and Ys had certain changes in different areas, years and cultivars, but these changes were not very large. Moreover, in practical application, k and b could be obtained and corrected by only two treatments (with known N application rate). Here, we would like to highlight the selection of Ys. In previous studies [15], the L1 we selected was semi-expanded at jointing. In this study, we changed this approach. As studied by Li [16], the heart leaf of the main stem with more than 80% extension was divided as L1, or when the heart leaf does not reach 80% extension, the next leaf was divided as L1. In addition, in this study, SPADL1 had higher R2 than other Yss, which was different from our previous studies. Some studies showed that the lower leaves were more effective than the upper leaves in diagnosis of rice plant N nutrition or in response to soil N supply [1619]. In studies on rice [20], wheat [21,22], corn [23], cotton [24], barley [25] and potato [26], the top 1 leaf was also used for the measurement of the SPAD value. Here, as our previous research point of view, when the leaf value model was applied, Ys was selected depending on the different purposes of users, such as convenience of measurement, size of correlation and N application rate [15]. Now, we prefer to use SPADL3.

Appraisals of variable N application based on the leaf value model

According to the spatial variation of soil nutrients, precision agriculture divides the planting area into several small plots and determines the fertilization rate according to the soil N of each plot, thereby maximising the utilisation of the cultivated land and fertilizer resources [27]. Variable rate fertilization technology is the core of precision agriculture, and to a certain extent, the fertilization model determines the application effect of variable rate fertilization technology. The effect of variable N application based on the leaf value model showed that this model can evaluate Nt by measuring the SPAD value to obtain Ys and can determine the N application rate according to the soil nutrient situation of the operating unit. Increasing the yield could also reduce the yield difference and increase the N agronomic efficiency, N partial productivity and N contribution rate.

The target yield fertilization model [28], the fertilizer effect function model [29] and the N nutrition index [30,31] are commonly used models that guide field crop N application. Compared with these models, the use of the leaf value model in the determination of Nz was similar to that of the fertilizer effect function model. The leaf value model also considers the soil nutrient status and other factors similar to those considered by the target yield model. Different from the diagnosis of plant N nutrition, the leaf value model is used to estimate the Nx by Ys. Our study showed that the leaf value model had stability and applicability.

In the construction and design of the leaf value model, the leaf value model at jointing was used for N topdressing in the season, and the leaf value model at heading was used to determine the base fertilizer and tiller fertilizer in the next season. For the application of the leaf value model at heading, further research is needed to show its application effect.

Conclusions

The leaf value model, a linear relationship N fertilization model based on SPAD value constructed in 2015–2016, is a simple, quick and dynamic N decision model. For the parameters of the leaf value model and the relationship constructed, the results of this study show that there is extremely significant curve fitting relationship between Nx and yield, and significant or extremely significant linear fitting relationship is observed between Nx and Ys. Combined with the data of three years (2015, 2016, 2019), leaf value model has certain stability and applicability in different areas, years and cultivars. This N application model is worth considering and recommending in rice planting N decision-making. In addition, considering the change of Ys effect, Nz and leaf value model coefficient, extensive research should be carried out under different environmental conditions.

References

  1. 1. Galloway JN, Townsend AR, Eeisman JW, Bekunda M, Cai Z, Freney JR, et al. Transformation of the nitrogen cycle: Recent trends, questions, and potential solutions. Science. 2008; 320: 889–892. pmid:18487183
  2. 2. Dong C, Zhao G, Su B, Chen X, Zhang S. Decision model of variable nitrogen fertilizer in winter wheat returning green stage based on UAV multi-spectral images. Spectrosc Spect Anal. 2019; 39: 3599–3605.
  3. 3. Bai Y. Review on research in plant nutrition and fertilizers. Sci Agric Sin. 2015; 48: 3477–3492.
  4. 4. Chen P, Haboudane D, Tremblay N, Wang J, Vigneault P, Li B. New spectral indicator assessing the efficiency of crop nitrogen treatment in corn and wheat. Remote Sens Environ. 2010; 114: 1987–1997.
  5. 5. Bullock DG, Anderson DS. Evaluation of the minolta SPAD-502 chlorophyll meter for nitrogen management in corn. J Plant Nutr. 1998; 21: 741–755.
  6. 6. Liu X, Zhang K, Zhang Z, Cao Q, Lv Z, Yuan Z, et al. Canopy chlorophyll density based index for estimating nitrogen status and predicting grain yield in rice. Front Plant Sci. 2017; 8.
  7. 7. Walsh OS, Klatt AR, Solie JB, Godsey CB, Raun WR. Use of soil moisture data for refined Green Seeker sensor based nitrogen recommendations in winter wheat (Triticum aestivum L.). Precis Agric. 2013; 14: 343–356.
  8. 8. Wang Y, Wang D, Shi P, Omasa K. Estimating rice chlorophyll content and leaf nitrogen concentration with a digital still color camera under natural light. Plant Methods. 2014; 10: 1–11.
  9. 9. Peng S, Garcia FV, Laza RC, Cassman KG. Adjustment for specific leaf weigF improves chlorophyll meter's estimate of rice leaf nitrogen concentration. Agron J. 1993; 85: 987–990.
  10. 10. Li G, Ding Y, Xue L, Wang S. Research progress on diagnosis of nitrogen nutrition and fertilization recommendation for rice by use chlorophyll meter. Plant Nutr Fert Sci. 2005; 11: 412–416.
  11. 11. Huang J, He F, Cui K, Buresh RJ, Xu B, Gong W, et al. Determination of optimal nitrogen rate for rice varieties using a chlorophyll meter. Field Crops Res. 2008; 105: 70–80.
  12. 12. Peng S, Garcia FV, Laza RC, Sanico AL, Visperas RM, Cassman kG. Increased N—use efficiency using a chlorophyll meter on high—yielding irrigated rice. Field Crops Res. 1996; 47: 243–252.
  13. 13. Dobermann A, Witt C, Dawe D, Abdulrachman S, Gines HC, Nagarajan R, et al. Site-specific nutrient management for intensive rice cropping systems in Asia. Field Crops Res. 2002; 74: 37–66.
  14. 14. Huang Z, Mi C, Hou Y, Liu S, Zheng H, Wang N, et al. Comparison of target yield method and ecological balanced fertilization model. J Agric Resou Envi. 2014; 31: 495–499.
  15. 15. Li J, Feng Y, Mou G, Xu G, Luo Q, Luo K, et al. Construction and application effect of the leaf value model based on SPAD value in rice. Sci Agric Sin. 2017; 50: 4714–4724.
  16. 16. Li G, Xue L, You J, Wang S, Ding Y, Wu H, et al. Spatial distribution of leaf N content and SPAD value and determination of the suitable leaf for N diagnosis in rice. Sci Agric Sin. 2007; 40: 1127–1134.
  17. 17. Wang S, Ji Z, Liu S, Ding Y, Cao W. Relationships between balance of nitrogen supply-demand and nitrogen translocation and senescence of leaves at different positions of rice. Sci Agric Sin. 2003; 36: 1261–1265.
  18. 18. Wang S, Zhu Y, Jiang H, Cao W. Positional differences in nitrogen and sugar concentrations of upper leaves relate to plant N status in rice under different N rates. Field Crops Res. 2006; 96: 224–234.
  19. 19. Zhou Q, Wang J. Comparison of upper leaf and lower leaf of rice plants in response to supplemental nitrogen levels. J Plant Nutr. 2003; 26: 607–617.
  20. 20. Balasubramanian A. Adaptation of the chlorophyll meter (SPAD) technology for real-time N management in rice: a review. Int Rice Res Notes. 2000; 25: 4–8.
  21. 21. Reeves D. Determination of wheat nitrogen status with a hand-held chlorophyll meter: Influence of management practices. J Plant Nutr. 1993; 16: 781–796.
  22. 22. Follett RH, Follett RF, Halvorson AD. Use of a chlorophyll meter to evaluate the nitrogen status of dryland winter wheat. Commun Soil Sci Plant Anal. 1992; 23: 687–697.
  23. 23. Wood CW, Reeves DW, Duffield RR, Edmisten KL. Field chlorophyll measurement for evaluation of corn nitrogen status. J Plant Nutr. 1992; 15: 487–500.
  24. 24. Wu F, Wu L, Xu F. Chlorophyll meter to predict nitrogen sidedress requirements for short-season cotton (Gossypium hirsutum L.). Field Crops Res. 1998; 56: 309–314.
  25. 25. Peltonen J, Virtanen A, Haggren E. Using a chlorophyll meter to optimize nitrogen fertilizer application for intensively-managed small cereals. J Agron Crop Sci. 1995; 174: 309–318.
  26. 26. Rodrigues M. Establishment of continuous critical levels for Indices of plant and presidedress soil nitrogen status in the potato crop. Commun Soil Sci Plant Anal. 2004; 35: 2067–2085.
  27. 27. Morala FJ, Terrónb JM, Rebollo FJ. Site-specific management zones based on the Rasch model and geostatistical techniques. Comput Electron Agr. 2011; 75: 223–230.
  28. 28. Miransari M, Mackenzie AF. Development of a soil N test for fertilizer requirement for wheat. J Plant Nutr. 2011; 34: 762–777.
  29. 29. Afzal S, Islamb M, Obaid-Ur Rehmanc. Application of mitscherlich-bray equation for fertilizer use on groundnut. Commun Soil Sci Plant Anal. 2014; 45: 1636–1645.
  30. 30. Wang Y, Shi P, Zhang G, Ran J, Shi W, Wang D. A critical nitrogen dilution curve for japonica rice based on canopy images. Field Crops Res. 2016; 198: 93–100.
  31. 31. Wang Y, Shi P, Ji R, Min J, Shi W, Wang D. Development of a model using the nitrogen nutrition index to estimate in-season rice nitrogen requirement. Field Crops Res. 2020; 245: 107664.