[go: up one dir, main page]
More Web Proxy on the site http://driver.im/
Next Article in Journal
Unmanned Aerial Vehicle Object Tracking by Correlation Filter with Adaptive Appearance Model
Previous Article in Journal
Research on a Seepage Monitoring Model of a High Core Rockfill Dam Based on Machine Learning
You seem to have javascript disabled. Please note that many of the page functionalities won't work as expected without javascript enabled.
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Reflectance-Based Vegetation Index Assessment of Four Plant Species Exposed to Lithium Chloride

Department of Environmental Engineering and Earth Sciences, Clemson University, Clemson, SC 29631, USA 2 Department of Statistics, Colorado State University, Fort Collins, CO 80523, USA
*
Author to whom correspondence should be addressed.
Sensors 2018, 18(9), 2750; https://doi.org/10.3390/s18092750
Submission received: 23 June 2018 / Revised: 2 August 2018 / Accepted: 17 August 2018 / Published: 21 August 2018
(This article belongs to the Section Remote Sensors)
Figure 1
<p>Relative reflectance of treatment to control for each of four species as acquired by CP.</p> ">
Figure 2
<p>Relative reflectance of treatment to control for each of four species as acquired by FOV.</p> ">
Figure 3
<p>Weeks 1, 2, and 3 of the first <span class="html-italic">A. thaliana</span> experiment; response was similar in the second experiment. (<b>a</b>) Week 1 control plants; (<b>b</b>) Week 1 treatment plants which exhibit slight chlorosis along the edges of older leaves; (<b>c</b>) Week 2 control plants; (<b>d</b>) Week 2 treatment plants which exhibit chlorosis, necrosis, and decreased biomass; (<b>e</b>) Week 3 control plants; (<b>f</b>) Week 3 plants which exhibit significant necrosis.</p> ">
Figure 4
<p>Example (<b>a</b>) <span class="html-italic">B. napus</span> and (<b>b</b>) <span class="html-italic">H. annuus</span> plants at week 3. Arrows in (<b>a</b>) are pointing out areas of chlorosis around leaf edges. The arrow in (<b>b</b>) is pointing out a necrotic spot.</p> ">
Figure 5
<p>Chlorophyll content (Chl <span class="html-italic">a + b</span>, µg/mL) by week and treatment level for each species considered. Lines (solid = control, dashed = treatment) indicate mean values at each week.</p> ">
Figure 6
<p>Relative water content (RWC) by week and treatment level for each species considered. Lines (solid = control, dashed = treatment) indicate mean values at each week.</p> ">
Figure 7
<p>Dry biomass (mg) by week and treatment level for each species considered. Lines (solid = control, dashed = treatment) indicate mean values at each week.</p> ">
Figure 8
<p>CLAI by week and treatment level for each species considered. Lines (solid = control, dashed = treatment) indicate mean values at each week.</p> ">
Figure 9
<p>Height by week and treatment level for <span class="html-italic">H. annuus</span>. Lines (solid = control, dashed = treatment) indicate mean values at each week.</p> ">
Figure A1
<p>Control and treatment spectra for <span class="html-italic">A. thaliana</span> as acquired by CP and FOV.</p> ">
Figure A2
<p>Control and treatment spectra for <span class="html-italic">Z. mays</span> as acquired by CP and FOV.</p> ">
Figure A3
<p>Control and treatment spectra for <span class="html-italic">H. annuus</span> as acquired by CP and FOV.</p> ">
Figure A4
<p>Control and treatment spectra for <span class="html-italic">B. napus</span> as acquired by CP and FOV.</p> ">
Versions Notes

Abstract

:
This study considers whether a relationship exists between response to lithium (Li) exposure and select vegetation indices (VI) determined from reflectance spectra in each of four plant species: Arabidopsis thaliana, Helianthus annuus (sunflower), Brassica napus (rape), and Zea mays (corn). Reflectance spectra were collected every week for three weeks using an ASD FieldSpec Pro spectroradiometer with both a contact probe (CP) and a field of view probe (FOV) for plants treated twice weekly in a laboratory setting with 0 mM (control) or 15 mM of lithium chloride (LiCl) solution. Plants were harvested each week after spectra collection for determination of relevant physical endpoints such as relative water content and chlorophyll content. Mixed effects analyses were conducted on selected endpoints and vegetation indices (VI) to determine the significance of the effects of treatment level and length of treatment as well as to determine which VI would be appropriate predictors of treatment-dependent endpoints. Of the species considered, A. thaliana exhibited the most significant effects and corresponding shifts in reflectance spectra. Depending on the species and endpoint, the most relevant VIs in this study were NDVI, PSND, YI, R1676/R1933, R750/R550, and R950/R750.

1. Introduction

Anthropogenic activities can result in the release of a wide array of contaminants, particularly metals, to the environment. Responsible environmental stewardship involves the management and remediation of such releases. Numerous remediation strategies exist, depending on the circumstance, but the technique considered here is reflectance spectroscopy. Reflectance spectroscopy has potential for use as a cost-effective, non-destructive analytical technique for detecting and assessing plant stress, specifically metal stress [1,2,3].
Chemometric mathematical methods are often used to analyze reflectance spectra when extensive samples are available; a few hundred are required to develop a robust model. For smaller sample sizes, as in this study, vegetation indices (VI; mathematical combinations of different reflectance spectral bands) can be used to provide rapid and convenient semi-analytical measures of vegetation activity, which in turn can provide indication of plant health [4].
The ability to detect metal exposure in plants is particularly relevant when employing phytoremediation strategies. Phytoremediation is defined, for our purposes, as the use of green plants for environmental clean-up, or the use of plants to remove or neutralize pollutants in the biosphere [5]. Specific applications for reflectance spectroscopy as a tool complimentary to phytoremediation could include: (1) early recognition of contaminated areas through identification of plant stress indicative of metal exposure, (2) surveillance of sites with existing contamination or sites with the potential to become contaminated, (3) assessment of phytoremediation efforts, (4) assessment of risk to human health and the environment through coarse quantification of contamination, etc.
The contaminant of concern in this study is lithium. Lithium (Li) is widely used in the USA, which is the leading producer and consumer of Li materials, finding utility in ceramics and glass, aluminum production, the medical industry, certain batteries and greases, nuclear reactor coolant, radiation dosimeters, and historically, in nuclear weapon development [6,7,8,9]. Although Li is not a radioactive concern, it is an anthropogenic contaminant related to the nuclear fuel cycle and to legacy waste and contamination from nuclear weapons development [8,10,11]. For example, historical waste-disposal activities at the Department of Energy’s Oak Ridge Y-12 Plant resulted in the release of Li to the groundwater [8]. Numerous studies have shown that there are shifts in plant reflectance spectra due to metal stress (or simulated metal stress) [2,12,13,14,15,16,17,18,19,20,21,22,23,24] but to the authors’ knowledge none have considered response to Li exposure.
To further investigate the use of reflectance spectroscopy as a useful method for assessing metal stress in plants, reflectance spectra for four species of plants were collected every week for three weeks using an ASD FieldSpec Pro spectroradiometer with both a contact probe and a field of view (FOV) probe. These plants were treated twice weekly in a laboratory setting with 0 mM (control) or 15 mM of lithium chloride (LiCl) solution and harvested weekly, immediately after spectra collection, for assessment of physical endpoints such as relative water content and chlorophyll content. The specific objectives of this exploratory study are to: (1) identify changes in plant status due to Li exposure in multiple plant species and (2) determine if reflectance is useful in identifying these changes through the utilization of both previously defined and newly determined vegetation indices.

1.1. Plant Species Considered

The species considered were Arabidopsis thaliana, Helianthus annuus (sunflower), Brassica napus (rape), and Zea mays (corn). A. thaliana is a member of the mustard family that is closely related to various crop plants. It has been the subject of intense study over the past several decades and is considered a model organism and ideal for use in the laboratory setting for biological research [25]. Several species of the Brassica family, which are vegetable and oilseed crops, have been identified as metal accumulators and are considered potential phytoremediation candidates [26,27,28]. H. annuus is an ornamental flower as well as an important environmental crop, and it has been shown to be an effective phytoremediation crop [29,30,31]. Z. mays has also been shown to have metal phytoremediation potential [31]. It is the major feed grain (>95%) in the United States and is also processed into a broad assortment of food stuffs, from cereals to sweeteners. Additionally, Z. mays has industrial utility as a component in the fabrication of fuel ethanol. The United States is currently the world’s largest producer and exporter of corn [32]. Several studies have considered the reflectance spectra of corn and sunflower, from assessing pigment concentrations to nutrient/water status and photosynthetic efficiency at both leaf and canopy scales [33,34,35,36,37,38,39,40,41,42,43,44,45].

1.2. Lithium in Plants

Lithium is the lightest metal, although as it is highly reactive, does not occur naturally in its elemental form. It occurs in various minerals and salts, and typically enters the environment through weathering processes. Once in the environment, Li is easily transported to above-ground plant biomass as it shares the potassium (K+) transport carrier; all plants will take up Li [7,8,9].
Lithium uptake, tolerance, and toxicity symptoms are all species specific, although stimulatory effects are commonly seen at low levels of Li [7,9,46,47]. Although symptoms of lithium toxicity are not distinct between species, general symptoms include chlorosis, necrotic spot development, leaf curling, and reduced biomass, all similar symptoms as exhibited in pathogen defense [48,49]. There is preferential uptake of Li in older leaves, so effects will be less pronounced on new growth [49,50,51]. LiCl has been used previously to consider Li uptake and stress, but although chlorine (Cl) is an essential micronutrient for higher plants, at high plant tissue concentrations Cl as the chloride ion (Cl) can be toxic [7,26,52,53]. However, the concentrations of LiCl used here are considered below levels of which Cl might be toxic; contribution of Cl to the effects in this case are considered negligible [52,53].
Several studies have considered Li toxicity in plants [7,26,47,48,49,50,51,54,55,56,57,58] although the mechanism is not fully understood; plants have each type of enzyme shown in animals and yeast, respectively, to be Li sensitive, namely inositol monophosphatase and HAL2 nucleotidase. It is fairly unknown which of these two enzymes is the major target of Li action [54]. Generally, at high concentrations, Li increases the production of ethylene, which is known to inhibit plant growth. The mechanism is not wholly understood, although the “inositol depletion hypothesis” is generally accepted. This theory holds that Li+ inhibits inositol monophosphatase, which ultimately triggers ACC (aminocyclopropane carboxylic acid) synthase, resulting in an increase in ethylene biosynthesis [7,49,54,55]. There is also evidence to suggest that Li toxicity can phenotypically resemble magnesium (Mg) deficiency and that Li+ can bind to proteins (such as chlorophylls) that normally bind with Mg2+ [56].

2. Materials and Methods

2.1. Plant Growth and Treatment

The soil mix used was a 4:1 mixture of PGX (Promix PGX, Premier Horticulture Inc., Quakertown, PA, USA) and perlite (Hoffman Horticultural Perlite, Good Earth, Lancaster, NY, USA). Soil was mixed and placed in square plastic grow pots (10.8 × 10.8 × 12.7 cm) with perforated bottoms to allow water seepage; soil was well-hydrated prior to sowing seeds. Nutrient solution used to water and treat plants was made with DI water, 1/32 strength Murashige and Skoog basal medium (137.5 mg L−1) (cat. no. M5519, Sigma-Aldrich, St. Louis, MO, USA), and 250 mg L−1 MES hydrate (cat. no. M2933, Sigma-Aldrich), using KOH to pH balance to 5.7. A. thaliana seeds (WT-02-41-01 Columbia [alias Col-0] Wildtype, LEHLE Seeds, Round Rock, TX, USA) were soaked in nutrient solution, exposed to red light for 30 min to synchronize germination, and pipetted onto the soil pots. Z. mays seeds (Burpee Sweet Corn Bi-Licious Hybrid, Burpee Garden Products Co, Warminster, PA, USA) H. annuus seeds (Snow Country Black Oil Sunflowers, Ridley Inc., Mankato, MN, USA), and B. napus seeds (Winfred Brassica Rape, outsidepride.com, lot: M31-9-2WIN, Independence, OR, USA) were planted in the soil pots by hand. Following the sowing of the seeds, the 1/32 nutrient solution was further diluted to 1/64 strength for subsequent treatments of all plant species.
After planting, arbitrary sets of 6 pots each were transferred to plastic tubs (40 × 31.75 × 15.24 cm) in 3 cm distilled water. Tubs were placed in rows of up to four on growth shelves, 42 cm beneath growth lights. Plants were on a 9 h light: 15 h dark cycle under ambient laboratory environmental conditions. At the seedling stage, plants were culled and/or redistributed to ensure an appropriate number of plants in each pot, based on the appearance of health: A. thaliana were culled to three seedlings per pot, B. napus were culled to 10 plants per pot, Z. mays were culled to three plants per pot, and H. annuus were culled to one plant per pot. Immediately prior to LiCl treatment, pots were randomly rearranged between tubs (six pots per tub, no longer in DI water) such that each tub, now serving as a treatment group, had similar size and quality plants. Each experiment had three treatment group tubs and three control tubs for a total of 36 samples (i.e., plants) per experiment. Each experiment was conducted twice such that there was a total of 72 plant samples per species, or 288 total samples.
Once pots were arranged for treatment, spike solution was evenly applied to the top of each pot as 100 mL (25 mL delivered to each quadrant) of 15 mM LiCl in 1/64 strength nutrient solution twice weekly, with control plants receiving 100 mL 1/64 nutrient solution only. Two pots were randomly selected from each treatment group for weekly spectra collection and harvest. After each application of nutrient solution, the plants were rotated within the tubs and the tubs were rearranged among the growth shelves to account for potential variation in lighting or other environmental conditions.

2.2. Equipment Setup and Spectra Collection

Reflectance spectra were collected using a FieldSpec Pro (FSP 350-2500P; Analytical Spectral Devices (ASD), Boulder, CO, USA) which is a full range (350 nm to 2500 nm) portable spectroradiometer (with sampling intervals/spectral resolutions of 1.4 nm/3 nm and 2 nm/10 nm for 350 to 1000 nm and 1000 to 2500 nm respectively) [59]. Contact probe (CP) spectra were collected using a leaf clip attachment on individual leaves. The CP provides light (3.825 V, 4.05 W low intensity bulb) and collects reflectance spectra. The leaf clip attachment has both a white (for white reference) and black (to minimize back scatter) background. Multiple reflectance readings were taken on the leaves of each species of plant to obtain an overall representation of reflectance.
FOV spectra were collected using an 8° probe (i.e., a viewing angle of 8°). Incident light was provided by two halogen lamps (Pro Lamp, 14.5 V, 50 W, P/N 145378, ASD) angled at 30 degrees from horizontal. The lights were 180° apart at 30.5 cm from the center of pot on the horizontal and 76.2 cm above the table surface. The fore optics probe was centered between the lights at 66.7 cm above the plane of the pot surface, resulting in a spot size diameter of 9.32 cm. Reflective surfaces were covered with light-absorbent material to minimize noise and thus variability in spectra, and dark room conditions were approximated by surrounding the lights and fore optics with a black felt canopy. Tripod surfaces were also wrapped in black felt and the table surface was lined with light-absorbent black rubber. The white reference was a calibrated Spectralon (25.4 × 25.4 cm, LabSphere, North Sutton, NH, USA) panel of 99% reflectance that was elevated to a height equivalent to a grow pot. Four spectra, each collected at a different arbitrary rotation of the pot, were acquired and then averaged to get an overall assessment of the reflectance of the sample. FOV spectra were always acquired prior to CP because it is possible for the CP to injure the plant and therefore affect subsequent FOV readings.

2.3. Collection of Physical Measures

As metal stress is known to mimic drought stress [3], plants were harvested after spectra collection each week to determine chlorophyll content and relative water content. The concentrations of chlorophyll a (Chl a) and chlorophyll b (Chl b) were determined for each replicate (i.e., pot) [26,60,61]. Four circular leaf subsamples were collected from representative leaves of the plants in a pot using a #3 cork borer (Fisher Scientific, Pittsburgh, PA, USA). Leaf samples were stored in the dark at 4 °C in capped 20 mL vials (KG-33 borosilicate glass; Kimble Chase, Vineland, NJ, USA) containing 2 mL 100% ethanol for three days before absorbance (A) at 665 nm, 649 nm, 629 nm, and 696 nm, with an offset at 750 nm, was determined for 1.5 mL subsamples for each vial using a NanoDrop 2000c UV-Vis spectrophotometer (Thermo Scientific, Wilmington, DE, USA). Disposable methacrylate cuvettes with transmission from 300 to 800 nm >80% were used with the 1.5 mL subsample (Cole Palmer, Vernon Hills, IL, USA). Chlorophyll content was determined using appropriate, previously published [61] equations where A x is absorbance at x nm:
Chl   a   ( μ g / mL ) = 5.2007 A 649 + 13.5275 A 665 ,
Chl   b   ( μ g / mL ) = 22.4327 A 649 7.0741 A 665 .
After leaves were sampled for chlorophyll content, additional sufficient leaves were removed to obtain between 1000 and 2000 mg of fresh mass for each replicate (i.e., pot) to determine relative water content. Samples were placed in weigh boats, fresh mass was obtained, samples were dried in an oven to a constant mass, and dry mass was obtained. A sample’s relative water content (RWC) was then expressed as Equation (3):
RWC = 1 dry   mass   fresh   mass
Lithium uptake has been shown to have a stimulatory effect at low levels, and potentially cause necrosis at high levels depending on the plant species [47,48], so remaining plant material from above was similarly retained, dried, and weighed to determine each sample’s overall dry biomass. Note that dry biomass was determined for all experiments except the first A. thaliana experiment, as these samples were initially used elsewhere for teaching purposes.
A visual assessment of the proportion of a pot covered by plant material was performed by overlaying a 6 × 6 (18 mm × 18 mm) grid on top-down photos of each group of plants at each of three timepoints, forming 36 squares with 49 evenly-spaced points (grid intersections). Photographs were taken immediately prior to spectra collection, directly above each six-pot treatment group and six-pot control group in the same manner each week. However, to account for any potential change in magnification or alignment, gridlines were laid based on pot dimensions, which were definitively consistent. Using the grid intersections, an additional endpoint, Coarse Leaf Area Index (CLAI), was defined as:
CLAI = N leaf   N total
where Ntotal is the total number of points in the grid and Nleaf is the number of points on leaf material. CLAI provides an approximate indication of how much of the pot surface is covered by plant material. Finally, heights of plants above the top of the pot were measured. The only plants with noticeable variation in heights were H. annuus and thus are the only results reported for plant height.

2.4. Data Analysis

Initially, fifteen VI were considered for applicable spectra acquisition technique(s) (i.e., FOV and/or CP), including indices from the literature [e.g., normalized difference vegetation index (NDVI), photochemical reflective index (PRI), normalized phaeophytinization index (NPQI), structural independent pigment index (SIPI), pigment specific normalized difference (PSND), normalized difference index (NDI), water index (WI), yellowness index (YI)] [15,19,39,62,63,64,65,66,67,68,69,70,71]. Six VI (R1390/R1454, PRI, SIPI, NDI, R1110/R810, and R725/R675) were removed from consideration as they were highly colinear with other VI. Pearson correlation coefficients (PCC) were utilized to determine collinearity, and VI were considered for exclusion if the absolute value of the PCC between two (or more) was greater than 0.75 for CP and 0.50 for FOV. Preference was given to VI previously associated with metal exposure.
For each species at each view, a mixed-effects model analysis was conducted for each vegetation index and each plant endpoint (RWC, Chl a + b, dry biomass, CLAI, and/or height) to consider the fixed effects of week (1, 2, 3), treatment group (control or 15 mM LiCl treatments), and week-by-treatment interaction with a random effect for sampling unit. The Kenward-Rogers approximation of degrees of freedom was used to account for variation among the week-by-treatment combinations. When a treatment effect or week-by-treatment interaction was significant for a plant stress endpoint (RWC, Chl a + b, biomass, CLAI, and/or height), subsequent linear mixed-effects analyses were conducted to consider the endpoint measure as a dependent variable with vegetation indices as predictors that also had significant treatment or week-by-treatment interaction effects from the mixed-effects model analysis as predictors. SAS Software v. 9.4 (SAS Institute Inc., Cary, NC, USA) was used for all statistical analyses and a significance level of 0.05 was used for all tests of significance.

3. Results

3.1. Spectra

Average relative reflectance spectra of treatment compared to control are shown by week and technique (CP or FOV) for each species in Figure 1 (CP) and Figure 2 (FOV). Average reflectance spectra for each experimental group (i.e., control and treatment) are shown in Appendix A (Figure A1, Figure A2, Figure A3 and Figure A4). Average reflectance data (from which the figures were developed) is included in the supplemental online material. Notice that the vertical axes in Figure 1 and Figure 2 are consistent between Z. mays, H. annuus, and B. napus, but needed an extended range for A. thaliana.

3.2. Phenotypic Observations

A. thaliana had visible symptoms of Li exposure starting in week 1, from a slight yellowing (chlorosis) at the leaf tips in week 1 to significant necrosis in week 3 (Figure 3).
B. napus plants showed some symptoms of toxicity around leaf edges in weeks 2 and 3 (Figure 4a). A few H. annuus plants began exhibiting slight symptoms of lithium toxicity between weeks 2 and 3, as slightly mottled leaves and occasional necrotic spots (Figure 4b). Z. mays showed no symptoms of toxicity. Observed symptoms are typical of lithium toxicity; different species of plants are known to have varying tolerances to lithium exposure [48,72].

3.3. Endpoints

Plots of chlorophyll content, relative water content, dry biomass, and CLAI by treatment level and time are shown in Figure 5, Figure 6, Figure 7 and Figure 8 respectively for each species. A plot of H. annuus height is shown in Figure 9. Summary statistics for these endpoints are provided in Appendix B. Notice in Figure 6 (RWC) that the vertical axis is different for A. thaliana compared to the other species as the treatment plants had a substantially greater range of values. Other figures depicting endpoint results have consistent scales on the vertical axes.
Results (p-values) of the statistical analysis for relative water content, chlorophyll content, dry biomass, CLAI, and height are shown in Table 1. F(n,d) in table headers indicates the degrees of freedom with which the F test statistic and associated p-value were calculated, where n is numerator degrees of freedom, and d denominator degrees of freedom. Note that denominator degrees of freedom vary due to differing numbers of observations and Kenward-Rogers degrees of freedom approximation. There are no significant time, treatment, or time by treatment interaction effects for B. napus chlorophyll content, nor for Z. mays CLAI. All other endpoints had significant differences in time for each species. There was a significant treatment effect on Chl a + b and a significant time by treatment interaction effect on height for H. annuus.

3.4. Vegetation Indices and Overall Response

Detailed results (p-values) from the initial mixed-effects analysis of the vegetation indices are contained in Appendix C. Significant p-values (p < 0.05; bolded) indicate a statistically significant difference between the treatment, week, or treatment-by-week interaction outcome means. As above, denominator degrees of freedom vary. From these values we see that methods of spectra acquisition as well as the results between species are generally varied. However, this is not wholly unexpected as four fairly different species were considered; these species vary in leaf size, height, structure, etc.
Nearly all A. thaliana VIs had significant treatment, time, and interaction effects between treatment and time, for both spectra acquisition techniques (CP and FOV). For other species, however, results differ between VI, spectra acquisition techniques, endpoints, and time/treatment. Every species exhibited significant time-dependent differences in VI. Corn only had significant treatment effects for YI acquired by FOV. H. annuus had significant treatment or interaction differences in PSND, YI, R750/R550, and R1636/R1933. B. napus had significant treatment or interaction differences in WI, YI, R750/R550, and R1636/1933. However, only A. thaliana and H. annuus had significant treatment effects with respect to measured endpoints (RWC, Chl a + b, biomass, CLAI, and/or height) (Table 1), so only these species were considered in secondary (i.e., follow-up) analysis to determine which VIs were significant predictors of these endpoints. Results of the corresponding mixed effects model are shown in Table 2 below. The only index with no predictive power was WI; all other VI were significant predictors of at least two endpoints. Interestingly, R750/R550 was a significant predictor in at least one view (CP or FOV) of all endpoints considered, as further discussed below. This index is therefore deemed the most promising of those considered for prediction of lithium exposure.

4. Discussion

In this study, four species of plants were treated with lithium chloride over three weeks to discern if symptoms of the exposure could be identified from their reflectance spectra through the use of VI. Plants were harvested each week, with reflectance spectra collected at the whole plant and leaf scales. Various physical endpoints were measured immediately after spectra acquisition. Controls were maintained in parallel with the treatment plants at the same sample size, as matching controls are important to be able to discern if temporal changes in reflectance result from contaminant exposure or from the natural growth pattern or inherent biological variability of the plants [17,22,69].
The progressive Li exposure had slightly different effects, as well as severity of effects depending on the species. All species considered exhibited significant changes in time, but only physical measures of H. annuus and A. thaliana exhibited significant treatment effects. Although Z. mays did not demonstrate any obvious symptoms of Li toxicity, B. napus did show slight symptoms of toxicity around the leaf edges starting in week 1. Uptake of Li is expected to be higher in dicots (e.g., A. thaliana, H. annuus, and B. napus) than in monocots (e.g., Z. mays) [7], which might explain why no symptoms were seen in Z. mays. Similarly, the low biomass of A. thaliana as compared to plants with much higher biomass for an equivalent soil concentration may have led to the more severe toxicity symptoms.
Observed symptoms in A. thaliana, H. annuus, and B. napus were consistent with Li toxicity, which in turn are similar to those resulting from Mg deficiency [56]. Mg is essential for photosynthesis as it is the central element of the chlorophyll molecule [73]. It is possible that the competitive binding of Li results in Mg deficiency, exacerbating the toxicity symptoms. Of the initial 15 VI considered, seven were found to be significant predictors of these treatment-dependent endpoints. Comparisons between species and the associated relationships with the acquired spectra and calculated VIs are discussed below.

4.1. Chlorophyll Content

Excess metal exposure negatively affects photosynthetic processes and typically induces a general stress reaction in plants [74]. Photosynthetic pigments typically decrease with metal exposure, which has obvious consequences for photosynthesis and plant growth. Inhibition of photosynthesis is one effect that most metals have in common when present at toxic concentrations; reduction in photosynthetic efficiency will be seen as an increase in reflectance in the visible range, as less light is being utilized for photosynthesis and chlorophyll production [3].
A. thaliana and H. annuus showed significant changes in chlorophyll content by treatment level, which corresponds to the increase in reflectance in the visible region (see Figure 1 and Figure 2; notice differences in H. annuus are subtle compared to A. thaliana). Interestingly, for H. annuus, chlorophyll content increased from week 1 to week 2, but remained essentially constant from week 2 to week 3, whereas Z. mays chlorophyll content remained the same from week 1 to week 2 and decreased from week 2 to week 3 (Figure 5). Also, although not statistically significant, chlorophyll content in treatment plants appears to increase slightly above that in control plants in Z. mays in week 3 (Figure 5), which corresponds to the increase in absorbance (and thus decrease in reflectance) seen in the Z. mays spectra in Figure 1 and Figure 2. There were no significant differences in B. napus chlorophyll content by week or treatment (Table 1), and B. napus also generally had lower chlorophyll content than the other species (Figure 5 and Table A1); Z. mays and H. annuus had the highest chlorophyll contents.
The vegetation indices that proved to be significant predictors of chlorophyll content for H. annuus and A. thaliana were PSND, YI, R750/R550, and R1636/R1933. PSND, which is determined by the normalized difference (R800 − R680)/(R800 + R680), has been shown to be correlated with pigment concentration per unit area [70]. Here, PSND was found to be a significant predictor of chlorophyll content only for spectra acquired by CP, not FOV (Table 2). However, for A. thaliana, PSND as determined using FOV was a significant predictor of CLAI. Taken together, this is consistent with previous findings. Note that CLAI did not have significant treatment differences for H. annuus, and PSND also did not have significant treatment differences for H. annuus spectra acquired via FOV. Similarly, R1636/R1933, which has been previously suggested to be associated with cesium (Cs, also an alkali metal) exposure in A. thaliana [69], was also a significant predictor of chlorophyll content when assessed using CP spectra, and CLAI when assessed utilizing FOV spectra.
YI is an approximated second derivative of the spectra at about 600 nm and is intended to provide indication of chlorosis, i.e., leaf yellowing [67]. In this study, YI was found to be a significant predictor of chlorophyll content with FOV spectra for H. annuus and with CP spectra for A. thaliana. As the YI was originally developed at the CP scale, it is likely that soil background will influence this VI, particularly in the instances like that of A. thaliana week 3 treatment plants where much of the soil is exposed to the fore optics. H. annuus leaves are large and broad and, being closer to the fore optics than A. thaliana plants (and much closer than the soil), leaf reflectance dominates the spectra. Also, FOV may potentially capture the mottled nature of the yellowing of H. annuus leaves better than CP. R750/R550 was originally developed for remote sensing of chlorophyll [68,71] but has also been found to be correlated with metal content [15]. Here, R750/R550 was found to be a significant predictor of chlorophyll content at both CP and FOV scales for both species considered.

4.2. Water Content

Metals can disrupt the plant-water balance [3], which can be seen as an increase in reflectance in the mid-infrared region, as water absorbs fairly strongly at 1450 nm, 1940 nm, and 2500 nm, with slight increase in the near infrared region, associated with slight water absorption [75,76]. A. thaliana was the only species that exhibited significant treatment effects associated with RWC (Table 1), although all species exhibited significant temporal differences, with RWC decreasing in time (Figure 6). The difference in treatment versus control reflectance for A. thaliana is clear in Figure 1 and Figure 2 as spikes in relative reflectance at the aforementioned water absorption wavelengths that get larger each week. Although difficult to see in the relative spectra, temporal differences in the remaining species can be seen in the average (i.e., raw) reflectance spectra as an increase in reflectance at these same wavelengths in both treatment and control plants (Figure A2, Figure A3 and Figure A4). Differences for each species except corn are greater between weeks 2 and 3 than between weeks 1 and 2, although the overall range of RWC is narrow for control plants, varying by only a few percent. It is interesting to note that the 1940 nm band in Figure 2 appears to give a qualitatively more consistent relative response between weeks than the other water bands. Although not statistically significant, B. napus RWC of treatment plants appears to decrease in time more than the control plants. This species, along with H. annuus had increasing variability (i.e., standard deviation) in RWC with time (Table A1), which explains why the apparent difference is not statistically significant.
The vegetation indices that proved to be significant predictors of RWC for A. thaliana were R750/R550 with FOV, R950/R750 with CP, and (R950 − R750)/(R950 + R750) with CP. The ratio (R950 − R750)/(R950 + R750) was considered because normalized differences are frequently utilized to improve upon a simple ratio by accounting for background reflectance [71]. In this instance, the normalized ratio has similar significance to its corresponding simple ratio, R950/R750, providing good indication that background reflectance is not a confounding factor. Similar to R1636/R1933 discussed above, R950/R750 was previously found to be associated with Cs exposure in A. thaliana, being correlated with RWC and Chl a + b at the CP scale and CLAI at the FOV scale in Cs contaminated plants [69]. Here, R950/R750 was not found to be significant for Chl a + b or CLAI (although p = 0.0661 [CP] and p = 0.0663 [FOV], which would be significant at a 90% confidence level), but it is found to be significant for RWC with CP. Also consistent with this study is that the Cs study also found R750/R550 at the FOV scale to be a significant predictor of RWC. This suggests that, as mentioned above, R750/R550 along with R950/R750, may be good predictors of plants exposed to alkali metals.

4.3. Remaining Endpoints: Size and Shape

Although there are no significant treatment effects, Z. mays, B. napus, and H. annuus share similar magnitudes of and temporal changes in dry biomass (Figure 7). By week 3, plants have more than doubled their week 1 biomass, and B. napus has the highest average mass followed by Z. mays, then H. annuus, although Z. mays has the largest variation in biomass; the maximum Z. mays biomass was about 1700 mg, followed by a maximum of about 1200 mg for both H. annuus and B. napus. The temporal increase in FOV acquired reflectance for these species (for both treatment and control groups) in the near-IR is partially due to the increase in biomass of the plants (Figure A2, Figure A3 and Figure A4). For H. annuus, not only do the plants get taller (Figure 7), and therefore become closer to the fore optics, but there is also an increase in number and size of leaves. H. annuus leaves are broad and grow parallel to horizontal, resulting in a large reflective surface perpendicular to the fore optics. There will be an increase in light scattering within a plant containing a greater proportion of cell surfaces exposed to intercellular air space, due to different indices of refraction of these materials. As near infrared light is not used for photosynthesis, increasing light scatter in the infrared region means less transmission of light through the plant and more reflection back to the fore optics. Differences in plant structure will therefore affect reflectance of light in the near-IR; larger leaf areas will result in higher reflectance, whereas cell degradation or reduction in leaf thickness will result in lower reflectance, as there will be less light scattering within the plant leaf [77].
Although the mean CLAI for H. annuus is comparable to that for B. napus (~50%–60% coverage), the variability is much greater with some plants at almost 100% coverage by week 3. Also, B. napus grew less than 10 cm above the top of the pot so are much further away from the fore optics than H. annuus. Although there were no differences in Z. mays heights between groups (details not shown), plants grew to ~45 cm in week 3, about double the H. annuus control plant mean (Table A2). However, the leaves are narrow and grow at an angle vice horizontal, meaning there is less leaf surface directly exposed to the fore optics compared to H. annuus. In fact, there are no significant differences in time or treatment for Z. mays CLAI, which remains essentially constant at ~35% coverage. These points together can be seen as differences between species in the magnitude of reflectance in the near-IR of FOV spectra in Figure A2, Figure A3 and Figure A4. For Z. mays, reflectance in the near-IR varies between ~0.18 and 0.25, for B. napus between ~0.32 to 0.50, and for H. annuus between ~0.48 to 0.80 from weeks 1 to 3.
Contrasting with the other three species, A. thaliana had significant treatment and temporal differences in biomass. A. thaliana plants also had much less biomass than the other species considered, as they are small, squatty plants prior to flowering. The dry biomass maximum for A. thaliana was less than 300 mg. Field-of-view control reflectance for A. thaliana increased in the near-IR each week (Figure A1), corresponding to the increase in biomass of the control plants. Treatment reflectance in the near-IR was lower than the control reflectance for both the CP and FOV in weeks 2 and 3 (Figure 1 and Figure 2); the lower CP reflectance indicates a difference in leaf structure. Lower FOV reflectance indicates a lower biomass comparatively, and a shift in the shape also indicates a structural difference. A. thaliana also has significant treatment differences for CLAI, with the control plant mean exceeding the H. annuus mean CLAI. However, prior to flowering A. thaliana does not grow vertical as H. annuus and Z. mays, so reflectance in the near-IR of ~0.3–0.4 in the FOV spectra relates more to biomass.
All species showed slightly increased reflectance in the FOV treatment group as compared to the control in week 1 (Figure 2). This implies a possible initial stimulatory effect at these time points, which is commonly seen response to low-level lithium exposure [7,47]. However, it also appears that there may be structural degradation at week 2 for all four species as treatment reflectance in FOV spectra is lower than control reflectance (Figure 2). In week 3 these differences remain in A. thaliana and Z. mays spectra but are no longer apparent in H. annuus or B. napus spectra. It has been suggested that plants could simultaneously experience both stimulatory effects and toxicity symptoms from Li exposure, due to different concentrations of Li within the plant [50,51]. A separate study saw stimulated growth along with slight chlorosis in leaves of snap bean (Phaseolus vulgaris) exposed to 4 ppm lithium nitrate (LiNO3) [47].
The vegetation indices that proved to be significant predictors of biomass for A. thaliana were YI with FOV and R750/R550 with both CP and FOV. Significant predictors of CLAI were NDVI, PSND, R950/R750, R1636/R1933, and (R950 − R750)/(R950 + R750) with FOV, and R750/R550 with both CP and FOV. NDVI has been well-associated with green biomass and leaf area [63]; here NDVI is not associated with A. thaliana biomass, but treatment-induced changes in biomass were also associated with necrosis and browning of the leaves. Biomass and CLAI are plant scale endpoints as opposed to leaf scale endpoints, however, R750/R550 was also a significant predictor of biomass at the leaf scale. One explanation for this is that biomass is related to leaf thickness (i.e., thinner leaves would result in less biomass), which would be seen in CP acquisition. Also, as statistically significant treatment effects were observed, A. thaliana biomass is related to lithium exposure. It is therefore highly likely that R750/R550 is responding to lithium uptake, particularly as it is a significant predictor of all other treatment-dependent endpoints considered. The remaining VI have been previously discussed.
The significant predictors of H. annuus height were PSND, YI, and R750/R550, which are the same indices as were significant for chlorophyll content in H. annuus (Table 2). YI was significant with FOV spectra; as H. annuus leaves mottle slightly with time, the leaves also become broader and closer to the fore optics which could also be an increase in yellowness. PSND was significant with CP spectra, which is at first counterintuitive. However, as the plant grows taller, the internal structure of the plant changes as well. As discussed above, PSND has been associated with areal pigment concentration which likely changes as the plant grows. R750/R550 was significant at both the leaf and plant scales. H. annuus height has significant interaction effects, which essentially indicates that height differences depend on treatment duration. Lithium uptake into plants follows the dose given, although the specific relationship is species dependent [7]. Thus, height is likely associated (negatively) with lithium uptake. R750/R550, as above, is likely more indicative of lithium uptake than height specifically, as also discussed above for other endpoints.

5. Conclusions

This study explored the feasibility of using reflectance spectroscopy as a quick and easy-to-use technique in supplementing phytoremediation efforts. Certain vegetation indices seem promising for selected endpoints and species but the variable responses of plants to similar Li concentrations makes applying VI across species less reliable. Treating species with various concentrations of Li to induce a similar level of toxicity may be a more appropriate assessment of vegetation indices than assessing a certain level of Li across all species. For A. thaliana, R750/R550 was the best indicator of plant response to Li exposure, as it was a significant predictor of all endpoints considered, and it also has been previously associated with alkali metal (i.e., Cs) exposure [15,69]. Similarly, although with less predictive power, PSND, R950/R750, and R1636/R1933 were also previously associated with metal exposure [69], and would be reasonable VI to consider for identifying metal stress in the future. It is likely that a combination of VI and spectra collection techniques (i.e., CP and FOV) could provide the overall best approximation of plant stress status by accounting for both whole plant and leaf optical properties; multi-index use should be given future consideration in studies utilizing both CP and FOV.
Although we found some success in the laboratory for identifying relationships between symptoms induced by low-level Li exposure and reflectance spectra, environmental and sampling conditions were controlled; therefore, care should be given if the intent is to extrapolate to field studies. Measurements taken in the field may not be as consistent or informative as measurements taken in the laboratory due to extraneous and potentially unknown environmental factors. However, this study did find that the normalized difference (R950 − R750)/(R950 + R750) provided comparable results to its corresponding simple ratio in a situation with low background reflection. In another setting, utilization of normalized differences where simple ratios were useful here may still be fitting.
Reflectance spectra provide useful information, but they can be expected to be different for distinct species. However, when treatment reflectance spectra are considered relative to a control, stress may be able to be quantified across species. The general difference in metal toxicity symptoms between all four species explains why reflectance spectra did not always shift in similar ways. Knowledge of the mechanisms involved in plant species uptake and response to a desired metal is necessary to appropriately apply vegetation indices as predictors of stress. Consideration of individual reflectance spectra as well as treatment spectra relative to control can provide additional insight into understanding results.

Supplementary Materials

The following are available online at https://www.mdpi.com/1424-8220/18/9/2750/s1, Table S1: Average reflectance data acquired by FOV and CP.

Author Contributions

M.C.D. and W.W.K. conceived the experiments and provided equipment and supplies; M.C.D., W.W.K., and N.E.M. designed the experiments; N.E.M. and C.T.S. performed the experiments; N.E.M. and J.L.S. analyzed the data; N.E.M. and J.L.S. wrote the paper; M.C.D., W.W.K., and T.E.J. provided critical review of the final manuscript.

Funding

This project was conducted in conjunction with work accomplished under Contract No. DE-AC09-08SR22470 with the U.S. Department of Energy. This work is supported by U.S. DOE—National Nuclear Security Administration through the Office of Defense Nuclear Non-proliferation Research and Development—NA-22 (Grant No. DE-FG52-05NA27035).

Conflicts of Interest

The authors declare no conflict of interest. The results were reviewed and approved for publication by funding sponsor, but the sponsor had no role in the design of the study; in the collection, analyses, or interpretation of data; or in the writing of the manuscript.

Appendix A

Figure A1. Control and treatment spectra for A. thaliana as acquired by CP and FOV.
Figure A1. Control and treatment spectra for A. thaliana as acquired by CP and FOV.
Sensors 18 02750 g0a1
Figure A2. Control and treatment spectra for Z. mays as acquired by CP and FOV.
Figure A2. Control and treatment spectra for Z. mays as acquired by CP and FOV.
Sensors 18 02750 g0a2
Figure A3. Control and treatment spectra for H. annuus as acquired by CP and FOV.
Figure A3. Control and treatment spectra for H. annuus as acquired by CP and FOV.
Sensors 18 02750 g0a3
Figure A4. Control and treatment spectra for B. napus as acquired by CP and FOV.
Figure A4. Control and treatment spectra for B. napus as acquired by CP and FOV.
Sensors 18 02750 g0a4

Appendix B

Table A1. Summary statistics for sample endpoints by week and treatment group. Cells contain the mean ± the standard deviation of the mean. For each cell n = 12, except dry biomass for A. thaliana for which n = 6 for each cell, and Chl a + b for A. thaliana in Week 3 for which n = 11 (insufficient plant material available for one sample).
Table A1. Summary statistics for sample endpoints by week and treatment group. Cells contain the mean ± the standard deviation of the mean. For each cell n = 12, except dry biomass for A. thaliana for which n = 6 for each cell, and Chl a + b for A. thaliana in Week 3 for which n = 11 (insufficient plant material available for one sample).
WeekGroupChl a + b (µg/mL)RWCDry Biomass (mg)CLAI
Zea mays
1Control24.83 ± 2.300.916 ± 0.005309.0 ± 102.10.371 ± 0.105
Treatment24.92 ± 2.760.915 ± 0.004334.7 ± 93.50.304 ± 0.086
2Control25.34 ± 2.280.912 ± 0.005502.4 ± 195.30.352 ± 0.084
Treatment26.79 ± 1.880.914 ± 0.005552.2 ± 180.60.349 ± 0.077
3Control20.84 ± 4.210.902 ± 0.009896.0 ± 368.00.378 ± 0.096
Treatment23.64 ± 2.120.905 ± 0.006818.4 ± 285.50.364 ± 0.093
Arabidopsis thaliana
1Control15.49 ± 2.180.928 ± 0.006139.1 ± 27.80.563 ± 0.065
Treatment15.42 ± 1.330.918 ± 0.010179.6 ± 77.70.505 ± 0.104
2Control18.35 ± 2.230.922 ± 0.010264.4 ± 31.60.704 ± 0.067
Treatment15.47 ± 4.340.882 ± 0.012183.4 ± 69.50.480 ± 0.109
3Control19.82 ± 2.930.916 ± 0.005279.3 ± 25.10.730 ± 0.107
Treatment6.06 ± 3.400.629 ± 0.137212.2 ± 69.10.461 ± 0.061
Brassica napus
1Control13.66 ± 1.890.946 ± 0.003271.0 ± 50.60.505 ± 0.068
Treatment13.67 ± 1.100.943 ± 0.004313.9 ± 91.50.493 ± 0.068
2Control14.41 ± 1.130.943 ± 0.008492.5 ± 146.40.617 ± 0.100
Treatment13.51 ± 2.130.934 ± 0.005520.4 ± 112.60.602 ± 0.115
3Control15.12 ± 2.840.917 ± 0.018924.7 ± 145.50.631 ± 0.080
Treatment14.45 ± 2.720.910 ± 0.019831.6 ± 165.10.672 ± 0.095
Helianthus annuus
1Control18.71 ± 2.350.932 ± 0.006262.9 ± 130.60.427 ± 0.195
Treatment17.93 ± 3.240.928 ± 0.008234.5 ± 113.20.398 ± 0.150
2Control23.92 ± 2.620.929 ± 0.008397.7 ± 109.90.565 ± 0.137
Treatment20.73 ± 3.120.931 ± 0.007345.8 ± 162.20.527 ± 0.152
3Control23.96 ± 2.880.911 ± 0.017805.8 ± 201.30.651 ± 0.180
Treatment20.79 ± 3.740.913 ± 0.013739.2 ± 218.20.556 ± 0.237
Table A2. Summary statistics for H. annuus height by week and treatment group. Cells contain the mean ± the standard deviation of the mean. For each cell, n = 12.
Table A2. Summary statistics for H. annuus height by week and treatment group. Cells contain the mean ± the standard deviation of the mean. For each cell, n = 12.
WeekGroupHeight (cm)
1Control7.76 ± 2.04
Treatment8.89 ± 2.87
2Control11.91 ± 2.80
Treatment12.30 ± 3.60
3Control24.32 ± 4.63
Treatment18.86 ± 4.54

Appendix C

Table A3. Results of initial mixed-effects analysis of VI. Significant p-values shown in bold.
Table A3. Results of initial mixed-effects analysis of VI. Significant p-values shown in bold.
SpeciesVICPFOV
Tmt (F(1,10))Week (F(2,56))Interaction (F(2,56))Tmt (F(1,10))Week (F(2,56))Interaction (F(2,56))
H. annuusNDVI------0.23850.02030.0813
WI0.50190.00690.40490.33170.03340.8654
PSND0.02950.01270.85550.31090.02170.0746
YI0.11770.61250.46280.0022<0.00010.0079
R950/R7500.3633<0.00010.74130.75590.15250.2403
R750/R5500.0071<0.00010.15730.0008<0.00010.0269
R1636/R19330.06530.88310.69970.8969<0.00010.0061
(R950 − R750)/(R950 + R750)0.3618<0.00010.73860.76190.14810.2351
Z. maysNDVI------0.75680.00060.6193
WI0.61740.08430.49980.6413<0.00010.9903
PSND0.21750.00030.07070.71040.00050.5628
YI0.49960.00300.17360.0197<0.00010.0729
R950/R7500.4220<0.00010.19040.9706<0.00010.2914
R750/R5500.85500.00020.48730.17120.00130.7066
R1636/R19330.35380.03880.25610.60070.19990.6718
(R950 − R750)/(R950 + R750)0.4221<0.00010.19100.9879<0.00010.2896
B. napusNDVI------0.50290.00040.8694
WI0.0103<0.00010.02500.11170.02130.6315
PSND0.1105<0.00010.39020.48650.00310.9486
YI0.65430.00270.02420.5517<0.00010.9030
R950/R7500.8266<0.00010.12670.1847<0.00010.1889
R750/R5500.32000.75600.04310.26390.00010.6164
R1636/R19330.0006<0.00010.00040.4581<0.00010.6011
(R950 − R750)/(R950 + R750)0.8263<0.00010.12760.1821<0.00010.1806
A. thalianaNDVI------<0.0001<0.0001<0.0001
WI0.0003<0.0001<0.00010.00340.0060<0.0001
PSND<0.0001<0.0001<0.0001<0.0001<0.0001<0.0001
YI<0.0001<0.0001<0.0001<0.0001<0.0001<0.0001
R950/R750<0.0001<0.0001<0.0001<0.0001<0.0001<0.0001
R750/R5500.0001<0.0001<0.0001<0.00010.6830<0.0001
R1636/R1933<0.0001<0.0001<0.0001<0.0001<0.0001<0.0001
(R950 − R750)/(R950 + R750)<0.0001<0.0001<0.0001<0.0001<0.0001<0.0001

References

  1. Chaerle, L.; Van Der Straeten, D. Imaging techniques and the early detection of plant stress. Trends Plant Sci. 2000, 5, 495–501. [Google Scholar] [CrossRef]
  2. Milton, N.M.; Ager, C.M.; Eiswerth, B.A.; Power, M.S. Arsenic-Induced and Selenium-Induced Changes in Spectral Reflectance and Morphology of Soybean Plants. Remote Sens. Environ. 1989, 30, 263–269. [Google Scholar] [CrossRef]
  3. Slonecker, E.T. Analysis of the effects of heavy metals on vegetation hyperspectral reflectance properties. In Hyperspectral Remote Sensing of Vegetation; Thenkabail, P.S., Lyon, J.G., Huete, A., Eds.; CRC Press: Boca Raton, FL, USA, 2012; p. 561. [Google Scholar]
  4. Agelet, L.E.; Hurburgh, C.R. A Tutorial on Near Infrared Spectroscopy and Its Calibration. Crit. Rev. Anal. Chem. 2010, 40, 246–260. [Google Scholar] [CrossRef]
  5. Salt, D.E.; Blaylock, M.; Kumar, N.P.B.A.; Dushenkov, V.; Ensley, B.D.; Chet, I.; Raskin, I. Phytoremediation: A Novel Strategy for the Removal of Toxic Metals from the Environment Using Plants. Nat. Biotechnol. 1995, 13, 468. [Google Scholar] [CrossRef]
  6. Cameron, J.R.; Kenney, G.; Johnson, N.; Daniels, F. Radiation Dosimeter Utilizing Thermoluminescence of Lithium Fluoride. Science 1961, 134, 333–334. [Google Scholar] [CrossRef] [PubMed]
  7. Hawrylak-Nowak, B.; Kalinowska, M.; Szymanska, M. A Study on Selected Physiological Parameters of Plants Grown Under Lithium Supplementation. Biol. Trace Elem. Res. 2012, 149, 425–430. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  8. Kszos, L.A.; Stewart, A.J. Review of lithium in the aquatic environment: Distribution in the united states, toxicity and case example of groundwater contamination. Ecotoxicology 2003, 12, 439–447. [Google Scholar] [CrossRef] [PubMed]
  9. Aral, H.; Vecchio-Sadus, A. Toxicity of lithium to humans and the environment—A literature review. Ecotoxicol. Environ. Saf. 2008, 70, 349–356. [Google Scholar] [CrossRef] [PubMed]
  10. Crowley, K.D.; Ahearne, J.F. Managing the environmental legacy of US nuclear-weapons production—Although the waste from America’s arms buildup will never be “cleaned up,” human and environmental risks can be reduced and managed. Am. Sci. 2002, 90, 514–523. [Google Scholar] [CrossRef]
  11. Balonov, M.; Dubourg, M.; Efremenkov, V.; Gilpin, J. Management of Waste Containing Tritium and Carbon-14; Technical Reports Series; IAEA: Vienna, Austria, 2004. [Google Scholar]
  12. Woodhouse, R.; Heeb, M.; Berry, W.; Hoshizaki, T.; Wood, M. Analysis of Remote Reflection Spectroscopy to Monitor Plant Health. Adv. Space Res. 1994, 14, 199–202. [Google Scholar] [CrossRef]
  13. Bandaru, V.; Hansen, D.J.; Codling, E.E.; Daughtry, C.S.; White-Hansen, S.; Green, C.E. Quantifying arsenic-induced morphological changes in spinach leaves: Implications for remote sensing. Int. J. Remote Sens. 2010, 31, 4163–4177. [Google Scholar] [CrossRef]
  14. Collins, W.; Chang, S.H.; Raines, G.; Canney, F.; Ashley, R. Airborne Biogeophysical Mapping of Hidden Mineral-Deposits. Econ. Geol. 1983, 78, 737–749. [Google Scholar] [CrossRef]
  15. Davids, C.; Tyler, A.N. Detecting contamination-induced tree stress within the Chernobyl exclusion zone. Remote Sens. Environ. 2003, 85, 30–38. [Google Scholar] [CrossRef]
  16. Dunagan, S.C.; Gilmore, M.S.; Varekamp, J.C. Effects of mercury on visible/near-infrared reflectance spectra of mustard spinach plants (Brassica rapa P.). Environ. Pollut. 2007, 148, 301–311. [Google Scholar] [CrossRef] [PubMed]
  17. Horler, D.N.H.; Dockray, M.; Barber, J. The red edge of plant leaf reflectance. Int. J. Remote Sens. 1983, 4, 273–288. [Google Scholar] [CrossRef]
  18. Kooistra, L.; Salas, E.A.L.; Clevers, J.G.P.W.; Wehrens, R.; Leuven, R.S.E.W.; Nienhuis, P.H.; Buydens, L.M.C. Exploring field vegetation reflectance as an indicator of soil contamination in river floodplains. Environ. Pollut. 2004, 127, 281–290. [Google Scholar] [CrossRef]
  19. Sridhar, B.B.M.; Han, F.X.; Diehl, S.V.; Monts, D.L.; Su, Y. Monitoring the effects of arsenic and chromium accumulation in Chinese brake fern (Pteris vittata). Int. J. Remote Sens. 2007, 28, 1055–1067. [Google Scholar] [CrossRef]
  20. Sridhar, B.B.M.; Han, F.X.; Diehl, S.V.; Monts, D.L.; Su, Y. Spectral reflectance and leaf internal structure changes of barley plants due to phytoextraction of zinc and cadmium. Int. J. Remote Sens. 2007, 28, 1041–1054. [Google Scholar] [CrossRef]
  21. Sridhar, B.B.M.; Vincent, R.K.; Roberts, S.J.; Czajkowski, K. Remote sensing of soybean stress as an indicator of chemical concentration of biosolid amended surface soils. Int. J. Appl. Earth Obs. Geoinf. 2011, 13, 676–681. [Google Scholar] [CrossRef]
  22. Milton, N.M.; Eiswerth, B.A.; Ager, C.M. Effect of Phosphorus Deficiency on Spectral Reflectance and Morphology of Soybean Plants. Remote Sens. Environ. 1991, 36, 121–127. [Google Scholar] [CrossRef]
  23. Schwaller, M.R.; Schnetzler, C.C.; Marshall, P.E. The changes in leaf reflectance of sugar maple (Acer-saccharum Marsh) seedlings in response to heavy-metal stress. Int. J. Remote Sens. 1983, 4, 93–100. [Google Scholar] [CrossRef]
  24. Su, Y.; Sridhar, B.B.M.; Han, F.X.; Diehl, S.V.; Monts, D.L. Effect of bioaccumulation of Cs and Sr natural isotopes on foliar structure and plant spectral reflectance of Indian mustard (Brassica juncea). Water Air Soil Pollut. 2007, 180, 65–74. [Google Scholar] [CrossRef]
  25. NSF. Arabidopsis: The Model Plant. Available online: http://www.nsf.gov/pubs/2002/bio0202/model.htm#content (accessed on 5 December 2017).
  26. Li, X.; Gao, P.; Gjetvaj, B.; Westcott, N.; Gruber, M.Y. Analysis of the metabolome and transcriptome of Brassica carinata seedlings after lithium chloride exposure. Plant Sci. 2009, 177, 68–80. [Google Scholar] [CrossRef]
  27. Marchiol, L.; Assolari, S.; Sacco, P.; Zerbi, G. Phytoextraction of heavy metals by canola (Brassica napus) and radish (Raphanus sativus) grown on multicontaminated soil. Environ. Pollut. 2004, 132, 21–27. [Google Scholar] [CrossRef] [PubMed]
  28. Memon, A.R.; Schroder, P. Implications of metal accumulation mechanisms to phytoremediation. Environ. Sci. Pollut. Res. 2009, 16, 162–175. [Google Scholar] [CrossRef] [PubMed]
  29. Vanaja, M.; Yadav, S.K.; Archana, G.; Lakshmi, N.J.; Reddy, P.R.R.; Vagheera, P.; Razak, S.K.A.; Maheswari, M.; Venkateswarlu, B. Response of C-4 (maize) and C-3 (sunflower) crop plants to drought stress and enhanced carbon dioxide concentration. Plant Soil Environ. 2011, 57, 207–215. [Google Scholar] [CrossRef]
  30. Prasad, M.N.V. Sunflower (Helianthus annuus L.)—A potential crop for environmental industry. Helia 2007, 30, 167–174. [Google Scholar]
  31. Wuana, R.A.; Okieimen, F.E. Phytoremediation potential of maize (Zea mays L.). A review. Afr. J. Gen. Agric. 2010, 6, 275–287. [Google Scholar]
  32. Capehart, T. Corn & Other Feedgrains. Available online: https://www.ers.usda.gov/topics/crops/corn-and-other-feedgrains/ (accessed on 7 June 2018).
  33. Grzesiak, M.T.; Filek, W.; Hura, T.; Kocurek, M.; Pilarski, J. Leaf optical properties during and after drought stress in triticale and maize genotypes differing in drought tolerance. Acta Physiol. Plant. 2010, 32, 433–442. [Google Scholar] [CrossRef]
  34. Peng, Y.; Gitelson, A.A. Remote estimation of gross primary productivity in soybean and maize based on total crop chlorophyll content. Remote Sens. Environ. 2012, 117, 440–448. [Google Scholar] [CrossRef]
  35. Al-Abbas, A.H.; Barr, R.; Hall, J.D.; Crane, F.L.; Baumgardner, M.F. Spectra of Normal and Nutrient-Deficient Maize Leaves. Agron. J. 1974, 66, 16–20. [Google Scholar] [CrossRef]
  36. Ercoli, L.; Mariotti, M.; Masoni, A.; Massantini, F. Relationship between nitrogen and chlorophyll content and spectral properties in maize leaves. Eur. J. Agron. 1993, 2, 113–117. [Google Scholar] [CrossRef]
  37. Blackmer, T.M.; Schepers, J.S.; Varvel, G.E. Light Reflectance Compared with Other Nitrogen Stress Measurements in Corn Leaves. Agron. J. 1994, 86, 934–938. [Google Scholar] [CrossRef]
  38. Penuelas, J.; Gamon, J.A.; Fredeen, A.L.; Merino, J.; Field, C.B. Reflectance Indexes Associated with Physiological-Changes in Nitrogen-Limited and Water-Limited Sunflower Leaves. Remote Sens. Environ. 1994, 48, 135–146. [Google Scholar] [CrossRef]
  39. Penuelas, J.; Baret, F.; Filella, I. Semiempirical Indexes to Assess Carotenoids Chlorophyll—A Ratio from Leaf Spectral Reflectance. Photosynthetica 1995, 31, 221–230. [Google Scholar]
  40. Masoni, A.; Ercoli, L.; Mariotti, M.; Barberi, P. Changes in Spectral Properties of Aging and Senescing Maize and Sunflower Leaves. Physiol. Plant. 1994, 91, 334–338. [Google Scholar] [CrossRef]
  41. Mariotti, M.; Ercoli, L.; Masoni, A. Spectral properties of iron-deficient corn and sunflower leaves. Remote Sens. Environ. 1996, 58, 282–288. [Google Scholar] [CrossRef]
  42. Daughtry, C.S.T.; Walthall, C.L.; Kim, M.S.; de Colstoun, E.B.; McMurtrey, J.E. Estimating corn leaf chlorophyll concentration from leaf and canopy reflectance. Remote Sens. Environ. 2000, 74, 229–239. [Google Scholar] [CrossRef]
  43. Reddy, G.S.; Rao, C.L.N.; Venkataratnam, L.; Rao, P.V.K. Influence of plant pigments on spectral reflectance of maize, groundnut and soybean grown in semi-arid environments. Int. J. Remote Sens. 2001, 22, 3373–3380. [Google Scholar] [CrossRef]
  44. Zhao, D.L.; Reddy, K.R.; Kakani, V.G.; Read, J.J.; Carter, G.A. Corn (Zea mays L.) growth, leaf pigment concentration, photosynthesis and leaf hyperspectral reflectance properties as affected by nitrogen supply. Plant Soil 2003, 257, 205–217. [Google Scholar] [CrossRef]
  45. Schlemmer, M.; Gitelson, A.; Schepers, J.; Ferguson, R.; Peng, Y.; Shanahan, J.; Rundquist, D. Remote estimation of nitrogen and chlorophyll contents in maize at leaf and canopy levels. Int. J. Appl. Earth Obs. Geoinf. 2013, 25, 47–54. [Google Scholar] [CrossRef] [Green Version]
  46. Schrauzer, G.N. Lithium: Occurrence, dietary intakes, nutritional essentiality. J. Am. Coll. Nutr. 2002, 21, 14–21. [Google Scholar] [CrossRef] [PubMed]
  47. Mcstay, N.G.; Rogers, H.H.; Anderson, C.E. Effects of Lithium on Phaseolus-vulgaris L. Sci. Total Environ. 1980, 16, 185–191. [Google Scholar] [CrossRef]
  48. Bingham, F.T.; Bradford, G.R.; Page, A.L. Toxicity of Lithium to Plants. Calif. Agric. 1964, 18, 6–7. [Google Scholar]
  49. Naranjo, M.A.; Romero, C.; Belles, J.M.; Montesinos, C.; Vicente, O.; Serrano, R. Lithium treatment induces a hypersensitive-like response in tobacco. Planta 2003, 217, 417–424. [Google Scholar] [CrossRef] [PubMed]
  50. Kent, N.L. The influence of lithium salts on certain cultivated plants and their parasitic diseases. Ann. Appl. Biol. 1941, 28, 189–209. [Google Scholar] [CrossRef]
  51. Kent, N.L. Absorption, translocation and ultimate fate of lithium in the wheat plant. New Phytol. 1941, 40, 291–298. [Google Scholar] [CrossRef]
  52. White, P.J.; Broadley, M.R. Chloride in soils and its uptake and movement within the plant: A review. Ann. Bot. 2001, 88, 967–988. [Google Scholar] [CrossRef]
  53. Boyko, A.; Golubov, A.; Bilichak, A.; Kovalchuk, I. Chlorine ions but not sodium ions alter genome stability of Arabidopsis thaliana. Plant Cell Physiol. 2010, 51, 1066–1078. [Google Scholar] [CrossRef] [PubMed]
  54. Bueso, E.; Alejandro, S.; Carbonell, P.; Perez-Amador, M.A.; Fayos, J.; Belles, J.M.; Rodriguez, P.L.; Serrano, R. The lithium tolerance of the Arabidopsis cat2 mutant reveals a cross-talk between oxidative stress and ethylene. Plant J. 2007, 52, 1052–1065. [Google Scholar] [CrossRef] [PubMed]
  55. Mulkey, T.J. Alteration of growth and gravitropic response of maize roots by lithium. Gravit. Space Res. 2005, 18, 119–120. [Google Scholar]
  56. Duff, M.C.; Kuhne, W.W.; Halverson, N.V.; Chang, C.-S.; Kitamura, E.; Hawthorn, L.; Martinez, N.E.; Stafford, C.; Milliken, C.E.; Caldwell, E.F.; et al. mRNA Transcript abundance during plant growth and the influence of Li+ exposure. Plant Sci. 2014, 229, 262–279. [Google Scholar] [CrossRef] [PubMed]
  57. Sneva, F.A. Lithium Toxicity in Seedlings of 3 Cool Season Grasses. Plant Soil 1979, 53, 219–224. [Google Scholar] [CrossRef]
  58. Bartolo, M.E.; Carter, J.V. Lithium Decreases Cold-Induced Microtubule Depolymerization in Mesophyll-Cells of Spinach. Plant Physiol. 1992, 99, 1716–1718. [Google Scholar] [CrossRef] [PubMed]
  59. ASD. FieldSpec® Pro: User’s Guide; Revision, C., Ed.; Analytical Spectral Devices: Boulder, CO, USA, 2002. [Google Scholar]
  60. Knudson, L.L.; Tibbitts, T.W.; Edwards, G.E. Measurement of Ozone Injury by Determination of Leaf Chlorophyll Concentration. Plant Physiol. 1977, 60, 606–608. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  61. Papista, E.; Acs, E.; Boddi, B. Chlorophyll-alpha determination with ethanol—A critical test. Hydrobiologia 2002, 485, 191–198. [Google Scholar] [CrossRef]
  62. Penuelas, J.; Filella, I.; Lloret, P.; Munoz, F.; Vilajeliu, M. Reflectance Assessment of Mite Effects on Apple-Trees. Int. J. Remote Sens. 1995, 16, 2727–2733. [Google Scholar] [CrossRef]
  63. Rouse, J.W., Jr.; Haas, R.H.; Deering, D.W.; Schell, J.A.; Harlan, J.C. Monitoring the Vernal Advancement and Retrogradation (Green Wave Effect) of Natural Vegetation; Remote Sensing Center, Texas A&M University: College Station, TX, USA, 1974. [Google Scholar]
  64. Gamon, J.A.; Penuelas, J.; Field, C.B. A narrow-waveband spectral index that tracks diurnal changes in photosynthetic efficiency. Remote Sens. Environ. 1992, 41, 35–44. [Google Scholar] [CrossRef]
  65. Serrano, L. Effects of leaf structure on reflectance estimates of chlorophyll content. Int. J. Remote Sens. 2008, 29, 5265–5274. [Google Scholar] [CrossRef]
  66. Penuelas, J.; Pinol, J.; Ogaya, R.; Filella, I. Estimation of plant water concentration by the reflectance water index WI (R900/R970). Int. J. Remote Sens. 1997, 18, 2869–2875. [Google Scholar] [CrossRef]
  67. Adams, M.L.; Philpot, W.D.; Norvell, W.A. Yellowness index: An application of spectral second derivatives to estimate chlorosis of leaves in stressed vegetation. Int. J. Remote Sens. 1999, 20, 3663–3675. [Google Scholar] [CrossRef]
  68. Gitelson, A.; Merzlyak, M.N. Spectral Reflectance Changes Associated with Autumn Senescence of Aesculus-hippocastanum L. and Acer-platanoides L. Leaves. Spectral Features and Relation to Chlorophyll Estimation. J. Plant Physiol. 1994, 143, 286–292. [Google Scholar] [CrossRef]
  69. Martinez, N.E.; Sharp, J.L.; Kuhne, W.W.; Johnson, T.E.; Stafford, C.T.; Duff, M.C. Assessing the use of reflectance spectroscopy in determining CsCl stress in the model species Arabidopsis thaliana. Int. J. Remote Sens. 2015, 36, 5887–5915. [Google Scholar] [CrossRef]
  70. Blackburn, G.A. Quantifying chlorophylls and caroteniods at leaf and canopy scales: An evaluation of some hyperspectral approaches. Remote Sens. Environ. 1998, 66, 273–285. [Google Scholar] [CrossRef]
  71. 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]
  72. Shahzad, B.; Tanveer, M.; Hassan, W.; Shah, A.N.; Anjum, S.A.; Cheema, S.A.; Ali, I. Lithium toxicity in plants: Reasons, mechanisms and remediation possibilities—A review. Plant Physiol. Biochem. 2016, 107, 104–115. [Google Scholar] [CrossRef] [PubMed]
  73. Hermans, C.; Verbruggen, N. Physiological characterization of Mg deficiency in Arabidopsis thaliana. J. Exp. Bot. 2005, 56, 2153–2161. [Google Scholar] [CrossRef] [PubMed]
  74. Appenroth, K. Definition of “heavy metals” and their role in biological systems. Soil Heavy Met. 2010, 19, 19–29. [Google Scholar]
  75. Kumar, L.; Schmidt, K.; Dury, S.; Skidmore, A. Imaging spectrometry and vegetation science. In Imaging Spectrometry; Van der Meer, F.D., de Jong, S.M., Eds.; Springer: Dordrecht, The Netherlands, 2006; p. 111. [Google Scholar]
  76. Knipling, E.B. Physical and physiological basis for the reflectance of visible and near-infrared radiation from vegetation. Remote Sens. Environ. 1970, 1, 155–159. [Google Scholar] [CrossRef]
  77. Slaton, M.R.; Hunt, E.R.; Smith, W.K. Estimating near-infrared leaf reflectance from leaf structural characteristics. Am. J. Bot. 2001, 88, 278–284. [Google Scholar] [CrossRef] [PubMed]
Figure 1. Relative reflectance of treatment to control for each of four species as acquired by CP.
Figure 1. Relative reflectance of treatment to control for each of four species as acquired by CP.
Sensors 18 02750 g001
Figure 2. Relative reflectance of treatment to control for each of four species as acquired by FOV.
Figure 2. Relative reflectance of treatment to control for each of four species as acquired by FOV.
Sensors 18 02750 g002
Figure 3. Weeks 1, 2, and 3 of the first A. thaliana experiment; response was similar in the second experiment. (a) Week 1 control plants; (b) Week 1 treatment plants which exhibit slight chlorosis along the edges of older leaves; (c) Week 2 control plants; (d) Week 2 treatment plants which exhibit chlorosis, necrosis, and decreased biomass; (e) Week 3 control plants; (f) Week 3 plants which exhibit significant necrosis.
Figure 3. Weeks 1, 2, and 3 of the first A. thaliana experiment; response was similar in the second experiment. (a) Week 1 control plants; (b) Week 1 treatment plants which exhibit slight chlorosis along the edges of older leaves; (c) Week 2 control plants; (d) Week 2 treatment plants which exhibit chlorosis, necrosis, and decreased biomass; (e) Week 3 control plants; (f) Week 3 plants which exhibit significant necrosis.
Sensors 18 02750 g003
Figure 4. Example (a) B. napus and (b) H. annuus plants at week 3. Arrows in (a) are pointing out areas of chlorosis around leaf edges. The arrow in (b) is pointing out a necrotic spot.
Figure 4. Example (a) B. napus and (b) H. annuus plants at week 3. Arrows in (a) are pointing out areas of chlorosis around leaf edges. The arrow in (b) is pointing out a necrotic spot.
Sensors 18 02750 g004
Figure 5. Chlorophyll content (Chl a + b, µg/mL) by week and treatment level for each species considered. Lines (solid = control, dashed = treatment) indicate mean values at each week.
Figure 5. Chlorophyll content (Chl a + b, µg/mL) by week and treatment level for each species considered. Lines (solid = control, dashed = treatment) indicate mean values at each week.
Sensors 18 02750 g005
Figure 6. Relative water content (RWC) by week and treatment level for each species considered. Lines (solid = control, dashed = treatment) indicate mean values at each week.
Figure 6. Relative water content (RWC) by week and treatment level for each species considered. Lines (solid = control, dashed = treatment) indicate mean values at each week.
Sensors 18 02750 g006
Figure 7. Dry biomass (mg) by week and treatment level for each species considered. Lines (solid = control, dashed = treatment) indicate mean values at each week.
Figure 7. Dry biomass (mg) by week and treatment level for each species considered. Lines (solid = control, dashed = treatment) indicate mean values at each week.
Sensors 18 02750 g007
Figure 8. CLAI by week and treatment level for each species considered. Lines (solid = control, dashed = treatment) indicate mean values at each week.
Figure 8. CLAI by week and treatment level for each species considered. Lines (solid = control, dashed = treatment) indicate mean values at each week.
Sensors 18 02750 g008
Figure 9. Height by week and treatment level for H. annuus. Lines (solid = control, dashed = treatment) indicate mean values at each week.
Figure 9. Height by week and treatment level for H. annuus. Lines (solid = control, dashed = treatment) indicate mean values at each week.
Sensors 18 02750 g009
Table 1. Results (p-values) of initial statistical analysis of measured endpoints. Bold numbers indicate significance at the 0.05 level.
Table 1. Results (p-values) of initial statistical analysis of measured endpoints. Bold numbers indicate significance at the 0.05 level.
SpeciesEndpointTmt
(F(1,10))
Week
(F(2,56))
Interaction (F(2,56))
A. thalianaChl a + b<0.0001<0.0001<0.0001
RWC<0.0001<0.0001<0.0001
Dry biomass0.0479<0.00010.0014
CLAI<0.00010.03090.0002
Z. maysChl a + b0.1067<0.00010.1922
RWC0.2580<0.00010.4192
Dry biomass0.9923<0.00010.5890
CLAI0.26520.43380.4271
H. annuusChl a + b0.0014<0.00010.2899
RWC0.9459<0.00010.5909
Dry biomass0.4353<0.00010.8903
CLAI0.36690.00060.7505
Height0.1703<0.00010.0031
B. napusChl a + b0.34590.15730.7325
RWC0.1600<0.00010.5512
Dry biomass0.8766<0.00010.0599
CLAI0.8302<0.00010.4791
Table 2. Results of secondary mixed-effects analysis. Significant p-values shown in bold.
Table 2. Results of secondary mixed-effects analysis. Significant p-values shown in bold.
EndpointVIH. annuusA. thaliana
CPFOVCPFOV
Chl a + bNDVI------0.3617
WI----0.81540.9284
PSND0.0063--0.02000.3836
YI--0.02020.00030.1850
R950/R750----0.06610.4488
R750/R550<0.00010.00300.02050.0235
R1636/R1933----0.02560.7458
(R950 − R750)/(R950 + R750)----0.09370.3938
RWCNDVI------0.5406
WI----0.99370.4283
PSND----0.53330.7267
YI----0.92890.5922
R950/R750----0.02570.0632
R750/R550----0.64460.0194
R1636/R1933----0.33560.6893
(R950 − R750)/(R950 + R750)----0.04020.0572
Dry biomassNDVI------0.5192
WI----0.15650.9351
PSND----0.69890.5206
YI----0.39000.0007
R950/R750----0.79970.1980
R750/R550----<0.0023<0.0001
R1636/R1933----0.55270.6381
(R950 − R750)/(R950 + R750)----0.76800.1066
CLAINDVI------0.0040
WI----0.56400.3781
PSND----0.12820.0058
YI----0.24480.9941
R950/R750----0.46030.0663
R750/R550----0.00030.0415
R1636/R1933----0.23420.0049
(R950 − R750)/(R950 + R750)----0.37390.0266
HeightNDVI--------
WI--------
PSND<0.0001------
YI--0.0229----
R950/R750--------
R750/R550<0.0001<0.0001----
R1636/R1933--------
(R950 − R750)/(R950 + R750)--------
-- Not included in secondary analysis.

Share and Cite

MDPI and ACS Style

Martinez, N.E.; Sharp, J.L.; Johnson, T.E.; Kuhne, W.W.; Stafford, C.T.; Duff, M.C. Reflectance-Based Vegetation Index Assessment of Four Plant Species Exposed to Lithium Chloride. Sensors 2018, 18, 2750. https://doi.org/10.3390/s18092750

AMA Style

Martinez NE, Sharp JL, Johnson TE, Kuhne WW, Stafford CT, Duff MC. Reflectance-Based Vegetation Index Assessment of Four Plant Species Exposed to Lithium Chloride. Sensors. 2018; 18(9):2750. https://doi.org/10.3390/s18092750

Chicago/Turabian Style

Martinez, Nicole E., Julia L. Sharp, Thomas E. Johnson, Wendy W. Kuhne, Clay T. Stafford, and Martine C. Duff. 2018. "Reflectance-Based Vegetation Index Assessment of Four Plant Species Exposed to Lithium Chloride" Sensors 18, no. 9: 2750. https://doi.org/10.3390/s18092750

APA Style

Martinez, N. E., Sharp, J. L., Johnson, T. E., Kuhne, W. W., Stafford, C. T., & Duff, M. C. (2018). Reflectance-Based Vegetation Index Assessment of Four Plant Species Exposed to Lithium Chloride. Sensors, 18(9), 2750. https://doi.org/10.3390/s18092750

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