Sources of variability in canopy reflectance and the convergent properties of plants

Discrimination of Canopy Structural Types in the Sierra Nevada Mountains in Central California

Accurate information about ecosystem structure and biogeochemical properties is essential to providing better estimates ecosystem functioning. Airborne LiDAR (light detection and ranging) is the most accurate way to retrieve canopy structure. However, accurately obtaining both biogeochemical traits and structure parameters requires concurrent measurements from imaging spectrometers and LiDARs. Our main objective was to evaluate the use of imaging spectroscopy (IS) to provide vegetation structural information. We developed models to estimate structural variables (i.e., biomass, height, vegetation heterogeneity and clumping) using IS data with a random forests model from three forest ecosystems (i.e., an oak-pine low elevation savanna, a mixed conifer/broadleaf mid-elevation forest, and a high-elevation montane conifer forest) in the Sierra Nevada Mountains, California. We developed and tested general models to estimate the four structural variables with accuracies greater than 75%, for the structurally and ecologically different forest sites, demonstrating their applicability to a diverse range of forest ecosystems. The model R2 for each structural variable was least in the conifer/broadleaf forest than either the low elevation savanna or the montane conifer forest. We then used the structural variables we derived to discriminate site-specific, ecologically meaningful descriptions of canopy structural types (CST). Our CST results demonstrate how IS data can be used to create comprehensive and easily interpretable maps of forest structural types that capture their major structural features and trends across different vegetation types in the Sierra Nevada Mountains. The mixed conifer/broadleaf forest and montane conifer forest had the most complex structures, containing six and five CSTs respectively. The identification of CSTs within a site allowed us to better identify the main drivers of structural variability in each ecosystem. CSTs in open savanna were driven mainly by differences in vegetation cover; in the mid-elevation mixed forest, by the combination of biomass and canopy height; and in the montane conifer forest, by vegetation heterogeneity and clumping.

Radiative transfer modelling reveals why canopy reflectance follows function

Optical remote sensing is potentially highly informative to track Earth's plant functional diversity. Yet, causal explanations of how and why plant functioning is expressed in canopy reflectance remain limited. Variation in canopy reflectance can be described by radiative transfer models (here PROSAIL) that incorporate plant traits affecting light transmission in canopies. To establish causal links between canopy reflectance and plant functioning, we investigate how two plant functional schemes, i.e. the Leaf Economic Spectrum (LES) and CSR plant strategies, are related to traits with relevance to reflectance. These traits indeed related to both functional schemes, whereas only traits describing leaf properties correlated with the LES. In contrast, traits related to canopy structure showed no correlation to the LES, but to CSR strategies, as the latter integrates both plant economics and size traits, rather than solely leaf economics. Multiple optically relevant traits featured comparable or higher correspondence to the CSR space than those traits originally used to allocate CSR scores. This evidences that plant functions and strategies are directly expressed in reflectance and entails that canopy 'reflectance follows function'. This opens up new possibilities to understand differences in plant functioning and to harness optical remote sensing data for monitoring Earth´s functional diversity.

A Random Forest Machine Learning Approach for the Retrieval of Leaf Chlorophyll Content in Wheat

Developing rapid and non-destructive methods for chlorophyll estimation over large spatial areas is a topic of much interest, as it would provide an indirect measure of plant photosynthetic response, be useful in monitoring soil nitrogen content, and offer the capacity to assess vegetation structural and functional dynamics. Traditional methods of direct tissue analysis or the use of handheld meters, are not able to capture chlorophyll variability at anything beyond point scales, so are not particularly useful for informing decisions on plant health and status at the field scale. Examining the spectral response of plants via remote sensing has shown much promise as a means to capture variations in vegetation properties, while offering a non-destructive and scalable approach to monitoring. However, determining the optimum combination of spectra or spectral indices to inform plant response remains an active area of investigation. Here, we explore the use of a machine learning approach to enhance the estimation of leaf chlorophyll (Chlt), defined as the sum of chlorophyll a and b, from spectral reflectance data. Using an ASD FieldSpec 4 Hi-Res spectroradiometer, 2700 individual leaf hyperspectral reflectance measurements were acquired from wheat plants grown across a gradient of soil salinity and nutrient levels in a greenhouse experiment. The extractable Chlt was determined from laboratory analysis of 270 collocated samples, each composed of three leaf discs. A random forest regression algorithm was trained against these data, with input predictors based upon (1) reflectance values from 2102 bands across the 400–2500 nm spectral range; and (2) 45 established vegetation indices. As a benchmark, a standard univariate regression analysis was performed to model the relationship between measured Chlt and the selected vegetation indices. Results show that the root mean square error (RMSE) was significantly reduced when using the machine learning approach compared to standard linear regression. When exploiting the entire spectral range of individual bands as input variables, the random forest estimated Chlt with an RMSE of 5.49 µg·cm−2 and an R2 of 0.89. Model accuracy was improved when using vegetation indices as input variables, producing an RMSE ranging from 3.62 to 3.91 µg·cm−2, depending on the particular combination of indices selected. In further analysis, input predictors were ranked according to their importance level, and a step-wise reduction in the number of input features (from 45 down to 7) was performed. Implementing this resulted in no significant effect on the RMSE, and showed that much the same prediction accuracy could be obtained by a smaller subset of indices. Importantly, the random forest regression approach identified many important variables that were not good predictors according to their linear regression statistics. Overall, the research illustrates the promise in using established vegetation indices as input variables in a machine learning approach for the enhanced estimation of Chlt from hyperspectral data.

Analysis of Sentinel-2 and RapidEye for Retrieval of Leaf Area Index in a Saltmarsh Using a Radiative Transfer Model

The Sentinel satellite fleet of the Copernicus Programme offers new potential to map and monitor plant traits at fine spatial and temporal resolutions. Among these traits, leaf area index (LAI) is a crucial indicator of vegetation growth and an essential variable in biodiversity studies. Numerous studies have shown that the radiative transfer approach has been a successful method to retrieve LAI from remote-sensing data. However, the suitability and adaptability of this approach largely depend on the type of remote-sensing data, vegetation cover and the ecosystem studied. Saltmarshes are important wetland ecosystems threatened by sea level rise among other human- and animal-induced changes. Therefore, monitoring their vegetation status is crucial for their conservation, yet few LAI assessments exist for these ecosystems. In this study, the retrieval of LAI in a saltmarsh ecosystem is examined using Sentinel-2 and RapidEye data through inversion of the PROSAIL radiative transfer model. Field measurements of LAI and some other plant traits were obtained during two succeeding field campaigns in July 2015 and 2016 on the saltmarsh of Schiermonnikoog, a barrier island of the Netherlands. RapidEye (2015) and Sentinel-2 (2016) data were acquired concurrent to the time of the field campaigns. The broadly employed PROSAIL model was inverted using two look-up tables (LUTs) generated in the spectral band’s settings of the two sensors and in which each contained 500,000 records. Different solutions from the LUTs, as well as, different Sentinel-2 spectral subsets were considered to examine the LAI retrieval. Our results showed that generally the LAI retrieved from Sentinel-2 had higher accuracy compared to RapidEye-retrieved LAI. Utilising the mean of the first 10 best solutions from the LUTs resulted in higher R2 (0.51 and 0.59) and lower normalised root means square error (NRMSE) (0.24 and 0.16) for both RapidEye and Sentinel-2 data respectively. Among different Sentinel-2 spectral subsets, the one comprised of the four near-infrared (NIR) and shortwave infrared (SWIR) spectral bands resulted in higher estimation accuracy (R2 = 0.44, NRMSE = 0.21) in comparison to using other studied spectral subsets. The results demonstrated the feasibility of broadband multispectral sensors, particularly Sentinel-2 for retrieval of LAI in the saltmarsh ecosystem via inversion of PROSAIL. Our results highlight the importance of proper parameterisation of radiative transfer models and capacity of Sentinel-2 spectral range and resolution, with impending high-quality global observation aptitude, for retrieval of plant traits at a global scale.

Assessing Across-Scale Optical Diversity and Productivity Relationships in Grasslands of the Italian Alps

The linearity and scale-dependency of ecosystem biodiversity and productivity relationships (BPRs) have been under intense debate. In a changing climate, monitoring BPRs within and across different ecosystem types is crucial, and novel remote sensing tools such as the Sentinel-2 (S2) may be adopted to retrieve ecosystem diversity information and to investigate optical diversity and productivity patterns. But are the S2 spectral and spatial resolutions suitable to detect relationships between optical diversity and productivity? In this study, we implemented an integrated analysis of spatial patterns of grassland productivity and optical diversity using optical remote sensing and Eddy Covariance data. Across-scale optical diversity and ecosystem productivity patterns were analyzed for different grassland associations with a wide range of productivity. Using airborne optical data to simulate S2, we provided empirical evidence that the best optical proxies of ecosystem productivity were linearly correlated with optical diversity. Correlation analysis at increasing pixel sizes proved an evident scale-dependency of the relationships between optical diversity and productivity. The results indicate the strong potential of S2 for future large-scale assessment of across-ecosystem dynamics at upper levels of observation.

Combining biophysical parameters, spectral indices and multivariate hyperspectral models for estimating yield and water productivity of spring wheat across different agronomic practices

Manipulating plant densities under different irrigation rates can have a significant impact on grain yield and water use efficiency by exerting positive or negative effects on ET. Whereas traditional spectral reflectance indices (SRIs) have been used to assess biophysical parameters and yield, the potential of multivariate models has little been investigated to estimate these parameters under multiple agronomic practices. Therefore, both simple indices and multivariate models (partial least square regression (PLSR) and support vector machines (SVR)) obtained from hyperspectral reflectance data were compared for their applicability for assessing the biophysical parameters in a field experiment involving different combinations of three irrigation rates (1.00, 0.75, and 0.50 ET) and five plant densities (D₁: 150, D₂: 250, D₃: 350, D₄: 450, and D₅: 550 seeds m^-2) in order to improve productivity and water use efficiency of wheat. Results show that the highest values for green leaf area, aboveground biomass, and grain yield were obtained from the combination of D₃ or D₄ with 1.00 ET, while the combination of 0.75 ET and D₃ was the best treatment for achieving the highest values for water use efficiency. Wheat yield response factor (ky) was acceptable when the 0.75 ET was combined with D₂, D₃, or D₄ or when the 0.50 ET was combined with D₂ or D₃, as the ky values of these combinations were less than or around one. The production function indicated that about 75% grain yield variation could be attributed to the variation in seasonal ET. Results also show that the performance of the SRIs fluctuated when regressions were analyzed for each irrigation rate or plant density specifically, or when the data of all irrigation rates or plant densities were combined. Most of the SRIs failed to assess biophysical parameters under specific irrigation rates and some specific plant densities, but performance improved substantially for combined data of irrigation rates and some specific plant densities. PLSR and SVR produced more accurate estimations of biophysical parameters than SRIs under specific irrigation rates and plant densities. In conclusion, hyperspectral data are useful for predicting and monitoring yield and water productivity of spring wheat across multiple agronomic practices.

Spectral Heterogeneity Predicts Local-Scale Gamma and Beta Diversity of Mesic Grasslands

Plant species diversity is an important metric of ecosystem functioning, but field assessments of diversity are constrained in number and spatial extent by labor and other expenses. We tested the utility of using spatial heterogeneity in the remotely-sensed reflectance spectrum of grassland canopies to model both spatial turnover in species composition and abundances (β diversity) and species diversity at aggregate spatial scales (γ diversity). Shannon indices of γ and β diversity were calculated from field measurements of the number and relative abundances of plant species at each of two spatial grains (0.45 m2 and 35.2 m2) in mesic grasslands in central Texas, USA. Spectral signatures of reflected radiation at each grain were measured from ground-level or an unmanned aerial vehicle (UAV). Partial least squares regression (PLSR) models explained 59–85% of variance in γ diversity and 68–79% of variance in β diversity using spatial heterogeneity in canopy optical properties. Variation in both γ and β diversity were associated most strongly with heterogeneity in reflectance in blue (350–370 nm), red (660–770 nm), and near infrared (810–1050 nm) wavebands. Modeled diversity was more sensitive by a factor of three to a given level of spectral heterogeneity when derived from data collected at the small than larger spatial grain. As estimated from calibrated PLSR models, β diversity was greater, but γ diversity was smaller for restored grassland on a lowland clay than upland silty clay soil. Both γ and β diversity of grassland can be modeled by using spatial heterogeneity in vegetation optical properties provided that the grain of reflectance measurements is conserved.

Multi-Crop Green LAI Estimation with a New Simple Sentinel-2 LAI Index (SeLI)

The spatial quantification of green leaf area index (LAIgreen), the total green photosynthetically active leaf area per ground area, is a crucial biophysical variable for agroecosystem monitoring. The Sentinel-2 mission is with (1) a temporal resolution lower than a week, (2) a spatial resolution of up to 10 m, and (3) narrow bands in the red and red-edge region, a highly promising mission for agricultural monitoring. The aim of this work is to define an easy implementable LAIgreen index for the Sentinel-2 mission. Two large and independent multi-crop datasets of in situ collected LAIgreen measurements were used. Commonly used LAIgreen indices applied on the Sentinel-2 10 m × 10 m pixel resulted in a validation R² lower than 0.6. By calculating all Sentinel-2 band combinations to identify high correlation and physical basis with LAIgreen, the new Sentinel-2 LAIgreen Index (SeLI) was defined. SeLI is a normalized index that uses the 705 nm and 865 nm centered bands, exploiting the red-edge region for low-saturating absorption sensitivity to photosynthetic vegetation. A R² of 0.708 (root mean squared error (RMSE) = 0.67) and a R² of 0.732 (RMSE = 0.69) were obtained with a linear fitting for the calibration and validation datasets, respectively, outperforming established indices. Sentinel-2 LAIgreen maps are presented.

Improving Forest Aboveground Biomass (AGB) Estimation by Incorporating Crown Density and Using Landsat 8 OLI Images of a Subtropical Forest in Western Hunan in Central China

Forest aboveground biomass (AGB) estimation modeling based on remote sensing is an important method for large-scale biomass estimation; the accuracy of the estimation models has been a topic of broad and current interest. In this study, we used permanent sample plot data and Landsat 8 Operational Land Imager (OLI) images of western Hunan. Remote-sensing-based models were developed for different vegetation types, and different crown density classes were incorporated. The linear model, linear dummy variable model, and linear mixed-effects model were used to determine the most effective and accurate method for remote-sensing-based AGB estimation. The results show that the adjusted coefficient of determination (R2adj) and root mean square error (RMSE) of the linear dummy model and linear mixed-effects model were significantly better than those of the linear model; the R2adj increased more than 0.16 and the RMSE decreased more than 2.12 for each vegetation type, and the F-test also showed significant differences between the linear model and linear dummy variable model and between the linear model and linear mixed-effects model. The accuracies of the AGB estimations of the linear dummy variable model and the linear mixed-effects model were significantly better than those of linear model in the thin and dense crown density classes. There were no significant differences in the AGB estimation performance between the linear dummy variable model and linear mixed-effects model; these two models were more flexible and more suitable than the linear model for remote-sensing-based AGB estimation. The results of this study provide a new approach for solving the low-accuracy estimations of linear models.

Reproductive success of soybean (Glycine max L. Merril) cultivars and exotic lines under high daytime temperature

The objectives were to (a) quantify the effects of high daytime temperature (HDT) from gametogenesis to full bloom on photosynthesis and pod set in soybean (Glycine max L. Merril) genotypes and (b) assess the relationships among photosynthesis, cardinal temperatures for pollen germination, in vitro pollen germination percentage, canopy reflectance, and pod-set percentage. Three field experiments were conducted, and Experiment I had HDT between gametogenesis and full bloom (36.5°C to 38.6°C) compared with Experiments II and III (29.5°C to 31.6°C; optimum temperature). HDT decreased photosynthesis (22%) and pod-set percent (11%) compared with Experiment III. Cultivars had higher photosynthesis and pod-set percent than plant introduction (PI) lines. The cultivars (i.e., IA3023 and KS4694) and PI lines (i.e., PI393540 and PI588026A) were HDT tolerant and susceptible, respectively. The decreased pod-set percentage in susceptible genotypes (PI lines) was associated with pollen characteristics. Significant positive (r2  ≥ 0.67) association between photosynthesis, cardinal temperatures for pollen germination (Topt and Tmax ) with pod-set percentage was observed. However, a negative (r2  ≥ -0.43) association between photosynthesis and pod set with canopy reflectance at visible spectrum was observed. In vitro pollen germination and canopy reflectance at visible spectrum can be used as a high-throughput phenotypic tool for breeding HDT-tolerant genotypes.

Quantifying Vegetation Biophysical Variables from Imaging Spectroscopy Data: A Review on Retrieval Methods

An unprecedented spectroscopic data stream will soon become available with forthcoming Earth-observing satellite missions equipped with imaging spectroradiometers. This data stream will open up a vast array of opportunities to quantify a diversity of biochemical and structural vegetation properties. The processing requirements for such large data streams require reliable retrieval techniques enabling the spatiotemporally explicit quantification of biophysical variables. With the aim of preparing for this new era of Earth observation, this review summarizes the state-of-the-art retrieval methods that have been applied in experimental imaging spectroscopy studies inferring all kinds of vegetation biophysical variables. Identified retrieval methods are categorized into: (1) parametric regression, including vegetation indices, shape indices and spectral transformations; (2) nonparametric regression, including linear and nonlinear machine learning regression algorithms; (3) physically based, including inversion of radiative transfer models (RTMs) using numerical optimization and look-up table approaches; and (4) hybrid regression methods, which combine RTM simulations with machine learning regression methods. For each of these categories, an overview of widely applied methods with application to mapping vegetation properties is given. In view of processing imaging spectroscopy data, a critical aspect involves the challenge of dealing with spectral multicollinearity. The ability to provide robust estimates, retrieval uncertainties and acceptable retrieval processing speed are other important aspects in view of operational processing. Recommendations towards new-generation spectroscopy-based processing chains for operational production of biophysical variables are given.

The CCB-ID approach to tree species mapping with airborne imaging spectroscopy

BACKGROUND

Biogeographers assess how species distributions and abundances affect the structure, function, and composition of ecosystems. Yet we face a major challenge: it is difficult to precisely map species across landscapes. Novel Earth observations could overcome this challenge for vegetation mapping. Airborne imaging spectrometers measure plant functional traits at high resolution, and these measurements can be used to identify tree species. In this paper, I describe a trait-based approach to species identification with imaging spectroscopy, the Center for Conservation Biology species identification (CCB-ID) method, which was developed as part of an ecological data science evaluation competition.

METHODS

These methods were developed using airborne imaging spectroscopy data from the National Ecological Observatory Network (NEON). CCB-ID classified tree species using trait-based reflectance variation and decision tree-based machine learning models, approximating a morphological trait and dichotomous key method inspired by botanical classification. First, outliers were removed using a spectral variance threshold. The remaining samples were transformed using principal components analysis (PCA) and resampled to reduce common species biases. Gradient boosting and random forest classifiers were trained using the transformed and resampled feature data. Prediction probabilities were calibrated using sigmoid regression, and sample-scale predictions were averaged to the crown scale.

RESULTS

CCB-ID received a rank-1 accuracy score of 0.919, and a cross-entropy cost score of 0.447 on the competition test data. Accuracy and specificity scores were high for all species, but precision and recall scores varied for rare species. PCA transformation improved accuracy scores compared to models trained using reflectance data, but outlier removal and data resampling exacerbated class imbalance problems.

DISCUSSION

CCB-ID accurately classified tree species using NEON data, reporting the best scores among participants. However, it failed to overcome several species mapping challenges like precisely identifying rare species. Key takeaways include (1) selecting models using metrics beyond accuracy (e.g., recall) could improve rare species predictions, (2) within-genus trait variation may drive spectral separability, precluding efforts to distinguish between functionally convergent species, (3) outlier removal and data resampling can exacerbate class imbalance problems, and should be carefully implemented, (4) PCA transformation greatly improved model results, and (5) targeted feature selection could further improve species classification models. CCB-ID is open source, designed for use with NEON data, and available to support species mapping efforts.

Landscape variation in canopy nitrogen and carbon assimilation in a temperate mixed forest

Canopy nitrogen (N) is a key factor regulating carbon cycling in forest ecosystems through linkages among foliar N and photosynthesis, decomposition, and N cycling. This analysis examined landscape variation in canopy nitrogen and carbon assimilation in a temperate mixed forest surrounding Harvard Forest in central Massachusetts, USA by integration of canopy nitrogen mapping with ecosystem modeling, and spatial data from soils, stand characteristics and disturbance history. Canopy %N was mapped using high spectral resolution remote sensing from NASA's AVIRIS (Airborne Visible/Infrared Imaging Spectrometer) instrument and linked to an ecosystem model, PnET-II, to estimate gross primary productivity (GPP). Predicted GPP was validated with estimates derived from eddy covariance towers. Estimated canopy %N ranged from 0.5 to 2.9% with a mean of 1.75% across the study region. Predicted GPP ranged from 797 to 1622 g C m^-2 year^-1 with a mean of 1324 g C m^-2 year^-1. The prediction that spatial patterns in forest growth are associated with spatial patterns in estimated canopy %N was supported by a strong, positive relationship between field-measured canopy %N and aboveground net primary production. Estimated canopy %N and GPP were related to forest composition, land-use history, and soil drainage. At the landscape scale, PnET-II GPP was compared with predicted GPP from the BigFoot project and from NASA's MODIS (Moderate Resolution Imaging Spectroradiometer) data products. Estimated canopy %N explained much of the difference between MODIS GPP and PnET-II GPP, suggesting that global MODIS GPP estimates may be improved if broad-scale estimates of foliar N were available.

Using imaging spectroscopy to detect variation in terrestrial ecosystem productivity across a water-stressed landscape

A central challenge to understanding how climate anomalies, such as drought and heatwaves, impact the terrestrial carbon cycle, is quantification and scaling of spatial and temporal variation in ecosystem gross primary productivity (GPP). Existing empirical and model-based satellite broadband spectra-based products have been shown to miss critical variation in GPP. Here, we evaluate the potential of high spectral resolution (10 nm) shortwave (400-2,500 nm) imagery to better detect spatial and temporal variations in GPP across a range of ecosystems, including forests, grassland-savannas, wetlands, and shrublands in a water-stressed region. Estimates of GPP from eddy covariance observations were compared against airborne hyperspectral imagery, collected across California during the 2013-2014 HyspIRI airborne preparatory campaign. Observations from 19 flux towers across 23 flight campaigns (102 total image-flux tower pairs) showed GPP to be strongly correlated to a suite of spectral wavelengths and band ratios associated with foliar physiology and chemistry. A partial least squares regression (PLSR) modeling approach was then used to predict GPP with higher validation accuracy (adjusted R²  = 0.71) and low bias (0.04) compared to existing broadband approaches (e.g., adjusted R²  = 0.68 and bias = -5.71 with the Sims et al. model). Significant wavelengths contributing to the PLSR include those previously shown to coincide with Rubisco (wavelengths 1,680, 1,740, and 2,290 nm) and V_cmax (wavelengths 1,680, 1,722, 1,732, 1,760, and 2,300 nm). These results provide strong evidence that advances in satellite spectral resolution offer significant promise for improved satellite-based monitoring of GPP variability across a diverse range of terrestrial ecosystems.

Proximal Optical Sensors for Nitrogen Management of Vegetable Crops: A Review

Optimal nitrogen (N) management is essential for profitable vegetable crop production and to minimize N losses to the environment that are a consequence of an excessive N supply. Proximal optical sensors placed in contact with or close to the crop can provide a rapid assessment of a crop N status. Three types of proximal optical sensors (chlorophyll meters, canopy reflectance sensors, and fluorescence-based flavonols meters) for monitoring the crop N status of vegetable crops are reviewed, addressing practical caveats and sampling considerations and evaluating the practical use of these sensors for crop N management. Research over recent decades has shown strong relationships between optical sensor measurements, and different measures of crop N status and of yield of vegetable species. However, the availability of both: (a) Sufficiency values to assess crop N status and (b) algorithms to translate sensor measurements into N fertilizer recommendations are limited for vegetable crops. Optical sensors have potential for N management of vegetable crops. However, research should go beyond merely diagnosing crop N status. Research should now focus on the determination of practical fertilization recommendations. It is envisaged that the increasing environmental and societal pressure on sustainable crop N management will stimulate progress in this area.

Leaf dorsoventrality as a paramount factor determining spectral performance in field-grown wheat under contrasting water regimes

The effects of leaf dorsoventrality and its interaction with environmentally induced changes in the leaf spectral response are still poorly understood, particularly for isobilateral leaves. We investigated the spectral performance of 24 genotypes of field-grown durum wheat at two locations under both rainfed and irrigated conditions. Flag leaf reflectance spectra in the VIS-NIR-SWIR (visible-near-infrared-short-wave infrared) regions were recorded in the adaxial and abaxial leaf sides and at the canopy level, while traits providing information on water status and grain yield were evaluated. Moreover, leaf anatomical parameters were measured in a subset of five genotypes. The spectral traits studied were more affected by the leaf side than by the water regime. Leaf dorsoventral differences suggested higher accessory pigment content in the abaxial leaf side, while water regime differences were related to increased chlorophyll, nitrogen, and water contents in the leaves in the irrigated treatment. These variations were associated with anatomical changes. Additionally, leaf dorsoventral differences were less in the rainfed treatment, suggesting the existence of leaf-side-specific responses at the anatomical and biochemical level. Finally, the accuracy in yield prediction was enhanced when abaxial leaf spectra were employed. We concluded that the importance of dorsoventrality in spectral traits is paramount, even in isobilateral leaves.

Retrieval of canopy water content of different crop types with two new hyperspectral indices: Water Absorption Area Index and Depth Water Index

Crop canopy water content (CWC) is an essential indicator of the crop's physiological state. While a diverse range of vegetation indices have earlier been developed for the remote estimation of CWC, most of them are defined for specific crop types and areas, making them less universally applicable. We propose two new water content indices applicable to a wide variety of crop types, allowing to derive CWC maps at a large spatial scale. These indices were developed based on PROSAIL simulations and then optimized with an experimental dataset (SPARC03; Barrax, Spain). This dataset consists of water content and other biophysical variables for five common crop types (lucerne, corn, potato, sugar beet and onion) and corresponding top-of-canopy (TOC) reflectance spectra acquired by the hyperspectral HyMap airborne sensor. First, commonly used water content index formulations were analysed and validated for the variety of crops, overall resulting in a R² lower than 0.6. In an attempt to move towards more generically applicable indices, the two new CWC indices exploit the principal water absorption features in the near-infrared by using multiple bands sensitive to water content. We propose the Water Absorption Area Index (WAAI) as the difference between the area under the null water content of TOC reflectance (reference line) simulated with PROSAIL and the area under measured TOC reflectance between 911 and 1271 nm. We also propose the Depth Water Index (DWI), a simplified four-band index based on the spectral depths produced by the water absorption at 970 and 1200 nm and two reference bands. Both the WAAI and DWI outperform established indices in predicting CWC when applied to heterogeneous croplands, with a R² of 0.8 and 0.7, respectively, using an exponential fit. However, these indices did not perform well for species with a low fractional vegetation cover (< 30%). HyMap CWC maps calculated with both indices are shown for the Barrax region. The results confirmed the potential of using generically applicable indices for calculating CWC over a great variety of crops.

Biological processes dominate seasonality of remotely sensed canopy greenness in an Amazon evergreen forest

Satellite observations of Amazon forests show seasonal and interannual variations, but the underlying biological processes remain debated. Here we combined radiative transfer models (RTMs) with field observations of Amazon forest leaf and canopy characteristics to test three hypotheses for satellite-observed canopy reflectance seasonality: seasonal changes in leaf area index, in canopy-surface leafless crown fraction and/or in leaf demography. Canopy RTMs (PROSAIL and FLiES), driven by these three factors combined, simulated satellite-observed seasonal patterns well, explaining c. 70% of the variability in a key reflectance-based vegetation index (MAIAC EVI, which removes artifacts that would otherwise arise from clouds/aerosols and sun-sensor geometry). Leaf area index, leafless crown fraction and leaf demography independently accounted for 1, 33 and 66% of FLiES-simulated EVI seasonality, respectively. These factors also strongly influenced modeled near-infrared (NIR) reflectance, explaining why both modeled and observed EVI, which is especially sensitive to NIR, captures canopy seasonal dynamics well. Our improved analysis of canopy-scale biophysics rules out satellite artifacts as significant causes of satellite-observed seasonal patterns at this site, implying that aggregated phenology explains the larger scale remotely observed patterns. This work significantly reconciles current controversies about satellite-detected Amazon phenology, and improves our use of satellite observations to study climate-phenology relationships in the tropics.

The spatial sensitivity of the spectral diversity-biodiversity relationship: an experimental test in a prairie grassland

Remote sensing has been used to detect plant biodiversity in a range of ecosystems based on the varying spectral properties of different species or functional groups. However, the most appropriate spatial resolution necessary to detect diversity remains unclear. At coarse resolution, differences among spectral patterns may be too weak to detect. In contrast, at fine resolution, redundant information may be introduced. To explore the effect of spatial resolution, we studied the scale dependence of spectral diversity in a prairie ecosystem experiment at Cedar Creek Ecosystem Science Reserve, Minnesota, USA. Our study involved a scaling exercise comparing synthetic pixels resampled from high-resolution images within manipulated diversity treatments. Hyperspectral data were collected using several instruments on both ground and airborne platforms. We used the coefficient of variation (CV) of spectral reflectance in space as the indicator of spectral diversity and then compared CV at different scales ranging from 1 mm² to 1 m² to conventional biodiversity metrics, including species richness, Shannon's index, Simpson's index, phylogenetic species variation, and phylogenetic species evenness. In this study, higher species richness plots generally had higher CV. CV showed higher correlations with Shannon's index and Simpson's index than did species richness alone, indicating evenness contributed to the spectral diversity. Correlations with species richness and Simpson's index were generally higher than with phylogenetic species variation and evenness measured at comparable spatial scales, indicating weaker relationships between spectral diversity and phylogenetic diversity metrics than with species diversity metrics. High resolution imaging spectrometer data (1 mm² pixels) showed the highest sensitivity to diversity level. With decreasing spatial resolution, the difference in CV between diversity levels decreased and greatly reduced the optical detectability of biodiversity. The optimal pixel size for distinguishing α diversity in these prairie plots appeared to be around 1 mm to 10 cm, a spatial scale similar to the size of an individual herbaceous plant. These results indicate a strong scale-dependence of the spectral diversity-biodiversity relationships, with spectral diversity best able to detect a combination of species richness and evenness, and more weakly detecting phylogenetic diversity. These findings can be used to guide airborne studies of biodiversity and develop more effective large-scale biodiversity sampling methods.

Are optical indices good proxies of seasonal changes in carbon fluxes and stress-related physiological status in a beech forest?

This study investigates the functionality of a Mediterranean-mountain beech forest in Central Italy using simultaneous determinations of optical measurements, carbon (C) fluxes, leaf eco-physiological and biochemical traits during two growing seasons (2014-2015). Meteorological variables showed significant differences between the two growing seasons, highlighting a heat stress coupled with a reduced water availability in mid-summer 2015. As a result, a different C sink capacity of the forest was observed between the two years of study, due to the differences in stressful conditions and the related plant physiological status. Spectral indices related to vegetation (VIs, classified in structural, chlorophyll and carotenoid indices) were computed at top canopy level and used to track CO₂ fluxes and physiological changes. Optical indices related to structure (EVI 2, RDVI, DVI and MCARI 1) were found to better track Net Ecosystem Exchange (NEE) variations for 2014, while indices related to chlorophylls (SR red edge, CL red edge, MTCI and DR) provided better results for 2015. This suggests that when environmental conditions are not limiting for forest sink capacity, structural parameters are more strictly connected to C uptake, while under stress conditions indices related to functional features (e.g., chlorophyll content) become more relevant. Chlorophyll indices calculated with red edge bands (SR red edge, NDVI red edge, DR, CL red edge) resulted to be highly correlated with leaf nitrogen content (R²>0.70), while weaker, although significant, correlations were found with chlorophyll content. Carotenoid indices (PRI and PSRI) were strongly correlated with both chlorophylls and carotenoids content, suggesting that these indices are good proxies of the shifting pigment composition related to changes in soil moisture, heat stress and senescence. Our work suggests the importance of integrating different methods as a successful approach to understand how changing climatic conditions in the Mediterranean mountain region will impact on forest conditions and functionality.

Genetic, morphological, and spectral characterization of relictual Niobrara River hybrid aspens (Populus × smithii)

PREMISE OF THE STUDY

Aspen groves along the Niobrara River in Nebraska have long been a biogeographic curiosity due to morphological differences from nearby remnant Populus tremuloides populations. Pleistocene hybridization between P. tremuloides and P. grandidentata has been proposed, but the nearest P. grandidentata populations are currently several hundred kilometers east. We tested the hybrid-origin hypothesis using genetic data and characterized putative hybrids phenotypically.

METHODS

We compared nuclear microsatellite loci and chloroplast sequences of Niobrara River aspens to their putative parental species. Parental species and putative hybrids were also grown in a common garden for phenotypic comparison. On the common garden plants, we measured leaf morphological traits and leaf-level spectral reflectance profiles, from which chemical traits were derived.

KEY RESULTS

The genetic composition of the three unique Niobrara aspen genotypes is consistent with the hybridization hypothesis and with maternal chloroplast inheritance from P. grandidentata. Leaf margin dentition and abaxial pubescence differentiated taxa, with the hybrids showing intermediate values. Spectral profiles allowed statistical separation of taxa in short-wave infrared wavelengths, with hybrids showing intermediate values, indicating that traits associated with internal structure of leaves and water absorption may vary among taxa. However, reflectance values in the visible region did not differentiate taxa, indicating that traits related to pigments are not differentiated.

CONCLUSIONS

Both genetic and phenotypic results support the hypothesis of a hybrid origin for these genetically unique aspens. However, low genetic diversity and ongoing ecological and climatic threats to the hybrid taxon present a challenge for conservation of these relictual boreal communities.

Artificial Neural Network to Predict Vine Water Status Spatial Variability Using Multispectral Information Obtained from an Unmanned Aerial Vehicle (UAV)

Water stress, which affects yield and wine quality, is often evaluated using the midday stem water potential (Ψ_stem). However, this measurement is acquired on a per plant basis and does not account for the assessment of vine water status spatial variability. The use of multispectral cameras mounted on unmanned aerial vehicle (UAV) is capable to capture the variability of vine water stress in a whole field scenario. It has been reported that conventional multispectral indices (CMI) that use information between 500–800 nm, do not accurately predict plant water status since they are not sensitive to water content. The objective of this study was to develop artificial neural network (ANN) models derived from multispectral images to predict the Ψ_stem spatial variability of a drip-irrigated Carménère vineyard in Talca, Maule Region, Chile. The coefficient of determination (R²) obtained between ANN outputs and ground-truth measurements of Ψ_stem were between 0.56–0.87, with the best performance observed for the model that included the bands 550, 570, 670, 700 and 800 nm. Validation analysis indicated that the ANN model could estimate Ψ_stem with a mean absolute error (MAE) of 0.1 MPa, root mean square error (RMSE) of 0.12 MPa, and relative error (RE) of −9.1%. For the validation of the CMI, the MAE, RMSE and RE values were between 0.26–0.27 MPa, 0.32–0.34 MPa and −24.2–25.6%, respectively.

Hyperspectral reflectance sensing to assess the growth and photosynthetic properties of wheat cultivars exposed to different irrigation rates in an irrigated arid region

Simultaneous indirect assessment of multiple and diverse plant parameters in an exact and expeditious manner is becoming imperative in irrigated arid regions, with a view toward creating drought-tolerant genotypes or for the management of precision irrigation. This study aimed to evaluate whether spectral reflectance indices (SRIs) in three parts of the electromagnetic spectrum ((visible-infrared (VIS), near-infrared (NIR)), and shortwave-infrared (SWIR)) could be used to track changes in morphophysiological parameters of wheat cultivars exposed to 1.00, 0.75, and 0.50 of the estimated evapotranspiration (ETc). Significant differences were found in the parameters of growth and photosynthetic efficiency, and canopy spectral reflectance among the three cultivars subjected to different irrigation rates. All parameters were highly and significantly correlated with each other particularly under the 0.50 ETc treatment. The VIS/VIS- and NIR/VIS-based indices were sufficient and suitable for assessing the growth and photosynthetic properties of wheat cultivars similar to those indices based on NIR/NIR, SWIR/NIR, or SWIR/SWIR. Almost all tested SRIs proved to assess growth and photosynthetic parameters, including transpiration rate, more efficiently when regressions were analyzed for each water irrigation rate individually. This study, the type of which has rarely been conducted in irrigated arid regions, indicates that spectral reflectance data can be used as a rapid and non-destructive alternative method for assessment of the growth and photosynthetic efficiency of wheat under a range of water irrigation rates.

Hyperspectral Monitoring of Green Roof Vegetation Health State in Sub-Mediterranean Climate: Preliminary Results

In urban and industrial environments, the constant increase of impermeable surfaces has produced drastic changes in the natural hydrological cycle. Decreasing green areas not only produce negative effects from a hydrological-hydraulic perspective, but also from an energy point of view, modifying the urban microclimate and generating, as shown in the literature, heat islands in our cities. In this context, green infrastructures may represent an environmental compensation action that can be used to re-equilibrate the hydrological and energy balance and reduce the impact of pollutant load on receiving water bodies. To ensure that a green infrastructure will work properly, vegetated areas have to be continuously monitored to verify their health state. This paper presents a ground spectroscopy monitoring survey of a green roof installed at the University of Calabria fulfilled via the acquisition and analysis of hyperspectral data. This study is part of a larger research project financed by European Structural funds aimed at understanding the influence of green roofs on rainwater management and energy consumption for air conditioning in the Mediterranean area. Reflectance values were acquired with a field-portable spectroradiometer that operates in the range of wavelengths 350–2500 nm. The survey was carried out during the time period November 2014–June 2015 and data were acquired weekly. Climatic, thermo-physical, hydrological and hydraulic quantities were acquired as well and related to spectral data. Broadband and narrowband spectral indices, related to chlorophyll content and to chlorophyll–carotenoid ratio, were computed. The two narrowband indices NDVI₇₀₅ and SIPI turned out to be the most representative indices to detect the plant health status.

Canopy near-infrared reflectance and terrestrial photosynthesis

Global estimates of terrestrial gross primary production (GPP) remain highly uncertain, despite decades of satellite measurements and intensive in situ monitoring. We report a new approach for quantifying the near-infrared reflectance of terrestrial vegetation (NIR_V). NIR_V provides a foundation for a new approach to estimate GPP that consistently untangles the confounding effects of background brightness, leaf area, and the distribution of photosynthetic capacity with depth in canopies using existing moderate spatial and spectral resolution satellite sensors. NIR_V is strongly correlated with solar-induced chlorophyll fluorescence, a direct index of photons intercepted by chlorophyll, and with site-level and globally gridded estimates of GPP. NIR_V makes it possible to use existing and future reflectance data as a starting point for accurately estimating GPP.

Recognition of Orobanche cumana Below-Ground Parasitism Through Physiological and Hyper Spectral Measurements in Sunflower (Helianthus annuus L.)

Broomrape (Orobanche and Phelipanche spp.) parasitism is a severe problem in many crops worldwide, including in the Mediterranean basin. Most of the damage occurs during the sub-soil developmental stage of the parasite, by the time the parasite emerges from the ground, damage to the crop has already been done. One feasible method for sensing early, below-ground parasitism is through physiological measurements, which provide preliminary indications of slight changes in plant vitality and productivity. However, a complete physiological field survey is slow, costly and requires skilled manpower. In recent decades, visible to-shortwave infrared (VIS-SWIR) hyperspectral tools have exhibited great potential for faster, cheaper, simpler and non-destructive tracking of physiological changes. The advantage of VIS-SWIR is even greater when narrow-band signatures are analyzed with an advanced statistical technique, like a partial least squares regression (PLS-R). The technique can pinpoint the most physiologically sensitive wavebands across an entire spectrum, even in the presence of high levels of noise and collinearity. The current study evaluated a method for early detection of Orobanche cumana parasitism in sunflower that combines plant physiology, hyperspectral readings and PLS-R. Seeds of susceptible and resistant O. cumana sunflower varieties were planted in infested (15 mg kg^-1 seeds) and non-infested soil. The plants were examined weekly to detect any physiological or structural changes; the examinations were accompanied by hyperspectral readings. During the early stage of the parasitism, significant differences between infected and non-infected sunflower plants were found in the reflectance of near and shortwave infrared areas. Physiological measurements revealed no differences between treatments until O. cumana inflorescences emerged. However, levels of several macro- and microelements tended to decrease during the early stage of O. cumana parasitism. Analysis of leaf cross-sections revealed differences in range and in mesophyll structure as a result of different levels of nutrients in sunflower plants, manifesting the presence of O. cumana infections. The findings of an advanced PLS-R analysis emphasized the correlation between specific reflectance changes in the SWIR range and levels of various nutrients in sunflower plants. This work demonstrates potential for the early detection of O. cumana parasitism on sunflower roots using hyperspectral tools.

Global leaf trait estimates biased due to plasticity in the shade

The study of leaf functional trait relationships, the so-called leaf economics spectrum^1,2, is based on the assumption of high-light conditions (as experienced by sunlit leaves). Owing to the exponential decrease of light availability through canopies, however, the vast majority of the world's vegetation exists in at least partial shade. Plant functional traits vary in direct dependence of light availability³, with different traits varying to different degrees, sometimes in conflict with expectations from the economic spectrum³. This means that the derived trait relationships of the global leaf economic spectrum are probably dependent on the extent to which observed data in existing large-scale plant databases represent high-light conditions. Here, using an extensive worldwide database of within-canopy gradients of key physiological, structural and chemical traits³, along with three different global trait databases^4,5, we show that: (1) accounting for light-driven trait plasticity can reveal novel trait relationships, particularly for highly plastic traits (for example, the relationship between net assimilation rate per area (A_a) and leaf mass per area (LMA)); and (2) a large proportion of leaf traits in current global plant databases reported as measured in full sun were probably measured in the shade. The results show that even though the majority of leaves exist in the shade, along with a large proportion of observations, our current understanding is too focused on conditions in the sun.

Do the energy fluxes and surface conductance of boreal coniferous forests in Europe scale with leaf area?

Earth observing systems are now routinely used to infer leaf area index (LAI) given its significance in spatial aggregation of land surface fluxes. Whether LAI is an appropriate scaling parameter for daytime growing season energy budget, surface conductance (G_s ), water- and light-use efficiency and surface-atmosphere coupling of European boreal coniferous forests was explored using eddy-covariance (EC) energy and CO₂ fluxes. The observed scaling relations were then explained using a biophysical multilayer soil-vegetation-atmosphere transfer model as well as by a bulk G_s representation. The LAI variations significantly alter radiation regime, within-canopy microclimate, sink/source distributions of CO₂ , H₂ O and heat, and forest floor fluxes. The contribution of forest floor to ecosystem-scale energy exchange is shown to decrease asymptotically with increased LAI, as expected. Compared with other energy budget components, dry-canopy evapotranspiration (ET) was reasonably 'conservative' over the studied LAI range 0.5-7.0 m² m^-2 . Both ET and G_s experienced a minimum in the LAI range 1-2 m² m^-2 caused by opposing nonproportional response of stomatally controlled transpiration and 'free' forest floor evaporation to changes in canopy density. The young forests had strongest coupling with the atmosphere while stomatal control of energy partitioning was strongest in relatively sparse (LAI ~2 m² m^-2 ) pine stands growing on mineral soils. The data analysis and model results suggest that LAI may be an effective scaling parameter for net radiation and its partitioning but only in sparse stands (LAI <3 m² m^-2 ). This finding emphasizes the significance of stand-replacing disturbances on the controls of surface energy exchange. In denser forests, any LAI dependency varies with physiological traits such as light-saturated water-use efficiency. The results suggest that incorporating species traits and site conditions are necessary when LAI is used in upscaling energy exchanges of boreal coniferous forests.

Effects of climate and fire on short-term vegetation recovery in the boreal larch forests of Northeastern China

Understanding the influence of climate variability and fire characteristics in shaping postfire vegetation recovery will help to predict future ecosystem trajectories in boreal forests. In this study, I asked: (1) which remotely-sensed vegetation index (VI) is a good proxy for vegetation recovery? and (2) what are the relative influences of climate and fire in controlling postfire vegetation recovery in a Siberian larch forest, a globally important but poorly understood ecosystem type? Analysis showed that the shortwave infrared (SWIR) VI is a good indicator of postfire vegetation recovery in boreal larch forests. A boosted regression tree analysis showed that postfire recovery was collectively controlled by processes that controlled seed availability, as well as by site conditions and climate variability. Fire severity and its spatial variability played a dominant role in determining vegetation recovery, indicating seed availability as the primary mechanism affecting postfire forest resilience. Environmental and immediate postfire climatic conditions appear to be less important, but interact strongly with fire severity to influence postfire recovery. If future warming and fire regimes manifest as expected in this region, seed limitation and climate-induced regeneration failure will become more prevalent and severe, which may cause forests to shift to alternative stable states.

Biodiversity and agriculture in dynamic landscapes: Integrating ground and remotely-sensed baseline surveys

Sustainable biodiversity and land management require a cost-effective means of forecasting landscape response to environmental change. Conventional species-based, regional biodiversity assessments are rarely adequate for policy planning and decision making. We show how new ground and remotely-sensed survey methods can be coordinated to help elucidate and predict relationships between biodiversity, land use and soil properties along complex biophysical gradients that typify many similar landscapes worldwide. In the lower Zambezi valley, Mozambique we used environmental, gradient-directed transects (gradsects) to sample vascular plant species, plant functional types, vegetation structure, soil properties and land-use characteristics. Soil fertility indices were derived using novel multidimensional scaling of soil properties. To facilitate spatial analysis, we applied a probabilistic remote sensing approach, analyzing Landsat 7 satellite imagery to map photosynthetically active and inactive vegetation and bare soil along each gradsect. Despite the relatively low sample number, we found highly significant correlations between single and combined sets of specific plant, soil and remotely sensed variables that permitted testable spatial projections of biodiversity and soil fertility across the regional land-use mosaic. This integrative and rapid approach provides a low-cost, high-return and readily transferable methodology that permits the ready identification of testable biodiversity indicators for adaptive management of biodiversity and potential agricultural productivity.