Airborne Imaging Spectroscopy Research Articles

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

BackgroundBiogeographers 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.MethodsThese 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.ResultsCCB-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.DiscussionCCB-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.

Invasive tree species detection in the Eastern Arc Mountains biodiversity hotspot using one class classification

Eucalyptus spp. and Acacia mearnsii are common exotic tree species in eastern Africa that have shown (strong) invasive behavior in some regions. Acacia mearnsii is considered a highly invasive species that is replacing native species and Eucalyptus spp. are known to consume high amounts of groundwater with suspected effects on native flora. Mapping the occurrence of these species in the Taita Hills, Kenya (part of the Eastern Arc Mountains Biodiversity Hotspot) is important as there is lack of knowledge on their occurrence and ecological impact in the area. Mapping methods that require a lot of fieldwork are impractical in areas like the Taita Hills, where the terrain is rugged and the infrastructure is poor. Our aim was hence to map the occurrence of these tree species in a 100 km2 area using airborne imaging spectroscopy and laser scanning. We used a one class biased support vector machine (BSVM) classifier as it needs labeled training data only for the positive classes (A. mearnsii and Eucalyptus spp.), which potentially reduces the amount of required fieldwork. We also introduce a new approach for parameterizing and setting the threshold level simultaneously for the BSVM classifier. The second aim was to link the occurrence of these species to selected environmental variables. The results showed that the BSVM classifier is suitable for mapping Acacia mearnsii and Eucalyptus spp., holding the potential to improve the efficiency of field data collection. The introduced parametrization/threshold selection method performed better than other commonly used approaches. The crown level F1-score was 0.76 for Eucalyptus spp. and 0.78 for A. mearnsii. We show that Eucalyptus spp. and A. mearnsii trees cover 0.8% and 1.6% of the study area, respectively. Both species are particularly located on steeper slopes and higher altitudes. Both species have significant occurrences in areas close to the biggest remaining native forest patch (Ngangao) in the study area. Nonetheless, follow-up studies are needed to evaluate their impact on the native flora and fauna, as well as their impact on the water resources. The maps created in this study in combination with such follow-up studies could serve as base data to generate guidelines that authorities can use to take action in handling the problems these species are causing.

Anthropogenic aerosol emissions mapping and characterization by imaging spectroscopy – application to a metallurgical industry and a petrochemical complex

ABSTRACTThis paper is focused on the retrieval of industrial aerosol optical thickness (AOT) and microphysical properties by means of airborne imaging spectroscopy. Industrial emissions generally lead to optically thin plumes requiring an adapted detection method taking into account the weak proportion of particles sought in the atmosphere. To this end, a semi-analytical model combined with the Cluster-Tuned Matched Filter (CTMF) algorithm is presented to characterize those plumes, requiring the knowledge of the soil under the plume. The model allows the direct computation of the at-sensor radiance when a plume is included in the radiative transfer. When applied to industrial aerosol classes as defined in this paper, simulated spectral radiances can be compared to ‘real’ MODTRAN (Moderate Resolution Atmospheric Transmission) radiances using the Spectral Angle Mapper (SAM). On the range from 0.4 to 0.7 µm, for three grounds (water, vegetation, and bright one), SAM scores are lower than 0.043 in the worst case (a both absorbing and scattering particle over a bright ground), and usually lower than 0.025. The darker the ground reflectance is, the more accurate the results are (typically for reflectance lower than 0.3). Concerning AOT retrieval capabilities, with a pre-calculated model for a reference optical thickness of 0.25, we are able to retrieve plume AOT at 550 nm in the range 0.0 to 0.4 with an error usually ranging between 9% and 13%. The first test case is a CASI (Compact Airborne Spectrographic Imager) image acquired over the metallurgical industry of Fos-sur-Mer (France). First results of the use of the model coupled with CTMF algorithm reveal a scattering aerosol plume with particle sizes increasing with the distance from the stack (from detection score of 54% near the stack for particles with a diameter of 0.1 µm, to 69% away from it for 1.0 µm particles). A refinement is made then to estimate more precisely aerosol plume properties, using a multimodal distribution based on the previous results. It leads to find a mixture of sulfate and brown carbon particles with a plume AOT ranging between 0.2 and 0.5. The second test case is an AHS (Airborne Hyperspectral Scanner) image acquired over the petrochemical site of Antwerp (Belgium). The first CTMF application results in detecting a brown carbon aerosol of 0.1 µm mode (detection score is 51%). Refined results show the evolution of the AOT decreasing from 0.15 to 0.05 along the plume for a mixture of brown carbon fine mode and 0.3 µm radius of sulfate aerosol.

Previsual symptoms of Xylella fastidiosa infection revealed in spectral plant-trait alterations

Plant pathogens cause significant losses to agricultural yields and increasingly threaten food security1, ecosystem integrity and societies in general2-5. Xylella fastidiosa is one of the most dangerous plant bacteria worldwide, causing several diseases with profound impacts on agriculture and the environment6. Primarily occurring in the Americas, its recent discovery in Asia and Europe demonstrates that X. fastidiosa's geographic range has broadened considerably, positioning it as a reemerging global threat that has caused socioeconomic and cultural damage7,8. X. fastidiosa can infect more than 350 plant species worldwide9, and early detection is critical for its eradication8. In this article, we show that changes in plant functional traits retrieved from airborne imaging spectroscopy and thermography can reveal X. fastidiosa infection in olive trees before symptoms are visible. We obtained accuracies of disease detection, confirmed by quantitative polymerase chain reaction, exceeding 80% when high-resolution fluorescence quantified by three-dimensional simulations and thermal stress indicators were coupled with photosynthetic traits sensitive to rapid pigment dynamics and degradation. Moreover, we found that the visually asymptomatic trees originally scored as affected by spectral plant-trait alterations, developed X. fastidiosa symptoms at almost double the rate of the asymptomatic trees classified as not affected by remote sensing. We demonstrate that spectral plant-trait alterations caused by X. fastidiosa infection are detectable previsually at the landscape scale, a critical requirement to help eradicate some of the most devastating plant diseases worldwide.

Mapping canopy defoliation by herbivorous insects at the individual tree level using bi-temporal airborne imaging spectroscopy and LiDAR measurements

Defoliation by herbivorous insects is a widespread forest disturbance driver, affecting global forest health and ecosystem dynamics. Compared with time- and labor-intensive field surveys, remote sensing provides the only realistic approach to mapping canopy defoliation by herbivorous insects over large spatial and temporal scales. However, the spectral and structural signatures of defoliation by insects at the individual tree level have not been well studied. Additionally, the predictive power of spectral and structural metrics for mapping canopy defoliation has seldom been compared. These critical knowledge gaps prevent us from consistently detecting and mapping canopy defoliation by herbivorous insects across multiple scales. During the peak of a gypsy moth outbreak in Long Island, New York in summer 2016, we leveraged bi-temporal airborne imaging spectroscopy (IS, i.e., hyperspectral imaging) and LiDAR measurements at 1 m spatial resolution to explore the spectral and structural signatures of canopy defoliation in a mixed oak-pine forest. We determined that red edge and near-infrared spectral regions within the IS data were most sensitive to crown-scale defoliation severity. LiDAR measurements including B70 (i.e., 70th bincentile height), intensity skewness, and kurtosis were effectively able to detect structural changes caused by herbivorous insects. In addition to canopy leaf loss, increased exposure of understory and non-photosynthetic materials contributed to the detected spectral and structural signatures. Comparing the ability of individual sensors to map canopy defoliation, the LiDAR-only Ordinary Least-Square (OLS) model performed better than the IS-only model (Adj. R-squared = 0.77, RMSE = 15.37% vs. Adj. R-squared = 0.63, RMSE = 19.11%). The IS + LiDAR model improved on performance of the individual sensors (Adj. R-squared = 0.81, RMSE = 14.46%). Our study improves our understanding of spectral and structural signatures of defoliation by herbivorous insects and presents a novel approach for mapping insect defoliation at the individual tree level. Additionally, with the current and next generation of spaceborne sensors (e.g., WorldView-3, Landsat, Sentinel-2, HyspIRI, and GEDI), higher accuracy and frequent monitoring of insect defoliation may become more feasible across a range of spatial scales, which are critical for ecological research and management of forest resources including the economic consequences of forest insect infestations (e.g., reduced growth and increased mortality), as well as for informing and testing of carbon cycle models.

Measuring short-term post-fire forest recovery across a burn severity gradient in a mixed pine-oak forest using multi-sensor remote sensing techniques

Understanding post-fire forest recovery is pivotal to the study of forest dynamics and global carbon cycle. Field-based studies indicated a convex response of forest recovery rate to burn severity at the individual tree level, related with fire-induced tree mortality; however, these findings were constrained in spatial/temporal extents, while not detectable by traditional optical remote sensing studies, largely attributing to the contaminated effect from understory recovery. Here, we examined whether the combined use of multi-sensor remote sensing techniques (i.e., 1 m simultaneous airborne imaging spectroscopy and LiDAR and 2 m satellite multi-spectral imagery) to separate canopy recovery from understory recovery would enable to quantify post-fire forest recovery rate spanning a large gradient in burn severity over large-scales. Our study was conducted in a mixed pine-oak forest in Long Island, NY, three years after a top-killing fire. Our studies remotely detected an initial increase and then decline of forest recovery rate to burn severity across the burned area, with a maximum canopy area-based recovery rate of 10% per year at moderate forest burn severity class. More intriguingly, such remotely detected convex relationships also held at species level, with pine trees being more resilient to high burn severity and having a higher maximum recovery rate (12% per year) than oak trees (4% per year). These results are one of the first quantitative evidences showing the effects of fire adaptive strategies on post-fire forest recovery, derived from relatively large spatial-temporal scales. Our study thus provides the methodological advance to link multi-sensor remote sensing techniques to monitor forest dynamics in a spatially explicit manner over large-scales, with important implications for fire-related forest management and constraining/benchmarking fire effect schemes in ecological process models.

Plastic adaptations of foraging strategies to variation in forage quality in Alpine chamois (Rupicapra rupicapra)

Foraging efficiency strongly affects individual fitness and is influenced by diverse factors such as food quality and quantity, as well as intra- and inter-specific interactions. We investigated whether Alpine chamois (Rupicapra rupicapra (Linnaeus, 1758)) in a protected area in the Swiss Alps adapted their foraging behaviour to forage availability and quality by modeling the bite and step rates of individuals on vegetation nitrogen content, relative plant cover, sex, daytime, air temperature, and slope. Vegetation characteristics were derived using remote sensing data from airborne imaging spectroscopy data sets and feeding locations determined using a theodolite. Chamois increased their bite rates with decreasing forage nitrogen content, decreasing slope, and increasing temperature. Step rates were higher at high temperatures and decreased with increasing relative plant cover. Males showed higher bite rates and lower step rates than females. Daytime had no influence on either bite or step rates. An increase in bite rate may represent a plastic adaptation of foraging behaviour to compensate for lower nutritional quality of the available vegetation. Our results show variability in foraging behaviour according to both vegetation characteristics and physical environment and emphasize the use of remote sensing data to investigate relationships between habitat and subtle behavioural adaptations in ungulates.

Transferability of species distribution models for the detection of an invasive alien bryophyte using imaging spectroscopy data

Remote sensing is a promising tool for detecting invasive alien plant species. Mapping and monitoring those species requires accurate detection. So far, most studies relied on models that are locally calibrated and validated against available field data. Consequently, detecting invasive alien species at new study areas requires the acquisition of additional field data which can be expensive and time-consuming. Model transfer might thus provide a viable alternative. Here, we mapped the distribution of the invasive alien bryophyte Campylopus introflexus to i) assess the feasibility of spatially transferring locally calibrated models for species detection between four different heathland areas in Germany and Belgium and ii) test the potential of combining calibration data from different sites in one species distribution model (SDM). In a first step, four different SDMs were locally calibrated and validated by combining field data and airborne imaging spectroscopy data with a spatial resolution ranging from 1.8 m to 4 m and a spectral resolution of about 10 nm (244 bands). A one-class classifier, Maxent, which is based on the comparison of probability densities, was used to generate all SDMs. In a second step, each model was transferred to the three other study areas and the performance of the models for predicting C. introflexus occurrences was assessed. Finally, models combining calibration data from three study areas were built and tested on the remaining fourth site. In this step, different combinations of Maxent modelling parameters were tested. For the local models, the area under the curve for a test dataset (test AUC) was between 0.57–0.78, while the test AUC for the single transfer models ranged between 0.45–0.89. For the combined models the test AUC was between 0.54–0.9. The success of transferring models calibrated in one site to another site highly depended on the respective study site; the combined models provided higher test AUC values than the locally calibrated models for three out of four study sites. Furthermore, we also demonstrated the importance of optimizing the Maxent modelling parameters. Overall, our results indicate the potential of a combined model to map C. introflexus without the need for new calibration data.

An Approach for Foliar Trait Retrieval from Airborne Imaging Spectroscopy of Tropical Forests

Spatial information on forest functional composition is needed to inform management and conservation efforts, yet this information is lacking, particularly in tropical regions. Canopy foliar traits underpin the functional biodiversity of forests, and have been shown to be remotely measurable using airborne 350–2510 nm imaging spectrometers. We used newly acquired imaging spectroscopy data constrained with concurrent light detection and ranging (LiDAR) measurements from the Carnegie Airborne Observatory (CAO), and field measurements, to test the performance of the Spectranomics approach for foliar trait retrieval. The method was previously developed in Neotropical forests, and was tested here in the humid tropical forests of Malaysian Borneo. Multiple foliar chemical traits, as well as leaf mass per area (LMA), were estimated with demonstrable precision and accuracy. The results were similar to those observed for Neotropical forests, suggesting a more general use of the Spectranomics approach for mapping canopy traits in tropical forests. Future mapping studies using this approach can advance scientific investigations and applications based on imaging spectroscopy.

Cast Shadow Detection to Quantify the Aerosol Optical Thickness for Atmospheric Correction of High Spatial Resolution Optical Imagery

The atmospheric correction of optical remote sensing data requires the determination of aerosol and gas optical properties. A method is presented which allows the detection of the aerosol scattering effects from optical remote sensing data at spatial sampling intervals below 5 m in cloud-free situations from cast shadow pixels. The derived aerosol optical thickness distribution is used for improved atmospheric compensation. In a first step, a novel spectral cast shadow detection algorithm determines the shadow areas using spectral indices. Evaluation of the cast shadow masks shows an overall classification accuracy on an 80% level. Using the such derived shadow map, the ATCOR atmospheric compensation method is iteratively applied on the shadow areas in order to find the optimum aerosol amount. The aerosol optical thickness is found by analyzing the physical atmospheric correction of fully shaded pixels in comparison to directly illuminated areas. The shadow based aerosol optical thickness estimation method (SHAOT) is tested on airborne imaging spectroscopy data as well as on photogrammetric data. The accuracy of the reflectance values from atmospheric correction using the such derived aerosol optical thickness could be improved from 3–4% to a level of better than 2% in reflectance for the investigated test cases.

Evaluation of digital soil mapping approaches with large sets of environmental covariates

Abstract. The spatial assessment of soil functions requires maps of basic soil properties. Unfortunately, these are either missing for many regions or are not available at the desired spatial resolution or down to the required soil depth. The field-based generation of large soil datasets and conventional soil maps remains costly. Meanwhile, legacy soil data and comprehensive sets of spatial environmental data are available for many regions.Digital soil mapping (DSM) approaches relating soil data (responses) to environmental data (covariates) face the challenge of building statistical models from large sets of covariates originating, for example, from airborne imaging spectroscopy or multi-scale terrain analysis. We evaluated six approaches for DSM in three study regions in Switzerland (Berne, Greifensee, ZH forest) by mapping the effective soil depth available to plants (SD), pH, soil organic matter (SOM), effective cation exchange capacity (ECEC), clay, silt, gravel content and fine fraction bulk density for four soil depths (totalling 48 responses). Models were built from 300–500 environmental covariates by selecting linear models through (1) grouped lasso and (2) an ad hoc stepwise procedure for robust external-drift kriging (georob). For (3) geoadditive models we selected penalized smoothing spline terms by component-wise gradient boosting (geoGAM). We further used two tree-based methods: (4) boosted regression trees (BRTs) and (5) random forest (RF). Lastly, we computed (6) weighted model averages (MAs) from the predictions obtained from methods 1–5.Lasso, georob and geoGAM successfully selected strongly reduced sets of covariates (subsets of 3–6 % of all covariates). Differences in predictive performance, tested on independent validation data, were mostly small and did not reveal a single best method for 48 responses. Nevertheless, RF was often the best among methods 1–5 (28 of 48 responses), but was outcompeted by MA for 14 of these 28 responses. RF tended to over-fit the data. The performance of BRT was slightly worse than RF. GeoGAM performed poorly on some responses and was the best only for 7 of 48 responses. The prediction accuracy of lasso was intermediate. All models generally had small bias. Only the computationally very efficient lasso had slightly larger bias because it tended to under-fit the data. Summarizing, although differences were small, the frequencies of the best and worst performance clearly favoured RF if a single method is applied and MA if multiple prediction models can be developed.

Quantification of Soil Properties with Hyperspectral Data: Selecting Spectral Variables with Different Methods to Improve Accuracies and Analyze Prediction Mechanisms

We explored the potentials of both non-imaging laboratory and airborne imaging spectroscopy to assess arable soil quality indicators. We focused on microbial biomass-C (MBC) and hot water-extractable C (HWEC), complemented by organic carbon (OC) and nitrogen (N) as well-studied spectrally active parameters. The aggregation of different spectral variable selection strategies was used to analyze benefits for reachable estimation accuracies and to explore spectral predictive mechanisms for MBC and HWEC. With selected variables, quantification accuracies improved markedly for MBC (laboratory: RPD = 2.32 instead of 1.33 with full spectra; airborne: 2.35 instead of 1.80) and OC (laboratory: RPD = 3.08 instead of 2.36; airborne: 2.20 instead of 1.94). Patterns of selected variables indicated similarities between HWEC and OC, but significant differences between all other soil variables. This agreed to our results of indirect approaches in which both (i) wet-chemical data of OC and N and (ii) spectra fitted to measured OC and N values were used to estimate MBC and HWEC. Compared to these approaches, we found marked benefits of laboratory and airborne data for a direct spectral quantification of MBC (but not for HWEC). This suggests specificity of spectra for MBC, usable for the determination of this important soil parameter.

Sensitivity of Common Vegetation Indices to the Canopy Structure of Field Crops

Leaf inclination angle distribution is an important canopy structure characteristic which directly impacts the fraction of the intercepted solar radiation. Together with the leaf area index (LAI) it determines the structure and fractional cover of a homogeneous crop canopy. Unfortunately, this key canopy parameter has usually been ignored when applying common multispectral vegetation indices to the mapping of LAI, although its impact is known from model simulations. Therefore, we measured leaf angles and determined their distribution (quantified using the leaf mean tilt angle, MTA) for six crop species with different structures growing on 162 plots with a broad range of LAI (1.1–5.0) and leaf chlorophyll content (26–94 μg cm−2). Next, we calculated six vegetation indices widely used for LAI monitoring—the normalized difference vegetation index (NDVI), the enhanced vegetation index (EVI), the two band enhanced vegetation index (EVI2), the modified triangular vegetation index (MTVI2), the optimized soil adjusted vegetation index (OSAVI) and the wide dynamic range vegetation index (WDRVI)—from airborne imaging spectroscopy data. We calculated the Spearman’s correlation coefficient R s , a non-parametric statistic chosen because of the non-normal distribution of canopy parameters. All studied indices depended on the LAI ( 0.50 ≤ R s ≤ 0.71 ) , but the dependence on the MTA was of similar magnitude ( − 0.83 ≤ R s ≤ − 0.53 ) with EVI, EVI2, OSAVI and MTVI2 depending more strongly on MTA than on LAI. All studied indices were good proxies ( 0.78 ≤ R s ≤ 0.88 ) for vegetation fractional cover (Fcover) which, for homogeneous crop canopies, is a nonlinear function of LAI and MTA. EVI2 and MTVI2 were the most strongly correlated with Fcover, although the difference to the other studied indices was small. This first study involving a large range of crop structures confirms the results from canopy reflectance simulations and emphasizes the necessity of leaf angle information for the successful mapping of LAI with Earth observation data.

Effect of Spatial Filtering on Characterizing Soil Properties From Imaging Spectrometer Data

Airborne imaging spectroscopy covering wavelength range of 0.35–2.5 $\mu$ m can be used to quantify soil textural properties and chemical constituents. In this paper, we evaluate the effects of spatial resolution on the quantification of soil constituents using a lasso algorithm-based ensemble bootstrapping framework. Airborne visible infrared imaging spectrometer data collected at 7.6 m resolution over Bird's Point New Madrid (BPNM) floodway in Missouri, USA, is upscaled using a spatial filter to simulate a satellite-based sensor and generate multiple coarser resolution datasets, including the originally proposed 60.8 m hyperspectral infrared imager like data. The simulated data at multiple spatial resolutions are used in an ensemble lasso algorithm-based modeling framework for developing quantitative prediction models and spatial mapping of the soil constituents. We outline an evaluation framework with a set of metrics that considers the point-scale model performance as well as the consistency of cross-scale spatial predictions. The model results demonstrate that the ensemble quantification method is scalable, and further the model structure indicates the persistence of important spectral features across spatial resolutions. The probability density functions of the constituents over the BPNM landscape show that it is similar for multiple spatial resolutions. Finally, a comparison of the model predictions with statistical central values together with the within pixel variance across fine to coarse resolutions indicate that the model accurately captures the median values of the fine subgrid that the coarse-resolution data is composed of. This study establishes the feasibility for quantifying soil constituents from space-borne hyperspectral sensors.

Classification of Tree Species in a Diverse African Agroforestry Landscape Using Imaging Spectroscopy and Laser Scanning

Airborne imaging spectroscopy (IS) and laser scanning (ALS) have been explored widely for tree species classification during the past decades. However, African agroforestry areas, where a few exotic tree species are dominant and many native species occur less frequently, have not yet been studied. Obtaining maps of tree species would provide useful information for the characterization of agroforestry systems and detecting invasive species. Our objective was to study tree species classification in a diverse tropical landscape using IS and ALS data at the tree crown level, with primary interest in the exotic tree species. We performed multiple analyses based on different IS and ALS feature sets, identified important features using feature selection, and evaluated the impact of combining the two data sources. Given that a high number of tree species with limited sample size (499 samples for 31 species) was expected to limit the classification accuracy, we tested different approaches to group the species based on the frequency of their occurrence and Jeffries–Matusita (JM) distance. Surface reflectance at wavelengths between 400–450 nm and 750–800 nm, and height to crown width ratio, were identified as important features. Nonetheless, a selection of minimum noise fraction (MNF) transformed reflectance bands showed superior performance. Support vector machine classifier performed slightly better than the random forest classifier, but the improvement was not statistically significant for the best performing feature set. The highest F1-scores were achieved when each of the species was classified separately against a mixed group of all other species, which makes this approach suitable for invasive species detection. Our results are valuable for organizations working on biodiversity conservation and improving agroforestry practices, as we showed how the non-native Eucalyptus spp., Acacia mearnsii and Grevillea robusta (mean F1-scores 76%, 79% and 89%, respectively) trees can be mapped with good accuracy. We also found a group of six fruit bearing trees using JM distance, which was classified with mean F1-score of 65%. This was a useful finding, as these species could not be classified with acceptable accuracy individually, while they all share common economic and ecological importance.

On the challenges of using field spectroscopy to measure the impact of soil type on leaf traits

Abstract. Understanding the causes of variation in functional plant traits is a central issue in ecology, particularly in the context of global change. Spectroscopy is increasingly used for rapid and non-destructive estimation of foliar traits, but few studies have evaluated its accuracy when assessing phenotypic variation in multiple traits. Working with 24 chemical and physical leaf traits of six European tree species growing on strongly contrasting soil types (i.e. deep alluvium versus nearby shallow chalk), we asked (i) whether variability in leaf traits is greater between tree species or soil type, and (ii) whether field spectroscopy is effective at predicting intraspecific variation in leaf traits as well as interspecific differences. Analysis of variance showed that interspecific differences in traits were generally much stronger than intraspecific differences related to soil type, accounting for 25 % versus 5 % of total trait variation, respectively. Structural traits, phenolic defences and pigments were barely affected by soil type. In contrast, foliar concentrations of rock-derived nutrients did vary: P and K concentrations were lower on chalk than alluvial soils, while Ca, Mg, B, Mn and Zn concentrations were all higher, consistent with the findings of previous ecological studies. Foliar traits were predicted from 400 to 2500 nm reflectance spectra collected by field spectroscopy using partial least square regression, a method that is commonly employed in chemometrics. Pigments were best modelled using reflectance data from the visible region (400–700 nm), while all other traits were best modelled using reflectance data from the shortwave infrared region (1100–2500 nm). Spectroscopy delivered accurate predictions of species-level variation in traits. However, it was ineffective at detecting intraspecific variation in rock-derived nutrients (with the notable exception of P). The explanation for this failure is that rock-derived elements do not have absorption features in the 400–2500 nm region, and their estimation is indirect, relying on elemental concentrations covarying with structural traits that do have absorption features in that spectral region (constellation effects). Since the structural traits did not vary with soil type, it was impossible for our regression models to predict intraspecific variation in rock-derived nutrients via constellation effects. This study demonstrates the value of spectroscopy for rapid, non-destructive estimation of foliar traits across species, but highlights problems with predicting intraspecific variation indirectly. We discuss the implications of these findings for mapping functional traits by airborne imaging spectroscopy.

Tree Density and Forest Productivity in a Heterogeneous Alpine Environment: Insights from Airborne Laser Scanning and Imaging Spectroscopy

We outline an approach combining airborne laser scanning (ALS) and imaging spectroscopy (IS) to quantify and assess patterns of tree density (TD) and forest productivity (FP) in a protected heterogeneous alpine forest in the Swiss National Park (SNP). We use ALS data and a local maxima (LM) approach to predict TD, as well as IS data (Airborne Prism Experiment—APEX) and an empirical model to estimate FP. We investigate the dependency of TD and FP on site related factors, in particular on surface exposition and elevation. Based on reference data (i.e., 1598 trees measured in 35 field plots), we observed an underestimation of ALS-based TD estimates of 40%. Our results suggest a limited sensitivity of the ALS approach to small trees as well as a dependency of TD estimates on canopy heterogeneity, structure, and species composition. We found a weak to moderate relationship between surface elevation and TD (R2 = 0.18–0.69) and a less pronounced trend with FP (R2 = 0.0–0.56), suggesting that both variables depend on gradients of resource availability. Further to the limitations faced in the sensitivity of the applied approaches, we conclude that the combined application of ALS and IS data was convenient for estimating tree density and mapping FP in north-facing forested areas, however, the accuracy was lower in south-facing forested areas covered with multi-stemmed trees.

Habitat Mapping and Quality Assessment of NATURA 2000 Heathland Using Airborne Imaging Spectroscopy

Appropriate management of (semi-)natural areas requires detailed knowledge of the ecosystems present and their status. Remote sensing can provide a systematic, synoptic view at regular time intervals, and is therefore often suggested as a powerful tool to assist with the mapping and monitoring of protected habitats and vegetation. In this study, we present a multi-step mapping framework that enables detailed NATURA 2000 (N2000) heathland habitat patch mapping and the assessment of their conservation status at patch level. The method comprises three consecutive steps: (1) a hierarchical land/vegetation type (LVT) classification using airborne AHS imaging spectroscopy and field reference data; (2) a spatial re-classification to convert the LVT map to a patch map based on life forms; and (3) identification of the N2000 habitat type and conservation status parameters for each of the patches. Based on a multivariate analysis of 1325 vegetation reference plots acquired in 2006–2007, 24 LVT classes were identified that were considered relevant for the assessment of heathland conservation status. These labelled data were then used as ground reference for the supervised classification of the AHS image data to an LVT classification map, using Linear Discriminant Analysis in combination with Sequential-Floating-Forward-Search feature selection. Overall classification accuracies for the LVT mapping varied from 83% to 92% (Kappa ≈ 0.82–0.91), depending on the level of detail in the hierarchical classification. After converting the LVT map to a N2000 habitat type patch map, an overall accuracy of 89% was obtained. By combining the N2000 habitat type patch map with the LVT map, two important conservation status parameters were directly deduced per patch: tree and shrub cover, and grass cover, showing a strong similarity to an independent dataset with estimates made in the field in 2009. The results of this study indicate the potential of imaging spectroscopy for detailed heathland habitat characterization of N2000 sites in a way that matches the current field-based workflows of the user.

Field and Airborne Spectroscopy Cross Validation—Some Considerations

Field spectroscopy is increasingly used in various fields of science: either as a research tool in its own right or in support of airborne- or space-based optical instruments for calibration or validation purposes. Yet, while the use of the instruments appears deceptively simple, the processes of light and surface interactions are complex to be measured in full and are further complicated by the multidimensionality of the measurement process. This study exemplifies the cross validation of in situ point spectroscopy and airborne imaging spectroscopy data across all processing stages within the spectroscopy information hierarchy using data from an experiment focused on vegetation. In support of this endeavor, this study compiles the fundamentals of spectroscopy, the challenges inherent to field and airborne spectroscopy, and the best practices proposed by the field spectroscopy community. This combination of theory and case study shall enable the reader to develop an understanding of 1) some of the commonly involved sources of errors and uncertainties, 2) the techniques to collect high-quality spectra under natural illumination conditions, and 3) the importance of appropriate metadata collection to increase the long-term usability and value of spectral data.

Landscape‐scale GPP and carbon density inform patterns and impacts of an invasive tree across wet forests of Hawaii

Plant invasion typically occurs within a landscape-scale framework of abiotic and biotic conditions, often resulting in emergent feedbacks among environment, ecosystem functions, and the dominance of invasive species. Understanding the mechanisms underlying successful invasions is an important component of conservation and management efforts, but this has been poorly investigated in a spatially explicit manner. Knowing where and why invasion patterns change throughout the landscape enables managers to use context-specific controls on the spread of invasive species. Using high-resolution airborne imaging spectroscopy, we studied plant performance in growth within and across landscapes to examine the dominance and spatial distribution of an invasive tree, Psidium cattleianum (strawberry guava), in heterogeneous environmental conditions of a submontane Hawaiian tropical forest. We assessed invader performance using the GPP ratio index, which is the relative difference in remotely sensed estimates of gross primary productivity between canopies of guava and canopies of the invaded plant community. In addition, we used airborne LiDAR data to evaluate the impacts of guava invasion on the forest aboveground carbon density in different environments. Structural equation modeling revealed that substrate type and elevation above sea level interact and amplify landscape-scale differences in productivity between the invasive species and the host plant community (GPP ratio); differences that ultimately control levels of dominance of guava. We found shifts in patterns of forest carbon storage based on both gradual increase of invader dominance and changes in environmental conditions. Overall, our results demonstrate that the remotely sensed index defined as the GPP ratio provided an innovative spatially explicit approach to track and predict the success of invasive plants based in their canopy productivity, particularly within a landscape-scale framework of varying environmental factors such as soils and elevation. This approach may help managers accurately predict where invaders of forests, scrublands, or grasslands are likely to exhibit high levels of dominance before the environment is fully invaded.