ORIGINAL RESEARCH article

Front. Plant Sci., 04 September 2017

Sec. Plant Breeding

Volume 8 - 2017 | https://doi.org/10.3389/fpls.2017.01517

Development of a Critical Nitrogen Dilution Curve Based on Leaf Area Duration in Wheat

  • 1. National Engineering and Technology Center for Information Agriculture, Key Laboratory for Crop System Analysis and Decision Making, Ministry of Agriculture, Jiangsu Key Laboratory for Information Agriculture, Jiangsu Collaborative Innovation Center for Modern Crop Production, Nanjing Agricultural University Nanjing, China

  • 2. Henan Institute of Science and Technology Xinxiang, China

  • 3. Key Laboratory of Soil Chemistry and Environmental Protection, Institute of Soil Science, Chinese Academy of Sciences Nanjing, China

Abstract

Precise quantification of plant nitrogen (N) nutrition status is essential for crop N management. The concept of critical N concentration (Nc) has been widely used for assessment of plant N status. This study aimed to develop a new winter wheat Nc dilution curve based on leaf area duration (LAD). Four field experiments were performed on different cultivars with different N fertilization modes in the Yangtze River basin and Yellow River basin in China. Results showed that the increase in LAD with increasing cumulative thermal time took the shape of an “S” type curve; whereas shoot N concentration decreased with increasing LAD, according to a power function. Both LAD and shoot N concentration increased with increasing N application. The new LAD based Nc dilution curve was determined and described as Nc = 1.6774 LAD−0.37 when LAD > 0.13. However, when LAD ≤ 0.13, Nc was constant and can be calculated by the equation when LAD = 0.13. The validation of Nc dilution curve with dataset acquired from independent experiments confirmed that N nutrition index (NNI) predictions based on the newly established Nc dilution curve could precisely diagnose N deficiency at different plant growth stages. The integrated N nutrition index (NNIinte), which was obtained by the weighted mean of NNI, was used to estimate shoot N concentration, shoot dry matter, LAD, and yield using regression functions. The linear relationships between NNIinte and these growth variables were well correlated. These results provided enough evidence that the new LAD–based Nc dilution curve could effectively and precisely diagnoses N deficiency in winter wheat crops.

Introduction

Nitrogen (N) is a major nutritional element required for crop growth and productivity (Eickhout et al., 2006). In-season N diagnosis and management in crop production is an important agronomic practice for ensuring food security and environmental sustainability (Ata-Ul-Karim et al., 2016a,b). N is an integral component of chlorophyll and protein related with plant vigor, yield and quality and appropriate N fertilizer application is not only imperative for improving crop production but also for increasing the processing quality of grain in cereals (Wang et al., 2002; Zhu et al., 2003). Previous studies have shown a strong relationship between N status and wheat yield and quality. Excessive N fertilizer application has negative impacts, such as decreased N use efficiency and deterioration of the ecological environment (Jeuffroy and Meynard, 1997; Zhao et al., 2016a). Therefore, accurate estimation of crop N status and precision N management has become pressing study topics for crop growers and scientists.

Reliable estimation of crop N requirement and diagnosing nutritional status at critical crop growth stages is imperative for optimizing the quantitative and qualitative aspects of crop production (Ata-Ul-Karim et al., 2016b). Non-destructive diagnostic tools such as chlorophyll meter, digital camera, active and passive sensors (Mistele and Schmidhalter, 2008; Lee and Lee, 2013; Yuan et al., 2016) have been widely used to assess the crop N status. Besides non-destructive diagnostic tools, the concept of critical N (Nc), which is crop specific, simple and biologically sound diagnostic tool, has also widely been used as a classical destructive diagnostic tool for diagnosing crop N status. The Nc dilution curves based on shoot biomass have been established in various agronomic crops, including corn (Plénet and Lemaire, 1999; Ziadi et al., 2008; Yue et al., 2014), rice (Sheehy et al., 1998; Ata-Ul-Karim et al., 2013; Huang et al., 2015; He et al., 2017), wheat (Justes et al., 1994; Ziadi et al., 2010; Yue et al., 2012; Zhao et al., 2012). Attempts have also been made to establish the Nc dilution curves on organ specific dry matter (leaf dry matter and panicle dry matter) in wheat (Yao et al., 2014a; Zhao et al., 2016b), on leaf dry matter and stem dry matter in rice ecotypes (Ata-Ul-Karim et al., 2014a; Yao et al., 2014b), and on leaf dry matter in maize (Zhao et al., 2017). The existing dilution curves based on whole plant and organ specific dry matter were developed for crops grown in a specific region, and the model parameters vary with climate. Additionally, dry matter is generally measured by destructive sampling. Leaf area index (LAI) and canopy coverage is a fundamental variable in agronomic and environmental research have been widely used as reference plant index for monitoring crop growth, predicting grain yield, and optimization of crop management practices (Ata-Ul-Karim et al., 2014b; Wang et al., 2016). With the advances in sensor technologies, it is convenient to measure LAI non-destructively. LAI measurements are more easily performed than quantification of shoot dry matter, providing an alternative method for evaluation of plant N status diagnosis. Previous studies have developed Nc dilution curves based on LAI, providing an alternative method for plant Nc assessment and N status diagnosis (Confalonieri et al., 2011; Ata-Ul-Karim et al., 2014b; Zhao B. et al., 2014). Crop dry matter accumulation occurs through photosynthetic production, which is depends heavily on the canopy leaves characteristics. Leaf area, its duration, and photosynthetic rate are the key factors driving matter accumulation and yield (Gastal and Bélanger, 1993; Hammer and Wright, 1994; Zhao et al., 2005). Previous studies have shown that N allocation to leaves for growth and photosynthesis to be prioritized within wheat plants (Jeuffroy et al., 2002). N nutritional status, leaf area, and leaf area duration affect both short-term growth status and future productivity. LAI is an instantaneous estimate of the photosynthetic area. However, LAI does not represent photosynthesis as a continuous process occurring within the photosynthetic area and does not precisely explain the formation of dry matter. Leaf area duration (LAD) is the integral of LAI and growth period, and comprehensively incorporates the size and duration of leaf area (Peltonen-Sainio et al., 2008). Therefore, compared with LAI, LAD better represents the process of dry matter accumulation and incorporates plant N dilution mechanisms.

The N diagnostic tools such as N nutrition index (NNI), accumulated N deficit (AND) and N requirement (NR) could be derived from the Nc dilution curve for quantifying the N status, developmental processes and radiation-use efficiency of the plant in response to N supply (Mills et al., 2009; Ata-Ul-Karim et al., 2017a,b) as well as for making decisions on N application (Lemaire et al., 2008). NNI is one of the most widely used Nc curve based N diagnostic tool for in-season plant N status diagnosis (Zhao et al., 2017). The relationships between NNI and relative yield have been previously implicated in assessing yield potential in wheat, corn (Ziadi et al., 2008, 2010), and rice (Ata-Ul-Karim et al., 2016a). Moreover, these relationships have also been used to estimate crop NR for a corrective N fertilization during crop growth period and for assessment of grain amylose and protein content in rice (Ata-Ul-Karim et al., 2017a,c). This concept can successfully differentiate the sub-optimal and supra-optimal N growth conditions in crop production, and the diagnostic tools derived from these curves could be integrated with crop simulation models to assist crop N management (Lemaire et al., 2008; Ata-Ul-Karim et al., 2017b).

The objectives of this study were to develop a Nc dilution curve based on LAD for winter wheat which combines the merits of shoot dry matter- and LAI-based Nc curves and to test the curve through evaluation of the effects of N deficiency on growth variables and yield. The projected results would provide an alternative diagnostic tool for assessing N deficiency limitations to winter wheat growth.

Materials and methods

Experiments design

Four experiments at two sites using different cultivars, N levels, split applications, and top-dressing stages were conducted, as summarized in Table 1. Experiment 4 is the historical data obtained from the study of Zhao B. et al. (2014). Urea used as the N fertilizer was applied before sowing (as base fertilizer) and from the jointing to booting stages (as top-dressing fertilizer), respectively. The experimental fields were plowed and subsequently harrowed before sowing. Pre-emergence herbicides were used to control weeds at early growth stages. Plots were also regularly hand-weeded until canopy was closed to prevent weed damage. Insecticides were used to prevent insect damage when necessary. No irrigation was in Experiments 1, 2 and 4. Irrigations (900 m3 ha−1 per time) at regrowth and jointing were performed in Experiment 3. All other agronomic management practices were used according to local management practices recommendations to ensure maximum potential productivity, i.e., no factor other than N was limiting.

Table 1

ExperimentSeasonSoil characteristicsCultivarN level (Kg ha−1)N Topdressing stage and ratioSampling stages
Exp. 1
(Yizheng, 32°16′ N, 119°10′ E)
2010-2011Type: clay soil
Organic matter: 14.26 g kg−1
Total N: 1.28 g kg−1
Available P: 45.46 mg kg−1
Available K: 87.23 mg kg−1
Yangmai16
(YM16)
0 (N0)
225 (N1)
300 N2
Top 3rd leaf emergence;70%(R1)
Top 3rd leaf emergence;60%(R2)
Top 3rd leaf emergence;50%(R3)
Top 3rd leaf emergence;40%(R4)
Top 3rd leaf emergence;30%(R5)
Spring
re-growth
Jointing
Booting
Heading
Filling
Exp. 2
(Yizheng, 32°16′ N, 119°10′ E)
2010-2011Type: clay soil
Organic matter: 13.54 g kg−1
Total N: 1.2 g kg–1 g kg−1
Available P: 43.01 mg kg−1
Available K: 83.25 mg kg−1
Yangmai16
(YM16)
0(N0)
225(N1)
300(N2);
Top 3rd leaf emergence(T1); 60%(R1)
Top 3rd leaf emergence(T1); 50%(R2)
Top 3rd leaf emergence(T1); 40%(R3)
Top 2rd leaf emergence(T2); 60%(R1)
Top 2rd leaf emergence(T2); 50%(R2)
Top 2rd leaf emergence(T2); 40%(R3) Flag leaf emergence(T3); 60%(R1) Flag leaf emergence(T3); 50%(R2) Flag leaf emergence(T3); 40%(R3)
Spring
re-growth
Jointing
Booting
Heading
Filling
Exp. 3 (Xinxiang, 35°11′ N, 113°48′ E)2011-2012Type: loam soil
Organic matter: 15.7 g kg−1
Total N: 1.45 g kg–1 g kg−1
Available P: 67.54 mg kg−1
Available K: 87.43 mg kg−1
Aikang58
(AK58)
0(N0)
75(N1)
150 N2
225(N3)
300 (N4)
Top 3rd leaf emergence;65%(R1)
Top 3rd leaf emergence;50%(R2)
Top 3rd leaf emergence;35%(R3)
Spring
re-growth
Jointing
Booting
Heading
Filling
Exp. 4
(Yizheng, 32°16′ N, 119°10′ E)
2010-2011Type: clay soil
Organic matter: 13.5 g kg−1
Total N: 1.1 g kg–1 g kg−1
Available P: 43 mg kg−1
Available K: 82 mg kg−1
Ningmai13
(NM13)
0(N0)
75(N1)
150 N2
225(N3)
300 (N4)
375 (N5)
Top 3rd leaf emergence;50%Spring
re-growth
Jointing
Booting
Heading

Basic information of the four experiments.

Experiment 1 followed a split plot design, with N level as the main plot factor and N splits as the subplot factor. The study was conducted in Yizheng, China (32°16′ N, 119°10′ E) during the 2010–2011 wheat growing season. Cultivar Yangmai 16 (YM16, grain protein 14.2%) was planted on October 18th by row seeding, with a row spacing of 25 cm. Planting density was 2.40 × 106 plants ha−1. P (135 kg ha−1 P2O5) and K (190 kg ha−1 K2O) fertilizers were incorporated into the soil as monocalcium phosphate (Ca(H2PO4)2) and potassium chloride (KCl) before planting.

Experiment 2 applied a split-split plot design with N level as main plot factor, N split applications as the subplot factor and N application stage as the sub-sub plot factor. The experiment was conducted in Yizheng, China (32°16′ N, 119°10′ E) during the 2010–2011 wheat season. Cultivar Yangmai 16 (YM16, grain protein 14.2%) was planted on November 7th by row seeding, with a row spacing of 25 cm. The area of each plot was 24 m2. Planting density was 2.40 × 106 plants ha−1. Phosphatic (120 kg ha−1 P2O5) and potassium (150 kg ha−1 K2O) fertilizers were incorporated into the soil as monocalcium phosphate (Ca(H2PO4)2) and potassium chloride (KCl) before planting.

Experiment 3 was a split plot design with N level as the main plot factor, and N split applications as the subplot factor. The experiment was conducted in Xinxiang, China (35°11′ N, 113°48′ E) during the 2011–2012 wheat growing season. Cultivar Aikang 58 (AK58, grain protein 14.48%) was planted on November 7th by row seeding, with a row spacing of 25 cm. The area of each plot was 24 m2. Planting density was 3.0 × 106 plants ha−1. Phosphatic (150 kg ha−1 P2O5) and potassium (150 kg ha−1 K2O) fertilizers were incorporated into the soil as monocalcium phosphate (Ca(H2PO4)2) and potassium chloride (KCl) before planting.

Experiment 4 was conducted involving N levels in Yizheng, China (32°16′ N, 119°10′ E) during the wheat growing season of 2010–2011. Cultivar Ningmai 13 (NM13, Grain protein 10.9%) was planted on 11th Nov. by row seeding, and row spacing was 25 cm. Random plot design was adopted with three replications. Six N levels including 0 (N0), 75(N1), 150(N2), 225(N3), 300(N4), and 375(N5) kg N ha−1, were tested, with 50% of the N applied at pre planting and 50% at jointing. The area of each plot was 30 m2. Planting density was 2.40 × 106 plants ha−1. Phosphatic (96 kg ha−1 P2O5) and potassium (120 kg ha−1 K2O) fertilizers were incorporated into the soil before planting as monocalcium phosphate (Ca(H2PO4)2) and potassium chloride (KCl).

Plant sampling, N determination, and daily temperature

Twenty representative plants were sampled from each plot at regrowth, jointing, booting, heading, anthesis, grain filling, and maturity for determination of LAI, shoot dry matter, and N concentration. For shoot dry matter measurements, samples were dried in a forced-draft oven at 80°C until constant weight was achieved. Shoot N was determined following the micro-Kjeldahl method (Bremner and Mulvancy, 1982). LAI was measured using an LI-3000 (Li-COR, Lincoln, US). Daily temperature measurements were acquired from a weather station near the experimental plots.

Data processing

Statistical analysis

The data for determination of critical N points were analyzed according to the methodology proposed by Justes et al. (1994). For each sampling date, experiment and wheat cultivar, the amounts of LAD and plant N concentration were subjected to analysis of variance (ANOVA) using GLM procedures in IBM SPSS Version19.0 (IBM Corporation, Armonk, New York). Significant effects of treatment on mean LAD were tested for using the least significant difference (LSD0.05) test, and the results were used to classify treatments as N-limiting or non-N-limiting. An allometric function was used to determine the relationship between the observed decreases in N concentration with increasing LAD using Microsoft Excel 2010 (Microsoft Corporation, Redmond, WA, USA).

Relative growing degree days

In order to account for variation in thermal time among different years and locations, relative growing degree days (RGDD, d°C) was expressed as the ratio of growing degree days (GDD, above 0°, d°C) at a sampling time point to the total GDD across the entire growth period (TGDD, d°C).

To perform real-time measurements, the average TGDD from previous years or a real-time TGDD estimated from real-time weather and historic weather data could be substituted for in-season TGDD (Bannayan and Hoogenboom, 2008).

Leaf area duration

Leaf area duration (LAD) was calculated using the LAI at each sampling date and the corresponding RGDD according to the equation:

The LAD between two sampling time points was approximated according to the equation:

where t1 and t2 are the sampling times across which RGDD was quantified. LAI1 and LAI2 are the LAI at the sampling time t1 and t2, respectively.

Establishment and validation of critical nitrogen dilution curve

Data from experiment 1 and 3 were used to develop the Nc dilution curve following the computation method of Justes et al. (1994). If at the same measurement date, statistical analysis distinguished at least one set of N limiting growth and non-N-limiting growth data points, these data were used to define the Nc dilution curve. For each measurement date, the variation in total N concentration and LAD was assimilated into a bilinear relation composed of (a) an oblique line showing the increase in joint LAD and N concentration and (b) a vertical line corresponding to the increase in N concentration and no variation in LAD. The Nc dilution curve was validated for N-limiting and non-N-limiting situations within the range for which it was established using the data collected from the independent experiments (Experiment 2 and 4).

Determination of NNI and integrated NNI

The NNI at each sampling date was calculated according to the equation (Justes et al., 1994):

where Nac and Nc, respectively, are the actual N concentration and critical N concentration of the shoot expressed in % dry matter. If NNI = 1, N nutritional status was considered to be optimum, while NNI > 1 and NNI < 1 indicated excess and deficient N nutrition, respectively.

An integrated NNI can be obtained by the weighted mean of NNI (Lemaire et al., 2008) as illustrated as follows.

where NNIinte is the integrated NNI, T represents the duration sampling occurred over (days or GDD), NNIi is the instantaneous NNI value for sampling time point i, and ti is the time interval between sampling time point i and sampling time point i-1.

Results

Leaf area duration and shoot N concentration in wheat under different N treatments

Table 2 showed the increase in LAD and the decreases in shoot N concentration (Nac) with the development stages in winter wheat cultivars. Although the four experiments were conducted under different cultivars, treatments, sites, years, the change trends in LAD and Nac are similar. LAD and Nac increased with increasing N application, the response of LAD and Nac to N treatment varied with N split applications and topdressing time in different experiments. Topdressing time also showed significant effects on LAD and Nac. The LAD ranged from 0.43 to 1.82, 0.46 to 1.824, 0.075 to 2.361, and 0.046 to 1.317 while Nac ranged from 0.67 to 3.7%, 0.83 to 3.87%, 0.95 to 3.35%, and 1.02 to 4.13% in Experiment 1–4, respectively. The LAD of cultivar Aikang58 in Xinxiang (Experiment 3) was higher than that of experiments in Yizheng (Experiments 1, 2 and 4). However, the variation of Nac in different experiments was minor.

Table 2

ExperimentTreatmentGrowth stage
Spring re-growthJointingBootingHeadingFilling
NacLADNacLADNacLADNacLADNacLAD
Exp. 1 (Yizheng, 2010–2011)N02.89c0.043d2.23d0.110d1.14d0.344d0.88d0.525d0.67c0.695d
N1R13.21b0.056c2.57bc0.152c1.57c0.509c1.65a0.886c1.09ab1.339c
N1R23.24b0.061b2.82b0.166c1.82b0.583bc1.46b0.992b1.06b1.497b
N1R33.34b0.058bc2.71b0.169bc1.74bc0.620b1.27c1.032b1.12a1.521b
N1R43.33b0.059b3.06a0.190b1.93b0.633b1.61a1.048b1.05b1.542b
N1R53.58a0.061b3.09a0.202a2.21a0.685b1.44b1.115b1.08ab1.619b
N2R13.64a0.059b2.60bc0.168bc1.88b0.596bc1.68a1.020b1.16a1.538b
N2R23.60a0.058bc2.54c0.189b2.21a0.642b1.53ab1.077b1.09a1.583b
N2R33.65a0.062ab3.17a0.205a2.18a0.696a1.61a1.179a1.05b1.706ab
N2R43.61a0.063a2.96ab0.217a2.01ab0.767a1.57a1.280a1.08ab1.844a
N2R53.70a0.068a3.14a0.223a2.08ab0.775a1.61a1.275a1.10a1.820a
F prob.
LSD
**********
0.1320.0050.2740.0210.2210.0780.1510.0910.0620.172
Exp. 2 (Yizheng, 2010–2011)N03.07c0.046c2.05c0.111c1.05d0.314c0.94c0.485c0.83c0.672c
N1R1T13.71ab0.056b2.58b0.176b1.82b0.607b1.47b1.003b1.16b1.472b
N1R1T23.48b0.045c2.50b0.158b1.79b0.576b1.46b0.944b1.10b1.348
N1R1T33.47b0.067a2.45b0.196a1.49c0.594b1.42b0.922b1.06b1.284b
N1R2T13.93a0.050bc3.06a0.172b2.12a0.650ab1.73a1.059b1.34a1.537b
N1R2T23.63b0.062a2.88a0.198a1.84b0.653ab1.64a1.058b1.29a1.561b
N1R2T33.84a0.057b3.04a0.192a1.54c0.625b1.51b0.979b1.07b1.425b
N2R1T13.81a0.068a2.83a0.196a2.39a0.713a1.55b1.224a1.37a1.813a
N2R1T23.56b0.062a2.89a0.192a1.79b0.673ab1.66a1.104ab1.36a1.573b
N2R1T34.03a0.057b2.79a0.195a1.59c0.639b1.51b1.009b1.27a1.464b
N2R2T13.80a0.067a2.89a0.216a2.13a0.768a1.62a1.297a1.38a1.966a
N2R2T23.89a0.061a2.84a0.205a1.82b0.716a1.65a1.212a1.39a1.854a
N2R2T33.87a0.063a2.96a0.224a1.73b0.732a1.66a1.218a1.27a1.824a
F prob.
LSD
**********
0.2410.0070.2240.0350.2360.0870.0710.1930.1720.284
Exp. 3 (Xinxiang, 2011–2012)N02.86d0.075c2.77c0.143c1.34d0.612c1.38d0.852c0.95c1.303d
N1R12.96c0.092b2.50cd0.175b1.71c0.780b1.47d1.093b1.09b1.657c
N1R22.96c0.097b2.37d0.178b1.47d0.801b1.49cd1.105b0.97c1.682bc
N1R33.01bc0.102b2.85b0.188ab1.63c0.857b1.82b1.201b1.04b1.872b
N2R13.03bc0.084c2.85b0.166b1.97b0.796b1.61c1.175b1.18a1.933b
N2R23.08b0.100b2.55cd0.191ab1.90b0.879b1.83b1.251b1.12a2.010ab
N2R33.24a0.105b3.12a0.202a1.77c0.909ab1.75b1.298ab1.00bc2.097ab
N3R12.92cd0.102b2.85b0.189ab2.22a0.929ab2.06a1.348a1.23a2.240a
N3R23.30a0.105b2.82b0.204a1.98b0.969a2.04a1.398a1.08b2.317a
N3R33.29a0.113a2.91ab0.201a1.98b1.016a1.93b1.452a1.20a2.292a
N4R13.24a0.115a2.65c0.222a2.06ab1.005a2.13a1.421a1.16a2.273a
N4R23.31a0.112a3.06a0.217a2.26a1.063a1.98ab1.523a1.12ab2.415a
N4R33.35a0.115a3.02a0.217a2.32a1.051a2.16a1.509a1.16a2.361a
F prob.
LSD
**********
0.0870.0080.1430.0290.2120.0960.1170.2720.0810.265
Exp. 4 (Yizheng, 2010–2011)N02.71d0.046d1.74c0.066d1.05d0.320d1.02d0.422d
N13.22c0.119c2.03b0.149c1.46bc0.473c1.10d0.614c
N23.80ab0.170b2.10b0.216b1.60b0.739b1.37c0.975b
N33.87ab0.210a2.18b0.277a1.83b0.966a1.75b1.276a
N44.13a0.196a2.60a0.262a2.47a0.962a1.91a1.284a
N53.65b0.228a2.95a0.295a2.51a0.991a2.01a1.317a
F prob.
LSD
********
0.3130.0340.3940.0490.3970.2210.1820.189

Nitrogen concentration at different growth stages and varied N rates.

LAD, leaf area duration; Nac, shoot nitrogen concentration; LSD, least significant difference.

*

F statistic significant at the 0.05 probability level. Data within a column followed by a different letter are significantly different (P < 0.05).

Determination of the LAD-based critical nitrogen dilution curve

Data from experiments 1 and 3 were used to develop the Nc dilution curve following the method of Justes et al. (1994). Growth rate and cultivar did not significantly affect Nc at both study sites, and the curves for two cultivars showed no significant differences. The unified dilution curve was determined as showed in Figure 1. The constant value of Nc (3.57) is calculated by the equation when LAD = 0.13 when LAD ≤ 0.13. When LAD > 0.13, Nc was calculated as follows.

Figure 1

Validation and application of the established critical nitrogen dilution curve

The LAD-based Nc dilution curve determined in this study was validated with data acquired from independent experiments (Experiment 2 and 4) using a standard method (Justes et al., 1994; Ziadi et al., 2008, 2010; Zhao B. et al., 2014). The data points under various N treatments from Experiment 2 and 4 were characterized by sub-optimal N or supra-optimal N growth conditions based on significant (P ≤ 0.001) differences in shoot biomass for each sampling date, site, and year. Treatments were considered N-limiting when shoot biomass increased significantly with increasing N supply, the N treatment is considered as sub-optimal N supply, while under supra-optimal N supply had no significant increase in biomass with increasing N supply (LSD < 0.05). Data points attained from sub-optimal treatments were located almost below the Nc curve while those of supra-optimal N supply was located close to or above the Nc curve (Figure 2). The newly established LAD-based Nc dilution curve well distinguished the sub-optimal and supra-optimal N supply conditions. Therefore, it was concluded that the curve could be used for in-season assessment of wheat N nutrition status. Data points were selected only from non-N-limiting treatments for the determination of upper limit curve (Nmax) while the data points from N-limiting treatments for which have no N application were used to determine lower limit curve (Nmin) (Figure 2).

Figure 2

Determination of plant N deficiency

Two N application levels in Experiment 2 were relatively high (225 and 300 kg ha−1). However, some N split applications still showed N deficiencies at a few growth stages (Figure 3A). Treatments N1R1T1, N1R1T2, N1R1T3, N1R2T3, and N2R1T3 were N deficient at few developmental stages when the base fertilizer level was relatively lower (R1 treatments), or the topdressing was applied later (T3 treatments). In Experiment 4, shoot N concentration under N4, and N5 treatments was significantly higher than shoot Nc. On the other hand, shoot N concentration under N3 was close to the Nc. However, N0, N1, and N2 showed significant N deficiencies due to low levels of N application (Figure 3B). The Nc fell approximately between N3 and N4 levels.

Figure 3

Relationships between NNIinte and crop growth variables

The effects of N deficiencies on growth variables, i.e., dry matter(DM), shoot nitrogen concentration (Nac), LAD and yield was expressed in relative values (DM/DMmax, Nac/Nacmax, LAD/LADmax, and yield/yieldmax) in order to account for variation in pedoclimatic conditions across the different experiments. The Experiment 1 and 3 were used to develop the linear relationship. DM/DMmax represents the ratio between DM accumulated by a crop receiving a limiting N supply and the maximum shoot biomass value measured at the same moment in time, assuming that the maximum biomass (DMmax) was accumulated under optimum N supply. DMmax was calculated as the mean of DM for the group of treatments giving the highest DM value, and not significantly different from each other, at a 5% significance level. The same methodology was applied to calculate DM/DMmax, Nac/Nacmax, LAD/LADmax, and yield/yieldmax. All of the relative growth variables strongly correlated with NNIinte (Figure 4). The data of experiment 2 was used to validate the relationships between the NNIinte and the relative values (DM/DMmax, Nac/Nacmax, LAD/LADmax, and yield/yieldmax). The results showed that RRMSE, agreement index, and R2 ranged from 6.61% to 11.99%, 0.9515 to 0.9829, and 0.7862 to 0.8968 (Table 3), which indicated that the established relationships are robust. The variation in NNIinte explained more than 70% of variability associated with any of these four variables. When NNIinte equaled 0.80, Nac, LAD, DM, and yield decreased by 42.17, 39.79, 30.44, and 29.95%, respectively.

Figure 4

Table 3

EquationRRMSE (%)Agreement indexPr > FR2
Nac/Nacmax = 2.0603 NNIinte −1.08319.830.9753<0.00010.7862
LAD/LADmax = 1.8813 NNIinte −0.935711.940.9672<0.00010.8193
DM/DMmax = 1.4546 NNIinte −0.498811.990.9515<0.00010.7862
Yield/Yieldmax = 1.4835 NNIinte −0.4936.610.9829<0.00010.8968

Validation of the relationships between the relative growth variables and NNIinte.

Nac, shoot nitrogen Concentration; LAD, leaf area duration; DM, shoot dry matter.

Discussion

Crop photosynthetic production capacity is limited by photosynthetic area (which can be quantified through LAI measurements), the duration of the photosynthetic area, and photosynthetic rate. Plant biomass production can be estimated by measuring LAD, which integrated the size and duration of the photosynthetic area (Gastal and Bélanger, 1993; Hammer and Wright, 1994; Zhao et al., 2005). In this study, LAD increased with the increase of N levels while shoot nitrogen concentration decreased with the N levels (Table 2). Similar results for shoot biomass responses to development time and N levels have been reported in previous studies on wheat (Justes et al., 1994; Zhao et al., 2012), rice (He et al., 2017), and maize (Plénet and Lemaire, 1999).

Most wheat Nc dilution curves have been developed based on dry matter accumulation (Justes et al., 1994; Yue et al., 2012; Zhao et al., 2012), while Zhao B. et al. (2014) developed a Nc dilution curve based on LAI in winter wheat (Figure 5). This study developed a new LAD-based Nc dilution curve by using the data collected through experiments within the Yangtze river and Yellow river basins of China. Comparing with dry matter-based Nc dilution curves (Justes et al., 1994; Yue et al., 2012; Zhao et al., 2012), LAD explains the process of dry matter accumulation, which can be effected by the size and duration of photosynthetic area. Therefore, the rationales for these two Nc dilution curves types are the same. To account for variations in thermal conditions among the different regions and seasons, thermal time was taken into account for calculating LAD in our study. Additionally, LAI is more easily obtained than shoot dry matter due to most of the instruments designed for measuring LAI. Compared to the LAI-based Nc curve, LAD-based curve can well describe the N dilution during the entire crop growth period, while LAI-based curve only quantifies the Nc dynamics till heading stage (when LAI reaches to its maximum value) (Ata-Ul-Karim et al., 2014b; Zhao B. et al., 2014). Due to the decline of LAI during post anthesis growth period, it does not take the shape of an “S” type curve, hence cannot be used as an appropriate plant index for the elucidation of plant N dilution. Therefore, it is not appropriate to integrate the Nc curves into the crop growth models, especially those which use the Nc concentration of entire crop growth period for simulation of plant N dynamics, such as CERES-wheat (Ritchie and Otter, 1985) and APSIM-wheat (Zhao Z. et al., 2014). Thus, LAD-based Nc curve combines the merits of LAI-based and dry matter-based Nc curve. The initial N concentration in the Nc dilution curve developed in present study was close to that determined by Zhao et al. (2012); yet, our results of N concentration (for a low to middle grain protein cultivar) showed significant differences in the N concentration of high grain protein cultivars used by Yue et al. (2014) and Justes et al. (1994). High grain protein cultivars tend to have better rooting capabilities, allowing uptake and assimilation of more N. Thus, more cultivars should be tested to validate this methodology.

Figure 5

The NNI at various wheat growth stages was calculated using the LAD-based Nc dilution curve. As N treatment levels were applied with different splits and top-dressing stages, there were multiple NNI for each N treatment level (Figure 3). NNI values, which can be used for real-time diagnosis of plant N nutrition status, change with time. Due to different N application modes and soil N mineralization, soil N supply capacity changes constantly, the value of instantaneous NNI was highly variable (Figure 3) and could not be applied to detect N deficiency throughout the entire growth period (Lemaire et al., 2008). NNIinte incorporates the knowledge of time course of the NNI during this period (Lemaire and Gastal, 1997; Colnenne et al., 2002). The linear correlations between NNIinte and DM, LAD, Nac, and yield, were very well (Figure 4). The slopes for the relationships between NNIinte and DM and yield were lower than those for the relationships between NNIinte and Nac and LAI, which suggests that the effects of NNIinte on Nac and LAD were greater and more sensitive to N deficiency than the effects on DM and yield. These relationships between NNIinte and DM, LAD, Nac, and yield were similar to those reported in a previous study (Colnenne et al., 2002; Lemaire et al., 2008).

Conclusions

This study has developed a new Nc dilution curve based on LAD in winter wheat, which combines the merits of shoot dry matter- and LAI-based Nc curve. The results showed that the allometric relationships between N uptake and LAD were robust across the different conditions and was described as Nc = 1.6774LAD−0.37 (LAD > 0.13), and Nc = 3.57, when LAD ≤ 0.13. NNI derived from the curve established here, effectively identified the limiting and non-limiting N nutrition at different growth stages. The NNIinte, which was based on Nc, was applied to regression functions to estimate Nc, shoot dry matter, actual shoot N content, LAD, and yield measured in independent experiments at different growth periods. These relationships showed strong linear correlations between NNIinte and the growth variables. Our results indicated that the developed Nc curve provide a new procedure and alternate method to assess crop N status for N diagnose and management.

Statements

Author contributions

XW and LT wrote the manuscript; XW analyzed the experiments data; SA, YZ and WC provided advice and edited the manuscript; LT, WC, and YZ planned experiments and XW, TY, and LL performed experiments. All authors read and approved the final manuscript.

Acknowledgments

This work was supported by grants from the National High-Tech Research and Development Program of China (2013AA102404), National Natural Science Foundation of China (31671582), Jiangsu Agriculture Science and Technology Innovation Fund (CX[14]2116), the Priority Academic Program Development of Jiangsu Higher Education Institutions (PAPD) of China, and The 111 Project (B16026).

Conflict of interest

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

References

  • 1

    Ata-Ul-KarimS. T.CaoQ.ZhuY.TangL.RehmaniM. I. A.CaoW. (2016a). Non-destructive assessment of plant nitrogen parameters using leaf chlorophyll measurements in rice. Front. Plant Sci. 7:1829. 10.3389/fpls.2016.01829

  • 2

    Ata-Ul-KarimS. T.LiuX.LuZ.YuanZ.ZhuY.CaoW. (2016b). In-season estimation of rice grain yield using critical nitrogen dilution curve. Field Crops Res. 195, 18. 10.1016/j.fcr.2016.04.027

  • 3

    Ata-Ul-KarimS. T.YaoX.LiuX.CaoW.ZhuY. (2013). Development of critical nitrogen dilution curve of Japonica rice in Yangtze river reaches. Field Crops Res. 149, 149158. 10.1016/j.fcr.2013.03.012

  • 4

    Ata-Ul-KarimS. T.ZhuY.CaoQ.RehmaniM. I. A.CaoW.TangL. (2017a). In-season assessment of grain protein and amylose content in rice using critical nitrogen dilution curve. Eur. J. Agron.90, 139151. 10.1016/j.eja.2017.08.001

  • 5

    Ata-Ul-KarimS. T.ZhuY.LiuX.CaoQ.TianY.CaoW. (2017b). Comparison of different critical nitrogen dilution curves for nitrogen assessment in rice. Sci. Rep.7:42679. 10.1038/srep42679

  • 6

    Ata-Ul-KarimS. T.ZhuY.LuX. J.CaoQ.TianY. C.CaoW. (2017c). Estimation of nitrogen fertilizer requirement for rice crop using critical nitrogen dilution curve. Field Crops Res. 2017, 3240. 10.1016/j.fcr.2016.10.009

  • 7

    Ata-Ul-KarimS. T.YaoX.LiuX.CaoW.ZhuY. (2014a). Determination of critical nitrogen dilution curve based on stem dry matter in rice. PLoS ONE9:e104540. 10.1371/journal.pone.0104540

  • 8

    Ata-Ul-KarimS. T.ZhuY.YaoX.CaoW. (2014b). Determination of critical nitrogen dilution curve based on leaf area index in rice. Field Crops Res.167, 7685. 10.1016/j.fcr.2014.07.010

  • 9

    BannayanM.HoogenboomG. (2008). Weather Analogue: a tool for lead time prediction of daily weather data realizations based on a modified k-Nearest Neighbor approach. Environ. Model.23, 703713. 10.1016/j.envsoft.2007.09.011

  • 10

    BremnerJ. M.MulvancyC. S. (1982). Nitrogen-total, in Methods of Soil Analysis, Part II, ed PageA. L. (Madison, WI: American Society of Agronomy), 595624.

  • 11

    ColnenneC.MeynardJ. M.RocheR.ReauR. (2002). Effects of nitrogen deficiencies on autumnal growth of oilseed rape, Eur. J. Agron.17, 1128. 10.1016/S1161-0301(01)00140-X

  • 12

    ConfalonieriR.DebelliniC.PirondiniM.PossentiP.BergaminiL.BarlassinaG.et al. (2011). A new approach for determining rice critical nitrogen concentration. J. Agric. Sci.149, 633638. 10.1017/S0021859611000177

  • 13

    EickhoutB.BouwmanA. F.ZeijtsH. V. (2006). The role of nitrogen in world food production and environmental sustainability. Agric. Ecosys. Environ.116, 414. 10.1016/j.agee.2006.03.009

  • 14

    GastalF.BélangerG. (1993). The effects of nitrogen fertilization and the growing season on photosynthesis of field-grown tall fescue (Festuca arundinacea schreb.) canopies. Ann. Bot.72, 401408. 10.1006/anbo.1993.1125

  • 15

    HammerG. L.WrightG. C. (1994). A theoretical analysis of nitrogen and radiation effects on radiation use efficiency in peanut. Aust. J. Agric. Res. 45, 575589. 10.1071/AR9940575

  • 16

    HeZ.QiuX.Ata-Ul-KarimS. T.LiY.LiuX.CaoQ.et al. (2017). Development of a critical nitrogen dilution curve of double cropping rice in south China. Front. Plant Sci.8:638. 10.3389/fpls.2017.00638

  • 17

    HuangS.MiaoY.ZhaoG.YuanF.MaX.TanC.et al. (2015). Satellite remote sensing-based in-season diagnosis of rice nitrogen status in northeast China. Remote Sens. 7, 1064610667. 10.3390/rs70810646

  • 18

    JeuffroyM. H.MeynardJ. M. (1997). Azote: production agricole et environement, in Assimilation de l'azote chez les Plantes, ed Morot-GaudryJ. F. (Paris: INRA), 369380.

  • 19

    JeuffroyM. H.NeyB.OuryA. (2002). Integrated physiological and agronomic modeling of N capture and use within the plant. J. Exp. Bot.53, 809823. 10.1093/jexbot/53.370.809

  • 20

    JustesE.MaryB.MeynardJ. M.MachetJ. M.Thelier-HucheL. (1994). Determination of a critical nitrogen dilution curve for winter wheat crops. Ann. Bot.74, 397407. 10.1006/anbo.1994.1133

  • 21

    LeeK. J.LeeB. W. (2013). Estimation of rice growth and nitrogen nutrition status using color digital camera image analysis. Eur. J. Agron. 48, 5765. 10.1016/j.eja.2013.02.011

  • 22

    LemaireG.GastalF. (1997). N uptake and distribution in plant canopies, in Diagnosis of the Nitrogen Status in Crops, ed LemaireG. (Heidelberg: Springer-Verlag), 343.

  • 23

    LemaireG.JeuffroyM. H.GastalF. (2008). Diagnosis tool for plant and crop N status in vegetative stage: theory and practices for crop N management. Eur. J. Agron. 28, 614624. 10.1016/j.eja.2008.01.005

  • 24

    MillsA.MootD. J.JamiesonP. D. (2009). Quantifying the effect of nitrogen on productivity of cocksfoot (Dactylis glomerata L.) pastures. Eur. J. Agron.30, 6369. 10.1016/j.eja.2008.07.008

  • 25

    MisteleB.SchmidhalterU. (2008). Estimating the nitrogen nutrition index using spectral canopy reflectance measurements. Eur. J. Agron. 29, 184190. 10.1016/j.eja.2008.05.007

  • 26

    Peltonen-SainioP.ForsmanK.PoutalaT. (2008). Crop management effects on pre-and post-anthesis changes in leaf area index and leaf area duration and their contribution to grain yield and yield components in spring cereals. J. Agron. Crop Science. 179, 4761. 10.1111/j.1439-037X.1997.tb01146.x

  • 27

    PlénetD.LemaireG. (1999). Relationships between dynamics of nitrogen uptake and dry matter accumulation in maize crops. determination critical N concentration. Plant Soil216, 6582. 10.1023/A:1004783431055

  • 28

    RitchieJ.OtterS. (1985). Description and performance of CERES-Wheat: a user-oriented wheat yield model. ARS U.S. Dep. Agric.38, 159175.

  • 29

    SheehyJ. E.DionoraM. J. A.MitchellP. L.PengS.CassmanK. G.LemaireG.et al. (1998). Critical nitrogen concentrations: implications for high-yielding rice (Oryza Sativa, L.) cultivars in the tropics. Field Crops Res.59, 3141. 10.1016/S0378-4290(98)00105-1

  • 30

    WangY. F.YuZ. W.LiS. X.YuS. L. (2002). Effects of nitrogen application amount on content of protein components and processing quality of wheat grain. Sci. Agric. Sin. 35, 10711078.

  • 31

    WangY.ShiP. H.ZhangG.RanJ.ShiW. M.WangD. J.et al. (2016). A critical nitrogen dilution curve for japonica rice based on canopy images, Field CropsResearch, 198, 93100. 10.1016/j.fcr.2016.08.032

  • 32

    YaoX.Ata-Ul-KarimS. T.ZhuY.TianY.LiuX.CaoW. (2014b). Development of critical nitrogen dilution curve in rice based on leaf dry matter. Eur. J. Agron. 55, 2028. 10.1016/j.eja.2013.12.004

  • 33

    YaoX.ZhaoB.TianY. C.LiuX. J.NiJ.CaoW. X.et al. (2014a). Using leaf dry matter to quantify the critical nitrogen dilution curve for winter wheat cultivated in eastern China. Field Crops Res. 159, 3342. 10.1016/j.fcr.2013.12.007

  • 34

    YuanZ.Ata-Ul-KarimS. T.CaoQ.LuZ.CaoW.ZhuY.et al. (2016). Indicators for diagnosing nitrogen status of rice based on chlorophyll meter readings. Field Crops Res. 185, 1220. 10.1016/j.fcr.2015.10.003

  • 35

    YueS.MengQ.ZhaoR.LiF.ChenX.ZhangF.et al. (2012). Critical nitrogen dilution curve for optimizing nitrogen management of winter wheat production in the North China plain. Agron. J. 104, 523529. 10.2134/agronj2011.0258

  • 36

    YueS.SunF.MengQ.ZhaoR.FeiL.ChenX.et al. (2014). Validation of a critical nitrogen curve for summer maize in the north China plain. Pedosphere24, 7683. 10.1016/S1002-0160(13)60082-X

  • 37

    ZhaoB.Ata-Ul-KarimS. T.LiuZ.XiaoJ.LiuZ.QinA.et al. (2017). Development of a critical nitrogen dilution curve based on leaf dry matter for summer maize. Field Crops Res.208, 6068. 10.1016/j.fcr.2017.03.010

  • 38

    ZhaoB.Ata-Ul-KarimS. T.YaoX.TianY.CaoW.ZhuY.et al. (2016a). A new curve of critical nitrogen concentration based on spike dry matter for winter wheat in eastern China. PLoS ONE11:e0164545. 10.1371/journal.pone.0164545

  • 39

    ZhaoB.LiuZ.Ata-Ul-KarimS. T.XiaoJ.LiuZ.QiA.et al. (2016b). Rapid and nondestructive estimation of the nitrogen nutrition index in winter barley using chlorophyll measurements. Field Crops Res.185, 5968. 10.1016/j.fcr.2015.10.021

  • 40

    ZhaoB.YaoX.TianY.LiuX.Ata-Ul-KarimS. T.NiJ.et al. (2014). New critical nitrogen curve based on leaf area index for winter wheat. Agron. J. 106, 379389. 10.2134/agronj2013.0213

  • 41

    ZhaoB.YaoX.TianY.LiuX.CaoW.ZhuY. (2012). Accumulative nitrogen deficit models of wheat aboveground part based on critical nitrogen concentration. Chin. J. Appl. Ecol. 23, 31413148.

  • 42

    ZhaoD.ReddyK. R.KakaniV. G.ReddyV. R. (2005). Nitrogen deficiency effects on plant growth, leaf photosynthesis, and hyperspectral reflectance properties of sorghum. Eur. J. Agron. 22, 391403. 10.1016/j.eja.2004.06.005

  • 43

    ZhaoZ.WangE.WangZ.ZangH.LiuY.AngusJ. F. (2014). A reappraisal of the critical nitrogen concentration of wheat and its implications on crop modeling. Field Crop. Res. 164, 6573. 10.1016/j.fcr.2014.05.004

  • 44

    ZhuX. K.GuoW. S.ZhouJ. L.HongH. U.YingZ.LiC. Y.et al. (2003). Effects of nitrogen on grain yield, nutritional and processing quality of wheat for different end uses. Agric. Sci. China2, 609616.

  • 45

    ZiadiN.BélangerG.ClaessensA.LefebvreL.CambourisA. N.TremblayN.et al. (2010). Determination of a critical nitrogen dilution curve for spring wheat. Agron. J. 102, 241250. 10.2134/agronj2009.0266

  • 46

    ZiadiN.BrassardM.BélangerG.CambourisA. N.TremblayN.NolinM. C.et al. (2008). Critical nitrogen curve and nitrogen nutrition index for corn in eastern Canada. Agron. J. 100, 271276. 10.2134/agrojnl2007.0059

Summary

Keywords

wheat, leaf area duration, critical nitrogen dilution curve, nitrogen nutrition index, yield

Citation

Wang X, Ye T, Ata-Ul-Karim ST, Zhu Y, Liu L, Cao W and Tang L (2017) Development of a Critical Nitrogen Dilution Curve Based on Leaf Area Duration in Wheat. Front. Plant Sci. 8:1517. doi: 10.3389/fpls.2017.01517

Received

12 May 2017

Accepted

17 August 2017

Published

04 September 2017

Volume

8 - 2017

Edited by

Jianjun Chen, University of Florida, United States

Reviewed by

Khawar Jabran, Düzce University, Turkey; Luis Morales-Salinas, University of Chile, Chile

Updates

Copyright

*Correspondence: Liang Tang

This article was submitted to Crop Science and Horticulture, a section of the journal Frontiers in Plant Science

Disclaimer

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

Outline

Figures

Cite article

Copy to clipboard


Export citation file


Share article

Article metrics