Estimating Plant Nitrogen Concentration of Maize Using a Leaf Fluorescence Sensor across Growth Stages

: Nitrogen (N) is one of the most essential nutrients that can signiﬁcantly a ﬀ ect crop grain yield and quality. The implementation of proximal and remote sensing technologies in precision agriculture has provided new opportunities for non-destructive and real-time diagnosis of crop N status and precision N management. Notably, leaf ﬂuorescence sensors have shown high potential in the accurate estimation of plant N status. However, most studies using leaf ﬂuorescence sensors have mainly focused on the estimation of leaf N concentration (LNC) rather than plant N concentration (PNC). The objectives of this study were to (1) determine the relationship of maize (Zea mays L.) LNC and PNC, (2) evaluate the main factors inﬂuencing the variations of leaf ﬂuorescence sensor parameters, and (3) establish a general model to estimate PNC directly across growth stages. A leaf ﬂuorescence sensor, Dualex 4, was used to test maize leaves with three di ﬀ erent positions across four growth stages in two ﬁelds with di ﬀ erent soil types, planting densities, and N application rates in Northeast China in 2016 and 2017. The results indicated that the total leaf N concentration (TLNC) and PNC had a strong correlation (R 2 = 0.91 to 0.98) with the single leaf N concentration (SLNC). The TLNC and PNC were a ﬀ ected by maize growth stage and N application rate but not the soil type. When used in combination with the days after sowing (DAS) parameter, modiﬁed Dualex 4 indices showed strong relationships with TLNC and PNC across growth stages. Both modiﬁed chlorophyll concentration (mChl) and modiﬁed N balance index (mNBI) were reliable predictors of PNC. Good results could be achieved by using information obtained only from the newly fully expanded leaves before the tasseling stage (VT) and the leaves above panicle at the VT stage to estimate PNC. It is concluded that when used together with DAS, the leaf ﬂuorescence sensor (Dualex 4) can be used to reliably estimate maize PNC across growth stages.


Introduction
Maize (Zea mays L.) is one of the three major grain crops in the world along with rice (Oryza sativa L.) and wheat (Triticum aestivum L.) [1]. Nitrogen (N) is one of the most essential nutrients that significantly affect maize yield, biomass, and grain quality. However, over-application of N is common in Chinese agriculture, resulting in many environmental problems [2,3]. Precision N management strategies aiming to apply the optimal amount of N fertilizer at the right time and place can help increase crop yield and N use efficiency, as well as reduce N surplus and environmental impacts [4][5][6].
The development and implementation of proximal and remote sensing technologies has provided new opportunities for non-destructive and real-time crop N status estimation on different scales [7][8][9]. Since leaf pigment concentrations, especially the chlorophyll concentration (Chl), can affect leaf reflectance properties and are highly correlated with N status, many optical sensors have been used to measure the canopy reflectance of a specific area or the entire crop field to estimate plant N status [7,[10][11][12]. However, reflectance sensors are usually influenced by light conditions, soil and water background, and often saturate under high biomass conditions. Proximal fluorescence sensing is insensitive to soil backgrounds, environmental light, or biomass conditions and may overcome some of the problems in N status diagnosis [7,13]. As a traditional and standard indicator of plant N status, plant N concentration (PNC) is closely related to yield performance [14][15][16]. Therefore, PNC has been widely estimated through various sensing methods, and used as a reference in different N diagnostic methods, with critical PNC values established for different crops and growth stages [7,17,18]. Most importantly, unlike reflectance indices, fluorescence signals have stronger relationships with PNC as they are mainly affected by leaf Chl concentration but not by biomass or leaf area index (LAI) [13,[19][20][21].
Dualex 4 (Force-A, Orsay, France) is a portable leaf fluorescence sensor and has been used in the past few years to monitor crop physiology and study N status diagnosis [22]. Apart from measuring Chl concentration through leaf transmittance, Dualex 4 can also measure leaf epidermal flavonoids (Flav) by comparing the Chl fluorescence induced by ultra-violet (UV) excitation at 375 nm to that induced by red light at 650 nm wavelength [23][24][25], and provides a new Chl/Flav ratio called N balance index (NBI). Numerous recent studies have focused on utilizing Dualex 4 to estimate N status in a variety of crops, and have revealed a significant relationship between Dualex 4 readings and N indicators. For example, Dualex 4-based Chl readings were found to be highly related to leaf Chl concentrations in four crops including corn, soybean (Glycine max L. Merr.), spring wheat (Triticum aestivum L.), and canola (Brassica napus L.) (R 2 = 0.69-0.90) [26]. Cartelat et al. [27] showed a strong linear relationship between phenolics (Phen) measured by Dualex sensor and leaf N concentration (LNC) (R 2 = 0.76) and further displayed the correlation between Phen and N nutrition index (NNI) (r = −0.60) for wheat. It was found that NBI calculated using a chlorophyll metercombined with an older version Dualex sensor could predict PNC and NNI most accurately for muskmelon (Cucumis melo L. cv. Tezac) (R 2 = 0.79-0.93 and 0.80-0.95) [28].
However, it is sometimes difficult to successfully use various proximal and remote sensors, including fluorescence sensors, because spectral data may vary due to different factors, such as soil conditions, crop growth stages, and leaf positions [7,9,21,[29][30][31][32]. The variations of different soils in terms of quality, water, nutrition, and temperature often affect crop growth and lead to changes in crop properties [33][34][35][36]. Stress events resulting from environmental situations may eventually change the content of some compounds like leaf pigments [37][38][39]. This could modify leaf optical and fluorescence properties and be monitored by proximal or remote sensing technologies.
It has been a great challenge to use proximal and remote sensing technologies to reliably estimate PNC across growth stages [40]. The active canopy sensor GreenSeeker-based vegetation indices (VIs) could be used to predict LAI and aboveground biomass well (R 2 = 0.83-0.89), but had a poor performance for PNC estimation (R 2 = 0.47) across V5-V10 growth stages for spring maize in Northeast China [41]. Different VIs and prediction models will be needed to estimate PNC at different growth stages [40,42]. PNC can be more reliably estimated at later growth stages, but at early growth stages before crop canopy closure, the performance of most prediction models has been quite poor, even with hyperspectral remote sensing data [40,42]. Gabriel et al. [43] compared two leaf chlorophyll sensors (SPAD-502 and Dualex) for estimating maize LNC on two different dates, and demonstrated similar performance of the two sensors, with R 2 = 0.43 and 0.62 for SPAD and R 2 = 0.42 and 0.68 for Dualex. Tremblay et al. [29] found that Dualex could be used to diagnose maize N status successfully within 21 days after topdressing but lost sensitivity at later stages. A hand-held canopy fluorescence sensor Multiplex has been used to detect N status in early growth stages of maize and it was proven that the Multiplex parameters were strongly influenced by N dose [20,44]. Moreover, strong relationships between fluorescence indices and N indicators of rice at different growth stages were revealed by Huang et al. [45]. To overcome the influence of growth stage or other factors, N sufficiency index (NSI) or response index (RI) is generally calculated by using a well-fertilized area as the reference to diagnose crop N status [41,45,46], but such approach will require a well-fertilized area or N rich plot, and generally cannot improve the prediction of PNC across growth stages [45]. Therefore, more studies are still needed to develop methods for reliable prediction of PNC across growth stages.
Several studies have reported the use of fluorescence sensing technology to evaluate crop N status. For example, Yang et al. [47] showed that there were consistent positive correlations between fluorescence parameters and LNC for different rice cultivars with R 2 varying from 0.70 to 0.90. Another study demonstrated that NBI and NBI1 based on the Multiplex 2 fluorescence sensor were linearly related to LNC with a high coefficient of determination for two turfgrass cultivars (R 2 = 0.85-0.87 and R 2 = 0.75-0.78, respectively) [48]. A study using a Dualex sensor calibrated three optical indices (Chl, Flav and NBI) against LNC of grapevine (Vitis vinifera L.) and verified that NBI was optimal for estimating LNC with a root-mean-square error (RMSE) smaller than 2 mg of N g −1 dry weight. The threshold values of NBI from 11 to 18 at flowering and 8 to 11 at bunch closure were proposed [17]. Recently, Zhang et al. [49] found Dualex 4 sensor parameters were significantly correlated with rice PNC at different growth stages (R 2 = 0.43-0.77) or across growth stages (R 2 = 0.52-0.69).
To date, few studies have reported on how to use leaf fluorescence sensors to accurately estimate crop PNC, especially developing strategies to effectively overcome the influence of different growth stages. Therefore, the objectives of this study were to (1) determine the relationship of maize LNC and PNC, (2) evaluate the main factors influencing the variations of leaf fluorescence sensor parameters, and (3) develop a practical strategy to reliably estimate maize PNC using leaf fluorescence sensor across different growth stages.

Study Site and Soil Description
The study was conducted in Lishu County (43 • 02 -43 • 46 N, 123 • 45 -124 • 53 E), Jilin Province in Northeast China from 2016 to 2017. Located in North Temperate Zone with four distinct seasons, this region has a semi-humid continental monsoon climate. The mean annual average temperature is 6.6 • C, and the annual average precipitation is 556 mm, about 80% of which occurs during the crop growing season from May to September. Figure 1 shows the precipitation distribution and mean temperature during the growing season in 2016-2017 in Lishu.

138
Two sites with different soil types were selected for this study. Site 1 has Aeolian sandy soil 139 (loamy sand), which is classified as typic Cryopsamments. Site 2 has Black soil (loamy clay), which is

143
All the treatments of the experiment were the same in both fields. Using a split-plot design with 144 three replications, a local maize cultivar Liangyu 66 was sown in early May each year with planting 145 densities of 55,000, 70,000 and 85,000 plant ha −1 as the main plots, and six N treatments were 146 established as the subplots, with total N doses of 0, 60, 120, 180, 240, and 300 kg ha -1 , respectively.

147
The N fertilizers were applied in two split applications: 1/3 was broadcasted and incorporated into 148 the soil with rotary tillage as basal N using ammonium sulfate before sowing, and 2/3 was banded as 149 side-dress N using urea at around V8-V9 growth stage in early July. The subplots in a wide-narrow 150 row planting pattern were 9 × 12 m with 1 m wide alley between the subplots. Sufficient phosphate 151 (90 kg P2O5 ha −1 ) and potash (90 kg K2O ha −1 ) fertilizers were applied before sowing to make sure P 152 and K nutrients were not limiting for each plot.

153
There was no irrigation at Site 2, while about 50 mm of water was irrigated into the soil around 154 mid-July at Site 1 in each year due to water stress in sandy soil. All plots were kept free of weeds,  Two sites with different soil types were selected for this study. Site 1 has Aeolian sandy soil (loamy sand), which is classified as typic Cryopsamments. Site 2 has Black soil (loamy clay), which is classified as typic Haploboroll in the United States Department of Agriculture (USDA) Soil Taxonomy.

Experimental Design
All the treatments of the experiment were the same in both fields. Using a split-plot design with three replications, a local maize cultivar Liangyu 66 was sown in early May each year with planting densities of 55,000, 70,000 and 85,000 plant ha −1 as the main plots, and six N treatments were established as the subplots, with total N doses of 0, 60, 120, 180, 240, and 300 kg ha −1 , respectively. The N fertilizers were applied in two split applications: 1/3 was broadcasted and incorporated into the soil with rotary tillage as basal N using ammonium sulfate before sowing, and 2/3 was banded as side-dress N using urea at around V8-V9 growth stage in early July. The subplots in a wide-narrow row planting pattern were 9 × 12 m with 1 m wide alley between the subplots. Sufficient phosphate (90 kg P 2 O 5 ha −1 ) and potash (90 kg K 2 O ha −1 ) fertilizers were applied before sowing to make sure P and K nutrients were not limiting for each plot.
There was no irrigation at Site 2, while about 50 mm of water was irrigated into the soil around mid-July at Site 1 in each year due to water stress in sandy soil. All plots were kept free of weeds, insects, and diseases with pesticides based on local standard practices. Detailed information about the experiments conducted in this study in 2016 and 2017 is listed in Table 1.

Dualex 4 Sensor Data Collection, Plant Sampling, and Measurements
Three representative plants located in the center rows of each plot of six N treatments at the 70,000 plants ha −1 density in the two fields were selected to be cut at the ground level at V8, V12, and VT in 2016, and V4, V11, and VT in 2017. Each single leaf was separated from the whole plant, and leaves in the same position of the three plants in each plot were mixed together. The stems were also mixed together.
The Dualex 4 sensor (Force-A, Orsay, France) was used in this study for proximal sensing. Three representative plants located in the inner rows of each plot of six N treatments at three densities were selected to be sampled and measured by the sensor. In particular, this sensor measures a leaf surface area of 20 mm 2 . The plant samples and sensor readings were obtained at V6, V8, V12, and VT growth stages in 2016 and at V6, V8, and VT growth stages in 2017 in each field. Dualex 4 values were measured at around the leaf blade midpoint to avoid midribs or physical damage on the adaxial side (upper side) of the uppermost, second, and third fully expanded leaves before VT stage, as well as the leaf above the panicle, panicle leaf, and the leaf below the panicle at the VT growth stage for each plant in each plot. After sensing, all plant samples in each plot were separated into leaves and stems except for the V6 growth stage.
All samples were oven-dried at 105°C for 30 min, then dried at 70°C to a constant weight, and ground into fine powders to determine N concentration using a modified Kjeldahl digestion method [50]. The total NLC (TLNC) is the sum of the product of the N concentration of each leaf and its proportional weight, while PNC is calculated by adding the product of the N concentration of each organ and its proportional weight. As the leaves were the main component of the plant at the V6 growth stage, the TLNC at this stage was also used as PNC in this study.
The data acquisition dates and days after sowing (DAS) are shown in Table 1. For the convenience of discussion, the first, second, and third leaves counted from the top of the maize plant selected for sensor measurement in this research were abbreviated as Leaf 1, Leaf 2, and Leaf 3, respectively. The three Dualex 4 parameters (Chl, Flav, and NBI) at the same leaf position of three plants in each plot were averaged and used as the mean reading of that leaf position for the plot. The details of the Dualex 4 parameters are shown in Table 2.

Statistical Analysis
The data collected were pooled together to establish the relationships between single LNC (SLNC) and TLNC and PNC. In addition, the data collected by fluorescence sensing for each site, year, growth stage, and N rate were divided into a calibration dataset (two replications of the observations) and a validation dataset (one replication of the observations). The mean, standard deviation (SD), coefficient of variation (CV, %) of SLNC, TLNC, and PNC were calculated using Microsoft Excel (Microsoft Corporation, Redmond, WA, USA).
Besides three raw Dualex 4 parameters, three modified parameters derived from them based on DAS were used in this study: modified Chl (mChl), modified Flav (mFlav), and modified NBI (mNBI) were computed as the ratio of Chl to DAS, the product of Flav and DAS, and the ratio of NBI to the square of DAS ( Table 2). The coefficient of determination (R 2 ) values for all relationships were calculated using SPSS 21.0 (SPSS Inc., Chicago, IL, USA), and the model with the highest R 2 between the Dualex parameters and the two N indices (TLNC and PNC) was selected and listed in this paper, and was used for further estimation of PNC. The performance of the established relationship models for predicting PNC was evaluated by comparing R 2 , RMSE, and relative error (RE). The higher the R 2 and the lower the RMSE and RE, the higher the accuracy of the prediction models. Agronomic and proximal sensing data were subjected to the least significant difference (LSD) test at a 5% significance level to assess differences between the means of treatments using SPSS 21.0.

Interrelationships of SLNC, TLNC, and PNC
We plotted SLNC against TLNC and PNC in Figure 2. The SLNC was highly related to both TLNC and PNC with R 2 varying from 0.91 to 0.98. The best relationships appeared in the results for Leaf 1 (R 2 = 0.98 and 0.96), followed by Leaf 2 and finally Leaf 3. As the sum of each single leaf, the TLNC was more related to SLNC than that of the whole plant, because the PNC was also affected by the stem N concentration. However, the strong correlations of TLNC and PNC to SLNC indicated the high probability to use information acquired from a single leaf, especially Leaf 1, to understand the N status of maize through the prediction of TLNC or PNC.   Figure 3 indicated that PNC was slightly lower than TLNC, and both were significantly affected by the growth stage and N rate. Different soil types only resulted in different TLNC values (mean = 22.91 and 23.93 g kg −1 for Aeolian sandy soil field (Site 1) and Black soil field (Site 2), respectively). For PNC, there was no significant difference between the two fields, with 20.07 and 20.69 g kg −1 for Aeolian sandy soil field (Site 1) and Black soil field (Site 2), respectively (Figure 3a). In contrast, TLNC of four growth stages differed significantly, from 17.06 g kg −1 at the V12 growth stage to 31.03 g kg −1 at the V6 growth stage, while PNC ranged from 12.10 g kg −1 at the VT growth stage to 31.28 g kg −1 at the V6 growth stage, without a significant difference between the V12 and VT growth stages (Figure 3b). The mean values of TLNC and PNC increased with N rates from the lowest under 0 kg ha −1 N rate treatment (18.10 g kg −1 for TLNC and 16.25 g kg −1 for PNC) to the highest under 300 kg ha −1 N rate treatment (27.09 g kg −1 for TLNC and 23.44 g kg −1 for PNC). The values of TLNC were significantly different under low N treatments (0 and 60 kg ha −1 ), but not under high N treatments (240 and 300 kg ha −1 ). However, the values of PNC were not significantly different under either low or high N treatments (Figure 3c).

Effects of Soil Type, Growth Stage, and N Rate on Dualex 4 Parameters
The differences in Chl, Flav, and NBI measured with a Dualex 4 sensor as influenced by soil types, growth stages, and N rates were analyzed and displayed in Figure 4. The Dualex 4 sensor-based fluorescence parameter values obtained from the three differently positioned leaves did not exhibit significant variability under the two soil types, except for Chl values of Leaf 2, with 41.70 for Aeolian sandy soil field (Site 1) and 43.14 for Black soil field (Site 2) (Figure 4a,d,g). Fluorescence parameter values changed less significantly with growth stages than TLNC and PNC (Figures 3b and 4b,e,h). For Chl, there was no consistent trend as maize grew, although there was a significant difference at late growth stages for Leaf 1 and at early stages for Leaf 2 and Leaf 3 (Figure 4b). Flav values showed a constant decreasing trend with the growth of maize, similar to the changes of TLNC and PNC, but with less variation. Comparatively, the variation trend of Leaf 1 was the same as that of PNC, while the values of the other two leaves changed significantly with growth stages (Figures 3b and 4e). NBI also showed significant differences among different growth stages, particularly for Leaf 2 and Leaf 3, but NBI and TLNC and PNC exhibited opposite trends from Flav, increasing with the increase of DAS (Figures 3b and 4h). Besides, the Chl and NBI values increased with N rate up to 240 kg ha −1 for all three differently positioned leaves (Figure 4c,i). There was a negative relationship between Flav and N rate, with higher levels of Flav found in low N treatments and lower levels of Flav found in high N treatments. However, Flav did not vary much from 180 to 300 kg ha −1 N rates (Figure 4f)

284
The results above indicated that growth stage was the main factor that influences the general  The results above indicated that growth stage was the main factor that influences the general model's estimation of maize N concentration. To overcome the problem of the high inconsistency between the Dualex 4 parameters and the TLNC and PNC across growth stages, the modified Dualex 4 parameters incorporating the information of DAS were computed and the variations of these new parameters are displayed in Figure 5. It was evident that changes in these new parameters of differently positioned leaves varied markedly among different growth stages after combining with DAS, although a reverse trend was observed in mFlav compared with the original Flav. The mNBI showed a similar trend as mChl with the growth stage.

300
All data acquired from the two study sites across growth stages were gathered together and 301 divided into calibration and validation datasets to establish a general model for PNC estimation 302 (Table 3).

303
The mean values of TLNC (23.08-23.59 g kg -1 ) were larger than those of PNC (19.91-20.61 g kg -304 1 ) in both datasets. For the calibration dataset, across all growth stages, the PNC was more variable 305 (CV = 42%) than the TLNC (CV = 31%), and the validation dataset had similar variability as the 306 calibration dataset with a CV of 44% for PNC and of 31% for TLNC.

Relationships between Dualex 4 Parameters and TLNC or PNC
All data acquired from the two study sites across growth stages were gathered together and divided into calibration and validation datasets to establish a general model for PNC estimation (Table 3). The mean values of TLNC (23.08-23.59 g kg −1 ) were larger than those of PNC (19.91-20.61 g kg −1 ) in both datasets. For the calibration dataset, across all growth stages, the PNC was more variable (CV = 42%) than the TLNC (CV = 31%), and the validation dataset had similar variability as the calibration dataset with a CV of 44% for PNC and of 31% for TLNC.
The original Dualex 4 readings, Chl, Flav, and NBI, demonstrated poor relationships with TLNC or PNC across growth stages, with R 2 ranging from 0.01 to 0.34 (Table 4). The mean values of Chl, Flav, and NBI for the three differently positioned leaves were calculated and their relationships with TLNC and PNC at four growth stages are displayed in Figure 6. The impact of growth stages on these relationships were significant.

325
The mChl was most related to TLNC and PNC with R 2 ranging from 0.75 to 0.84 for single-326 position leaves. mNBI also performed very well, with R 2 ranging from 0.74 to 0.80. Nevertheless, the 327 result for mFlav (R 2 = 0.49-0.56) was not as good as the above-mentioned two parameters, but it was 328 much better than that of the original Flav (Table 4). A power function relationship was found for the 329 models established based on the values measured on every single leaf for mChl, mFlav, and mNBI.

330
Leaf 1 showed the best potential to estimate PNC with the highest R 2 values (Table 4) Table 4). Compared with TLNC, PNC was more related to the modified 334 Dualex 4 parameters (mChl, mFlav, and mNBI). This indicated that the modified Dualex 4 parameters 335 had greater potential for direct estimation of PNC accurately (Figure 7 and Table 4). However, after combining with DAS, which was directly connected to the growth stage, the modified Dualex 4 parameters (mChl, mFlav, and mNBI) were highly related to TLNC and PNC (Table 4 and Figure 7). As maize grew, the values of mChl and mNBI decreased with the increase of TLNC and PNC (Figure 7a,c,d,f), while the values of mFlav increased with the growth stages and exhibited negative relationships with both TLNC and PNC accordingly (Figure 7b,e).

337 338
The mChl was most related to TLNC and PNC with R 2 ranging from 0.75 to 0.84 for single-position leaves. mNBI also performed very well, with R 2 ranging from 0.74 to 0.80. Nevertheless, the result for mFlav (R 2 = 0.49-0.56) was not as good as the above-mentioned two parameters, but it was much better than that of the original Flav (Table 4). A power function relationship was found for the models established based on the values measured on every single leaf for mChl, mFlav, and mNBI. Leaf 1 showed the best potential to estimate PNC with the highest R 2 values (Table 4). Moreover, by calculating the mean values of each parameter of the three leaves, the correlations between the modified Dualex 4 parameters and TLNC or PNC were improved to a certain extent, especially for mFlav (Figure 7b,e and Table 4). Compared with TLNC, PNC was more related to the modified Dualex 4 parameters (mChl, mFlav, and mNBI). This indicated that the modified Dualex 4 parameters had greater potential for direct estimation of PNC accurately (Figure 7 and Table 4). Remote Sens. 2020, 1, x FOR PEER REVIEW 12 of 22

355
showed the highest accuracy among the three positioned leaves (Table 5)

The Estimation of PNC Using Different Modified Dualex 4 Parameters
In order to establish good relationships between modified Dualex 4 parameters and PNC as well as TLNC, the above models describing their relationships were further evaluated with the validation dataset. The validation results of mChl for estimating TLNC and PNC were the best with the highest R 2 (0.71-0.79 and 0.73-0.84) and the lowest RMSE (3.27-0.81 and 3.46-4.48) and RE (14%-16% and 17-23%) ( Table 5). Although the mFlav showed a better relationship with TLNC and PNC than the original Flav, it still did not perform as well as the other two modified Dualex parameters in predicting TLNC and PNC, with the lowest R 2 (0.42-0.53 and 0.40-0.53) and the highest RMSE (4.83-5.38 and 5.96-6.73) and RE (21%-23% and 30%-34%) ( Table 5). For single leaf, mChl obtained from Leaf 1 predicted TLNC and PNC most accurately, while mFlav and mNBI obtained from Leaf 3 showed the highest accuracy among the three positioned leaves (Table 5). Besides, the results were further improved by averaging the mFlav and mNBI readings measured from three differently positioned leaves (Table 5 and Figure 8b,c,e,f). The data distribution for the V6 growth stage deviated from the 1:1 line more significantly than the other growth stages, which would influence the estimation results to some extent (Figure 8). Above all, the modified Dualex 4 parameters showed a better performance for direct PNC estimation across growth stages, but not for TLNC estimation, which implied that it was needless to estimate maize PNC by estimating TLNC using sensor readings acquired primarily from single leaves. acquired primarily from single leaves.

Feasibility of Estimating Maize N Status Using Single Leaf-based Dualex 4 Parameters
The SLNC of the fully expanded three leaves selected in our study was linearly related to both TLNC and PNC, with high R 2 values. Considering the good relationships between SLNC and TLNC and PNC (Figure 2), and the significant change trend of Dualex 4 parameters with the increase of N rate (Figures 3c and 4c,f,i), this study verified that it would be possible to predict maize N status using the information from a single leaf of maize and be possible to predict PNC directly or indirectly by predicting TLNC first.
The Chl values elevated with an increasing N rate (Figure 4c), which conforms to the finding that leaf chlorophyll is strongly affected by leaf N [51][52][53]. Leaf N contributed a large proportion to plant photosynthesis [54]. However, as an indicator of epidermal flavonoids, a carbon-based secondary metabolite, Flav exhibited a reverse trend with the changes of N rate (Figure 4f) because the carbon-based metabolite would be produced under low N level and would consequently cause the enhanced synthesis of flavonoids [55,56]. As the ratio of Chl to Flav, NBI increased with N rate ( Figure 4i) and has been considered to be a sensitive indicator of crop N status [13,28,44,57]. The results of this study indeed showed that the range of NBI variation was greater than that of Chl and Flav under different N rates (Figure 4c,f,i).
The Dualex 4 parameters' good discrimination ability of different N rates and the saturation at high N levels (240 and 300 kg ha −1 ) (Figure 4) was in agreement with other studies using the fluorescence-based method [58]. The suggested optimal N rate for the experimental sites in our study was approximately 180-200 kg ha −1 [59], so the N rates of 240 and 300 kg ha −1 exceeded the optimal level. Some previous studies have concluded that fluorescence indices were insensitive to high crop N content in the range of optimal to excessive N levels [58,60]. In this study, as shown in Figure 3c, both TLNC and PNC varied clearly with the changes of N rate from 0 through 240 kg ha −1 , but there was no significant difference between 240 and 300 kg ha −1 N rates. Thus, in this study, the saturation of Dualex 4 parameters on high N rates was closely related to the N status of maize, and we could not draw a similar conclusion that the fluorescence indices were insensitive to high N content here.

Main Factor(s) Affecting the Establishment of the General Model and the Best Parameter(s) for PNC Estimation
Compared with stems, leaves are more important for photosynthesis. However, maize plants grown in the Aeolian sandy soil field (Site 1) were much more vulnerable to water deficiency than maize plants grown in the Black soil field (Site 2). Due to the N accumulation in stems, the PNC obtained in the Aeolian sandy soil field increased when the N concentration of stems was taken into account. Therefore, the difference in PNC between the two fields vanished (Figure 3a). The significant variation observed in TLNC values but the similarity observed in PNC values in the two study sites (Figure 3a) implies that it may be more difficult to evaluate TLNC from different soil fields, and it would be more feasible to establish a PNC predicting model neglecting the effect of soil types. Nevertheless, it is necessary to distinguish soil types when estimating PNC indirectly by estimating TLNC first (Figure 4a,d,g).
The changes of TLNC or PNC as affected by growth stages (Figure 3b) have been verified in diverse crops in other studies [28,61,62]. This is why growth stages must be strictly defined and distinguished for N status evaluation in most studies. Thus, the growth stage plays an important role in establishing general models for PNC estimation according to the findings of this study.
Chl measured by Dualex 4 has been calibrated in µg cm −2 units [22], which means Chl is a proxy of surface-based N. However, the TLNC and PNC mentioned in this study were both calculated based on mass. The area-based LNC changes very little during the growing season due to the expression of area-based LNC forces changes in specific leaf areas [63]. The time we selected to conduct the sensing was before the reproductive phase, and there was no obvious redistribution of N from the leaves to the ears that would result in the significant reduction of LNC. Furthermore, the Dualex 4 only measures a green leaf surface area of 20 mm 2 that would not exhibit too much variation of N concentration. Thus, the surface-based parameter Chl is relatively more stable than the mass-based TLNC and PNC (Figures 3b and 4b).
Different from Chl, the epidermal Flav estimated by Dualex 4, which is the representative of total leaf Flav [27,64], is considered as a surrogate of dry leaf mass per area with no units [22]. Meyer et al. [23] reported the positive correlations between mass-based polyphenols (mainly Flav) and leaf mass per area for woody plants, which revealed the regularity of the decrease in Flav was caused by the increase in leaf dry mass. Thicker leaves had a larger leaf mass per area, and the accumulation of organic matter in leaves played the role of Flav dilution. In addition, a previous study has reported that the accumulation of Flav was highly sensitive to light intensity [65]. As the maize grows, the leaf area and canopy closure of the uppermost leaves gradually increase, resulting in increased degree of shading. This may limit the expression of Flav synthetic genes and lead to a reduction in flavonoid content in these leaves during later growth stages. More studies are still needed on the periodic changes of Flav.
As the ratio of Chl to Flav, NBI is equivalent to mass-based Chl (Chl g −1 dry mass). Since Chl is relatively stable and Flav showed a downward trend as the maize grew, it is easy to understand the reverse trend of NBI to Flav (Figure 4h). The inconsistent relationships between the original Dualex 4 parameters (Chl, Flav, and NBI) and TLNC or PNC at different growth stages led to poor model performance across the growth stages (Table 4 and Figure 6). This makes the original Dualex 4 parameters unfeasible for estimating TLNC and PNC across growth stages.
Varvel et al. [46] pointed out that it is effective to evaluate crop N status by calculating growing degree days (GDD) when a specific growth stage cannot be determined. Following the idea of Teal et al. [66], who predicted maize yield using normalized difference vegetation index (NDVI) and GDD, DAS was used together with Dualex 4 parameters in this study to explore the potential of real-time PNC estimation across growth stages. The results of this study indicated the relationships between the DAS-based modified Dualex 4 parameters and TLNC and PNC across growth stages were significantly improved and were simulated by power function models. Furthermore, this study revealed a stronger relationship between the modified parameters and PNC than TLNC. In addition, the good relations between SLNC and PNC implied that PNC could be directly estimated well using a leaf sensor instead of indirect estimation through TLNC (Table 5 and Figure 8). This indicates the great potential of leaf-based Dualex 4 for assessing plant N status, not just maize leaf N status.
It has been reported that polyphenols (Phen, including Flav) may be a more specific indicator of crop N status than Chl because of the sensitivity of Chl to sulfur stress [67]. Padilla et al. [58] demonstrated the strong linear relationships between Flav, measured by a Multiplex sensor, and cucumber (Cucumis sativus 'Strategos') N concentration (R 2 = 0.87-0.95). Besides, some studies pointed out the superiority of NBI (or Chl/Phen) in N status estimation in a variety of crops, including maize [28,29,31,44,48,57,68]. However, mChl (Chl) performed best while mFlav (Flav) performed worst in this study.
The relatively poor performance of mFlav indicated the difficulty of estimating PNC accurately using Flav, while the difficulty for estimating Phen in leaves has been attributed to different solubility, distribution locations, and species varieties [23]. Flav has previously been shown to perform well in N status estimation for other crops such as vegetable, rice, and sweet cherry (Prunus avium L.) [49,52,68,69]. However, the relatively thick maize leaves could affect the transmittance of the light emitted by the sensor and then lead to the change of the parameter [22]. Measurements taken from both the adaxial (top) and abaxial (bottom) sides of the leaves may help improve maize PNC estimation by Flav, although the flavonoid contents in the adaxial and abaxial sides of the leaves are highly correlated in some crops [27,57]. Moreover, the change of Flav content is influenced by various factors in addition to N such as leaf thickness, light condition, water stress, pathogen attack, low temperature, and the onset of senescence [17,30,65,70]. Although affected by the poor performance of mFlav, the estimation results of mNBI-based models showed similar R 2 values with mChl, indicating the good performance of NBI for N indices estimation. These results conform to the findings of some previous studies [27][28][29]31,44,49,57].

The Most Suitable Leaf Position for Sensing Measurements and PNC Estimation
We tested three differently positioned leaves at each growth stage in this study in order to understand their effects on evaluating maize N status. The changes of the Dualex 4-based parameters with leaf positions were relatively smaller than that with the growth stage and N rate. Nevertheless, there were slight differences among the three leaves, which consequently displayed different performance in PNC prediction.
It was noted that leaf Chl and Flav contents depended on leaf age (leaf position) induced by light conditions [71,72]. The photosynthetic photon flux density would decrease with increasing canopy depth during crop growth and canopy development [73]. The priority supplement of N in younger leaves (upper leaves) re-translocated from older leaves (lower leaves) leads to the preferential distribution of N to the upper leaves [74] and consequently leads to higher rates of canopy photosynthesis, which needs a higher content of Chl. The finding of higher Chl values in Leaf 1 (the uppermost leaf) than in Leaf 3 (the lowest leaf) was similar to the study results for rice [74,75]. Higher Flav contents in older leaves have been reported [71,76,77], similar to the results of this study. Besides, the decrease of leaf mass per area along with leaf thickness as leaf position declines, the limited sunlight in the lower leaves mentioned above caused the change of Flav with the leaf position. However, there was a slightly increasing tendency of Chl measurement in Leaf 2, which may be caused by the close locations of the three leaves. As the uppermost leaf (Leaf 1) was newly expanded soon after Leaf 2, the N concentration of these two leaves would not differ too much. In addition, the average values used in this analysis may help reduce some variations. The Flav is an N-free compound, so its continuous change with leaf positions was more distinct.
To diagnose crop N status quickly and effectively, the test frequency of a leaf sensor must be reduced. Thus, selecting the most appropriate leaf rather than several leaves of a plant for sensing measurement with high accuracy and reliability is important. The measurements taken from Leaf 1 for all three modified Dualex parameters provided the highest model R 2 in TLNC and PNC predictions. The best performance of mFlav appeared in Leaf 3 in PNC estimation, and the mean values of three leaves improved the results significantly, but the results were not as good as the other two parameters. For mNBI, although the best result appeared in Leaf 3 too, it did not show a significant superiority over the results of the other two leaves. Thus, in conclusion, the measurement of Leaf 1 may be sufficient to estimate maize N status using Dualex 4-based mChl and mNBI whereas Leaf 3 is preferred to accurately estimate PNC using mFlav.

Implications for Practical Application and Future Research Needs
This study found that combining with DAS, the Dualex 4 leaf fluorescence sensor could be used to reliably estimate maize PNC across growth stages, even if only one leaf (uppermost) was sampled. The estimated PNC can then be compared with threshold or critical PNC values to diagnose maize N status by calculating N nutrition index (NNI) [41]. For this purpose, maize aboveground biomass will be needed and can be reliably estimated using active canopy sensors, like GreenSeeker or Crop Circle sensors [41,78]. Research is needed to develop methods to use Dualex 4 sensor to estimate NNI directly, as demonstrated for SPAD chlorophyll meter [79]. Based on the relationship between SPAD chlorophyll meter readings and PNC and the critical N concentration curve for maize, Yang et al. [80] established the critical or optimal chlorophyll meter (SPAD) reading curve. Measured chlorophyll meter readings can then be compared with critical chlorophyll meter readings at specific biomass to calculate NNI, without the need to estimate PNC. This is an innovative idea and should be tested with the Dualex sensor.
The relationship between leaf and canopy fluorescence parameters could be further investigated through the use of a simple model taking into account leaf fluorescence profile inside the canopy, structural variations of the canopy, and background reflection [81]. Such information can be useful to evaluate the possibility to use the proposed method in conjunction with remote sensing fluorescence measurements obtained from an unmanned aerial vehicle (UAV) [82], aerial or satellite remote sensing [83,84].
In theory, NBI combining both Chl and Flav should be more sensitive to crop N status than using Chl or Flav alone, as found by previous research [28,29,31,44,48,57,68]. However, NBI did not perform better than Chl in this study. Those previous studies either used an earlier model of Dualex sensor (cannot measure Chl) and chlorophyll meter to calculate NBI or used canopy fluorescence sensor Multiplex or studied other crops. Therefore, more studies are needed to confirm our results with maize and further evaluate the potential of improving corn N status diagnosis using both Chl and Flav information, as well as other related soil, weather, and management data.

Conclusions
This study demonstrated the reliability of maize PNC estimation by fluorescence parameters obtained from single leaves using a leaf sensor Dualex 4. The fluorescence parameter values did not exhibit significant variability under the two different soil conditions, which indicated the needlessness of establishing a soil-specific predicting model in this study, while the variations between the PNC and Dualex 4 parameters at different growth stages make it difficult to estimate PNC across all growth stages using the original Dualex 4 parameters. Nevertheless, the modified Dualex 4 parameters using the information of DAS overcame the problem caused by growth stage changes and allowed accurate estimation of PNC using a general model across growth stages. Among the three parameters obtained by Dualex 4, mChl and mNBI were more reliable indicators for PNC estimation. It was sufficient to take the fluorescence measurement from the uppermost leaf of maize. Further research needs to focus on improving the accuracy of PNC estimation in early growth stages and practical methods for maize N status diagnosis using proximal leaf fluorescence sensors.