Evaluating suitability of MODIS-Terra images for reproducing historic sediment concentrations in water bodies: Lake Tana,Ethiopia

Government and NGO funded conservation programs are being implemented in developing countries with the potential benefit of reduced sediment inflow into fresh water lakes. However, these claims are difficult to verify due to limited historical sediment concentration data in lakes and rivers. Remote sensing can potentially aid in monitoring sediment concentration. With almost daily availability over the past ten years and consistent atmospheric correction applied to the images, Moderate Resolution Imaging Spectroradiometer (MODIS) 250 meter images are potential resources capable of monitoring future concentrations and reconstructing historical sediment concentration records. In this paper, site-specific relationships are developed between reflectance in near-infrared (NIR) images and three factors: total suspended solids (TSS), turbidity and Secchi depth for Lake Tana near the mouth of the Gumara River. The first two sampling campaigns on November 27, 2010 and May 13, 2011 are used in calibration. Reflectance in the NIR varies linearly with turbidity (R² = 0.89) and TSS (R² = 0.95). Secchi depth fit best to an exponential relation with R² of 0.74. The relationships are validated using a third sample set collected on November 7, 2011 with RMSE of 11 Nephelometric Turbidity Units (NTU) for Turbidity, 16.5 mg l⁻¹ for TSS and 0.12 meters for Secchi depth. The MAE was 10% for TSS, 14% for turbidity and 0.1% for Secchi depth. Using the relationship for TSS, a 10-year time series of sediment concentration in Lake Tana near the Gumara River was plotted. It was found that after the severe drought of 2002 and 2003 the concentration in the lake increased significantly. The results showed that MODIS images are potential cost effective tools to monitor suspended sediment concentration and obtain a past history of concentration for evaluating the effect of best management practices.     

Predicting stem borer density in maize using RapidEye data and generalized linear models

Average maize yield in eastern Africa is 2.03 t ha⁻¹ as compared to global average of 6.06 t ha⁻¹ due to biotic and abiotic constraints. Amongst the biotic production constraints in Africa, stem borers are the most injurious. In eastern Africa, maize yield losses due to stem borers are currently estimated between 12% and 21% of the total production. The objective of the present study was to explore the possibility of RapidEye spectral data to assess stem borer larva densities in maize fields in two study sites in Kenya. RapidEye images were acquired for the Bomet (western Kenya) test site on the 9th of December 2014 and on 27th of January 2015, and for Machakos (eastern Kenya) a RapidEye image was acquired on the 3rd of January 2015. Five RapidEye spectral bands as well as 30 spectral vegetation indices (SVIs) were utilized to predict per field maize stem borer larva densities using generalized linear models (GLMs), assuming Poisson (‘Po’) and negative binomial (‘NB’) distributions. Root mean square error (RMSE) and ratio prediction to deviation (RPD) statistics were used to assess the models performance using a leave-one-out cross-validation approach. The Zero-inflated NB (‘ZINB’) models outperformed the ‘NB’ models and stem borer larva densities could only be predicted during the mid growing season in December and early January in both study sites, respectively (RMSE = 0.69–1.06 and RPD = 8.25–19.57). Overall, all models performed similar when all the 30 SVIs (non-nested) and only the significant (nested) SVIs were used. The models developed could improve decision making regarding controlling maize stem borers within integrated pest management (IPM) interventions.     

Visualizing the ill-posedness of the inversion of a canopy radiative transfer model: A case study for Sentinel-2

Monitoring biophysical and biochemical vegetation variables in space and time is key to understand the earth system. Operational approaches using remote sensing imagery rely on the inversion of radiative transfer models, which describe the interactions between light and vegetation canopies. The inversion required to estimate vegetation variables is, however, an ill-posed problem because of variable compensation effects that can cause different combinations of soil and canopy variables to yield extremely similar spectral responses. In this contribution, we present a novel approach to visualise the ill-posed problem using self-organizing maps (SOM), which are a type of unsupervised neural network. The approach is demonstrated with simulations for Sentinel-2 data (13 bands) made with the Soil-Leaf-Canopy (SLC) radiative transfer model. A look-up table of 100,000 entries was built by randomly sampling 14 SLC model input variables between their minimum and maximum allowed values while using both a dark and a bright soil. The Sentinel-2 spectral simulations were used to train a SOM of 200 × 125 neurons. The training projected similar spectral signatures onto either the same, or contiguous, neuron(s). Tracing back the inputs that generated each spectral signature, we created a 200 × 125 map for each of the SLC variables. The lack of spatial patterns and the variability in these maps indicate ill-posed situations, where similar spectral signatures correspond to different canopy variables. For Sentinel-2, our results showed that leaf area index, crown cover and leaf chlorophyll, water and brown pigment content are less confused in the inversion than variables with noisier maps like fraction of brown canopy area, leaf dry matter content and the PROSPECT mesophyll parameter. This study supports both educational and on-going research activities on inversion algorithms and might be useful to evaluate the uncertainties of retrieved canopy biophysical and biochemical state variables.     

Assessing the potential of MODIS/Terra version 5 images to improve near shore lake bathymetric surveys

Estimation of storage volume of inland water bodies is made using bathymetric surveys. Such surveys are often terminated far inward from shores due to operational limitations. The gap in depth data for the shore area is often overcome by extrapolating the established volume-area relationship beyond the survey boundary. This introduces large errors in volume estimation. In this study, the potential of the Moderate Resolution Imaging Spectroradiometer MODIS version-5 images to improve an existing bathymetric model is assessed over Lake Tana, Ethiopia. MODIS images are used to generate a new bathymetric model of the unsurveyed shore. The new model is validated through its ability to reproduce measured lake levels. The result is compared with lake level estimates made from an existing bathymetric survey. Results showed that MODIS improved bathymetric model reproduced lake water level with a RMSE of 0.14 m as compared to 0.82 m using the existing bathymetric model. The Mean Absolute Error (MAE) has also improved from 0.78 m to 0.13 m. These results suggest that despite their unavailability on cloudy days, MODIS images provide significant advantages for improving bathymetric models.     

Comparative analysis of different retrieval methods for mapping grassland leaf area index using airborne imaging spectroscopy

Fine scale maps of vegetation biophysical variables are useful status indicators for monitoring and managing national parks and endangered habitats. Here, we assess in a comparative way four different retrieval methods for estimating leaf area index (LAI) in grassland: two radiative transfer model (RTM) inversion methods (one based on look-up-tables (LUT) and one based on predictive equations) and two statistical modelling methods (one partly, the other entirely based on in situ data). For prediction, spectral data were used that had been acquired over Majella National Park in Italy by the airborne hyperspectral HyMap instrument. To assess the performance of the four investigated models, the normalized root mean squared error (nRMSE) and coefficient of determination (R²) between estimates and in situ LAI measurements are reported (n = 41). Using a jackknife approach, we also quantified the accuracy and robustness of empirical models as a function of the size of the available calibration data set. The results of the study demonstrate that the LUT-based RTM inversion yields higher accuracies for LAI estimation (R² = 0.91, nRMSE = 0.18) as compared to RTM inversions based on predictive equations (R² = 0.79, nRMSE = 0.38). The two statistical methods yield accuracies similar to the LUT method. However, as expected, the accuracy and robustness of the statistical models decrease when the size of the calibration database is reduced to fewer samples. The results of this study are of interest for the remote sensing community developing improved inversion schemes for spaceborne hyperspectral sensors applicable to different vegetation types. The examples provided in this paper may also serve as illustrations for the drawbacks and advantages of physical and empirical models.     

Aboveground biomass mapping in French Guiana by combining remote sensing,forest inventories and environmental data

Mapping forest aboveground biomass (AGB) has become an important task, particularly for the reporting of carbon stocks and changes. AGB can be mapped using synthetic aperture radar data (SAR) or passive optical data. However, these data are insensitive to high AGB levels (>150 Mg/ha, and >300 Mg/ha for P-band), which are commonly found in tropical forests. Studies have mapped the rough variations in AGB by combining optical and environmental data at regional and global scales. Nevertheless, these maps cannot represent local variations in AGB in tropical forests. In this paper, we hypothesize that the problem of misrepresenting local variations in AGB and AGB estimation with good precision occurs because of both methodological limits (signal saturation or dilution bias) and a lack of adequate calibration data in this range of AGB values. We test this hypothesis by developing a calibrated regression model to predict variations in high AGB values (mean >300 Mg/ha) in French Guiana by a methodological approach for spatial extrapolation with data from the optical geoscience laser altimeter system (GLAS), forest inventories, radar, optics, and environmental variables for spatial inter- and extrapolation. Given their higher point count, GLAS data allow a wider coverage of AGB values. We find that the metrics from GLAS footprints are correlated with field AGB estimations (R² = 0.54, RMSE = 48.3 Mg/ha) with no bias for high values. First, predictive models, including remote-sensing, environmental variables and spatial correlation functions, allow us to obtain “wall-to-wall” AGB maps over French Guiana with an RMSE for the in situ AGB estimates of ∼50 Mg/ha and R² = 0.66 at a 1-km grid size. We conclude that a calibrated regression model based on GLAS with dependent environmental data can produce good AGB predictions even for high AGB values if the calibration data fit the AGB range. We also demonstrate that small temporal and spatial mismatches between field data and GLAS footprints are not a problem for regional and global calibrated regression models because field data aim to predict large and deep tendencies in AGB variations from environmental gradients and do not aim to represent high but stochastic and temporally limited variations from forest dynamics. Thus, we advocate including a greater variety of data, even if less precise and shifted, to better represent high AGB values in global models and to improve the fitting of these models for high values.     

Using multi-polarization C- and L-band synthetic aperture radar to estimate biomass and soil moisture of wheat fields

Biomass and soil moisture are two important parameters for agricultural crop monitoring and yield estimation. In this study, the Water Cloud Model (WCM) was coupled with the Ulaby soil moisture model to estimate both biomass and soil moisture for spring wheat fields in a test site in western Canada. This study exploited both C-band (RADARSAT-2) and L-band (UAVSAR) Synthetic Aperture Radars (SARs) for this purpose. The WCM-Ulaby model was calibrated for three polarizations (HH, VV and HV). Subsequently two of these three polarizations were used as inputs to an inversion procedure, to retrieve either soil moisture or biomass without the need for any ancillary data. The model was calibrated for total canopy biomass, the biomass of only the wheat heads, as well as for different wheat growth stages. This resulted in a calibrated WCM-Ulaby model for each sensor-polarization-phenology-biomass combination. Validation of model retrievals led to promising results. RADARSAT-2 (HH-HV) estimated total wheat biomass with root mean square (RMSE) and mean average (MAE) errors of 78.834 g/m² and 58.438 g/m²; soil moisture with errors of 0.078 m³/m³ (RMSE) and 0.065 m³/m³ (MAE) are reported. During the period of crop ripening, L-band estimates of soil moisture had accuracies of 0.064 m³/m³ (RMSE) and 0.057 m³/m³ (MAE). RADARSAT-2 (VV-HV) produced interesting results for retrieval of the biomass of the wheat heads. In this particular case, the biomass of the heads was estimated with accuracies of 38.757 g/m² (RSME) and 33.152 g/m² (MAE). For wider implementation this model will require additional data to strengthen the model accuracy and confirm estimation performance. Nevertheless this study encourages further research given the importance of wheat as a global commodity, the challenge of cloud cover in optical monitoring and the potential of direct estimation of the weight of heads where wheat production lies.     

Synergy of in situ and space borne observation for snow depth mapping in the Swiss Alps

The Swiss Federal Institute for Snow and Avalanche Research in Davos (SLF) provides snow depth maps for Switzerland on a spatial resolution of 1 km × 1 km. These snow depth maps are derived from snow station measurements using a spatial interpolation method based on the dependency of snow depth and altitude. During a winter season the number of operating snow stations varies and the area-wide snow depth interpolation becomes increasingly difficult in spring. The objective of the study is to develop an operational and near-real time method to calculate snow depth maps using a combination of in situ snow depth measurements and the snow cover extent provided from space borne observations. The operational daily snow cover product obtained from the polar-orbiting NOAA-AVHRR satellite is used to gain an additional set of virtual snow stations to densify the in situ measurements for an improved spatial interpolation. The capacity of this method is demonstrated on selected days during winter 2005. Cross-validation tests are conducted to examine the quantitative accuracy of the synergetic interpolation method. The error estimators prove the decrease in error variance and increase of overall accuracy pointing out the high capacity of this new interpolation method that can be run in near real-time over a large horizontal domain at high horizontal resolution. A solid method for snow–no snow classification in the processing of the satellite data is essential to the quality of the snow depth maps.     

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.     

Comparative assessment of thematic accuracy of GLC maps for specific applications using existing reference data

Inputs to various applications and models, current global land cover (GLC) maps are based on different data sources and methods. Therefore, comparing GLC maps is challenging. Statistical comparison of GLC maps is further complicated by the lack of a reference dataset that is suitable for validating multiple maps. This study utilizes the existing Globcover-2005 reference dataset to compare thematic accuracies of three GLC maps for the year 2005 (Globcover, LC-CCI and MODIS). We translated and reinterpreted the LCCS (land cover classification system) classifier information of the reference dataset into the different map legends. The three maps were evaluated for a variety of applications, i.e., general circulation models, dynamic global vegetation models, agriculture assessments, carbon estimation and biodiversity assessments, using weighted accuracy assessment. Based on the impact of land cover confusions on the overall weighted accuracy of the GLC maps, we identified map improvement priorities. Overall accuracies were 70.8 ± 1.4%, 71.4 ± 1.3%, and 61.3 ± 1.5% for LC-CCI, MODIS, and Globcover, respectively. Weighted accuracy assessments produced increased overall accuracies (80–93%) since not all class confusion errors are important for specific applications. As a common denominator for all applications, the classes mixed trees, shrubs, grasses, and cropland were identified as improvement priorities. The results demonstrate the necessity of accounting for dissimilarities in the importance of map classification errors for different user application. To determine the fitness of use of GLC maps, accuracy of GLC maps should be assessed per application; there is no single-figure accuracy estimate expressing map fitness for all purposes.     

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.     

Mapping seagrass coverage and spatial patterns with high spatial resolution IKONOS imagery

Seagrass habitats in subtidal coastal waters provide a variety of ecosystem functions and services and there is an increasing need to acquire information on spatial and temporal dynamics of this resource. Here, we explored the capability of IKONOS (IKO) data of high resolution (4 m) for mapping seagrass cover [submerged aquatic vegetation (%SAV) cover] along the mid-western coast of Florida, USA. We also compared seagrass maps produced with IKO data with that obtained using the Landsat TM sensor with lower resolution (30 m). Both IKO and TM data, collected in October 2009, were preprocessed to calculate water depth invariant bands to normalize the effect of varying depth on bottom spectra recorded by the two satellite sensors and further the textural information was extracted from IKO data. Our results demonstrate that the high resolution IKO sensor produced a higher accuracy than the TM sensor in a three-class % SAV cover classification. Of note is that the OA of %SAV cover mapping at our study area created with IKO data was 5–20% higher than that from other studies published. We also examined the spatial distribution of seagrass over a spatial range of 4–240 m using the Ripley’s K function [L(d)] and IKO data that represented four different grain sizes [4 m (one IKO pixel), 8 m (2 × 2 IKO pixels), 12 m (3 × 3 IKO pixels), and 16 m (4 × 4 IKO pixels)] from moderate-dense seagrass cover along a set of six transects. The Ripley’s K metric repeatedly indicated that seagrass cover representing 4 m × 4 m pixels displayed a dispersed (or slightly dispersed) pattern over distances of <4–8 m, and a random or slightly clustered pattern of cover over 9–240 m. The spatial pattern of seagrass cover created with the three additional grain sizes (i.e., 2 × 24 m IKO pixels, 3 × 34 m IKO pixels, and 4 × 4 m IKO pixels) show a dispersed (or slightly dispersed) pattern across 4–32 m and a random or slightly clustered pattern across 33–240 m. Given the first report on using satellite observations to quantify seagrass spatial patterns at a spatial scale from 4 m to 240 m, our novel analyses of moderate-dense SAV cover utilizing Ripley’s K function illustrate how data obtained from the IKO sensor revealed seagrass spatial information that would be undetected by the TM sensor with a 30 m pixel size. Use of the seagrass classification scheme here, along with data from the IKO sensor with enhanced resolution, offers an opportunity to synoptically record seagrass cover dynamics at both small and large spatial scales.     

Defining optimal DEM resolutions and point densities for modelling hydrologically sensitive areas in agricultural catchments dominated by microtopography

Defining critical source areas (CSAs) of diffuse pollution in agricultural catchments depends upon the accurate delineation of hydrologically sensitive areas (HSAs) at highest risk of generating surface runoff pathways. In topographically complex landscapes, this delineation is constrained by digital elevation model (DEM) resolution and the influence of microtopographic features. To address this, optimal DEM resolutions and point densities for spatially modelling HSAs were investigated, for onward use in delineating CSAs. The surface runoff framework was modelled using the Topographic Wetness Index (TWI) and maps were derived from 0.25 m LiDAR DEMs (40 bare-earth points m⁻²), resampled 1 m and 2 m LiDAR DEMs, and a radar generated 5 m DEM. Furthermore, the resampled 1 m and 2 m LiDAR DEMs were regenerated with reduced bare-earth point densities (5, 2, 1, 0.5, 0.25 and 0.125 points m⁻²) to analyse effects on elevation accuracy and important microtopographic features. Results were compared to surface runoff field observations in two 10 km² agricultural catchments for evaluation. Analysis showed that the accuracy of modelled HSAs using different thresholds (5%, 10% and 15% of the catchment area with the highest TWI values) was much higher using LiDAR data compared to the 5 m DEM (70–100% and 10–84%, respectively). This was attributed to the DEM capturing microtopographic features such as hedgerow banks, roads, tramlines and open agricultural drains, which acted as topographic barriers or channels that diverted runoff away from the hillslope scale flow direction. Furthermore, the identification of ‘breakthrough’ and ‘delivery’ points along runoff pathways where runoff and mobilised pollutants could be potentially transported between fields or delivered to the drainage channel network was much higher using LiDAR data compared to the 5 m DEM (75–100% and 0–100%, respectively). Optimal DEM resolutions of 1–2 m were identified for modelling HSAs, which balanced the need for microtopographic detail as well as surface generalisations required to model the natural hillslope scale movement of flow. Little loss of vertical accuracy was observed in 1–2 m LiDAR DEMs with reduced bare-earth point densities of 2–5 points m⁻², even at hedgerows. Further improvements in HSA models could be achieved if soil hydrological properties and the effects of flow sinks (filtered out in TWI models) on hydrological connectivity are also considered.     

Fire danger assessment in Iran based on geospatial information

Fire danger assessment is a vital issue to alleviate the impacts of wildland fires. In this study, a fire danger assessment system is proposed, which extensively uses geographical databases to characterize the spatial variations of fire danger conditions in Iran. This assessment requires three steps: (i) generation of the required input variables, (ii) methods to integrate those variables for creating synthetic indices and (iii) validation of those indices versus fire occurrence data. This fire danger model is based on previous works but adapted to Iranian conditions. It includes an estimation of the fire ignition potential (both considering human and climatic factors) and fire propagation potential. The former was generated from a logistic regression approach based on a wide range of input variables. The fire propagation probability was estimated from the Flammap fire behavior model. A first stage for validation of our fire danger system was based on comparing the estimated danger values to actual fire occurrence, based on satellite detected active fires and burned areas. The logistic regression model for fire ignition probability estimated 72.7% of true ignitions. Detected hotspots occurred more frequently in areas with higher fire ignition probability (average value: 0.65) than non hotspots (average value: 0.4). Propagation probability showed higher values for areas with higher proportion of burned area (r = 0.68, p < 0.001).     

Testing estimation of water surface in Italian rice district from MODIS satellite data

Recent changes in rice crop management within Northern Italy rice district led to a reduction of seeding in flooding condition, which may have an impact on reservoir water management and on the animal and plant communities that depend on the flooded paddies. Therefore, monitoring and quantifying the spatial and temporal variability of water presence in paddy fields is becoming important. In this study we present a method to estimate dynamics of presence of standing water (i.e. fraction of flooded area) in rice fields using MODIS data. First, we produced high resolution water presence maps from Landsat by thresholding the Normalised Difference Flood Index (NDFI) made: we made it by comparing five Landsat 8 images with field-obtained information about rice field status and water presence. Using these data we developed an empirical model to estimate the flooding fraction of each MODIS cell. Finally we validated the MODIS-based flooding maps with both Landsat and ground information. Results showed a good predictability of water surface from Landsat (OA = 92%) and a robust usability of MODIS data to predict water fraction (R² = 0.73, EF = 0.57, RMSE = 0.13 at 1 × 1 km resolution). Analysis showed that the predictive ability of the model decreases with the greening up of rice, so we used NDVI to automatically discriminate estimations for inaccurate cells in order to provide the water maps with a reliability flag. Results demonstrate that it is possible to monitor water dynamics in rice paddies using moderate resolution multispectral satellite data. The achievement is a proof of concept for the analysis of MODIS archives to investigate irrigation dynamics in the last 15 years to retrieve information for ecological and hydrological studies.     

Temporal analysis of remotely sensed turbidity in a coastal archipelago

A topographically fragmental archipelago with dynamic waters set the preconditions for assessing coherent remotely sensed information. We generated a turbidity dataset for an archipelago coast in the Baltic Sea from MERIS data (FSG L1b), using CoastColour L1P, L2R and L2W processors. We excluded land and mixed pixels by masking the imagery with accurate (1:10 000) shoreline data. Using temporal linear averaging (TLA), we produced satellite-imagery datasets applicable to temporal composites for the summer seasons of three years. The turbidity assessments and temporally averaged data were compared to in situ observations obtained with coastal monitoring programs. The ability of TLA to estimate missing pixel values was further assessed by cross-validation with the leave-one-out method. The correspondence between L2W turbidity and in situ observations was good (r = 0.89), and even after applying TLA the correspondence remained acceptable (r = 0.78). The datasets revealed spatially divergent temporal water characteristics, which may be relevant to the management, design of monitoring and habitat models. Monitoring observations may be spatially biased if the temporal succession of water properties is not taken into account in coastal areas with anisotropic dispersion of waters and asynchronous annual cycles. Accordingly, areas of varying turbidity may offer a different habitat for aquatic biota than areas of static turbidity, even though they may appear similar if water properties are measured for short annual periods.     

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.     

Retrieving aboveground biomass of wetland Phragmites australis (common reed) using a combination of airborne discrete-return LiDAR and hyperspectral data

Wetland biomass is essential for monitoring the stability and productivity of wetland ecosystems. Conventional field methods to measure or estimate wetland biomass are accurate and reliable, but expensive, time consuming and labor intensive. This research explored the potential for estimating wetland reed biomass using a combination of airborne discrete-return Light Detection and Ranging (LiDAR) and hyperspectral data. To derive the optimal predictor variables of reed biomass, a range of LiDAR and hyperspectral metrics at different spatial scales were regressed against the field-observed biomasses. The results showed that the LiDAR-derived H_p99 (99th percentile of the LiDAR height) and hyperspectral-calculated modified soil-adjusted vegetation index (MSAVI) were the best metrics for estimating reed biomass using the single regression model. Although the LiDAR data yielded a higher estimation accuracy compared to the hyperspectral data, the combination of LiDAR and hyperspectral data produced a more accurate prediction model for reed biomass (R² = 0.648, RMSE = 167.546 g/m², RMSE_r = 20.71%) than LiDAR data alone. Thus, combining LiDAR data with hyperspectral data has a great potential for improving the accuracy of aboveground biomass estimation.     

Determination of the factors governing soil erodibility using hyperspectral visible and near-infrared reflectance spectroscopy

Soil erodibility, which is difficult to estimate and upscaling, was determined in this study using multiple spectral models of soil properties (soil organic matter (SOM), water-stable aggregates (WSA) > 0.25 mm, the geometric mean radius (D_g)). Herein, the soil erodibility indicators were calculated, and soil properties were quantitatively analyzed based on laboratory simulation experiments involving two selected contrasting soils. In addition, continuous wavelet transformation was applied to the reflectance spectra (350–2500 nm) of 65 soil samples from the study area. To build the relationship, the soil properties that control erodibility were identified prior to the spectral analysis. In this study, the SOM, D_g and WSA >0.25 mm were selected to represent the most significant soil properties controlling erodibility and describe the erodibility indicator based on a logarithmic regression model as a function of SOM or WSA > 0.25 mm. Five, six and three wavelet features were observed to calibrate the estimated soil properties model, and the best performance was obtained with a combination feature regression model for SOM (R² = 0.86, p < 0.01), D_g (R² = 0.79, p < 0.01) and WSA >0.25 mm (R² = 0.61, p < 0.01), respectively. One part of the wavelet features captured amplitude variations in the broad shape of the reflectance spectra, and another part captured variations in the shape and depth of the soil dry substances. The wavelet features for the validated dataset used to predict the SOM, WSA >0.25 mm and D_g were not significantly different compared with the calibrated dataset. The synthesized spectral models of soil properties, and the formation of a new equation for soil erodibility transformed from the spectral models of soil properties are presented in this study. These results show that a spectral analytical approach can be applied to complex datasets and provide new insights into emerging dynamic variation with erodibility estimation.     

Estimate of heavy metals in soil and streams using combined geochemistry and field spectroscopy in Wan-sheng mining area,Chongqing, China

Heavy metals contaminated soils and water will become a major environmental issue in the mining areas. This paper intends to use field hyper-spectra to estimate the heavy metals in the soil and water in Wan-sheng mining area in Chongqing. With analyzing the spectra of soil and water, the spectral features deriving from the spectral of the soils and water can be found to build the models between these features and the contents of Al, Cu and Cr in the soil and water by using the Stepwise Multiple Linear Regression (SMLR). The spectral features of Al are: 480 nm, 500 nm, 565 nm, 610 nm, 680 nm, 750 nm, 1000 nm, 1430 nm, 1755 nm, 1887 nm, 1920 nm, 1950 nm, 2210 nm, 2260 nm; The spectral features of Cu are: 480 nm, 500 nm, 610 nm, 750 nm, 860 nm, 1300 nm, 1430 nm, 1920 nm, 2150 nm, 2260 nm; And the spectral features of Cr are: 480 nm, 500 nm, 610 nm, 715 nm, 750 nm, 860 nm, 1300 nm, 1430 nm, 1755 nm, 1920 nm, 1950 nm. With these features, the best models to estimate the heavy metals in the study area were built according to the maximal R². The R² of the models of estimating Al, Cu and Cr in the soil and water are 0.813, 0.638, 0.604 and 0.742, 0.584, 0.513 respectively. And the gradient maps of these three types of heavy metals’ concentrations can be created by using the Inverse distance weighted (IDW).The gradient maps indicate that the heavy metals in the soil have similar patterns, but in the North-west of the streams in the study area, the contents are of great differences. These results show that it is feasible to predict contaminated heavy metals in the soils and streams due to mining activities by using the rapid and cost-effective field spectroscopy.