Relationships between polarimetric SAR backscatter and forest canopy and sub-canopy biophysical properties

This paper highlights the potential of multiwaveband polarimetric SAR data for the estimation of both canopy (percentage canopy closure) and sub-canopy (stem biomass) biophysical variables of a Sitka spruce forest in upland Wales. Stand stem biomass was estimated using forest survey data on diameter at breast height (DBH) and tree height from 0.01 ha plots. Photographs of the forest canopy were taken using a camera fitted with a wide-angle fisheye lens from a number of locations within a stand. The photographs were later digitized and estimates of stand percentage canopy closure were derived using image processing software. It was found that C-band HV and VV, and L-band HV and VV polarization backscatter were significantly related to stem biomass. There was no sensitivity to percentage canopy closure using single polarization backscatter but highly significant relationships were obtained using ratios of single polarization backscatter and variables derived from the polarization signatures. The strong correlations between C-band backscatter and stem biomass indicated a relationship between the structure of the top crown layer and sub-canopy biomass.     

Integration of radar and Landsat-derived foliage projected cover for woody regrowth mapping, Queensland, Australia

In Queensland, Australia, forest areas are discriminated from non-forest by applying a threshold (∼ 12%) to Landsat-derived Foliage Projected Cover (FPC) layers (equating to ∼ 20% canopy cover), which are produced routinely for the State. However, separation of woody regrowth following agricultural clearing cannot be undertaken with confidence, and is therefore not mapped routinely by State Agencies. Using fully polarimetric C-, L- and P-band NASA AIRSAR and Landsat FPC data for forests and agricultural land near Injune, central Queensland, we corroborate that woody regrowth dominated by Brigalow (Acacia harpophylla) cannot be discriminated using either FPC or indeed C-band data alone, because the rapid attainment of a canopy cover leads to similarities in both reflectance and backscatter with remnant forest. We also show that regrowth cannot be discriminated from non-forest areas using either L-band or P-band data alone. However, mapping can be achieved by thresholding and intersecting these layers, as regrowth is unique in supporting both a high FPC (> ∼ 12%) and C-band SAR backscatter (> ~ − 18 dB at HV polarisation) and low L-band and P-band SAR backscatter (e.g. < =∼ 14 dB at L-band HH polarisation). To provide a theoretical explanation, a wave scattering model based on that of Durden et al. [Durden, S.L., Van Zyl, J.J. & Zebker, H.A. (1989). Modelling and observation of radar polarization signature of forested areas. IEEE Trans. Geoscience and Remote Sensing, 27, 290-301.] was used to demonstrate that volume scattering from leaves and small branches in the upper canopy leads to increases in C-band backscattering (particularly HV polarisations) from regrowth, which increases proportionally with FPC. By contrast, low L-band and P-band backscatter occurs because of the lack of double bounce interactions at co-polarisations (particularly HH) and volume scattering at HV polarisation from the stems and branches, respectively, when their dimensions are smaller than the wavelength. Regrowth maps generated by applying simple thresholds to both FPC and AIRSAR L-band data showed a very close correspondence with those mapped using same-date 2.5 m Hymap data and an average 73.7% overlap with those mapped through time-series comparison of Landsat-derived land cover classifications. Regrowth mapped using Landsat-derived FPC from 1995 and JER-1 SAR data from 1994-1995 also corresponded with areas identified within the time-series classification and true colour stereo photographs for the same period. The integration of Landsat FPC and L-band SAR data is therefore expected to facilitate regrowth mapping across Queensland and other regions of Australia, particularly as Japan's Advanced Land Observing System (ALOS) Phase Arrayed L-band SAR (PALSAR), to be launched in 2006, will observe at both L-band HH and HV polarisations.     

Decomposition of polarimetric synthetic aperture radar backscatter from upland and flooded forests

The goal of this research was to decompose polarimetric Synthetic Aperture Radar (SAR) imagery of upland and flooded forests into three backscatter types: single reflection, double reflection, and cross-polarized backscatter. We used a decomposition method that exploits the covariance matrix of backscatter terms. First we applied this method to SAR imagery of dihedral and trihedral corner reflectors positioned on a smooth, dry lake bed, and verified that it accurately isolated the different backscatter types. We then applied the method to decompose multi-frequency Jet Propulsion Laboratory (JPL) airborne SAR (AIRSAR) backscatter from upland and flooded forests to explain scattering components in SAR imagery from forested surfaces. For upland ponderosa pine forest in California, as SAR wavelength increased from C-band to P-band, scattering with an odd number of reflections decreased and scattering with an even number of reflections increased. There was no obvious trend with wavelength for cross-polarized scattering. For a bald cypress-tupelo floodplain forest in Georgia, scattering with an odd number of reflections dominated at C-band. Scattering power with an even number of reflections from the flooded forest was strong at L-band and strongest at P-band. Cross-polarized scattering may not be a major component of total backscatter at all three wavelengths. Various forest structural classes and land cover types were readily distinguishable in the imagery derived by the decomposition method. More importantly, the decomposition method provided a means of unraveling complex interactions between radar signals and vegetated surfaces in terms of scattering mechanisms from targets. The decomposed scattering components were additions to the traditional HH and V V backscatter. One cautionary note: the method was not well suited to targets with low backscatter and a low signal-to-noise ratio.     

Impact of spatial variability of tropical forest structure on radar estimation of aboveground biomass

Understanding the spatial variability of tropical forest structure and its impact on the radar estimation of aboveground biomass (AGB) is important to assess the scale and accuracy of mapping AGB with future low frequency radar missions. We used forest inventory plots in old growth, secondary succession, and forest plantations at the La Selva Biological Station in Costa Rica to examine the spatial variability of AGB and its impact on the L-band and P-band polarimetric radar estimation of AGB at multiple spatial scales. Field estimation of AGB was determined from tree size measurements and an allometric equation developed for tropical wet forests. The field data showed very high spatial variability of forest structure with no spatial dependence at a scale above 11 m in old-growth forest. Plot sizes of greater than 0.25 ha reduced the coefficients of variation in AGB to below 20% and yielded a stationary and normal distribution of AGB over the landscape. Radar backscatter measurements at all polarization channels were strongly positively correlated with AGB at three scales of 0.25 ha, 0.5 ha, and 1.0 ha. Among these measurements, PHV and LHV showed strong sensitivity to AGB < 300 Mg ha^− 1 and AGB < 150 Mg ha^− 1 respectively at the 1.0 ha scale. The sensitivity varied across forest types because of differences in the effects of forest canopy and gap structure on radar attenuation and scattering. Spatial variability of structure and speckle noise in radar measurements contributed equally to degrading the sensitivity of the radar measurements to AGB at spatial scales less than 1.0 ha. By using algorithms based on polarized radar backscatter, we estimated AGB with RMSE = 22.6 Mg ha^− 1 for AGB < 300 Mg ha^− 1 at P-band and RMSE = 23.8 Mg ha^− 1 for AGB < 150 Mg ha^− 1 at L-band and with the accuracy optimized at 1-ha scale within 95% confidence interval. By adding the forest height, estimated from the C-band Interferometry data as an independent variable to the algorithm, the AGB estimation improved beyond the backscatter sensitivity by 20% at P-band and 40% at L-band. The results suggested the estimation of AGB can be improved substantially from the fusion of lidar or InSAR derived forest height with the polarimetric backscatter.     

JERS-1/SAR backscatter and its relationship with biomass of regenerating forests

It has been hypothesized that regenerating tropical forests are large atmospheric carbon sinks. Accurate estimates of the location, extent and biomass of regenerating tropical forests are needed in order to quantify their contribution to global carbon budgets. Synthetic Aperture Radar (SAR) data are independent of near-constant tropical cloud cover and have proved useful for locating and mapping the extent of regenerating tropical forests. To estimate the biomass of regenerating tropical forests we need to determine the nature and strength of the relationship between radar backscatter and biomass for different types of regenerating forest. To further investigate this, two extreme forms of regenerating forest were considered; they were block-logged (clear-cut) forest in the Tapajós area of Pará State, Brazil and selectively-logged forest in Southern Cameroon. Biomass was estimated alometrically for 15 plots in Tapajós and 34 plots in Cameroon and was related to L-band backscatter derived from the JERS-1 SAR. The relationship between backscatter and biomass was strong for the Tapajós study area and weak for the Cameroonian study area. It was concluded that there is potential for the use of JERS-1/SAR to locate, map and estimate biomass for young regenerating forests following block-logging rather than selective-logging.     

Airborne P-band SAR applied to the aboveground biomass studies in the Brazilian tropical rainforest

The objective of this study is to evaluate the relationship between the response (σ°) of airborne P-band Synthetic Aperture Radar (SAR) polarimetric data versus biomass values of primary forest and secondary succession. To ensure that different landscapes of “Terra firme” tropical forest of the Brazilian Amazon were represented, a test-site was selected in the lower Tapajós river region (Pará State). The microwave signals from the P-band polarimetric images were related to the aboveground biomass data by statistical regression models (logarithmic and polynomial functions). In the field survey, physiognomic and structural aspects of primary forest and regrowth were collected and afterwards the biomass was estimated using allometric equations based on dendrometric parameters. As an example of the potential use of P-band polarimetric images, they were classified by a contextual classifier (ICM), whose thematic stratification of land use/land cover was associated with biomass class intervals for mapping purposes. The main objective of this P-band experiment is to improve this tool for regional mapping of Amazon landscape changes, due to the growing rate of land use occupation.     

Soil moisture limitations on monitoring boreal forest regrowth using spaceborne L-band SAR data

A study was carried out to investigate the utility of L-band SAR data for estimating aboveground biomass in sites with low levels of vegetation regrowth. Data to estimate biomass were collected from 59 sites located in fire-disturbed black spruce forests in interior Alaska. PALSAR L-band data (HH and HV polarizations) collected on two dates in the summer/fall of 2007 and one date in the summer of 2009 were used. Significant linear correlations were found between the log of aboveground biomass (range of 0.02 to 22.2 t ha^-1) and σ° (L-HH) and σ° (L-HV) for the data collected on each of the three dates, with the highest correlation found using the L-HV data collected when soil moisture was highest. Soil moisture, however, did change the correlations between L-band σ° and aboveground biomass, and the analyses suggest that the influence of soil moisture is biomass dependent. The results indicate that to use L-band SAR data for mapping aboveground biomass and monitoring forest regrowth will require development of approaches to account for the influence that variations in soil moisture have on L-band microwave backscatter, which can be particularly strong when low levels of aboveground biomass occur.     

The BIOMASS mission: Mapping global forest biomass to better understand the terrestrial carbon cycle

In response to the urgent need for improved mapping of global biomass and the lack of any current space systems capable of addressing this need, the BIOMASS mission was proposed to the European Space Agency for the third cycle of Earth Explorer Core missions and was selected for Feasibility Study (Phase A) in March 2009. The objectives of the mission are 1) to quantify the magnitude and distribution of forest biomass globally to improve resource assessment, carbon accounting and carbon models, and 2) to monitor and quantify changes in terrestrial forest biomass globally, on an annual basis or better, leading to improved estimates of terrestrial carbon sources (primarily from deforestation); and terrestrial carbon sinks due to forest regrowth and afforestation. These science objectives require the mission to measure above-ground forest biomass from 70° N to 56° S at spatial scale of 100-200 m, with error not exceeding ± 20% or ± 10 t ha^− 1 and forest height with error of ± 4 m. To meet the measurement requirements, the mission will carry a P-Band polarimetric SAR (centre frequency 435 MHz with 6 MHz bandwidth) with interferometric capability, operating in a dawn-dusk orbit with a constant incidence angle (in the range of 25°-35°) and a 25-45 day repeat cycle. During its 5-year lifetime, the mission will be capable of providing both direct measurements of biomass derived from intensity data and measurements of forest height derived from polarimetric interferometry. The design of the BIOMASS mission spins together two main observational strands: (1) the long heritage of airborne observations in tropical, temperate and boreal forest that have demonstrated the capabilities of P-band SAR for measuring forest biomass; (2) new developments in recovery of forest structure including forest height from Pol-InSAR, and, crucially, the resistance of P-band to temporal decorrelation, which makes this frequency uniquely suitable for biomass measurements with a single repeat-pass satellite. These two complementary measurement approaches are combined in the single BIOMASS sensor, and have the satisfying property that increasing biomass reduces the sensitivity of the former approach while increasing the sensitivity of the latter. This paper surveys the body of evidence built up over the last decade, from a wide range of airborne experiments, which illustrates the ability of such a sensor to provide the required measurements.At present, the BIOMASS P-band radar appears to be the only sensor capable of providing the necessary global knowledge about the world's forest biomass and its changes. In addition, this first chance to explore the Earth's environment with a long wavelength satellite SAR is expected to make yield new information in a range of geoscience areas, including subsurface structure in arid lands and polar ice, and forest inundation dynamics.     

Sensitivity of space-borne SAR data to forest parameters over sloping terrain. Theory and experiment

Recently, SAR data proved to be useful for the retrieval of forest biomass. However, the effects of terrain slope must be addressed towards the generalization of biomass retrieval for varied forest and environmental conditions. To this aim, we developed experimental and theoretical approaches allowing the study of multi-frequency/multi-polarization forest backscatter of a given forest type, as a function of forest parameters and SAR local incidence angle over the relief. The experimental results showed that the sensitivity of SAR data to biomass was similar to that obtained over a flat terrain, only if the backscatter data were calibrated for slope effects. Moreover, the backscatter must also be corrected for its angular decrease, which can be removed using a simple angular model developed under assumptions of theoretical equations. The highest correlation of corrected backscatter with forest parameters related to aboveground biomass (such as stand age and bole volume) was achieved at L-HV 55° (R ²     

Exploring LiDAR-RaDAR synergy—predicting aboveground biomass in a southwestern ponderosa pine forest using LiDAR, SAR and InSAR

Scanning Light Detecting and Ranging (LiDAR), Synthetic Aperture Radar (SAR) and Interferometric SAR (InSAR) were analyzed to determine (1) which of the three sensor systems most accurately predicted forest biomass, and (2) if LiDAR and SAR/InSAR data sets, jointly considered, produced more accurate, precise results relative to those same data sets considered separately. LiDAR ranging measurements, VHF-SAR cross-sectional returns, and X- and P-band cross-sectional returns and interferometric ranges were regressed with ground-estimated (from dbh) forest biomass in ponderosa pine forests in the southwestern United States. All models were cross-validated. Results indicated that the average canopy height measured by the scanning LiDAR produced the best predictive equation. The simple linear LiDAR equation explained 83% of the biomass variability (n = 52 plots) with a cross-validated root mean square error of 26.0 t/ha. Additional LiDAR metrics were not significant to the model. The GeoSAR P-band (λ = 86 cm) cross-sectional return and the GeoSAR/InSAR canopy height (X-P) captured 30% of the forest biomass variation with an average predictive error of 52.5 t/ha. A second RaDAR-FOPEN collected VHF (λ ∼ 7.8 m) and cross-polarized P-band (λ = 88 cm) cross-sectional returns, none of which proved useful for forest biomass estimation (cross-validated R² = 0.09, RMSE = 63.7 t/ha). Joint consideration of LiDAR and RaDAR measurements produced a statistically significant, albeit small improvement in biomass estimation precision. The cross-validated R² increased from 83% to 84% and the prediction error decreased from 26.0 t/ha to 24.9 t/ha when the GeoSAR X-P interferometric height is considered along with the average LiDAR canopy height. Inclusion of a third LiDAR metric, the 60th decile height, further increased the R² to 85% and decreased the RMSE to 24.1 t/ha. On this 11 km² ponderosa pine study area, LiDAR data proved most useful for predicting forest biomass. RaDAR ranging measurements did not improve the LiDAR estimates.     

On the relationship between the circumference of rubber trees and L-band waves

Despite substantial research conducted within the forestry domain, detailed assessments to monitor plantations and support their sustainable management have been understudied. This article attempts to fill this gap through coupling fully polarimetric L-band data and contemporary data mining methods for the estimation of tree circumference as: (1) a primary dataset for biomass accumulation studies; and, (2) critical information for operational management in rubber plantations. We used two rubber plantation sites in Subang (West Java) and Jember (East Java), Indonesia, to evaluate the capability of L-band radar data. Although polarimetric features derived from polarimetric decomposition theorems have been advocated by others, we show that backscatter coefficients, especially HV polarization, remain an important dataset for this research domain. Using Subang data to build the model, we found that modern machine learning methods do not always deliver the best performance. It appears that the data being ingested plays a significant role in obtaining a good model, hence careful selection of datasets from multiple forms of polarimetric SAR data needs to be further considered. The highest coefficient of determination (R² = 0.79) was achieved by Yamaguchi decomposition features with the aid of partial least squares regression. Nonetheless, we note that the R² gap was insignificant to the backscatter coefficient when random forests regression was used (R² = 0.78). Overall, only the backscatter coefficient dataset delivered fairly consistent results with any regression model, with the average R² being about 0.67. When tuning parameters were not assessed, random forests consistently outweighed support vector regressions in all forms of datasets. The latter generated a substantial increase in R² when a linear kernel was used instead of the popular radial basis function. The issue of transferability of the model is also addressed in this article. It appears that similarity of terrain characteristics substantially influences the model’s performance. Models developed in Subang, which has gentle slopes, seem valid only in plantations with similar terrain. Validation attempts in very flat terrain within two plantation sectors in Jember delivered a poor result, although they have similar elevations to the Subang site. In contrast, validation in a plantation sector with similar, gently sloping terrain achieved an R² of about 0.6 using some datasets.     

Estimation of tree growth in a conifer plantation over 19 years from multi-satellite L-band SAR

The general capability of synthetic aperture radar (SAR) for monitoring forest ecosystems is well documented. However, the majority of SAR studies of forest dynamics use only imagery acquired by one SAR system and are thus limited to the lifecycle of a particular satellite. The synergistic analysis of SAR data from one of the earliest spaceborne SAR missions, the SEASAT mission, with the Japanese JERS-1 satellite-borne SAR is presented. Biophysical parameters frequently retrieved from SAR are tree biomass using backscatter and tree height from the interferometric phase. One potential application that has not been thoroughly examined is mapping of incremental tree growth from SAR backscatter changes. Tree growth measures biomass changes over time, and is correlated to the amount of carbon sequestered by the trees. This paper examines the retrieval of tree growth from multitemporal spaceborne L-band SAR. A SEASAT SAR image from 1978 and a JERS-1 SAR image from 1997 over Thetford forest, UK are used to retrieve tree growth of Corsican Pine stands. Incremental growth was estimated from the changes in backscatter coefficient, and compared to the expected tree growth from general yield class models used by the UK Forestry Commission. The accuracy of the retrieval algorithm depends on the minimum forest stand size included in the analysis. For managed forest plantations, multitemporal L-band SAR has some potential for detecting incremental biomass to support sustainable forest management.     

Forest canopy height and carbon estimation at Monks Wood National Nature Reserve, UK, using dual-wavelength SAR interferometry

Forest canopy height is a critical parameter in better quantifying the terrestrial carbon cycle. It can be used to estimate aboveground biomass and carbon pools stored in the vegetation, and predict timber yield for forest management. Polarimetric SAR interferometry (PolInSAR) uses polarimetric separation of scattering phase centers derived from interferometry to estimate canopy height. A limitation of PolInSAR is that it relies on sufficient scattering phase center separation at each pixel to be able to derive accurate forest canopy height estimates. The effect of wavelength-dependent penetration depth into the canopy is known to be strong, and could potentially lead to a better height separation than relying on polarization combinations at one wavelength alone. Here we present a new method for canopy height mapping using dual-wavelength SAR interferometry (InSAR) at X- and L-band. The method is based on the scattering phase center separation at different wavelengths. It involves the generation of a smoothed interpolated terrain elevation model underneath the forest canopy from repeat-pass L-band InSAR data. The terrain model is then used to remove the terrain component from the single-pass X-band interferometric surface height to estimate forest canopy height. The ability of L-band to map terrain height under vegetation relies on sufficient spatial heterogeneity of the density of scattering elements that scatter L-band electromagnetic waves within each resolution cell. The method is demonstrated with airborne X-band VV polarized single-pass and L-band HH polarized repeat-pass SAR interferometry using data acquired by the E-SAR sensor over Monks Wood National Nature Reserve, UK. This is one of the first radar studies of a semi-natural deciduous woodland that exhibits considerable spatial heterogeneity of vegetation type and density. The canopy height model is validated using airborne imaging LIDAR data acquired by the Environment Agency. The rmse of the LIDAR canopy height estimates compared to theodolite data is 2.15 m (relative error 17.6%). The rmse of the dual-wavelength InSAR-derived canopy height model compared to LIDAR is 3.49 m (relative error 28.5%). From the canopy height maps carbon pools are estimated using allometric equations. The results are compared to a field survey of carbon pools and rmse values are presented. The dual-wavelength InSAR method could potentially be delivered from a spaceborne constellation similar to the TerraSAR system.     

Empirical relationships between AIRSAR backscatter and LiDAR-derived forest biomass, Queensland, Australia

To evaluate the use of multi-frequency, polarimetric Synthetic Aperture Radar (SAR) data for quantifying the above ground biomass (AGB) of open forests and woodlands, NASA JPL AIRSAR (POLSAR) data were acquired over a 37 × 60 km area west of Injune, central Queensland, Australia. From field measurements recorded within 32 50 × 50 m plots, AGB was estimated by applying species-specific allometric equations to stand measurements. AGB was then scaled-up to the larger area using relationships established with Light Detection and Ranging (LiDAR) data acquired over 150 (10 columns, 15 rows) 500 × 150 m cells (or Primary Sampling Units, PSUs) spaced 4 × 4 km apart in the north- and east-west directions. Large-scale (1 : 4000) stereo aerial photographs were also acquired for each PSU to assess species composition. Based on the LiDAR extrapolations, the median AGB for the PSU grid was 82 Mg ha^− 1 (maximum 164 Mg ha^− 1), with the higher levels associated with forests containing a high proportion of Angophora and Callitris species. Empirical relationships between AGB and SAR backscatter confirmed that C-, L- and P-band saturated at different levels and revealed a greater strength in the relationship at higher incidence angles and a larger dynamic range and consistency of relationships at HV polarizations. A higher level of saturation (above ∼50 Mg ha^− 1) was observed at C-band HV compared to that reported for closed forests which was attributable to a link between foliage projected cover (FPC) and AGB. The study concludes that L-band HV backscatter data acquired at incidence angles approaching or exceeding 45° are best suited for estimating the AGB up to the saturation level of ∼80-85 Mg ha^− 1. For regional mapping of biomass below the level of saturation, the use of the Japanese Space Exploration Agency (JAXA) Advanced Land Observing Satellite (ALOS) Phase Arrayed L-band SAR (PALSAR) is advocated.     

Aboveground biomass retrieval in tropical forests — The potential of combined X- and L-band SAR data use

In the context of reducing emissions from deforestation and forest degradation (REDD) and the international effort to reduce anthropogenic greenhouse gas emissions, a reliable assessment of aboveground forest biomass is a major requirement. Especially in tropical forests which store huge amounts of carbon, a precise quantification of aboveground biomass is of high relevance for REDD activities. This study investigates the potential of X- and L-band SAR data to estimate aboveground biomass (AGB) in intact and degraded tropical forests in Central Kalimantan, Borneo, Indonesia. Based on forest inventory data, aboveground biomass was first estimated using LiDAR data. These results were then used to calibrate SAR backscatter images and to upscale the biomass estimates across large areas and ecosystems. This upscaling approach not only provided aboveground biomass estimates over the whole biomass range from woody regrowth to mature pristine forest but also revealed a spatial variation due to varying growth condition within specific forest types. Single and combined frequencies, as well as mono- and multi-temporal TerraSAR-X and ALOS PALSAR biomass estimation models were analyzed for the development of accurate biomass estimations. Regarding the single frequency analysis overall ALOS PALSAR backscatter is more sensitive to AGB than TerraSAR-X, especially in the higher biomass range (> 100 t/ha). However, ALOS PALSAR results were less accurate in low biomass ranges due to a higher variance. The multi-temporal L- and X-band combined model achieved the best result and was therefore tested for its temporal and spatial transferability. The achieved accuracy for this model using nearly 400 independent validation points was r² = 0.53 with an RMSE of 79 t/ha. The model is valid up to 307 t/ha with an accuracy requirement of 50 t/ha and up to 614 t/ha with an accuracy requirement of 100 t/ha in flat terrain. The results demonstrate that direct biomass measurements based on the synergistic use of L- and X-band SAR can provide large-scale AGB estimations for tropical forests. In the context of REDD monitoring the results can be used for the assessment of the spatial distribution of the biomass, also indicating trends in high biomass ranges and the characterization of the spatial patterns in different forest types.     

Correlation Analysis between L-band SAR and ForestStem Volume in Northeast China

Because Synthetic Aperture Radar(SAR)can penetrate into forest canopy and interact with the primary stem volume contents of the trees (trunk and branch),SAR data are widely used for forest stem volume estimation.This paper investigated the correlation between SAR data and forest stem volume in Xunke,Heilongjiang using the stand-wise forest inventory data in 2003 and ALOS PALSAR data for five dates in 2007.The influences of season and polarizations on the relationship between stem volume and SAR data were studied by analyzing the scatterplots;that was followed by interpretation of the mechanisms primarily based on a forest radar backscattering model-water cloud model.The results showed that the relationship between HV polarization backscatter and stem volume is better than HH polarization,and SAR data in summer dry conditions are more correlated to stem volume than the data acquired in other conditions.The interferometric coherence with 46-day temporal baseline is negatively correlated to the stem volume.The correlation coefficients from winter coherence are higher than those from summer coherence and backscatter.The study results suggest using the interferometric coherence in winter as the best choice for forest stem volume estimation with L-band SAR data.     

Sensitivity of SAR data to post-fire forest regrowth in Mediterranean and boreal forests

Disturbed forests may need decades to reach a mature stage and optically-based vegetation indices are usually poorly suited for monitoring purposes due to the rapid saturation of the signal with increasing canopy cover. Spaceborne synthetic aperture radar (SAR) data provide an alternate monitoring approach since the backscattered microwave energy is sensitive to the vegetation structure. Images from two regions in Spain and Alaska were used to analyze SAR metrics (cross-polarized backscatter and co-polarized interferometric coherence) from regrowing forests previously affected by fire. TerraSAR-X X-band backscatter showed the lowest sensitivity to forest regrowth, with the average backscatter increasing by 1-2 dB between the most recent fire scar and the unburned forest. Increased sensitivity (around 3-4 dB) was observed for C-band Envisat Advanced Synthetic Aperture (ASAR) backscatter. The Advanced Land Observing Satellite (ALOS) Phased Array-type L-band Synthetic Aperture Radar (PALSAR) L-band backscatter presented the highest dynamic range from unburned to recently burned forests (approximately 8 dB). The interferometric coherence showed low sensitivity to forest regrowth at all SAR frequencies. For Mediterranean forests, five phases of forest regrowth were discerned whereas for boreal forest, up to four different regrowth phases could be discerned with L-band SAR data. In comparison, the Normalized Difference Vegetation Index (NDVI) provided reliable differentiation only for the most recent development stages. The results obtained were consistent in both environments.     

Retrieval of growing stock volume in boreal forest using hyper-temporal series of Envisat ASAR ScanSAR backscatter measurements

Methods for the estimation of forest growing stock volume (GSV) are a major topic of investigation in the remote sensing community. The boreal zone contains almost 30% of global forest by area but measurements of forest resources are often outdated. Although past and current spaceborne synthetic aperture radar (SAR) backscatter data are not optimal for forest-related studies, a multi-temporal combination of individual GSV estimates can improve the retrieval as compared to the single-image case. This feature has been included in a novel GSV retrieval approach, hereafter referred to as the BIOMASAR algorithm. One innovative aspect of the algorithm is its independence from in situ measurements for model training. Model parameter estimates are obtained from central tendency statistics of the backscatter measurements for unvegetated and dense forest areas, which can be selected by means of a continuous tree canopy cover product, such as the MODIS Vegetation Continuous Fields product. In this paper, the performance of the algorithm has been evaluated using hyper-temporal series of C-band Envisat Advanced SAR (ASAR) images acquired in ScanSAR mode at 100 m and 1 km pixel size. To assess the robustness of the retrieval approach, study areas in Central Siberia (Russia), Sweden and Québec (Canada) have been considered. The algorithm validation activities demonstrated that the automatic approach implemented in the BIOMASAR algorithm performed similarly to traditional approaches based on in situ data. The retrieved GSV showed no saturation up to 300 m³/ha, which represented almost the entire range of GSV at the study areas. The relative root mean square error (RMSE) was between 34.2% and 48.1% at 1 km pixel size. Larger errors were obtained at 100 m because of local errors in the reference datasets. Averaging GSV estimates over neighboring pixels improved the retrieval statistics substantially. For an aggregation factor of 10 × 10 pixels, the relative RMSE was below 25%, regardless of the original resolution of the SAR data.     

On the influence of canopy structure on the radar backscattering of mangrove forests

This paper is the third of a series which aims to evaluate the effects of canopy structure on the polarimetric radar response of mangrove forests. It complements the experimental and theoretical study of closed canopies presented in the previous papers by analysing two different mangrove stands of equal biomass but which greatly differ in their structure. For the three considered frequencies (C-, L- and P-band), experimental observations show that the back-scattering from the open declining stand is higher than that of the closed forest. The corresponding enhancement factor increases with wavelength and shows maximum values for the HH polarization. The identification of the scattering mechanisms occurring between the incident radar wave and the forest components was performed with the assistance of a polarimetric scattering model based on a radiative transfer approach. For the co-polarizations, results of the simulation study confirm that the backscatter enhancement is mainly due to an increase of either the surface scattering or the interaction component. For the cross-polarization HV at L- and P-bands, the increase of the volume component, originating from a stronger interaction with bigger branches, is found to be responsible for the observed enhancement. These findings confirm the large effect of the canopy structure on the forest backscatter and give rise to two important applications. First, the mapping of open declining mangrove stands appears feasible by using either the backscattering coefficient values, especially at P-HH and P-HV, or the HH-VV phase difference at P-band. Second, the use of the σ °-biomass statistical relationships must be restricted to homogeneous closed canopies.     

Assessment of polarimetric SAR data for discrimination between wet versus dry soil moisture conditions

Radarsat-2 imagery from extreme dry versus wet conditions are compared in an effort to determine the value of using polarimetric synthetic aperture radar (SAR) data for improving estimation of fuel moisture in a chronosequence of Alaskan boreal black spruce ecosystems (recent burns, regenerating forests dominated by shrubs, open canopied forests, moderately dense forest cover). Results show strong distinction between wet and dry conditions for C-HH and C-LR polarized backscatter, and Freeman–Durden and van Zyl surface bounce decomposition parameters (35–65% change for all but the dense spruce site). These four SAR variables have high potential for evaluation of within site surface soil moisture, as well as for relative distinction between wet and dry conditions across sites for lower biomass and sparse canopy forested sites. However, for any given test site except the shrubby regrowth site, van Zyl volume, surface, and double bounce scattering all result in similar percentage increases from dry to wet soil condition. This indicates that for most of these test sites/cases moisture enhances the magnitude of the return for all scattering mechanisms evaluated. Thus, differences in scattering from the interaction of biomass, surface roughness, and moisture condition across sites remains an issue and backscatter due to surface roughness or biomass cannot be uniquely estimated. In contrast, the Cloude–Pottier C-band decomposition variables appear invariant to soil moisture, but may instead account for variations in ecosystem structure and biomass. Further investigation is needed, as results warrant future research focused on evaluation of multiple polarimetric parameters in algorithm development.