Improving the spectral unmixing algorithm to map water turbidity Distributions

https://doi.org/10.1016/j.envsoft.2009.02.013Get rights and content

Abstract

In this paper we evaluate the suitability of the spectral unmixing algorithm to map the turbidity in the Curuai floodplain lake and enhance its applicability using autocorrelation modelling. The Spectral Unmixing Model (SMM) was applied to a Moderate Resolution Imaging Spectroradiometer (MODIS) surface reflectance (MOD09) image, taking in-situ measurements close to the acquisition date. Fraction images of inorganic matter-laden water, dissolved organic matter-laden water, and phytoplankton-laden water were generated by SMM, using 4 MODIS spectral bands (blue, green, red, and near infrared). These endmembers were selected based on the dominance of these components, which affect water turbidity. These fraction images allowed assessing the turbidity distribution in the study area but showing only places with high or low turbidity. The kernel estimation algorithm was then used to verify the spatial correlation among the in-situ measurement data. The occurrence of clusters suggests that there are different spatial water regimes. One spatial regression model was then compiled for each water regime, each of which presented a better turbidity estimation as opposed to the one derived from the Ordinary Least Square (OLS). The methodology applied was hence useful to analyze the spatial distribution of turbidity in the Curuai floodplain lake.

Introduction

Water quality modelling is the linkage between the sources of pollution and the in-stream water quality of a given water body (Freni et al., 2009). Different models have already been applied to a variety of studies including the assessment of relocation impact (Dias et al., 2009), trans-boundary conflicts (Zhao, 2009). Many of those models, however, use in-situ water quality data as input (Wu et al., 2009).

Conventional water quality monitoring is expensive and time consuming. This is particularly problematic if the water bodies to be examined are large. Conventional techniques also bring about a high probability of undersampling (Hadjimitsis et al., 2006). Conversely, remote sensing is a powerful tool to assess aquatic systems and is particularly useful in remote areas such as the Amazon lakes (Alcântara et al., 2008).

The remote sensing technique has been used not only for water surface modelling but also for modelling land surface temperature (Sheng et al., 2009), annual variability of nitrate concentrations (Montzka et al., 2008), spatially distributed sediment budgets (Wilkinson et al., 2009), crop yield (Pan et al., 2009).

Data collected using this technique can provide a synoptic overview of such large aquatic environments, which could otherwise not be observed at a glance (Dekker et al., 1995). However, remote sensing is not easily applied to aquatic environment monitoring mainly because of concerns related to the mixture of the optically active substances (OAS) in the water. Several approaches have been proposed to cope with this issue such as derivative analysis (Goodin et al., 1993), the continuum removal (Kruse et al., 1993), and spectral mixture analysis (Novo and Shimabukuro, 1994, Oyama et al., 2009).

The two first approaches are more suitable for hyperspectral images, whereas the spectral mixture analysis can be used for both hyperspectral and multispectral images. The SMM has largely been used for spectral mixture analysis, uncoupling the reflectance of each image pixel (Tyler et al., 2006) into the proportion of each water component contributing to the signal. The result of a spectral mixture analysis is a set of fraction images representing the proportion of each water component per image pixel. This technique has been applied to TM/Landsat images to determine the concentration of suspended particles (Mertes et al., 1993), chlorophyll-a concentration (Novo and Shimabukuro, 1994); as well as to MODIS images, to determine the chlorophyll-a concentration in the Amazon floodplain (Novo et al., 2006), to characterize the composition of optically complex waters in the Amazon (Rudorff et al., 2007), and to study turbidity distribution in the Amazon floodplain (Alcântara et al., 2008).

Remote sensing data has been extensively used to detect and to quantify water quality variables in lakes and reservoirs (Kloiber et al., 2002). One of the most important variables to monitor water quality is turbidity, because it reveals information about the availability of light, as determined by the concentration of inorganic and organic (alive and dead) components in the water column (Goodin et al., 1993). Although turbidity is brought about by organic and inorganic particles, one unresolved issue is to distinguish between them using remote sensing (Wetzel, 2001). The Spectral Unmixing Model (SMM) can, however, be useful to analyze the turbidity caused by inorganic particles and by phytoplankton cell scattering.

A previous study used the SMM to decompose the optical water components and map the turbidity distribution in the Curuai floodplain at high water level (Alcântara et al., 2008). However, it is uncertain if this approach will work when applied to rising water conditions.

We believe that changes in turbidity can be detected in remote sensing images because they bring about changes in the upward irradiance leaving the water surface. This paper describes the experiment carried out to test this hypothesis and to assess the Amazon floodplain lake turbidity using MODIS imagery.

Section snippets

Study area

The Amazon River basin drains an area of approximately 6 × 106 km2, which represents 5% of the Earth surface. The Central Amazon has large floodplains covering around 300,900 km2 (Hess et al., 2003), including 110,000 km2 of the main part ‘Varzeas’ (floodplains of the white water river; Junk, 1997). At high water, the Amazon River flows into the floodplains, and fills both temporary and permanent lakes which might merge to each other. The ‘Lago Curuai’ floodplain (Fig. 1) covers an area varying

In-situ measurements and remote sensing data

The turbidity ground data acquisition was carried out from February 1st to February 14th 2004, during the rising water period. Turbidity measurements were taken at 215 sampling stations using the HORIBA U-10 multi-sensor. This equipment provides turbidity measurements in NTU (Nephelometric Turbidity Unit) with a resolution of 1 NTU. The locations of the sampling stations were determined with the aid of spectral analyses of Landsat/TM images taken at similar water level (Barbosa, 2005). These

Methodological approach

The turbidity distribution was assessed using fraction images derived from the Linear Spectral Mixing Model, using four MODIS spectral bands (3 – blue, 4 – green, 1 – red and 2 – near infrared) with a spatial resolution of 250 m. In order to evaluate the turbidity distribution observed in the MODIS fraction images, in-situ measurements acquired during in February 2004 (a few days apart of the MODIS acquisition) were used to apply the Ordinary Least Square (OLS), spatial lag, and spatial error

Results and discussion

Fig. 7a shows that the water, of the entire Curuai floodplain lake was rich in inorganic matter (Iss), with a particularly high proportion in the Poção lake (Fig. 7b). The images illustrating the distribution of dissolved organic matter (Dom) revealed that this was particularly apparent in the Salé lake, and in the border region (i.e. the region between water and forest). This is mainly due to the fact that some organic matter in decomposition is transported into the floodplain by the water

Conclusions

This present work evaluates the suitability of the Spectral Unmixing Model to map the turbidity distribution in the Curuai floodplain. The main conclusions are:

  • 1.

    The fraction images for the endmembers selected directly from the MODIS image based on dominance of water components allowed assessing the turbidity in the Curuai floodplain lake.

  • 2.

    Owing to non-linearity in the Amazon floodplain waters, the unmixing model does not work in an optimal way. This is also due to autocorrelation presented in the

Acknowledgements

This work was supported by the Brazilian funding agency FAPESP under Grant 02/09911-1. The work of E.H. de Alcântara was supported by the Brazilian Council for Scientific and Technological Development (CNPq). The authors would like to thank Carsten Reinhard from the School of Biological Sciences, University of Bristol, for the English review and the anonymous reviews for their constructive comments and suggestions, which helped improve this paper.

References (41)

  • G. Pan et al.

    Using QuickBird imagery and a production efficiency model to improve crop yield estimation in the semi-arid hilly Loess Plateau, China

    Environ. Model. Softw.

    (2009)
  • J. Sheng et al.

    Comparison of land surface temperature (LST) modeled with a spatially-distributed solar radiation model (SRAD) and remote sensing data

    Environ. Model. Softw.

    (2009)
  • S.N. Wilkinson et al.

    Modelling and testing spatially distributed sediment budgets to relate erosion processes to sediment yields

    Environ. Model. Softw.

    (2009)
  • J. Wu et al.

    Application of chaos and fractal models to water quality time series prediction

    Environ. Model. Softw.

    (2009)
  • L. Zhao

    Model of collective cooperation and reallocation of benefits related to conflicts over water pollution across regional boundaries in a Chinese river basin

    Environ. Model. Softw.

    (2009)
  • E.H. Alcântara et al.

    Turbidity in the Amazon floodplain assessed through a spatial regression model applied to fraction images derived from MODIS/Terra

    IEEE Trans. Geosci. Remote Sens.

    (2008)
  • L. Anselin

    Spatial Econometrics: Methods and Models

    (1988)
  • T.C. Bailey et al.

    Interactive Spatial Data Analysis

    (1995)
  • B.H. Baltaqi

    A Companion to Theoretical Econometrics

    (2003)
  • Barbosa, C.C.F., 2005. Sensoriamento remoto da dinâmica de circulação da água do sistema planície de Curuai/ Rio...
  • Cited by (0)

    View full text