Estimating vegetation parameter for soil erosion assessment in an alpine catchment by means of QuickBird imagery

Soil erosion rates in alpine regions are related to high spatial variability complicating assessment of risk and damages. A crucial parameter triggering soil erosion that can be derived from satellite imagery is fractional vegetation cover (FVC). The objective of this study is to assess the applicability of normalized differenced vegetation index (NDVI), linear spectral unmixing (LSU) and mixture tuned matched filtering (MTMF) in estimating abundance of vegetation cover in alpine terrain. To account for the small scale heterogeneity of the alpine landscape we used high resolved multispectral QuickBird imagery (pixel resolution = 2.4 m) of a site in the Urseren Valley, Central Swiss Alps (67 km²). A supervised land-cover classification was applied (total accuracy 93.3%) prior to the analysis in order to stratify the image. The regression between ground truth FVC assessment and NDVI as well as MTMF-derived vegetation abundance was significant (r² = 0.64, r² = 0.71, respectively). Best results were achieved for LSU (r² = 0.85). For both spectral unmixing approaches failed to estimate bare soil abundance (r² = 0.39 for LSU, r² = 0.28 for MTMF) due to the high spectral variability of bare soil at the study site and the low spectral resolution of the QuickBird imagery. The LSU-derived FVC map successfully identified erosion features (e.g. landslides) and areas prone to soil erosion. FVC represents an important but often neglected parameter for soil erosion risk assessment in alpine grasslands.     

Spatial dependency of cholera prevalence on potential cholera reservoirs in an urban area,Kumasi, Ghana

Cholera has been a public health burden in Ghana since the early 1970s. Between 1999 and 2005, a total of 25,636 cases and 620 deaths were officially reported to the WHO. In one of the worst affected urban cities, fecal contamination of surface water is extremely high, and the disease is reported to be prevalent among inhabitants living in close proximity to surface water bodies. Surface runoff from dump sites is a major source of fecal and bacterial contamination of rivers and streams in the study area. This study aims to determine (a) the impacts of surface water contamination on cholera infection and (b) detect and map arbitrary shaped clusters of cholera. A Geographic Information System (GIS) based spatial analysis is used to delineate potential reservoirs of the cholera vibrios; possibly contaminated by surface runoff from open space refuse dumps. Statistical modeling using OLS model reveals a significant negative association between (a) cholera prevalence and proximity to all the potential cholera reservoirs (R² = 0.18, p < 0.001) and (b) cholera prevalence and proximity to upstream potential cholera reservoirs (R² = 0.25, p < 0.001). The inclusion of spatial autoregressive coefficients in the OLS model reveals the dependency of the spatial distribution of cholera prevalence on the spatial neighbors of the communities. A flexible scan statistic identifies a most likely cluster with a higher relative risk (RR = 2.04, p < 0.01) compared with the cluster detected by circular scan statistic (RR = 1.60, p < 0.01). We conclude that surface water pollution through runoff from waste dump sites play a significant role in cholera infection.     

Mapping beech (Fagus sylvatica L.) forest structure with airborne hyperspectral imagery

Estimating forest structural attributes using multispectral remote sensing is challenging because of the saturation of multispectral indices at high canopy cover. The objective of this study was to assess the utility of hyperspectral data in estimating and mapping forest structural parameters including mean diameter-at-breast height (DBH), mean tree height and tree density of a closed canopy beech forest (Fagus sylvatica L.). Airborne HyMap images and data on forest structural attributes were collected from the Majella National Park, Italy in July 2004. The predictive performances of vegetation indices (VI) derived from all possible two-band combinations (VI_(i,j) = (R_i − R_j)/(R_i + R_j), where R_i and R_j = reflectance in any two bands) were evaluated using calibration (n = 33) and test (n = 20) data sets. The potential of partial least squares (PLS) regression, a multivariate technique involving several bands was also assessed. New VIs based on the contrast between reflectance in the red-edge shoulder (756–820 nm) and the water absorption feature centred at 1200 nm (1172–1320 nm) were found to show higher correlations with the forest structural parameters than standard VIs derived from NIR and visible reflectance (i.e. the normalised difference vegetation index, NDVI). PLS regression showed a slight improvement in estimating the beech forest structural attributes (prediction errors of 27.6%, 32.6% and 46.4% for mean DBH, height and tree density, respectively) compared to VIs using linear regression models (prediction errors of 27.8%, 35.8% and 48.3% for mean DBH, height and tree density, respectively). Mean DBH was the best predicted variable among the stand parameters (calibration R² = 0.62 for an exponential model fit and standard error of prediction = 5.12 cm, i.e. 25% of the mean). The predicted map of mean DBH revealed high heterogeneity in the beech forest structure in the study area. The spatial variability of mean DBH occurs at less than 450 m. The DBH map could be useful to forest management in many ways, e.g. thinning of coppice to promote diameter growth, to assess the effects of management on forest structure or to detect changes in the forest structure caused by anthropogenic and natural factors.     

Terrestrial laser scanning to estimate plot-level forest canopy fuel properties

This paper evaluates the potential of a terrestrial laser scanner (TLS) to characterize forest canopy fuel characteristics at plot level. Several canopy properties, namely canopy height, canopy cover, canopy base height and fuel strata gap were estimated. Different approaches were tested to avoid the effect of canopy shadowing on canopy height estimation caused by deployment of the TLS below the canopy. Estimation of canopy height using a grid approach provided a coefficient of determination of R² = 0.81 and an RMSE of 2.47 m. A similar RMSE was obtained using the 99th percentile of the height distribution of the highest points, representing the 1% of the data, although the coefficient of determination was lower (R² = 0.70). Canopy cover (CC) was estimated as a function of the occupied cells of a grid superimposed upon the TLS point clouds. It was found that CC estimates were dependent on the cell size selected, with 3 cm being the optimum resolution for this study. The effect of the zenith view angle on CC estimates was also analyzed. A simple method was developed to estimate canopy base height from the vegetation vertical profiles derived from an occupied/non-occupied voxels approach. Canopy base height was estimated with an RMSE of 3.09 m and an R² = 0.86. Terrestrial laser scanning also provides a unique opportunity to estimate the fuel strata gap (FSG), which has not been previously derived from remotely sensed data. The FSG was also derived from the vegetation vertical profile with an RMSE of 1.53 m and an R² = 0.87.     

Exploring hyperspectral bands and estimation indices for leaf nitrogen accumulation in wheat

Hyperspectral sensing can provide an effective means for fast and non-destructive estimation of leaf nitrogen (N) status in crop plants. The objectives of this study were to design a new method to extract hyperspectral spectrum information, to explore sensitive spectral bands, suitable bandwidth and best vegetation indices based on precise analysis of ground-based hyperspectral information, and to develop regression models for estimating leaf N accumulation per unit soil area (LNA, g N m⁻²) in winter wheat (Triticum aestivum L.). Three field experiments were conducted with different N rates and cultivar types in three consecutive growing seasons, and time-course measurements were taken on canopy hyperspectral reflectance and LNA under the various treatments. Then, normalized difference spectral indices (NDSI) and ratio spectral indices (RSI) based on the original spectrum and the first derivative spectrum were constructed within the range of 350–2500 nm, and their relationships with LNA were quantified. The results showed that both LNA and canopy hyperspectral reflectance in wheat changed with varied N rates, with consistent patterns across different cultivars and seasons. The sensitive spectral bands for LNA existed mainly within visible and near infrared regions. The best spectral indices for estimating LNA in wheat were found to be NDSI (R₈₆₀, R₇₂₀), RSI (R₉₉₀, R₇₂₀), NDSI (FD₇₃₆, FD₅₂₆) and RSI (FD₇₂₅, FD₅₁₆), and the regression models based on the above four spectral indices were formulated as Y = 26.34x^1.887, Y = 5.095x − 6.040, Y = 0.609 e^3.008x and Y = 0.388x^1.260, respectively, with R² greater than 0.81. Furthermore, expanding the bandwidth of NDSI (R₈₆₀, R₇₂₀) and RSI (R₉₉₀, R₇₂₀) from 1 nm to 100 nm at 1 nm interval produced the LNA monitoring models with similar performance within about 33 nm and 23 nm bandwidth, respectively, over which the statistical parameters of the models became less stable. From testing of the derived equations, the model for LNA estimation on NDSI (R₈₆₀, R₇₂₀), RSI (R₉₉₀, R₇₂₀), NDSI (FD₇₃₆, FD₅₂₆) and RSI (FD₇₂₅, FD₅₁₆) gave R² over 0.79 with more satisfactory performance than previously reported models and physical models in wheat. It can be concluded that the present hyperspectral parameters of NDSI (R₈₆₀, R₇₂₀), RSI (R₉₉₀, R₇₂₀), NDSI (FD₇₃₆, FD₅₂₆) and RSI (FD₇₂₅, FD₅₁₆) can be reliably used for estimating LNA in winter wheat.     

Retrieval of chlorophyll and nitrogen in Norway spruce (Picea abies L. Karst.) using imaging spectroscopy

The research evaluated the information content of spectral reflectance (laboratory and airborne data) for the estimation of needle chlorophyll (C_AB) and nitrogen (C_N) concentration in Norway spruce (Picea abies L. Karst.) needles. To identify reliable predictive models different types of spectral transformations were systematically compared regarding the accuracy of prediction. The results of the cross-validated analysis showed that C_AB can be well estimated from laboratory and canopy reflectance data. The best predictive model to estimate C_AB was achieved from laboratory spectra using continuum-removal transformed data (R²cv = 0.83 and a relative RMSEcv of 8.1%, n = 78) and from hyperspectral HyMap data using band-depth normalised spectra (R²cv = 0.90, relative RMSEcv = 2.8%, n = 13). Concerning the nitrogen concentration, we observed somewhat weaker relations, with however still acceptable accuracies (at canopy level: R²cv = 0.57, relative RMSEcv = 4.6%). The wavebands selected in the regression models to estimate C_AB were typically located in the red edge region and near the green reflectance peak. For C_N, additional wavebands related to a known protein absorption feature at 2350 nm were selected. The portion of selected wavebands attributable to known absorption features strongly depends on the type of spectral transformation applied. A method called “water removal” (WR) produced for canopy spectra the largest percentage of wavebands directly or indirectly related to known absorption features. The derived chlorophyll and nitrogen maps may support the detection and the monitoring of environmental stressors and are also important inputs to many bio-geochemical process models.     

Spectral mixture analysis to monitor defoliation in mixed-aged Eucalyptus globulus Labill plantations in southern Australia using Landsat5-TM and EO-1Hyperion data

Defoliation is a key parameter of forest health and is associated with reduced productivity and tree mortality. Assessing the health of forests requires regular observations over large areas. Satellite remote sensing provides a cost-effective alternative to traditional ground-based assessment of forest health, but assessing defoliation can be difficult due to mixed pixels where vegetation cover is low or fragmented. In this study we apply a novel spectral unmixing technique, referred to as weighted Multiple Endmember Spectral Mixture Analysis (wMESMA), to Landsat 5-TM and EO-1 Hyperion data acquired over a Eucalyptus globulus (Labill.) plantation in southern Australia. This technique combines an iterative mixture analysis cycle allowing endmembers to vary on a per pixel basis (MESMA) and a weighting algorithm that prioritizes wavebands based on their robustness against endmember variability. Spectral mixture analysis provides an estimate of the physically interpretable canopy cover, which is not necessarily correlated with defoliation in mixed-aged plantations due to natural variation in canopy cover as stands age. There is considerable variability in the degree of defoliation as well as in stand age among sites and in this study we found that results were significantly improved by the inclusion of an age correction algorithm for both the multi-spectral (R²_{no age correction} = 0.55 vs R²_{age correction} = 0.73 for Landsat) and hyperspectral (R²_{no age correction} = 0.12 vs R²_{age correction} = 0.50 for Hyperion) image data. The improved accuracy obtained from Landsat compared to the Hyperion data illustrates the potential of applying SMA techniques for analysis of multi-spectral datasets such as MODIS and SPOT-VEGETATION.     

Predicting forest structural parameters using the image texture derived from WorldView-2 multispectral imagery in a dryland forest,Israel

Estimation of forest structural parameters by field-based data collection methods is both expensive and time consuming. Satellite remote sensing is a low-cost alternative in modeling and mapping structural parameters in large forest areas. The current study investigates the potential of using WordView-2 multispectral satellite imagery for predicting forest structural parameters in a dryland plantation forest in Israel. The relationships between image texture features and the several structural parameters such as Number of Trees (NT), Basal Area (BA), Stem Volume (SV), Clark-Evans Index (CEI), Diameter Differentiation Index (DDI), Contagion Index (CI), Gini Coefficient (GC), and Standard Deviation of Diameters at Breast Heights (SDDBH) were examined using correlation analyses. These variables were obtained from 30 m × 30 m square-shaped plots. The Standard Deviation of Gray Levels (SDGL) as a first order texture feature and the second order texture variables based on Gray Level Co-occurrence Matrix (GLCM) were calculated for the pixels that corresponds to field plots. The results of the correlation analysis indicate that the forest structural parameters are significantly correlated with the image texture features. The highest correlation coefficients were calculated for the relationships between the SDDBH and the contrast of red band (r = 0.75, p < 0.01), the BA and the entropy of blue band (r = 0.73, p < 0.01), and the GC and the contrast of blue band (r = 0.71, p < 0.01). Each forest structural parameter was modeled as a function of texture measures derived from the satellite image using stepwise multi linear regression analyses. The determination coefficient (R²) and root mean square error (RMSE) values of the best fitting models, respectively, are 0.38 and 109.56 ha⁻¹ for the NT; 0.54 and 1.79 m² ha⁻¹ for the BA; 0.42 and 27.18 m³ ha⁻¹ for the SV; 0.23 and 0.16 for the CEI; 0.32 and 0.05 for the DDI; 0.25 and 0.06 for the CI; 0.50 and 0.05 for the GC; and 0.67 and 0.70 for the SDDBH. The leave-one-out cross-validation technique was applied for validation of the best-fitted models (R² > 0.50). In conclusion, cross-validated statistics confirmed that the structural parameters including the BA, SDDBH, and GC can be predicted and mapped with a reasonable accuracy using the texture features extracted from the spectral bands of WorldView-2 image.     

Spectral mixture analysis to assess post-fire vegetation regeneration using Landsat Thematic Mapper imagery: Accounting for soil brightness variation

Post-fire vegetation cover is a crucial parameter in rangeland management. This study aims to assess the post-fire vegetation recovery 3 years after the large 2007 Peloponnese (Greece) wildfires. Post-fire recovery landscapes typically are mixed vegetation-substrate environments which makes spectral mixture analysis (SMA) a very effective tool to derive fractional vegetation cover maps. Using a combination of field and simulation techniques this study aimed to account for the impact of background brightness variability on SMA model performance. The field data consisted out of a spectral library of in situ measured reflectance signals of vegetation and substrate and 78 line transect plots. In addition, a Landsat Thematic Mapper (TM) scene was employed in the study. A simple SMA, in which each constituting terrain feature is represented by its mean spectral signature, a multiple endmember SMA (MESMA) and a segmented SMA, which accounts for soil brightness variations by forcing the substrate endmember choice based on ancillary data (lithological map), were applied. In the study area two main spectrally different lithological units were present: relatively bright limestone and relatively dark flysch (sand-siltstone). Although the simple SMA model resulted in reasonable regression fits for the flysch and limestones subsets separately (coefficient of determination R² of respectively 0.67 and 0.72 between field and TM data), the performance of the regression model on the pooled dataset was considerably weaker (R² = 0.65). Moreover, the regression lines significantly diverged among the different subsets leading to systematic over-or underestimations of the vegetative fraction depending on the substrate type. MESMA did not solve the endmember variability issue. The MESMA model did not manage to select the proper substrate spectrum on a reliable basis due to the lack of shape differences between the flysch and limestone spectra,. The segmented SMA model which accounts for soil brightness variations minimized the variability problems. Compared to the simple SMA and MESMA models, the segmented SMA resulted in a higher overall correlation (R² = 0.70), its regression slope and intercept were more similar among the different substrate types and its resulting regression lines more closely resembled the expected one-one line. This paper demonstrates the improvement of a segmented approach in accounting for soil brightness variations in estimating vegetative cover using SMA. However, further research is required to evaluate the model's performance for other soil types, with other image data and at different post-fire timings.     

Generation and evaluation of gross primary productivity using Landsat data through blending with MODIS data

Spatial and Temporal Adaptive Reflectance Fusion Model (STARFM) has been used for the blending of Landsat and MODIS data. Specifically, the 30 m Landsat-7 ETM+ (Enhanced Thematic Mapper plus) surface reflectance was predicted for a period of 10 years (2000–2009) as the product of observed ETM+ and MODIS surface reflectance (MOD09A1) on the predicted and observed ETM+ dates. A pixel based analysis for six observed ETM+ dates covering winter and summer crops showed that the prediction method was more accurate for NIR band (mean r² = 0.71, p ≤ 0.01) compared to green band (mean r² = 0.53; p ≤ 0.01). A recently proposed chlorophyll index (CI), which involves NIR and green spectral bands, was used to retrieve gross primary productivity (GPP) as the product of CI and photosynthetic active radiation (PAR). The regression analysis of GPP derived from closet observed and synthetic ETM+ showed a good agreement (r² = 0.85, p ≤ 0.01 and r² = 0.86, p ≤ 0.01) for wheat and sugarcane crops, respectively. The difference between the GPP derived from synthetic and observed ETM+ (prediction residual) was compared with the difference in GPP values from observed ETM+ on the two dates (temporal residual). The prediction residuals (mean value of 1.97 g C/m² in 8 days) was found to be significantly lower than the temporal residuals (mean value of 4.46 g C/m² in 8 days) that correspondence to 12% and 27%, respectively, of GPP values (mean value of 16.53 g C/m² in 8 days) from observed ETM+ data, implying that the prediction method was better than temporal pixel substitution. Investigating the trend in synthetic ETM+ GPP values over a growing season revealed that phenological patterns were well captured for wheat and sugarcane crops. A direct comparison between the GPP values derived from MODIS and synthetic ETM+ data showed a good consistency of the temporal dynamics but a systematic error that can be read as bias (MODIS GPP over estimation). Further, the regression analysis between observed evapotranspiration and synthetic ETM+ GPP showed good agreement (r² = 0.66, p ≤ 0.01).     

Estimation of gross primary production in wheat from in situ measurements

Gross primary production (GPP) is a parameter of significant importance for carbon cycle and climate change research. Remote sensing combined with other climate and meteorological data offers a convenient tool for large scale GPP estimation. This paper presents a study of GPP estimation using three methods with in situ measurements of canopy reflectance, LAI, and the photosynthetically active radiation (PAR). First, because LAI is considered as an indicator of the factor of absorbed PAR (f_APAR), it provides reasonable estimates of GPP for all types of wheat with coefficient of determination R² of 0.7353. The second method uses four kinds of vegetation indices (VIs) to estimate GPP because these indices are suggested to be reliable candidates in the estimation of light use efficiency (LUE). Good determination coefficients were acquired in estimating GPP with R² ranging from the lowest of 0.7604 for NDVI to the highest of 0.8505 for EVI. A new method was proposed for the estimation of GPP following the Monteith logic, which considering GPP as a product of VI × VI × PAR. Results indicated that this method can provide the best estimates of GPP as determination coefficient R² increased largely compared to the other two methods. EVI × EVI × PAR was demonstrated to be the most suitable for the estimation of GPP with the highest R² of 0.9207, which was about 10% larger as compared to GPP estimated from the single EVI. These results will be helpful for the development of new models of GPP estimation with all remote sensing inputs.     

Estimating sagebrush cover in semi-arid environments using Landsat Thematic Mapper data

Sagebrush ecosystems of the western US provide important habitat for several ungulate and vertebrate species. As a consequence of energy development, these ecosystems in Wyoming have been subjected to a variety of anthropogenic disturbances. Land managers require methodology that will allow them to consistently catalog sagebrush ecosystems and evaluate potential impact of proposed anthropogenic activities. This study addresses the utility of remotely sensed and ancillary geospatial data to estimate sagebrush cover using ordinal logistic regression. We demonstrate statistically significant prediction of ordinal sagebrush cover categories using spectral (χ² = 113; p < 0.0001) and transformed indices (χ² = 117; p < 0.0001). Both Landsat spectral bands (c-value = 0.88) and transformed indices (c-value = 0.89) can distinguish sites with closed, moderate and open cover sagebrush cover categories from no cover. The techniques described in this study can be used for estimating categories of sagebrush cover in arid ecosystems.     

Soil erosion modelled with USLE and PESERA using QuickBird derived vegetation parameters in an alpine catchment

The focus of soil erosion research in the Alps has been in two categories: (i) on-site measurements, which are rather small scale point measurements on selected plots often constrained to irrigation experiments or (ii) off-site quantification of sediment delivery at the outlet of the catchment. Results of both categories pointed towards the importance of an intact vegetation cover to prevent soil loss. With the recent availability of high-resolution satellites such as IKONOS and QuickBird options for detecting and monitoring vegetation parameters in heterogeneous terrain have increased. The aim of this study is to evaluate the usefulness of QuickBird derived vegetation parameters in soil erosion models for alpine sites by comparison to Cesium-137 (Cs-137) derived soil erosion estimates. The study site (67 km²) is located in the Central Swiss Alps (Urseren Valley) and is characterised by scarce forest cover and strong anthropogenic influences due to grassland farming for centuries. A fractional vegetation cover (FVC) map for grassland and detailed land-cover maps are available from linear spectral unmixing and supervised classification of QuickBird imagery. The maps were introduced to the Pan-European Soil Erosion Risk Assessment (PESERA) model as well as to the Universal Soil Loss Equation (USLE). Regarding the latter model, the FVC was indirectly incorporated by adapting the C factor. Both models show an increase in absolute soil erosion values when FVC is considered. In contrast to USLE and the Cs-137 soil erosion rates, PESERA estimates are low. For the USLE model also the spatial patterns improved and showed “hotspots” of high erosion of up to 16 t ha⁻¹ a⁻¹. In conclusion field measurements of Cs-137 confirmed the improvement of soil erosion estimates using the satellite-derived vegetation data.     

Estimating the parameters of forest inventory using machine learning and the reduction of remote sensing features

Locally computed statistics of image texture and a case-based reasoning (CBR) system were evaluated for mapping of forest attributes. Cluster analysis was preferred to regression models, as a pre-selection method of features. The best stand-based accuracy using satellite sensor images was 74.64 m⁻³ ha⁻¹ (36%) RMSE for stand volume, 1.98 m⁻³ ha⁻¹ a⁻¹ (49%) for annual increase in stand volume, where κ = 0.23 for stand growth classes and κ = 0.41 for dominant tree species in stands. The top pixel-based accuracy using orthophotos was 76.54 m⁻³ ha⁻¹ (41%) RMSE for stand volume, 1.87 m⁻³ ha⁻¹ a⁻¹ (44%) for annual increase in stand volume, where κ = 0.24 for stand growth classes and κ = 0.38 for dominant tree species in stands. Mean saturation in 30 m radius was the most useful feature when orthophotos were used, and standard deviation of Landsat ETM 6.2 values in 80 m radius was the best when satellite sensor images were used. The most valuable feature components (radii, channels and local statistics) for orthophotos were: 30 m kernel radius, lightness and the mean of pixel values; for satellite sensor images: 80 m kernel radius, near-infrared channel (ETM 4) and the mean of pixel values. Locally computed statistics.     

A time-integrated MODIS burn severity assessment using the multi-temporal differenced normalized burn ratio (dNBRMT)

Burn severity is an important parameter in post-fire management. It incorporates both the direct fire impact (vegetation depletion) and ecosystem responses (vegetation regeneration). From a remote sensing perspective, burn severity is traditionally estimated using Landsat's differenced normalized burn ratio (dNBR). In this case study of the large 2007 Peloponnese (Greece) wildfires, Landsat dNBR estimates correlated reasonably well with Geo composite burn index (GeoCBI) field data of severity (R² = 0.56). The usage of Landsat imagery is, however, restricted by cloud cover and image-to-image normalization constraints. Therefore a multi-temporal burn severity approach based on coarse spatial, high temporal resolution moderate resolution imaging spectroradiometer (MODIS) imagery is presented in this study. The multi-temporal dNBR (dNBR_MT) is defined as the 1-year integrated difference between burned pixels and their unique control pixels. These control pixels were selected based on time series similarity and spatial context and reflect how burned pixels would have behaved in the case no fire had occurred. Linear regression between downsampled Landsat dNBR and dNBR_MT estimates resulted in a moderate-high coefficient of determination R² = 0.54. dNBR_MT estimates are indicative for the change in vegetation productivity due to the fire. This change is considerably higher for forests than for more sparsely vegetated areas like shrub lands. Although Landsat dNBR is superior for spatial detail, MODIS-derived dNBR_MT estimates present a valuable alternative for burn severity mapping at continental to global scale without image availability constraints. This is beneficial to compare trends in burn severity across regions and time. Moreover, thanks to MODIS's repeated temporal sampling, the dNBR_MT accounts for both first- and second-order fire effects.     

A global map of rainfed cropland areas (GMRCA) at the end of last millennium using remote sensing

The overarching goal of this study was to produce a global map of rainfed cropland areas (GMRCA) and calculate country-by-country rainfed area statistics using remote sensing data. A suite of spatial datasets, methods and protocols for mapping GMRCA were described. These consist of: (a) data fusion and composition of multi-resolution time-series mega-file data-cube (MFDC), (b) image segmentation based on precipitation, temperature, and elevation zones, (c) spectral correlation similarity (SCS), (d) protocols for class identification and labeling through uses of SCS R²-values, bi-spectral plots, space-time spiral curves (ST-SCs), rich source of field-plot data, and zoom-in-views of Google Earth (GE), and (e) techniques for resolving mixed classes by decision tree algorithms, and spatial modeling. The outcome was a 9-class GMRCA from which country-by-country rainfed area statistics were computed for the end of the last millennium. The global rainfed cropland area estimate from the GMRCA 9-class map was 1.13 billion hectares (Bha). The total global cropland areas (rainfed plus irrigated) was 1.53 Bha which was close to national statistics compiled by FAOSTAT (1.51 Bha). The accuracies and errors of GMRCA were assessed using field-plot and Google Earth data points. The accuracy varied between 92 and 98% with kappa value of about 0.76, errors of omission of 2–8%, and the errors of commission of 19–36%.     

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.     

Water-removed spectra increase the retrieval accuracy when estimating savanna grass nitrogen and phosphorus concentrations

Information about the distribution of grass foliar nitrogen (N) and phosphorus (P) is important for understanding rangeland vitality and for facilitating the effective management of wildlife and livestock. Water absorption effects in the near-infrared (NIR) and shortwave-infrared (SWIR) regions pose a challenge for nutrient estimation using remote sensing. The aim of this study was to test the utility of water-removed (WR) spectra in combination with partial least-squares regression (PLSR) and stepwise multiple linear regression (SMLR) to estimate foliar N and P, compared to spectral transformation techniques such as first derivative, continuum removal and log-transformed (Log(1/R)) spectra. The study was based on a greenhouse experiment with a savanna grass species (Digitaria eriantha). Spectral measurements were made using a spectrometer. The D. eriantha was cut, dried and chemically analyzed for foliar N and P concentrations. WR spectra were determined by calculating the residual from the modelled leaf water spectra using a nonlinear spectral matching technique and observed leaf spectra. Results indicated that the WR spectra yielded a higher N retrieval accuracy than a traditional first derivative transformation (R²=0.84, RMSE = 0.28) compared to R²=0.59, RMSE = 0.45 for PLSR. Similar trends were observed for SMLR. The highest P retrieval accuracy was derived from WR spectra using SMLR (R²=0.64, RMSE = 0.067), while the traditional first derivative and continuum removal resulted in lower accuracy. Only when using PLSR did the first derivative result in a higher P retrieval accuracy (R²=0.47, RMSE = 0.07) than the WR spectra (R²=0.43, RMSE = 0.070). It was concluded that the water removal technique is a promising technique to minimize the perturbing effect of foliar water content when estimating grass nutrient concentrations.     

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.