Next Article in Journal
Satellite Image Time Series Clustering via Time Adaptive Optimal Transport
Next Article in Special Issue
Predicting Equivalent Water Thickness in Wheat Using UAV Mounted Multispectral Sensor through Deep Learning Techniques
Previous Article in Journal
Integrated Preprocessing of Multitemporal Very-High-Resolution Satellite Images via Conjugate Points-Based Pseudo-Invariant Feature Extraction
Previous Article in Special Issue
Ir-UNet: Irregular Segmentation U-Shape Network for Wheat Yellow Rust Detection by UAV Multispectral Imagery
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Incorporating Multi-Scale, Spectrally Detected Nitrogen Concentrations into Assessing Nitrogen Use Efficiency for Winter Wheat Breeding Populations

by
Raquel Peron-Danaher
1,2,3,
Blake Russell
4,
Lorenzo Cotrozzi
1,5,†,
Mohsen Mohammadi
4 and
John J. Couture
1,3,5,*
1
Department of Entomology, Purdue University, West Lafayette, IN 47906, USA
2
Interdisciplinary Life Science Education Program, Purdue University, West Lafayette, IN 47906, USA
3
Center for Plant Biology, Purdue University, West Lafayette, IN 47906, USA
4
Department of Agronomy, Purdue University, West Lafayette, IN 47906, USA
5
Department of Forestry and Natural Resources, Purdue University, West Lafayette, IN 47906, USA
*
Author to whom correspondence should be addressed.
Current Address: Department of Agriculture, Forestry, and the Environment, University of Pisa, 56124 Pisa, Italy.
Remote Sens. 2021, 13(19), 3991; https://0-doi-org.brum.beds.ac.uk/10.3390/rs13193991
Submission received: 1 September 2021 / Revised: 24 September 2021 / Accepted: 26 September 2021 / Published: 6 October 2021
(This article belongs to the Special Issue Proximal and Remote Sensing for Precision Crop Management)

Abstract

:
Annually, over 100 million tons of nitrogen fertilizer are applied in wheat fields to ensure maximum productivity. This amount is often more than needed for optimal yield and can potentially have negative economic and environmental consequences. Monitoring crop nitrogen levels can inform managers of input requirements and potentially avoid excessive fertilization. Standard methods assessing plant nitrogen content, however, are time-consuming, destructive, and expensive. Therefore, the development of approaches estimating leaf nitrogen content in vivo and in situ could benefit fertilization management programs as well as breeding programs for nitrogen use efficiency (NUE). This study examined the ability of hyperspectral data to estimate leaf nitrogen concentrations and nitrogen uptake efficiency (NUpE) at the leaf and canopy levels in multiple winter wheat lines across two seasons. We collected spectral profiles of wheat foliage and canopies using full-range (350–2500 nm) spectroradiometers in combination with leaf tissue collection for standard analytical determination of nitrogen. We then applied partial least-squares regression, using spectral and reference nitrogen measurements, to build predictive models of leaf and canopy nitrogen concentrations. External validation of data from a multi-year model demonstrated effective nitrogen estimation at leaf and canopy level (R2 = 0.72, 0.67; root-mean-square error (RMSE) = 0.42, 0.46; normalized RMSE = 12, 13; bias = −0.06, 0.04, respectively). While NUpE was not directly well predicted using spectral data, NUpE values calculated from predicted leaf and canopy nitrogen levels were well correlated with NUpE determined using traditional methods, suggesting the potential of the approach in possibly replacing standard determination of plant nitrogen in assessing NUE. The results of our research reinforce the ability of hyperspectral data for the retrieval of nitrogen status and expand the utility of hyperspectral data in winter wheat lines to the application of nitrogen management practices and breeding programs.

1. Introduction

Wheat is one of the most widely cultivated cereal crops in the world. Annually, 220 million ha of wheat are planted, producing over 670 million tons of grain, which is mostly used as a food source for humans and animals [1]. Because of the widespread consumption, this cereal plays an essential role in global food security, representing ~20% of total caloric intake and protein ingestion for many [1]. Recent studies have shown a continued increase in wheat consumption by humans, which established a new record in 2018–2019 of approximately 602 million metric tons (MMT) of wheat consumed worldwide [2]. Additionally, projections estimate that the demand for wheat will increase by nearly 60% within the next 30 years due to projected human population growth [3]. This predicted increase in consumption poses a significant challenge for wheat breeders and producers, especially considering the estimated decrease in available areas for cultivation and an increase of extreme weather events [1].
The USA predominantly produces winter wheat (70–80% of total wheat planted) due to the cultivars’ high yield potential [4]. Winter wheat is planted in late summer and early fall, then germinates and establishes roots, and has a vernalization period in the winter, where a temperature below 8 °C is necessary to induce flowering [5]. Wheat grain is composed of approximately 12–15% of protein [6], and consequently wheat needs a high rate of nitrogen uptake that in turn requires an elevated rate of nitrogen fertilization for optimal productivity [7]. The high fertilization rate is due to the low efficiency of nitrogen uptake and remobilization to grain during filling [8,9]. Consequently, a large amount of nitrogen is applied in fields annually, over 100 million tons globally, to ensure optimal wheat productivity [10], and wheat fertilization represents nearly 20% of the total fertilizer input used in crop plants globally [11]. The high amount of fertilizer input required for optimal wheat production, however, can harm the environment and represents a substantial production cost for farmers [1,9,10,12,13]. Developing approaches to quantify wheat nitrogen status in the field can aid in determining optimal fertilizer inputs which will help mitigate the adverse effects of under- and over-fertilization.
In addition to improving the management of field-based fertilization, selecting lines with enhanced nitrogen use efficiency (NUE) can help mitigate the high demand by wheat for nitrogen fertilizer inputs. NUE is defined as the ability of a plant to increase grain dry-weight yield per unit of nitrogen uptake from the soil [14,15] and is composed of two key elements: nitrogen uptake efficiency (NUpE) and nitrogen utilization efficiency (NUtE) [16]. Monitoring nitrogen fertilizer input, as well as quantifying NUE, requires accurate nitrogen determination in plant tissue, especially at the scale for in-field monitoring and assessing agronomic breeding populations. Traditional methods are labor-intensive and time-consuming, requiring tissue collection and chemical analysis. These analyses also provide temporally limited information, and because of the costs, measurements are collected infrequently, limiting information about nutrient dynamics throughout a growing season. Thus, there is a critical need to develop more rapid and efficient methods for assessing NUE, in vivo and in situ, to assist breeding programs [17].
Recent advances in sensor-based approaches for monitoring plant status have revealed promising outcomes [18], and hyperspectral reflectance has gained considerable attention by demonstrating the ability to estimate plant chemical composition. The ability of hyperspectral data to confirm plant chemical composition relies on variation in the absorption of light due to differential patterns of vibrational excitation of molecular bonds, primarily C-H, N-H, and O-H bonds, at specific wavelengths in the visible (VIS; 400–700 nm), near-infrared (NIR; 700–1100 nm), and short-wave infrared (SWIR; 1100–2500 nm) ranges. Thus, variations in plant spectral profiles are determined by the chemical and physical composition of the vegetation [19,20], and the use of spectral data from vegetation, combined with mathematical methods, has been demonstrated to be an efficient approach to accurately predict plant traits, such as water status, physiological processes, and primary and secondary metabolites [21,22,23,24,25,26,27,28,29,30,31,32,33].
While previous studies showed that hyperspectral data applied to nitrogen assessment in wheat is an efficient method for estimating nitrogen concentrations in wheat [22,23,29,32,34,35,36], those studies were mostly conducted using a single measurement approach and often did not explore either the scalability of hyperspectral data over measurement scales (e.g., leaf vs. canopy) or the transferability of models across seasons [31]. Furthermore, the ability of spectral data to scale across different genotypes within species and years is not well understood and needs further study to identify the range of applications for the use of hyperspectral data in agronomic practices [30,31,37].
To address these knowledge gaps, we examined the ability of hyperspectral data to identify traits related to nitrogen dynamics in different wheat lines, across two growing seasons, and assessed the ability of hyperspectral data in estimating NUpE in a wheat breeding population. Specifically, we (1) developed predictive models using hyperspectral data to estimate nitrogen concentration at leaf and canopy levels, (2) assessed model performance to estimate nitrogen concentration across different wheat lines and growing seasons, (3) estimated NUE directly from spectral data, and (4) used predicted nitrogen values to calculate NUE components and test the accuracy of this method with conventional NUE calculations.

2. Materials and Methods

2.1. Experimental Design

This project was conducted at the Agronomy Center for Research and Education at Purdue University (ACRE), West Lafayette, IN, USA (40°47′10”N, −87°01′04”W, annual precipitation: 1047 nn/year, average annual temperature: 11.5 °C). The soil type at ACRE is described as a combination of: Rockfield silt loam, Fincastle silt loam, and Toronto silt loam (USDA Web Soil Survey). Additional information about detailed soil characteristics can be found in Appendix A, Table A1. The 2017 experiment was planted in September 2016 and included four different wheat lines, showing a range of yield performance under nitrogen fertilization regimes, selected from the Purdue soft-red winter wheat breeding program, listed in Appendix B, Table A2. The 2018 experiment was planted in October 2017 and included 32 different lines, also selected from the Purdue soft-red winter wheat breeding program, showing a range of yield performance under nitrogen fertilization regimes, and listed in Appendix B, Table A2. The goal of this paper was to determine the ability of hyperspectral data to estimate foliar nitrogen levels and to assess the utility of the approach in breeding wheat lines for NUE and not to assess responses of individual wheat lines to different fertilizer levels. Specific information about the responses of different wheat lines to variable fertilization and the resulting NUE can be found in [38].
In both years, all wheat lines were planted in eight replicate plots for each genotype: four plots with baseline fertilization and four plots with baseline and supplemental fertilization. Plot dimensions were 1.22 m × 3.05 m and contained seven rows with 15 cm intervals; seeds were planted at a density of 370 seeds m−2. In both years, prior to planting, all plots were fertilized with monoammonium phosphate ((NH4)H2PO4 (N, P, K: 11-52-0)) at a rate of 24.6 kg ha−1. In the spring following planting, four replicates of each line received supplemental fertilizer treatments with urea (N, P, K: 46-0-0) treated with Limus (BASF, Germany) during the stem elongation growth stage (Zadoks 30) at a rate of 112.09 kg ha−1 establishing a low, consisting of lines that were fertilized only prior to planting, and a high, consisting of lines fertilized both prior and after planting, fertilization treatments. Detailed soil composition is described in Appendix A, Table A1.

2.2. Spectral and Reference Tissue Collections

2.2.1. Data Collection 2017

Hyperspectral data were collected from leaf and canopy levels using full-range spectroradiometers (HR-1024i, Spectra Vista Corporation, Poughkeepsie, NY, USA). Collections were carried out on all four wheat lines at six different time points (TPs) over one growing season. Each TP represented a different physiological plant stage: booting (TP1-26/04/2017; Zadoks 45), onset of heading (TP2-08/05/2017; Zadoks 50), heading (TP3-15/05/2017; Zadoks 58), flowering (TP4-22/05/2017; Zadoks 60), post-anthesis (TP5-06/06/2017; Zadoks 68), and maturity (TP6-12/06/2017; Zadoks 91). Each plot was subdivided into three subplots, and spectral and foliar collections were carried out in one randomly selected subplot within each plot. A total of 288 samples, with paired spectral and foliar measurements, were collected and distributed as follows: 32 samples from TP1 and TP2 (4 lines × 2 fertilization treatments × 4 plots = 32) and 64 samples from TP3, TP4, TP5, and TP6 each (4 lines × 2 fertilization treatments × 8 plots × 4 time points = 256). During the first two time points (TP1 and TP2), weather and logistical constraints limited us to sampling only four of the eight total plots. Canopy reflectance measurements were conducted between 11:00 and 13:00 h., under clear sky conditions, using an eight-degree foreoptic (representing approximately a 12 cm−2 area) mounted to a tripod at one meter above the canopy, with a foreoptic angle approximately at nadir of the target. Following the canopy measurements, four flag leaves were collected and stored in separate paper bags on ice. After returning to the laboratory (<2 h), leaf reflectance was measured on the adaxial leaf surface, using the same spectroradiometer but fitted with a fiber optic cable attached to a plant probe with a leaf clip, containing an internal halogen light source fitted in a lens tube adapter for small spot (ideal for viewing samples with narrow area, such as wheat leaves). Relative reflectance for both leaf and canopy measurements was determined by dividing vegetation radiance by the radiance of a 99% reflective white panel (Spectralon, Labsphere, North Sutton, NH, USA) that was collected either just above the canopy, for canopy measurements, or mounted within the internal leaf clip, for leaf measurements, measured every 15 or 12 spectral collections for leaf and canopy measurements, respectively. After measurements, leaf samples were flash-frozen in liquid nitrogen and stored in a freezer at −20 °C until they were oven-dried at 90 °C for seven days. After reaching a constant dry mass, samples were ground for chemical analysis using a ball mill.

2.2.2. Data Collection 2018

Hyperspectral data were collected from leaf and canopy levels from all 32 wheat lines at a single time point: heading and the beginning of flowering (from 24/05/2018 to 06/06/2018; Zadoks 58–60). Each plot was subdivided into two subplots, and spectral and foliar collections were carried out in one randomly selected subplot within each plot for spectral and leaf collection. The reason for this reduction in subplot sampling was to reduce sampling time so samples were of comparable growth stages. A total of 128 samples with paired spectral and foliar measurements were collected (32 lines × 2 treatment × 2 plots = 128). The sample selection for nitrogen analysis was carried out assigning random samples from subplots collected across all genotypes and both treatments, and the spectral data collection and tissue sampling and processing for leaf and canopy measurements followed the same protocol as described in the 2017 data collection.

2.3. Standard Nitrogen Analysis

Analysis of nitrogen concentrations from foliar tissue collected from both years was performed using a Thermo Finnigan Flash 1112 elemental analyzer (San Jose, CA, USA). Atropine (CE Elantech, Lakewood, NJ, USA) was used as a standard [39,40]. The total leaf nitrogen concentration, in percentage dry mass, for the 2017 and 2018 datasets, ranged from 0.64 to 4.85 and 1.24 to 4.88, respectively.

2.4. Chemometric Modeling

We generated chemometric models to predict nitrogen concentration in wheat from spectral data using partial least-squares regression (PLSR) [41,42]. In cases where predictor variables are highly correlated and outnumber response variables, such as often the case with hyperspectral data, classical regression techniques can produce unreliable beta coefficient estimates [42]. In contrast to standard regression approaches, PLSR reduces many inter-correlated predictor variables into relatively few, uncorrelated latent variables and is a preferred method in chemometrics [29,43,44,45,46,47,48]. The number of latent variables included in each model was determined by the reduction of the predicted residual error sum of squares (PRESS) statistic [49] following leave-one-out (LOO) cross-validation. The latent variables were then combined into a linear model to predict nitrogen concentrations.
To avoid the influence of canopy geometry and solar intensity during canopy measurements and to optimize PLSR models, reflectance data were vector-normalized [50]. The first derivative of spectral data was then calculated on the normalized reflectance data to highlight spectral differences for both leaf and canopy measurements. The wavelength range 400–2400 nm was used in both leaf and canopy models to capture absorption features related to both protein and chlorophyll content [19]. When using canopy spectral data, the wavelength region between 1800 and 2000 nm was omitted to remove the noise caused by atmospheric water absorption. Preliminary models identified poorly predicted outliers, detected by elevated reflectance in the visible wavelength region, spectral jumps that occur when the leaf clip is not fully engaged, or inconsistencies or errors in the reference data [29,48,50,51,52,53,54,55,56,57,58]. Outliers removed accounted for approximately 10 and 15% of the dataset for leaf and canopy collections, respectively.
For the 2017 data, we evaluated the model performance by conducting 500 randomized permutations of the entire 2017 dataset using 70% of the data, divided into 80:20 split for calibration and cross-validation, respectively, with the remaining 30% held out for external validation. For each permutation, we tracked the model goodness of fit (R2), overall error rate (root-mean-square error, RMSE), and normalized root-mean-square error (NRMSE) over the data range together with a bias to assess model performance when applied to the validation dataset. These randomized analyses generated a distribution of fit statistics allowing for the assessment of model stability as well as uncertainty in model predictions. We also determined the contribution of individual wavelength for the PLSR using the variable importance for the projection (VIP) selection statistic and standardized (i.e., centered and scaled) coefficients selection statistics, which explain the contribution of a single wavelength for the variance in the latent variables [44,58]. The modeling approach was performed using the ‘pls’ package in R (www.r-project.org accessed on 2 October 2019).
To assess cross-year transferability of spectral models, we applied models that were developed to predict nitrogen concentrations in 2017 to the spectral data collected in 2018. To determine if we could improve the model fit by adding data collection in 2018, we included 50% of the reference data (spectral data that had a paired reference nitrogen measurement) collected in 2018 (64 samples) and created a cross-year model, following the same modeling approach described above, and used the remaining 50% of the reference data (64 samples) collected in 2018 as external validation.

Categorization of Standardized Coefficients and VIP within Regions of Chlorophyll and Protein Absorption Features

To align spectral features with chlorophyll and protein absorption features, we categorized wavelength coefficients and VIP statistics (Appendix D, Table A4 and Table A5) having a higher contribution to model performance into general regions of literature-published chlorophyll and protein absorption features. First, we averaged standardized coefficients and VIP values every 30 nm starting at 400 nm to account for multicollinearity. Then, we ranked the top 30 averaged standardized coefficient values. For the VIP statistic, we excluded all values that were below 0.8, as suggested by [41], then averaged the VIP values as we did for coefficients, and then ranked the top 30 VIP values. Next, we selected literature-published chlorophyll and protein absorption features (Appendix C, Table A3) and binned wavelengths 40 nm on center of the published wavelength peak (i.e., 20 nm before and 20 nm after the described peak). This process was carried out to capture minor absorption features associated with chlorophyll and protein outside of the singular wavelength absorption feature. Finally, we identified averaged wavelength coefficients and VIP values that fell within the binned values for each chlorophyll and protein absorption feature range (Appendix D, Table A6, Table A7 and Table A8).

2.5. NUE Quantification

We used two approaches to investigate the potential of spectral data to benefit NUE estimation. First, we tested the ability of spectral data to directly predict NUpE, following a similar chemometric modeling protocol described above, with the exception that we used only LOO cross-validation. We used this approach because the logistical constraints of collecting the biomass and chemical data to calculate NUpE, a component of NUE, limited the number of samples we were able to collect, and the total number of samples used in this analysis was considerably smaller (n = 39) than the number of samples collected for previous chemometric modeling of nitrogen concentrations. Second, we calculated NUpE using spectrally predicted nitrogen to test the accuracy of predicting component processes of NUE using spectral data. Reference NUpE was calculated using the nitrogen straw values, and estimated NUpE was calculated by substituting mid-season nitrogen values predicted by spectral models for straw values. NUpE was calculated as the following:
N U p E = G r a i n   N + S t r a w   N N   S u p p l i e d
where G r a i n   N is the nitrogen concentration in the grain, S t r a w   N is the whole-plant (leaves, stems, tillers, organs) nitrogen level, and N   S u p p l i e d is the quantity of nitrogen applied via fertilizer (kg N ha−1) [38].

3. Results

3.1. Chemometric Models

3.1.1. Modeling and Cross-Year Validation 2017

The 2017 model accurately predicted nitrogen at both leaf and canopy levels, although there was a decline in external validation performance (Table 1, Figure 1), and these declines were comparable.
Standardized coefficients and VIP profiles for nitrogen prediction were more prominent in regions of known chlorophyll and protein absorption features. Overall standardized coefficient and VIP peaks aligned with wavelengths at approximately 430, 460, 550, 700, 910, 1000, 1400, 1500, 1700, 1940, 2050 and 2300 nm (Figure 2, Appendix D, Table A4 and Table A5). The top 30 standardized coefficients in the leaf-level model were found more in the SWIR region (1100-2500 nm) than the canopy models, which leveraged more coefficients in the VNIR region (Appendix D, Table A6).
The 2017 nitrogen prediction models applied to the 2018 dataset revealed a skewed nitrogen prediction (Figure 3), with a decrease in R2 (0.71 to 0.69 and 0.73 to 0.45 for leaf and canopy models, respectively), an overall change in RMSE (0.45 to 0.41 and 0.42 to 0.52 for leaf and canopy models, respectively) and NRMSE (14 to 11 and 13 to 12 for leaf and canopy models, respectively), and an increase in bias (0.009 to −0.34 and 0.003 to −1.17 for leaf and canopy estimates, respectively).

3.1.2. Multi-Year Models: Combining 2017 and 2018 Datasets

The multi-year model showed improved results compared with the cross-year predictions. The multi-year model yielded an accurate predictive model for nitrogen for leaf and canopy levels (Table 2) with good external validation (Figure 4).
Similar to the 2017 models, the multi-year models presented peaks for standardized coefficients and VIP profiles for nitrogen prediction in spectral regions with known chlorophyll and protein absorption features (Figure 5, Appendix D Table A7 and Table A8). Multi-year coefficients and VIP values were highly correlated with 2017 outputs for leaf (Pearson’s correlation r = 0.64, p < 0.001, and 0.94, p <0.001, for coefficients and VIP values, respectively, n = 2001) and canopy (Pearson’s correlation r = 0.68, p < 0.001, and 0.98, p < 0.001, for coefficients and VIP values, respectively, n = 1801) models. Standardized coefficients and VIP statistics generally aligned with wavelengths at approximately 430, 500, 630 700, 790, 900, 1000, 1300, 1400, 1900, 2000, 2100, 2300, 2350 and 2390 nm (Figure 5 and Appendix D Table A7 and Table A8). Again, similar to the 2017 models, the multi-year leaf-level models leveraged relatively more coefficients in the SWIR region compared to the canopy model, which had more coefficients in the VNIR region (Appendix D Table A8).

3.2. NUE Prediction

Models directly predicting NUpE were less accurate than the models predicting nitrogen. NUpE models presented an average R2 of 0.52 and 0.37, an average RMSE of 0.30 and 0.28, an average NRMSE of 8 and 9, and an average bias of 0.01 and 0.01 for leaf and canopy levels, respectively (Figure 6). NUpE calculations were strongly correlated at both leaf and canopy levels (Pearson’s correlation r = 0.86, p <0 .001, n = 18, and 0.76, p = 0.004, n = 12, respectively) when we compared NUpE using reference nitrogen or nitrogen levels predicted by the multi-year model (Figure 7). We also found that using spectrally predicted nitrogen for NUpE calculations resulted in an overestimation of NUpE values. Normalizing NUpE for both estimation and reference data, however, decreased bias by 95 and 93% and revealed a higher correlation with r = 0.93, p = 0.001, and bias = 0.11 for leaf level and r = 0.87, p = 0.0002, and bias = 0.17 for canopy level, respectively.

4. Discussion

Wheat production has doubled in the last four decades, and the need to increase production of this globally important crop, while minimizing fertilization inputs, is critical. In this paper, we show that (1) hyperspectral data, collected at multiple different scales, can be used for estimation of foliar nitrogen levels in wheat, (2) updating of models to include a greater genetic variability and nitrogen range improves prediction accuracy, and (3) spectral estimates of foliar nitrogen, at different scales, can potentially be used as a replacement for reference data to identify NUE components in a breeding population.
The use of hyperspectral data to estimate plant traits, including nitrogen, is becoming more common with the expansion of precision and digital agriculture, and recent studies have demonstrated positive results for monitoring nitrogen status in several crops [27,31,34,39,53,54,55,56,57,58,59,60]. While several of these studies explored the ability of hyperspectral data to access nitrogen content in crops, including winter wheat, most studies have focused on the VIS and NIR wavelength ranges, relying primarily on relationships with chlorophyll content to interpret nitrogen content [60]. It is possible that using multiple wheat stages over the season in the models from 2017 influenced the estimation outcomes because of changes in leaf color. However, the linkage between chlorophyll and nitrogen can vary among different ecosystems, environmental conditions, and biotic and abiotic stressors [59,60]. Methods that rely solely on chlorophyll content for nitrogen assessment can be misleading, especially in the context of nitrogen fertilization management, once several other factors lead to chlorophyll reduction in leaves, such as sulfur deficiency, vegetative stage of the plants, and disease [60], and thus measurements are likely more sensitive to changes in leaf color (i.e., green vs. senescent vegetation). Chlorophyll accounts for a smaller portion of nitrogen contained within leaves compared with protein (e.g., Rubisco), which can account for over 50% of nitrogen content in plant tissue [60,61,62,63,64]. Rubisco also is a principal source for nitrogen remobilization as well and is a key enzyme driving photosynthesis [61]. Full-spectrum data can leverage both chlorophyll and protein absorption features that are closely related to leaf nitrogen content. Our findings suggest that spectral data can accurately predict nitrogen at both leaf- and canopy-level measurements, corroborating with the findings of other studies [27,29,31,60,61,64].
Our assessment of a single-season model revealed a high prediction accuracy for nitrogen in the 2017 leaf-level model. This response was likely bolstered by data collected across the growing season in 2017 that captured a wide range of nitrogen concentrations in leaves, especially at lower levels of nitrogen in mature leaves due to nitrogen remobilization [64]. The estimation of nitrogen from samples collected in 2018 using the model built solely from the model developed in 2017, however, revealed a decline in model performance when compared with the 2017-only model. We anticipated this decline because of the introduction of increasing genetic variability and a greater range of nitrogen and because we included a wider range of lines from the NUE breeding population used in 2017. A decline in model performance when transferring models across years to novel genetic material has been previously demonstrated [65]. We found that when the model was updated to capture the new cultivars and increased range of nitrogen levels (i.e., 2017 and 2018 multi-year model) performance returned to near-original performance. These outcomes highlight the importance of updating predictive spectral models as novel genetic material and values outside the original model range become available and increase transferability across years and genotypes but highlight that developmental stage and environment also likely influence these results.
We found that the wavelengths having significant contributions to models estimating nitrogen were related to chlorophyll and protein absorption features, which we expected given the contribution of these compounds to nitrogen levels in vegetative tissue [59]. We also found the wavelength region near 1400 nm was an important feature for predicting nitrogen at both leaf and canopy levels. In fresh leaves, this feature is found in a wavelength range strongly influenced by water absorption yet has previously been shown to contribute to spectral models estimating nitrogen in fresh and dry samples [51,66,67,68,69,70]. This outcome, coupled with those of other studies, suggests that, even while contained in a region strongly affected by water absorption, this feature plays an important role in nitrogen retrievals using spectral data.
Spectral data have been directly, through statistical approaches, or indirectly, via indices, used to estimate plant physiology, chemical composition, and water status [21,24,28,29,32,33,37,54,58]. Spectral data have also been used to assess within-plant nitrogen dynamics, including NUE and the sub-component NUpE [71]. NUpE is affected by several factors, including the vegetative stage of the plant, soil nitrogen content, and above-ground biomass [39]. While we found that spectral data did not directly estimate NUpE well, we could use spectrally predicted nitrogen as a reasonable surrogate for reference nitrogen when calculating NUpE. Although we found that bias in the relationship between NUpE calculated using either reference or spectrally predicted nitrogen was high, this can possibly be accounted for because the reference nitrogen for the traditionally calculated NUpE values was end-of-season nitrogen and had much lower nitrogen levels than the mid-season spectrally predicted nitrogen, which was much higher, explaining the overestimation of predicted NUpE values. One potential issue with modeling NUpE directly using mid-season foliar spectral data is the potential temporal disconnection with the nitrogen levels estimated during the mid-season and the yield data needed to calculate NUpE. One logistical drawback we encountered when trying to model NUpE using spectral data was a small sample size, which potentially did not provide a sufficient range of NUpE to characterize a relationship between spectral data and NUpE. The lack of a direct spectral relationship with NUpE contrasts with recent studies relating spectral data with NUE parameters [72,73]. Previous work has shown a significant correlation of canopy reflectance and the Maccioni index in the grain filling phase but highlighted the importance of large-plot research to study the relationship of vegetation indexes and nitrogen use dynamics [71]. It has also been suggested that heading is the optimal time to predict NUE using canopy reflectance, but the correlation between reflectance and NUE during the heading phase declined over a second season [63]. Regardless of the lack of a direct relationship between spectral data and NUpE, the ability of spectrally estimated mid-season nitrogen to calculate NUpE shows that spectral data have considerable potential in monitoring nitrogen dynamics in wheat, especially in a breeding context.

5. Conclusions

Our results demonstrate that hyperspectral data, combined with chemometric modeling, is an effective method to monitor nitrogen dynamics in a diverse breeding population of wheat. We found reasonable transferability of models across years using full-range spectral data and improved model performance as we included greater genetic, developmental, and environmental variation of spectral and reference data. We also found promising model performance for canopy-based models, suggesting the potential for scaling measurements from leaves to fields using unmanned or manned aerial vehicles and highlighting the ability for whole-field monitoring of nitrogen management as a potential replacement for classical phenotyping approaches. Our research builds on the importance of including genotypic, phenotypic, and temporal variation in developing hyperspectral models estimating plant functional traits [30,31]. Moreover, we feel an important outcome of this research is the incorporation of hyperspectral data into crop nitrogen management and plant breeding for enhanced NUE. Outcomes of this work highlight the potential of hyperspectral data for precision nitrogen assessment in crops in both management and breeding contexts.

Author Contributions

Conceptualization, R.P.-D., L.C. and J.J.C.; methodology, R.P.-D., B.R., L.C., M.M. and J.J.C.; validation, R.P.-D. and J.J.C.; formal analysis, R.P.-D. and J.J.C.; investigation, R.P.-D. and J.J.C.; resources, J.J.C. and M.M.; data curation, R.P.-D.; writing—original draft preparation, R.P.-D.; writing—review and editing, R.P.-D., B.R., L.C., M.M. and J.J.C.; supervision, J.J.C.; project administration, J.J.C.; funding acquisition, J.J.C. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by USDA NIFA Hatch awards IND011490 to J.J.C. and 1013073 to M.M.

Acknowledgments

We would like to thank Allie Duffy, Aimee Jasmine Stramowski, Anupreet Saini, Hanna Quellhorst, Kirsten Brichler, Marguerite Bolt, and Taylor Nelson for assistance with field collections.

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, or in the decision to publish the results.

Appendix A

Table A1. Detailed soil characteristics of the experimental plots.
Table A1. Detailed soil characteristics of the experimental plots.
ComponentsPlot 5 NorthPlot 5 South
Phosphorus (Bray P1lbs/acre)4353
Potassium (pounds/acre)195188
Calcium (pounds/acre)23262214
Magnesium (pounds/acre)715621
Base Sat’n (%)97100
Calcium Sat’n (%)6266
Ca/Mg Ratio22.2
CEC (meq/100 g)11.099.95
Potassium Sat’n (%)2.522.71
Magnesium Sat’n (%)3130
Mg/K Ratio12.511.2
pH6.937
Lime Index69.7370

Appendix B

Table A2. List of wheat winter lines used in both years of experiments with respective yield indication. Wheat lines: 2017—2017 dataset; 2018—2018 dataset. Specific information on cultivar yield responses to variable fertilizer levels can be found in [38]. Copyright © 2020 Russell, Guzman and Mohammadi. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
Table A2. List of wheat winter lines used in both years of experiments with respective yield indication. Wheat lines: 2017—2017 dataset; 2018—2018 dataset. Specific information on cultivar yield responses to variable fertilizer levels can be found in [38]. Copyright © 2020 Russell, Guzman and Mohammadi. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
DatasetWheat Line ID (Purdue ID)
201704606RA1-7-1-4 (PU01)
053A1-2-5-3-5 (PU02)
10101RA1-6-2 (PU26)
11407A1-6 (PU24)
20180175A1-31-4-1 (PU17)
03549A1-18-25-4 (PU05)
04606RA1-7-1 (PU07)
04606RA1-7-1-4 (PU01)
04606RA1-7-1-6 (PU23)
04719A1-16-1-1-47-4 (PU19)
05247A1-7-3-29 (PU10)
05247A1-7-3-98 (PU08)
05247A-7-7-3-1 (PU25)
05251A1-1-77-16-3 (PU22)
0537A1-3-12-1 (PU14)
053A1-2-5-3-5-3 (PU02)
0566A1-3-1-48 (PU09)
0570A1-2-39-2-4 (PU06)
0570A1-8-5-1 (PU03)
057RA1-8-5-3 (PU13)
057RA1-8-5-33 (PU11)
06497A1-7-3 (PU21)
07117B1-29-7-9-9-4-3-6-3 (PU29)
0722A1-1-1-7 (PU04)
07469A1-6-1-1 (PU18)
0762A1-2-8 (PU15)
08334A1-31 (PU12)
10101RA1-6-2 (PU26)
10221A1-8-1 (PU27)
10222A1-9-2 (PU30)
1041RB1-10 (PU16)
10447A1-5 (PU28)
10565C1-1 (PU20)
11407A1-6 (PU24)
CHECK—P25R40
CHECK—P25R62

Appendix C

Table A3. Specific absorption bands related to plant chlorophyll a and b combined, chlorophyll a, or chlorophyll b and protein and the respective absorption mechanisms described in literature (Curran et al., 1989; Gitelson et al., 1995; Datt et al., 1998; Curran et al., 2001; Kumar et al., 2005; Berger et al., 2020) [19,59,73,74,75,76].
Table A3. Specific absorption bands related to plant chlorophyll a and b combined, chlorophyll a, or chlorophyll b and protein and the respective absorption mechanisms described in literature (Curran et al., 1989; Gitelson et al., 1995; Datt et al., 1998; Curran et al., 2001; Kumar et al., 2005; Berger et al., 2020) [19,59,73,74,75,76].
ChemicalWL Highest PeakAbsorption MechanismsLiterature Described the WL
Chlorophyll a420Electron transitionKumar et al., 2005
Chlorophyll a430Electron transitionCurran et al., 1989
Chlorophyll b435Electron transitionKumar et al., 2005
Chlorophyll b460Electron transitionCurran et al., 1989
Chlorophyll a490Electron transitionKumar et al., 2005
Chlorophyll530Electron transitionCurran et al., 2001
Chlorophyll a 550Electron transitionDatt et al., 1998 and Gitelson et al., 1995
Chlorophyll630Electron transitionCurran et al., 2001
Chlorophyll b640Electron transitionCurran et al., 1989
Chlorophyll b643Electron transitionKumar et al.,
Chlorophyll a660Electron transitionCurran et al., 1989; Kumar et al., 2005
Chlorophyll675Electron transitionDatt et al., 1998
Chlorophyll700Not describedCurran et al., 2001; Gitelson et al., 1995
Protein910C-H stretch, 3rd overtoneCurran et al., 1989
Protein1020N-H stretchCurran et al., 1989
Protein1500Not describedKumar et al., 2005
Protein1510N-H stretch, 1st overtoneCurran et al., 1989
Protein1520N-H stretch, 1st overtoneBerger et al., 2020
Protein 1680C-H strecth, 1st overtoneKumar et al., 2005
Protein1690C-H strecth, 1st overtoneBerger et al., 2020
Protein1730C-H stretch Kumar et al., 2005
Protein 1940O-H strech, O-H deformationCurran et al., 1989; Kumar et al., 2005
Protein1960N-H assymmetryBerger et al., 2020
Protein1980N-H assymmetryCurran et al., 1989
Protein2050N-H stretch, N=H rotationKumar et al., 2005
Protein2060N-H stretch, N=H rotationCurran et al., 1989
Protein2130N-H stretchCurran et al., 1989
Protein 2170Not describedKumar et al., 2005
Protein2180N-H rotation, C-H stretch, C-O stretch, C=O stretchCurran et al., 1989
Protein2200N-H rotation, C-H stretch, C-O stretch, C=O stretchBerger et al., 2020
Protein2240C-H stretchCurran et al., 1989
Protein2270C-H stretchBerger et al., 2020
Protein2290C-H rotation, C=O stretch, N-H stretchKumar et al., 2005
Protein2300C-H rotation, C=O stretch, N-H stretchCurran et al., 1989
Protein2350CH2 rotation, C-H deformationCurran et al., 1989

Appendix D

Table A4. Standardized coefficient values of wavelengths associated with chlorophyll a and b combined, chlorophyll a, or chlorophyll b, and protein absorption feature for 2017 and multi-year models at both leaf and canopy levels. Bold values, standardized coefficient value within the top 30 highest standardized coefficients. NA, wavelength removed from the canopy data.
Table A4. Standardized coefficient values of wavelengths associated with chlorophyll a and b combined, chlorophyll a, or chlorophyll b, and protein absorption feature for 2017 and multi-year models at both leaf and canopy levels. Bold values, standardized coefficient value within the top 30 highest standardized coefficients. NA, wavelength removed from the canopy data.
2017 ModelsMulti-Year Models
WL Highest PeakChemicalLeaf LevelCanopy LevelLeaf LevelCanopy Level
420Chlorophyll a0.420.090.690.09
430Chlorophyll a0.290.050.590.08
435Chlorophyll b0.260.050.570.09
460Chlorophyll b0.20.050.330.09
490Chlorophyll a0.190.060.210.10
530Chlorophyll0.370.080.580.13
550Chlorophyll a 0.310.120.380.20
630Chlorophyll0.050.060.100.13
640Chlorophyll b0.060.050.130.11
643Chlorophyll b0.060.040.130.10
660Chlorophyll a0.060.030.150.07
675Chlorophyll0.140.140.310.21
700Chlorophyll0.380.250.580.32
910Protein0.10.050.150.11
1020Protein0.050.040.130.06
1500Protein 0.050.010.050.02
1510Protein 0.030.010.050.02
1520Protein0.030.010.040.01
1680Protein 0.020.010.050.01
1690Protein0.020.010.050.01
1730Protein 0.040.010.060.01
1940Protein0.05NA0.08NA
1960Protein0.06NA0.10NA
1980Protein0.07NA0.11NA
2050Protein 0.160.030.170.05
2060Protein 0.170.030.200.05
2130Protein0.070.010.180.02
2170Protein 0.050.010.120.04
2180Protein0.030.010.100.03
2200Protein0.030.010.070.03
2240Protein0.080.030.140.05
2270Protein0.080.030.140.04
2290Protein 0.080.030.150.04
2300Protein 0.090.030.140.04
2350Protein 0.130.030.150.05
Table A5. Variable importance for the projection (VIP) values of wavelengths associated with chlorophyll a and b combined, chlorophyll a, or chlorophyll b, and protein absorption feature for 2017 and multi-year models at both leaf and canopy levels. Bold values, VIP ≥ 0.8. NA, wavelength removed from the canopy data.
Table A5. Variable importance for the projection (VIP) values of wavelengths associated with chlorophyll a and b combined, chlorophyll a, or chlorophyll b, and protein absorption feature for 2017 and multi-year models at both leaf and canopy levels. Bold values, VIP ≥ 0.8. NA, wavelength removed from the canopy data.
2017 ModelsMulti-Year Models
Wavelength (nm)ChemicallLeaf LevelCanopy LevelLeaf LevelCanopy Level
420Chlorophyll a0.81.60.81.6
430Chlorophyll a0.91.60.91.6
435Chlorophyll b1.01.71.01.7
460Chlorophyll b1.11.71.11.7
490Chlorophyll a1.51.71.51.7
530Chlorophyll1.91.41.91.4
550Chlorophyll a 1.61.71.61.7
630Chlorophyll1.11.41.41.4
640Chlorophyll b1.11.21.21.2
643Chlorophyll b1.21.21.21.2
660Chlorophyll a1.31.01.30.9
675Chlorophyll1.51.21.51.2
700Chlorophyll1.61.41.61.4
910Protein0.70.60.70.7
1020Protein0.80.40.90.4
1500Protein1.00.61.00.7
1510Protein1.00.61.00.7
1520Protein1.00.61.00.6
1680Protein0.91.10.91.1
1690Protein0.91.20.91.1
1730Protein1.00.91.00.9
1940Protein0.9NA0.9NA
1960Protein1.1NA1.0NA
1980Protein1.1NA1.0NA
2050Protein0.90.70.90.7
2060Protein0.90.71.00.7
2130Protein0.80.30.80.4
2170Protein0.60.20.70.3
2180Protein0.60.20.60.3
2200Protein0.50.20.60.3
2240Protein0.90.70.90.8
2270Protein0.70.40.70.4
2290Protein0.60.30.60.3
2300Protein0.70.30.60.3
2350Protein0.70.30.60.3
Table A6. Top 30 standardized coefficient values and their respective wavelength ranges from 2017 leaf and canopy models. The coefficients values displayed in this table represent the average standardized coefficient value every 30nm, as described by the wavelength range (WL Range (nm)). Colors highlighting the three main areas of the full-spectral range: green, visible (VIS); orange, near-infrared (NIR); and red, short-wave infrared (SWIR).
Table A6. Top 30 standardized coefficient values and their respective wavelength ranges from 2017 leaf and canopy models. The coefficients values displayed in this table represent the average standardized coefficient value every 30nm, as described by the wavelength range (WL Range (nm)). Colors highlighting the three main areas of the full-spectral range: green, visible (VIS); orange, near-infrared (NIR); and red, short-wave infrared (SWIR).
2017 Models
Leaf LevelCanopy Level
WL Range (nm)CoefficientWL Range (nm)Coefficient
400-4290.51400-4290.10
430-4590.21430-4590.06
460-4890.18460-4890.06
490-5190.23490-5190.07
520-5490.40520-5490.08
550-5790.21550-5790.16
670-6990.28580-6090.12
700-7290.33610-6390.08
760-7890.12670-6990.22
790-8190.12700-7290.34
910-9390.13730-7590.57
940-9690.14760-7890.07
970-9990.41790-8190.11
1120-11490.09820-8490.09
1300-13290.11850-8790.08
1330-13590.15880-9090.05
1360-13890.16910-9390.08
1390-14190.13940-9690.08
1420-14490.13970-9990.06
1450-14790.111000-10290.04
1870-18990.281030-10590.05
1900-19290.111060-10890.04
1990-20190.111120-11490.04
2020-20490.141150-11790.04
2050-20790.181240-12690.04
2230-22590.091330-13590.05
2290-23190.091360-13890.08
2320-23490.121390-14190.14
2350-23790.111420-14490.07
2380-24000.132370-24000.05
Table A7. Variable importance for the projection (VIP) values ≥ 0.8 2017 and multi-year models at both leaf and canopy levels. The VIP values displayed in this table represent the average VIP value after averaging VIP data every 30nm, as described by the wavelength range (WL Range (nm)). Colors highlighting the three main areas of the full-spectral range: green, visible (VIS); orange, near-infrared (NIR); and red, short-wave infrared (SWIR).
Table A7. Variable importance for the projection (VIP) values ≥ 0.8 2017 and multi-year models at both leaf and canopy levels. The VIP values displayed in this table represent the average VIP value after averaging VIP data every 30nm, as described by the wavelength range (WL Range (nm)). Colors highlighting the three main areas of the full-spectral range: green, visible (VIS); orange, near-infrared (NIR); and red, short-wave infrared (SWIR).
2017 ModelMulti-Year Models
Leaf LevelCanopy LevelLeaf LevelCanopy Level
WL Range (nm)VIP ≥ 0.8WL Range (nm)VIP ≥ 0.8WL Range (nm)VIP ≥ 0.8WL Range (nm)VIP ≥ 0.8
430-4591.1400-4291.6430-4591.1400-4291.5
460-4891.2430-4591.7460-4891.2430-4591.7
490-5191.9460-4891.7490-5191.9460-4891.7
520-5491.9490-5191.6520-5491.9490-5191.6
550-5791.4520-5491.4550-5791.3520-5491.4
580-6091.2550-5791.7580-6091.2550-5791.7
610-6391.1580-6091.6610-6391.1580-6091.7
640-6691.3610-6391.5640-6691.3610-6391.5
670-6991.6640-6690.9670-6991.6640-6690.8
700-7291.7670-6991.4700-7291.7670-6991.4
730-7591.8700-7291.4730-7592.0700-7291.4
760-7891.1730-7591.7760-7891.1730-7591.7
790-8191.1760-7891.0790-8191.1760-7891.0
820-8491.0790-8191.4820-8491.0790-8191.4
850-8790.9820-8491.5850-8790.9820-8491.5
1000-10290.8850-8791.41000-10290.9850-8791.3
1030-10590.8880-9090.81030-10590.8880-9090.8
1120-11490.91060-10891.41090-11190.81060-10891.4
1270-12990.81210-12391.41120-11490.91210-12391.4
1300-13290.91240-12691.31210-12390.81240-12691.2
1330-13590.91270-12991.11270-12990.81270-12991.2
1360-13890.91420-14490.81300-13290.91420-14490.8
1390-14191.11600-16291.21330-13590.91600-16291.2
1420-14491.01660-16891.21360-13890.91660-16891.1
1450-14791.01690-17191.11390-14191.11690-17191.1
1480-15091.01720-17490.91420-14491.01720-17490.9
1510-15391.0 1450-14790.92040-20690.8
1540-15690.9 1480-15091.02220-22490.8
1570-15990.9 1510-15390.9
1600-16290.8 1540-15690.9
1660-16890.9 1570-15990.9
1690-17191.1 1600-16290.9
1720-17490.9 1660-16890.9
1750-17790.9 1690-17191.0
1780-18090.8 1720-17490.9
1810-18390.9 1750-17790.9
1840-18691.0 1780-18090.8
1870-18991.1 1810-18390.9
1900-19290.9 1840-18690.9
1930-19591.0 1870-18991.1
1960-19891.1 1900-19290.9
1990-20191.1 1930-19591.0
2020-20490.9 1960-19891.0
2050-20790.9 1990-20191.0
2080-21090.8 2020-20490.9
2110-21390.8 2050-20791.0
2230-22591.0 2080-21090.8
2110-21390.8
2140-21690.8
2230-22590.9
Table A8. Top 30 standardized coefficient values and their respective wavelength ranges from multi-year leaf and canopy models. The coefficient values displayed in this table represent the average standardized coefficient value every 30 nm, as described by the wavelength range (WL Range (nm)). Colors highlighting the three main areas of the full-spectral range: green, visible (VIS); orange, near-infrared (NIR); and red, short-wave infrared (SWIR).
Table A8. Top 30 standardized coefficient values and their respective wavelength ranges from multi-year leaf and canopy models. The coefficient values displayed in this table represent the average standardized coefficient value every 30 nm, as described by the wavelength range (WL Range (nm)). Colors highlighting the three main areas of the full-spectral range: green, visible (VIS); orange, near-infrared (NIR); and red, short-wave infrared (SWIR).
Multi-Year Models
Leaf LevelCanopy Level
WL Range (nm)CoefficientWL Range (nm)Coefficient
400-4290.77400-4290.10
430-4590.45430-4590.10
460-4890.23460-4890.09
490-5190.25490-5190.11
520-5490.64520-5490.14
550-5790.23550-5790.24
640-6690.16580-6090.20
670-6990.52610-6390.15
700-7290.40670-6990.31
730-7590.26700-7290.35
760-7890.15730-7590.54
790-8190.17760-7890.18
910-9390.18790-8190.11
940-9690.24820-8490.13
970-9990.51850-8790.08
1000-10290.15880-9090.07
1300-13290.16910-9390.17
1330-13590.18940-9690.18
1360-13890.30970-9990.12
1390-14190.331030-10590.07
1420-14490.271090-11190.08
1450-14790.151120-11490.07
1870-18990.151150-11790.07
1900-19290.281330-13590.20
2050-20790.221360-13890.14
2110-21390.171390-14190.21
2140-21690.191420-14490.08
2230-22590.161780-20090.08
2350-23790.162340-23690.07
2380-24000.192370-24000.09

References

  1. Shiferaw, B.; Smale, M.; Braun, H.-J.; Duveiller, E.; Reynolds, M.; Muricho, G. Crops that feed the world Past successes and future challenges to the role played by wheat in global food security. Food Secur. 2013, 5, 291–317. [Google Scholar] [CrossRef] [Green Version]
  2. FAO. GIEWS Crop Prospects and Food Situation. 2 July 2019. Available online: http://www.fao.org/documents/card/en/c/ca5327en (accessed on 15 January 2020).
  3. Phillips, B.S.; Norton, R. Global Wheat Production and Fertilizer Use. Better Crop. 2012, 96, 4–6. [Google Scholar]
  4. Rose, G.J. Crop Production. Soil Sci. 1956, 81, 152. [Google Scholar] [CrossRef]
  5. Li, G.; Yu, M.; Fang, T.; Cao, S.; Carver, B.F.; Yan, L. Vernalization requirement duration in winter wheat is controlled by T a VRN—A 1 at the protein level. Plant J. 2013, 76, 742–753. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  6. Dong, K.; Zhen, S.; Cheng, Z.; Cao, H.; Ge, P.; Yan, Y. Proteomic Analysis Reveals Key Proteins and Phosphoproteins upon Seed Germination of Wheat (Triticum aestivum L.). Front. Plant Sci. 2015, 6. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  7. Delogu, G.; Cattivelli, L.; Pecchioni, N.; De Falcis, D.; Maggiore, T.; Stanca, A. Uptake and agronomic efficiency of nitrogen in winter barley and winter wheat. Eur. J. Agron. 1998, 9, 11–20. [Google Scholar] [CrossRef]
  8. Barbottin, A.; LeComte, C.; Bouchard, C.; Jeuffroy, M.-H. Nitrogen Remobilization during Grain Filling in Wheat: Genotypic and Environmental Effects. Crop. Sci. 2005, 45, 1141–1150. [Google Scholar] [CrossRef]
  9. Robertson, G.; Vitousek, P.M. Nitrogen in Agriculture: Balancing the Cost of an Essential Resource. Annu. Rev. Environ. Resour. 2009, 34, 97–125. [Google Scholar] [CrossRef] [Green Version]
  10. Garnett, T.; Plett, D.; Heuer, S.; Okamoto, M. Genetic approaches to enhancing nitrogen-use efficiency (NUE) in cereals: Challenges and future directions. Funct. Plant Biol. 2015, 42, 921. [Google Scholar] [CrossRef]
  11. Heisey, P.W.; Norton, G. Fertilizers and other farm chemicals. In Handbook of Agricultural Economics; Evenson, R., Pingali, P., Eds.; Elsevier North Holland: Amsterdam, The Netherlands, 2007; Volume 3, pp. 2741–2777. [Google Scholar]
  12. Wei, F.; Yan, Z.; Yongchao, T.; Weixing, C.; Xia, Y.; Yingxue, L. Monitoring leaf nitrogen accumulation in wheat with hyper-spectral remote sensing. Acta Ecol. Sin. 2008, 28, 23–32. [Google Scholar] [CrossRef]
  13. Hitz, K.; Clark, A.J.; Van Sanford, D.A. Identifying nitrogen-use efficient soft red winter wheat lines in high and low nitrogen environments. Field Crop. Res. 2017, 200, 1–9. [Google Scholar] [CrossRef] [Green Version]
  14. Moll, R.H.; Kamprath, E.J.; Jackson, W.A. Analysis and interpretation of factors which contribute to efficiency of nitrogen-utilization. Agron. J. 1982, 74, 562–564. [Google Scholar] [CrossRef]
  15. Kong, L.; Xie, Y.; Hu, L.; Feng, B.; Li, S. Remobilization of vegetative nitrogen to developing grain in wheat (Triticum aestivum L.). Field Crop. Res. 2016, 196, 134–144. [Google Scholar] [CrossRef]
  16. Han, M.; Okamoto, M.; Beatty, P.H.; Rothstein, S.J.; Good, A.G. The Genetics of Nitrogen Use Efficiency in Crop Plants. Annu. Rev. Genet. 2015, 49, 269–289. [Google Scholar] [CrossRef]
  17. Araus, J.L.; Kefauver, S.C.; Zaman-Allah, M.; Olsen, M.S.; Cairns, J. Translating High-Throughput Phenotyping into Genetic Gain. Trends Plant Sci. 2018, 23, 451–466. [Google Scholar] [CrossRef] [Green Version]
  18. Roitsch, T.; Cabrera-Bosquet, L.; Fournier, A.; Ghamkhar, K.; Jiménez-Berni, J.; Pinto, F.; Ober, E.S. Review: New sensors and data-driven approaches—A path to next generation phenomics. Plant Sci. 2019, 282, 2–10. [Google Scholar] [CrossRef]
  19. Curran, P.J. Remote sensing of foliar chemistry. Remote Sens. Environ. 1989, 30, 271–278. [Google Scholar] [CrossRef]
  20. Cotrozzi, L.; Townsend, P.A.; Pellegrini, E.; Nali, C.; Couture, J.J. Reflectance spectroscopy: A novel approach to better understand and monitor the impact of air pollution on Mediterranean plants. Environ. Sci. Pollut. Res. 2017, 25, 8249–8267. [Google Scholar] [CrossRef]
  21. Serbin, S.P.; Dillaway, D.N.; Kruger, E.L.; Townsend, P.A. Leaf optical properties reflect variation in photosynthetic metabolism and its sensitivity to temperature. J. Exp. Bot. 2011, 63, 489–502. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  22. Li, F.; Mistele, B.; Hu, Y.; Chen, X.; Schmidhalter, U. Reflectance estimation of canopy nitrogen content in winter wheat using optimised hyperspectral spectral indices and partial least squares regression. Eur. J. Agron. 2014, 52, 198–209. [Google Scholar] [CrossRef]
  23. Féret, J.-B.; Berger, K.; de Boissieu, F.; Malenovský, Z. PROSPECT-PRO for estimating content of nitrogen-containing leaf proteins and other carbon-based constituents. Remote Sens. Environ. 2020, 252, 112173. [Google Scholar] [CrossRef]
  24. Ollinger, S.V. Sources of variability in canopy reflectance and the convergent properties of plants. New Phytol. 2010, 189, 375–394. [Google Scholar] [CrossRef]
  25. Couture, J.J.; Serbin, S.P.; Townsend, P.A. Spectroscopic sensitivity of real-time, rapidly induced phytochemical change in response to damage. New Phytol. 2013, 198, 311–319. [Google Scholar] [CrossRef]
  26. Serbin, S.P.; Singh, A.; Desai, A.R.; Dubois, S.G.; Jablonski, A.D.; Kingdon, C.; Kruger, E.L.; Townsend, P. Remotely estimating photosynthetic capacity, and its response to temperature, in vegetation canopies using imaging spectroscopy. Remote Sens. Environ. 2015, 167, 78–87. [Google Scholar] [CrossRef]
  27. Yuan, M.; Couture, J.J.; Townsend, P.A.; Ruark, M.D.; Bland, W.L. Spectroscopic Determination of Leaf Nitrogen Concentration and Mass Per Area in Sweet Corn and Snap Bean. Agron. J. 2016, 108, 2519–2526. [Google Scholar] [CrossRef]
  28. Cotrozzi, L.; Couture, J.J.; Cavender-Bares, J.; Kingdon, C.; Fallon, B.; Pilz, G.; Pellegrini, E.; Nali, C.; A Townsend, P. Using foliar spectral properties to assess the effects of drought on plant water potential. Tree Physiol. 2017, 37, 1582–1591. [Google Scholar] [CrossRef] [PubMed]
  29. Silva-Perez, V.; Molero, G.; Serbin, S.P.; Condon, A.G.; Reynolds, M.P.; Furbank, R.T.; Evans, J.R. Hyperspectral reflectance as a tool to measure biochemical and physiological traits in wheat. J. Exp. Bot. 2017, 69, 483–496. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  30. Meacham-Hensold, K.; Montes, C.; Wu, J.; Guan, K.; Fu, P.; Ainsworth, E.A.; Pederson, T.; Moore, C.; Brown, K.L.; Raines, C.; et al. High-throughput field phenotyping using hyperspectral reflectance and partial least squares regression (PLSR) reveals genetic modifications to photosynthetic capacity. Remote Sens. Environ. 2019, 231, 111176. [Google Scholar] [CrossRef] [PubMed]
  31. Bruning, B.; Berger, B.; Lewis, M.; Liu, H.; Garnett, T. Approaches, applications, and future directions for hyperspectral vegetation studies: An emphasis on yield-limiting factors in wheat. Plant Phenome J. 2020, 3, 1–22. [Google Scholar] [CrossRef]
  32. Cotrozzi, L.; Peron, R.; Tuinstra, M.R.; Mickelbart, M.V.; Couture, J.J. Spectral Phenotyping of Physiological and Anatomical Leaf Traits Related with Maize Water Status. Plant Physiol. 2020, 184, 1363–1377. [Google Scholar] [CrossRef]
  33. Campos-Medina, V.A.; Cotrozzi, L.; Stuart, J.J.; Couture, J.J. Spectral characterization of wheat functional trait responses to Hessian fly: Mechanisms for trait-based resistance. PLoS ONE 2019, 14, e0219431. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  34. Tan, C.; Du, Y.; Zhou, J.; Wang, D.; Luo, M.; Zhang, Y.; Guo, W. Analysis of Different Hyperspectral Variables for Diagnosing Leaf Nitrogen Accumulation in Wheat. Front. Plant Sci. 2018, 9, 674. [Google Scholar] [CrossRef]
  35. Liang, L.; Di, L.; Huang, T.; Wang, J.; Lin, L.; Wang, L.; Yang, M. Estimation of Leaf Nitrogen Content in Wheat Using New Hyperspectral Indices and a Random Forest Regression Algorithm. Remote Sens. 2018, 10, 1940. [Google Scholar] [CrossRef] [Green Version]
  36. Duan, D.-D.; Zhao, C.-J.; Li, Z.-H.; Yang, G.-J.; Yang, W.-D. Estimating total leaf nitrogen concentration in winter wheat by canopy hyperspectral data and nitrogen vertical distribution. J. Integr. Agric. 2019, 18, 1562–1570. [Google Scholar] [CrossRef]
  37. Yendrek, C.R.; Tomaz, T.; Montes, C.; Cao, Y.; Morse, A.M.; Brown, P.J.; McIntyre, L.; Leakey, A.D.; Ainsworth, E.A. High-Throughput Phenotyping of Maize Leaf Physiological and Biochemical Traits Using Hyperspectral Reflectance. Plant Physiol. 2016, 173, 614–626. [Google Scholar] [CrossRef]
  38. Russell, B.; Guzman, C.; Mohammadi, M. Cultivar, Trait and Management System Selection to Improve Soft-Red Winter Wheat Productivity in the Eastern United States. Front. Plant Sci. 2020, 11. [Google Scholar] [CrossRef] [PubMed]
  39. Muñoz-Huerta, R.F.; Guevara-Gonzalez, R.G.; Contreras-Medina, L.M.; Torres-Pacheco, I.; Prado-Olivarez, J.; Ocampo-Velazquez, R.V. A Review of Methods for Sensing the Nitrogen Status in Plants: Advantages, Disadvantages and Recent Advances. Sensors 2013, 13, 10823–10843. [Google Scholar] [CrossRef] [PubMed]
  40. Dumas, J.B.A. Procedes de l’analyse Organic. Annales de Chimie et de Physique. Ann. Chem. Phys. 1831, 247, 198–213. [Google Scholar]
  41. Wold, S.; Ruhe, A.; Wold, H.; Dunn, W. The collinearity problem in linear regression. The partial least squares (PLS) approach to generalized inverses. SIAM J. Sci. Stat. Comput. 1984, 5, 735–743. [Google Scholar] [CrossRef] [Green Version]
  42. Wold, S.; Sjöström, M.; Eriksson, L. PLS-regression: A basic tool of chemometrics. Chemom. Intell. Lab. Syst. 2001, 58, 109–130. [Google Scholar] [CrossRef]
  43. Grossman, Y.; Ustin, S.; Jacquemoud, S.; Sanderson, E.; Schmuck, G.; Verdebout, J. Critique of stepwise multiple linear regression for the extraction of leaf biochemistry information from leaf reflectance data. Remote Sens. Environ. 1996, 56, 182–193. [Google Scholar] [CrossRef]
  44. Cotrozzi, L.; Lorenzini, G.; Nali, C.; Pellegrini, E.; Saponaro, V.; Hoshika, Y.; Arab, L.; Rennenberg, H.; Paoletti, E. Hyperspectral Reflectance of Light-Adapted Leaves Can Predict Both Dark- and Light-Adapted Chl Fluorescence Parameters, and the Effects of Chronic Ozone Exposure on Date Palm (Phoenix dactylifera). Int. J. Mol. Sci. 2020, 21, 6441. [Google Scholar] [CrossRef] [PubMed]
  45. Bolster, K.L.; Martin, M.E.; Aber, J.D. Determination of carbon fraction and nitrogen concentration in tree foliage by near infrared reflectances: A comparison of statistical methods. Can. J. For. Res. 1996, 26, 590–600. [Google Scholar] [CrossRef] [Green Version]
  46. Atzberger, C.; Guérif, M.; Baret, F.; Werner, W. Comparative analysis of three chemometric techniques for the spectroradiometric assessment of canopy chlorophyll content in winter wheat. Comput. Electron. Agric. 2010, 73, 165–173. [Google Scholar] [CrossRef]
  47. Serbin, S.P.; Singh, A.; McNeil, B.E.; Kingdon, C.; Townsend, P. Spectroscopic determination of leaf morphological and biochemical traits for northern temperate and boreal tree species. Ecol. Appl. 2014, 24, 1651–1669. [Google Scholar] [CrossRef] [Green Version]
  48. Couture, J.J.; Singh, A.; Rubert-Nason, K.F.; Serbin, S.; Lindroth, R.L.; Townsend, P. Spectroscopic determination of ecologically relevant plant secondary metabolites. Methods Ecol. Evol. 2016, 7, 1402–1412. [Google Scholar] [CrossRef] [Green Version]
  49. Chen, S.; Hong, X.; Harris, C.J.; Sharkey, P.M. Sparse Modeling Using Orthogonal Forward Regression With PRESS Statistic and Regularization. IEEE Trans. Syst. Man Cybern. Part B (Cybern.) 2004, 34, 898–911. [Google Scholar] [CrossRef]
  50. Wang, Z.; Chlus, A.; Geygan, R.; Ye, Z.; Zheng, T.; Singh, A.; Couture, J.J.; Cavender-Bares, J.; Kruger, E.L.; Townsend, P.A. Foliar functional traits from imaging spectroscopy across biomes in eastern North America. New Phytol. 2020, 228, 494–511. [Google Scholar] [CrossRef]
  51. Cotrozzi, L.; Couture, J.J. Hyperspectral assessment of plant responses to multi-stress environments: Prospects for managing protected agrosystems. Plants People Planet 2019, 2, 244–258. [Google Scholar] [CrossRef]
  52. Marchica, A.; Loré, S.; Cotrozzi, L.; Lorenzini, G.; Nali, C.; Pellegrini, E.; Remorini, D. Early Detection of Sage (Salvia officinalis L.) Responses to Ozone Using Reflectance Spectroscopy. Plants 2019, 8, 346. [Google Scholar] [CrossRef] [Green Version]
  53. Yu, K.-Q.; Zhao, Y.-R.; Li, X.-L.; Shao, Y.-N.; Liu, F.; He, Y. Hyperspectral Imaging for Mapping of Total Nitrogen Spatial Distribution in Pepper Plant. PLoS ONE 2014, 9, e116205. [Google Scholar] [CrossRef] [PubMed]
  54. Chong, I.-G.; Jun, C.-H. Performance of some variable selection methods when multicollinearity is present. Chemom. Intell. Lab. Syst. 2005, 78, 103–112. [Google Scholar] [CrossRef]
  55. Acevedo, M.; Zurn, J.D.; Molero, G.; Singh, P.; He, X.; Aoun, M.; McCandless, L. The role of wheat in global food security. In Agricultural Development and Sustainable Intensification: Technology and Policy Challenges in the Face of Climate Change; Routledge: Abingdon, UK, 2018; pp. 81–110. [Google Scholar] [CrossRef]
  56. Hawkesford, M.J. Reducing the reliance on nitrogen fertilizer for wheat production. J. Cereal Sci. 2013, 59, 276–283. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  57. Yousfi, S.; Marin Peira, J.F.; De la Horra, G.R.; Ablanque, P.V.M. Remote Sensing: Useful Approach for Crop Nitrogen Management and Sustainable Agriculture. In Soil Managment and Plant Nutrition for susteinable Crop Production; IntechOpen: London, UK, 2019; pp. 1–12. [Google Scholar] [CrossRef]
  58. Calzone, A.; Cotrozzi, L.; Remorini, D.; Lorenzini, G.; Nali, C.; Pellegrini, E. Oxidative stress assessment by a spectroscopic approach in pomegranate plants under a gradient of ozone concentrations. Environ. Exp. Bot. 2020, 182, 104309. [Google Scholar] [CrossRef]
  59. Berger, K.; Verrelst, J.; Féret, J.-B.; Wang, Z.; Wocher, M.; Strathmann, M.; Danner, M.; Mauser, W.; Hank, T. Crop nitrogen monitoring: Recent progress and principal developments in the context of imaging spectroscopy missions. Remote Sens. Environ. 2020, 242, 111758. [Google Scholar] [CrossRef]
  60. Li, H.; Zhang, Y.; Lei, Y.; Antoniuk, V.; Hu, C. Evaluating Different Non-Destructive Estimation Methods for Winter Wheat (Triticum aestivum L.) Nitrogen Status Based on Canopy Spectrum. Remote Sens. 2019, 12, 95. [Google Scholar] [CrossRef] [Green Version]
  61. Homolová, L.; Malenovský, Z.; Clevers, J.G.; García-Santos, G.; Schaepman, M.E. Review of optical-based remote sensing for plant trait mapping. Ecol. Complex. 2013, 15, 1–16. [Google Scholar] [CrossRef] [Green Version]
  62. Kokaly, R.F.; Asner, G.; Ollinger, S.; Martin, M.E.; Wessman, C.A. Characterizing canopy biochemistry from imaging spectroscopy and its application to ecosystem studies. Remote Sens. Environ. 2009, 113, S78–S91. [Google Scholar] [CrossRef]
  63. Chapin, F.S.; Bloom, A.J.; Field, C.B.; Waring, R.H. Plant Responses to Multiple Environmental Factors. BioScience 1987, 37, 49–57. [Google Scholar] [CrossRef]
  64. Ellis, R. The most abundant protein in the world. Trends Biochem. Sci. 1979, 4, 241–244. [Google Scholar] [CrossRef]
  65. Ecarnot, M.; Compan, F.; Roumet, P. Assessing leaf nitrogen content and leaf mass per unit area of wheat in the field throughout plant cycle with a portable spectrometer. Field Crop. Res. 2013, 140, 44–50. [Google Scholar] [CrossRef]
  66. Ge, Y.; Atefi, A.; Zhang, H.; Miao, C.; Ramamurthy, R.K.; Sigmon, B.; Yang, J.; Schnable, J.C. High-throughput analysis of leaf physiological and chemical traits with VIS–NIR–SWIR spectroscopy: A case study with a maize diversity panel. Plant Methods 2019, 15, 1–12. [Google Scholar] [CrossRef] [Green Version]
  67. Mahajan, G.R.; Sahoo, R.N.; Pandey, R.N.; Gupta, V.K.; Kumar, D. Using hyperspectral remote sensing techniques to monitor nitrogen, phosphorus, sulphur and potassium in wheat (Triticum aestivum L.). Precis. Agric. 2014, 15, 499–522. [Google Scholar] [CrossRef]
  68. Fourty, T.; Baret, F.; Jacquemoud, S.; Schmuck, G.; Verdebout, J. Leaf optical properties with explicit description of its biochemical composition: Direct and inverse problems. Remote Sens. Environ. 1996, 56, 104–117. [Google Scholar] [CrossRef]
  69. Wang, J.; Chen, J.M.; Ju, W.; Qiu, F.; Zhang, Q.; Fang, M.; Chen, F. Limited Effects of Water Absorption on Reducing the Accuracy of Leaf Nitrogen Estimation. Remote Sens. 2017, 9, 291. [Google Scholar] [CrossRef] [Green Version]
  70. Frels, K.; Guttieri, M.; Joyce, B.; Leavitt, B.; Baenziger, P.S. Evaluating canopy spectral reflectance vegetation indices to estimate nitrogen use traits in hard winter wheat. Field Crop. Res. 2018, 217, 82–92. [Google Scholar] [CrossRef]
  71. Pavuluri, K.; Chim, B.K.; Griffey, C.A.; Reiter, M.S.; Balota, M.; Thomason, W.E. Canopy spectral reflectance can predict grain nitrogen use efficiency in soft red winter wheat. Precis. Agric. 2014, 16, 405–424. [Google Scholar] [CrossRef]
  72. Nigon, T.; Yang, C.; Paiao, G.D.; Mulla, D.; Knight, J.; Fernández, F. Prediction of Early Season Nitrogen Uptake in Maize Using High-Resolution Aerial Hyperspectral Imagery. Remote Sens. 2020, 12, 1234. [Google Scholar] [CrossRef] [Green Version]
  73. Gitelson, A.A.; Merzlyak, M.N. Signature Analysis of Leaf Reflectance Spectra: Algorithm Development for Remote Sensing of Chlorophyll. J. Plant Physiol. 1996, 148, 494–500. [Google Scholar] [CrossRef]
  74. Datt, B. Remote Sensing of Chlorophyll a, Chlorophyll b, Chlorophyll a+b, and Total Carotenoid Content in Eucalyptus Leaves. Remote Sens. Environ. 1998, 66, 111–121. [Google Scholar] [CrossRef]
  75. Curran, P.J.; Dungan, J.L.; Peterson, D.L. Estimating the foliar biochemical concentration of leaves with reflectance spectrometry: Testing the Kokaly and Clark methodologies. Remote Sens. Environ. 2001, 76, 349–359. [Google Scholar] [CrossRef]
  76. Kumar, L.; Schmidt, K.; Dury, S.; Skidmore, A. Imaging Spectrometry and Vegetation Science. In Imaging Spectrom; Springer: Dordrecht, The Netherlands, 2006; Volume 4, pp. 111–155. [Google Scholar] [CrossRef]
Figure 1. Models for 2017: observed vs. predicted external validation values of nitrogen levels from the 2017 dataset. (a) External validation of nitrogen estimates using leaf-level spectral measurements. (b) External validation of nitrogen estimates using canopy-level spectral measurements. Error bars for estimated values represent the standard deviation generated from 500 model permutations using an 80:20 calibration:cross-validation approach. Dashed line: 1:1 relationship. R2: model goodness of fit; RMSE: root-mean-square error (% dry mass); NRMSE: normalized root-mean-square error (%). TP1, time point 1 = booting (Zadoks 45); TP2, time point 2 = onset of heading (Zadoks 50); TP3, time point 3 = heading (Zadoks 58); TP4, time point 4 = flowering (Zadoks 60); TP5, time point 5 = post-anthesis (Zadoks 68); TP6, time point 6 = maturity (Zadoks 91).
Figure 1. Models for 2017: observed vs. predicted external validation values of nitrogen levels from the 2017 dataset. (a) External validation of nitrogen estimates using leaf-level spectral measurements. (b) External validation of nitrogen estimates using canopy-level spectral measurements. Error bars for estimated values represent the standard deviation generated from 500 model permutations using an 80:20 calibration:cross-validation approach. Dashed line: 1:1 relationship. R2: model goodness of fit; RMSE: root-mean-square error (% dry mass); NRMSE: normalized root-mean-square error (%). TP1, time point 1 = booting (Zadoks 45); TP2, time point 2 = onset of heading (Zadoks 50); TP3, time point 3 = heading (Zadoks 58); TP4, time point 4 = flowering (Zadoks 60); TP5, time point 5 = post-anthesis (Zadoks 68); TP6, time point 6 = maturity (Zadoks 91).
Remotesensing 13 03991 g001
Figure 2. Standardized coefficients and variable importance for the projection (VIP) selection for 2017 models. (a,c) Standardized coefficients of 2017 models for leaf-level and canopy-level measurements, respectively; (b,d) VIP of 2017 models for leaf-level and canopy-level measurements, respectively. Shaded areas: regions of absorption feature present between 400 and 2400 nm for chlorophyll (green) and protein (red). Specific absorption features within these regions are listed in Appendix C Table A3.
Figure 2. Standardized coefficients and variable importance for the projection (VIP) selection for 2017 models. (a,c) Standardized coefficients of 2017 models for leaf-level and canopy-level measurements, respectively; (b,d) VIP of 2017 models for leaf-level and canopy-level measurements, respectively. Shaded areas: regions of absorption feature present between 400 and 2400 nm for chlorophyll (green) and protein (red). Specific absorption features within these regions are listed in Appendix C Table A3.
Remotesensing 13 03991 g002
Figure 3. Observed vs. predicted values of nitrogen level estimates for 2018 collections using coefficients from the 2017 model. (a) Estimates generated using leaf-level spectral measurements; (b) estimates generated using canopy-level spectral measurements. Error bars for estimated values represent the standard deviation generated from 500 model permutations using an 80:20 calibration:cross-validation approach. Dashed line: 1:1 relationship. R2: model goodness of fit; RMSE: root-mean-square error; NRMSE: normalized root-mean-square error.
Figure 3. Observed vs. predicted values of nitrogen level estimates for 2018 collections using coefficients from the 2017 model. (a) Estimates generated using leaf-level spectral measurements; (b) estimates generated using canopy-level spectral measurements. Error bars for estimated values represent the standard deviation generated from 500 model permutations using an 80:20 calibration:cross-validation approach. Dashed line: 1:1 relationship. R2: model goodness of fit; RMSE: root-mean-square error; NRMSE: normalized root-mean-square error.
Remotesensing 13 03991 g003
Figure 4. Multi-year models: observed vs. predicted external validation values of nitrogen levels derived from the multi-year dataset. (a) External validation of nitrogen estimates generated using leaf-level spectral measurements. (b) External validation of nitrogen estimates generated using canopy-level spectral measurements. 2017: 2017 dataset, 2018: 2018 dataset. Error bars for estimated values represent the standard deviation generated from 500 model permutations using an 80:20 calibration:cross-validation approach. Dashed line: 1:1 relationship. R2: model goodness of fit; RMSE: root-mean-square error; NRMSE: normalized root-mean-square error.
Figure 4. Multi-year models: observed vs. predicted external validation values of nitrogen levels derived from the multi-year dataset. (a) External validation of nitrogen estimates generated using leaf-level spectral measurements. (b) External validation of nitrogen estimates generated using canopy-level spectral measurements. 2017: 2017 dataset, 2018: 2018 dataset. Error bars for estimated values represent the standard deviation generated from 500 model permutations using an 80:20 calibration:cross-validation approach. Dashed line: 1:1 relationship. R2: model goodness of fit; RMSE: root-mean-square error; NRMSE: normalized root-mean-square error.
Remotesensing 13 03991 g004
Figure 5. Multi-year models: standardized coefficients and variable importance for the projection (VIP) selection for multi-year models. (a,c) Standardized coefficients of the multi-year models for leaf-level and canopy-level measurements, respectively; (b,d) VIP of the multi-year models for leaf-level and canopy-level measurements, respectively. Shaded areas: absorption features present between 400 and 2400 nm for chlorophyll (green) and protein (red). Specific absorption features within these regions listed in Appendix C Table A3.
Figure 5. Multi-year models: standardized coefficients and variable importance for the projection (VIP) selection for multi-year models. (a,c) Standardized coefficients of the multi-year models for leaf-level and canopy-level measurements, respectively; (b,d) VIP of the multi-year models for leaf-level and canopy-level measurements, respectively. Shaded areas: absorption features present between 400 and 2400 nm for chlorophyll (green) and protein (red). Specific absorption features within these regions listed in Appendix C Table A3.
Remotesensing 13 03991 g005
Figure 6. Observed vs. predicted values of NUpE estimated using spectral data. (a) Cross-validated NUpE estimates determined using leaf-level spectral measurements. (b) Cross-validated NUpE estimates determined using canopy-level spectral measurements. Error bars for estimated values represent the standard deviation generated from 500 model permutations using an 80:20 calibration:cross-validation approach. Dashed line: 1:1 relationship. R2: model goodness of fit; RMSE: root-mean-square error; NRMSE: normalized root-mean-square error.
Figure 6. Observed vs. predicted values of NUpE estimated using spectral data. (a) Cross-validated NUpE estimates determined using leaf-level spectral measurements. (b) Cross-validated NUpE estimates determined using canopy-level spectral measurements. Error bars for estimated values represent the standard deviation generated from 500 model permutations using an 80:20 calibration:cross-validation approach. Dashed line: 1:1 relationship. R2: model goodness of fit; RMSE: root-mean-square error; NRMSE: normalized root-mean-square error.
Remotesensing 13 03991 g006
Figure 7. Relationships between nitrogen uptakes (NUpEs) calculated using either reference or estimated nitrogen values. (a) Leaf-level measurements, (b) canopy-level measurements. R2: model goodness of fit; RMSE: root-mean-square error; NRMSE: normalized root-mean-square error. NUpE reference: standard NUpE measurement. NUpE estimated: NUpE calculated using nitrogen predicted values as the input. Dashed line: 1:1 relationship.
Figure 7. Relationships between nitrogen uptakes (NUpEs) calculated using either reference or estimated nitrogen values. (a) Leaf-level measurements, (b) canopy-level measurements. R2: model goodness of fit; RMSE: root-mean-square error; NRMSE: normalized root-mean-square error. NUpE reference: standard NUpE measurement. NUpE estimated: NUpE calculated using nitrogen predicted values as the input. Dashed line: 1:1 relationship.
Remotesensing 13 03991 g007
Table 1. Performance metrics of the 2017 model: number of latent variables (LVs), goodness of fit (R2), root-mean-square error (RMSE), normalized root-mean-square error (NRMSE), and bias for calibration (C), cross-validation (CV), and external validation (EV). Data generated via cross-validation using 500 random permutations of the dataset with 70% of the total data divided into 80:20 split for C and CV, respectively, for models predicting nitrogen concentrations from wheat foliage using hyperspectral data. EV was generated by applying coefficients generated in the C:CV approach to the 30% of the data excluded from the C:CV approach. Data shown as mean ± standard deviation.
Table 1. Performance metrics of the 2017 model: number of latent variables (LVs), goodness of fit (R2), root-mean-square error (RMSE), normalized root-mean-square error (NRMSE), and bias for calibration (C), cross-validation (CV), and external validation (EV). Data generated via cross-validation using 500 random permutations of the dataset with 70% of the total data divided into 80:20 split for C and CV, respectively, for models predicting nitrogen concentrations from wheat foliage using hyperspectral data. EV was generated by applying coefficients generated in the C:CV approach to the 30% of the data excluded from the C:CV approach. Data shown as mean ± standard deviation.
Data LevelLVR2RMSENRMSE (%)Bias
CCVEVCCVEVCCVEVCCVEV
Leaf50.90 ± 0.000.84 ± 0.020.710.24 ± 0.010.31 ± 0.020.45810144.10 × 10−18 ± 0.02−0.0002 ± 0.05−0.009
Canopy60.87 ± 0.010.85 ± 0.020.730.27 ± 0.010.3 ± 0.020.4299134.89 × 10−18 ± 0.02−0.0015 ± 0.040.003
Table 2. Multi-year models’ performance: number of latent variables (LVs), model goodness of fit (R2), root-mean-square error (RMSE), normalized root-mean-square error (NRMSE), and bias for calibration (C), cross-validation (CV), and external validation (EV). Data generated via cross-validation using 500 random permutations of the dataset with 70% of the total data divided into 80:20 split for C and CV, respectively, for models predicting nitrogen concentrations from wheat foliage using hyperspectral data. EV was generated by applying coefficients generated from the C:CV approach to the 30% of the data excluded from the C:CV approach. Data shown as mean ± standard deviation.
Table 2. Multi-year models’ performance: number of latent variables (LVs), model goodness of fit (R2), root-mean-square error (RMSE), normalized root-mean-square error (NRMSE), and bias for calibration (C), cross-validation (CV), and external validation (EV). Data generated via cross-validation using 500 random permutations of the dataset with 70% of the total data divided into 80:20 split for C and CV, respectively, for models predicting nitrogen concentrations from wheat foliage using hyperspectral data. EV was generated by applying coefficients generated from the C:CV approach to the 30% of the data excluded from the C:CV approach. Data shown as mean ± standard deviation.
Data LevelLVR2RMSENRMSE (%)Bias
CCVEVCCVEVCCVEVCCVEV
Leaf60.89 ± 0.010.84 ± 0.010.720.01 ± 0.010.32 ± 0.010.4271012−1.2 × 10−17 ± 0.00−0.0043−0.06
Canopy50.87 ± 0.010.84 ± 0.020.670.26 ± 0.010.29 ± 0.020.46910135.98 × 1018 ± 0.00−0.00160.04
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Peron-Danaher, R.; Russell, B.; Cotrozzi, L.; Mohammadi, M.; Couture, J.J. Incorporating Multi-Scale, Spectrally Detected Nitrogen Concentrations into Assessing Nitrogen Use Efficiency for Winter Wheat Breeding Populations. Remote Sens. 2021, 13, 3991. https://0-doi-org.brum.beds.ac.uk/10.3390/rs13193991

AMA Style

Peron-Danaher R, Russell B, Cotrozzi L, Mohammadi M, Couture JJ. Incorporating Multi-Scale, Spectrally Detected Nitrogen Concentrations into Assessing Nitrogen Use Efficiency for Winter Wheat Breeding Populations. Remote Sensing. 2021; 13(19):3991. https://0-doi-org.brum.beds.ac.uk/10.3390/rs13193991

Chicago/Turabian Style

Peron-Danaher, Raquel, Blake Russell, Lorenzo Cotrozzi, Mohsen Mohammadi, and John J. Couture. 2021. "Incorporating Multi-Scale, Spectrally Detected Nitrogen Concentrations into Assessing Nitrogen Use Efficiency for Winter Wheat Breeding Populations" Remote Sensing 13, no. 19: 3991. https://0-doi-org.brum.beds.ac.uk/10.3390/rs13193991

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