Monitoring of rain water storage in forests with satellite radar

The sensitivity of radar backscatter to the amount of intercepted rain in temperate deciduous forests is analyzed to determine the feasibility of retrieval of this parameter from satellite radar data. A backscatter model is validated with X-band radar measurements of a single tree exposed to rain. A good agreement between simulation and measurements is observed and this demonstrates the ability of radar to measure the amount of intercepted rain. The backscatter model is next applied to simulate different satellite radar configurations. To account for forest variability, the backscatter difference between a wet and dry forest canopy is calculated for four deciduous tree species, above a wet and a dry soil. On average, the simulated backscatter of a wet forest canopy is 1 dB higher than the backscatter of a dry forest canopy at co-polarized L-band and 2 dB at co-polarized C and X-band. The simulated sensitivity is in agreement with observations. It is argued that current satellites can retrieve the amount of intercepted rain at best with a reliability of 50%, due to the variability in soil moisture, species composition and system noise. The authors expect that the reliability will improve with the launch of the next generation radar satellites. The results of this analysis may also be used to assess the influence of rain, fog or dew upon other radar applications for temperate deciduous forests     

Backscattering properties of boreal forests at the C- and X-bands

The backscattering properties of boreal forests are studied using empirical airborne and spaceborne radar data from Finland. Airborne measurements were carried out in the summer of 1992 by the HUTSCAT scatterometer at the Teijo test area in southern Finland. The HUTSCAT scatterometer is an eight-channel helicopter-borne profiling radar operating at the C- and X-bands. The ranging capability of the HUTSCAT scatterometer was employed in the semiempirical modeling of forest backscatter. The backscatter profile information was used in the analysis of the canopy transmissivity and the canopy backscattering coefficient by distinguishing backscattering contributions from the canopy and the ground. Additionally, ERS-1 C-band satellite SAR measurements were obtained for the Teijo test area and for the reference test area in Sodankyla in northern Finland. The radar results were compared with operational ground-based forest assessment data on forest compartments (stands) of the area. The key parameter investigated was the stem (bole) volume per hectare. The results obtained show the behavior of the canopy transmissivity and the canopy backscatter as a function of stem volume (directly related to the forest biomass). The influence of seasonal and diurnal changes on, and the effects of the changes in soil moisture to the backscattering coefficient were also investigated     

Seasonal dynamics of C-band backscatter of boreal forests withapplications to biomass and soil moisture estimation

The seasonal changes of the C-band backscattering properties of boreal forests are investigated by applying 1) a semiempirical forest backscattering model and 2) multitemporal ERS-1 SAR data from two test areas in Finland. The semiempirical modeling of forest canopy volume backscattering and extinction properties is based on high-resolution data from the authors' ranging scatterometer HUTSCAT. The response of ERS-1 SAR to forest stem volume (biomass) and other forest characteristics is investigated by employing the National Forest Inventory sample plots, stand-wise forest inventory data and LANDSAT- and SPOT-based digital land use maps. The results show that the correlation between the backscattering coefficient and forest stem volume (biomass) varies from positive to negative depending on canopy and soil moisture. Additionally, the seasonal snow cover and soil freezing/thawing effects cause drastic changes in the radar response. A novel method for the estimation of forest stem volume (biomass) is introduced. This technique is based on the use of: 1) multitemporal ERS-1 SAR data; 2) reference sample plot information; and 3) the semiempirical backscattering model. It is shown that the multitemporal ERS-1 SAR images can be successfully used for estimating the forest stem volume. The effects of soil moisture variations to ERS-1 SAR results have been analyzed using hydrological soil moisture model and in situ data. The results indicate that the semiempirical model can he used for predicting the soil and canopy moisture variations in ERS-1 images     

A three-component scattering model for polarimetric SAR data

An approach has been developed that involves the fit of a combination of three simple scattering mechanisms to polarimetric SAR observations. The mechanisms are canopy scatter from a cloud of randomly oriented dipoles, evenor double-bounce scatter from a pair of orthogonal surfaces with different dielectric constants and Bragg scatter from a moderately rough surface. This composite scattering model is used to describe the polarimetric backscatter from naturally occurring scatterers. The model is shown to describe the behavior of polarimetric backscatter from tropical rain forests quite well by applying it to data from NASA/Jet Propulsion Laboratory's (JPLs) airborne polarimetric synthetic aperture radar (AIRSAR) system. The model fit allows clear discrimination between flooded and nonflooded forest and between forested and deforested areas, for example. The model is also shown to be usable as a predictive tool to estimate the effects of forest inundation and disturbance on the fully polarimetric radar signature. An advantage of this model fit approach is that the scattering contributions from the three basic scattering mechanisms can be estimated for clusters of pixels in polarimetric SAR images. Furthermore, it is shown that the contributions of the three scattering mechanisms to the HH, HV, and VV backscatter can be calculated from the model fit. Finally, this model fit approach is justified as a simplification of more complicated scattering models, which require many inputs to solve the forward scattering problem     

Measurements of the propagation parameters of tree canopies at MMWfrequencies

The presence of trees in a given scene can hamper detection of nearby targets by millimeter-wave (MMW) radars especially at near grazing incidence. Proper characterization of scattering and attenuation in tree canopies is important for optimal detection algorithms. In this paper, a new technique for determining the extinction and volume backscattering coefficients in tree canopies using the measured radar backscatter response is proposed and verified experimentally. The technique, which can be applied to already available wideband radar backscatter data, is used to compute the extinction and volume backscattering coefficients of different tree canopies under various physical conditions. The dynamic range of these coefficients are presented and results at 35 GHz are compared with results at 95 GHz     

Calibration of spaceborne scatterometers using tropical rainforests

Wind scatterometers are radar systems designed specifically to measure the normalized radar backscatter coefficient (σ°) of the ocean's surface in order to determine the near-surface wind vector. Postlaunch calibration of a wind scatterometer can be performed with an extended-area natural target such as the Amazon tropical rain forest. Rain forests exhibit a remarkably high degree of homogeneity in their radar response over a very large area though some spatial and temporal variability exist. The authors present a simple technique for calibrating scatterometer data using tropical rain forests, Using a polynomial model for the rolloff of σ° with incidence angle, the technique determines gain corrections to ensure consistency between different antennas and processing channels. Corrections for the time varying instrument gain are made consistent with a seasonally fixed rain forest response; however, without ground stations or aircraft flights, it is difficult to uniquely distinguish between seasonal variations in the rain forest and slow variations of the system gain. Applying the corrections, the intrinsic variability of the σ° of the rain forest is estimated to be ±0.15 dB, which is the limit of the accuracy of calibration using the rain forest. The technique is illustrated with Seasat scatterometer (SASS) data and applied to ERS-1 Active Microwave Instrument scatterometer (Escat) data. Gain corrections of up to several tenths of a decibel are estimated for SASS. Corrections for Escat data are found to be very small, suggesting that Escat data is well calibrated     

Radar estimates of aboveground biomass in boreal forests ofinterior Alaska

Airborne SAR data gathered by the NASA/JPL three-frequency, polarimetric, radar system in winter, spring, and summer over the Bonanza Creek Experimental Forest, near Fairbanks, AK, are compared to estimates of whole-tree aboveground dry biomass from 21 forest stands and two clear-cuts. While C-band radar backscatter shows little sensitivity to biomass, L- and P-band radar backscatter increase by more than 6 dB when biomass increases from 5 to 200 tons/ha. Using second-order polynomial regressions, biomass values are predicted from the radar at L- and P-band and compared to actual biomass values. At P-band HV-polarization, the error in predicted biomass is about 30% of the actual biomass. When HV-, HH-, and VV-polarization are used together in the regression, the error in predicted biomass is about 20%. Errors obtained using L-band data are a few percents larger. These errors are caused by uncertainties in actual stand biomass estimates, significant inner-stand spatial variations in biomass, unusual conditions of forest stands following natural disturbances, along with interactions of the radar signals with a complex three-dimensional structure of the canopy. Multiple incidence angle data reveal that the incidence angle &thetas; _i of the radar illumination is also a factor influencing the retrieval of biomass, even at HV-polarization, when &thetas;_i>50° or &thetas;_i<25°. Finally, the radar response of the forest-and thereby the regression curves for biomass retrieval-are dependent on the seasonal and environmental conditions     

The use of a helicopter mounted ranging scatterometer forestimation of extinction and scattering properties of forest canopies.II. Experimental results for high-density aspen

For pt.I see ibid., vol.26, no.2, p.140-3 (1988). An analytic expression is derived that describes the backscatter power from a semi-infinite plane parallel homogeneous canopy as a function of distance from an airborne radar. This model is fitted to observed data for a high-density aspen canopy by a modification of a technique developed by A.A. Tyapkin (1960). This inversion of the model provides unbiased estimates of the canopy extinction and backscattering parameters. An active radar calibrator located underneath the canopy provides an independent method of determining the volume extinction coefficient. The results reported indicate that the coefficients change throughout the year. A comparison of these coefficients with H.J. Eom and A.K. Funk's (1984) disk model, using measured canopy properties, shows that at C-band frequency, only a part of the scattering and absorption can be attributed to the canopy leaves     

Multiple Incidence Angle SIR-B Experiment Over Argentina: Mapping of Forest Units

Multiple incidence angle SIR-B data of the Cordón la Grasa region of the Chubut Province of Argentina are used to discriminate nate various forest types by their relative brightness versus incidence angle signatures. The region consists of several species of Nothofagas which change in canopy structure with elevation, slope, and exposure. In general, the factors that appear to impact the radar response most are canopy structure, density, and ground cover (presence or absence of dead trunks and branches in particular). The results of this work indicate that 1) different forest species, and structures of a singlee-species, may be discriminated using multiple incidence angle radar-imagery and 2) it is essential to consider the variation in backscatter due to incidence angle when analyzing and comparing data collected at varying frequencies and polarizations.     

Radiative transfer model for microwave bistatic scattering from forest canopies

A bistatic forest scattering model is developed to simulate scattering coefficients from forest canopies. The model is based on the Michigan Microwave Canopy Scattering (MIMICS) model (hence called Bi-MIMICS) and uses radiative transfer theory, where the first-order fully polarimetric transformation matrix is used. Bistatic radar systems offer advantages over monostatic radar systems because of the additional information provided by the diversity of the geometry. By simulating the forest canopy scattering from multiple viewpoints, we can better understand how the forest scatterers' shape, orientation, density, and permittivity affect the canopy scattering. Bi-MIMICS is parametrized using selected forest stands with different canopy compositions and structure. The simulation results show that bistatic scattering is more sensitive to forest biomass changes than backscattering. Analyzing scattering contributions from different parts of the canopy gives us a better understanding of the microwave's interaction with the tree components. The ground effects can also be studied. Knowledge of the canopy's bistatic scattering behavior combined with additional synthetic aperture radar measurements can be used to improve forest parameter retrievals. The simulation results of the model provide the required information for the design of future bistatic radar systems for forest sensing applications.     

Millimeter wave backscatter from deciduous trees

Measurements of the millimeter-wave backscatter from deciduous (leafed) trees are explained using a simple scattering model. The backscatter from individual leaves was measured in the laboratory and used to justify the use of an average leaf radar cross section when computing the normalized radar cross section (NRCS) of tree canopies. NRCS measurements of canopies show that the direction in which incident radiation impinges on the canopy is an important factor in characterizing radar backscatter. Comparisons of measured NRCS values demonstrate that planophil and erectophil trees can be distinguished based on their backscatter at 215 GHz     

A Physical-Optics Model for Double-Bounce Scattering From Tree Stems Standing on an Undulating Ground Surface

In this paper, a model for prediction of radar backscatter from coniferous forests in the VHF and UHF band is proposed. The model includes the double-bounce scattering originating from vertical stems standing on an undulating ground surface and is based on a physical-optics approach. The model can be used to assess the importance of ground topography in synthetic aperture radar (SAR) imagery of forests, and it is applicable to SAR systems using horizontally transmit and receive polarization (HH). The model was validated against data from the airborne SAR systems CARABAS-II and LORA. Precision measurements of ground topography and forest characterization at a single tree level were used as model input to simulate SAR images. The simulated images were compared to radar data in the frequency bands 22–82 and 225–470 MHz, and it was found that the model could predict much of the variation in backscatter observed in images ($R^{2} = hbox {0.44}$ and 0.65 at best, for the lower and higher frequency band, respectively), which should be compared to $R^{2} approx hbox {0.1}$ if the same model, but assuming a flat ground, was used. The results thus indicate that ground topography must be considered when predicting the variations in backscatter in the SAR images studied. The model did, however, fail to predict the absolute values of the backscattered intensity. The reason for the discrepancy is believed to be the value chosen for stem dielectric constant and unmodeled effects due to wave attenuation, tilting stems, and small-scale surface roughness.      

A three-dimensional radar backscatter model of forest canopies

A three-dimensional forest backscatter model, which takes full account of spatial position of trees in a forest stand is described. A forest stand was divided into cells according to arbitrary spatial resolution. The cells may include “crown”, “trunk”, and “gap” components, determined by the shape, size and position of the trees. The forest floor is represented by a layer of “ground” cells. A ray tracing method was used to calculate backscattering components of 1) direct crown backscatter, 2) direct backscattering from ground, 3) direct backscattering from trunk, 4) crown-ground scattering, and 5) trunk-ground scattering. Both the attenuation and time-delay of microwave signals within cells other than “gap” were also calculated from ray tracing. The backscattering Mueller matrices of these components within the same range intervals were incoherently added to yield the total backscattering of an image pixel. By assuming a zero-mean, multiplicative Gaussian noise for image speckle, the high-resolution images were aggregated to simulate a SAR image with a given spatial resolution and number of independent samples (looks). A well-characterized 150 m×200 m forest stand in Maine, USA, was used to parameterize the model. The simulated radar backscatter coefficients were compared with actual JPL SAR data. The model gives reasonable prediction of backscattering coefficients averaged over the entire stand with agreement between model and data within 1.35 dB for all channels. The correlations between simulated images and SAR data (10 by 15 pixels) were positive and significant at the 0.001 level for all frequencies (P, L, and C bands) and polarizations (HH, HV, and VV)     

Discrete scatter model for microwave radar and radiometer responseto corn: comparison of theory and data

As part of the Multisensor Aircraft Campaign, MACHYDRO, two microwave sensors, NASA's Airborne Synthetic Aperture Radar (AIRSAR) and Pushbroom Microwave Radiometer (PBMR) collected data over the same corn fields during the summer of 1990. During these flights, measurements were made on the ground of soil moisture and plant parameters. In this paper the measured canopy and soil parameters are used in a discrete scatter model to predict the response of both sensors (radar and radiometer). A distorted Born approximation is used to compute the scattering coefficient for the corn canopy. The backscatter coefficient gives the radar response and the radiometer response is obtained by integrating the bistatic coefficient over all scattering angles above ground. The objective of this analysis is to test the model and, in particular, to determine how well a single set of plant parameters and single model can yield agreement with both the radar and radiometer measurements. The model values are in reasonably good agreement with the measurements at horizontal polarization and reflect observed changes in soil moisture     

A radar cross-section model for power lines at millimeter-wave frequencies

The knowledge of radar backscatter characteristics of high-voltage power lines is of great importance in the development of a millimeter-wave wire detection system. In this paper, a very high-frequency technique based on an iterative physical optics approach is developed for predicting polarimetric radar backscattering behavior of power lines of arbitrary strand arrangement. In the proposed scattering model the induced surface current is obtained using the tangent plane approximation in an iterative manner where the first-order current, obtained from the incident wave, is used as the source for the second-order current and so on. The approximation is valid for frequencies where the cable strand diameter is on the order of or larger than the wavelength. It is shown that the copolarized backscatter is dominated by the contribution from the first-order PO currents, whereas the cross-polarized backscatter is generated by the second- and higher order PO currents. Using this model, the effects of radar antenna footprint, surface irregularities, and cable sag (when suspended between towers) on radar backscatter are studied. To verify the validity of the proposed model, theoretical results are compared at 94 GHz with experimental results and are found to be in good agreement.     

Low VHF-band backscatter from coniferous forests on sloping terrain

Low-frequency synthetic aperture radar (SAR) is a promising technique for stem volume retrieval, particularly for dense forests, due to the good penetration of forest canopies. However, it is well known that the dominant scattering mechanism, the trunk-ground dihedral interaction, decreases rapidly on sloping terrain. In this paper, we use low VHF-band SAR data, collected with CARABAS over dense coniferous forests in Sweden, to examine the effect of topography. Using flight passes with different headings, the effect of slope and aspect angle on backscatter is characterized. For tall trees (/spl sim/30 m), on the steepest slopes in the test-site (up to /spl sim/12/spl deg/), differences of up to 8 dB are observed between images acquired with different look directions relative to the slope. A physical model is developed to investigate the different scattering mechanisms and their sensitivity to terrain slopes. The model shows that the trunk-ground scattering still dominates the response for large trees on moderate slopes, and a semiempirical model for the effect of topography on backscatter is proposed. The model shows good agreement with measurements, indicating the possibility of using it to compensate for the effects of sloping terrain when retrieving stem volume in coniferous forest.     

Backward and forward scattering by the melting layer composed ofspheroidal hydrometeors at 5-100 GHz

This paper addresses the behavior of the differential reflectivity, specific attenuation, and specific phase shift due to a melting layer composed of oblate-spheroidal hydrometeors. The results are based on a melting layer model and scattering computations derived from the point-matching technique with the truncation and recurrence adjusted. Computations at 5-100 GHz for five raindrop size distributions at rain rates below 12.5 mm/h are presented. In general, the reflectivity factor and differential reflectivity features with height at centimeter wavelengths agree with available radar measurements. At millimeter wavelengths, contributions to the radar backscatter from smaller hydrometeors become more and more important as the frequency increases and approaches 100 GHz. This should be instructive for utilizing millimeter wavelength radar techniques in radar remote sensing studies of the melting layer. Corresponding vertical profiles of the specific attenuation and phase shift are also presented at 5-100 GHz. The differential attenuation and phase shift indicate the particle shape effects. These attenuation and phase shift become more and more considerable as the frequency increases. Such forward scattering calculations should prove useful for studying propagation effects caused by the melting layer for satellite-earth communications, including depolarizations     

Simulation of interferometric SAR response for characterizing thescattering phase center statistics of forest canopies

A coherent scattering model for tree canopies is employed in order to characterize the sensitivity of an interferometric SAR (INSAR) response to the physical parameters of forest stands. The concept of an equivalent scatterer for a collection of scatterers within a pixel, representing the vegetation particles of tree structures, is used for identifying the scattering phase center of the pixel whose height is measured by an INSAR. Combining the recently developed coherent scattering model for tree canopies and the INSAR Δk-radar-equivalence algorithm, accurate statistics of the scattering phase-center location of forest stands are obtained numerically for the first time. The scattering model is based on a Monte Carlo simulation of scattering from fractal-generated tree structures, and therefore is capable of preserving the absolute phase of the backscatter. The model can also account for coherent effects due to the relative position of individual scatterers and the inhomogeneous extinction experienced by a coherent wave propagating through the random collection of vegetation particles. The location of the scattering phase center and the correlation coefficient are computed using the Δk-radar equivalence simply by simulating the backscatter response at two slightly different frequencies. The model is successfully validated using the measured data acquired by JPL TOPSAR over a selected pine stand in Raco, MI. A sensitivity analysis is performed to characterize the response of coniferous and deciduous forest stands to a multifrequency and multipolarization INSAR in order to determine an optimum system configuration for remote sensing of forest parameters     

Rain backscatter measurements at millimeter wavelengths

The US Army Ballistic Research Laboratory (BRL) conducted an experiment in 1973 to measure the properties of radar backscatter from rain at millimeter wavelengths. Rain backscatter and attenuation were measured with pulse radars operating simultaneously at 9.375, 35, 70, and 95 GHz over a wide range of rain intensities while continuous measurements of raindrop size and rainfall rate were made. This report describes the measurement technique, details of the instrumentation, the data analysis procedure, and the rain backscatter data obtained from A-scope photographs and video tapes