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.     

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.     

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.     

Predicting grain protein content of winter wheat using remote sensing data based on nitrogen status and water stress

Advanced site-specific knowledge of grain protein content of winter wheat from remote sensing data would provide opportunities to manage grain harvest differently, and to maximize output by adjusting input in fields. In this study, remote sensing data were utilized to predict grain protein content. Firstly, the leaf nitrogen content at winter wheat anthesis stage was proved to be significantly correlated with grain protein content (R² = 0.36), and spectral indices significantly correlated to leaf nitrogen content at anthesis stage were potential indicators for grain protein content. The vegetation index, VI_green, derived from the canopy spectral reflectance at green and red bands, was significantly correlated to the leaf nitrogen content at anthesis stage, and also highly significantly correlated to the final grain protein content (R² = 0.46). Secondly, the external conditions, such as irrigation, fertilization and temperature, had important influence on grain quality. Water stress at grain filling stage can increase grain protein content, and leaf water content is closely related to irrigation levels, therefore, the spectral indices correlated to leaf water content can be potential indicators for grain protein content. The spectral reflectance of TM channel 5 derived from canopy spectra or image data at grain filling stage was all significantly correlated to grain protein content (R² = 0.31 and 0.37, respectively). Finally, not only this study proved the feasibility of using remote sensing data to predict grain protein content, but it also provided a tentative prediction of the grain protein content in Beijing area using the reflectance image of TM channel 5.     

Analysis of common canopy vegetation indices for indicating leaf nitrogen accumulations in wheat and rice

The common spectra wavebands and vegetation indices (VI) were identified for indicating leaf nitrogen accumulation (LNA), and the quantitative relationships of LNA to canopy reflectance spectra were determined in both wheat (Triticum aestivum L.) and rice (Oryza sativa L.). The 810 and 870 nm are two common spectral wavebands indicating LNA in both wheat and rice. Among all ratio vegetation indices (RVI), difference vegetation indices (DVI) and normalized difference vegetation indices (NDVI) of 16 wavebands from the MSR16 radiometer, RVI (870, 660) and RVI (810, 660) were most highly correlated to LNA in both wheat and rice. In addition, the relations between VIs and LNA gave better results than relations between single wavebands and LNA in both wheat and rice. Thus LNA in both wheat and rice could be indicated with common VIs, but separate regression equations are better for LNA monitoring.     

Applying different inversion techniques to retrieve stand variables of summer barley with PROSPECT + SAIL

This study describes the retrieval of state variables (LAI, canopy chlorophyll, water and dry matter contents) for summer barley from airborne HyMap data by means of a canopy reflectance model (PROSPECT + SAIL). Three different inversion techniques were applied to explore the impact of the employed method on estimation accuracies: numerical optimization (downhill simplex method), a look-up table (LUT) and an artificial neural network (ANN) approach. By numerical optimization (Num Opt), reliable estimates were obtained for LAI and canopy chlorophyll contents (LAI × C_ab) with r² of 0.85 and 0.94 and RDP values of 1.81 and 2.65, respectively. Accuracies dropped for canopy water (LAI × C_w) and dry matter contents (LAI × C_m). Nevertheless, the range of leaf water contents (C_w) was very narrow in the studied plant material. Prediction accuracies generally decreased in the order Num Opt > LUT > ANN. This decrease in accuracy mainly resulted from an increase in offset in the obtained values, as the retrievals from the different approaches were highly correlated. The same decreasing order in accuracy was found for the difference between the measured spectra and those reconstructed from the retrieved variable values. The parallel application of the different inversion techniques to one collective data set was helpful to identify modelling uncertainties, as shortcomings of the retrieval algorithms themselves could be separated from uncertainties in model structure and parameterisation schemes.     

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.     

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.     

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.     

Canopy foliar nitrogen retrieved from airborne hyperspectral imagery by correcting for canopy structure effects

A statistical relationship between canopy mass-based foliar nitrogen concentration (%N) and canopy bidirectional reflectance factor (BRF) has been repeatedly demonstrated. However, the interaction between leaf properties and canopy structure confounds the estimation of foliar nitrogen. The canopy scattering coefficient (the ratio of BRF and the directional area scattering factor, DASF) has recently been suggested for estimating %N as it suppresses the canopy structural effects on BRF. However, estimation of %N using the scattering coefficient has not yet been investigated for longer spectral wavelengths (>855 nm). We retrieved the canopy scattering coefficient for wavelengths between 400 and 2500 nm from airborne hyperspectral imagery, and then applied a continuous wavelet analysis (CWA) to the scattering coefficient in order to estimate %N. Predictions of %N were also made using partial least squares regression (PLSR). We found that %N can be accurately retrieved using CWA (R² = 0.65, RMSE = 0.33) when four wavelet features are combined, with CWA yielding a more accurate estimation than PLSR (R² = 0.47, RMSE = 0.41). We also found that the wavelet features most sensitive to %N variation in the visible region relate to chlorophyll absorption, while wavelet features in the shortwave infrared regions relate to protein and dry matter absorption. Our results confirm that %N can be retrieved using the scattering coefficient after correcting for canopy structural effect. With the aid of high-fidelity airborne or upcoming space-borne hyperspectral imagery, large-scale foliar nitrogen maps can be generated to improve the modeling of ecosystem processes as well as ecosystem-climate feedbacks.     

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.     

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.     

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.     

A visible band index for remote sensing leaf chlorophyll content at the canopy scale

Leaf chlorophyll content is an important variable for agricultural remote sensing because of its close relationship to leaf nitrogen content. The triangular greenness index (TGI) was developed based on the area of a triangle surrounding the spectral features of chlorophyll with points at (670 nm, R₆₇₀), (550 nm, R₅₅₀), and (480 nm, R₄₈₀), where R_λ is the spectral reflectance at wavelengths of 670, 550 and 480, respectively. The equation is TGI = −0.5[(670 − 480)(R₆₇₀ − R₅₅₀) − (670 − 550)(R₆₇₀ − R₄₈₀)]. In 1999, investigators funded by NASA's Earth Observations Commercialization and Applications Program collaborated on a nitrogen fertilization experiment with irrigated maize in Nebraska. Airborne Visible/Infrared Imaging Spectrometer (AVIRIS) data and Landsat 5 Thematic Mapper (TM) data were acquired along with leaf chlorophyll meter and other data on three dates in July during late vegetative growth and early reproductive growth. TGI was consistently correlated with plot-averaged chlorophyll-meter values at the spectral resolutions of AVIRIS, Landsat TM, and digital cameras. Simulations using the Scattering by Arbitrarily Inclined Leaves (SAIL) canopy model indicate an interaction among TGI, leaf area index (LAI) and soil type at low crop LAI, whereas at high LAI and canopy closure, TGI was only affected by leaf chlorophyll content. Therefore, TGI may be the best spectral index to detect crop nitrogen requirements with low-cost digital cameras mounted on low-altitude airborne platforms.     

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).     

Deriving leaf mass per area (LMA) from foliar reflectance across a variety of plant species using continuous wavelet analysis

Leaf mass per area (LMA), the ratio of leaf dry mass to leaf area, is a trait of central importance to the understanding of plant light capture and carbon gain. It can be estimated from leaf reflectance spectroscopy in the infrared region, by making use of information about the absorption features of dry matter. This study reports on the application of continuous wavelet analysis (CWA) to the estimation of LMA across a wide range of plant species. We compiled a large database of leaf reflectance spectra acquired within the framework of three independent measurement campaigns (ANGERS, LOPEX and PANAMA) and generated a simulated database using the PROSPECT leaf optical properties model. CWA was applied to the measured and simulated databases to extract wavelet features that correlate with LMA. These features were assessed in terms of predictive capability and robustness while transferring predictive models from the simulated database to the measured database. The assessment was also conducted with two existing spectral indices, namely the Normalized Dry Matter Index (NDMI) and the Normalized Difference index for LMA (ND_LMA).Five common wavelet features were determined from the two databases, which showed significant correlations with LMA (R²: 0.51–0.82, p < 0.0001). The best robustness (R² = 0.74, RMSE = 18.97 g/m² and Bias = 0.12 g/m²) was obtained using a combination of two low-scale features (1639 nm, scale 4) and (2133 nm, scale 5), the first being predominantly important. The transferability of the wavelet-based predictive model to the whole measured database was either better than or comparable to those based on spectral indices. Additionally, only the wavelet-based model showed consistent predictive capabilities among the three measured data sets. In comparison, the models based on spectral indices were sensitive to site-specific data sets. Integrating the ND_LMA spectral index and the two robust wavelet features improved the LMA prediction. One of the bands used by this spectral index, 1368 nm, was located in a strong atmospheric water absorption region and replacing it with the next available band (1340 nm) led to lower predictive accuracies. However, the two wavelet features were not affected by data quality in the atmospheric absorption regions and therefore showed potential for canopy-level investigations. The wavelet approach provides a different perspective into spectral responses to LMA variation than the traditional spectral indices and holds greater promise for implementation with airborne or spaceborne imaging spectroscopy data for mapping canopy foliar dry biomass.     

Applications of AnnAGNPS model for soil loss estimation and nutrient loading for Malaysian conditions

The study evaluated the performance and suitability of AnnAGNPS model in assessing runoff, sediment loading and nutrient loading under Malaysian conditions. The watershed of River Kuala Tasik in Malaysia, a combination of two sub-watersheds, was selected as the area of study. The data for the year 2004 was used to calibrate the model and the data for the year 2005 was used for validation purposes. Several input parameters were computed using methods suggested by other researchers and studies carried out in Malaysia. The study shows that runoff was predicted well with an overall R² value of 0.90 and E value of 0.70. Sediment loading was able to produce a moderate result of R² = 0.66 and E = 0.49, nitrogen loading predictions were slightly better with R² = 0.68 and E = 0.53, and phosphorus loading performance was slightly poor with an R² = 0.63 and E = 0.33. The erosion map developed was in agreement with the erosion risk map produced by the Department of Agriculture, Malaysia. Rubber estates and urban areas were found to be the main contributors to soil erosion. The simulation results showed that AnnAGNPS has the potential to be used as a valuable tool for planning and management of watersheds under Malaysian conditions.     

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.     

Using spectral information from the NIR water absorption features for the retrieval of canopy water content

Canopy water content (CWC) is important for mapping and monitoring the condition of the terrestrial ecosystem. Spectral information related to the water absorption features at 970 nm and 1200 nm offers possibilities for deriving information on CWC. In this study, we compare the use of derivative spectra, spectral indices and continuum removal techniques for these regions. Hyperspectral reflectance data representing a range of canopies were simulated using the combined PROSPECT + SAILH model. Best results in estimating CWC were obtained by using spectral derivatives at the slopes of the 970 nm and 1200 nm water absorption features. Real data from two different test sites were analysed. Spectral information at both test sites was obtained with an ASD FieldSpec spectrometer, whereas at the second site HyMap airborne imaging spectrometer data were also acquired. Best results were obtained for the derivative spectra. In order to avoid the potential influence of atmospheric water vapour absorption bands the derivative of the reflectance on the right slope of the canopy water absorption feature at 970 nm can best be used for estimating CWC.