Commercial tree species discrimination using airborne AISA Eagle hyperspectral imagery and partial least squares discriminant analysis (PLS-DA) in KwaZulu–Natal,South Africa

Discriminating commercial tree species using hyperspectral remote sensing techniques is critical in monitoring the spatial distributions and compositions of commercial forests. However, issues related to data dimensionality and multicollinearity limit the successful application of the technology. The aim of this study was to examine the utility of the partial least squares discriminant analysis (PLS-DA) technique in accurately classifying six exotic commercial forest species (Eucalyptus grandis, Eucalyptus nitens, Eucalyptus smithii, Pinus patula, Pinus elliotii and Acacia mearnsii) using airborne AISA Eagle hyperspectral imagery (393–900 nm). Additionally, the variable importance in the projection (VIP) method was used to identify subsets of bands that could successfully discriminate the forest species. Results indicated that the PLS-DA model that used all the AISA Eagle bands (n = 230) produced an overall accuracy of 80.61% and a kappa value of 0.77, with user’s and producer’s accuracies ranging from 50% to 100%. In comparison, incorporating the optimal subset of VIP selected wavebands (n = 78) in the PLS-DA model resulted in an improved overall accuracy of 88.78% and a kappa value of 0.87, with user’s and producer’s accuracies ranging from 70% to 100%. Bands located predominantly within the visible region of the electromagnetic spectrum (393–723 nm) showed the most capability in terms of discriminating between the six commercial forest species. Overall, the research has demonstrated the potential of using PLS-DA for reducing the dimensionality of hyperspectral datasets as well as determining the optimal subset of bands to produce the highest classification accuracies.     

Spectral resampling based on user-defined inter-band correlation filter: C3 and C4 grass species classification

In this paper, a user-defined inter-band correlation filter function was used to resample hyperspectral data and thereby mitigate the problem of multicollinearity in classification analysis. The proposed resampling technique convolves the spectral dependence information between a chosen band-centre and its shorter and longer wavelength neighbours. Weighting threshold of inter-band correlation (WTC, Pearson's r) was calculated, whereby r = 1 at the band-centre. Various WTC (r = 0.99, r = 0.95 and r = 0.90) were assessed, and bands with coefficients beyond a chosen threshold were assigned r = 0. The resultant data were used in the random forest analysis to classify in situ C₃ and C₄ grass canopy reflectance. The respective WTC datasets yielded improved classification accuracies (kappa = 0.82, 0.79 and 0.76) with less correlated wavebands when compared to resampled Hyperion bands (kappa = 0.76). Overall, the results obtained from this study suggested that resampling of hyperspectral data should account for the spectral dependence information to improve overall classification accuracy as well as reducing the problem of multicollinearity.     

Evaluating the robustness of models developed from field spectral data in predicting African grass foliar nitrogen concentration using WorldView-2 image as an independent test dataset

In this paper, we evaluate the extent to which the resampled field spectra compare with the actual image spectra of the new generation multispectral WorldView-2 (WV-2) satellite. This was achieved by developing models from resampled field spectra data and testing them on an actual WV-2 image of the study area. We evaluated the performance of reflectance ratios (RI), normalized difference indices (NDI) and random forest (RF) regression model in predicting foliar nitrogen concentration in a grassland environment. The field measured spectra were used to calibrate the RF model using a randomly selected training (n = 70%) nitrogen data set. The model developed from the field spectra resampled to WV-2 wavebands was validated on an independent field spectral test dataset as well as on the actual WV-2 image of the same area (n = 30%, bootstrapped a 100 times). The results show that the model developed using RI could predict nitrogen with a mean R² of 0.74 and 0.65 on an independent field spectral test data set and on the actual WV-2 image, respectively. The root mean square error of prediction (RMSE %) was 0.17 and 0.22 for the field test data set and the WV-2 image, respectively. Results provide an insight on the magnitude of errors that are expected when up-scaling field spectral models to airborne or satellite image data. The prediction also indicates the unceasing relevance of field spectroscopy studies to better understand the spectral models critical for vegetation quality assessment.     

Modeling the spatial distribution of above-ground carbon in Mexican coniferous forests using remote sensing and a geostatistical approach

Forest conservation is considered an option for mitigating the effect of greenhouse gases on global climate, hence monitoring forest carbon pools at global and local levels is important. The present study explores the capability of remote-sensing variables (vegetation indices and textures derived from SPOT-5; backscattering coefficient and interferometric coherence of ALOS PALSAR images) for modeling the spatial distribution of above-ground biomass in the Environmental Conservation Zone of Mexico City. Correlation and spatial autocorrelation coefficients were used to select significant explanatory variables in fir and pine forests. The correlation for interferometric coherence in HV polarization was negative, with correlations coefficients r = −0.83 for the fir and r = −0.75 for the pine forests. Regression-kriging showed the least root mean square error among the spatial interpolation methods used, with 37.75 tC/ha for fir forests and 29.15 tC/ha for pine forests. The results showed that a hybrid geospatial method, based on interferometric coherence data and a regression-kriging interpolator, has good potential for estimating above-ground biomass carbon.     

Assessing plantation canopy condition from airborne imagery using spectral mixture analysis and fractional abundances

Pine plantations in Australia are subject to a range of abiotic and biotic damaging agents that affect tree health and productivity. In order to optimise management decisions, plantation managers require regular intelligence relating to the status and trends in the health and condition of trees within individual compartments. Remote sensing technology offers an alternative to traditional ground-based assessment of these plantations. Automated estimation of foliar crown health, especially in degraded crowns, can be difficult due to mixed pixels when there is low or fragmented vegetation cover. In this study we apply a linear spectral unmixing approach to high spatial resolution (50 cm) multispectral imagery to quantify the fractional abundances of the key image endmembers: sunlit canopy, shadow, and soil. A number of Pinus radiata tree crown attributes were modelled using multiple linear regression and endmember fraction images. We found high levels of significance (r² = 0.80) for the overall crown colour and colour of the crown leader (r² = 0.79) in tree crowns affected by the fungal pathogen Sphaeropsis sapinea, which produces both needle necrosis and chlorosis. Results for stands associated with defoliation and chlorosis through infestation by the aphid Essigella californica were lower with an r² = 0.33 for crown transparency and r² = 0.31 for proportion of crown affected. Similar analysis of data from a nitrogen deficient site produced an outcome somewhat in between the other two damaging agents. Overall the sunlit canopy image fraction has been the most important variable used in the modelling of forest condition for all damaging agents.     

Augmenting forest inventory attributes with geometric optical modelling in support of regional susceptibility assessments to bark beetle infestations

Assessment of the susceptibility of forests to mountain pine beetle (Dendroctonus ponderosae Hopkins) infestation is based upon an understanding of the characteristics that predispose the stands to attack. These assessments are typically derived from conventional forest inventory data; however, this information often represents only managed forest areas. It does not cover areas such as forest parks or conservation regions and is often not regularly updated resulting in an inability to assess forest susceptibility. To address these shortcomings, we demonstrate how a geometric optical model (GOM) can be applied to Landsat-5 Thematic Mapper (TM) imagery (30 m spatial resolution) to estimate stand-level susceptibility to mountain pine beetle attack. Spectral mixture analysis was used to determine the proportion of sunlit canopy and background, and shadow of each Landsat pixel enabling per pixel estimates of attributes required for model inversion. Stand structural attributes were then derived from inversion of the geometric optical model and used as basis for susceptibility mapping. Mean stand density estimated by the geometric optical model was 2753 (standard deviation ± 308) stems per hectare and mean horizontal crown radius was 2.09 (standard deviation ± 0.11) metres. When compared to equivalent forest inventory attributes, model predictions of stems per hectare and crown radius were shown to be reasonably estimated using a Kruskal–Wallis ANOVA (p < 0.001). These predictions were then used to create a large area map that provided an assessment of the forest area susceptible to mountain pine beetle damage.     

Retrieval of leaf water content spanning the visible to thermal infrared spectra

The objective of this study was to investigate the entire spectra (from visible to the thermal infrared; 0.390–14.0 μm) to retrieve leaf water content in a consistent manner. Narrow-band spectral indices (calculated from all possible two band combinations) and a partial least square regression (PLSR) were used to assess the strength of each spectral region. The coefficient of determination (R²) and root mean square error (RMSE) were used to report the prediction accuracy of spectral indices and PLSR models. In the visible-near infrared and shortwave infrared (VNIR–SWIR), the most accurate spectral index yielded R² of 0.89 and RMSE of 7.60%, whereas in the mid infrared (MIR) the highest R² was 0.93 and RMSE of 5.97%. Leaf water content was poorly predicted using two-band indices developed from the thermal infrared (R² = 0.33). The most accurate PLSR model resulted from MIR reflectance spectra (R² = 0.96, RMSE = 4.74% and RMSE cross validation RMSE_CV = 6.17%) followed by VNIR–SWIR reflectance spectra (R² = 0.91, RMSE = 6.90% and RMSE_CV = 7.32%). Using thermal infrared (TIR) spectra, the PLSR model yielded a moderate retrieval accuracy (R² = 0.67, RMSE = 13.27% and RMSE_CV = 16.39%). This study demonstrated that the mid infrared (MIR) and shortwave infrared (SWIR) domains were the most sensitive spectral region for the retrieval of leaf water content.     

Modelling the standing timber volume of Baden-Württemberg—A large-scale approach using a fusion of Landsat,airborne LiDAR and National Forest Inventory data

Remote sensing-based timber volume estimation is key for modelling the regional potential, accessibility and price of lignocellulosic raw material for an emerging bioeconomy. We used a unique wall-to-wall airborne LiDAR dataset and Landsat 7 satellite images in combination with terrestrial inventory data derived from the National Forest Inventory (NFI), and applied generalized additive models (GAM) to estimate spatially explicit timber distribution and volume in forested areas. Since the NFI data showed an underlying structure regarding size and ownership, we additionally constructed a socio-economic predictor to enhance the accuracy of the analysis. Furthermore, we balanced the training dataset with a bootstrap method to achieve unbiased regression weights for interpolating timber volume. Finally, we compared and discussed the model performance of the original approach (r² = 0.56, NRMSE = 9.65%), the approach with balanced training data (r² = 0.69, NRMSE = 12.43%) and the final approach with balanced training data and the additional socio-economic predictor (r² = 0.72, NRMSE = 12.17%). The results demonstrate the usefulness of remote sensing techniques for mapping timber volume for a future lignocellulose-based bioeconomy.     

Remote quantification of phycocyanin in potable water sources through an adaptive model

Cyanobacterial blooms in water supply sources in both central Indiana USA (CIN) and South Australia (SA) are a cause of great concerns for toxin production and water quality deterioration. Remote sensing provides an effective approach for quick assessment of cyanobacteria through quantification of phycocyanin (PC) concentration. In total, 363 samples spanning a large variation of optically active constituents (OACs) in CIN and SA waters were collected during 24 field surveys. Concurrently, remote sensing reflectance spectra (R_rs) were measured. A partial least squares–artificial neural network (PLS–ANN) model, artificial neural network (ANN) and three-band model (TBM) were developed or tuned by relating the R_rs with PC concentration. Our results indicate that the PLS–ANN model outperformed the ANN and TBM with both the original spectra and simulated ESA/Sentinel-3/Ocean and Land Color Instrument (OLCI) and EO-1/Hyperion spectra. The PLS–ANN model resulted in a high coefficient of determination (R²) for CIN dataset (R² = 0.92, R: 0.3–220.7 μg/L) and SA (R² = 0.98, R: 0.2–13.2 μg/L). In comparison, the TBM model yielded an R² = 0.77 and 0.94 for the CIN and SA datasets, respectively; while the ANN obtained an intermediate modeling accuracy (CIN: R² = 0.86; SA: R² = 0.95). Applying the simulated OLCI and Hyperion aggregated datasets, the PLS–ANN model still achieved good performance (OLCI: R² = 0.84; Hyperion: R² = 0.90); the TBM also presented acceptable performance for PC estimations (OLCI: R² = 0.65, Hyperion: R² = 0.70). Based on the results, the PLS–ANN is an effective modeling approach for the quantification of PC in productive water supplies based on its effectiveness in solving the non-linearity of PC with other OACs. Furthermore, our investigation indicates that the ratio of inorganic suspended matter (ISM) to PC concentration has close relationship to modeling relative errors (CIN: R² = 0.81; SA: R² = 0.92), indicating that ISM concentration exert significant impact on PC estimation accuracy.     

PTrees: A point-based approach to forest tree extraction from lidar data

This paper introduces PTrees, a multi-scale dynamic point cloud segmentation dedicated to forest tree extraction from lidar point clouds. The method process the point data using the raw elevation values (Z) and compute height (H = Z − ground elevation) during post-processing using an innovative procedure allowing to preserve the geometry of crown points. Multiple segmentations are done at different scales. Segmentation criteria are then applied to dynamically select the best set of apices from the tree segments extracted at the various scales. The selected set of apices is then used to generate a final segmentation. PTrees has been tested in 3 different forest types, allowing to detect 82% of the trees with under 10% of false detection rate. Future development will integrate crown profile estimation during the segmentation process in order to both maximize the detection of suppressed trees and minimize false detections.     

Improving the prediction of arsenic contents in agricultural soils by combining the reflectance spectroscopy of soils and rice plants

Visible and near-infrared reflectance spectroscopy provides a beneficial tool for investigating soil heavy metal contamination. This study aimed to investigate mechanisms of soil arsenic prediction using laboratory based soil and leaf spectra, compare the prediction of arsenic content using soil spectra with that using rice plant spectra, and determine whether the combination of both could improve the prediction of soil arsenic content. A total of 100 samples were collected and the reflectance spectra of soils and rice plants were measured using a FieldSpec3 portable spectroradiometer (350–2500 nm). After eliminating spectral outliers, the reflectance spectra were divided into calibration (n = 62) and validation (n = 32) data sets using the Kennard-Stone algorithm. Genetic algorithm (GA) was used to select useful spectral variables for soil arsenic prediction. Thereafter, the GA-selected spectral variables of the soil and leaf spectra were individually and jointly employed to calibrate the partial least squares regression (PLSR) models using the calibration data set. The regression models were validated and compared using independent validation data set. Furthermore, the correlation coefficients of soil arsenic against soil organic matter, leaf arsenic and leaf chlorophyll were calculated, and the important wavelengths for PLSR modeling were extracted. Results showed that arsenic prediction using the leaf spectra (coefficient of determination in validation, R_v² = 0.54; root mean square error in validation, RMSE_v = 12.99 mg kg⁻¹; and residual prediction deviation in validation, RPD_v = 1.35) was slightly better than using the soil spectra (R_v² = 0.42, RMSE_v = 13.35 mg kg⁻¹, and RPD_v = 1.31). However, results also showed that the combinational use of soil and leaf spectra resulted in higher arsenic prediction (R_v² = 0.63, RMSE_v = 11.94 mg kg⁻¹, RPD_v = 1.47) compared with either soil or leaf spectra alone. Soil spectral bands near 480, 600, 670, 810, 1980, 2050 and 2290 nm, leaf spectral bands near 700, 890 and 900 nm in PLSR models were important wavelengths for soil arsenic prediction. Moreover, soil arsenic showed significantly positive correlations with soil organic matter (r = 0.62, p < 0.01) and leaf arsenic (r = 0.77, p < 0.01), and a significantly negative correlation with leaf chlorophyll (r = −0.67, p < 0.01). The results showed that the prediction of arsenic contents using soil and leaf spectra may be based on their relationships with soil organic matter and leaf chlorophyll contents, respectively. Although RPD of 1.47 was below the recommended RPD of >2 for soil analysis, arsenic prediction in agricultural soils can be improved by combining the leaf and soil spectra.     

Plant phenolics and absorption features in vegetation reflectance spectra near 1.66 μm

Past laboratory and field studies have quantified phenolic substances in vegetative matter from reflectance measurements for understanding plant response to herbivores and insect predation. Past remote sensing studies on phenolics have evaluated crop quality and vegetation patterns caused by bedrock geology and associated variations in soil geochemistry. We examined spectra of pure phenolic compounds, common plant biochemical constituents, dry leaves, fresh leaves, and plant canopies for direct evidence of absorption features attributable to plant phenolics. Using spectral feature analysis with continuum removal, we observed that a narrow feature at 1.66 μm is persistent in spectra of manzanita, sumac, red maple, sugar maple, tea, and other species. This feature was consistent with absorption caused by aromatic CH bonds in the chemical structure of phenolic compounds and non-hydroxylated aromatics. Because of overlapping absorption by water, the feature was weaker in fresh leaf and canopy spectra compared to dry leaf measurements. Simple linear regressions of feature depth and feature area with polyphenol concentration in tea resulted in high correlations and low errors (% phenol by dry weight) at the dry leaf (r² = 0.95, RMSE = 1.0%, n = 56), fresh leaf (r² = 0.79, RMSE = 2.1%, n = 56), and canopy (r² = 0.78, RMSE = 1.0%, n = 13) levels of measurement. Spectra of leaves, needles, and canopies of big sagebrush and evergreens exhibited a weak absorption feature centered near 1.63 μm, short ward of the phenolic compounds, possibly consistent with terpenes. This study demonstrates that subtle variation in vegetation spectra in the shortwave infrared can directly indicate biochemical constituents and be used to quantify them. Phenolics are of lesser abundance compared to the major plant constituents but, nonetheless, have important plant functions and ecological significance. Additional research is needed to advance our understanding of the spectral influences of plant phenolics and terpenes relative to dominant leaf biochemistry (water, chlorophyll, protein/nitrogen, cellulose, and lignin).     

A new calibration of the effective scattering albedo and soil roughness parameters in the SMOS SM retrieval algorithm

This study focuses on the calibration of the effective vegetation scattering albedo (ω) and surface soil roughness parameters (H_R, and N_Rp, p = H,V) in the Soil Moisture (SM) retrieval from L-band passive microwave observations using the L-band Microwave Emission of the Biosphere (L-MEB) model. In the current Soil Moisture and Ocean Salinity (SMOS) Level 2 (L2), v620, and Level 3 (L3), v300, SM retrieval algorithms, low vegetated areas are parameterized by ω = 0 and H_R = 0.1, whereas values of ω = 0.06 − 0.08 and H_R = 0.3 are used for forests. Several parameterizations of the vegetation and soil roughness parameters (ω, H_R and N_Rp, p = H,V) were tested in this study, treating SMOS SM retrievals as homogeneous over each pixel instead of retrieving SM over a representative fraction of the pixel, as implemented in the operational SMOS L2 and L3 algorithms. Globally-constant values of ω = 0.10, H_R = 0.4 and N_Rp = −1 (p = H,V) were found to yield SM retrievals that compared best with in situ SM data measured at many sites worldwide from the International Soil Moisture Network (ISMN). The calibration was repeated for collections of in situ sites classified in different land cover categories based on the International Geosphere-Biosphere Programme (IGBP) scheme. Depending on the IGBP land cover class, values of ω and H_R varied, respectively, in the range 0.08–0.12 and 0.1–0.5. A validation exercise based on in situ measurements confirmed that using either a global or an IGBP-based calibration, there was an improvement in the accuracy of the SM retrievals compared to the SMOS L3 SM product considering all statistical metrics (R = 0.61, bias = −0.019 m³ m⁻³, ubRMSE = 0.062 m³ m⁻³ for the IGBP-based calibration; against R = 0.54, bias = −0.034 m³ m⁻³ and ubRMSE = 0.070 m³ m⁻³ for the SMOS L3 SM product). This result is a key step in the calibration of the roughness and vegetation parameters in the operational SMOS retrieval algorithm. The approach presented here is the core of a new forthcoming SMOS optimized SM product.     

Classifying individual tree species under leaf-off and leaf-on conditions using airborne lidar

In this paper, a methodology for individual tree-based species classification using high sampling density and small footprint lidar data is clarified, corrected and improved. For this purpose, a well-defined directed graph (digraph) is introduced and it plays a fundamental role in the approach. It is argued that there exists one and only one such unique digraph that describes all four pure events and resulting disjoint sets of laser points associated with a single tree in data from a two-return lidar system. However, the digraph is extendable so that it fits an n-return lidar system (n > 2) with higher logical resolution. Furthermore, a mathematical notation for different types of groupings of the laser points is defined, and a new terminology for various types of individual tree-based concepts defined by the digraph is proposed. A novel calibration technique for estimating individual tree heights is evaluated. The approach replaces the unreliable maximum single laser point height of each tree with a more reliable prediction based on shape characteristics of a marginal height distribution of the whole first-return point cloud of each tree. The result shows a reduction of the RMSE of the tree heights of about 20% (stddev = 1.1 m reduced to stddev = 0.92 m). The method improves the species classification accuracy markedly, but it could also be used for reducing the sampling density at the time of data acquisition. Using the calibrated tree heights, a scale-invariant rescaled space for the universal set of points for each tree is defined, in which all individual tree-based geometric measurements are conducted. With the corrected and improved classification methodology the total accuracy raises from 60% to 64% for classifying three leaf-off individual tree deciduous species (N = 200 each) in West Virginia, USA: oaks (Quercus spp.), red maple (Acer rubrum), and yellow poplar (Liriodendron tuliperifera).     

Aboveground total and green biomass of dryland shrub derived from terrestrial laser scanning

Sagebrush (Artemisia tridentata), a dominant shrub species in the sagebrush-steppe ecosystem of the western US, is declining from its historical distribution due to feedbacks between climate and land use change, fire, and invasive species. Quantifying aboveground biomass of sagebrush is important for assessing carbon storage and monitoring the presence and distribution of this rapidly changing dryland ecosystem. Models of shrub canopy volume, derived from terrestrial laser scanning (TLS) point clouds, were used to accurately estimate aboveground sagebrush biomass. Ninety-one sagebrush plants were scanned and sampled across three study sites in the Great Basin, USA. Half of the plants were scanned and destructively sampled in the spring (n = 46), while the other half were scanned again in the fall before destructive sampling (n = 45). The latter set of sagebrush plants was scanned during both spring and fall to further test the ability of the TLS to quantify seasonal changes in green biomass. Sagebrush biomass was estimated using both a voxel and a 3-D convex hull approach applied to TLS point cloud data. The 3-D convex hull model estimated total and green biomass more accurately (R² = 0.92 and R² = 0.83, respectively) than the voxel-based method (R² = 0.86 and R² = 0.73, respectively). Seasonal differences in TLS-predicted green biomass were detected at two of the sites (p < 0.001 and p = 0.029), elucidating the amount of ephemeral leaf loss in the face of summer drought. The methods presented herein are directly transferable to other dryland shrubs, and implementation of the convex hull model with similar sagebrush species is straightforward.     

Testing the discrimination and detection limits of WorldView-2 imagery on a challenging invasive plant target

Invasive plants pose significant threats to biodiversity and ecosystem function globally, leading to costly monitoring and management effort. While remote sensing promises cost-effective, robust and repeatable monitoring tools to support intervention, it has been largely restricted to airborne platforms that have higher spatial and spectral resolutions, but which lack the coverage and versatility of satellite-based platforms. This study tests the ability of the WorldView-2 (WV2) eight-band satellite sensor for detecting the invasive shrub mesquite (Prosopis spp.) in the north-west Pilbara region of Australia. Detectability was challenged by the target taxa being largely defoliated by a leaf-tying biological control agent (Gelechiidae: Evippe sp. #1) and the presence of other shrubs and trees. Variable importance in the projection (VIP) scores identified bands offering greatest capacity for discrimination were those covering the near-infrared, red, and red-edge wavelengths. Wavelengths between 400 nm and 630 nm (coastal blue, blue, green, yellow) were not useful for species level discrimination in this case. Classification accuracy was tested on three band sets (simulated standard multispectral, all bands, and bands with VIP scores ≥1). Overall accuracies were comparable amongst all band-sets (Kappa = 0.71–0.77). However, mesquite omission rates were unacceptably high (21.3%) when using all eight bands relative to the simulated standard multispectral band-set (9.5%) and the band-set informed by VIP scores (11.9%). An incremental cover evaluation on the latter identified most omissions to be for objects <16 m². Mesquite omissions reduced to 2.6% and overall accuracy significantly improved (Kappa = 0.88) when these objects were left out of the confusion matrix calculations. Very high mapping accuracy of objects >16 m² allows application for mapping mesquite shrubs and coalesced stands, the former not previously possible, even with 3 m resolution hyperspectral imagery. WV2 imagery offers excellent portability potential for detecting other species where spectral/spatial resolution or coverage has been an impediment. New generation satellite sensors are removing barriers previously preventing widespread adoption of remote sensing technologies in natural resource management.     

Retrieval of leaf area index in alpine wetlands using a two-layer canopy reflectance model

In this paper, we focused on the retrieval of the LAI in an alpine wetland located in western part of China in late August and early July 2011. A two-layer canopy reflectance model (ACRM) was used to establish the relationships between the LAI and the reflectance of near-infrared (NIR) and red (RED) wavebands. The reflectance data were derived from Landsat TM L1T product and the Terra and Aqua MODIS 16-day and 8-day composite reflectance products (MOD/MYD09) at 250 m resolution. Due to the lack of the information about some major input parameters for ACRM, which are sensitive to model outputs in the reflectance of NIR and RED wavebands, the inverse problem was ill-posed. To overcome this problem, a method of increasing the sensitivity of the LAI while reducing the influence of other model free parameters based on the study of free parameters’ sensitivity to the ACRM outputs and the region’s features was studied. The area of interest was divided into two parts using the approximately statistic normalized difference vegetation index (NDVI) value around 0.5. One part was sparse vegetation (0.1 < NDVI < 0.5), which is more sensitive to soil background effects and less sensitive to the canopy biophysical and biochemical variables. The other part was dense vegetation (0.5 ≤ NDVI < 1.0), which is less sensitive to soil background effects and more sensitive to plant canopies and leaf parameters. Then, the relationships of ρ_nir–LAI and ρ_red–LAI were established using a look-up table algorithm for the two parts. Furthermore, a regularization technique for fast pixel-wise retrieval was introduced to reduce the elements of LUT sets while maintaining a relatively high accuracy. The results were very promising compared to the field measured LAI values that the correlation (R²) of the measured LAI values and retrieved LAI values reached 0.95, and the root-mean-square deviation (RMSD) was 0.33 for late August, 2011, while the R² reached 0.82 and RMSD was 0.25 for early July 2011.     

Very high spatial resolution optical and radar imagery in tracking water level fluctuations of a small inland reservoir

Tracking water level fluctuations in small lakes and reservoirs is important in order to better understand and manage these ecosystems. A geographic object-based image analysis (GEOBIA) method using very high spatial and temporal resolution optical (Pléiades) and radar (COSMO-SkyMed and TerraSAR-X) remote sensing imagery is presented here which (1) tracks water level fluctuations via variations in water surface area and (2) avoids common difficulties found in using single-band radar images for water-land image classification. Results are robust, with over 98% of image surface area correctly classified into land or water, R² = 0.963 and RMSE = 0.42 m for a total water level fluctuation range of 5.94 m. Multispectral optical imagery is found to be more straightforward in producing results than single-band radar imagery, but the latter crucially increase temporal resolution to the point where fluctuations can be satisfactorily tracked in time. Moreover, an analysis suggests that high and medium spatial resolution imagery is sufficient, in at least some cases, in tracking the water level fluctuations of small inland reservoirs. Finally, limitations of the methodology presented here are briefly discussed along with potential solutions to overcome them.     

Modeling soil parameters using hyperspectral image reflectance in subtropical coastal wetlands

Developing spectral models of soil properties is an important frontier in remote sensing and soil science. Several studies have focused on modeling soil properties such as total pools of soil organic matter and carbon in bare soils. We extended this effort to model soil parameters in areas densely covered with coastal vegetation. Moreover, we investigated soil properties indicative of soil functions such as nutrient and organic matter turnover and storage. These properties include the partitioning of mineral and organic soil between particulate (>53 μm) and fine size classes, and the partitioning of soil carbon and nitrogen pools between stable and labile fractions. Soil samples were obtained from Avicennia germinans mangrove forest and Juncus roemerianus salt marsh plots on the west coast of central Florida. Spectra corresponding to field plot locations from Hyperion hyperspectral image were extracted and analyzed. The spectral information was regressed against the soil variables to determine the best single bands and optimal band combinations for the simple ratio (SR) and normalized difference index (NDI) indices. The regression analysis yielded levels of correlation for soil variables with R² values ranging from 0.21 to 0.47 for best individual bands, 0.28 to 0.81 for two-band indices, and 0.53 to 0.96 for partial least-squares (PLS) regressions for the Hyperion image data. Spectral models using Hyperion data adequately (RPD > 1.4) predicted particulate organic matter (POM), silt + clay, labile carbon (C), and labile nitrogen (N) (where RPD = ratio of standard deviation to root mean square error of cross-validation [RMSECV]). The SR (0.53 μm, 2.11 μm) model of labile N with R² = 0.81, RMSECV= 0.28, and RPD = 1.94 produced the best results in this study. Our results provide optimism that remote-sensing spectral models can successfully predict soil properties indicative of ecosystem nutrient and organic matter turnover and storage, and do so in areas with dense canopy cover.     

Using remote sensing to monitor the influence of river discharge on watershed outlets and adjacent coral Reefs: Magdalena River and Rosario Islands,Colombia

Worldwide, coral reef ecosystems are being increasingly threatened by sediments loads from river discharges, which in turn are influenced by changing rainfall patterns due to climate change and by growing human activity in their watersheds. In this case study, we explored the applicability of using remote sensing (RS) technology to estimate and monitor the relationship between water quality at the coral reefs around the Rosario Islands, in the Caribbean Sea, and the rainfall patterns in the Magdalena River watershed. From the Moderate Resolution Imaging Spectroradiometer (MODIS), this study used the water surface reflectance product (MOD09GQ) to estimate water surface reflectance as a proxy for sediment concentration and the land cover product (MCD12Q1 V51) to characterize land cover of the watershed. Rainfall was estimated by using the 3B43 V7 product from the Tropical Rainforest Measuring Mission (TRMM). For the first trimester of each year, we investigated the inter-annual temporal variation in water surface reflectance at the Rosario Islands and at the three main mouths of the Magdalena River watershed. No increasing or decreasing trends of water surface reflectance were detected for any of the sites for the study period 2001–2014 (p > 0.05) but significant correlations were detected among the trends of each site at the watershed mouths (r = 0.57–0.90, p < 0.05) and between them and the inter-annual variation in rainfall on the watershed (r = 0.63–0.67, p < 0.05). Those trimesters with above-normal water surface reflectance at the mouths and above-normal rainfall at the watershed coincided with La Niña conditions while the opposite was the case during El Niño conditions. Although, a preliminary analysis of inter-annual land cover trends found only cropland cover in the watershed to be significantly correlated with water surface reflectance at two of the watershed mouths (r = 0.58 and 0.63, p < 0.05), the validation analysis draw only a 40.7% of accuracy in this land cover classification. This requires further analysis to confirm the impact of the cropland on the water quality at the watershed outlets. Spatial analysis with MOD09GQ imagery detected the overpass of river plumes from Barbacoas Bay over the Rosario Islands waters.