4 March 2015 Poyang Lake wetland vegetation biomass inversion using polarimetric RADARSAT-2 synthetic aperture radar data
Author Affiliations +
J. of Applied Remote Sensing, 9(1), 096077 (2015). doi:10.1117/1.JRS.9.096077
Poyang Lake is the largest freshwater lake in China and one of the most important wetlands in the world. Vegetation, an important component of wetland ecosystems, is one of the main sources of the carbon in the atmosphere. Biomass can quantify the contribution of wetland vegetation to carbon sinks and carbon sources. Synthetic aperture radar (SAR), which can operate in all day and weather conditions and penetrate vegetation to some extent, can be used to retrieve information about vegetation structure and the aboveground biomass. In this study, RADARSAT-2 polarimetric SAR data were used to retrieve aboveground vegetation biomass in the Poyang Lake wetland. Based on the canopy backscatter model, the vegetation backscatter characteristics in the C-band were studied, and a good relation between simulated backscatter and backscatter in the RADARSAT-2 imagery was achieved. Using the backscatter model, pairs of training data were built and used to train the back propagation artificial neural network. The biomass was retrieved using this ANN and compared with the field survey results. The root-mean-square error in the biomass estimation was 45.57  g/m2. This shows that the combination of the model and polarimetric decomposition components can efficiently improve the inversion precision.
Shen, Liao, Guo, and Liu: Poyang Lake wetland vegetation biomass inversion using polarimetric RADARSAT-2 synthetic aperture radar data



Biomass, in ecology, is the mass of living organisms in a given area or ecosystem at a given time.1 Here, as in most of the literature on the carbon cycle and remote sensing, biomass refers to the total dry weight of organic aboveground plant matter in a unit surface area. Wetlands, an important component of terrestrial ecosystems, play a key role in global climate change because one of their components, wetland vegetation, is the main source of the carbon in the atmosphere.2 Wetland vegetation biomass is a key indicator for evaluating the carbon sequestration capacity of wetlands. Biomass governs the potential amount of carbon that could be released into the atmosphere because of deforestation, and regional biomass changes have been associated with important outcomes in ecosystem functional characteristics and climate change.3 Therefore, biomass estimation of the wetland vegetation plays a key role in the understanding of dynamic changes of wetland ecosystems. Generally, conventional field methods of biomass estimation are time consuming and labor intensive and cannot provide details on the spatial distribution of a biomass. The advantages of remote sensing data, especially synthetic aperture radar (SAR) data, are that it can be collected repeatedly and processed conveniently, making it an efficient data source for large-area biomass retrieval, especially in areas of difficult access. Lu pointed out that remote sensing data have been the main data source for biomass estimation and are attracting increased attention from scientific researchers.3

Optical data can be used to extract spectral information in the form of the vegetation index, which is related to the biomass and can be used to retrieve biomass by application of a regression model.4,5 However, in areas of dense vegetation, saturation in optical bands influences the precision in this biomass inversion.5 Because SAR is able to penetrate vegetation and to operate at all times and in all types of weather and because it saturates only at very high levels of energy, it can be used to retrieve vegetation biophysical parameters (VBPs).6

Many studies have shown that SAR can be used to retrieve vegetation biomass. Henderson and Lewis reviewed the recent research from five major remote sensing journals, spanning the years 1965 to 2007, and found that the C-band is necessary for herbaceous wetland detection. HH polarization is preferable to VV, but cross-polarized data seem to contribute more than was thought earlier in some cases.7 Kasischke et al. concluded that imaging radar was suitable for ecological applications, especially the delineation of wetland inundation, vegetative cover, and the measurement of aboveground plant biomass.8 Lu reviewed radar data for biomass estimation and discussed the common saturation problem in radar data.3 Baghdadi et al. have evaluated the use of C-band polarimetric SAR (PolSAR) data for wetland mapping,9 and the results have shown that the HH and cross-polarizations are better than VV polarization and the cross-polarization data can provide the best separation. Srivastava et al. explored the potential of multiparametric (multi-incidence angle, multipolarized, temporal, high-resolution, and multifrequency) SAR data in wetland inventories and suggested that SAR data have a number of unique features that can be exploited in the assessment, monitoring, and management of a wetland ecosystem. Srivastava et al. also addressed the interaction of radar signals (frequency, incidence angle, and polarization) with various components of a wetland ecosystem (moisture content of vegetation or soil, surface roughness, structure of the vegetation, textural properties of the target, orientation of the target, etc.).6

Liao et al. and Shen et al. have retrieved the Poyang Lake wetland vegetation biomass using ENVISAT ASAR data with an artificial neural network (ANN).1011.12 Shen et al. and Yang et al. found that the fresh biomass of rice has a good relation with the backscatter coefficient of ENVISAT ASAR data with HH and HV polarizations.13,14 Because ASAR data are nonpolarimetric, the polarized difference and polarization ratio were used to improve the precision of the inversion procedure.11 To optimize livestock management in a water-saturated Andean wetland, Moreau and Toan used ERS SAR data to retrieve spatio-temporal information on the biomass of totora reeds and bofedal, which are a critical foraging resource for smallholders.15 The radar backscattering coefficient (σ0 measured in dB), measured by ERS, was found to be sensitive to both the humid and the dry biomass of totora reeds and bofedal grasslands. The sensitivity of the signal-to-biomass variation is high for dry biomass, less than 1kg/m2 for totora, and less than 2kg/m2 for bofedal.15 Yang et al. used relational models based on the analysis of the correlation between backscatter coefficients and field measurements to retrieve rice growth parameters from ASAR images.14

Polarimetric information can significantly improve wetland vegetation characterization for wetland mapping.16,17 Slatton et al. modeled wetland vegetation by using airborne PolSAR in herbaceous coastal wetlands in the C-, L-, and P-bands and showed that PolSAR has the potential for mapping coastal herbaceous wetlands.18 The short-wavelength SAR (X- and C-bands) interacts with the upper part of the crop canopy, thus offering the potential for retrieving crop biophysical parameters.19 The use of polarimetric RADARSAT-2 data has been investigated for operational mapping and for monitoring Canadian wetlands16,20 and the Poyang Lake wetland.21,22 Inoue et al. found that the backscatters from a polarimetric scatterometer are sensitive to rice biomass and have some connection with the angular distribution and size of the plant elements (leaves, stems, and heads).23 Target scattering decomposition has become a popular method for extracting the geophysical parameters of a natural target from PolSAR data,24 which can lead to unique parameters that optimize wetland classification.2021.22 Li et al. proposed a new inversion method for vegetation parameters based on the Freeman polarization decomposition.25

Retrieving wetland vegetation biomass from SAR data has often been carried out by implementation of linear and nonlinear multiple regressions. However, the performance is influenced by the highly complex nonlinear relation between vegetation biomass and the SAR backscatter. An alternative approach, based on ANN algorithms, which are becoming commonly available and readily usable, has been devised to retrieve environmental parameters from SAR data. ANNs are composed of many nonlinear computational elements such as neurons, operating in parallel and linked with each other through connections characterized by multiplying factors. This structure makes ANNs inherently suitable for addressing nonlinear problems.26 The use of these algorithms in remote sensing has often been found effective because they can simultaneously handle nonlinear mapping of a multidimensional input space onto the output space and also cope with complex statistical behavior. Jin and Liu used the experimental data obtained from active/passive microwave remote sensing and from ground-truth measurements to train the ANN to retrieve biomass parameters.27 The ANN model was first trained by using wheat canopy data from 1988, and was then used to retrieve biomass parameters of wheat canopy during growth in l989, which was compared with the ground-truth data. Frate and Wang applied an electromagnetic model and neural network algorithms to retrieve the biomass for sunflower fields using radar backscattering data. The electromagnetic model was used to generate the scattering coefficients for training and testing the network. The inversion results showed that the neural network was capable of performing the retrieval with good accuracy. By optimizing the structural complexity of the net, these researchers obtained a better inversion result.26 Frate and Wang used ANN to retrieve forest biomass from multifrequency (L- and P-bands), multipolarization (HH, VV, and HV) backscattering, after the training and pruning procedures had been made to the net, and compared the ANN retrieval accuracy with that yielded by linear and nonlinear regressions and by a model-based technique.28

In this study, polarimetric RADARSAT-2 data were used to retrieve vegetation biomass in the Poyang Lake wetland. First, the backscatter from vegetation was simulated on the basis of the canopy backscatter model. Then, the Freeman polarimetric decomposition components and HH, VV, and HV polarization backscatters from RADARSAT-2 data were used as the inputs to an ANN to retrieve the VBPs. Finally, the retrieved VBPs and Freeman polarimetric decomposition components were used as the inputs to another ANN to retrieve the biomass.


Study Area and Data


Study Area

The test site is located in the Poyang Lake wetland (inland lake wetland), Jiangxi Province, China in the lower Yangtze River Basin. The latitude and longitude are 28°22′–29°15′ N and 115°47′–116°45′ E (Fig. 1). The climate is characterized as a subtropical, humid, monsoon climate with a 1620 mm mean annual precipitation, and an annual average temperature of approximately 17°C. In summer, Poyang Lake is the largest freshwater body in China. By the end of the rainy season, the lake can extend up to 3500km2 in area, but during the dry season, it may shrink to less than 1000km2 and become a system of sublakes interspersed with mudflats, sediment beds, and vegetation. The variation in the water level of Poyang Lake during the year is very large. In 2007, the absolute fluctuation of the water level was up to 8.37 m. The area of the lake greatly varies with the fluctuation of its water level. The water level rises during the flood season, and the water surface suddenly expands. The water level drops during the dry season, and the lakebed appears and only a few wandering watercourses remain. In summary, the test site has dry and wet seasons. In the dry season (November to April), the wetland vegetation is above the water. It starts to rapidly grow from early spring, with the biomass reaching its highest level in April.29 In the wet season (June to September), the wetland vegetation is flooded and hardly grows.

Fig. 1

Spatial location and extent of the study area of Poyang Lake, Jiangxi Province, People’s Republic of China. The red polygon shows the extent of Poyang Lake wetland, the black polygon shows the extent of the RADARSAT-2 imagery acquired on April 8, 2011, and the green polygon shows the Nanjishan National Nature Reserve. (a) The Pauli decomposition image, (b) and (c) Landsat ETM+ acquired on February 4, 2003.


The predominant vegetation in Poyang Lake is carex and reed [Fig. 4(a)], which together cover more than 90% of the area.29,30 The height, water content, and plant density of the vegetation vary with the different hydrological stages from April to November. In April, the carex rapidly grows with green leaves of approximately 30 to 70cm2. In July, the plants are submerged in the water, and in November they turn yellow and the water content of the plants decreases.

The study area (the RADARSAT-2’s coverage) is mainly located in the southwest part of Poyang Lake wetland known as the Nanjishan National Nature Reserve (the green polygon in Fig. 1), covering an area of 33,300hm2. The land cover of the Nanjishan Wetland Nature Reserve consists of a profundal zone, swamp, marsh, sandy land, and grassland, where the predominant vegetation is also carex and reed.31



RADARSAT-2, carrying a fully polarimetric SAR, has opened up new opportunities for wetland classification by the use of PolSAR images.20 In this study, one scene of RADARSAT-2 data was acquired on April 8, 2011, which was during the dry season at Poyang Lake. In fact, at this time, the area was experiencing the most severe drought of the past 50 years. The main system parameters are listed in Table 1. In terms of pixels, the data set had the dimensions of 3080×5982 in range and azimuth, respectively.

Table 1

Main RADARSAT-2 system parameters.

RADARSAT-2 parameterValue
Acquisition dateApril 8, 2011
Frequency bandC-band
Central frequency5.405 GHz
ModeFine Quad Pol (HH, HV, VH, and VV)
Pass directionAscending
Range and azimuth pixel spacing4.73 m, 4.96 m
Resolution8 m
Incidence angle31.31°–33.00°


Field Survey

Because ground-truth data collection in a wetland is difficult to access, time consuming, and tedious,32 a quasisynchronous field survey was conducted in Poyang Lake wetland from April 7, 2011 to April 10, 2011. Within the area covered by the RADARSAT-2 imagery, 2 to 3 sample points (0.5×0.5m2 or 1×1m2) were randomly selected for each sample area, and 54 sample points within different sample areas (with intersite spacing >50m, 10 times the pixel size) were collected (Fig. 1).32 At each sample point, the fresh biomass, total height, stem height, stem density, stem diameter, leaf length, leaf width, leaf thickness, leaf number per plant, and the longitude and latitude of the sample point were recorded. After the grass was harvested, it was immediately weighed to determine the fresh biomass. The grass was then dried in an oven for 30 min at 120°C, and then dried to constant weight at 80°C to get the dry biomass and the moisture content.33 Because of the severe drought in April 2011,34 the dry biomass was a little lighter than the sample weight in the survey conducted from March 29 to April 3, 2007.11

According to ground-truth data, the dry biomass was well correlated with the vegetation height, stem-layer height, leaf length, leaf width, and stem diameter [Figs. 2(a) to 2(e)], but not so well with leaf number density [Fig. 2(f)]. The dry biomass negatively correlated with the vegetation moisture content [Fig. 2(g)]. Therefore, these parameters (the vegetation height, stem-layer height, leaf length, leaf width, stem diameter, and vegetation moisture) can be used to retrieve the dry aboveground biomass of the vegetation (The following discussion will be based on the dry aboveground biomass).

Fig. 2

Relation between the dry biomass and various vegetation biophysical parameters, (a) for vegetation height, (b) for stem layer height, (c) for leaf length, (d) for leaf width, (e) for stem diameter, (f) for leaf number density, and (g) for vegetation moisture content.





RADARSAT-2 Data Processing

First, the RADARSAT-2 image was calibrated to the Sinclair matrix by use of the sigma look-up table.35 Next, the backscatter coefficients (σ0) for HH, HV/VH, and VV polarizations, as well as the coherency elements, T3, were calculated. To reduce the speckle, the data were filtered through the Lee refined filter36 with a window size of 3×3. The Freeman polarimetric decomposition assumes that the covariance matrix of the backscattering data can be decomposed into three covariance matrices, corresponding to scattering mechanisms of the following types: direct surface (odd-bounce), dihedral-type (double and even-bounce), and a random volume.37 Thus, the Freeman three-component decomposition method was used to extract the components representing the volume, double-bounce, and surface scattering (Fig. 3).37,38 From the top left of Fig. 3, it can be seen that the predominant backscattering for wetland vegetation is volume scattering (Region A). However, at the lakeshore (Region B), the vegetation displays as red or yellow because of the high double-bounce scattering between the sparse vegetation and the ground.39 The predominant scattering for the mudbank (Region C) near the open water is surface scattering. Surface scattering is dominant for some of the water (Region D) because of the flotsam, but in some ponds (Region E), surface scattering (specular reflection) is dominant and results in the weak backscatter.

Fig. 3

Composite map of Freeman polarimetric decomposition.


All the data (the σ0 under HH, HV/VH, VV polarizations; the three components under Freeman decomposition) were projected to the same projection, which could be easily used to extract the sample points’ backscatter coefficients and the polarimetric decomposition components by use of their coordinates. To determine the vegetation extent, the unsupervised H-A-alpha-Wishart classification method40,41 was used to classify the extent of the water body and lakeshore. The classification result can be derived from Fig. 8, which divided the land cover into four types: grassland with different vegetation biomass level, water, lakeshore, and others (NULL in Fig. 8).


Backscatter Simulation

The canopy backscatter model used in this work42 was modified on the basis of Karam’s forest canopy scattering model,43 which fits the scattering by rice at different growth stages. The model was used to simulate the wetland vegetation’s backscatter. This simulation could be done because the model has a similar structure to rice but has a greater density (Fig. 4). Zhang used this modified model to retrieve the rice parameters.44 Here, the backscatter model was used to analyze the radar backscatter sensitivity to the wetland VBPs.

Fig. 4

(a) Poyang Lake wetland vegetation and (b) the diagram used in the backscatter model.


Because of the high vegetation coverage, it was assumed that the only vegetation within the study area was carex. There were three layers: (1) the leaf layer, (2) the trunk layer, and (3) the ground. There was no ear layer. Thus, the total backscatter can be expressed as


where σtotal is the total backscatter, σleaf and σstem are the volume backscatters from leaf and stem, σleaf-ground and σstem-ground are the double-bounce backscatters introduced by the interaction between the leaf, stem, and the ground, and σground is the direct backscatter from the ground. In this model, the Debye−Cole model was used to calculate the dielectric constant of the vegetation canopy.45 The inputs to the backscatter model are listed in Table 2.

Table 2

Statistics for the wetland vegetation parameters (inputs for the backscatter model).

NumberParameterMeanMaxMinStd. Dev.
1Leaf half length (cm)13.3220.427.582.78
2Leaf half width (mm)2.052.751.330.28
3Leaf half depth (thickness, mm)
4Leaf number density (11,80118,92845363688
5Vegetation moisture content (%)79.5183.9175.161.82
6Vegetation height (cm)48.9567.6737.676.01
7Stem-layer height (cm)15.3929.179.84.08
8Stem diameter (mm)2.222.881.410.33
9Stem density (20963548648669
10Wet aboveground biomass (1650.12363.6870.3368.5
11Dry aboveground biomass (348.8571.2181.080.97

Because of the wetland environment with high soil moisture (nearly saturated), the ground could be considered flooded in most cases. The flooded ground results in weak direct backscatter from the ground (water) because of specular reflection, but the ground can still contribute to the double-bounce between it and the leaf and stem, whereas the nonflooded ground can also contribute to the two-way attenuation (SAR signal coming to ground and then going back to the sensor through the vegetation).46,47 In addition, because of the same growth stage when the field survey was conducted, the vegetation moisture scatters closely around 80%, as did the vegetation moisture from the former field surveys.11 However, the different canopy height and density resulting from the different fertilizer of underlayed soil can affect the backscatter to some extent.

Shen et al. pointed out that the incidence angle has little effect on the backscatter in ENVISAT ASAR (C-band) data.11 Because the wavelength in this study was the same and because of the small range of incidence angles (31.31°N to 33.00°N), the effect of the incidence angle was ignored here.

The simulated backscatter coefficients were calculated by use of the field-sample parameters. Compared with the corresponding backscatter coefficients extracted from the RADARSAT-2 image, the simulation error was less than 1 dB, mostly ranging from 0.2 to 0.5 dB, and the R2 was 0.93, 0.94, and 0.95 for HH, VV, and HV polarizations, respectively (Fig. 5). The span of the backscatter coefficients was approximately 6 dB, which can show the different statuses of land covers.

Fig. 5

Comparison of backscatter coefficients simulated by the model and those extracted from the RADARSAT-2 image, (a) for HH polarization, (b) for VV polarization, and (c) for HV polarization.



Sensitivity Analysis

After simulation of the backscatter by use of the backscatter model, the backscatter data were used to analyze the relation between the VBPs and the scatter components, volume scattering, double-bounce scattering, and surface (single) scattering (Fig. 6). The sum of double-bounces (σleaf-ground,σstem-ground) between leaf, stem, and ground in Eq. (1) was taken as the double-bounce scatter in the Freeman polarimetric decomposition, the sum of σleaf and σstem was taken as the volume scatter, and the σground was taken as the surface scatter. The supposed flooded ground results in weak direct backscatter (surface scatter) from the ground, so only the sensitivities of volume scatter and double-bounce scatter (simulated using the backscatter model) to various vegetation biophysical parameters were analyzed, and then used to retrieve the vegetation biomass. When the backscatter was calculated, one of the VBPs was taken as variable (ranging from its minimum to maximum in Table 2), and the others were taken as invariable with their means (Table 2) as the value.

Fig. 6

The sensitivity of scatter components (simulated by use of the backscatter model) to various vegetation biophysical parameters, volume scatter to (a) leaf half length, (b) leaf half width, (c) leaf half depth, (d) leaf number density, (f) leaf layer height, and (g) stem radius; (e) double bounce scatter to leaf number density, and (h) stem radius.


From Fig. 6, it can be seen that with increasing leaf length, leaf width, leaf depth, leaf density, leaf-layer height, and stem radius, the volume scatter increases for HH, VV, and HV polarizations. The backscatter for HH and VV polarizations is similar and is stronger than for HV polarization. Because of the depolarization of the vegetation canopy, the HV backscatter increases with increasing volume scatter.

As leaf number density increases, the double-bounce scatter decreases because of the increasing attenuation caused by the vegetation leaf layer. The decreasing trends of the double-bounce scatter for HH, VV, and HV polarizations are different, decreasing for VV and HV polarizations but slightly stable for HH polarization. VV polarization is sensitive to double-bounce scatter. For double-bounce scatter, the VV polarization has higher backscatter than HH polarization, which shows that VV polarization can penetrate the canopy more deeply.48



Section 2.3 shows that the biophysical parameters, including leaf length, leaf width, leaf-layer height, leaf number density, and the stem diameter, are sensitive to both the dry biomass and the polarimetric decomposed components (the volume scattering and double-bounce scattering). As a final step, these parameters were chosen as the outputs and inputs for the BP ANNs to retrieve the biomass. The flowchart for the biomass inversion is shown in (Fig. 7). The flow can be divided into two parts, one for the inversion of the VBPs and the other for the biomass.

Fig. 7

The flowchart for biomass inversion using the polarimetric RADARSAT-2 image and the artificial neural network (ANN) method. The left-hand part is the ANN for vegetation parameter inversion, including training and inversion; the right-hand part is the ANN for vegetation biomass inversion, including training and inversion.


Here, the BP ANNs, standard feedforward multilayer perceptrons,26 were used to retrieve the vegetation biomass (Fig. 7). The BP ANNs were both designed with two hidden layers. The number of neurons in the hidden layers was calculated according to n=i×j+k/2,49 where i and j are the numbers of neurons in the input and output layers, and k is the number of training samples. The transfer functions (activation function26) for BP ANNs are the logarithm S transfer function {a=logsig(n)=1/[1+exp(n)]} and the linear transfer function [a=purelin(n)=n]. The training function for BP ANNs is the “traingd” training function, which updates weight and bias values according to gradient descent.


Biomass Inversion

First, the simulated volume scatter, double-bounce scatter, and backscatter for HH, VV, and HV polarizations obtained from the backscatter model were taken as the inputs, and the five surveyed VBPs were taken as the outputs to train the ANN. After training of this ANN, the Freeman polarimetric decomposition components (volume scatter, double-bounce scatter) and the backscatters for HH, VV, and HV polarizations from the RADARSAT-2 data were used to retrieve the biophysical parameters for the whole study area.

Second, the ANN for the biomass inversion was trained by use of the surveyed VBPs, the simulated volume, and double-bounce scatter as inputs and the biomass as output. By taking the retrieved VBPs and the Freeman polarimetric decomposition components (volume scatter, double-bounce scatter) from the RADARSAT-2 data as the inputs, the biomass for the whole study area was retrieved (Fig. 8).

Fig. 8

Biomass map retrieved from the RADARSAT-2 image.


The vegetation is mainly distributed along the shorelines on higher ground. The higher the elevation, the more biomass there is.


Accuracy Assessment

The survey data consisting of the 54 samples were divided into two groups, one for training the ANN and the other for accuracy assessment. The accuracy of the wetland vegetation estimate was measured by use of the root-mean-square error (RMSE):


where Bg is the biomass according to the ground truth, Be is the biomass according to the RADARSAT-2 PolSAR-based estimate, and n is the number of observations.

The RMSE for polarimetric decomposition-based biomass inversion was 45.57g/m2, which is 13.06% of the average biomass (348.8g/m2) obtained from the 54 surveyed data. The coefficient of determination, R2, between the surveyed biomass and the retrieved biomass is 0.87 (Fig. 9).

Fig. 9

Comparison of model-simulated and field-measured biomass values.




For the Poyang Lake vegetated (mainly carex) wetland, the total backscattering is dependent on the interaction of the microwave energy with both the canopy and the canopy-ground. Both the characteristics of the canopy, such as density, distribution, orientation, shape of the foliage, dielectric constant, height, and branches and the characteristics of the sensor, such as polarization, incidence angle, and wavelength are important in determining the amount of radiation backscattered toward the radar antenna.50

RADARSAT-2 PolSAR data appear to be well suited in mapping biomass in this area because the C-band (5.4 GHz) radar data are particularly sensitive to vegetation biomass when the canopy is above an underlying water surface or a water-saturated soil. This relation occurs because the dominant scattering mechanisms involve vegetation-water surface interactions.15 PolSAR data can provide more information for land cover classification and biomass inversion. The dependency of the polarimetric scattering features on the VBPs was investigated. Because of the distributed vegetation in the Poyang Lake wetland and the statistics feature around the field survey points, the incoherent Freeman polarimetric decomposition method was used to extract the polarimetric scattering features. The results confirm that the Freeman polarimetric decomposition components can be used to estimate the VBPs because of their sensitivity to the vegetation structure parameters.

Here, the backscatter model was used to analyze the relation between the backscatter and various VBPs. The results show that the backscatter model can be used to simulate the backscatter for wetland vegetation and to generate the training data for the ANN because wetland vegetation has the same structure as rice, for which the coefficients of determination are 0.93, 0.94, and 0.95 for HH, VV, and HV polarizations, respectively. However, SAR backscatter from vegetation-covered fields is a strong function of dielectric properties of the vegetation, vegetation canopy structure, vegetation volume along with moisture, and surface roughness of the underlying soil.46,51,52 Sometimes SAR sensors may receive strong backscatter because of the double-bounce between the water surface and the vegetation stem.53 Here, the soil properties were not considered in this study, because most of the sample locations were either flooded or oversaturated during image acquisition. Therefore, when the backscatter model is used to simulate the backscatter from vegetation, the backscatter will be underestimated because of the absence of the soil surface. The biomass will be wrongly estimated because of the assumption that the ground is covered by water in the backscatter model. If the ground soil was not flooded or oversaturated, the soil properties should be taken into account. If soil properties had been considered in the modeling process, the accuracy of the backscatter simulation and biomass inversion would have been expected to be improved. When the backscatter model is used, another factor that must be taken into consideration is the geometry (shape, size, geometric distribution, etc.) of the leaf. Here, the probability of geometric distribution function used in the model was adopted from Liu et al.,54 where the vegetation and study area are the same types as in this study.

Because of the complex nonlinear relation between radar backscatter and wetland VBPs, the responses of scatter components to every single VBP show poor linear relation (Fig. 6). The occupied volume and moisture of each layer must both be taken into account for each layer when wetland VBPs and vegetation biomass are retrieved by use of SAR data.51 Here, the BP ANN was used to retrieve wetland VBPs and vegetation biomass because of the ability to express the nonlinear relation, which the traditional empirical and semi-empirical methods cannot do. When biomass was retrieved, the VBPs and the surveyed biomass were taken into account together. As a result, the ANN achieved good results, with R2 on the order of 0.88 between the surveyed biomass and the retrieved biomass. However, the ANN was not optimized to the optimal solution, hence more work had to be done for the simulation and the optimization of the ANN. When remote sensing technology is used to retrieve vegetation biomass, the data source and methods influence the inversion precision to a great degree.55 During training of the ANN, the upper and lower limits came from the survey data. Because they were not the actual values for the wetland vegetation, their use influenced the inversion accuracy.

The maximum values of vegetation biomass obtained during the field survey were 1650.1g/m2 (wet) and 348.8g/m2 (dry), which do not reach the saturation level for C-band SAR.3 The RMSE of overall biomass estimation for this area was better than that obtained by use of ENVISAT ASAR data in other studies.11,54 The RMSE for polarimetric decomposition-based biomass inversion was 45.57g/m2, which is 13.06% of the average biomass (348.8g/m2) obtained from the survey data. The coefficient of determination between the surveyed biomass and the retrieved biomass was 0.87 (Fig. 9).



This study focused on the application of ANNs combined with a backscatter model to retrieve wetland vegetation biomass with RADARSAT-2 polarimetric SAR data. The backscatter model was used to analyze the relation between the backscatter and various VBPs. The results show that the Freeman polarimetric decomposition components are sensitive to the wetland VBPs and can be used to retrieve the VBPs for the Poyang Lake wetland. One trained ANN, taking the simulated backscatter, volume scatter, and double-bounce scatter as inputs and five VBPs as outputs, was used to retrieve the wetland VBPs. The retrieved VBPs, combined with polarimetric decomposition components from RADARSAT-2 data, were used to retrieve wetland vegetation biomass by use of another trained ANN. The RMSE for polarimetric decomposition-based biomass inversion was 45.57g/m2, which is 13.06% of the average biomass. The coefficient of determination between the surveyed biomass and the retrieved biomass was 0.87. The results indicate the potential for biomass estimation in wetland environments by use of a combination of PolSAR data, backscatter models, and ANNs.

A general problem with ANNs is the optimization, which influences the inversion accuracy, hence more work remains to be done. As far as the backscatter model is concerned, the soil status should be taken into consideration in subsequent work.


The work is supported by the Natural Science Foundation of China (Grant No. 41401483), the Key National Natural Science Foundation of China (Grant No. 61132006), the Director Foundation of CEODE, CAS (Grant No. Y1ZZ05101B) and the Open Fund of State Key Laboratory of Remote Sensing Science (Grant No. OFSLRSS201205). The authors thank PolSARPro software for the open source programs and the anonymous reviewers for their insightful and helpful comments.



S. Wang, X. Li and Y. Zhou, “Progress of method for wetland vegetation biomass (in Chinese),” Geog. Geo-Inform. Sci. 20(5), 104–109 (2004), http://dx.doi.org/10.3969/j.issn.1672-0504.2004.05.025Google Scholar


G. Berndes, M. Hoogwijk and R. van den Broek, “The contribution of biomass in the future global energy supply: a review of 17 studies,” Biomass Bioenergy. 25(1), 1–28 (2003).BMSBEO0961-9534http://dx.doi.org/10.1016/S0961-9534(02)00185-XGoogle Scholar


D. Lu, “The potential and challenge of remote sensing-based biomass estimation,” Int. J. Remote Sens. 27(7), 1297–1328 (2006).IJSEDK0143-1161http://dx.doi.org/10.1080/01431160500486732Google Scholar


P. Chen et al., “New index for crop canopy fresh biomass estimation (in Chinese),” Spectrosc. Spectral Anal. 30(2), 512–517 (2010).http://dx.doi.org/10.3964/j.issn.1000-0593(2010)02-0512-06Google Scholar


R. Li and J. Liu, “An estimation of wetland vegetation biomass in the Poyang Lake using Landsat ETM data (in Chinese),” Acta Geogr. Sin. 56(5), 532–540 (2001).http://dx.doi.org/10.3321/j.issn:0375-5444.2001.05.004Google Scholar


H. S. Srivastava et al., “Potential applications of multi-parametric synthetic aperture radar (SAR) data in wetland inventory: a case study of Keoladeo National Park (a World Heritage and Ramsar site),” in Proc. of Taal 2007: The 12th World Lake Conf., M. Sengupta and R. Dalwani, Eds., pp. 1862–1879, Bharatpur, India (2008).Google Scholar


F. M. Henderson and A. J. Lewis, “Radar detection of wetland ecosystems: a review,” Int. J. Remote Sens. 29(20), 5809–5835 (2008).IJSEDK0143-1161http://dx.doi.org/10.1080/01431160801958405Google Scholar


E. S. Kasischke, J. M. Melack and M. C. Dobson, “The use of imaging radars for ecological applications—a review,” Remote Sens. Environ. 59(2), 141–156 (1997).RSEEA70034-4257http://dx.doi.org/10.1016/S0034-4257(96)00148-4Google Scholar


N. Baghdadi et al., “Evaluation of C-band SAR data for wetlands mapping,” Int. J. Remote Sens. 22(1), 71–88 (2001).IJSEDK0143-1161http://dx.doi.org/10.1080/014311601750038857Google Scholar


J. Liao, L. Dong and G. Shen, “Neural network algorithm and backscattering model for biomass estimation of wetland vegetation in Poyang Lake area using Envisat ASAR data,” in Proc. IEEE Int. Geosci. Remote Sens. Symp. 2009, Cape Town, Vol. 4, pp. IV-180–IV-183, IEEE (2009), http://dx.doi.org/10.1109/IGARSS.2009.5417344Google Scholar


G. Shen et al., “Retrieval of biomass in Poyang Lake wetland from ENVISAT ASAR data (in Chinese),” Chinese High Tech. Lett. 19(6), 644–649 (2009).http://dx.doi.org/10.3772/j.issn.1002-0470.2009.06.016Google Scholar


J. Liao, G. Shen and L. Dong, “Biomass estimation of wetland vegetation in Poyang Lake area using ENVISAT advanced synthetic aperture radar data,” J. Appl. Remote Sens. 7(1), 073579 (2013).1931-3195http://dx.doi.org/10.1117/1.JRS.7.073579Google Scholar


S. Shen et al., “A scheme for regional rice yield estimation using ENVISAT ASAR data,” Sci. China Ser. D. 52(8), 1183–1194 (2009).SCDEF81006-9313http://dx.doi.org/10.1007/s11430-009-0094-zGoogle Scholar


S. Yang et al., “Rice mapping and monitoring using ENVISAT ASAR data,” IEEE Trans. Geosci. Remote Sens. 5(1), 108–112 (2008).IGRSD20196-2892http://dx.doi.org/10.1109/LGRS.2007.912089Google Scholar


S. Moreau and T. L. Toan, “Biomass quantification of Andean wetland forages using ERS satellite SAR data for optimizing livestock management,” Remote Sens. Environ. 84, 477–492 (2003).RSEEA70034-4257http://dx.doi.org/10.1016/S0034-4257(02)00111-6Google Scholar


R. Touzi, A. Deschamps and G. Rother, “Wetland characterization using polarimetric RADARSAT-2 capability,” Can. J. Remote Sens. 33(sup1), S56–S67 (2007).CJRSDP0703-8992http://dx.doi.org/10.5589/m07-047Google Scholar


B. Brisco et al., “SAR polarimetric change detection for flooded vegetation,” Int. J. Digit. Earth 6(2), 1–12 (2013).1753-8947http://dx.doi.org/10.1080/17538947.2011.608813Google Scholar


K. C. Slatton et al., “Modeling wetland vegetation using polarimetric SAR,” in Proc. IEEE Int. Geosci. Remote Sens. Symp. 1996, Vol. 1, pp. 263–265, IEEE, Lincoln, Nebraska (1996).http://dx.doi.org/10.1109/IGARSS.1996.516311Google Scholar


F. T. Ulaby et al., “Relating the microwave backscattering coefficient to leaf area index,” Remote Sens. Environ. 14(1–3), 113–133 (1984).RSEEA70034-4257http://dx.doi.org/10.1016/0034-4257(84)90010-5Google Scholar


B. Brisco et al., “Water resource applications with RADARSAT-2-a preview,” Int. J. Digit. Earth 1(1), 130–147 (2008).1753-8947http://dx.doi.org/10.1080/17538940701782577Google Scholar


J. Liao and Q. Wang, “Wetland characterization and classification using polarimetric RADARSAT-2 data (in Chinese),” Remote Sens. Land Resour. 20(3), 70–73 (2009).Google Scholar


Q. Wang, J. Liao and B. Tian, “Polarimetric SAR data classification with eigenvalue-based decompositions (in Chinese),” Remote Sens. Technol. Appl. 25(01), 31–37 (2010).Google Scholar


Y. Inoue et al., “Season-long daily measurements of multifrequency (Ka, Ku, X, C, and L) and full-polarization backscatter signatures over paddy rice field and their relationship with biological variables,” Remote Sens. Environ. 81(2–3), 194–204 (2002).RSEEA70034-4257http://dx.doi.org/10.1016/S0034-4257(01)00343-1Google Scholar


S. R. Cloude and E. Pottier, “An entropy based classification scheme for land applications of polarimetric SAR,” IEEE Trans. Geosci. Remote Sens. 35(1), 68–78 (1997).IGRSD20196-2892http://dx.doi.org/10.1109/36.551935Google Scholar


T. Li et al., “A novel vegetation parameters inversion method based on the Freeman decomposition (in Chinese),” J. Electron. Inf. Technol. 33(4), 781–786 (2011).http://dx.doi.org/10.3724/SP.J.1146.2010.00763Google Scholar


F. Del Frate and L. F. Wang, “Sunflower biomass estimation using a scattering model and a neural network algorithm,” Int. J. Remote Sens. 22(7), 1235–1244 (2001).IJSEDK0143-1161http://dx.doi.org/10.1080/01431160151144323Google Scholar


Y. Jin and C. Liu, “Biomass parameter retrieval by using ANN model (in Chinese),” J. Remote Sens. 1(2), 83–87 (1997).Google Scholar


F. Del Frate and D. Solimini, “On neural network algorithms for retrieving forest biomass from SAR data,” IEEE Trans. Geosci. Remote Sens. 42(1), 24–34 (2004).IGRSD20196-2892http://dx.doi.org/10.1109/TGRS.2003.817220Google Scholar


Y. Peng, Y. Jian and R. Li, “Community diversity of aquatic plants in the lakes of Poyang Plain District of China (in Chinese),” J. Central South Forest. Univ. 23(4), 22–27 (2003).http://dx.doi.org/10.3969/j.issn.1673-923X.2003.04.014Google Scholar


Q. Hu et al., “Methane emission from a Carex-dominated wetland in Poyang Lake (in Chinese),” Acta Ecol. Sin. 31(17), 4851–4857 (2011).1872-2032Google Scholar


X. Liu, S. Fan and B. Hu, Comprehensive and Scientific Survey of Jiangxi Nanjishan Wetland Nature Reserve (in Chinese), China Forestry Publishing House, Beijing, China (2005).Google Scholar


P. Patel and H. S. Srivastava, “Ground truth planning for synthetic aperture radar (SAR): addressing various challenges using statistical approach,” Int. J. Adv. Remote Sens. GIS Geogr. 1(2), 1–17 (2013).Google Scholar


G. Hua et al., “Study on the maize backscatter signatures based on polarimetric SAR data (in Chinese),” Jiangsu Agr. Sci. 39(3), 562–565 (2011).http://dx.doi.org/10.3969/j.issn.1002-1302.2011.03.217Google Scholar


Jiangxi Hydrographic Office, “The overview of Jiangxi rainfall and river regime in April,,” The letter of Jiangxi river regime, http://www.jxsw.cn/Item/17931.aspx (2011).Google Scholar


Macdonald, Dettwiler and Associates Ltd. (MDA), RADARSAT-2 Application Look-Up Tables (LUTs), MDA, Richmond, British Colombia (2009).Google Scholar


J. Lee, “Refined filtering of image noise using local statistics,” Comput. Vis. Graph. 15(4), 380–389 (1981).CVGPDB0734-189Xhttp://dx.doi.org/10.1016/S0146-664X(81)80018-4Google Scholar


A. Freeman and S. L. Durden, “A three-component scattering model for polarimetric SAR data,” IEEE Trans. Geosci. Remote Sens. 36(3), 963–973 (1998).IGRSD20196-2892http://dx.doi.org/10.1109/36.673687Google Scholar


Y. Yamaguchi et al., “Four-component scattering model for polarimetric SAR image decomposition,” IEEE Trans. Geosci. Remote Sens. 43(8), 1699–1706 (2005).IGRSD20196-2892http://dx.doi.org/10.1109/TGRS.2005.852084Google Scholar


C. Xie et al., “Analysis of ALOS PALSAR InSAR data for mapping water level changes in Yellow River Delta wetlands,” Int. J. Remote Sens. 34(6), 2047–2056 (2013).IJSEDK0143-1161http://dx.doi.org/10.1080/01431161.2012.731541Google Scholar


L. Ferro-Famil, E. Pottier and J.-S. Lee, “Unsupervised classification of multifrequency and fully polarimetric SAR images based on the H/A/Alpha-Wishart classifier,” IEEE Trans. Geosci. Remote Sens. 39(11), 2332–2342 (2001).IGRSD20196-2892http://dx.doi.org/10.1109/36.964969Google Scholar


J.-S. Lee et al., “Unsupervised classification using polarimetric decomposition and the complex Wishart classifier,” IEEE Trans. Geosci. Remote Sens. 37(5), 2249–2258 (1999).IGRSD20196-2892http://dx.doi.org/10.1109/36.789621Google Scholar


C. Wang et al., “Characterizing L-Band scattering of paddy rice in Southeast China with radiative transfer model and multitemporal ALOS/PALSAR imagery,” IEEE Trans. Geosci. Remote Sens. 47(4), 988–998 (2009).IGRSD20196-2892http://dx.doi.org/10.1109/TGRS.2008.2008309Google Scholar


M. A. Karam et al., “A microwave polarimetric scattering model for forest canopies based on vector radiative transfer theory,” Remote Sens. Environ. 53(1), 16–30 (1995).RSEEA70034-4257http://dx.doi.org/10.1016/0034-4257(95)00048-6Google Scholar


Y. Zhang, Acreage Extraction and Biomass Estimation of Paddy Rice Based on Microwave Remote Sensing and Methane Emissions Simulation from Paddy Fields, Zhejiang University, Hangzhou, China (2009).Google Scholar


J. Liao et al., “Study of microwave dielectric properties of rice growth stages (in Chinese),” J. Remote Sens. 6(1), 19–23 (2002).http://dx.doi.org/10.3321/j.issn:1007-4619.2002.01.004Google Scholar


H. S. Srivastava et al., “A semi-empirical modeling approach to calculate two-way attenuation in radar backscatter from soil due to crop cover,” Curr. Sci. 100(12), 1871–1874 (2011).CUSCAM0011-3891Google Scholar


P. Patel, H. S. Srivastava and R. R. Navalgund, “Use of synthetic aperture radar polarimetry to characterize wetland targets of Keoladeo National Park, Bharatpur, India,” Curr. Sci. 97(4), 529–537 (2009).CUSCAM0011-3891Google Scholar


Y. Wang et al., “Understanding the radar backscattering from flooded and nonflooded Amazonian forests: results from canopy backscatter modeling,” Remote Sens. Environ. 54(3), 324–332 (1995).RSEEA70034-4257http://dx.doi.org/10.1016/0034-4257(95)00140-9Google Scholar


X. Luo, Artificial Neural Network Theory Model Algorithm and Application (in Chinese), Guangxi Normal University Press, Guilin, Guangxi, China (2005).Google Scholar


P. J. van Oevelen and D. H. Hoekman, “Radar backscatter inversion techniques for estimation of surface soil moisture: EFEDA-Spain and HAPEX-Sahel case studies,” IEEE Trans. Geosci. Remote Sens. 37(1), 113–123 (1999).IGRSD20196-2892http://dx.doi.org/10.1109/36.739141Google Scholar


P. Patel, H. S. Srivastava and R. R. Navalgund, “Estimating wheat yield: an approach for estimating number of grains using cross-polarised ENVISAT-1 ASAR data,” Proc. SPIE 6410, 641009 (2006).PSISDG0277-786Xhttp://dx.doi.org/10.1117/12.693930Google Scholar


H. S. Srivastava et al., “Multi-frequency and multi-polarized SAR response to thin vegetation and scattered trees,” Curr. Sci. 97(3), 425–429 (2009).CUSCAM0011-3891Google Scholar


H. S. Srivastava, P. Patel and R. R. Navalgund, “Application potentials of synthetic aperture radar interferometry for land-cover mapping and crop-height estimation,” Curr. Sci. 91(6), 783–788 (2006).CUSCAM0011-3891Google Scholar


J. Liu, J. Liao and G. Shen, “Retrieval of wetland vegetation biomass in Poyang Lake based on quad-polarization image (in Chinese),” Remote Sens. Land Resour. 94(3), 38–43 (2012).http://dx.doi.org/10.6046/gtzyyg.2012.03.08Google Scholar


R. J. Daoust and D. L. Childers, “Quantifying aboveground biomass and estimating net aboveground primary production for wetland macrophytes using a non-destructive phenometric technique,” Aquat. Bot. 62(2), 115–133 (1998).AQBODShttp://dx.doi.org/10.1016/S0304-3770(98)00078-3Google Scholar


Guozhuang Shen received his BS degree from Zhejiang University in 2003 and his PhD degree from the Institute of Remote Sensing Applications, Chinese Academy of Sciences, in 2008. Since 2008, he has been working at the Institute of Remote Sensing and Digital Earth, Chinese Academy of Sciences. Now he mainly does research on flood inundation information extraction, lakes’ variations responding to the global change, and the wetland ecosystem.

Jingjuan Liao received her BS and MS degrees from Nanjing University in 1987 and 1990, respectively, and her PhD degree from the Institute of Geophysics, Chinese Academy of Sciences, in 1993. Since 1993, she has been working on radar remote sensing applications at the Institute of Remote Sensing and Digital Earth, Chinese Academy of Sciences. Her current research interests include microwave remote sensing for surface parameters estimation and the integration of remote sensing observations in wetlands and lakes.

Huadong Guo graduated from the Geology Department at Nanjing University in 1977 and received his MSc degree from the Graduate University of the Chinese Academy of Science (CAS) in 1981. He is a guest professor at eight universities in China and has written more than 200 publications, including 16 books. His current research includes radar remote sensing, applications of Earth observing technologies to global change, and Digital Earth.

Ju Liu received her MS degree from Center for Earth Observation and Digital Earth, Chinese Academy of Sciences, in 2012. Currently, she works at Construction Information Center of Hebei Province.

Guozhuang Shen, Jingjuan Liao, Huadong Guo, Ju Liu, "Poyang Lake wetland vegetation biomass inversion using polarimetric RADARSAT-2 synthetic aperture radar data," Journal of Applied Remote Sensing 9(1), 096077 (4 March 2015). http://dx.doi.org/10.1117/1.JRS.9.096077




Synthetic aperture radar

Data modeling



Back to Top