Next Article in Journal
On the Segmentation of the Cephalonia–Lefkada Transform Fault Zone (Greece) from an InSAR Multi-Mode Dataset of the Lefkada 2015 Sequence
Next Article in Special Issue
Estimating Nitrogen from Structural Crop Traits at Field Scale—A Novel Approach Versus Spectral Vegetation Indices
Previous Article in Journal
Diagnosis of Xinmo (China) Landslide Based on Interferometric Synthetic Aperture Radar Observation and Modeling
Previous Article in Special Issue
Modeling Mid-Season Rice Nitrogen Uptake Using Multispectral Satellite Data
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

In-Season Diagnosis of Rice Nitrogen Status Using Proximal Fluorescence Canopy Sensor at Different Growth Stages

1
International Center for Agro-Informatics and Sustainable Development (ICASD), College of Resources and Environmental Sciences, China Agricultural University, Beijing 100193, China
2
Institute of Geography, University of Cologne, 50923 Köln, Germany
3
Chinese Academy of Agricultural Engineering Planning & Design, Beijing 100125, China
4
Precision Agriculture Center, Department of Soil, Water and Climate, University of Minnesota, St. Paul, MN 55108, USA
5
Department of Geography, Minnesota State University, Mankato, MN 56001, USA
6
National Engineering and Technology Center for Information Agriculture, Nanjing Agricultural University, Nanjing 210095, China
7
Key Laboratory of Digital Earth Science, Institute of Remote Sensing and Digital Earth, Chinese Academy of Sciences, Beijing 100094, China
*
Author to whom correspondence should be addressed.
Remote Sens. 2019, 11(16), 1847; https://doi.org/10.3390/rs11161847
Submission received: 7 July 2019 / Revised: 1 August 2019 / Accepted: 6 August 2019 / Published: 8 August 2019
(This article belongs to the Special Issue Remote Sensing for Precision Nitrogen Management)

Abstract

:
Precision nitrogen (N) management requires an accurate and timely in-season assessment of crop N status. The proximal fluorescence sensor Multiplex®3 is a promising tool for monitoring crop N status. It performs a non-destructive estimation of plant chlorophyll, flavonol, and anthocyanin contents, which are related to plant N status. The objective of this study was to evaluate the potential of proximal fluorescence sensing for N status estimation at different growth stages for rice in cold regions. In 2012 and 2013, paddy rice field experiments with five N supply rates and two varieties were conducted in northeast China. Field samples and fluorescence data were collected in the leaf scale (LS), on-the-go (OG), and above the canopy (AC) modes using Multiplex®3 at the panicle initiation (PI), stem elongation (SE), and heading (HE) stages. The relationships between the Multiplex indices or normalized N sufficient indices (NSI) and five N status indicators (above-ground biomass (AGB), leaf N concentration (LNC), plant N concentration (PNC), plant N uptake (PNU), and N nutrition index (NNI)) were evaluated. Results showed that Multiplex measurements taken using the OG mode were more sensitive to rice N status than those made in the other two modes in this study. Most of the measured fluorescence indices, especially the N balance index (NBI), simple fluorescence ratios (SFR), blue–green to far-red fluorescence ratio (BRR_FRF), and flavonol (FLAV) were highly sensitive to N status. Strong relationships between these fluorescence indices and N indicators, especially the LNC, PNC, and NNI were revealed, with coefficients of determination (R2) ranging from 0.40 to 0.78. The N diagnostic results indicated that the normalized N sufficiency index based on NBI under red illumination (NBI_RNSI) and FLAV achieved the highest diagnostic accuracy rate (90%) at the SE and HE stages, respectively, while NBI_RNSI showed the highest diagnostic consistency across growth stages. The study concluded that the Multiplex sensor could be used to reliably estimate N nutritional status for rice in cold regions, especially for the estimation of LNC, PNC, and NNI. The normalized N sufficiency indices based on the Multiplex indices could further improve the accuracy of N nutrition diagnosis by reducing the influences of inter-annual variations and different varieties, as compared with the original Multiplex indices.

1. Introduction

Nitrogen (N) is an essential nutrient for plant growth and development. However, excessive N fertilizer applications have led to severe environmental impacts in China [1,2]. Therefore, there has been a growing interest in developing precision N management strategies in agricultural research for many years. This requires the development of efficient and timely crop N status diagnosis strategies and technologies [3].
Plant or leaf N concentration is an indication of N nutritional status. The traditional N testing method in the laboratory is time-consuming and impractical for characterizing spatial and temporal variability in crop N status in precision N management. Alternatively, it is known that there is a strong relationship between plant chlorophyll content and N content, although other stress factors, such as water deficit, cold, heat, excess light, leaf or root pathogens, etc., can also influence chlorophyll content directly or indirectly [4,5]. Therefore, various instruments based on measuring chlorophyll have been developed to provide indirect, nondestructive, and real-time estimations of leaf N content [3,6]. For example, the Soil Plant Analysis Development chlorophyll meter (SPAD) developed by Minolta Camera Co. (Osaka, Japan) is a widely used portable instrument for measuring chlorophyll in leaves. The SPAD meter measures the difference in absorption between the red (660 nm) and near-infrared (940 nm) wavelengths [7,8]. Leaf chlorophyll absorbs red light but not infrared, therefore, the SPAD readings indicate plant chlorophyll concentration and N content [9,10,11]. However, the reliability of SPAD results is affected by factors such as growth stage, irradiance, water status, and leaf thickness [3,5,9].
Optical non-destructive remote sensing methods based on canopy reflectance measurements have also been widely used [12,13]. The high measuring efficiency of reflectance spectroscopy sensors and the strong correlation between their measurements and crop physiological and biochemical parameters offer a high potential for N management [14,15]. Proximal active sensors, such as GreenSeeker (NTech Industries, Inc., Ukiah, CA) and Crop Circle (Holland Scientific, NE, USA), have been used to diagnose N nutritional status in real-time and to guide in-season precision management for rice N fertilization [16,17]. However, the results based on the canopy reflectance are affected by various factors, such as soil characteristics, crop growth stages, and saturation under high biomass conditions [16,18]. In addition, it is more difficult to estimate chlorophyll or N status using optical remote sensing methods as the contribution of leaf area index and biomass to canopy reflectance is much greater than that of chlorophyll or N concentration [19,20].
Unlike reflectance indices, the fluorescence spectra are less affected by biomass or leaf area index [15,21]. At different N nutritional levels, the fluorescence intensities of leaves are significantly different near the 440 nm (Blue, B), 525 nm (Green, G), 685–690 nm (Red, R), and 735–740 nm (Near-infrared, NIR) wavelengths [22,23]. Studies have shown that the fluorescence ratio of NIR and R bands is highly correlated with chlorophyll concentration [24,25] and leaf N concentration [26]. Because the fluorescence ratio is only related to chlorophyll concentration or photosynthetic activities, soil background does not affect the vegetation spectra at the early growth stages. Longchamps and Khosla [27] observed that N supply levels in corn could be differentiated as early as the V5 phenological stage using a proximal fluorescence sensor. In contrast, the test results only became reliable starting from the V8 growth stage based on reflective sensors [28,29]. Therefore, chlorophyll fluorescence sensing is a powerful solution to the shortcomings of proximal reflectance sensors in crop N status monitoring.
In addition to chlorophyll content, polyphenols (mainly flavonols) can also be used to estimate plant N status. When N is deficient, polyphenols increase significantly due to carbon and N balance regulation mechanisms [30], although an increase in polyphenol concentration may also be related to their photoprotective roles [31]. These compounds are mainly concentrated in epidermal cells and have typical absorption peaks in the ultraviolet region [31,32,33]. Thus, N status diagnosis is improved by combining the polyphenol and chlorophyll fluorescence [34,35]. Lejealle et al. [36] demonstrated that the N balance index (NBI), the ratio of chlorophyll to flavonol, had a better and more stable correlation with leaf N concentration. Leaf fluorescence sensor Dualex (FORCE-A, Orsay, Paris, France) and canopy fluorescence sensor Multiplex (FORCE-A, Orsay, Paris, France) can be used to estimate plant polyphenol contents as well as chlorophyll content. Thus, in addition to N, these fluorescence sensors can detect physiological and biochemical plant parameters such as anthocyanins [37], or plant diseases [38]. Furthermore, Multiplex is a canopy sensor that is more efficient than a leaf sensor such as Dualex. It allows rapid large-area measurements with simultaneous GPS data recording for field map generation [39,40]. However, studies based on Multiplex are still limited, especially for rice N status diagnosis and precision N management. Zhang et al. [41] identified FERARI (fluorescence excitation ratio anthocyanin relative index), SFR_G (simple fluorescence ratio under green excitation), SPAD/DUAD (Dualex reading), and SPAD reading as valuable indicators for monitoring corn N status at early stages. In addition, they found leaf-scale (LS) Multiplex measurements were better related to N treatments than the Multiplex readings made from a small distance above the canopy (AC). Li et al. [42] reported that both leaf-based SPAD and canopy Multiplex indices could be used to predict rice leaf N contents. The on-the-go (OG) measurement mode means placing the sensor probe in the canopy close to the blade and continuously collecting data during the move. This mode is most efficient in data collection. In a recent study, Diago et al. [39] confirmed the capability of the Multiplex sensor using the OG mode to estimate key nutritional parameters in grapevine leaves in motion by calibrating Multiplex against the leaf-clip Dualex sensor. However, further systematic and comprehensive study is necessary to investigate the application potential of canopy fluorescence sensing for rice N status diagnosis in motion.
Hence, the main objectives of this study were to: (1) Determine the feasibility of using canopy multispectral fluorescence sensing system such as the Multiplex sensor to estimate N status in paddy rice by comparing the results of three measurement modes (LS, AC, and OG); and (2) establish and validate the estimation models for N indicators based on the optimal Multiplex indices. In addition, to reduce the influences of varieties, years, sites, and other factors, the normalized N sufficiency index (NSI) was calculated and included in the analysis of the fluorescence indices. Well-fertilized reference plots were used to normalize reflectance measurements as more stable rice N diagnostic results might be obtained when calculating the NSI [43].

2. Materials and Methods

2.1. Experimental Design

Two field trials were conducted at the Jiansanjiang Experiment Station of China Agricultural University (47°15′N, 132°39′E), Sanjiang Plain, Heilongjiang Province, northeast China. The field experiments in 2012 and 2013 included five different N rates (0, 70, 100, 130, 160 kg N ha−1) and two Japonica rice varieties, Kongyu 131 (KY 131) and Longjing 21 (LJ 21). These two represent the main varieties in this region: KY131 has 11 leaves, four elongation nodes and about 127 maturity days, while LJ 21 has 12 leaves and needs 133 maturity days. Planting density was approximately 30 hills m−2 for KY 131 and 28 hills m−2 for LJ 21, with an identical row spacing of 0.3 m. The plot size was 4.5 m × 9 m. The N fertilizer applications were split into 40%, 30%, 30%, and applied before the transplanting, at the active tillering stage, and at the stem elongation (SE) stage, respectively. Phosphate (50 kg P2O5 ha−1) fertilizers were applied before transplanting, and potash (100 kg K2O ha−1) fertilizers were applied as two splits, 50% as basal fertilizer and 50% as panicle fertilizer at the SE stage. The two experiments were carried out in a randomized complete block design with three replicates. Table 1 lists the details of the two experiments.

2.2. Fluorescence Measurements

The portable fluorescence sensor Multiplex®3 was used in this study. It is an active sensor involving four emission light sources (UV_A, green, red or blue) to excite the fluorescence in plant tissues. Generally, the UV_A (375 nm), green (530 nm), and red (630 nm) emission light sources were used for plant monitoring while the blue (470 nm) emission light source was used for calibration. The sensor has three filtered detectors for fluorescence recording including blue–green fluorescence (447 nm) (BGF), red fluorescence (665 nm) (RF) and far-red fluorescence (735 nm) (FRF). The fluorescence measurements were not influenced by ambient light conditions [41].
All the variables provided by the Multiplex sensor and their explanations are listed in Table 2. There are nine measured single fluorescence variables under three excitations and ten calculated indices. The Simple Fluorescence Ratio (SFR) index is the ratio of the FRF and RF emission under red (SFR_R) or green (SFR_G) illumination. SFR is related to leaf chlorophyll content. Due to the chlorophyll absorption waveband overlapping with its fluorescence emission red band, the chlorophyll re-absorption occurs at the shorter red wavelength rather than at the far-red wavelength. Therefore, using the FRF as a reference, the absorption of the RF reflects the content of chlorophyll [24,44]. Accordingly, SFR increases as chlorophyll content increases. The Flavonols (FLAV) index compares the fluorescence emission density of the far-red fluorescence under ultraviolet (FRF_UV) and red excitation (FRF_R). It is related to the flavonoid concentration of the epidermal layer [45,46]. The N Balance Index (NBI) is defined as the ratio of SFR and FLAV [47]. Therefore, NBI is proportional to both chlorophyll and flavonoid concentrations. Blue–green to far-red fluorescence ratio (BRR_FRF) index is the ratio of BGF and FRF under UV excitation. The ratio of fluorescence emission at 440 nm and 740 nm (F440/F740) wavelength has been shown to be sensitive to environmental changes and growth conditions, and can detect plant stress before visible symptoms occur [18,48]. The Anthocyanins (ANTH) index and FERARI are both proven to be correlated with skin anthocyanin content [49].
Measurements were taken in motion with the sensor bottom along and just touching the crop leaves on two representative rows in the center of each plot. The data collected from the two rows were averaged to represent the plot [50]. This approach was named “measuring in motion” or “on-the-go” (OG) mode. Notably, in the OG mode, the Multiplex sensor was placed right above the crop leaves manually while walking forward as it is impractical to use a motorized vehicle or to mount the sensor on a vehicle in paddy fields. For comparison purposes, data were also measured in the LS and AC modes. The leaf-borne measurements in the LS mode were collected in the laboratory by taking ten leaves in the second position from the top. In the AC mode, the measurements were collected in the field by randomly selecting ten representative hills and placing the Multiplex sensor at a distance of approximately 10 cm above each selected hill. The average reading was then used to represent the plot.

2.3. Plant Sampling and Measurements

At the PI, SE, and HE stages, for each plot three to six plant samples were collected in each plot from the same rows where fluorescence sensor measurements were acquired. Various N status indicators, including above-ground biomass (AGB), plant N concentration (PNC), leaf N concentration (LNC), plant N uptake (PNU), and N nutrition index (NNI), were determined. The detailed sampling dates and related information are listed in Table 1. Roots from all the plant samples were removed and the samples were cleaned with water and then separated into leaves, stems, and panicles (e.g., HE stage). All samples were dried at 105 °C for half an hour to reduce plant metabolic activities. After being dried in an oven at 70–80 °C to constant weight, the samples were weighed. N concentrations for leaves, stems, and panicles were determined using the standard Kjeldahl method. The PNC is the sum of the products of each organ’s N concentration and its proportional weight. The PNU was determined by multiplying PNC by AGB. The NNI is defined as the ratio of the actual PNC ( N a ) and the critical N concentration ( N c ), with N c being calculated following the equation developed for rice in this region [51].

2.4. Statistical Analysis

The Multiplex data of the three measurement modes at each sampling stage, year, and cultivar obtained under the varied N supply were subjected to analysis of variance (ANOVA) using SAS software (SAS Institute, Cary, NC, USA). Moreover, the means for each treatment were compared using the least significant difference (LSD) test at the 95% level of significance. Relationships between the Multiplex indices and N status indicators were determined. All of the in-situ samples were divided into two groups by a stratified random sampling method, with approximately 2/3 of the data used for model calibration and the remaining for model validation. Simple linear regression analyses were performed with SPSS 20.0 (SPSS Inc., Chicago, Illinois, USA). The coefficient of determination (R2) was calculated for comparison. The relationships between Multiplex indices and N status indicators established at different growth stages were validated. The root mean square error (RMSE) and the relative error (RE), shown in Equations (1) and (2), between the predicted and observed values were used for evaluation.
RMSE = 1 n i = 1 n ( y i y ^ i ) 2
RE ( % ) = R M S E y ¯ × 100
where y i , y ^ i , and y ¯ were the observed, predicted, and mean value, respectively.
In order to evaluate normalized vegetation indices for improving the estimation of N nutrition indicators, the well-fertilized plots were used as N sufficient reference to calculate the NSI index. The NSI index equals to the ratio of Multiplex indices of the plots receiving normal N rates and the well-fertilized plots. In this study, the treatment with the largest shoot dry matter was defined as the well-fertilized plot, corresponding to the treatment of 130 or 160 kg N ha−1.
Finally, the NNI diagnostic results of validation data using Multiplex indices were compared to the observed NNI by areal agreement and Kappa statistics [52]. Both used the same diagnostic criteria: N was deficient when NNI < 0.95, N was optimal when NNI is between 0.95–1.05, and N was in surplus when NNI > 1.05. The areal agreement (%) and Kappa statistics were used to determine the accuracy of the diagnostic results. The areal agreement indicates the percentage of two groups sharing a common category or diagnostic class [53]. The Kappa statistic is a more robust measure of the agreement of two classifications by correcting the agreement that occurs by chance. When Kappa equals 1, it indicates that the two categorization systems are identical. Kappa ≥ 0.60 indicates a satisfactory agreement, while the Kappa < 0.4 indicates weak agreement [54].

3. Results

3.1. Comparison of the Three Measurement Modes

The three aforementioned measurement modes were applied for each treatment plot. The results were compared to determine the best measurement mode. The abbreviation of the measurement mode is added to the variable as a prefix. For example, AC_SFR_G represents the Multiplex index SFR_G obtained from above the canopy.
Figure 1 shows box plots of Multiplex indices obtained from the three measuring modes at two phenological stages. Since FLAV and ANTH are the log transformation values of FER_RUV and FER_RG, only FLAV and ANTH were selected for the analysis to avoid duplication. During both the PI and SE growth stages, the mean values of the Multiplex indices (except for the BRR_FRF) obtained using the OG mode were significantly higher (P ≤ 0.05) than the leaf-borne indices measured in the LS mode, whereas the latter were significantly higher than the measurements obtained in the AC mode. In addition, substantial large data ranges were shown in the NBI_G and NBI_R indices measured in motion and the leaf-borne BRR_FRF value.
The ANOVA results are listed in Table 3, indicating the sensitivities of the Multiplex indices to five N rates (0, 70, 100, 130, and 160 kg N ha−1). N fertilization treatment effects were significant for most of the Multiplex variables measured in the OG mode whereas the opposite was true for the leaf-borne measurements. Comparatively, the N treatment effects were more significant for the readings obtained in the OG mode than those collected in the AC mode, while the measurements collected using the LS mode showed the minimum sensitivity to N rate changes. As an example, Figure 2 further demonstrates the comparison of the selected Multiplex indices at each N application rate for cultivar KY 131 in 2013. From Figure 2 we can also see that the indices measured using the OG mode could better distinguish the effects of N supply compared to the other two modes. Especially, the Multiplex indices measured in the OG mode performed consistently better than the others for differentiating high N-application treatment (≥100 kg ha−1) effects. Similar results were achieved for the variety LJ 21, or for the experiments conducted in 2012.
In conclusion, the ANOVA analysis results showed that most of the Multiplex variables obtained in the OG mode were more sensitive to N supply, followed by the AC modes, while measurements made in the LS mode were the least influenced by N supply (Table 3, Figure 2). Thus, only the results of the OG mode were used for further analysis and discussion.

3.2. Changes in Multiplex Indices (“On-The-Go” Mode) over Growth Stages under Different N Supplies

The SFR_G, SFR_R, NBI_G, and NBI_R indices demonstrated an increasing trend as N rate increased, while a decreasing trend was shown for FLAV (Figure 2). Comparatively, the ANTH and BRR_FRF values were less sensitive to the changes in N rates. The values of these SFR_G, SFR_R, NBI_G, and NBI_R indices increased from the PI to SE stage but decreased slightly from the SE to HE stage, because the panicle formation decreased the chlorophyll/N concentration in the upper layer at the HE stage. The opposite was true for the BRR_FRF, FLAV, and ANTH. NBI_G and NBI_R could differentiate different N application rates the best regardless of the growth stages, followed by SFR_G, SFR_R, BRR_FRF, and FLAV. The performance of ANTH was the worst (Figure 2).

3.3. Correlations between Multiplex Indices (“On-The-Go” Mode) and N Status Indicators

The linear regression results of the seven Multiplex indices and the five N status indicators at three growth stages across the two rice varieties are shown in Table 4. The SFR_G, SFR_R, NBI_G, and NBI_R indices were positively correlated with the N indicators whereas the BRR_FRF, ANTH, and FLAV were inversely associated with them. The R2 of the regression models based on these indices varied from 0.03 to 0.78. The best performing index varied at different stages, but NBI_G and NBI_R showed steady high correlations with all five N status indicators. The second-best performing indices were BRR_FRF and FLAV. The SFR_G and SFR_R indices displayed high or moderate correlations with the N indicators during the PI or HE stage, respectively.
Compared to the counterpart of the standard indices, the normalized sufficiency indices SFR_GNSI, SFR_RNSI, and ANTHNSI exhibited better linear relationships with LNC, PNC, and NNI in most of the cases, especially at the SE and HE stages. The NBI_GNSI and NBI_RNSI displayed enhanced relationships with the LNC and PNC at the PI and HE stages, and with the NNI at the HE stage. The BRR_FRFNSI, FLAVNSI showed improved associations with PNC at the PI and HE stages, and with NNI at the SE and HE stages. All the standard indices showed moderate–high relationships with the AGB and PNU during the PI and HE stages, while at the HE stage, greatly improved R2 values were obtained using the normalized indices.

3.4. Validation of the Estimation Models for N Status Indicators

In order to diagnose rice N status, linear regression models between the Multiplex indices and N indicators were established. The regression models varied across growth stages. Table 5 lists the best performing models at the PI, SE, and HE stages. The best performing indices differed across the stages. However, the relationships of NBI_G and NBI_R with N indicators were relatively more stable. After normalization, the NBI_RNSI showed an absolute advantage for N status estimation at the PI and SE growth stages, while the FLAVNSI demonstrated to be optimal for estimating most of the N indicators.
Figure 3 shows the RE values of the validation models for six Multiplex indices (SFR_G, SFR_R, BRR_FRF, FLAV, NBI_G, and NBI_R) and the N status indicators. The RE values for AGB and PNU estimations based on these indices decreased steadily with advancing growth stages, while a slightly increasing trend was observed for the LNC and PNC estimation models from the SE to HE stage. The RE values for LNC (4.50%–10.24%) and PNC (5.87%–10.87%) models were much smaller than those for AGB (15.49%–30.18%) and PNU (19.31%–31.25%), while the REs of NNI remained similar during the three growth stages. At the earlier to middle growth stages, NBI_R and NBI_G presented a lower RE than the other four indices for all the five N indicators. At the HE stage, however, the prediction accuracies of the six indices were similar.

3.5. Rice N Status Diagnosis

The best performing indices including SFR_G, BRR_FRF, NBI_G, NBI_R, NBI_GNSI, NBI_RNSI and FLAVNSI were validated using independent data sets (Table 5). Moderate–high model performance with R2 ranging from 0.34 to 0.82 was observed, especially for the LNC, PNC, and NNI estimations. The areal agreement and Kappa statistics were compared at the critical N fertilizer application stages (SE and HE) to evaluate the N diagnostic accuracies of the indices. Results confirmed that the NNI models based on NBI_R and NBI_G performed consistently well at the SE and HE growth stages, and their corresponding NSI indices further improved the results (Table 6). At the SE stage, the NBI_RNSI achieved the highest diagnostic accuracy (areal agreement = 90%; Kappa = 0.84), while the best accuracy was achieved by FLAV at the HE stage (areal agreement = 90%; Kappa = 0.76). In addition, across the two growth stages, the NBI_RNSI showed the highest diagnostic consistency, followed by the BRR_FRFNSI.

4. Discussion

4.1. Multiplex Measurement Modes and Estimation of Crop N Indicators by Fluorescence Indices

The N treatment effects were more significant for the readings obtained in the OG mode than those collected in the AC mode (Table 3, Figure 2), which is different from the finding by Diago et al. [39] who reported a 20% loss of information occurred when using the Multiplex on-the-go (compared to the AC mode) for N assessment of grapevine. This is because the OG measurements in this study were taken manually by placing the Multiplex sensor right on the top of the rice rather than a small distance above the rice canopy while passing through the rice paddy. In contrast, in the study by Diago et al. [39], the Multiplex sensor was mounted onto an all-terrain vehicle and placed 1.5 m above the ground so that the leaves on the mid-part of the canopy were automatically measured at a 20 cm distance, the same measuring distance as their AC mode. In addition, this study revealed that measurements made using the LS mode were the least sensitive to N supply, contrasting to the result by Zhang et al. [41] who found Multiplex measurements made from corn leaves were more capable of distinguishing plant N status than those made from above the plants. While the leaf scale measurements made in this study were collected in the laboratory by taking ten leaves in the second position from the top, the leaf-borne measurements by Zhang et al. [41] were made on 20 representative plants in the center two rows of each plot in the field, which is more similar to the OG rather than the LS method in this study. Another possible reason for their better results with the LS method is that the OG measurements do not give much time to choose leaves and result in more random leaf choosing than the LS method, which may have an unwanted tendency to choose “good” leaves. This is particularly true for maize, because systematic use of a representative leaf is easier, as it is well known which is the most representative leaf for each growth stage, given its determined growth. For rice, the individual leaves are quite small and the signal obtained during the measurements taken on a leaf is relatively weak and can be easily affected by other factors. One advantage the “measuring in motion” or OG mode has is efficiency, especially when the sensor is mounted on a vehicle or other automatic devices, which might make practical applications of such non-destructive technology over large areas possible. Bringing the sensor close to or even touching the leaves of the crop in OG mode may help reduce information loss. However, further well-designed studies are needed to confirm this finding.
Strong relationships between the Multiplex indices (SFR_G, SFR_R, BRR_FRF, FLAV, NBI_G, and NBI_R) measured in the OG mode and the five N indicators were achieved with low RE and high R2 values (Table 4, Figure 3). This finding conforms to previous research results in this field [26,42,55,56]. Many studies confirmed SFR was a good fluorescence index for chlorophyll content monitoring [24,38,57]. However, in this study, it was found that the R2 of the SFR_G, SFR_R for LNC, PNC, and NNI estimation decreased steadily from the early-stage to later stages, while an opposite trend was observed for FLAV. Padilla et al. [55] found that the relationships between the NNI and SFR_G changed with the phenological stages of cucumber (Cucumis sativus L.). Firstly, the consistency of the relationship between chlorophyll content and N concentration varied with crop development, leading to different performances of SFR for N concentration estimation. For example, the linear correlation between LNC and chlorophyll meter readings of rice was weaker at the SE stage than at other growth stages [58]. Secondly, the weaker differentiation ability of SFR under the unlimited N conditions may also be a reason [55]. The performance of FLAV increased from the PI to HE stages, which was confirmed by Padilla et al. [59], who found the relationship of FLAV and NNI increased at the middle to late growth stages. The better performance of FLAV at the later stage may be attributed to the accumulation of the flavonols content in leaves under light radiation [44,60]. The NBI_G and NBI_R indices were shown as the best indices for estimating the N indicators (Table 5). Many studies have proven that the NBI indices appeared to be the most efficient in estimating the N status [47,56,59,61]. This is because the NBI is a ratio of SFR and FLAV, which makes it more robust than using FLAV or SFR alone to reduce the effects of leaf age or other factors [34,36,47]. The NBI_G and NBI_R had similar performance in this study, as demonstrated by Longchamps and Khosla [27]. Moreover, Longchamps and Khosl [27] found that SFR was less sensitive to N application than NBI, which conforms to our results, as shown in Figure 2. In most cases, the SFR_G and SFR_R indices could not distinguish between the 100 and 130 kg N ha−1 treatments, but NBI could. The BRR_FRF index was significantly correlated with the N nutritional status and was especially sensitive to N deficiency in this research. When there is N stress, the fluorescence ratio of blue–green/far-red will increase after exposure to elevated UV radiation to avoid or alleviate the damage of the photosynthetic apparatus [62]. Generally, the UV-protection response takes place before the chlorophyll damage can be seen, so the BRR_FRF can also be considered as a potential index that can realize early N deficiency detection [63]. The BRR_FRF was also very sensitive to environmental stresses, such as disease and drought [38,63,64]. The ANTH index provided by Multiplex is commonly used to reflect anthocyanin content, which corresponds to the maturation degree of fruit [49,65]. In this research, the low values of ANTH were due to the low anthocyanin content in the rapid vegetative growth phase for rice [66]. Nevertheless, ANTH was also found to be closely related to the leaf chlorophyll concentration in some studies [38,41]. This study revealed that ANTH was significantly related to N status indicators in PI and HE growth stages with moderate R2 values (Table 4).

4.2. Normalized Nitrogen Sufficiency Fluorescence Indices

Our research involved two years and two varieties of experiments. In these experiments, N fertilization rate is the main variable, and the variation is so high that it will probably override any other source of variation. In a commercial field, many factors can influence N availability, N and Chlorophyll relationship, or Chlorophyll (+FLAV and ANTH) relations to fluorescence indices, including biotic or abiotic stresses at the moment of measurement or in the history of the crop or even the field. The normalized N sufficiency index approach has been suggested to reduce the influence of the varieties, developmental stages, and other variables on SPAD values or spectral data [3,11,67]. From the results of this study, in most cases, the normalized NSIs were better associated with the LNC, PNC, and NNI (Table 4). The R2 of the ANTHNSI was improved the most, followed by the NBI_GNSI and NBI_RNSI. However, the improvement in R2 for BRR_FRFNSI was minimal. The variance analysis of this study showed consistent results, which demonstrated that NSI indices could reduce the influence of inter-annual and growth stage differences. Since NNI itself is a diagnostic criterion, it represents an optimal N status when NNI is equal to one [68]. Most of the NSI indices greatly improved the NNI diagnostic accuracy at the critical topdressing (SE and HE) stages (Table 6). Similarly, Lu et al. [43] observed that the NNI inversion through the normalized vegetation indices further improved the N nutrition diagnostic results of rice.
Hussain et al. [69] proposed a critical NSI value of 0.90 for rice. However, in this study, when the NSI indices were 0.90, different optimal NNI values, ranging from 0.85 to 2.14, were derived by different indices at different developmental stages. Only the corresponding optimal NNI values for the NBI_GNSI and NBI_RNSI indices were close to one (ranging from 0.91 to 1.19). Therefore, to avoid the risk of misdiagnosis, the NSI threshold was not used as a diagnostic criterion directly. A possible reason for this is that the N fertilizer application rate in this study was only 1.3–1.6 times higher than the optimal amount instead of 1.8–2.0 times higher than recommended for the well-fertilized N plot as Hussain et al. [69] suggested. Furthermore, all of the Multiplex indices were divided by the readings of the N rate with the largest shoot dry matter at each sampling date to obtain a sufficiency index. However, Varvel et al. [70] suggested that the maximum readings within each cropping system, variety, and year should be considered as the normalized criterion. Obviously, with different normalization criteria, different sufficiency indices will be obtained, which will affect the corresponding analysis results. Another limitation of the NSI approach is that well-fertilized reference plots need to be established in each farmer’s field for practical application purposes, and some farmers may not be willing to do this. More in-depth and systematic research is expected in the future.

4.3. The Application Potential of the Multiplex Sensor

The Multiplex indices presented good R2 values for LNC and PNC estimations at the earlier growth stages (Table 4). In particular, the validation data showed that the RE values for LNC and PNC estimations were as low as 6%–7% (Figure 3c,d). This is consistent with the results of Cerovic et al. [71] and Agati et al. [56], who have shown a high correlation between the fluorescence index and LNC. NBI and LNC had a fairly linear relationship. Therefore, the NBI indices can be used to more accurately estimate a wider range of LNC. Agati et al. [56] also found the results based on reflectance imaging (camera picture) are less sensitive to N application than fluorescence-based indices. Research by Stroppiana et al. [72] and Yu et al. [19] on rice showed unsatisfactory results for the estimation of LNC and PNC based on reflectance spectroscopy. This is possibly due to the fact that the effect of N on the leaf area index and biomass is much greater than its effect on chlorophyll content. Second, near-infrared radiation is hardly absorbed in the canopy and is highly transmissive, so its correlation with leaf area index or biomass is extremely high; while visible light, especially the blue and red radiation, is easily absorbed by chlorophyll and its transmittance is low, so it is highly correlated with chlorophyll content [15,72]. On the other hand, changes in plant metabolism indicators are fast or slow due to changes in response to the environment. However, the sensitivity of reflectance-based parameters does not always provide satisfactory monitoring results [73]. Demotes-Mainard et al. [74] observed that changes in N concentration took precedence over changes in biomass. Thus, fluorescence-based techniques that are highly sensitive to plant N status information may address the limitation of reflectance-based methods [27,73]. Similarly, the Multiplex indices, especially the NBI_G and NBI_R, presented accurate estimation for NNI, with R2 reaching a maximum of 0.72–0.78, and the validation results also showed a low inversion error for NNI (RMSE ≤ 0.16, RE ≤ 15%) (Table 4). Many studies have confirmed that NBI has a strong estimation potential for NNI [47,55,59]. This is because NBI is the ratio of SFR to FLAV. The SFR index was considered to be an important parameter for estimating chlorophyll concentration, which was often used as an index of surface-based N [75], while the FLAV parameter directly reflects flavonol content, which is controlled by light as well as leaf mass per area, and has a strong correlation with leaf mass [76]. Therefore, NBI as the SFR/FLAV ratio is the best N nutrition diagnostic index.
The fluorescence-based indices are more sensitive to chlorophyll or N content than the reflectance-based indices and can detect the difference in N nutrition status earlier. However, the difference of the stage-based models between the indices and the N nutrition indicators based on the canopy reflectance instrument is smaller than that based on fluorescence [59]. The surface area of the crop involved in each test when using the canopy reflectance spectroscopy sensors is larger than the fluorescence sensors [59]. Therefore, canopy reflectance measurements are more representative, while fluorescence instruments require increasing the number of tests to obtain sufficiently representative data. Although the performance of the fluorescence sensor was quite good for estimating LNC, PNC, and NNI, 22%–60% of their variability was still not explained. In addition, the fluorescence sensor did not perform very well for estimating plant biomass in the middle to late stages (Table 4). It has been suggested to combine the fluorescence and reflectance data to improve the estimation of plant N status [3,73]. This may be one of the important research directions in the future.

5. Conclusions

This research compared the LS, OG, and AC measurement modes of the fluorescence instrument Multiplex®3 and determined that the OG mode was best suited for this rice N status study. Using the OG mode, stable test results and crop growth information were derived. The results revealed that the fluorescence indices of NBI, SFR, BRR_FRF, and FLAV were significantly correlated to all five N status indicators from the PI through HE growth stages. Among them, NBI_G and NBI_R were the best performing indices and highly correlated to LNC (R2 = 0.52–0.68), PNC (R2 = 0.52–0.71), NNI (R2 = 0.69–0.78), AGB (R2 = 0.47–0.64), and PNU (R2 = 0.68–0.72) at the three growth stages. The normalized sufficiency indices of the Multiplex parameters could greatly improve the LNC, PNC, and NNI estimation ability, especially at the HE stage. The N diagnostic results indicated that the NBI_RNSI and FLAV achieved the highest diagnostic accuracy rate (90%) at the SE and HE stage, respectively, while NBI_RNSI showed the highest analytical consistency across growth stages. The results suggest that the Multiplex sensor can be used to reliably estimate N nutritional status for rice in cold regions, especially for the estimation of LNC, PNC, and NNI. The normalized sufficiency indices based on Multiplex indices may further improve the accuracy of N nutrition diagnosis by reducing the differences between years and varieties.

Author Contributions

Y.M. and G.B. conceived and guide the study. S.H. and Q.C. conducted the field experiments. S.H. and H.Y. performed the data analysis. S.H. wrote the original manuscript. Y.M., F.Y. and V.I.S.L.-W. revised the manuscript. G.B. and H.Y. reviewed and edited the manuscript. All authors read and approved the final manuscript.

Funding

This research was financially supported by the National Key Research and Development Program of China (2016YFD0200600, 2016YFD0200602), National Basic Research Program (2015CB150405), and the Norwegian Ministry of Foreign Affairs (SINOGRAIN II, CHN-17/0019).

Acknowledgments

We would like to thank the supports by Wen Yang, Huamin Zhu, and Fengyan Liu at the Jiansanjiang Institute of Agricultural Sciences. We also would like to thank Jianning Shen, Weifeng Yu, and Shanshan Cheng for their fieldwork and contributions in spectral data collection.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Guo, J.H.; Liu, X.J.; Zhang, Y.; Shen, J.L.; Han, W.X.; Zhang, W.F.; Christie, P.; Goulding, K.W.; Vitousek, P.M.; Zhang, F.S. Significant acidification in major Chinese croplands. Science 2010, 327, 1008–1010. [Google Scholar] [CrossRef]
  2. Miao, Y.; Stewart, B.A.; Zhang, F. Long-term experiments for sustainable nutrient management in China. A review. Agron. Sustain. Dev. 2011, 31, 397–414. [Google Scholar] [CrossRef]
  3. Samborski, S.M.; Tremblay, N.; Fallon, E. Strategies to make use of plant sensors-based diagnostic information for nitrogen recommendations. Agron. J. 2009, 101, 800–816. [Google Scholar] [CrossRef]
  4. Evans, J.R. Nitrogen and photosynthesis in the flag leaf of wheat (Triticum aestivum L.). Plant Physiol. 1983, 72, 297–302. [Google Scholar] [CrossRef] [PubMed]
  5. Blackmer, T.M.; Schepers, J.S. Use of a chlorophyll meter to monitor nitrogen status and schedule fertigation for corn. J. Prod. Agric. 1995, 8, 56–60. [Google Scholar] [CrossRef]
  6. Schlemmer, M.R.; Francis, D.D.; Shanahan, J.F. Remotely measuring chlorophyll content in corn leaves with differing nitrogen levels and relative water content. Agron. J. 2005, 97, 106–112. [Google Scholar] [CrossRef]
  7. Schröder, J.J.; Neeteson, J.J.; Oenema, O.; Struik, P.C. Does the crop or the soil indicate how to save nitrogen in maize production? Reviewing the state of the art. Field Crop Res. 2000, 66, 151–164. [Google Scholar] [CrossRef]
  8. Schepers, J.S.; Blackmer, T.M.; Francis, D.D. Predicting N fertilizer needs for corn in humid regions: Using chlorophyll meters. In Predicting N Fertilizer Needs for Corn in Humid Regions; Bock, B.R., Kelley, K.R., Eds.; National Fertilizer and Environmental Research Center: Muscle Shoals, AL, USA, 1992; pp. 105–114. [Google Scholar]
  9. Schepers, J.S.; Francis, D.D.; Vigil, M.; Below, F.E. Comparison of corn leaf nitrogen concentration and chlorophyll meter readings. Commun. Soil Sci. Plan. 1992, 23, 2173–2187. [Google Scholar] [CrossRef]
  10. Markwell, J.; Osterman, J.C.; Mitchell, J.L. Calibration of the Minolta SPAD-502 leaf chlorophyll meter. Photosynth. Res. 1995, 46, 467–472. [Google Scholar] [CrossRef]
  11. Lin, F.F.; Qiu, L.F.; Deng, J.S.; Shi, Y.Y.; Chen, L.S.; Wang, K. Investigation of spad meter-based indices for estimating rice nitrogen status. Comput. Electron. Agric. 2010, 71, S60–S65. [Google Scholar] [CrossRef]
  12. Ali, M.M.; Al-Ani, A.; Eamus, D.; Tan, D.K.Y. Leaf nitrogen determination using non-destructive techniques—A review. J. Soil Sci. Plant Nutr. 2017, 40, 928–953. [Google Scholar]
  13. Mulla, D.J.; Miao, Y. Precision Farming. In Land Resources Monitoring, Modeling, and Mapping with Remote Sensing; Thenkabail, P.S., Ed.; CRC Press: Boca Raton, FL, USA, 2016. [Google Scholar]
  14. Gitelson, A.A.; Gritz, Y.; Merzlyak, M.N. Relationships between leaf chlorophyll content and spectral reflectance and algorithms for non-destructive chlorophyll assessment in higher plant leaves. J. Plant Physiol. 2003, 160, 271–282. [Google Scholar] [CrossRef] [PubMed]
  15. Heege, H.J.; Reusch, S.; Thiessen, E. Prospects and results for optical systems for site-specific on-the-go control of nitrogen-top-dressing in Germany. Precis. Agric. 2008, 9, 115–131. [Google Scholar] [CrossRef]
  16. Yao, Y.; Miao, Y.; Huang, S.; Gao, L.; Ma, X.; Zhao, G.; Jiang, R.; Chen, X.; Zhang, F.; Yu, K.; et al. Active canopy sensor-based precision n management strategy for rice. Agron. Sustain. Dev. 2012, 32, 925–933. [Google Scholar] [CrossRef]
  17. Cao, Q.; Miao, Y.; Shen, J.; Yu, W.; Yuan, F.; Cheng, S.; Huang, S.; Wang, H.; Yang, W.; Liu, F. Improving in-season estimation of rice yield potential and responsiveness to topdressing nitrogen application with Crop Circle active crop canopy sensor. Precis. Agric. 2016, 17, 136–154. [Google Scholar] [CrossRef]
  18. Olfs, H.W.; Blankenau, K.; Brentrup, F.; Jasper, J.; Link, A.; Lammel, J. Soil-and plant-based nitrogen-fertilizer recommendations in arable farming. J. Soil Sci. Plant Nutr. 2005, 168, 414–431. [Google Scholar] [CrossRef]
  19. Yu, K.; Li, F.; Gnyp, M.L.; Miao, Y.; Bareth, G.; Chen, X. Remotely detecting canopy nitrogen concentration and uptake of paddy rice in the Northeast China Plain. ISPRS J. Photogramm. Remote Sens. 2013, 78, 102–115. [Google Scholar] [CrossRef]
  20. Huang, S.; Miao, Y.; Zhao, G.; Yuan, F.; Ma, X.; Tan, C.; Yu, W.; Gnyp, M.; Lenz-Wiedemann, V.; Rascher, U. Satellite remote sensing-based in-season diagnosis of rice nitrogen status in Northeast China. Remote Sens. 2015, 7, 10646–10667. [Google Scholar] [CrossRef]
  21. Bredemeier, C.; Schmidhalter, U. Laser-induced chlorophyll fluorescence sensing to determine biomass and nitrogen uptake of winter wheat under controlled environment and field condition. In Proceedings of the 5th European Conference on Precision Agriculture, Uppsala, Sweden, 9–12 June 2005; Wageningen Academic Publishers: Wageningen, The Netherlands, 2005; pp. 273–280. [Google Scholar]
  22. McMurtrey, J.E., III; Chappelle, E.W.; Kim, M.S.; Meisinger, J.J.; Corp, L.A. Distinguishing nitrogen fertilization levels in field corn (Zea mays L.) with actively induced fluorescence and passive reflectance measurements. Remote Sens. Environ. 1994, 47, 36–44. [Google Scholar] [CrossRef]
  23. Langsdorf, G.; Buschmann, C.; Sowinska, M.; Babani, F.; Mokry, M.; Timmermann, F.; Lichtenthaler, H.K. Multicolour fluorescence imaging of sugar beet leaves with different nitrogen status by flash lamp UV-excitation. Photosynthetica 2000, 38, 539–551. [Google Scholar] [CrossRef]
  24. Gitelson, A.A.; Buschmann, C.; Lichtenthaler, H.K. The chlorophyll fluorescence ratio F-735/F-700 as an accurate measure of the chlorophyll content in plants. Remote Sens. Environ. 1999, 69, 296–302. [Google Scholar] [CrossRef]
  25. Cerovic, Z.G.; Goutouly, J.P.; Hilbert, G.; Destrac-Irvine, A.; Martinon, V.; Moise, N. Mapping winegrape quality attributes using portable fluorescence-based sensors. Frutic 2009, 9, 301–310. [Google Scholar]
  26. Yang, J.; Wei, G.; Shi, S.; Lin, D.; Sun, J.; Song, S.; Chen, B.; Zhang, Z. Analyzing the performance of fluorescence parameters in the monitoring of leaf nitrogen content of paddy rice. Sci. Rep. 2016, 6, 28787. [Google Scholar] [CrossRef] [PubMed]
  27. Longchamps, L.; Khosla, R. Early detection of nitrogen variability in maize using fluorescence. Agron. J. 2014, 106, 511–518. [Google Scholar] [CrossRef]
  28. Teal, R.K.; Tubana, B.; Girma, K.; Freeman, K.W.; Arnall, D.B.; Walsh, O.; Raun, W.R. In-season prediction of corn grain yield potential using normalized difference vegetation index. Agron. J. 2006, 98, 1488–1494. [Google Scholar] [CrossRef]
  29. Martin, K.L.; Girmaa, K.; Freemana, K.W.; Teala, R.K.; Tubańaa, B.; Arnalla, D.B.; Chunga, B.; Walsha, O.; Solieb, J.B.; Stoneb, M.L.; et al. Expression of variability in corn as influenced by growth stage using optical sensor measurements. Agron. J. 2007, 99, 384–389. [Google Scholar] [CrossRef]
  30. Jones, C.G.; Hartley, S.E. A protein competition model of phenolic allocation. Oikos 1999, 86, 27–44. [Google Scholar] [CrossRef]
  31. Burchard, P.; Bilger, W.; Weissenbock, G. Contribution of hydroxycinnamates and flavonoids to epidermal shielding of UV-A and UV-B radiation in developing rye primary leaves as assessed by ultraviolet-induced chlorophyll fluorescence measurements. Plant Cell Environ. 2000, 23, 1373–1380. [Google Scholar] [CrossRef]
  32. Knogge, W.; Weissenböck, G. Tissue-distribution of secondary phenolic biosynthesis in developing primary leaves of Avena sativa L. Planta 1986, 167, 196–205. [Google Scholar] [CrossRef]
  33. Cerovic, Z.G.; Ounis, A.; Cartelat, A.; Latouche, G.; Goulas, Y.; Meyer, S.; Moya, I. The use of chlorophyll fluorescence excitation spectra for the non-destructive in situ assessment of UV-absorbing compounds in leaves. Plant Cell Environ. 2002, 25, 1663–1676. [Google Scholar] [CrossRef]
  34. Tremblay, N.; Wang, Z.; Bélec, C. Evaluation of the Dualex for the assessment of corn nitrogen status. J. Soil Sci. Plant Nutr. 2007, 30, 1355–1369. [Google Scholar] [CrossRef]
  35. Tremblay, N.; Wang, Z.; Belec, C. Performance of Dualex in spring wheat for crop nitrogen status assessment, yield prediction and estimation of soil nitrate content. J. Soil Sci. Plant Nutr. 2009, 33, 57–70. [Google Scholar] [CrossRef]
  36. Lejealle, S.; Evain, S. Multiplex: A new diagnostic tool for management of nitrogen fertilization of turgrass. In Proceedings of the 10th International Conference on Precision Agriculture, Denver, CO, USA, 18–21 July 2010. [Google Scholar]
  37. Cerovic, Z.G.; Moise, N.; Agati, G.; Latouche, G.; Ghozlen, N.B.; Meyer, S. New portable optical sensors for the assessment of winegrape phenolic maturity based on berry fluorescence. J. Food Compos. Anal. 2008, 21, 650–654. [Google Scholar] [CrossRef]
  38. Yu, K.; Leufen, G.; Hunsche, M.; Noga, G.; Chen, X.; Bareth, G. Investigation of leaf diseases and estimation of chlorophyll concentration in seven barley varieties using fluorescence and hyperspectral indices. Remote Sens. 2013, 6, 64–86. [Google Scholar] [CrossRef]
  39. Diago, M.P.; Rey-Carames, C.; Le Moigne, M.; Fadaili, E.M.; Tardaguila, J.; Cerovic, Z.G. Calibration of non-invasive fluorescence-based sensors for the manual and on-the-go assessment of grapevine vegetative status in the field. Aust. J. Grape Wine R. 2016, 22, 438–449. [Google Scholar] [CrossRef]
  40. Song, X.; Yang, G.; Yang, C.; Wang, J.; Cui, B. Spatial variability analysis of within-field winter wheat nitrogen and grain quality using canopy fluorescence sensor measurements. Remote Sens. 2017, 9, 237. [Google Scholar] [CrossRef]
  41. Zhang, Y.P.; Tremblay, N.; Zhu, J.J. A first comparison of Multiplex® for the assessment of corn nitrogen status. J. Food Agric. Environ. 2012, 10, 1008–1016. [Google Scholar]
  42. Li, J.W.; Zhang, J.X.; Zhao, Z.; Lei, X.D.; Xu, X.L.; Lu, X.X.; Weng, D.L.; Gao, Y.; Cao, L.K. Use of fluorescence-based sensors to determine the nitrogen status of paddy rice. J. Agric. Sci. 2013, 151, 862–871. [Google Scholar] [CrossRef]
  43. Lu, J.; Miao, Y.; Wei, S.; Li, J.; Yuan, F. Evaluating different approaches to non-destructive nitrogen status diagnosis of rice using portable RapidSCAN active canopy sensor. Sci. Rep. 2017, 7, 14073. [Google Scholar] [CrossRef]
  44. Pedrós, R.; Goulas, Y.; Jacquemoud, S.; Louis, J.; Moya, I. FluorMODleaf: A new leaf fluorescence emission model based on the PROSPECT model. Remote Sens. Environ. 2010, 114, 155–167. [Google Scholar] [CrossRef] [Green Version]
  45. Ounis, A.; Cerovic, Z.G.; Briantais, J.M.; Moya, I. Dual-excitation FLIDAR for the estimation of epidermal UV absorption in leaves and canopies. Remote Sens. Environ. 2001, 76, 33–48. [Google Scholar] [CrossRef]
  46. Agati, G.; Cerovic, Z.G.; Pinelli, P.; Tattini, M. Light-induced accumulation of ortho-dihydroxylated flavonoids as non-destructively monitored by chlorophyll fluorescence excitation techniques. Environ. Exp. Bot. 2011, 73, 3–9. [Google Scholar] [CrossRef]
  47. Cartelat, A.; Cerovic, Z.G.; Goulas, Y.; Meyer, S.; Lelarge, C.; Prioul, J.L.; Barbottin, A.; Jeuffroy, M.H.; Gate, P.; Agati, G. Optically assessed contents of leaf polyphenolics and chlorophyll as indicators of nitrogen deficiency in wheat (Triticum aestivum L.). Field Crop. Res. 2005, 91, 35–49. [Google Scholar] [CrossRef]
  48. Lichtenthaler, H.K. Vegetation stress: An introduction to the stress concept in plants. J Plant Physiol. 1996, 148, 4–14. [Google Scholar] [CrossRef]
  49. Ghozlen, N.B.; Cerovic, Z.G.; Germain, C.; Toutain, S.; Latouche, G. Non-destructive optical monitoring of grape maturation by proximal sensing. Sensors 2010, 10, 10040–10068. [Google Scholar] [CrossRef] [PubMed]
  50. Huang, S.; Miao, Y.; Yuan, F.; Cao, Q.; Ye, H.; Lenz-Wiedemann, V.; Khosla, R.; Bareth, G. Proximal fluorescence sensing for in-season diagnosis of rice nitrogen status. Adv. Anim. Biosci. 2017, 8, 343–348. [Google Scholar] [CrossRef]
  51. Huang, S.; Miao, Y.; Cao, Q.; Yao, Y.; Zhao, G.; Yu, W.; Shen, J.; Yu, K.; Bareth, G. Critical nitrogen dilution curve for rice nitrogen status diagnosis in Northeast China. Pedosphere 2018, 28, 814–822. [Google Scholar] [CrossRef]
  52. Xia, T.; Miao, Y.; Wu, D.; Shao, H.; Khosla, R.; Mi, G. Active optical sensing of spring maize for in-season diagnosis of nitrogen status based on nitrogen nutrition index. Remote Sens. 2016, 8, 605. [Google Scholar] [CrossRef]
  53. Bausch, W.C.; Khosla, R. QuickBird satellite versus ground-based multi-spectral data for estimating nitrogen status of irrigated maize. Precis. Agric. 2010, 11, 274–290. [Google Scholar] [CrossRef]
  54. Landis, R.J.; Koch, G.G. The measurement of observer agreement for categorical data. Biometrics 1977, 33, 159–174. [Google Scholar] [CrossRef]
  55. Padilla, F.M.; Peña-Fleitas, M.T.; Gallardo, M.; Thompson, R.B. Proximal optical sensing of cucumber crop N status using chlorophyll fluorescence indices. Eur. J. Agron. 2016, 73, 83–97. [Google Scholar] [CrossRef]
  56. Agati, G.; Foschi, L.; Grossi, N.; Volterrani, M. In field non-invasive sensing of the nitrogen status in hybrid bermudagrass (Cynodon dactylon, × C. transvaalensis, Burtt Davy) by a fluorescence-based method. Eur. J. Agron. 2015, 63, 89–96. [Google Scholar] [CrossRef]
  57. Buschmann, C. Variability and application of the chlorophyll fluorescence emission ratio red/far-red of leaves. Photosynth. Res. 2007, 92, 261–271. [Google Scholar] [CrossRef] [PubMed]
  58. Yuan, Z.; Ata-Ul-Karim, S.T.; Cao, Q.; Lu, Z.; Cao, W.; Zhu, Y.; Liu, X. Indicators for diagnosing nitrogen status of rice based on chlorophyll meter readings. Field Crop. Res. 2016, 185, 12–20. [Google Scholar] [CrossRef]
  59. Padilla, F.M.; Peña-Fleitas, M.T.; Gallardo, M.; Thompson, R.B. Evaluation of optical sensor measurements of canopy reflectance and of leaf flavonols and chlorophyll contents to assess crop nitrogen status of muskmelon. Eur. J. Agron. 2014, 58, 39–52. [Google Scholar] [CrossRef]
  60. Barthod, S.; Cerovic, Z.; Epron, D. Can dual chlorophyll fluorescence excitation be used to assess the variation in the content of UV-absorbing phenolic compounds in leaves of temperate tree species along a light gradient? J. Exp. Bot. 2007, 58, 1753–1760. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  61. Goulas, Y.; Cerovic, Z.G.; Cartelat, A.; Moya, I. Dualex: A new instrument for field measurements of epidermal ultraviolet absorbance by chlorophyll fluorescence. Appl. Opt. 2004, 43, 4488–4496. [Google Scholar] [CrossRef]
  62. Lichtenthaler, H.K.; Schweiger, J. Cell wall bound ferulic acid, the major substance of the blue-green fluorescence emission of plants. J. Plant Physiol. 1998, 152, 272–282. [Google Scholar] [CrossRef]
  63. Buschmann, C.; Lichtenthaler, H.K. Principles and characteristics of multi-colour fluorescence imaging of plants. J. Plant Physiol. 1998, 152, 297–314. [Google Scholar] [CrossRef]
  64. Bürling, K.; Cerovic, Z.G.; Cornic, G.; Ducruet, J.M.; Noga, G.; Hunsche, M. Fluorescence-based sensing of drought-induced stress in the vegetative phase of four contrasting wheat genotypes. Environ. Exp. Bot. 2013, 89, 51–59. [Google Scholar] [CrossRef]
  65. Agati, G.; Pinelli, P.; Cortés, E.S.; Romani, A.; Cartelat, A.; Cerovic, Z.G. Nondestructive evaluation of anthocyanins in olive (Olea europaea) fruits by in situ chlorophyll fluorescence spectroscopy. J. Agric. Food Chem. 2005, 53, 1354–1363. [Google Scholar] [CrossRef] [PubMed]
  66. Chen, Y.; Chen, J.; Yan, B.; Cui, L.; Pan, J.; Kai, G. Agronomic traits of a new characteristic rice line ‘Huxuan 102’. Acta Agric. Shanghai 2015, 2, 60–64. [Google Scholar]
  67. Zubillaga, M.; Urricariet, S. Assessment of nitrogen status in wheat using aerial photography. Commun. Soil Sci. Plan. 2005, 36, 1787–1798. [Google Scholar] [CrossRef]
  68. Lemaire, G.; Jeuffroy, M.H.; Gastal, F. Diagnosis tool for plant and crop N status in vegetative stage: Theory and practices for crop N management. Eur. J. Agron. 2008, 28, 614–624. [Google Scholar] [CrossRef]
  69. Hussain, F.; Bronson, K.F.; Peng, S. Use of chlorophyll meter sufficiency indices for nitrogen management of irrigated rice in Asia. Agron. J. 2000, 92, 875–879. [Google Scholar]
  70. Varvel, G.E.; Wilhelm, W.W.; Shanahan, J.F.; Schepers, J.S. An algorithm for corn nitrogen recommendations using a chlorophyll meter based sufficiency index. Agron. J. 2007, 99, 701–706. [Google Scholar] [CrossRef]
  71. Cerovic, Z.G.; Ghozlen, N.B.; Milhade, C.; Obert, M.; Debuisson, S.; Le Moigne, M. Nondestructive diagnostic test for nitrogen nutrition of grapevine (Vitis vinifera L.) based on Dualex leaf-clip measurements in the field. J. Agr. Food Chem. 2015, 63, 3669–3680. [Google Scholar] [CrossRef] [PubMed]
  72. Stroppiana, D.; Boschetti, M.; Brivio, P.A.; Bocchi, S. Plant nitrogen concentration in paddy rice from field canopy hyperspectral radiometry. Field Crop. Res. 2009, 111, 119–129. [Google Scholar] [CrossRef]
  73. Tremblay, N.; Wang, Z.; Cerovic, Z.G. Sensing crop nitrogen status with fluorescence indicators. a review. Agron. Sustain. Dev. 2012, 32, 451–464. [Google Scholar] [CrossRef]
  74. Demotes-Mainard, S.; Boumaza, R.; Meyer, S.; Cerovic, Z.G. Indicators of nitrogen status for ornamental woody plants based on optical measurements of leaf epidermal polyphenol and chlorophyll contents. Sci. Hortic. 2008, 115, 377–385. [Google Scholar] [CrossRef]
  75. Agati, G.; Foschi, L.; Grossi, N.; Guglielminetti, L.; Cerovic, Z.G.; Volterrani, M. Fluorescence-based versus reflectance proximal sensing of nitrogen content in Paspalum vaginatum and Zoysia matrella turfgrasses. Eur. J. Agron. 2013, 45, 39–51. [Google Scholar] [CrossRef]
  76. Meyer, S.; Cerovic, Z.G.; Goulas, Y.; Montpied, P.; Demotesmainard, S.; Bidel, L.P.; Moya, I.; Dreyer, E. Relationships between optically assessed polyphenols and chlorophyll contents, and leaf mass per area ratio in woody plants: A signature of the carbon-nitrogen balance within leaves? Plant Cell Environ. 2006, 29, 1338–1348. [Google Scholar] [CrossRef] [PubMed]
Figure 1. Box plots of selected Multiplex index values for the above canopy, on-the-go, and leaf scale measuring modes at the panicle initiation and stem elongation stages in 2013: Simple fluorescence ratio under green excitation (SFR_G) (a), simple fluorescence ratio under red excitation (SFR_R) (b), flavonol (FLAV) (c), blue–green to far-red fluorescence ratio (BRR_FRF) (d), anthocyanins (ANTH) (e), nitrogen balance index under UV and red excitation (NBI_R) (f), and nitrogen balance index under UV and green excitation (NBI_G) (g). Within the same growth stage, different lowercase (panicle initiation stage) or uppercase letters (stem elongation stage) above or below the boxes indicate that the Multiplex index values differed significantly according to the least significant difference test at P ≤ 0.05.
Figure 1. Box plots of selected Multiplex index values for the above canopy, on-the-go, and leaf scale measuring modes at the panicle initiation and stem elongation stages in 2013: Simple fluorescence ratio under green excitation (SFR_G) (a), simple fluorescence ratio under red excitation (SFR_R) (b), flavonol (FLAV) (c), blue–green to far-red fluorescence ratio (BRR_FRF) (d), anthocyanins (ANTH) (e), nitrogen balance index under UV and red excitation (NBI_R) (f), and nitrogen balance index under UV and green excitation (NBI_G) (g). Within the same growth stage, different lowercase (panicle initiation stage) or uppercase letters (stem elongation stage) above or below the boxes indicate that the Multiplex index values differed significantly according to the least significant difference test at P ≤ 0.05.
Remotesensing 11 01847 g001
Figure 2. Mean value comparisons for each Multiplex index at the three growth stages and each N application rate (kg N ha−1) for variety KY 131 in 2013. Means and standard errors are shown in each cell for N rate treatments with each measurement mode (Above Canopy, AC; On-the-go, OG; Leaf Scale, LS). Different lowercase letters at the bottom of the plot at each growth stage indicate significant differences according to the least significant difference test at P ≤ 0.05.
Figure 2. Mean value comparisons for each Multiplex index at the three growth stages and each N application rate (kg N ha−1) for variety KY 131 in 2013. Means and standard errors are shown in each cell for N rate treatments with each measurement mode (Above Canopy, AC; On-the-go, OG; Leaf Scale, LS). Different lowercase letters at the bottom of the plot at each growth stage indicate significant differences according to the least significant difference test at P ≤ 0.05.
Remotesensing 11 01847 g002aRemotesensing 11 01847 g002b
Figure 3. The relative error (RE) values of the validation analysis based on the regression models of the six Multiplex indices and the N status indicators for (a) above-ground biomass (AGB), (b) plant N uptake (PNU), (c) leaf N concentration (LNC), (d) plant N concentration (PNC), and (e) nitrogen nutrition index (NNI) at the panicle initiation (PI), stem elongation (SE), and heading (HE) stages.
Figure 3. The relative error (RE) values of the validation analysis based on the regression models of the six Multiplex indices and the N status indicators for (a) above-ground biomass (AGB), (b) plant N uptake (PNU), (c) leaf N concentration (LNC), (d) plant N concentration (PNC), and (e) nitrogen nutrition index (NNI) at the panicle initiation (PI), stem elongation (SE), and heading (HE) stages.
Remotesensing 11 01847 g003
Table 1. Details of nitrogen (N) rate experiments with two rice cultivars conducted during 2012–2013.
Table 1. Details of nitrogen (N) rate experiments with two rice cultivars conducted during 2012–2013.
ExperimentYearCultivarTransplanting DateSampling Date and Stage
12012KY 13118 May21 June (PI), 29 June (SE), 23 July (HE)
12012LJ 2118 May25 June (PI), 2 July (SE), 23 July (HE)
22013KY 13117 May23 June (PI), 2 July (SE), 22 July (HE)
22013LJ 2117 May28 June (PI), 6 July (SE), 27 July (HE)
PI: Panicle Initiation; SE: Stem Elongation; HE: Heading.
Table 2. Description of the variables and indices recorded by the Multiplex sensor (modified from Table 1 by Zhang et al. [41]).
Table 2. Description of the variables and indices recorded by the Multiplex sensor (modified from Table 1 by Zhang et al. [41]).
VariablesFormulaExplanation
BGF_UV/Blue–green Fluorescence under UV excitation
RF_UV/Red Fluorescence under UV excitation
FRF_UV/Far-Red Fluorescence under UV excitation
BGF_G/Reflected Blue–Green light under Green excitation
RF_G/Red Fluorescence under Green excitation
FRF_G/Far-Red Fluorescence under Green excitation
RF_R/Red Fluorescence under Red excitation
FRF_R/Far-Red Fluorescence under Red excitation
SFR_GFRF_G/RF_GSimple Fluorescence Ratio under Green excitation
SFR_RFRF_R/RF_RSimple Fluorescence Ratio under Red excitation
BRR_FRFBGF_UV/FRF_UVBlue–green to Far-Red Fluorescence Ratio under UV excitation
FER_RUVFRF_R/FRF_UVFlavonols under Red and UV excitation
FLAVLog (FRF_R/FRF_UV)Flavonols under Red and UV excitation
FER_RGFRF_R/FRF_GAnthocyanins under Red and Green excitation
ANTHLog (FRF_R/FRF_G)Anthocyanins under Red and Green excitation
NBI_GFRF_UV/RF_GNitrogen Balance Index under UV and Green excitation
NBI_RFRF_UV/RF_RNitrogen Balance Index under UV and Red excitation
FERARI#Log (5000/FRF_R)Fluorescence Excitation Ratio Anthocyanin Relative Index
# the variable is not measured in the “on-the-go” mode.
Table 3. Significance test (ANOVA) of Multiplex variables measured in three modes across 5 N rates: Above canopy (AC), “on-the-go” (OG), and leaf scale (LS) at the panicle initiation (PI), stem elongation (SE), and heading (HE) growth stages for the rice varieties Kongyu 131 (KY 131) and Longjing 21 (LJ 21) in 2013.
Table 3. Significance test (ANOVA) of Multiplex variables measured in three modes across 5 N rates: Above canopy (AC), “on-the-go” (OG), and leaf scale (LS) at the panicle initiation (PI), stem elongation (SE), and heading (HE) growth stages for the rice varieties Kongyu 131 (KY 131) and Longjing 21 (LJ 21) in 2013.
VarietyStageBGF_U VRF_UVFRF_UVBGF_GRF_GFRF_GRF_RFRF_RSFR_GSFR_RBRR_FRFFLAVANTHNBI_GNBI_RFERARI
KY 131AC measurement mode
PINS**NS********NSNSNSNSNS**
SENSNSNSNSNSNSNSNSNS*NSNSNS*NSNS
OG measurement mode
PI***********************************
SE********NS***********************
HENS******NS***********************
LS measurement mode
PINSNSNSNSNSNSNS***NSNSNSNSNSNS
SENSNSNSNSNS*NS******NSNSNS****
HENSNSNS**NS*NS*****NSNS***
LJ 21AC measurement mode
PINS**NS*****************
SENS**NS********NSNSNSNS*
OG measurement mode
PINS**NSNS*NS************
SENS****NSNS**NS***************
HENS******NS*****NS***************NS******
LS measurement mode
PINSNSNSNSNSNSNSNS***NSNSNS*NSNS
SENS****NSNS*NSNS*NS**NSNS***
HENSNSNSNSNSNSNSNSNSNSNSNSNSNSNSNS
*** Significant at the 0.001 level; ** Significant at the 0.01 level; * Significant at the 0.05 level; NS Not significant.
Table 4. The coefficients of determination (R2) for the linear relationships between standard and normalized Multiplex indices and N status indicators (leaf N concentration (LNC), plant N concentration (PNC), N nutrition index (NNI), aboveground biomass (AGB) and plant N uptake (PNU)) for two varieties at the panicle initiation (PI), stem elongation (SE), and heading (HE) growth stages.
Table 4. The coefficients of determination (R2) for the linear relationships between standard and normalized Multiplex indices and N status indicators (leaf N concentration (LNC), plant N concentration (PNC), N nutrition index (NNI), aboveground biomass (AGB) and plant N uptake (PNU)) for two varieties at the panicle initiation (PI), stem elongation (SE), and heading (HE) growth stages.
Multiplex IndicesLNC (g kg−1)PNC (g kg−1)NNIAGB (t ha−1)PNU (kg ha−1)
PISEHEPISEHEPISEHEPISEHEPISEHE
Standard indices
SFR_G0.63 **0.30 **0.49 **0.64 **0.34 **0.46 **0.72 **0.31 **0.59 **0.60 **0.14 *0.41 **0.66 **0.21 **0.58 **
SFR_R0.59 **0.28 **0.42 **0.58 **0.34 **0.38 **0.66 **0.29 **0.54 **0.56 **0.13 *0.45 **0.61 **0.19 **0.57 **
BRR_FRF0.53 **0.52 **0.67 **0.47 **0.48 **0.66 **0.57 **0.48 **0.72 **0.50 **0.26 **0.33 **0.54 **0.39 **0.59 **
FLAV0.40 **0.64 **0.55 **0.39 **0.64 **0.55 **0.58 **0.73 **0.67 **0.55 **0.50 **0.38 **0.59 **0.68 **0.59 **
ANTH0.38 **0.12 *0.27 **0.41 **0.14 *0.33 **0.60 **0.10 *0.47 **0.60 **0.03NS0.36 **0.61 **0.06NS0.48 **
NBI_G0.54 **0.68 **0.62 **0.52 **0.71 **0.61 **0.69 **0.78 **0.76 **0.63 **0.50 **0.47 **0.68 **0.71 **0.71 **
NBI_R0.52 **0.67 **0.58 **0.52 **0.71 **0.56 **0.70 **0.77 **0.74 **0.64 **0.47 **0.51 **0.69 **0.68 **0.72 **
Normalized indices
SFR_GNSI0.58 **0.39 **0.67 **0.65 **0.42 **0.70 **0.55 **0.54 **0.69 **0.35 **0.45 **0.24 **0.43 **0.50 **0.52 **
SFR_RNSI0.57 **0.42 **0.62 **0.61 **0.46 **0.67 **0.52 **0.57 **0.68 **0.33 **0.45 **0.25 **0.40 **0.52 **0.52 **
BRR_FRFNSI0.49 **0.34 **0.63 **0.48 **0.41 **0.74 **0.41 **0.56 **0.76 **0.26 **0.50 **0.28 **0.33 **0.56 **0.58 **
FLAVNSI0.42 **0.51 **0.70 **0.44 **0.60 **0.76 **0.41 **0.74 **0.82 **0.26 **0.55 **0.34 **0.33 **0.70 **0.64 **
ANTHNSI0.51 **0.40 **0.57 **0.64 **0.40 **0.65 **0.54 **0.57 **0.56 **0.34 **0.52 **0.11 *0.42 **0.58 **0.35 **
NBI_GNSI0.59 **0.53 **0.69 **0.61 **0.63 **0.75 **0.55 **0.76 **0.78 **0.35 **0.55 **0.31 **0.43 **0.71 **0.61 **
NBI_RNSI0.60 **0.55 **0.69 **0.65 **0.64 **0.75 **0.58 **0.77 **0.79 **0.36 **0.56 **0.31 **0.46 **0.72 **0.61 **
** Correlation is significant at the 0.01 level; * Correlation is significant at the 0.05 level; NS Correlation is not significant.
Table 5. Equations and coefficients of determination of linear regression models (n = 40) at different growth stages based on the best performing Multiplex index and crop N indicators (LNC, PNC, NNI, PNU, and AGB).
Table 5. Equations and coefficients of determination of linear regression models (n = 40) at different growth stages based on the best performing Multiplex index and crop N indicators (LNC, PNC, NNI, PNU, and AGB).
Growth StageStandard IndicesModelR2Normalized IndicesModelR2
PISFR_GLNC = 4.468x + 5.9320.63NBI_RNSILNC = 23.918x + 10.4130.60
PISFR_GPNC = 2.912x + 4.9610.64NBI_RNSIPNC = 15.323x + 8.2470.65
PISFR_GNNI = 0.2442x-0.51880.72NBI_RNSINNI = 1.1412x − 0.11160.58
PINBI_RPNU = 88.184x-33.560.69NBI_RNSIPNU = 85.908x − 43.670.46
PINBI_GAGB = 1.5268x-1.15650.64NBI_RNSIAGB = 2.905x − 1.11840.36
SENBI_GLNC = 8.707x + 14.3520.68NBI_RNSILNC = 17.96x + 16.2790.55
SENBI_GPNC = 5.544x + 9.0820.71NBI_RNSIPNC = 12.317x + 9.5420.64
SENBI_GNNI = 0.5003x + 0.05820.78NBI_RNSINNI = 1.1571x + 0.06010.77
SENBI_GPNU = 51.494x-40.8730.71NBI_RNSIPNU = 120.8x − 42.1570.72
SENBI_GAGB = 1.7391x-0.45020.50NBI_RNSIAGB = 4.1975x − 0.59610.56
HEBRR_FRFLNC = -210.31x + 47.4520.67NBI_GNSILNC = 21.646x + 15.4730.69
HEBRR_FRFPNC = -131.79x + 24.3130.66FLAVNSIPNC = -31.59x + 49.5910.76
HENBI_GNNI = 0.5942x-0.0540.76FLAVNSINNI = -3.0631x + 4.39560.82
HENBI_RPNU = 213.07x-67.6230.72FLAVNSIPNU = -462.81x + 612.190.64
HENBI_RAGB = 8.3363x + 0.06090.51FLAVNSIAGB = -15.729x + 24.1310.34
Table 6. Agreement and Kappa statistics for different indices (SFR_G, SFR_R, BRR_FRF, FLAV, ANTH, NBI_G, and NBI_R) and corresponding normalized indices (SFR_GNSI, SFR_RNSI, BRR_FRFNSI, FLAVNSI, ANTHNSI, NBI_GNSI, and NBI_RNSI) regarding diagnostic results (Nitrogen Nutrition Index) at different growth stages.
Table 6. Agreement and Kappa statistics for different indices (SFR_G, SFR_R, BRR_FRF, FLAV, ANTH, NBI_G, and NBI_R) and corresponding normalized indices (SFR_GNSI, SFR_RNSI, BRR_FRFNSI, FLAVNSI, ANTHNSI, NBI_GNSI, and NBI_RNSI) regarding diagnostic results (Nitrogen Nutrition Index) at different growth stages.
ComparisonAgreement (%)Kappa statisticsComparisonAgreement (%)Kappa statistics
SEHESEHESEHESEHE
SFR_G and NNI75650.554 ***0.310 *SFR_GNSI and NNI70850.494 **0.661 ***
SFR_R and NNI70700.510 ***0.322 *SFR_RNSI and NNI75850.583 ***0.661 ***
BRR_FRF and NNI60700.363 *0.409 **BRR_FRFNSI and NNI80800.655 ***0.538 ***
FLAV and NNI75900.605 ***0.763 ***FLAVNSI and NNI75800.558 ***0.570 ***
ANTH and NNI55750.283NS0.355 *ANTHNSI and NNI80650.669 ***0.227NS
NBI_G and NNI75800.595 ***0.590 ***NBI_GNSI and NNI75850.673 ***0.698 ***
NBI_R and NNI75800.595 ***0.590 ***NBI_RNSI and NNI90850.840 ***0.698 ***
*** Significant at the 0.001 level; ** Significant at the 0.01 level; * Significant at the 0.05 level; NS Not significant.

Share and Cite

MDPI and ACS Style

Huang, S.; Miao, Y.; Yuan, F.; Cao, Q.; Ye, H.; Lenz-Wiedemann, V.I.S.; Bareth, G. In-Season Diagnosis of Rice Nitrogen Status Using Proximal Fluorescence Canopy Sensor at Different Growth Stages. Remote Sens. 2019, 11, 1847. https://doi.org/10.3390/rs11161847

AMA Style

Huang S, Miao Y, Yuan F, Cao Q, Ye H, Lenz-Wiedemann VIS, Bareth G. In-Season Diagnosis of Rice Nitrogen Status Using Proximal Fluorescence Canopy Sensor at Different Growth Stages. Remote Sensing. 2019; 11(16):1847. https://doi.org/10.3390/rs11161847

Chicago/Turabian Style

Huang, Shanyu, Yuxin Miao, Fei Yuan, Qiang Cao, Huichun Ye, Victoria I.S. Lenz-Wiedemann, and Georg Bareth. 2019. "In-Season Diagnosis of Rice Nitrogen Status Using Proximal Fluorescence Canopy Sensor at Different Growth Stages" Remote Sensing 11, no. 16: 1847. https://doi.org/10.3390/rs11161847

APA Style

Huang, S., Miao, Y., Yuan, F., Cao, Q., Ye, H., Lenz-Wiedemann, V. I. S., & Bareth, G. (2019). In-Season Diagnosis of Rice Nitrogen Status Using Proximal Fluorescence Canopy Sensor at Different Growth Stages. Remote Sensing, 11(16), 1847. https://doi.org/10.3390/rs11161847

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop