Next Article in Journal
Chemical Characterization of Dew Water Collected in Different Geographic Regions of Poland
Previous Article in Journal
Spatio-Temporal Analysis of Forest Fire Risk and Danger Using LANDSAT Imagery
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Detecting Aquatic Vegetation Changes in Taihu Lake, China Using Multi-temporal Satellite Imagery

State Key Laboratory of Lake Science and Environment, Nanjing Institute of Geography and Limnology, Chinese Academy of Sciences, Nanjing, PO BOX 210008, China
*
Author to whom correspondence should be addressed.
Sensors 2008, 8(6), 3988-4005; https://doi.org/10.3390/s8063988
Submission received: 6 May 2008 / Accepted: 19 June 2008 / Published: 25 June 2008

Abstract

:
We have measured the water quality and bio-optical parameters of 94 samples from Taihu Lake in situ and/or in the lab between June 10-18, 2007. A transparency-assisted decision tree was developed to more accurately divide the aquatic vegetation zone into a floating vegetation-dominated zone and a submerged vegetation-dominated zone, whose respective present biomass retrieval models were easily developed with an empirical approach because of the quasi-concurrence of ground field investigations with the satellite sensor flight over the lake. The significant quantitative relationships between the vegetation index NDVI (Normalized Difference Vegetation Index) of different images at different times were used to help develop the past biomass retrieval model on the basis of the present developed model. In Taihu Lake, the total covering area of aquatic vegetations decreased from 454.6 km2 in 2001 to 364.1 km2 in 2007. Correspondingly, the total biomass decreased from 489,000 tons in 2001 to 406,000 tons in 2007, suggesting that a great change in the ecological environment has been taking place in Taihu Lake over this period.

1. Introduction

Aquatic vegetation, generally existing in the shallow near-shore area, is a key component of lake ecosystems. This vegetation provides food, shelter and breeding habitats for aquatic animals like invertebrates, fish and wading birds, and helps maintain the balance of the lake ecosystem. In addition, it also plays an important role in maintaining a clean lake water quality by stabilizing sediments and providing a substrate for periphyton that actively removes nitrogen and phosphorus from the water column. At times and locations where submerged vegetation is very abundant, water is clear, and phytoplankton blooms are rare. It almost becomes a token indicator to determine whether the water quality can be expected to be good or not. The spatial extent and biomass of aquatic vegetation, especially submerged aquatic vegetation, are expected to increase significantly as a result of regional restoration programs, so the past and present status of a regional biomass and distribution of aquatic vegetation are necessary for the monitoring, management and understanding of shallow aquatic ecosystems, which need baseline measurements to document their extent and condition. Determining the vegetation cover and biomass is difficult, however, both on small and large spatial scales because of the spatial heterogeneity of these communities [1]. It has been indicated that aquatic vegetation yields spectrally distinct signals governed by the density of the vegetation, the openness of the canopy and the amounts, forms and orientations of the leaves [2-5]. Since traditional quantitative ground investigations on the scale of a whole lake are laborious, remote sensing methods are increasingly being used for mapping aquatic vegetation and estimating their distribution and biomass [6-9].
Multi-temporal remote sensing data can give valuable information about changes that have taken place in a given area [8]. It is also well suited to identify both emergent and submerged vegetation [1, 10]. Aerial photographs were commonly used for mapping aquatic vegetation in the earlier period [3, 5, 7], but with rapid development of satellite sensors, satellite multi-spectral scanner data are widely used at present. Landsat TM (Thematic Mapper) has been proven very effective for aquatic vegetation distribution and biomass mapping applications [11, 12], however, more detailed spatial monitoring is still not possible due to its somewhat coarse spatial resolution [13]. IKONOS image, with enough small pixel size, may be used to make fine-level habitat discrimination by making full use of within-habitat textural information in a supervised classification, which significantly improved the thematic map accuracy compared with Landsat TM image; however, it cannot spectrally resolve changes in community structure [14]. The use of Hymap images with hyper-spectral resolution has been attempted to provide a spectral discrimination of submerged macrophytes [15, 16]. Supervised classification [11, 14] or visual interpretation [1] is a common method for discrimination of emerged and submerged vegetation, both of which have significantly different spectral features in the visible and near-infrared wavelengths. Ackleson and Klemas [17] demonstrated that the accuracy of results could be increased significantly when a prior knowledge of water depth was added to the classification crteria because of the close correlation between the distribution of aquatic vegetation and water depth [18]. The spectral confusion between submerged vegetation canopy density and deep water could be solved significantly. Malthus and George [19] demonstrated that a combination of band 3 (520-600 nm), band 7 (760-900 nm), and band 8 (910-1050 nm) data from the Daedalus Airborne Thematic Mapper could discriminate between different macrophyte growth forms. The DN (Digital Number) value, a sum of contributions due to atmosphere, water column and bottom, has been most commonly used to estimate vegetation biomass with an empirical linear or non-linear fitting model [9]. If the influence of the water column could be removed from the remotely sensed images, the potential of optical remote sensing in littoral applications would be extended and the accuracy of classification and biomass retrieval would be improved to a certain extent [20, 21]. As optical field equipments have become available, the bio-optical model has appeared and is being widely used, which make it possible to accurately remove the contribution due to water column [22, 23] and then to separate between the signatures coming from water column and bottom.
This paper aims to present an approach to estimate the past biomass of aquatic vegetation in a shallow inland lake using the historical satellite imageries without any corresponding field investigations, based on the present satellite imagery provided along with concurrent or quasi-concurrent ground field investigations, and then to determine the historical change of the aquatic vegetation in Taihu Lake. We first present some general information about Taihu Lake and the field investigation, including instruments and sampling methods, and satellite imageries together with its preprocessing method. Then we describe the method for developing the quantitative estimation model for aquatic vegetation in Taihu Lake. Finally, some conclusions are drawn on the basis of the full discussion.

2. Study area

Taihu Lake, a large shallow lake with an average depth of 1.9 m (max. 2.6 m), and covering an area of 2427.8 km2 (including 51 islands), is one of the five largest fresh-water lakes of China and is located in the core of the Yangtze Delta in the lower reaches of Yangtze River, in East China, an area that has a developed economy (Figure 1). Unfortunately, in recent years, water quality has deteriorated and cyanobacteria blooms have covered large areas of the lake every summer since the 1990s. There exist four types of aquatic vegetation: (a) emerged vegetation, including Phragnites communis, Typha angustifolia, Zizania latifolia and so on, (b) submerged vegetation, including Ceratophyllum demersum, Vallisneria spiralis, P. malaianus and so on, (c) floating vegetation, including Eichhornia crassipes, Lemna minor and so on, and (d) floating-leaved vegetation, including Nyphar pumilum, Nymphaea tetragona and so on. They are mainly distributed in the several bays of the east and southeast and in the south littoral zones, where the water is clear and cyanobacteria blooms are rare.

3. Material

3.1 Field data

The campaign was carried out along the preset transects for 94 samples on 10-18 June 2007 (Figure 2), which were positioned by differential GPS (Global Positioning System, Trimble Ltd., USA) with a positioning precision of 2-5 m. Before the field investigations, the sampling locations were preset according to the prior knowledge of the aquatic vegetation distribution with help of a 1:50,000 scale digital topographic map. This prior knowledge comes from the previous six field investigations of the spatial distribution of aquatic vegetation for 95 samples in October 2006, 133 samples in May 2006, 27 samples in October 2005, 80 samples in June 2005, 25 samples in July 2003, and 78 samples in September 2002, respectively.
The remote sensing reflectance and the backscattering coefficient, respectively, were measured in situ with a dual channel spectrometer FieldSpec 931 (ASD Ltd.) and a HydroScat-6 Spectral Backscattering Sensor (HS-6, HOBI Lab Inc.) mounted at six wavelengths (centered at 442, 488, 532, 589, 676 and 852 nm, respectively). The instruments, methods of measurement and data processing are the same as in Ma et al. [22, 23]. Water samples were collected from the surface to about 30 cm below in the vertical direction with a standard 21 polyethylene water-fetching instrument immediately after measuring the spectra. They were then held in a freezer half filled with ice bags for preservation for approximately 4 h every afternoon, and then returned to the laboratory for concentration and absorption measurements. The transparency was measured with a Secchi disk at the shade side of the boat. We employed the boat-based rake method as described by Rodusky et al. [19], with an oyster-tongs-like rake apparatus operated from a boat, to collect aquatic vegetation, which were put on the deck or into a plastic basin for flushing with local lake waters in order to remove the attached mud. We also measured the reflectance of aquatic vegetation, mud and the mix of vegetation and mud collected from the bottom using the same spectrometer in the wavelength range of 350-1050 nm. All the clean vegetation at one site, including leaves brushed off, was put into a plastic bag with a serial number. The plastic bag was then held in a freezer half filled with ice bags for preservation and together with water samples, returned to the laboratory for species classification and biomass measurements. In addition, relative species depth (surface, sub-surface, and deep) within the water column was measured at each sample site, and the sampling time was also recorded point by point.
In the lab, the concentrations of DOC (Dissolved Organic Color, CDOC), chlorophyll-a (CCHL), and SPM (Suspended Particulate Matter, CSPM) were measured according to the investigation criteria for Chinese lakes [24]. We also measured the absorptions of phytoplankton pigment (aph), non-algal particulate (ad), and CDOM (Colored Dissolved Organic Matter, ag) according to the NASA SeaWiFS protocols described by Mueller et al. [25]. Details about the concentration and absorption measurements can be found in Ma et al. [22, 23]. The absorption and backscattering due to pure water, aw and bbw, refers to the definitions of Pope and Fry [26].

3.2 Satellite data

IRS P6 (Indian Remote Sensing Satellite Resourcesat-1) LISS-3 (Linear Imaging Self-Scanning Sensor-3) image with a spatial resolution of 23.5 m, mounted at three bands in the range from 520 nm to 860 nm and one band from 1550 nm to 1700 nm was first used to develop the basic and primary model for biomass estimation of aquatic vegetation in Taihu Lake. The imaging date selected for this research was 15 June 2007. In Taihu Lake, based on our experience, aquatic vegetation grows slowly before July every year, so we may reasonably assume that the vegetation changes little in the short time range of our field investigation and we may consider that the ground field investigation and the LISS-3 sensor flight over the lake are almost concurrent. In order to establish the annual change of aquatic vegetation in Taihu Lake, we also employed some Landsat TM images, with a spatial resolution of 30 m, with hardly any cloud cover and almost the same date, respectively, acquired on 26 July 2001, 13 July 2002, and 26 July 2004.

4. Methods

4.1 Image preprocessing

Two methods were employed to remove the contributions due to the atmosphere in order to develop the two different biomass estimation models, viz. an empirical model and an analytical model. One is referred to an assumption that the water-leaving radiance in the near-infrared wavelength is zero, which is commonly used in the case one water color remote sensing. This method is simply approximate and aims to empirically estimate the aquatic vegetation biomass. A pixel with the lowest DN value exists in each of bands, which is assumed to be fully due to atmosphere. The pixel is first identified if it meets the following two requirements: (a) it is in optically deep clear water, or in the shallow clear water but the bottom is mud without concomitant vegetation, and (b) its DN value is the lowest in the corresponding band mounted on the LISS-3 or TM sensor. Then we may subtract the contribution due to atmosphere from total signals to leave only the signal due to water column and bottom. The steps mentioned above are repeated until all the bands mounted at the sensor are finished. The other is complexly approximate, referring to the methods for terrestrial remote sensing. Here, FLAASH (Fast Line-of-sight Atmospheric Analysis of Spectral Hypercubes), an atmospheric correction code based on the MODTRAN (MODerate resolution atmospheric TRANsmission) 4 radiative transfer model, is used to convert LISS-3 and TM sensor radiance data to apparent reflectance, aiming to develop an analytical approach to estimate aquatic biomass. Subsequently, a geometrical correction was done using a 1:50,000 scale topographic map (spheroid: Krasovsky 1940, projection: Gauss-Krüger). The RMSE (Root Mean Square Error) for each of the images was maintained at less than 0.5 pixels. Finally, the Taihu Lake boundary from the same topographic map is used as a binary mask to create a water body-only image.

4.2 Image classification

The boundary between various classes of emergent, floating, floating-leaved and submerged vegetation, whose biomass in unit area have significant differences, is not distinct, as these communities merge into one another. To make it suitable for the remote sensing application and then to make the biomass estimation precision as high as possible, the aquatic vegetation zone is divided into two types: (a) Type one, with a characteristic of the leaf standing above the water surface, including the floating-leaved and emergent vegetation, named a floating vegetation-dominated zone, and (b) Type two, with the characteristic feature of the whole plant from canopy to root being submerged beneath the water surface, including submerged and floating vegetation, named as the submerged vegetation-dominated zone. Spectral features due to Type one in the range of red and near-infrared wavelengths are greatly different from the ones due to Type two (Figure 3). The two types of vegetations are easily separable in optical satellite remote sensing images.
Prior to separating the aquatic vegetation into the two types, it is of vital importance to first separate them from lake water. The submerged vegetation-dominated hydrophytes, however, commonly show a similar spectral feature with lake waters containing many phytoplankton or suspended particulate matter (Figure 3).
It is not easy to separate them only depending on their spectral characteristics. So we add the prior knowledge, including water depth and transparency, into the classification process to develop a decision tree as the following steps: (a) separate the vegetation and quasi-vegetation from lake water with help of the prior knowledge about relationships between depth, transparency and vegetation distribution, (b) to calculate vegetation index RVI (Ratio Vegetation Index, a ratio of reflectance in the near-infrared waveband to in the red wavelength) and NDVI (Normalized Difference Vegetation Index, a ratio of the reflectance difference between in the near-infrared and red wavebands to the sum of the reflectance in the two wavebands), and then to remove Type one according to the correct threshold values of RVI, NDVI and SDT (Secchi Disk Transparency), (c) to change the threshold values of RVI, NDVI and SDT to remove the non-vegetation from vegetation and quasi-vegetation zone, and then only leave the Type two.
The depth comes from a bottom topographic map with a scale of 1:50,000 mapped in 1998, and corrected in situ during this field investigation. The transparency may be retrieved from an empirical model given by Equation (1). There are good relationships, with determination coefficients (R2) of more than 0.68, between the Secchi disk transparencies measured in situ and the atmosphere-corrected DN values of the first (B1) and the second band (B2) mounted at LISS-3 sensor, respectively (Figure 4). However, there is a better relationship between the Secchi disk transparencies measured in situ and both of the atmosphere-corrected DN values of the former two bands mounted at LISS-3 sensor given as:
SDT = 279.330 0.708 × B 1 1.033 × B 2 ( R 2 = 0.82 , N = 70 )
where SDT (in cm), R2, B1 and B2 are the same as mentioned above; N is the total sampling number, not including samples in the floating vegetation-dominated zone.

4.3 Biomass estimation

The submerged vegetation zone is optically shallow in Taihu Lake [22], so [27]:
r r s = r r s c + r r s B r r s d p ( 1 exp { [ 1 cos ( θ w ) + D u c ] κ H } ) + 1 π ρ exp { [ 1 cos ( θ w ) + D u B ] κ H }
where
r r s d p ( 0.084 + 0.170 u ) u
D u c 1.03 ( 1 + 2.4 u ) 0.5
D u B 1.04 ( 1 + 5.4 u ) 0.5
u = b b / ( a + b b )
κ = a + b b
R r s 0.5 r r s 1 1.5 r r s
where Rrs is remote sensing reflectance; rrs is the remote sensing reflectance just below the surface; rrsc and rrsB are the signals, respectively, due to water column and bottom; θw is the subsurface solar zenith angle; H is the bottom depth; a is the total absorption coefficient, a sum of aph, ad, ag and aw; bb is the total backscattering coefficient measured by HS-6.
The FLAASH atmosphere-corrected signal was first converted to the signal just below the surface, which might be decomposed into the contributions due to water column and due to lake bottom using Equation (2) on the basis of the hypothesis conditions of Figure 5: (a) infinitely horizontal bottom with a homogenous mix of vegetation and mush with a Lambert property, and (b) homogenous water quality in vertical, so the contribution due to water column can be approximately removed by Equation (2) to leave the remote sensing reflectance only due to bottom. Then we developed the biomass estimation model by the relationship between the vegetation biomass and the remote sensing reflectance only due to bottom, however, which always was bad no matter what efforts we paid out. So here the semi-analytical approach is not suitable for aquatic vegetation biomass estimation in Taihu Lake. Subsequently, we have to attempt to develop an empirical model between the atmosphere-corrected DN value and the biomass measured in the field investigation.
In the LISS-3 image on 15 June 2007 good empirical relationships between the biomass in unit area of the two types of aquatic vegetation mentioned above, respectively, and NDVI and SDT, existed which can be described as:
{ B M 1 LISS 3 = 2.760 + 4.081 × NDV I LISS 3 ( R 2 = 0.623 , N = 14 ) B M 2 LISS 3 = 1.260 0.011 × SDT 2.415 × NDV I LISS 3 ( R 2 = 0.782 , N = 28 )
or
{ B M 1 LISS 3 = 2.760 + 4.081 × NDV I LISS 3 ( R 2 = 0.623 , N = 14 ) B M 2 LISS 3 = 0.137 3.542 × NDV I LISS 3 ( R 2 = 0.697 , N = 28 )
where BM1LISS3 (in kg/m2), BM2LISS3 (in kg/m2) and NDVILISS3, respectively, are the vegetation biomass in unit area of Type one and Type two and the vegetation index NDVI from LISS-3 image on 15 June 2007; SDT (in cm) is transparency coming from Equation (1); R2 is the determination coefficient and N is the total samples in the floating vegetation-dominated or submerged vegetation-dominated zone. There was no corresponding biomass data from field investigations on 26 July 2001, 13 July 2002, and 26 July 2004 when the Landsat TM sensor was flying over the Taihu Lake. In order to approximately acquire the historical aquatic vegetation biomass at the corresponding time, these Landsat TM image pixels are firstly re-sampled into the same spatial resolution of 23.5 m as LISS-3 image pixel. And then we use the pure pixels with growing vegetations in the land, whose canopy is almost invariable between TM and LISS-3 images, to develop the quantitative NDVI relationships between the two images, so:
NDV I LISS 3 = 0.114 + 0.224 × NDV I T M 2001 ( R 2 = 0.660 , N = 40 )
NDV I LISS 3 = 0.083 + 2.269 × NDV I T M 2002 ( R 2 = 0.561 , N = 40 )
NDV I LISS 3 = 0.104 + 0.310 × NDV I T M 2004 ( R 2 = 0.670 , N = 40 )
where NDVITM2001, NDVITM2002 and NDVITM2004 are the vegetation index NDVI from Landsat TM images, respectively, on 26 July 2001, 13 July 2002, 26 July 2004, R2 is the determination coefficient and N is the total sample number. It's a pity that the relationships of transparency between TM and LISS-3 images cannot be set up because of no transparency data measured in situ at the time when TM images was acquired. However, the water quality in those vegetation zones is hardly variational from 2001 to 2007 [28]. So we could, respectively, substitute Equation (11)-(13) into Equation (10), so:
{ B M 1 T M 2001 = 2.295 + 0.914 × NDV I T M 2001 B M 2 T M 2001 = 0.541 0.793 × NDV I T M 2001
{ B M 1 T M 2002 = 2.421 + 1.098 × NDV I T M 2002 B M 2 T M 2002 = 0.431 0.953 × NDV I T M 2002
{ B M 1 T M 2004 = 2.336 + 1.265 × NDV I T M 2004 B M 2 T M 2004 = 0.505 1.098 × NDV I T M 2004
where BM1 and BM2 (in kg/m2), respectively, denote the vegetation biomass of Type one and Type two; the subscripts represent the sensor and the date when the image was acquired by the sensor (here 26 July 2001, 13 July 2002, 26 July 2004). Finally, the error matrix from the powerful Accuracy Assessment tool in the ERDAS 9.0 software was used to evaluate the remote sensing classification accuracy and then the estimation accuracy for vegetation biomass in 2007 was assessed by comparison of the estimated and the measured ones, based on which, the accuracies in 2001, 2002 and 2004 were estimated by the error transfer law commonly used in the field of surveying and mapping sciences [29].

5. Results

5.1 Field investigation

This investigation showed that aquatic vegetations with a wet biomass range from 0.11 to 3.98 kg/m2 and a standard deviation of 0.98 kg/m2 existed only at 42 sampling sites. The emerged vegetation mainly is Phragmites communi, located within a near-shore range of about 100 m. It is outside this biomass estimation. The remainder is growing in six clusters: Potamogeton malaianus cluster, Nymphoides peltatum cluster, Nymphoides peltatum-Potamogeton malaianus cluster, Vallisneria natanus cluster, Vallisneria natanus-Elodea nuttalli cluster, and Najas marina cluster, as in Table 1. They can be classified into two classes: (a) the floating vegetation-dominated hydrophyte, growing at 14 sampling sites, mainly including Nymphoides peltatums, and (b) the submerged vegetation-dominated hydrophyte, growing at 28 sampling sites, mainly including Potamogeton malaianu.

5.1 Remote sensing classification and estimation

The satellite image classification, illustrated in Figure 6 and Table 2, shows in Taihu Lake that the aquatic vegetation has undergone a great spatial change since 2001. The coverage area in both 2001 and 2002 were basically constant around 450 km2, and in 2004 increased up to 482 km2.
However in 2007, it decreased rapidly to 361 km2. The vegetation in Gongshan Bay, Zhenhu Bay, Guangfu Bay and Xukou Bay were spreading around in 2001-2002; however in 2002-2004, they began to move back towards to the bay interior or one side of the shore. Noteworthily, the vegetation in Gongshan Bay almost completely vanished according to the remote sensing investigations in 2007. In fact, the field investigation showed that in the north of Gongshan Bay there is only a little vegetation with a very thinly scattered distribution. The same situation as in Gongshan Bay in 2007 also happened in Zhushan Bay in 2002. Some degradation happened to a certain extent along the route from Xiaoleishan island to the mouth of East Taihu Bay in 2002. On the contrary, since 2004, some vegetation was growing near the shore at the south of Xishan island. In addition, we obviously observe that some vegetation has appeared along the southern littoral zone from Qidou to Huzhou since 2004, and that it was spreading towards to the north at the same time of extending towards to the west along the shore.
The estimated vegetation biomass following the decision tree classification, illustrated as Figure 7 and Table 2, shows that the total biomass was basically stable about 480 thousand ton in 2001 and 2002, then increased up to 528 thousand ton in 2004, and then decreased rapidly up to 406 thousand ton in 2007 because of its covering area decreasing rapidly. However, in some special lake zone such as Zhenhu Bay and Guangfu Bay, the vegetation biomass increased significantly but its covering area only has a little difference, which shows that the vegetation has a more flourishing growth in these zones.
The classification precisions are up to 90.8%, 88.9%, 90.2% and 89.2%, respectively, for LISS-3 image in 2007, TM images in 2001, 2002 and 2004. For the LISS-3 image in 2007, the maximum and minimum absolute errors of biomass estimation, respectively, are 2.60 kg/m2 and 0.001 kg/m2, and the maximum and minimum relative errors, respectively, are 287.03% and 0.10%, which were acquired by analysis and comparison between 46 sampling sits growing aquatic vegetations (Figure 8). The average absolute errors of Equation (11)-(13) are, respectively, 0.049, 0.048, and 0.044, and the corresponding average relative errors are, respectively, 34.64%, 35.49%, and 29.55%. The total biomass estimation error in 2007 is about 46.3 thousand ton, and the errors in 2004, 2002 and 2001 are, respectively, about 70.2, 67.5, and 64.4 thousand ton according to the error transfer law from Equation (11)-(13), through Equation (10), to Equation (14)-(16).

6. Discussion

It is difficult or impossible to accurately estimate the historical biomass of aquatic vegetation and then to discover its change between the past and the present if there was no field investigations and also no corresponding remote sensing image obtained at the same time. The present empirical estimation model could be developed on the basis of the field investigation and the concurrent or quasi-concurrent satellite imagery at the corresponding time, and then may be used to work out the past estimation model by rightly fusing the relationship between vegetation indexes from the present and past satellite imageries. Another approach is firstly to remove the contribution due to water column from the atmosphere-corrected signals through the radiative transfer in natural waters on the basis of several perfect assumptions. And then the remote sensing signal only due to the bottom was acquired and then used to develop the retrieval model of the submerged vegetation or submerged vegetation-dominated hydrophyte biomass. The actual situation, however, is a little far away from the perfect assumptions, which may result in a big error for removing the contribution due to water column. So it is not applicable and practicable very well.
Accurate separation of the aquatic vegetation from lake area is the first step to accurately estimate the biomass using remote sensing images. In this paper, the aquatic vegetation zone in Taihu Lake was divided into the floating vegetation-dominated and submerged vegetation-dominated hydrophyte zones, which only is a way to improve biomass estimation precision based on the characteristic of different species of aquatic vegetation growing in a cluster from a viewpoint of remote sensing. In Taihu Lake, both of the water depth and transparency in aquatic vegetation zone have a distinct difference with in other zones. The prior knowledge on depth and transparency can help discriminate the confusion between vegetation and non-vegetation zones and then increase their separability and classification accuracy to a certain ex0tent; in the same way, the estimated biomass, especially for Type one by Equation (9) is more accurate than by Equation (10). It's a pity, however, that no such prior knowledge about transparency was provided in the corresponding date when the Landsat TM images used in this study was acquired. Although the historical transparency could be estimated using the empirical model as Equation (1) by the Landsat TM images, the estimated result could not be verified because of no data measured in situ. So we left off the transparency data in the historical biomass estimation using the Landsat TM images.
Maybe it is more perfect and more accurate for the total biomass if we may implement biomass estimation of each species of aquatic vegetations, which, however, is very difficult in Taihu Lake. Firstly, the aquatic vegetation grows in a cluster including several species, so there hardly is a pure vegetation pixel in the satellite imagery with a coarse spatial resolution of about 20-30 m, where only a single species of aquatic vegetation grows. The spectral signal of vegetation in one pixel of satellite imagery comes from all the aquatic vegetations in the range of field-of-view angle. Additionally, most of multi-spectral sensors have a coarse spectral resolution of about 60-90 nm in the visible wavelength, which results in a difficulty to differentiate vegetation species from a viewpoint of spectral characteristics. So the capability to differentiate the synthetical signal is essential here. Hyperspectral remote sensing is a good tool to implement the differentiation because of its narrow band width. However at present, there only is an on-orbit hyperspectral satellite sensor EO-1 Hyperion with the same spatial resolution of 30 m as Landsat TM, but a fine spectral resolution of about 10 nm. The Hyperion image with a swath of about 7 km needs to be ordered in a schedule time range. Covering all over the Taihu Lake needs ten Hyperion images, which, however, are in different paths with different acquiring time. So it is obvious up to now that the hyperspectral image acquisition is inconvenient and the sensor capability can not fully conform to the requirements of applications.
Phragmites communi is mainly regular distributed in the western and southern littoral zone with a range of about 100 m and its biomass is far more than other hydrophytes. According to our experience and previous investigation along the lake shore, there has only been a small change in its spatial distribution since 2001, so in this campaign, Phragmites communi was not investigated and the total biomass estimated in this paper doesn't contain of Phragmites communi.
It is very interesting to take notice of the phenomenon that the aquatic vegetation is disappearing gradually in some lake zones, such as Gongshan Bay, Xukou Bay and Dongshan Bay; however in other zones it is spreading. The reason why there is an opposite trend for the aquatic vegetation in the same lake is not clear now, which, however, is a very interesting and important problem for scientists in the field to discover change of the Taihu Lake ecosystem.

7. Conclusions

It is applicable and practical in technology and methodology for satellite remote sensing to monitor the spatial distribution of aquatic vegetations and their biomass in shallow inland lake. A decision tree with prior knowledge on water transparency may help improve the satellite remote sensing classification accuracy of aquatic vegetation to a certain extent. The empirical model for estimating aquatic vegetation biomass is more accurate and more practical than the semi-analytical model at present. The present empirical model could be conditionally transferred to develop the past model through the relationship between the vegetation index NDVI, respectively, coming from the present and the past satellite imageries.
In Taihu Lake, the total covering area of aquatic vegetation has decreased from 454.6 km2 in 2001 to 364.1 km2 in 2007; and vegetation in Gongshan Bay has vanished almost entirely, and a great decrease also took place in Xukou Bay and Dongshan Bay. However it is notewrthy that a large spread toward the west along the shore and concurrently toward the north was happening at the same time as the rapid total decrease. At present, the total biomass of aquatic vegetations, excluding Phragmites communi, is about 406 (±46.3) thousand tons, mainly distributed in Zhunhu Bay, Guangfu Bay, East Taihu Bay and its mouth zone, those zones joining Xishan island with Dongshan island, and some littoral zones along Qidou-Huzhou. However, it was about 528 (±70.2) thousand tons in 2004, 482 (±67.5) and 489 (±64.4) thousand tons in 2002 and 2001.

Acknowledgments

This research was jointly supported by grants from the National Natural Science Foundation of China (Grant No.40671138), the Natural Science Foundation of Jiangsu Province of China (Grant No. BK2007261), National Project of Scientific and Technical Supporting Programs Funded by Ministry of Science and Technology of China (Grant No. 2007BAC26B01). Qingjun Song, Jianjun Zhu, Xiaomei Wang, Dayong Bi, Jun He and Xiaozhou Sun are acknowledged for participating in the field campaigns and experiment analysis. And we also thank the three anonymous referees for their comments and suggestions to clarify a few points in this manuscript.

References and Notes

  1. Vis, C.; Hudon, C.; Carignan, R. An evaluation of approaches used todetermine the distribution and biomass of emergent and submerged aquatic macrophytes over large spatial scales. Aquat. Bot. 2003, 77, 187–201. [Google Scholar]
  2. Malthus, T.J.; George, D.G. Airborne remote sensing of macrophytes in Cefni Reservoir, Anglesey, UK. Aquat. Bot. 1997, 58, 317–332. [Google Scholar]
  3. Marshall, T.R.; Lee, P.F. Mapping aquatic macrophytes through digital image analysis of aerial photographs: an assessment. J. Aquat. Plant Manag. 1994, 32, 61–66. [Google Scholar]
  4. Penuelas, J.; Gamon, J.A.; Griffin, K.L.; Field, C.B. Assessing community type, plant biomass, pigment composition, and photosynthetic efficiency of aquatic vegetation from spectral reflectance. Remote Sens. Environ. 1993, 46, 110–118. [Google Scholar]
  5. Valta-Hulkkonen, K.; Pellikka, P.; Tanskanen, H.; Ustinov, A.; Sandman, O. Digital false colour aerial photographs for discrimination of aquatic macrophyte species. Aquat. Bot. 2003, 75, 71–88. [Google Scholar]
  6. Lehmann, A.; Jaquet, J.M.; Lachavanne, J.B. Contribution of GIS to submerged macrophyte biomass estimation and community structure modeling, Lake Geneva, Switzerland. Aquat. Bot. 1993, 47, 99–117. [Google Scholar]
  7. Lehmann, A.; Lachavanne, J.B. Geographic information systems and remote sensing in aquatic botany. Aquat. Bot. 1997, 58, 195–207. [Google Scholar]
  8. Valta-Hulkkonen, K.; Kanninen, A.; Pellikka, P. Remote sensing and GIS for detecting changes in the aquatic vegetation of a rehabilitated lake. Int. J. Remote Sens. 2004, 25, 5745–5758. [Google Scholar]
  9. Valta-Hulkkonen, K.; Pellikka, P.; Peltoniemi, J. Assessment of bidirectional effects over aquatic macrophyte vegetation in CIR aerial photographs. Photogramm. Eng. Remote Sens. 2004, 70, 581–587. [Google Scholar]
  10. Lee, K.H.; Lunetta, R.S. Wetland Detection Methods. In Wetland and Environmental Applications of GIS; Lyon, J.G., McCarthy, J., Eds.; CRC Press Inc.: Boca Raton, FL, USA, 1995; p. 373. [Google Scholar]
  11. Schweizer, D.; Armstrong, A.; Posada, J. Remote sensing characterization of benthic habitats and submerged vegetation biomass in Los Roques Archipelago National Park, Venezuela. Int. J. Remote Sens. 2005, 26, 2657–2667. [Google Scholar]
  12. Zhang, X. On the stimation of biomass of submerged vegetation using Landsat Thematic Mapper (TM) imagery: a case study of the Honghu Lake, PR China. Int. J. Remote Sens. 1998, 19, 11–20. [Google Scholar]
  13. Dottavio, C.L.; Dottavio, F.D. Potential benefits of new satellite sensors to wetland mapping. J. Photogramm. Eng. Remote Sens. 1984, 50, 599–606. [Google Scholar]
  14. Mumby, P.J.; Edwards, A.J. Mapping marine environments with IKONOS imagery: enhanced spatial resolution can deliver greater thematic accuracy. Remote Sens. Environ. 2002, 82, 248–257. [Google Scholar]
  15. Pinnel, N.; Heege, T.; Zimmermann, S. Spectral discrimination of submerged macrophytes in Lakes using hyperspectral remote sensing data. http://www.limno.biologie.tumuenchen.de/forschung/publika/Pinnel2004OO189.pdf 2006-11-20.
  16. Williams, D.J.; Rybicki, N.B.; Lombana, A.V.; O'Brien, T.M.; Gomez, R.B. Preliminary investigation of submerged aquatic vegetation mapping using hyperspectral remote sensing. Environ. Monit. Assess. 2003, 81, 383–392. [Google Scholar]
  17. Ackleson, S.G.; Klemas, V. Remote sensing of submerged aquatic vegetation in lower Chesapeake Bay: A comparison of Landsat MSS to TM imagery. Remote Sens. Environ. 1987, 22, 235–248. [Google Scholar]
  18. Havens, K.E. Submerged aquatic vegetation correlations with depth and light attenuating materials in a shallow subtropical lake. Hydrobiologia 2003, 493, 173–186. [Google Scholar]
  19. Matlthus, T.J.; George, D.G. Airborne remote sensing of macrophytes in Cefni Reservoir, Anglesey, UK. Aquat. Bot. 1997, 58, 317–332. [Google Scholar]
  20. Karpouzli, E.; Malthus, T.; Place, C.; Chui, A.M.; Garcia, M. I.; Mair, J. Underwater light characterization for correction of remotely sensed images. Int. J. Remote Sens. 2003, 24, 2683–2702. [Google Scholar]
  21. Mumby, P.J.; Clark, C.D.; Green, E.P.; Edwards, A.J. Benefits of water column correction and contextual editing for mapping coral reefs. Int. J. Remote Sens. 1998, 19, 203–210. [Google Scholar]
  22. Ma, R.; Tang, J.; Dai, J. Bio-optical model with optimal parameter suitable for Taihu Lake in water colour remote sensing. Int. J. Remote Sens. 2006, 27, 4303–4326. [Google Scholar]
  23. Ma, R.; Tang, J.; Dai, J.; Zhang, Y.; Song, Q. Absorption and scattering properties of water body in Taihu Lake, China: absorption. Int. J. Remote Sens. 2006, 27, 4275–4302. [Google Scholar]
  24. Huang, X. Eco-Investigation, Observation and Analysis of Lakes; Standard Press of China: Beijing, China, 1999; pp. 77–99. [Google Scholar]
  25. Mueller, J.L.; Fargion, G.S.; Mcclain, C.R. Ocean Optics Protocols for Satellite Ocean Color Sensor Validation; Revision 4; Greenbelt: Maryland, USA, 2003. [Google Scholar]
  26. Pope, R.M.; Fry, E.S. Absorption spectrum (380–700nm) of pure water II. Integrating cavity measurements. Appl. Opt. 1997, 36, 8710–8722. [Google Scholar]
  27. Lee, Z.; Carder, K.L.; Mobley, C.D.; Steward, R.G.; Path, J.S. Hyperspectral remote sensing for shallow waters: 2. Deriving bottom depths and water properties by optimization. Appl. Opt. 1999, 38, 3831–3843. [Google Scholar]
  28. Qin, B.; Hu, W.; Chen, W. Process and Mechanism of Environmental Changes of the Taihu Lake; Science Press: Beijing, China, 2004; pp. 10–13. [Google Scholar]
  29. SASG (Survey Adjustment Subject Group of Survey College in Wuhan University). Error theorem and foundations of survey adjustment; Wuhan University Press: Wuhan, China, 2003; pp. 89–102. [Google Scholar]
Figure 1. Location of Taihu Lake in China.
Figure 1. Location of Taihu Lake in China.
Sensors 08 03988f1
Figure 2. Locations of field samples investigated on 10-19 June 2007; there is aquatic vegetation in the black circle samples and there isn't in the cross samples.
Figure 2. Locations of field samples investigated on 10-19 June 2007; there is aquatic vegetation in the black circle samples and there isn't in the cross samples.
Sensors 08 03988f2
Figure 3. Comparisons between remote sensing reflectance only due to water (the thin gray curves), reflectance due to Type one (the thick black curves), reflectance due to Type two (the thin black curves), all of which were measured in situ. The x-axes is wavelength (in nm), and the y-axes is reflectance (dimensionless).
Figure 3. Comparisons between remote sensing reflectance only due to water (the thin gray curves), reflectance due to Type one (the thick black curves), reflectance due to Type two (the thin black curves), all of which were measured in situ. The x-axes is wavelength (in nm), and the y-axes is reflectance (dimensionless).
Sensors 08 03988f3
Figure 4. Relationships between the Secchi disk transparencies measured in situ and the atmosphere-corrected digital number values of: (a) the first band (B1) mounted at LISS-3 sensor, and (b) the second band (B2). Both of the y-axes are the Secchi disk transparencies (in cm), x-axes in (a) and (b), respectively, are B1 and B2 (in dimensionless).
Figure 4. Relationships between the Secchi disk transparencies measured in situ and the atmosphere-corrected digital number values of: (a) the first band (B1) mounted at LISS-3 sensor, and (b) the second band (B2). Both of the y-axes are the Secchi disk transparencies (in cm), x-axes in (a) and (b), respectively, are B1 and B2 (in dimensionless).
Sensors 08 03988f4
Figure 5. Sketch map of the hypothesis condition.
Figure 5. Sketch map of the hypothesis condition.
Sensors 08 03988f5
Figure 6. Spatial distribution of the aquatic vegetation except Phragmites communi on 26 July 2001, 13 July 2002, 26 July 2004, and 16 June 2007; Type one and Type two, respectively, are the floating vegetation-dominated hydrophyte, and the submerged vegetation-dominated hydrophyte.
Figure 6. Spatial distribution of the aquatic vegetation except Phragmites communi on 26 July 2001, 13 July 2002, 26 July 2004, and 16 June 2007; Type one and Type two, respectively, are the floating vegetation-dominated hydrophyte, and the submerged vegetation-dominated hydrophyte.
Sensors 08 03988f6
Figure 7. The satellite-estimated biomass of aquatic vegetation on 26 July 2001, 13 July 2002, 26 July 2004, and 16 June 2007.
Figure 7. The satellite-estimated biomass of aquatic vegetation on 26 July 2001, 13 July 2002, 26 July 2004, and 16 June 2007.
Sensors 08 03988f7
Figure 8. Scatter-map of biomass between the measured in situ (x-axis, in kg/m2) and the estimated by remote sensing (y-axis, in kg/m2). The thick line is the fitting curve by a least square method, and the thin is the 1:1 line.
Figure 8. Scatter-map of biomass between the measured in situ (x-axis, in kg/m2) and the estimated by remote sensing (y-axis, in kg/m2). The thick line is the fitting curve by a least square method, and the thin is the 1:1 line.
Sensors 08 03988f8
Table 1. Description of the six clusters and their concomitances.
Table 1. Description of the six clusters and their concomitances.
Vegetation clusterMain concomitant vegetation
Potamogeton malaianus clusterVallisneria natanus, Ceratophyllum denersum, Hydrilla verticillata
Nymphoides peltatum clusterPotamogeton malaianus, Trapa incise, Vallisneria natanus
Nymphoides peltatum-Potamogeton malaianus clusterTrapa incise, Ceratophyllum demersum, Vallisneria natanus
Vallisneria natanus clusterCeratophyllum demersum, Elodea nuttalli, Hydrilla verticillata
Vallisneria natanus-Elodea nuttalli clusterNymphoides peltatum and Ceratophyllum demersum
Najas marina clusterPotamogeton malaianus, Vallisneria natanus
Table 2. The covering area of aquatic vegetation and its biomass in Taihu Lake, no containing Phragmites communi; Type one and Type two, respectively, are the floating vegetation-dominated hydrophyte, and the submerged vegetation-dominated hydrophyte.
Table 2. The covering area of aquatic vegetation and its biomass in Taihu Lake, no containing Phragmites communi; Type one and Type two, respectively, are the floating vegetation-dominated hydrophyte, and the submerged vegetation-dominated hydrophyte.
DateType one (km2)Type two (km2)Total area (km2)Total biomass (thousand ton)
15 June 200772.4291.7364.1406
26 July 200489.1393.1482.2528
13 July 200271.7380.0451.7482
26 July 200175.8378.8454.6489

Share and Cite

MDPI and ACS Style

Ma, R.; Duan, H.; Gu, X.; Zhang, S. Detecting Aquatic Vegetation Changes in Taihu Lake, China Using Multi-temporal Satellite Imagery. Sensors 2008, 8, 3988-4005. https://doi.org/10.3390/s8063988

AMA Style

Ma R, Duan H, Gu X, Zhang S. Detecting Aquatic Vegetation Changes in Taihu Lake, China Using Multi-temporal Satellite Imagery. Sensors. 2008; 8(6):3988-4005. https://doi.org/10.3390/s8063988

Chicago/Turabian Style

Ma, Ronghua, Hongtao Duan, Xiaohong Gu, and Shouxuan Zhang. 2008. "Detecting Aquatic Vegetation Changes in Taihu Lake, China Using Multi-temporal Satellite Imagery" Sensors 8, no. 6: 3988-4005. https://doi.org/10.3390/s8063988

APA Style

Ma, R., Duan, H., Gu, X., & Zhang, S. (2008). Detecting Aquatic Vegetation Changes in Taihu Lake, China Using Multi-temporal Satellite Imagery. Sensors, 8(6), 3988-4005. https://doi.org/10.3390/s8063988

Article Metrics

Back to TopTop