Radiative transfer modeling of cross-polarized backscatter from a pine forest using the discrete ordinate and eigenvalue method

Radiative transfer models have been widely used to interpret the radar backscatter from forested areas. Most of these models are based on an iterative solution of the radiative transfer equation, usually solved up to first or second order, thus taking into account single and double scattering. Although this method leads to results agreeing well with copolarized backscatter measurements, it produces less accurate estimates for horizontal-vertical (HV) polarization. This paper presents a radiative transfer backscatter model that accounts for multiple scattering by using the discrete ordinate and eigenvalue method applied to a layered medium. Using parameters derived from an architectural tree model, calculations at C- and L-band are compared with HV data acquired for a maritime pine forest in the southwest of France during the Spaceborne Imaging Radar-C missions. Good agreement is found at C-band for all values of forest biomass, and reasonable agreement at L-band for high biomass, when the soil backscatter plays a minor role. For low biomass, the L-band modeling is inadequate because of difficulties in estimating the soil backscatter. Comparison with calculations from a first-order radiative transfer model shows that multiple scattering is significant, especially at C-band.     

Microwave scattering from mixed-species forests, Queensland, Australia

The potential of synthetic aperture radar (SAR) data for retrieving the above-ground and component (e.g., branch, trunk) biomass of mixed-species forests (including woodlands) typical to subtropical Queensland, Australia, was evaluated using a wave scattering model based on that of Durden et al. (1989). The model was parameterized using field data collected for nine forest types, which were selected through combined analysis of 1 : 4000 aerial photographs and light detection and ranging data. The simulated SAR backscatter data demonstrated a good correspondence at most frequencies and polarizations with Airborne SAR data. Analysis of scattering mechanisms revealed dominance of C-band horizontal-vertical (HV) volume scattering and increases with small-branch/foliage biomass, dominance of L- and P-band HH trunk-ground scattering and increases with trunk biomass, and dominance of L-band HV volume (branch) scattering and increases with large-branch biomass. The study concluded that above-ground biomass estimated using empirical relationships with selected SAR channels will be more reliable for forests of similar structural form due to dominance of microwave interaction with particular biomass components and the strength and consistency of relationships between these and the affiliated components that represent the total. In mixed-species forests, retrieval will be compromised by interaction with a greater diversity of structures and variability in relationships between structural components. Although empirical relationships with selected combinations of channels (e.g., L-band HH/HV) might allow retrieval of component and total biomass of forests containing trees of similar form (e.g., as mapped using Landsat sensor data), the use of SAR inversion models was considered a more appropriate route for retrieving the biomass of forests containing a mix of structural forms.     

Characterizing L-Band Scattering of Paddy Rice in Southeast China With Radiative Transfer Model and Multitemporal ALOS/PALSAR Imagery

Rice is a major food supply in southeast China. With increased population and urbanization, reliable rice mapping is critical in this region. Because of frequent cloud cover and precipitation during the rice-growing season, it is difficult to conduct large-area rice monitoring with optical remote sensing techniques. L-band synthetic aperture radar (SAR), with its all-weather day and night imaging and canopy penetration capabilities, provides a unique alternative. In this study, a first-order radiative transfer model was developed to simulate L-band scattering properties of paddy rice. Three Advanced Land Observing Satellite (ALOS)/Phased-Array-Type L-band Synthetic Aperture Radar (PALSAR) images in dual-polarization mode (HH and HV) acquired in early tillering (June 28, 2007), tillering (August 13, 2007), and heading (September 28, 2007) stages were processed to test the temporal variation of rice backscatter. It was found that plant height and leaf mass amount were the two major structural parameters that contributed to rice backscatter in PALSAR images. The variation of the simulated HH backscatter matched with PALSAR observations in sample fields, although the simulated backscatter coefficients were around 3 dB lower than image-extracted values. Leaf volume scattering and leaf–ground double bounce were found as the two major scattering components in L-band HH polarization and increased with leaf layer height and density. This paper demonstrated that L-band HH backscatter was more sensitive to rice's structural variation than the VV backscatter and may therefore be more useful in rice mapping and modeling studies.      

Sea ice monitoring by L-band SAR: an assessment based on literature and comparisons of JERS-1 and ERS-1 imagery

Spaceborne single-polarization C-band synthetic aperture radar (SAR) imagery is widely used to gather information about the state of the sea ice cover in the polar regions. C-band is regarded as a reasonable choice for all-season monitoring capabilities. For specific mapping tasks, however, other frequency bands can be more suitable. In the first part of this paper, the summary of a literature study dealing with the utilization of L-band SAR imagery for sea ice monitoring is presented. Investigations reveal that if deformation features such as ice ridges, rubble fields, and brash ice are to be mapped, L-band radar is superior in a number of cases. The second part of this paper addresses the comparison of JERS-1 and ERS-1 SAR images that were acquired over sea ice east of Svalbard and along the east coast of Greenland. The effects of the different frequencies, polarizations, and incidence angles of the two SAR systems are discussed. It is demonstrated that the images of both sensors complement one another in the analysis of ice conditions, resulting in a more detailed view of the sea ice cover state.     

Retrieval of ocean wave spectra and RAR MTF's fromdual-polarization SAR data

A method for the retrieval of the real aperture radar (RAR) modulation transfer function (MTF) and ocean wave spectra from dual-polarization (i.e., simultaneously acquired HH and VV polarizations) synthetic aperture radar (SAR) image data is described. The RAR MTF is estimated by applying empirical MTF estimation methodologies to inter-look cross spectra between various combinations of individual looks and available polarizations for a given radar frequency. The concept behind the nonlinear inversion is that any combination of like- and cross-polarization image spectra should return the same wave spectrum, in agreement with in situ and model wave spectra. This permits estimation of the RAR MTF on a case-by-case basis. The results are compared with theoretical treatments of the RAR MTF, which are shown to be inadequate for the range of conditions encountered in their data set. However, the theory and measurements fit well in describing the polarization dependence of the RAR MTF. The data set consists of SIR-C/X-SAR L-band and CCRS CV-580 C-band SAR data, in situ buoy measurements, and model data from field programs in Canadian waters in October and December 1994     

Multipolarization Radar Images for Geologic Mapping and Vegetation Discrimination

The NASA/JPL airborne synthetic aperture radar system produces radar image data simultaneously in four linear polarizations (HH, VV, VH, HV) at 24.6-cm wavelength (L-band), with 10-m resolution, across a swath width of approximately 10 km. The signal data are recorded optically and digitally and annotated in each of the channels to facilitate a completely automated digital correlation. Both standard amplitude, and also phase difference images are produced in the correlation process. Individual polarization and range-dependent gain functions improve the effective dynamic range, but as yet do not permit absolute quantitative measurements of the scattering coefficients. However, comparison of the relative intensities of the different polarizations in individual black-and-white and color composite images provides discriminatory mapping information. In the Death Valley, California, area, rough surfaces of young alluvial deposits produce strong responses at all polarizations. Smoother surfaces of older alluvial deposits show significantly lower responses. Evaporite deposits of different types and moisture contents have distinct polarization signatures. In the Wind River Basin, Wyoming, sedimentary rock units show polarization responses that relate to differences in weathering. Local intensity variations in like-polarization images result from topographic effects; strong cross-polarization responses denote the effects of vegetation cover and, in some cases, possible scattering from the subsurface. In the Savannah River Plant, South Carolina, forest cover characteristics are discriminated by polarization responses that reflect the density and structure of the canopy, and the presence or absence of standing water beneath the canopy.     

Mapping of forest types in Alaskan boreal forests using SAR imagery

Mapping of forest types in the Tanana river flood-plain, interior Alaska, is performed using a maximum-a-posteriori Bayesian classifier applied on SAR data acquired by the NASA/JPL three-frequency polarimetric AIRSAR system on several dates. Five vegetation types are separated, dominated by 1) white spruce, 2) balsam poplar, 3) black spruce, 4) alder/willow shrubs, and 5) bog/fen/nonforest vegetation. Open water of rivers and lakes is also separated. Accuracy of forest classification is investigated as a function of frequency and polarization of the radar, as well as the forest seasonal state, which includes winter/frozen, winter/thawed, spring/flooded, spring/unflooded, and summer/dry conditions. Classifications indicate that C-band is a more useful frequency for separating forest types than L- or P-bands, and HV polarization is the most useful polarization at all frequencies. The highest classification accuracy, with 90 percent of forest pixels classified correctly, is obtained by combining L-band HV and C-band HV data acquired in spring as seasonal river flooding recedes and before deciduous tree species have leaves. In 17 forest stands for which actual percentages of each tree species are known, the same radar data are capable of predicting tree species composition with less than 10 percent error. From the authors' classification, they predict that current and future spaceborne SAR systems will have limited mapping capabilities when used alone. Yet, RADARSAT combined with J-ERS-1 and ERS-1 could resolve forest types with 80 percent accuracy, separate nonforest areas resulting from commercial logging or forest wildfire, and map river edges     

Dual-band dual-polarized perforated microstrip antennas for SARapplications

For dual-band dual-polarized synthetic aperture radar (SAR) applications a compact low-profile design is investigated. The operating frequencies are in the L and C-bands, centered about 1.275 and 5.3 GHz, respectively. Since the C-band frequency is larger by a factor of four, its array elements and inter-element separations are smaller by the same ratio. Thus, to allow similar scan ranges for both bands, the L-band elements are selected as perforated patches to enable the placement of C-band elements within them. Stacked-patch configurations were used to meet the bandwidth requirements, especially in the L-band. The C-band element was designed numerically, but the perforated L-band one required final experimental optimization. Also, in the latter case of L-band, a balanced transmission line feed was used to minimize cross polarization. For the C-band elements, slot coupling was used and, to simplify the feed, symmetric parasitic slots were incorporated to minimize cross polarization. No vertical connections were utilized, and electromagnetic couplings resulted in a compact low-profile design, with an electrically and thermally symmetric geometry     

Stand age retrieval in production forest stands in New Zealand using C- and L-band polarimetric Radar

Regressions of single-, dual-, quad-, and full-polarization L- and C-band synthetic aperture radar (SAR) against stand age from 403 radiata pine stands in Kaingaroa Forest, New Zealand have been carried out, using the National Aeronautics and Space Administration's Airborne SAR instrument. The regressions attempted to find products suitable for the separation of young (two years or less) from old stands (25 years or older), and for the estimation of stand age. Local incidence angle had no significant effect for C-band, but was always significant for L-band, giving a standard error reduction of 13% to 32% for log stand age. Stand density was highly significant for both bands, giving standard error reductions of 7% to 47%. Single- and dual-polarization products were severely biased, and it was impossible to separate young and old stands, except L-band horizontal-(HH)-plus-horizontal-vertical (HV). C-band quad-polarization gave less bias and lower error than for that L-band, when local incidence angle and stand density were excluded. C-band full-polarization using covariance magnitudes gave no improvement over C-band quad-polarization, but L-band did give a significant improvement. The C- and L-band full-polarization products with six polarimetric indices gave significant improvements in the standard error. The results show that regressions of SAR data with stand age are possible with full- and quad-polarization L- and C-band datasets, although the prediction limits increase rapidly with stand age. The smallest error in estimated stand age, with an RMS of 3.22 years, was for L-band full-polarization with six polarimetric indices, calculated from a validation dataset. Separation of young and old forest stands was only possible for full-, quad-polarization, and the L-band HH-plus-HV products.     

Use of radar and optical remotely sensed data for soil moisture retrieval over vegetated areas

This work assesses the possibility of obtaining soil moisture maps of vegetated fields using information derived from radar and optical images. The sensor and field data were acquired during the SMEX'02 experiment. The retrieval was obtained by using a Bayesian approach, where the key point is the evaluation of probability density functions (pdfs) based on the knowledge of soil parameter measurements and of the corresponding remotely sensing data. The purpose is to determine a useful parameterization of vegetation backscattering effects through suitable pdfs to be later used in the inversion algorithm. The correlation coefficients between measured and extracted soil moisture values are R=0.68 for C-band and R=0.60 for L-band. The pdf parameters have been found to be correlated to the vegetation water content estimated from a Landsat image with correlation coefficients of R=0.65 and 0.91 for C- and L-bands, respectively. In consideration of these correlations, a second run of the Bayesian procedure has been performed where the pdf parameters are variable with vegetation water content. This second procedure allows the improvement of inversion results for the L-band. The results derived from the Bayesian approach have also been compared with a classical inversion method that is based on a linear relationship between soil moisture and the backscattering coefficients for horizontal and vertical polarizations.     

Fully polarimetric bistatic radar scattering behavior of forestedhills

The bistatic radar scattering measurements of forested hills were performed at grazing incidence and at azimuth scattering angles from 28° to 66° from the forward scatter plane. Using pulse-to-pulse switching between orthogonal transmitted polarizations, the radar simultaneously measures two orthogonally polarized components of the scattered wave to obtain full polarimetric information about the scattering process. These are the first fully polarimetric terrain clutter measurements to be conducted at large bistatic angles. The complete Stokes matrix, computed by averaging successive realizations of the polarization scattering matrix, is used to examine the polarization sensitivity of the bistatic clutter. It is found that the polarization state of the EM wave scattered out of the plane of incidence strongly depends on the polarization orientation of the incident electric field. Unlike the monostatic case, these two incident wave polarization states are found to produce substantially different scattered wave behavior when trees are viewed at large bistatic angles. Scattered fields resulting from vertically oriented incident fields are found to be highly polarized and to produce bistatic clutter power levels that are strongly dependent on the polarization of the receiving antenna. In contrast, horizontally oriented incident fields are found to produce weakly polarized scattered waves with bistatic clutter power levels that are insensitive to the polarization of the receiving antenna     

Compact polarimetry based on symmetry properties of geophysical media: the /spl pi//4 mode

We assess the performance of synthetic aperture radar (SAR) compact polarimetry architectures based on mixed basis measurements, where the transmitter polarization is either circular or orientated at 45/spl deg/(/spl pi//4), and the receivers are at horizontal and vertical polarizations with respect to the radar line of sight. An original algorithm is proposed to reconstruct the full polarimetric (FP) information from this architecture. The performance assessment is twofold: it first concerns the level of information preserved in comparison with FP, both for point target analysis and crop fields classification, using L-band SIRC/XSAR images acquired over Landes forest and Jet Propulsion Laboratory AIRSAR images acquired over Flevoland. Then, it addresses the space implementation complexity, in terms of processed swath, downloading features, power budget, calibration, and ionospheric effects. The polarization uniqueness in transmission of this mixed basis mode, hereafter referred to as the /spl pi//4 mode, maintains the standard lower pulse repetition frequency operation and hence maximizes the coverage of the sensor. Because of the mismatch between transmitter and receiver basis, the power budget is deteriorated by a factor of 3 dB, but it can partly be compensated.     

Specular null polarization theory: applications to radarmeteorology

Specular null polarization theory (SNPT) has been recently introduced for the case of coherent scattering where a 2×2 scattering matrix is sufficient to describe the scattering process. In this paper, SNPT is extended to the case of incoherent scattering. Optimum polarization states are derived and the results are discussed in relation to the classic radar optimum polarizations. In traditional radar polarimetry, modeling of the radar receive/transmit network is included in the radar voltage equation and consequently this affects the optimum polarizations and polarization responses of scatterers. SNPT eliminates this effect and therefore allows for a more direct analysis of scatterers. Modeling of ensembles of precipitation particles is used to illustrate the results of the analysis     

Improvement of Moisture Estimation Accuracy of Vegetation-Covered Soil by Combined Active/Passive Microwave Remote Sensing

The measured effects of vegetation canopies on radar and radiometric sensitivity to soil moisture are compared to first-order emission and scattering models. The models are found to predict the measured emission and backscattering with reasonable accuracy for various crop canopies at frequencies between 1.4 and 5.0 GHz, especially at angles of incidence less than 30°. The vegetation loss factor L (?) increases with frequency and is found to be dependent upon canopy type and water content. In addition, the effective radiometric power absorption coefficient of a mature corn canopy is roughly 1.75 times that calculated for the radar at the same frequency. Comparison of an L-band radiometer with a C-band radar shows the two systems to be complementary in terms of accurate soil moisture sensing over the extreme range of naturally occurring soil-moisture conditions. The combination of both an L-band radiometer and a C-band radar is expected to yield soil-moisture estimates that are accurate to better than +/-30 percent of true soil moisture, even for a soil under a lossy crop canopy such as mature corn. This is true even without any other ancillary information.     

Preliminary Report on Measurements of Forest Canopies with C-Band Radar Scatterometer at NASA/NSTL

This paper presents preliminary results of C-band radar scatterometer measurements of forest canopies of southeastern forests in the vicinity of NASA/NSTL. The results are as follows: 1) radar backscattering coefficients (BSC's) of deciduous forests such as oak, maple, blackgum, and cypress are higher than those of coniferous forests such as slash pine plantation and natural pine; 2) at a large incidence angle, where polarization effect is significant, and by ranging measurement, the VV polarization BSC's obtain peak value at the first few meters from the canopy top and decrease rather quickly, while the HH polarization BSC's obtain peak value at longer distances from the canopy top and decrease rather slowly through the canopy; 3) using the active radar calibrator for tree canopy attenuation measurement of a dense and a sparse live oak, it is found that the tree canopies with higher attenuations have higher BSC's for all three polarizations, with VV polarization containing the largest differential (2.2 dB).     

Airborne C-band SAR measurements of wet snow-covered areas

Using polarimetric synthetic aperture radar (SAR) data at C-band, the optimum polarization and range incidence angles for the classification of land covers found in the James Bay area, P.Q., Canada-open areas, lakes ice, and forests, all covered with wet snow-have been examined. The separability between classes shows that they can be classified by a single polarization. For &thetas;<30°, the overall classification accuracies were 97.1, 98, and 90.8% from HH-, cross, and VV-polarizations alone, respectively. They were greater than 99% for all polarizations at &thetas;>30°. However, the polarimetric parameters were not suitable for classification     

Delineation of inundated area and vegetation along the Amazonfloodplain with the SIR-C synthetic aperture radar

Floodplain inundation and vegetation along the Negro and Amazon rivers near Manaus, Brazil were accurately delineated using multi-frequency, polarimetric synthetic aperture radar (SAR) data from the April and October 1994 SIR-C missions. A decision-tree model was used to formulate rules for a supervised classification into five categories: water, clearing (pasture), aquatic macrophyte (floating meadow), nonflooded forest, and flooded forest. Classified images were produced and tested within three days of SIR-C data acquisition. Both C-band (5.7 cm) and L-band (24 cm) wavelengths were necessary to distinguish the cover types. HH polarization was most useful for distinguishing flooded from nonflooded vegetation (C-HH for macrophyte versus pasture, and L-HH for flooded versus nonflooded forest), and cross-polarized L-band data provided the best separation between woody and nonwoody vegetation. Between the April and October missions, the Amazon River level fell about 3.6 m and the portion of the study area covered by flooded forest decreased from 23% to 12%. This study demonstrates the ability of multifrequency SAR to quantify in near realtime the extent of inundation on forested floodplains, and its potential application for timely monitoring of flood events     

Change detection for thematic mapping by means of airbornemultitemporal polarimetric SAR imagery

The paper addresses the detection of changes in multitemporal polarimetric radar images, focusing on small objects and narrow linear features. The images were acquired at C- and L-band by the airborne EMISAR system. It is found that the radar intensities are better suited for change detection than the correlation coefficient and the phase difference between the co-polarized channels. In the case of linear features, there is no obvious difference between the C- and L-bands , and slight variations of the flight tracks are acceptable at look angles larger than 35 degrees. Theoretical detection thresholds are evaluated from the statistical distribution of the intensity ratio due to speckle. For the linear features and for urban environments, the observed thresholds are larger than the theoretical predictions. This is interpreted as an effect of radar intensity variations on length scales smaller than the spatial image resolution. The signature of urban areas is very sensitive to deviations between the flight tracks, and the sensitivity is larger at C-band than at L-band. On the other hand, the intensity contrast between buildings and the urban background is smaller at L-band and larger at C-band. For change detection, thresholds may have to be chosen separately for each object class because the intensity ratios of different object classes vary differently as a function of time     

Absolute radiometric calibration of the CCRS SAR

Determining the radar scattering coefficients from SAR (synthetic aperture radar) image data requires absolute radiometric calibration of the SAR system. The authors describe an internal calibration methodology for the airborne Canada Centre for Remote Sensing (CCRS) SAR system, based on radar theory, a detailed model of the radar system, and measurements of system parameters. The methodology is verified by analyzing external calibration data acquired over a six-month period in 1988 by the C-band radar using HH polarization. The results indicate that the overall error is ±0.8 dB (1σ) for incidence angles ±20° from antenna boresight. The dominant error contributions are due to the antenna radome and uncertainties in the elevation angle relative to the antenna boresight     

Extension of Kennaugh's optimal polarization concept to the asymmetric scattering matrix case

The polarization scattering matrix measured by a bistatic radar system generally will be asymmetric. The scattering matrix will also be asymmetric when the radar system is monostatic, but the intervening propagation medium is anisotropic. Kennaugh's optimal polarizations theory is generalized to the case of the asymmetric scattering matrix. The radar antenna polarizations to be used for maximum and zero power receptions are defined and geometrically interpreted on the Poincare sphere. These polarizations, termed optimal polarizations, may be used to enhance the level of target echo return, to discriminate against undesired interference sources, or to classify radar targets.