JEQ Grow Your Career With ASA
HOME HELP FEEDBACK SUBSCRIPTIONS ARCHIVE SEARCH TABLE OF CONTENTS
 QUICK SEARCH:   [advanced]


     


This Article
Right arrow Abstract Freely available
Right arrow Figures Only
Right arrow Full Text (PDF)
Right arrow Alert me when this article is cited
Right arrow Alert me if a correction is posted
Services
Right arrow Similar articles in this journal
Right arrow Similar articles in ISI Web of Science
Right arrow Similar articles in PubMed
Right arrow Alert me to new issues of the journal
Right arrow Download to citation manager
Citing Articles
Right arrow Citing Articles via HighWire
Right arrow Citing Articles via ISI Web of Science (17)
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Read, J. J.
Right arrow Articles by Reddy, K. R.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by Read, J. J.
Right arrow Articles by Reddy, K. R.
Agricola
Right arrow Articles by Read, J. J.
Right arrow Articles by Reddy, K. R.
Related Collections
Right arrow Crop Physiology & Metabolism
Right arrow Cotton
Right arrow Remote Sensing
Right arrow Plant Nutrition
Journal of Environmental Quality 31:1442-1452 (2002)
© 2002 American Society of Agronomy, Crop Science Society of America, and Soil Science Society of America

Narrow-Waveband Reflectance Ratios for Remote Estimation of Nitrogen Status in Cotton

J. J. Read*,a, L. Tarpleyb, J. M. McKiniona and K. R. Reddyc

a USDA-Agricultural Research Service, Crop Science Research Laboratory, P.O. Box 5367, Mississippi State, MS 39762
b Texas A&M Research and Extension Center, 1509 Aggie Dr., Beaumont, TX 77713
c Plant and Soil Science Dep., Box 9555, Mississippi State, MS 39762

* Corresponding author (jjread{at}msa-msstate.ars.usda.gov)

Received for publication January 17, 2001.

    ABSTRACT
 TOP
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 NOTES
 RESULTS AND DISCUSSION
 CONCLUSIONS
 REFERENCES
 
Tailoring nitrogen (N) fertilizer applications to cotton (Gossypium hirsutum L.) in response to leaf N status may optimize N use efficiency and reduce off-site effects of excessive fertilizer use. This study compared leaf and canopy reflectance within the 350 to 950 nm range in order to identify reflectance ratios sensitive to leaf chlorophyll (Chl), and hence N status, in cotton. Plants were grown outdoors in large pots using half-strength Hoagland's (control) solution until some three-row plots received a restricted supply of N. Treatments comprised control, 20% of control N at first flower bud (square) onward; 0 and 20% of control N at first flower onward; and 0% of control N at fruit-filling onward. Despite leaf N values ranging from 51 to 19 g kg-1 across treatments and sampling dates, a weak correlation was obtained between Chl and N (r2 = 0.32, df = 70). In general, N stress led to increased reflectance at 695 ± 2.5 nm (R695) and decreased reflectance at R410, and changes in leaf N were best correlated with either R695 or R755 in leaves and either R410 or R700 in canopies. The strongest associations between leaf constituent and canopy reflectance ratio were Chl vs. R415/R695 (r2 = 0.72), carotenoids vs. R415/R685 (r2 = 0.79), and N vs. R415/R710 (r2 = 0.70). The R415 measure appears to be a more stable spectral feature under N stress, as compared with more pronounced changes along the reflectance red edge (690–730 nm). Multiple regression identified a three-waveband canopy reflectance model that explained 80% of the variability in leaf N. Results indicate that remote sensing of N status in cotton is feasible using narrow-waveband reflectance ratios that involve the violet or blue region of the spectrum (400 to 450 nm) and the more commonly featured red-edge region.

Abbreviations: Chl, chlorophyll • NIR, near infrared range of the spectrum • R, leaf or canopy reflectance


    INTRODUCTION
 TOP
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 NOTES
 RESULTS AND DISCUSSION
 CONCLUSIONS
 REFERENCES
 
NITROGEN (N) FERTILITY is an important determinant in cotton productivity, quality, and profitability (Gerik et al., 1998). Crop N requirements in southeastern USA are relatively high and variable, ranging between 67 and 225 kg N ha-1 (or about 0.12 kg per kg lint) (Boquet et al., 1993). Meredith et al. (1997) reported cotton growers in the Mississippi Delta region apply about 112 kg ha-1, and found modern, high-yielding varieties may require larger amounts of nutrients in a shorter time period than older varieties. Moreover, because N supply influences vegetative and reproductive growth, there is a tendency by some growers to boost yield potentials by applying higher than recommended N rates (Boquet and Breitenbeck, 2000). Consequently, modifications to current N fertilization and irrigation practices have potential for reducing runoff and leaching of mobile forms of nitrogen (nitrite and nitrate), and thus minimizing environmental effects of cotton production on surface and ground water quality (Pan et al., 1997; Power et al., 2000). Some approaches to this problem include the use of crop simulation models (Hunt et al., 1998; Jones and Barnes, 2000; McKinion et al., 2001) and remote sensing for precision crop management (PCM) in cotton production (Barnes et al., 2000).

Remote sensing in the visible and near-infrared (NIR) wavelengths (between about 400 and 1000 nm) is based on the principle that changes in plant light interception and utilization influence the reflectance properties of vegetation (Carter, 1998; Jensen, 2000). Nitrogen deficiency in cotton decreases chlorophyll (Chl) and soluble protein content (Longstreth and Nobel, 1980), as well as the rate of leaf expansion and canopy development (Reddy et al., 1997; Gerik et al., 1998). Remote sensing of Chl has potential for rapidly estimating cotton N status, and hence crop productivity. A critical need in PCM is the unambiguous detection of plant N status, preferably at an early stage of crop development, which requires identification of spectral wavebands or indices in which vegetation reflectance is most responsive to unfavorable growth conditions (Carter, 1994).

The dominant biochemical factor affecting canopy reflectance in the visible wavelengths is the concentration and total amount of leaf chlorophylls a and b (Chappelle et al., 1992; Filella et al., 1995). Using a priori knowledge of leaf physiology, some authors have related reflectance at 675 nm (R675), where Chl a absorption is maximal, to the concentration of total Chl (Jacquemoud and Baret, 1990; Filella et al., 1995). Chappelle et al. (1992) found R675/R700 in ratio to the value obtained from well-fertilized control plots yielded a precise relationship with Chl content in N-deficient soybean [Glycine max (L.) Merr.]. Others have proposed using R550 (the green peak where Chl a and b absorption is minimal), because relatively low leaf reflectance at 675 nm decreases sensitivity of the 675-nm band to changes in Chl a. Thomas and Gausman (1977) examined single leaf reflectance in three principle wavebands, and reported 550 nm was superior to either 450 or 675 nm for relating leaf reflectance to pigment concentration in cotton. Reflectance at 550 nm or its ratio with a waveband in the near-infrared was associated with corn (Zea mays L.) grain yield, an estimator of N status (Blackmer et al., 1996), and with chlorophyll loss during leaf senescence in horse chestnut (Aesculus hippocastanum L.) and Norway maple (Acer platanoides L.) (Gitelson and Merzlyak, 1996).

The position of the red edge, which is the sharp increase in spectral reflectance between 690 and 730 nm, may shift to shorter wavelengths as the concentration of Chl a decreases (Horler et al., 1983; Carter, 1994). This stress-sensitive region in digital imagery is often referred to as the "blue shift" (Jensen, 2000). Because the reflectance red edge in cotton leaves is known to increase under N stress (Tarpley et al., 2000), this spectral feature has potential as an index for the remote sensing of crop canopies deficient in N. Increased reflectance of approximately R700 results from strong absorption by Chl in the visible region, near 670 nm, and relatively high reflectance by leaf structural features in the near-infrared region, between 720 and 920 nm. Several studies have related the red edge, or single-waveband ratios involving the red edge, to changes in Chl concentration (Gitelson and Merzlyak, 1996), photosynthetic activity (Carter et al., 1996; Carter, 1998), or micronutrient content (Adams et al., 2000). The use of simple reflectance ratios may correct for variations in irradiance, leaf orientation, irradiance angles, and shading (Carter, 1994). In leaves of cotton, single-waveband ratios that combined a low reflectance, red-edge measure with a waveband of high reflectance in the NIR region (755 to 920 and 1000 nm) were very sensitive indicators of N concentration (Tarpley et al., 2000).

Leaf N may also be associated with the concentration of other pigments. Carotenoids are a group of yellow to orange pigments involved in light harvesting for photosynthesis, and may also serve to protect chlorophyll from O2 damage when irradiance levels or other stressors are excessively high. The ratio of carotenoids to Chl a typically increases under stress (Thomas and Gausman, 1977) or in senescing leaves (Gitelson et al., 1996). Normalized difference pigment index [NPCI, (R680 - R430)/(R680 + R430)] and pigment simple ratio (PSR, R430/R680) appear to be sensitive reflectance indices of environmental and phenological changes in the carotenoid to Chl a ratio (Filella et al., 1995).

Besides leaf constituent absorption, vegetation reflectance is affected by the contribution of stems and leaf orientation to canopy reflectance (Carter and Miller, 1994; Jensen, 2000). Therefore, it is generally believed that in situ measurements with a hand-held instrument may provide more direct and accurate assessments of leaf Chl. Moreover, leaf spectra are useful in explaining how canopy spectra are influenced by changes in leaf reflectance as compared with changes in canopy structure. Radiative transfer models, which are beyond the scope of this paper, can be used to simulate the spectral contribution of leaf biochemical constituents to the canopy reflectance (e.g., Jacquemoud and Baret, 1990).

The present study determined optical properties in full canopies of cotton (i.e., without soil background) that best correlate with leaf Chl and N concentration. Because N deficiency alters cotton growth and development, reflectance ratios were compared between measurements obtained at the level of single leaves and full canopies. Measurements were conducted under natural illumination conditions to directly relate vegetation reflectance with remote (e.g., airborne) observations. Simple waveband ratios were evaluated as indicators of Chl content, and hence the activity of green vegetation. Linear response functions were sought across a wide range of vegetation and phenological conditions, so the reflectance ratios would be more generally applicable, and to facilitate calibration and validation of the index. Plant N status, and hence N stress, was ascertained from total leaf N values, because they provide a reliable estimate of N accumulated by the leaf prior to sampling (Gerik et al., 1998). We expect reflectance-based algorithms that are sensitive to N deficiency symptoms would provide functionality in image processing of a cotton plantation, watershed, or an entire county, which along with ground truth data on soil chemical and physical properties can be used to generate a site-specific soil fertility map. The intended outcome is to apply the necessary amounts of fertilizer N when and where it is needed to optimize cotton N management and yield.


    MATERIALS AND METHODS
 TOP
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 NOTES
 RESULTS AND DISCUSSION
 CONCLUSIONS
 REFERENCES
 
Experimental Facility and Treatment Design
Individual cotton (cv. NuCOTN 33B) plants were grown outdoors in large, free-draining polyvinyl chloride (PVC) pots (15-cm diam. x 65-cm depth) filled with sand and supplied water and nutrients via plastic pipe and a dripper system (Netafim,1 Fresno, CA). The pots were set side-by-side in 6.5-m-long racks (or rows) spaced 1 m apart and arranged in an east–west direction. Seeds were sown 17 May 1999, and plants were thinned to one per pot 11 June 1999. Each treatment comprised three 1-m-wide rows with plants spaced 0.15 m on center, giving 40 plants per row. This study was part of a larger study that comprised a total of 40 rows.

All control and border rows were fed a favorable supply of water and nutrients using half-strength Hoagland's nutrient solution via one Netafim dripper per pot rated at 1 L of liquid per hour. Plants were irrigated three times per day. If most plants were visibly wilted at midday, the duration of subsequent daily irrigations was increased to maintain soil water above the wilting point. A weather station adjacent to the field site continuously monitored radiation, precipitation, air temperature, and relative humidity, which aided in determining plant water requirements. The different N solutions for each treatment were contained in large mixing tanks and pumped under pressure by chemical metering pumps through plastic lines to the three-row treatment plots. All timing and duration of flows to individual plots was under computer-controlled switches and solenoid valves. At various times (listed below), the plants in some three-row plots were supplied nutrients with an osmotically balanced Hoagland's solution containing either 0 or 20% of control N. The five N treatments were: (i) half-strength Hoagland's supplied either up to and beyond the time N treatment was imposed (control), (ii) 20% of control N commenced 31 d after emergence at first floral bud (square) stage, (iii) 20% of control N commenced 52 d after emergence at flowering stage, (iv) 0% of control N commenced 52 d after emergence at flowering stage, and (v) 0% of control N commenced 89 d after emergence at fruit-filling stage.

Spectroradiometric and Physiological Measurements
Hyperspectral reflectance data from 350 to 950 nm with a resolution of about 1.5 nm were obtained from single leaves and full canopies (i.e., no soil background) with a GER 1500 spectroradiometer (Geophysical and Environmental Research Corp., Millbrook, NY). Measurements were taken at midday (1100 to 1300 h) under incident solar radiation from an angle perpendicular to the target. Four scans were averaged in every measurement. Immediately (2 min) prior to each set of reflectance measurements, a reference spectra was collected from a high reflectance, white reference panel (Spectralon; Labsphere, North Sutton, NH). Percent reflectance was computed by dividing observed radiance data by radiance measured from a white reference panel under approximately the same sun conditions, and multiplying by 100.

Leaf reflectance was measured from the uppermost, fully expanded leaves of five plants selected at random from each three-row treatment. Measurements were conducted with a 4° field of view foreoptic from a height of 0.3 m above the leaf. Thus, spectral reflectance was detected from an approximately 9.5-cm2 area (3.45 cm long x 2.75 cm wide) of the middle lobe. Because a fully absorptive (black) background was not placed behind a leaf during the measurement, we expect near-infrared reflectance is artificially enhanced, as compared with measurements obtained using an integrating sphere. This additive reflectance is due to the fact that healthy green leaves in the vicinity of our field of view may reflect approximately 80% of the incident near-infrared energy at 900 nm (Jensen, 2000). We minimized certain limitations to measuring reflectance under natural illumination conditions, such as specular reflectance (e.g., glare) and variability in irradiance angle and intensity by (i) positioning the radiometer between the sun and the target of the upper leaf (adaxial) surface, (ii) selecting the best angle with minimal specular reflectance, (iii) avoiding shadows, and (iv) completing measurements in <1 h.

Canopy reflectance was measured with a 23° field of view foreoptic from a height of 1 m above the central portion of a canopy (40-cm width). Thus, reflectance typically involved fully developed plant crowns arising from three consecutive pots, and virtually no soil background. A fiber optic cable was run to the radiometer along an extensible boom that attached to the top of a tripod, which also had a platform for attaching the spectroradiometer. Three separate rows were measured in each treatment by quickly positioning the boom as needed above each row. A nadir view (single viewing geometry) may be insufficient to fully characterize N stress in the cotton canopy. This is because measurements do not reveal different aspect of vertical structure, leaf area index, leaf orientation, and shading that have a strong influence on overall characteristics of the reflectance spectrum (Jensen, 2000). We expect these effects may largely cancel out as the number of measurements is increased, as in the present study. Further, our main intent was to compare reflectance indices between single leaves and canopies that correspond most strongly with Chl concentration. Further studies may require that the comparison of leaf and canopy reflectance be evaluated with radiative transfer models (Jacquemoud and Baret, 1990).

Immediately following leaf reflectance measurements, five uppermost, fully expanded leaves were harvested from each treatment, placed in a cooler over ice, and transported to the laboratory. A cork borer was used to remove five discs from each leaf, avoiding areas with large veins. Leaf discs were placed into vials containing 4.0 mL dimethyl sulfoxide (DMSO), and kept at room temperature overnight in the dark to extract pigments. Absorbances of the extract at 470, 648, and 664 nm were recorded and concentrations of carotenoid pigments and chlorophylls a, b and a + b were computed following the formulae of Chappelle et al. (1992). Chlorophyll concentration of the extract and the total disk surface area of 0.98 cm2 were used to compute total leaf Chl concentrations per unit projected leaf area. Fully expanded uppermost leaves were sampled for Chl because interaction with the incident light and proximity to the spectroradiometer would strongly affect canopy reflectance, and because these leaves express early symptoms of N stress, such as chlorosis. The five leaves used for pigment analysis were pooled and dried at 70°C for 72 h. Leaf tissue N concentration was subsequently determined on dried, ground samples according to standard micro-Kjeldahl methods (Nelson and Sommers, 1972), and expressed as either g kg-1 dry weight or g cm-2 leaf area.

Single-leaf gas exchange was determined during the boll filling period on 30 July and 5 August with a LI-COR (Lincoln, NE) Model LI-6400 photosynthesis system. Uppermost, fully expanded leaves were measured from five plants selected at random within each treatment, except 0% N at boll-filling onward, giving 20 observations per sampling date. Irradiance from a red–blue LED was set at 1500 µmol m-2 s-1 and ambient CO2 concentration was controlled at 350 µmol mol-1. Environmental conditions inside the leaf chamber were typically 27°C air temperature, 55% relative humidity, with a leaf-to-air vapor pressure deficit of about 1.4 kPa.

Data Analysis
A total of 73 observations for leaf pigments and 72 observations for leaf N concentration were generated from samples collected at approximately biweekly intervals. These data were analyzed for significant treatment effects (P < 0.05) across sampling dates with the proc MIXED procedure in SAS (SAS Institute, 1999). A balanced set of leaf constituent data (n = 69) is presented in Fig. 1 . This figure also shows the sampling dates when both leaf and canopy reflectance were measured, which represent leaf Chl and N data used in the correlation and regression analyses (described below).



View larger version (29K):
[in this window]
[in a new window]
 
Fig. 1. Evolution of (a) total chlorophyll, (b) carotenoids, and (c) nitrogen in uppermost fully expanded leaves of cotton, n = 69. Vertical arrows indicate sampling dates on which both leaf and canopy reflectance were measured (n = 22); these data were subsequently used in correlation analysis (see Fig. 3, 4, 5, and 6). Leaf and canopy reflectance measurements were obtained only once, on 23 August, in the 0% of control N at boll-filling onward treatment.

 


View larger version (23K):
[in this window]
[in a new window]
 
Fig. 3. Coefficients of determination (r2) for the linear relationship between (a) leaf reflectance and leaf chlorophyll concentration, (b) reflectance at 755 nm divided by reflectance at wavelength indicated by the horizontal axis, and (c) reflectance at 705 nm divided by each wavelength. All plot means of spectral reflectance and chlorophyll collected during the two-month study (n = 60) were included in correlation analyses.

 


View larger version (21K):
[in this window]
[in a new window]
 
Fig. 4. Coefficients of determination (r2) for the linear relationship between (a) leaf reflectance and total leaf nitrogen concentration, (b) reflectance at 755 nm divided by reflectance at wavelength indicated by the horizontal axis, and (c) reflectance at 700 nm divided by each wavelength. All plot means of spectral reflectance and chlorophyll collected during the two-month study (n = 60) were included in the correlation analyses.

 


View larger version (21K):
[in this window]
[in a new window]
 
Fig. 5. Coefficients of determination (r2) for the linear relationship between (a) canopy reflectance and total leaf chlorophyll concentration, (b) reflectance at 700 nm divided by reflectance at wavelength indicated by the horizontal axis, and (c) reflectance at 415 nm divided by reflectance at each wavelength. All plot means of spectral reflectance and chlorophyll collected during the two-month study (n = 22) were included in correlation analyses.

 


View larger version (21K):
[in this window]
[in a new window]
 
Fig. 6. Coefficients of determination (r2) for the linear relationship between (a) canopy reflectance and total leaf nitrogen, (b) reflectance at 415 nm divided by reflectance at wavelength indicated by the horizontal axis, and (c) reflectance at 710 nm divided by reflectance at each wavelength. All plot means of spectral reflectance and chlorophyll collected during the two-month study (n = 22) were included in the correlation analyses.

 
We developed a SAS (SAS Institute, 1999) program to determine single-waveband ratios most strongly associated with N status. Spectra were originally sampled in approximately 1.5-nm intervals. First, reflectance values were averaged across non-overlapping, 5-nm intervals to reduce noise, and still retain potentially important narrow-waveband features. Then, the spectra were averaged across replicate leaf (n = 5) or canopy (n = 3) measurements, each represented by a three-row plot in each treatment. Pearson's correlation coefficient (r) was calculated across all treatments and sampling dates to determine associations between reflectance (or reflectance ratios) measured in narrow wavebands and the concentration of Chl and N in leaves on a plot mean basis. We first plotted the coefficients of determination (r2) against individual wavebands (5-nm resolution) to determine spectral regions of greatest sensitivity to leaf N status (e.g., Fig. 3a). Then, the maximum r2 value and, as needed, a priori knowledge of leaf physiology, was used to determine wavebands from which a numerator reflectance should be selected for computing simple ratios that used reflectance values at 5-nm intervals between 350 and 950 nm as the denominator (e.g., Fig. 3b). Similarly, a final iteration involved all possible reflectance ratios formed when the denominator waveband best correlated with leaf N status was used as the numerator in ratio with the reflectance spectra at 5-nm intervals (e.g., Fig. 3c). Values for r2 from the linear relationship between a reflectance ratio and leaf constituent were plotted against the individual denominator wavebands. The r2 value indicates the fraction (0–1) of variability in leaf N status that can be accounted for by plant reflectance features. Linear regression analysis was used to develop relationships between the waveband ratio with maximum r2 (or potential to maximize r2) and changes in leaf N concentration across the different treatments and sampling dates.

We also used multiple linear regression to determine the combination of spectrally narrow canopy reflectance wavebands (1.5 nm resolution) that best explains variation in leaf N concentration. Because of the potential for inflation of r2 values (over-fitting data) when independent variables are autocorrelated, as is the case for spectral reflectance data (Thenkabail et al., 2000; Grossman et al., 1996), we used a partial least squares (PLS) procedure in SAS as a selection methodology (SAS Institute, 1999). Because PLS provides a list of all independent variables in the model, the first step involved ranking the 360 single wavebands based on two SAS outputs, the standardized regression coefficients and the Variable Importance for Projection (VIP) (Lindgren et al., 1994). Whereas the regression coefficients represent the importance each predictor has in the prediction of just the response, the VIP represents the value of each predictor in fitting the PLS model for both predictors and response. Sixty of the highest-ranked wavebands were then selected as inputs in the SAS-MAXR procedure, so the model would have fewer predictor variables (single wavebands) than dependent variables (66 observations for leaf N), and thus minimize the potential of over-fitting. We then selected the best 20-term model based on MAXR regression analysis. Then, the next 40 wavebands identified in the PLS ranking were included in the MAXR model statement, and the regression analysis was repeated, keeping the best 20-term model each time. This step was repeated until the PLS ranking list was exhausted. The best combination of wavebands was obtained from the final MAXR iterative analysis.


    RESULTS AND DISCUSSION
 TOP
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 NOTES
 RESULTS AND DISCUSSION
 CONCLUSIONS
 REFERENCES
 
Treatment and Seasonal Effects on Leaf Nitrogen Status and Photosynthesis
In plants grown under N restriction from either first square onward or first flower onward, N deficiency decreased the Chl, carotenoid, and N concentrations in uppermost fully expanded leaves (Fig. 1). Analysis of variance across all sampling dates (n = 73) indicated significant (P < 0.001) differences in these physiological traits between control and N-deficient plants. Averaged across all sampling dates, mean values for total Chl were 53.3, 47.6, 45.4, 37.6, and 36.9 µg cm-2 in controls, 0% N at boll-filling, 20% N at first flower, 20% N at first square, and 0% N at first flower, respectively. Mean values for carotenoids in this set of treatments were 12.5, 11.0, 10.9, 9.5, and 9.6 µg cm-2, respectively. A strong (P < 0.001) positive association was obtained between Chl and carotenoid concentration (r2 = 0.80, df = 72). This indicates cotton response to N stress in this study involved a general loss of photosynthetic pigments, which would decrease maximum light absorption in uppermost, fully expanded leaves. We further expect the observed stress-induced changes in spectral reflectance (discussed below) were largely influenced by the 20% N at first square and 0% N at first flower onward treatments, as these plants had relatively rapid and large reductions in leaf Chl and carotenoid concentrations (Fig. 1a,b).

Consistent with other N-deficit studies in cotton (Reddy et al., 1997), leaf N reached a minimum of about 18.7 g kg-1 dry weight on 9 July, and then remained level in plants grown at 20% control N from first square stage onward. Averaged across sampling dates, leaf N values were 41.2, 35.3, 31.3, 26.2, and 26.3 g kg-1 in controls, 0% N at boll-filling, 20% N at first flower, 20% N at first square, and 0% N at first flower, respectively (Fig. 1c). Gerik et al. (1998) reported that cotton leaf N levels > 45 g kg-1, 30 to 45 g kg-1, 25 to 30 g kg-1, and <25 g kg-1, respectively, are considered excessive, sufficient, low, and deficient. Deficient levels of N were recorded on most sampling dates among plants supplied 20% of control N at first square stage onward; whereas leaves from the 0% N at first flower treatment ranged from low to deficient in N across sampling dates (Fig. 1c).

Leaf N was weakly associated with either Chl a (r2 = 0.37, df = 70) or total Chl (r2 = 0.32); however, the relationship with Chl improved somewhat when N was expressed as g m-2 leaf area (r2 = 0.42; df = 49; some samples were ground prior to weighing). Luxury consumption of N may weaken the relationship between Chl and N, but was of little consequence in the present study because relatively few values for leaf N exceeded 45 g kg-1 (Gerik et al., 1998). Further, leaf N values ranged from 18.7 to 56.0 g kg-1, which is very similar to the range reported by Tarpley et al. (2000) for either growth-chamber plants with restricted N supply or field-grown plants with favorable nutrient conditions. Earlier studies indicate cotton photosynthesis would change about 50% across this range of deficient to excessive leaf N values, with both linear and quadratic responses being observed (Reddy et al., 1997; Gerik et al., 1998).

Leaf photosynthetic rates measured during the boll-filling period increased as Chl concentration increased across the four N treatments (Fig. 2) . Photosynthesis was lower on 5 August than 30 July in the present study, and was associated with a 50% reduction in leaf conductance to water vapor due in part to somewhat greater leaf-to-air vapor pressure deficit (1.3 vs. 1.6 kPa). Longstreth and Nobel (1980) demonstrated that photosynthetic potential in cotton relates strongly to Chl concentration, which was directly correlated with plant mineral status, notably N and potassium, and to a lesser extent phosphorus.



View larger version (25K):
[in this window]
[in a new window]
 
Fig. 2. Relationship between net CO2 exchange and leaf chlorophyll concentration in uppermost leaves of cotton grown under N-limiting conditions. Values represent single observations obtained from five leaves in each treatment and sampling date.

 
Final dry matter accumulation and lint yield were least in plants grown under N restriction from either first square or first flower onward (data not presented). These N-deficit treatments were affected most by N stress because it is at these developmental stages that cotton N resources are being shifted to increased growth and reproductive functions, respectively (Reddy et al., 1997; Boquet et al., 1993). This seasonal change in cotton N requirement is well known, and has prompted many growers to make single or multiple split applications of N based on crop progress and yield goals, and thereby maintain sufficient levels of N in leaves (Boquet and Breitenbeck, 2000; Hunt et al., 1998; Gerik et al., 1998). We expect remote sensing of Chl has potential as an early indicator of these N-deficient zones within the field prior to visible N-stress symptoms.

Spectral Reflectance Properties in Single Leaves and Full Canopies
When plot means were analyzed across all N treatments and sampling dates, changes in Chl were best correlated with single-leaf reflectance values centered at 755 nm (r2 = 0.15), and to a lesser extent 695 nm (Fig. 3a) . Because R755 is near a sun-induced fluorescence peak at approximately 760 nm, stress-induced changes in R755 may be influenced by factors other than Chl content. Nevertheless, Theisen (2000) reported that ratios utilizing far-red fluorescence at 755.56 nm consistently yielded strong correlation with leaf Chl in N-stressed bean (Phaseolus vulgaris L.). Sensitivity of a red edge feature at 695 nm is less tenuous, and represents the transition (or degree of slope) between maximal Chl absorption at approximately R675 and minimal absorption of light at approximately R700.

Because no further absorption occurs beyond approximately 750 nm, NIR reflectance becomes primarily a function of leaf structure and scatter. This suggests N-deficit treatments in the present study led to changes in both the pigmentation and internal structure of leaves. However, caution must be used when interpreting results for any single-waveband reflectance feature measured under naturally sunlit conditions. This is because reflectance at one determined wavelength is sensitive to the geometric arrangement between the leaf and canopy surface and sensor (the viewing geometry), variable irradiance, background effects, and green biomass. Carter (1994) proposed computing single-waveband reflectance ratios between a "stress sensitive" and "stress-insensitive" band as a way to correct for these variable effects, and subsequent studies have confirmed the validity of this method (Carter and Spiering, 2000, and references therein). Single-waveband ratios are mathematically similar to commonly used vegetation indices (Jensen, 2000), but can be optimized for detection of leaf N status (Carter and Spiering, 2000).

The association between Chl and leaf reflectance improved substantially when reflectance at 755 nm, a relatively Chl-insensitive waveband, was divided by reflectance values centered at either 705 or 555 nm (Fig. 3b). Reflectance at 705 nm is inversely related to changes in Chl a content, whereas R750 serves as an internal reference with respect to spectral reflectance due to leaf internal structure. In previous studies, R700/R750 was linearly related to Chl content in tobacco (Nicotiana tabacum L.) genotypes (Lichtenthaler et al., 1996), and in senescing leaves of two deciduous tress species (Gitelson et al., 1996). Sensitivity at R700 is a red-edge measure that typically occurs in the 680 to 760 nm range, and appears to be more strongly dependent on Chl content than reflectance near 550 nm (Gitelson and Merzlyak, 1996; Carter et al., 1996; Barnes et al., 2000; Carter and Spiering, 2000). Reflectance at about 550 nm or its ratio with a NIR waveband was closely associated with plant N status in corn (Blackmer et al., 1996), as well as two deciduous tree species with relatively low N content (Gitelson and Merzlyak, 1996). Overall, the best correlation with Chl in the present study (approximate r2 = 0.56) was obtained for leaf reflectance ratios R705/R930 and R705/R715 (Fig. 3c).

Similar to changes in Chl in our study, leaf N was best correlated with leaf reflectance values centered at 755 nm (r2 = 0.13), and to a much lesser extent at 695 nm (Fig. 4a) . The association with N improved substantially when R755 was divided by reflectance at either 700 or 590 nm (Fig. 4b), but the best correlation was obtained for reflectance ratios R700/R710 (r2 = 0.38) and R700/R755 (r2 = 0.34) (Fig. 4c). Because R700/R710 is essentially a feature of the red edge, we expect both wavebands are sensitive to Chl a, and sensitivity of the ratio is a result of variability in the reflectance response to early loss of Chl (chlorosis) in young, green leaves (see Fig. 2 in Gitelson et al., 1996). The proximity of these two wavebands would present difficulties in radiometer construction. Moreover, the commonly reported ratio of red edge to NIR reflectance, either R705/R930 for total Chl (Fig. 3c) or R700/R755 for N (Fig. 4b), was consistently identified as a sensitive spectral feature of N status.

Our success in identifying spectral indices sensitive to leaf N in cotton follows from the strong association between leaf Chl and N content observed in this and other nutrient stress studies (Longstreth and Nobel, 1980). For instance, the predominance of spectral features along the red edge in the present study is due to a combination of strong absorption by Chl near 670 nm, and high reflectance by leaf structural features between 720 and 920 nm (Jensen, 2000; Gitelson et al., 1996; Chappelle et al., 1992). In previous studies, single-leaf reflectance near 700 nm or in ratio with wavebands in the near-infrared region were closely associated with stress-induced changes in leaf N (Tarpley et al., 2000). They found that about half of the single-band reflectance ratios related to N concentration had r2 values of 0.64 or greater, and ranged up to 0.92. The best leaf reflectance ratios identified in the present study explained only about 35 to 40% of the variation in leaf N (Fig. 4).

Unexpectedly, treatment and seasonal changes in leaf Chl and N concentration were more closely associated with spectral reflectance from canopies (Fig. 5 and 6) than from leaves (Fig. 3 and 4). Further, leaf-level and canopy-level reflectance differed in the spectral bands with maximal sensitivity to plant N status. Presumably, these differences are due to additional factors, such as the influence of shadows, leaf density and orientation, and stems and/or flowers on canopy optical properties under N-limiting conditions (Carter et al., 1996; Barnes et al., 2000). In regard to single wavebands, the strong correlation observed between leaf N status and near-infrared reflectance (approximately 755 nm) in leaves was absent in canopies (Fig. 5 and 6). Multiple leaf layers reflect strongly in the near-infrared region due to additive reflectance (Jensen, 2000), but our results indicate N deficiency had little affect on either the density of leaves or number of leaf layers in cotton. The association in leaves between about 560 and 590 nm apparently shifted to higher wavebands between about 605 and 625 nm in canopies (Fig. 3a, 4a, 5a, and 6a), a region closely related to Chl a absorption (Chappelle et al., 1992). As with leaf reflectance data, canopy reflectance was best correlated with leaf Chl concentration at 700 nm (r2 = 0.52), but also at 625 nm (r2 = 0.43) (Fig. 5).

The canopy spectral association with Chl improved somewhat when reflectance at 700 nm was divided by reflectance at 415 nm (r2 = 0.67) or 675 nm (r2 = 0.50) (Fig. 5b). Because Chl a and carotenoid pigments absorb strongly at 420 nm (Chappelle et al., 1992), a relatively large amount of Chl must be lost before reflectance increases significantly at this wavelength. This suggests canopy ratios involving R415 are ratioing to more stable regions of the spectrum, as compared with more pronounced changes along the red edge under N stress. At 605 nm, and particularly at 695 and 710 nm, the absorptivity of Chl a and b is relatively weak (Lichtenthaler et al., 1996), approaching zero at wavelengths near 720 nm when Chl remains associated with chloroplast membranes. Thus, as Chl begins to decrease under N stress, leaf reflectance increases first and then most dramatically in the red edge due to further narrowing of the major Chl a absorption peak at 680 nm (Carter et al., 1996). Because canopy reflectance ratios (or indices) include the effect of both biomass and pigments, the significance of R675 is likely a function of strong absorption by leaves between 660 and 690 nm due to changes in Chl (Chappelle et al., 1992). Overall, the best coefficient of determination was obtained from the reflectance ratio R415/R695 (r2 = 0.72) (Fig. 5c). This resulted because reflectance typically decreased at 415 nm and increased at 695 nm as Chl decreased under N stress. Stress-induced loss of Chl also was associated with changes in R415/R605 (r2 = 0.68) (Fig. 5c), due mostly to increased reflectance at 605 nm by N-deficient leaves. Due in part to the close association between Chl and carotenoids in the present study, canopy reflectance was best correlated with carotenoid pigments at R415/R685 (r2 = 0.79), R415/R625 (r2 = 0.78), and R415/R505 (r2 = 0.56) (data not presented).

Unlike leaf reflectance data, leaf N was best correlated with canopy reflectance wavebands centered in the blue region of the spectrum (maximum r2 = 0.31 at 410 nm), and to a lesser extent in the far-red region, at about 700 nm (Fig. 6a). Strong correlation was obtained between leaf N and canopy reflectance when R415 was divided by either R710 (r2 = 0.70) or R585 (r2 = 0.67) (Fig. 6b). A slightly lower coefficient of determination was obtained for R710/R405 (r2 = 0.63), which is nearly the inverse of R415/R710 (Fig. 6c).

Among the various canopy reflectance ratios identified from correlation analysis, the most sensitive indicator of N stress in cotton appeared to be R360/R710, R415/R585, R415/R710, and R415/R695 (Fig. 7) . Somewhat better relationships with leaf N were obtained with quadratic rather than simple linear functions for each of these single waveband ratios. A weak correlation was obtained between changes in leaf N and R755/R700 in canopies, but this relationship was significant (r2 = 0.82) when N values > 40 g kg-1 were excluded (Fig. 7f, dashed line). Although close associations were obtained between R755/R700 and Chl concentration in the present study (canopy: r2 = 0.56, Fig. 5; leaf: r2 = 0.51, Fig. 3), variability in R755/R700 at high Chl values suggests this canopy reflectance ratio may not be a sensitive indicator of N stress in cotton (Fig. 7f).



View larger version (31K):
[in this window]
[in a new window]
 
Fig. 7. Relationships between leaf nitrogen (N) concentration and (af) various single-waveband canopy reflectance ratios across different N treatments and sampling dates. Best-fit regression lines are represented by solid lines, with the exception of R755/R700 for which a dashed line is used to represent the regression function when leaf N values > 40 g kg-1 are excluded (df = 17, r2 = 0.82**). **Significant at P < 0.01; otherwise, not significant (ns).

 
In plants irrigated with 20% control N at first floral bud (square) stage onward, values for R415/R695 decreased 27% between 1 and 8 July (from about 0.28 to 0.22) as leaf N decreased 56% (from about 42.8 to 18.7 g kg-1) (Fig. 7e). Carter et al. (1996) reported values for R420/R695 were greater in control leaves, as compared with leaves stressed by various biological and physiochemical agents. Sensitivity to leaf N status by R415/R695 is very similar to that of the normalized difference pigment index (NPCI), an index formed by R430 and R680 (Jensen, 2000). Filella et al. (1995) reported NPCI measurements from canopies of field-grown wheat (Triticum aestivum L.) closely followed phenological and N treatment difference in leaf Chl. Canopy NPCI values in the present study were closely associated with changes in Chl (r2 = 0.65) and N (r2 = 0.49) concentrations, but somewhat larger correlation was obtained with R415/R695 (Fig. 5 and 6). This indicates that the best-fit reflectance ratios are specific for N deficit treatments in cotton. Changes in Chl or carotenoid concentration were more closely associated with canopy-level than with leaf-level NPCI. The light regime is expected to differ between full canopies and single leaves due to multiple light scattering, leaf layering, shading effects, and active photosynthesis–temperature effects (Jensen, 2000; Carter et al., 1996).

Canopy Reflectance Optimization with Multiple Linear Regression
The best three-waveband model for predicting leaf N from canopy spectral reflectance had an r2 value of 0.80, as compared with r2 of 0.92 for the best seven-waveband model (Table 1). The three-waveband model identified R437, R610, and R689 spectral regions within 25 nm of the canopy reflectance ratios closely correlated with leaf Chl (R605, Fig. 5) and N (either R415, R585, or approximately R700; Fig. 6). Most regression models shared R437 and R610. The best six- and seven-waveband models identified R755±5 (Table 1), a reflectance region closely associated with N status in plants under N-limiting conditions (Fig. 3, 4, and 7). None of the models identified R700±10, which was a significant spectral region in the present study and a measure commonly associated with changes in leaf Chl and photosynthesis in other hyperspectral reflectance studies (Chappelle et al., 1992; Carter et al., 1996). Strong colinearity of this red-edge reflectance feature with the other spectral variables may have decreased the contribution of R700 to the model (Gitelson and Merzlyak, 1996). Models with more than four variables indicated that R912, R890, and R852 were important in explaining observed variation in leaf N of cotton. Although NIR reflectance in canopies (R710<={lambda}<=755) was insensitive to changes in N status, its association improved when expressed in ratio to 415 nm (Fig. 7). In contrast, NIR measures in single leaves were very sensitive to changes in Chl and N under restricted N supply (Fig. 3 and 4). Unfortunately, these single-waveband regression models may not be easily transferred to other situations or crops in the field, because as mentioned previously, viewing geometry and canopy structure have a strong effect on overall characteristics of the reflectance spectrum (Jensen, 2000). Such effects on single waveband reflectance are typically modeled using radiative transfer approaches (Jacquemoud and Baret, 1990).


View this table:
[in this window]
[in a new window]
 
Table 1. Narrow-waveband (1.5-nm bandwidth) canopy reflectance predictors of leaf N in cotton grown under restricted N supply from first square onward and from first flower stage onward. Multiple linear regression procedures involved 66 values for N concentration (g kg-1) as the dependent variable and 60 single wavebands as the independent variable across all sampling dates and N-deficit treatments (see vertical arrows in Fig. 1).

 
Somewhat stronger associations with leaf N were obtained using multiple linear regression procedures, as compared with correlation procedures involving simple reflectance ratio (Table 1 and Fig. 7). Nevertheless, both procedures identified similar reflectance features: (i) a red-edge measure, (ii) a violet–blue reflectance region (approximately 420 nm), (iii) a yellow reflectance region (approximately 610 nm), and (iv) near-infrared reflectance (755–920 nm). In previous studies of cotton leaf reflectance, good precision in estimating foliar N was obtained from ratios using 460 or 495 nm (blue), and to a lesser extent 416 nm (violet), with wavebands centered at either 540, 516, 590, or 685 nm (Tarpley et al., 2000). In the present study, leaf N was closely associated with R437/R610 (Fig. 7), a canopy reflectance ratio formed by discrete wavebands consistently identified using regression procedures (Table 1); however, somewhat stronger associations were obtained from the canopy indices identified with simple correlation analysis (Fig. 7).


    CONCLUSIONS
 TOP
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 NOTES
 RESULTS AND DISCUSSION
 CONCLUSIONS
 REFERENCES
 
Hyperspectral reflectance was determined from leaves and canopies of cotton grown outdoors in large PVC pots, having a range of variability in N status commonly found in the field. Best-fit, linear relationships were determined between spectral reflectance indices and leaf N status. Results indicate that specific pigment and optical properties of photosynthetic tissues can be exploited to detect N deficiency at the canopy level. Providing Chl content is closely coupled to cotton nutrient status, specific features of cotton reflectance can provide an effective tool for managing crop N inputs, either on a site-by-site basis or prior to visible stress symptoms (Carter et al., 1996; Pan et al., 1997; Jones and Barnes, 2000). We expect that proper and timely application of nutrients through remote sensing assessment of crop N status may reduce application cost to farmers, and decrease leaching and runoff of mobile forms of fertilizer N (nitrite and nitrate).

Our analysis to discern spectral reflectance ratios sensitive to leaf N status in cotton support results of Tarpley et al. (2000). They found that leaf reflectance ratios between wavebands in the red edge (700–716 nm) and a waveband in the very near infrared region (755–920 nm) provided good prediction of leaf N concentration in cotton. A noticeable difference in the present study is the prevalence of changes in the 400- to 450-nm spectral range in canopies. The best correlation with Chl was obtained with R415/R695 in canopies. Multiple linear regression found that R437, R610, and R689 in canopies were important in explaining the observed N-induced and seasonal variation in leaf N concentration.

Somewhat larger correlation was obtained between plant N status and canopy reflectance than leaf reflectance. Further, the spectral signature(s) differed between the two sightings, with reflectance measurements at either 410 ± 5 or 700 ± 10 nm in canopies and in the 695- to 755-nm range in leaves being the most sensitive indicators of plant N status.


    ACKNOWLEDGMENTS
 
This study was in part supported by The National Aeronautical and Space Administration–funded Remote Sensing Technology Center at Mississippi State University (NASA Grant no. NCC13-99001). Technical help from Mr. Kim Gourley and Mr. Wendell Ladner is greatly appreciated.


    NOTES
 TOP
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 NOTES
 RESULTS AND DISCUSSION
 CONCLUSIONS
 REFERENCES
 
1 Mention of a trademark, proprietary product, or vendor does not constitute a guarantee or warranty of the product by the USDA and does not imply its approval to the exclusion of other products or vendors that also may be suitable. Back


    REFERENCES
 TOP
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 NOTES
 RESULTS AND DISCUSSION
 CONCLUSIONS
 REFERENCES
 




This article has been cited by other articles:


Home page
J. Environ. Qual.Home page
D. R. Tilley, M. Ahmed, J. H. Son, and H. Badrinarayanan
Hyperspectral Reflectance Response of Freshwater Macrophytes to Salinity in a Brackish Subtropical Marsh
J. Environ. Qual., April 5, 2007; 36(3): 780 - 789.
[Abstract] [Full Text] [PDF]


Home page
J. Environ. Qual.Home page
J. A. Poss, W. B. Russell, and C. M. Grieve
Estimating Yields of Salt- and Water-Stressed Forages with Remote Sensing in the Visible and Near Infrared
J. Environ. Qual., May 31, 2006; 35(4): 1060 - 1071.
[Abstract] [Full Text] [PDF]


Home page
Crop Sci.Home page
P. J. Starks, D. Zhao, W. A. Phillips, and S. W. Coleman
Development of Canopy Reflectance Algorithms for Real-Time Prediction of Bermudagrass Pasture Biomass and Nutritive Values
Crop Sci., February 24, 2006; 46(2): 927 - 934.
[Abstract] [Full Text] [PDF]


Home page
Agron. J.Home page
M-C. Belanger, A. A. Viau, G. Samson, and M. Chamberland
Determination of a Multivariate Indicator of Nitrogen Imbalance (MINI) in Potato Using Reflectance and Fluorescence Spectroscopy
Agron. J., October 19, 2005; 97(6): 1515 - 1523.
[Abstract] [Full Text] [PDF]


Home page
Agron. J.Home page
D. Zhao, K. R. Reddy, V. G. Kakani, J. J. Read, and S. Koti
Selection of Optimum Reflectance Ratios for Estimating Leaf Nitrogen and Chlorophyll Concentrations of Field-Grown Cotton
Agron. J., January 1, 2005; 97(1): 89 - 98.
[Abstract] [Full Text] [PDF]


This Article
Right arrow Abstract Freely available
Right arrow Figures Only
Right arrow Full Text (PDF)
Right arrow Alert me when this article is cited
Right arrow Alert me if a correction is posted
Services
Right arrow Similar articles in this journal
Right arrow Similar articles in ISI Web of Science
Right arrow