Remote Sensing Retrieval of Leaf Area Index (lai) and Soil Water Content (wc) Using Hyperspectral Remote Sensing under Controlled Glass House Conditions for Spring Barley and Sugar Beet

Leaf area index (LAI) and water content (WC) in the root zone are two major hydro-meteorological parameters that exhibit a dominant control on water, energy and carbon fluxes, and are therefore important for any regional eco-hydrological or climatological study. To investigate the potential for retrieving these parameter from hyperspectral remote sensing, we have investigated plant spectral reflectance (400–2,500 nm, ASD FieldSpec3) for two major agricultural crops (sugar beet and spring barley) in the mid-latitudes, treated under different water and nitrogen (N) conditions in a greenhouse experiment over the growing period of 2008. Along with the spectral response, we have measured soil water content and LAI for 15 intensive measurement campaigns spread over the growing season and could demonstrate a significant response of plant reflectance characteristics to variations in water content and nutrient conditions. Linear and non-linear dimensionality analysis suggests that the full band reflectance information is well represented by the set of 28 vegetation spectral indices (SI) and most of the variance is explained by three to a maximum of eight variables. Investigation of linear dependencies between LAI and soil WC and pre-selected SI's indicate that: (1) linear regression using single SI is not sufficient to describe plant/soil variables over the range of experimental conditions, however, some improvement can be seen knowing crop species beforehand; (2) the improvement is superior when applying multiple linear regression using three 1703 explanatory SI's approach. In addition to linear investigations, we applied the non-linear CART (Classification and Regression Trees) technique, which finally did not show the potential for any improvement in the retrieval process.

explanatory SI's approach.In addition to linear investigations, we applied the non-linear CART (Classification and Regression Trees) technique, which finally did not show the potential for any improvement in the retrieval process.

Introduction
The quantitative description of water, carbon and energy fluxes between the land surface and the atmosphere plays an important role in current hydrological and climate research studies, especially at larger spatial scales [1].An effective tool to derive relevant land surface characteristics over the last decades has been the development and use of satellite remote sensing techniques in the visible (VIS: 350-700 nm) and near infrared (NIR: 700-1,300 nm), short wave infrared (SWIR: 1,300-2,500 nm), microwave and thermal infrared spectral region of the electromagnetic spectrum (see e.g., [2] for an extensive review of methods in the field of eco-hydrology).In particular, the monitoring of narrow band spectral reflectance from plant surfaces using multi-and hyper-spectral sensors has proven very promising in gathering a large variety of eco-physiological information including photosynthetic activity (absorbed photosynthetically active radiation, APAR) [3], vegetation structure (leaf area index, LAI) [3,4], biomass and ecosystem productivity (net ecosystem productivity/gross primary productivity, NEP/GPP) [5], plant/leaf water content [6] or vegetation light use efficiency (LUE) [7], which is the ratio of GPP to APAR and controlled by the phenological stage, solar radiation conditions and the availability of soil moisture and nutrients [8].
The relationship between spectral reflectance and ecophysiological characteristics is often modeled via the use of spectral indices (SI's) that are calculated from spectral reflectance of multiple bands.There have been a large number of investigations over the last years and a summary of developed and applied SI's is given in Table 1.Spectral indices are here grouped into three major classes: i.The first group, including the normalized differential vegetation index (NDVI) as the most prominent representative, shows response to the general greenness, biomass and structure of the vegetation.The NDVI and its variants have successfully been related to properties such as the leaf area index (LAI), the fraction of absorbed photosynthetically active radiation (FPAR) or the biomass of many different ecosystems and environments [9,10].ii.The second group of indices is related to the leaf pigment activity (such as xanthophyll or carotenoid) that shows sensitivity to plant physiological processes and in particular to the photosynthetic radiation use efficiency, a major component in many current eco-climatic models and analysis.Specifically, the photochemical reflectance index (PRI, the normalized difference of the 531 and 570 nm bands) has received large attention over the last years allowing relating spectral response to carbon fluxes and GPP [7,[11][12][13].As the PRI is sensitive to the plant xanthophyll cycle, active when dissipating excess energy under intensive radiation, and also dependent on water and nutrient availability, it will be a very useful indicator of the soil moisture and nutrient regime at a particular location [6,14].iii.The third group of indices mainly uses water absorption bands in the near-and mid-infrared region.
They are sensitive to the leaf/plant water concentration that is controlled by the soil water availability and climatic conditions.A detailed discussion about pros and cons of individual spectral bands/indices can be found in [6,15,16].In general, the indices of this group are applied in areas such as drought assessment, irrigation practice or wild fire risk [17].

Spectral Index Index Name Equation Reference
Greenness/Biomass/Canopy structure:
While both, LAI and root zone WC, are increasingly derived using the inversion of canopy reflectance models in the visible and infrared (LAI, see [36] for an overview) and in the microwave spectral region (WC, e.g., [37][38][39]), we here want to restrict our analysis to the use of hyperspectral reflectance data and in particular to the use of individual narrow band ratios (Table 1).There are two main reasons for this: First, the use of simple SI's will be significantly less demanding in terms of computational resources when compared to the inversion of radiation transfer models such as PROSAIL [40] and might therefore be necessary for product development, e.g., with the launch of EnMAP (Environmental Mapping and Analysis Program) envisaged for 2013.Also, the use of simple indices would allow, e.g., the applications of low cost camera systems including narrow band optical filters serving as an alternative for cheap and widespread application.
Therefore, one major aim of this study is to find the most significant and lowest dimensional spectral band combinations for extracting the LAI and WC.We have investigated two major agricultural crops in the mid-latitudes, namely sugar beet and spring barley, under different water and N-treatments in a greenhouse environment over the vegetation period of 2008.Different water and N-treatments will allow a broad spectrum of environmental conditions to be explored.Section 2 will describe the experimental setup, the used instruments, and the measured variables in more detail.Section 3 presenting the experimental results gives an extensive analysis of the temporal dynamics of differences in spectral reflectance for different experimental conditions.Here, we also try to extract the dominant spectral ranges and vegetation indices that allow soil moisture content in the root zone and N conditions to be derived from remote sensing.In particular we would like to derive and test an empirical relationship between spectral reflectance information and the soil WC and LAI.This is followed by a discussion and conclusion in Section 4.

Study Site and Experimental Setup
Experimental data have been obtained during the growing season 2008 at the field research station of the UFZ-Helmholtz Centre for Environmental Research in Bad Lauchstädt (BL), Saxony-Anhalt, Germany (51°23′32″N 11°52′33″E, altitude 113 m a.s.l., mean annual precipitation 484 mm, mean annual temperature 8.7 °C; [41]).Two main German agricultural crop species, spring barley (Hordeum vulgare L.) and sugar beet (Beta vulgaris ssp.), were chosen as objects of investigation.A total of 24 pots were planted using black soil in circular-shaped pots, 0.6 m in diameter and 0.8 m high, 12 for spring barley and 12 for sugar beet.There were four sugar beet plants per pot and spring barley was planted in rows, allowing five rows per pot separated by 0.125 m as illustrated in Figure 1.Selected crops were seeded and grown in the winter greenhouse equipped with removable walls and roof that were left open during the day and closed during the night and rainfall events.This allowed plant growth under most possible natural environment conditions while being able to examine crop response to artificially produced water stress under different nutrient conditions.After seeding at the beginning/mid of April, all pots were initially cultivated under normal, well watered conditions (defined by 60% of soil field capacity that has been derived from 24 h capillary saturation of 100 cm 3 soil cores, and 2 h ex-filtration on a 6 cm sandbed).The spectral and other measurements under controlled conditions started in the first week of June when soil/plot covering by plants was almost complete.Spring barley was in the flowering stage roughly 6-7 weeks after first emergence.For sugar beet the measurements started when 4-6 leaves were unrolled and fully developed.Table 2 gives a detailed overview about phenological stages and the time schedule of the experiment.At the start of the experiment, pots were treated with different water and nutrient conditions.Two variants concerning nutrient treatment and three different modes of irrigation were considered.The Bad Lauchstädt soil used in this experiment is Haplic Chernozem (FAO classification) developed from loess [42] with 22% clay [41] characterized by a very high organic carbon (C org = 19.8g kg −1 ) and total nitrogen (N tot = 1.6 g kg −1 ) contents [43].Before seeding, all sample-pots were filled with 150 kg of soil material, 12 pots (six for sugar beet and six for spring barley) received an addition of mineral fertilizer, while the other 12 pots were kept as a reference without fertilizer.The amount of mineral fertilizer applied was defined in accordance with the "Standard Bad Lauchstädt" treatment praxis (15 g P and N as well as 20 g K per pot).The two fertilizer treatments are called "N-low" or "N-high" in the following.Three different water treatments have been applied defined as follows: a "high" WC was defined by controlling WC at a value of 60% of the soil field capacity, for a "medium" variant this value was set to 30% of field capacity, and "low" variant was defined by no additional water supply.The water status/content was verified every second-third day by weighting of the pots, followed by an irrigation to compensate for the evapotranspiration losses.In total, 24 pots including two different crop species under two nutrient and three irrigation treatments are considered and investigated, and each combination has two replicates.

Instrumentation and Measurements
Top of canopy spectral signatures were taken using the FieldSpec Pro 3 (Analytical Spectral Devices, Boulder, CO, USA) field portable spectrometer with a full spectral range located in the VIS, NIR and SWIR region (350-2,500 nm).Instrument sampling intervals (1.4 nm-VIS and NIR; 2 nm-SWIR) undergo an automatic internal interpolation to 1 nm spacing using linear and polynomial equations [44].All spectral surveys were taken under two conditions: First, reflectance was recorded from nadir position about 0.45-0.55m above the top of the plant under sun illumination in the winter greenhouse with all walls and the roof removed.Second, two tungsten halogen quartz lamps: THQ (Dörr, Neu-Ulm, Germany) with 1,000 watt, 3,200 K color temperature, 26,000 Lumens of brightness each (http://www.doerrfoto.de/)under closed room conditions were used.Both lamps were installed in a closed greenhouse room with walls covered with black mat paint.Lamps were pointed from two opposite directions and mounted at a 30° angle according to the floor plane.For reflectance scanning under THQ illumination, every sample pot had to be carried into the dark room using special lift truck designed to carry heavy load.Reference measurements were taken for both modes using a Spectralon (LabSphere, North Sutton, NH, USA) white reference panel (Effective Spectral Range: 250-2,500 nm, Reflectance within 400-1,500 nm: 99% and within 250-2,500 nm: more than 95%, Thermal Stability: up to 400 °C) at the beginning of each time series and after changes of sun illumination conditions under mode 1. Spectral data were collected throughout the experiment between 11:00 a.m. and 1:00 p.m. in order to avoid large impacts of bidirectional effects.Each reflectance spectra was calculated as an average over four individual spectral measurements (each averaged from 25 internal surveys) slightly varied around the pot center to avoid boundary effects per pot.
LAI was measured using a LAI 2,000 plant canopy analyzer (LI-COR, NE, USA) directly after the spectral measurement under outside but shadowed conditions.Measurements were carried out considering individual calibration of the instrument in a way that the hemispherical field of view was least influenced by the experimental setup and conditions.Again, each data set was averaged over four individual measurements per pot.Soil moisture was measured as volumetric soil WC using a ThetaProbe ML2x FD-probe (Delta-T Devices Ltd, Cambridge, UK) for every sampling day.A probe length of 60 mm and a set of 4-10 measurements per pot (with regards to possible soil penetration dependent on soil softness/dryness) guaranteed an average water content in the upper root zone of the pot.In addition, vegetation height, and in case of sugar beet: leaves length, were monitored.
The glasshouse experiment was started in the first week of June 2008 and was continued over a six weeks period until the second week of July.Measurements were taken at 15 individual dates however, as spectral reflectance measurements under outside conditions were limited due to unfavorable weather conditions, and as the reflectance spectra and derived SI's did not show any significant differences between both illumination conditions, we have limited our further analysis to the use of data obtained under indoor conditions.Besides direct reflectance measurements, the 28 spectral indices as defined in Table 1 were calculated on the basis of averaged spectra.

Dynamics of Plant Spectral and Biological Properties
Figure 2a and b show the temporal dynamics of reflectance properties for Spring Barley (a) and Sugar Beet (b) separately for all combinations of water (WC) and nitrogen (N) treatments.Spectra for each crop type/N-/WC-treatment combination vary not only between single time steps but also in between different N-and WC-treatments due to the changes in a canopy structure and differences in pigment and water content, both caused by different availability of water and nutrients.It is obvious from Figure 2 that spectral response significantly changes over time, however, to a different degree for different treatments and for different crop types.As pointed out by [45,46], spectral domains of strong chlorophyll (blue and red) and water (water bands) absorption are more sensitive to low concentration of chlorophyll and water, while others (green and regions between water bands) are more sensitive to high water and N-concentration.From the reflectance spectra in Figure 2, where we see both, a clear temporal development that seems to be dominated by the plant growth and therefore LAI dynamics, and (to a much smaller extent) differences under different irrigation conditions, we would expect at least some potential for the extraction of LAI and WC from narrow band SI's.In summary, plant spectral reflectance as well as the temporal development of LAI and crop height responds significantly to differences in water and nutrient availability.Even though reflectance sensitivity is not wavelength independent and varies within the entire spectral range, this behavior suggests the possibility of an inverse retrieval of crop characteristics and controlling water and nutrient conditions.This potential will be explored in the next sub-sections.

Dimensionality of the Spectral Data
The spectral reflectance properties presented in Figure 2 also illustrate-while containing some information on LAI and WC-that large spectral ranges vary similarly with time and in response to nutrient and water conditions [47].Thus, many of the spectral bands are inadequate or redundant in their information content when retrieving plant and soil characteristics such as LAI, chlorophyll or leaf and soil water content [48].In order to get estimates of the dimensionality of our experimental data set, which is an important guidance on the minimum band/index number necessary for any retrieval algorithm-both linear and non-linear, local and non-local methods are used in the following.

Linear Methods:
First, the complete spectral information with a number of 2,151 bands was used in a correlation analysis.Correlation coefficients of each band combination were calculated and the frequency distribution of their absolute values is presented in Figure 4 (left).It can be seen that the majority of values (~80%) are correlated to each other with absolute values greater than 0.95, indicating the redundancy in the information content.This is also highlighted by the results of a principal component analysis (PCA) over the whole set of bands, giving a first principal component (PC) with a loading of 99.4% and a second and third PC explaining no more than 0.38% and 0.21% of total variance, respectively (without figure).
Analyzing the correlation between selected spectral indices (see Table 1), it can be seen, that the proportion of pairs with correlation showing absolute values close to 1 is reduced, however the PCA analysis (without figure) shows an even higher loading of the first PC (99.99%).The correlation and the PCA-analysis for both, the individual spectral bands and the pre-selected spectral indices, suggest only a very low number of variables (here 1-3, dependent on the threshold value) needed to represent the spectral database information.This, however, only holds as long as there is no significant non-linearity present in the data that is not implicitly included by the use of complex spectral indices and that cannot be recovered by linear methods.

Non-linear Methods:
In contrast to the traditional linear techniques, nonlinear methods have the ability to deal with complex nonlinear data.Previous studies have shown that nonlinear techniques outperform their linear counterparts on complex artificial tasks.Examples are for instance the Swiss Roll dataset comprising a set of points that lie on a spiral-like two-dimensional manifold within a three-dimensional space [49,50].These methods aim at finding the intrinsic dimensionality (d) of a data set, being the minimum number of parameters that is necessary in order to account for all information contained within.
Techniques for intrinsic dimensionality estimation can be subdivided into two main groups: (1) estimators based on the analysis of local properties of the data and (2) estimators based on the analysis of global properties of the data.Local methods used in this analysis comprise the correlation dimension estimator (CorrDim), the nearest neighbor dimension estimator (NNDim), and the maximum likelihood estimator (MaxLike).Global estimators additional to the PCA considered in this analysis are the packing number estimator (PackNum), and the geodesic minimum spanning tree estimator (GMST).These methods are briefly introduced in the Appendix, but readers are referred to the "Matlab Toolbox for Dimensionality Reduction drtoolbox" [51][52][53], for a more comprehensive discussion and comparison.
Table 3 compares these different methods for both the complete spectral dataset including all spectral bands and the dataset given by the spectral indices listed in Table 1.The results show similar values when comparing the estimated intrinsic dimensionalities for the two data sets, suggesting that the much smaller (N = 28 vs. N = 2151) set of indices might well be able to represent the entire information content of the original dataset of all spectral bands.Comparing d estimations from different techniques, they show some more pronounced deviations indicating the different foci of the methods with respect to the non-linear characteristics of the data.However, the estimated intrinsic dimensionality of both datasets ranges around 3-6 (PCA, GMST) and does not exceed a value of 8 (MAXLike).Table 3.Estimated intrinsic dimensionalities (d) of the complete spectral dataset including all spectral bands and of the dataset given by the pre-selected spectral indices.(CorrDim-correlation dimension estimator, NNDim-nearest neighbor dimension estimator, MaxLike-maximum likelihood estimator, PackNum-packing number estimator, GMSTgeodesic minimum spanning tree estimator, PCA-principal component analysis).These values give some guidance on the maximum number of SI's as independent variable that will be needed in any type of regression or non-linear analysis.

Retrieval of LAI and Soil Moisture by Linear Regression
The intrinsic dimensionality estimation as described in the previous section, suggest that most of the information content of the original spectral dataset can be well described by a set of 28 SI's (Table 1).We therefore restrict any further analysis to this reduced data set.In a first step, the simple linear dependencies between LAI and soil WC and single SI's are investigated.Standard linear regression was applied to extract the single SI with the maximum explanatory potential.Calculations were done for three variants: (1) using both crop types in one dataset and (2, 3) treating each crop species spring barley (2) and sugar beet (3) separately.Table 4 lists the coefficients of determination (R 2 ) when predicting LAI and WC from a single SI for the three variants (1-3), as described above.Also given for each of the variants are the SI's (according to Table 1) for which the maximum R 2 is achieved.
The results in Table 4 present a clear pattern: First, the R 2 values are relatively low in general with predictions being generally better for WC compared to LAI (maximum R 2 of 0.60 vs. 0.42) indicating that a single spectral index will not be sufficient to describe plant/soil variables over the range of conditions presented here by the individual nutrients and water treatments.Second, simple linear regression shows a better performance for both LAI and WC when applied to individual crop types rather than treating crops in a common application.In practice, this would mean that knowing crop species beforehand will significantly improve such retrieval, but will make the application of land use classification schemes necessary beforehand, thereby introducing potential misclassification error.The results from simple linear regressions (Table 4) as well as intrinsic dimensionality analysis of the spectral dataset suggest that multiple explanatory variables (3)(4) should be able to more adequately explain the dynamic dependencies of WC, LAI and the vegetation spectral responses (SI's).Table 4. Coefficients of determination (R 2 ) and root mean square error (rmse, in brackets) of simple univariate and multiple (three-fold) linear regressions between LAI/WC and pre-selected SI's (see Table 1).Values listed are the maximum R² values and the corresponding root mean square error (rmse) when varying through all 28 indices for univariate linear regression and through all combinations of SI's (one from each of the three functional groups, see Table 1) for multiple linear regression.This analysis has been performed using data for "both crops" (Spring Barley and Sugar Beet) together as well as separately ("only barley" and "only sugar").The columns "SI" refers to the spectral index (see Table 1) achieving the best fit.In a next step, we therefore tested multiple linear regression using three explanatory SI's-one from each functional group (see Table 1)-leading to a scenario of 560 possible combinations of three-parameter explanatory inputs.Calculations were done in a similar way using three variants (1) for the entire dataset incorporating both crops together and for each species separately: Spring Barley (2) and Sugar Beet (3).Results are shown in Table 4 (right), demonstrating a significant improvement in the retrieval of WC and LAI from a combination of SI's for both the calculated R 2 -values as well as rmse.Interestingly, in four of the six applications, the SI's that produced the best fit within the univariate regression analysis was not among the best combination of SI's in the multiple regression approach, indicating that some of the individual SI's might be better able to represented the combined plant response to water and nutrient stress, but a combination of each multiple SI's are more selective to a certain aspect of plant stress response and perform better.In particular the retrieval of LAI (for both crops) has gained drastically from inclusion of additional SI's confirming findings of e.g., [19] that LAI is not only influences by leaf chlorophyll concentration but also by other biochemical and crop/leaf structural factors that are possibly reflected in other SI's.

Multivariate Linear Regression
While the results from dimensionality analysis suggested the use of 3-4 independent variable SI's best representing the spectral information, we have investigated the consideration of further SI's into the modeling procedure as shown in Figure 5.While considerable improvements in the prediction of LAI is achieved by an increase to three variables, only small increases in the model performance are observed by including additional SI's as explanatory variables.For the WC retrieval, already one SI achieves relatively high R 2 -values that only slightly improved by adding further SI's.
The overall performance of the multiple-(3-SI)-regression-retrieval approach with R²-values of 0.65 for WC (both crops), and 0.67 and 0.33 for LAI (barley and sugar beet, respectively) compare very much to the upper limit of results that have been summarized by [19] showing a range of 0.05-0.66for the retrieval of LAI.However, it needs to be considered, that our experiments have been performed under local, partly idealized glasshouse conditions, with direct control on illumination conditions and without the disturbance of atmospheric effects.

Retrieval of LAI and Soil Moisture by Non-Linear Regression Trees
In addition to the multiple-linear regression approach, we have investigated the possible improvement and benefit of using non-linear techniques in the retrieval algorithm.A large variety of different techniques have been developed over the past decades, many of them emerging from the area of artificial intelligence, data mining and machine learning.Among those neuronal networks [54], fuzzy rules based approaches [55], support vector machines [56], random trees [57] and classification and regression trees (CART) [58] are the most popular.
The CART method was chosen here as it represents a simple, comprehensible, easily to implement and widely used method.It is a nonparametric modeling approach that can describe the responses of a dependent from a set of independent variables.In our application, CART recursively partition the set of SI's to find increasingly homogeneous subsets based on SI's splitting criteria using variance minimizing algorithms.The (dependent) LAI and WC data are partitioned into a series of descending left and right child nodes derived from parent nodes.Once the partitioning has ceased, the child nodes are designated as terminal nodes.For CART calibration, a cost function is defined as the product of the estimated probability of a node times the cost of a node, that is the average squared error over the observations in that node, summed over all terminal nodes [58].The "optimal" size of the tree has been extracted using 10-fold cross-validation and choosing the smallest tree that is within one standard error of the minimum-cost sub-tree.As CART, as well as all the other non-linear techniques mentioned above, show the ability to "overfit" the available data, we have compared multiple-linear regression (as presented in Section 3.3) and CART in a cross-validation framework using a leave-one-out strategy.
Results shown in Table 5 using the combination of SI's achieving the best results in the multiple-linear-regression analysis presented before (Table 4) clearly demonstrate the difficulties of using non-linear techniques.Using the standard settings of the CART tree building procedure, only the retrieval of LAI for both crops shows superior results to multiple-linear regression, while "fine-tuning" of individual settings in the tree-building has been shown to improve the validation results for one of the dependent variable and crop combinations, the decrease at the same time for the others.Table 5. Coefficients of determination (R 2 ) and root mean square errors (rmse) for multiple linear regression and CART between LAI and WC using the best combinations of pre-selected SI's from the analysis in Section 3.3 within a cross-validation (leave-one-out) framework.Results are representing cross-validation results for both crops (Spring Barley and Sugar Beet) together and separately for each one.While we do not want to generalize these finding to all other non-linear techniques listed above, we here want state that at least using the CART approach for the retrieval of LAI and WC for our experimental data did not show the potential for any improvement in the retrieval process.In consequence, this also means that most of the non-linearity that might be inherent in the LAI/WC to spectral reflectance dependency is already "included" in the non-linear combination of individual and selective narrow spectral bands of the SI's.

Summary and Conclusions
Time series of hyper-spectral reflectance data (400-2,500 nm) for two crop species-sugar beet and spring barley-over a vegetation period clearly demonstrated a significant response of plant reflectance characteristics to variations in WC and nutrient conditions.Based on these findings, we investigated the potential of retrieving plant LAI and soil WC from spectral reflectance data.First, the full range of spectral bands (N = 2151), as well as a reduced spectral data-base consisting of 28 spectral indices extracted from literature review (see Table 1; N = 28), have been analyzed for their information in a system analysis sense.Using linear and non-linear techniques, the intrinsic dimensionality of both datasets suggested a number of 3-4 but not more than 8 bands/indices to fully represent the information contained in the spectral data.Results also suggest, that the full band information is well represented by the set of 28 SI's that have been divided into three functional groups: (i) greenness/biomass/vegetation structure, (ii) pigment or light use efficiency, and (iii) water related indices (see Table 1).
Restricting the further analysis to the use of the set of SI's, we subsequently analyzed the potential for LAI and WC retrieval from the SI's using single and multiple linear regression.It was shown that including three indices in the retrieval algorithm-one index from each functional group-could significantly improve the results.In general, the retrieval of LAI and WC was significantly higher when only considering single species.R² values range from 0.65 for WC retrieval for both sugar beet and spring barley when treated separately, to 0.67 and 0.33 for LAI retrieval.These results are very much in line with a summary statistics of LAI retrieval results reviewed by [19] given values in the range of 0.05-0.66,however we have to consider that here we have been working under an almost ideal glasshouse situation with cultivated plants under optimal and controlled irrigation, fertilization and illumination conditions and were also able to achieve constant light source sensor configuration thereby avoiding bidirectional effects.
While the results achieved are promising, we still need to extend our analysis considering the following aspects: While the best retrieval results have been achieved considering only single crop species, this would require the application of classification schemes beforehand when applying glasshouse calibrated algorithms to "real world" field conditions.This is currently under investigation using field spectral measurements at selected experimental test plot of the Bad Lauchstädt long-term fertilization experiment and will be presented in near future.
The potential of non-linear retrieval algorithms needs to be extended to a more thorough analysis considering a larger variety of methods in addition to the CART approach presented here.However, first results presented here within a cross-validation framework suggest that most of the non-linearity in the LAI/WC to spectral reflectance relationship is already covered by the non-linearity of the SI's applied, but needs to be further tested.Our aim in publishing the here presented data set to the remote sensing community will hopefully support this requirement.
While our analysis has mainly built on the use of SI's, at least for the retrieval of LAI, the inversion of a leaf/canopy reflectance model such as PROSPECT/SAIL [36,46,59] might have the potential to improve the results, but again while being beyond the scope of this paper, this aspect needs to be investigated in more detail and is part of our current research activities that will certainly be presented in the near future.neighborhood graph defined on the dataset.The length function of the GMST is the sum of the Euclidean distances corresponding to all edges in the geodesic minimum spanning tree.

Figure 2 .
Figure 2. Temporal dynamics of spectral properties of Spring Barley (a) and Sugar Beet (b), plotted with regards to applied WC and N treatments (indicated in an upper right corner of all plots), during the BL 2008 experiment.

Figure 3 .
Figure 3. Temporal development of leaf area index (a) and of vegetation height (b) for Spring Barley (o) and Sugar Beet (+), plotted with regards to applied WC and N treatments (indicated in an upper right corner of all plots), BL 2008 Experiment.

Figure 4 .
Figure 4. Histogram of the absolute values of pairwise correlation between individual bands (left) and vegetation indices (right) that are listed in Table 1.N indicates the number of spectral bands or vegetation indices considered.The total number of pairs can be calculated as N * (N − 1) /2.

Figure 5 .
Figure 5. R 2 improvement within multiple linear regression analysis due to increase of explanatory variables for LAI (continuous line) and theta (discontinuous line) for both crops species together (top), (2) for barley (middle), and (3) for sugar beet (bottom).

Table 1 .
List of SI's to describe plant ecophysiological parameters reviewed from the literature and used throughout this paper.SI's are grouped into three general classes showing response to (a) greenness, biomass and structure of the vegetation, (b) light use efficiency, and (c) vegetation and leaf water content.R xxx indicates the reflectance at a specific wavelength (nm).

Table 2 .
Time schedule of the Bad Lauchstädt Experiment 2008.Black vertical pointers indicate in which part of the month particular events took place.