Multiple scattering in the high-frequency limit with second-order shadowing function from 2D anisotropic rough dielectric surfaces: I. Theoretical study

Abstract

In this paper the first- and second-order Kirchhoff approximation is applied to study the backscattering enhancement phenomenon, which appears when the surface rms slope is greater than 0.5. The formulation is reduced to the geometric optics approximation in which the second-order illumination function is taken into account. This study is developed for a two-dimensional (2D) anisotropic stationary rough dielectric surface and for any surface slope and height distributions assumed to be statistically even. Using the Weyl representation of the Green function (which introduces an absolute value over the surface elevation in the phase term), the incoherent scattering coefficient under the stationary phase assumption is expressed as the sum of three terms. The incoherent scattering coefficient then requires the numerical computation of a ten- dimensional integral. To reduce the number of numerical integrations, the geometric optics approximation is applied, which assumes that the correlation between two adjacent points is very strong. The model is then proportional to two surface slope probabilities, for which the slopes would specularly reflect the beams in the double scattering process. In addition, the slope distributions are related with each other by a propagating function, which accounts for the second-order illumination function. The companion paper is devoted to the simulation of this model and comparisons with an ‘exact’ numerical method.     

Multiple scattering in the high-frequency limit with second-order shadowing function from 2D anisotropic rough dielectric surfaces: II. comparison with numerical results

This second part presents illustrative examples of the model developed in the companion paper, which is based on the first- and second-order optics approximation. The surface is assumed to be Gaussian and the correlation height is chosen as anisotropic Gaussian. The incoherent scattering coefficient is computed for a height rms range from 0.5λ 1λwhere λ is the electromagnetic wavelength), for a slope rms range from 0.5 to 1 and for an incidence angle range from 0 to 70°. In addition, simulations are presented for an anisotropic Gaussian surface and when the receiver is not located in the plane of incidence. For a metallic and dielectric isotropic Gaussian surfaces, the cross- and co-polarizations are also compared with a numerical approach obtained from the forward.backward method with a novel spectral acceleration algorithm developed by Torrungrueng and Johnson (2001, JOSA A 18). (Some figures in this article are in colour only in the electronic version)     

Multiple scattering in the high-frequency limit with second-order shadowing function from 2D anisotropic rough dielectric surfaces: II. comparison with numerical results

Abstract

This second part presents illustrative examples of the model developed in the companion paper, which is based on the first- and second-order optics approximation. The surface is assumed to be Gaussian and the correlation height is chosen as anisotropic Gaussian. The incoherent scattering coefficient is computed for a height rms range from 0.5λ 1λwhere λ is the electromagnetic wavelength), for a slope rms range from 0.5 to 1 and for an incidence angle range from 0 to 70°. In addition, simulations are presented for an anisotropic Gaussian surface and when the receiver is not located in the plane of incidence. For a metallic and dielectric isotropic Gaussian surfaces, the cross- and co-polarizations are also compared with a numerical approach obtained from the forward.backward method with a novel spectral acceleration algorithm developed by Torrungrueng and Johnson (2001, JOSA A 18). (Some figures in this article are in colour only in the electronic version)     

Bistatic scattering from one-dimensional random rough homogeneous layers in the high-frequency limit with shadowing effect

Many fast asymptotic models of electromagnetic scattering from a single rough interface have been developed over the last few years, but only a few have been developed on stacks of rough interfaces. The specific case of very rough surfaces, compared to the incident wavelength, has not been treated before, which is the context of this paper. The model starts from the iteration of the Kirchhoff approximation to calculate the fields scattered by a rough layer, and is reduced to the high-frequency limit in order to rapidly obtain numerical results. The shadowing effect, important under grazing angles, is taken into account. The model can be applied to any given slope statistics. Then, the model is compared with a reference numerical method based on the method of moments, which validates the model in the high-frequency limit for lossless and lossy inner media.     

Multiple light scattering from metal and dielectric rough surfaces

Abstract

A numerical study is done on light scattering from one-dimensional random rough surfaces supporting both dielectrics and metals. The influence of the corrugation on the Brewster angle as well as on the angular distribution of transmitted intensity into the dielectric is investigated. The authors also obtain enhanced backscattering in the light reflected from metallic surfaces. Finally they discuss the influence of roughness on the drop of reflectance due to polariton absorption.     

Multiple light scattering from metal and dielectric rough surfaces

A numerical study is done on light scattering from one-dimensional random rough surfaces supporting both dielectrics and metals. The influence of the corrugation on the Brewster angle as well as on the angular distribution of transmitted intensity into the dielectric is investigated. The authors also obtain enhanced backscattering in the light reflected from metallic surfaces. Finally they discuss the influence of roughness on the drop of reflectance due to polariton absorption.     

Energy conservation of the scattering from one-dimensional random rough surfaces in the high-frequency limit

Energy conservation of the scattering from one-dimensional strongly rough dielectric surfaces is investigated using the Kirchhoff approximation with single reflection and by taking the shadowing phenomenon into account, both in reflection and transmission. In addition, because no shadowing function in transmission exists in the literature, this function is presented here in detail. The model is reduced to the high-frequency limit (or geometric optics). The energy conservation criterion is investigated versus the incidence angle, the permittivity of the lower medium, and the surface rms slope.     

Theoretical and computational aspects of scattering from rough surfaces: one-dimensional perfectly reflecting surfaces

We discuss the scattering of acoustic or electromagnetic waves from one-dimensional rough surfaces. We restrict the discussion in this report to perfectly reflecting Dirichlet surfaces (TE polarization). The theoretical development is for both infinite and periodic surfaces, the latter equations being derived from the former. We include both derivations for completeness of notation. Several theoretical developments are presented. They are characterized by integral equation solutions for the surface current or normal derivative of the total field. All the equations are discretized to a matrix system and further characterized by the sampling of the rows and columns of the matrix which is accomplished in either coordinate space (C) or spectral space (S). The standard equations are referred to here as CC equations of either the first (CC1) or second kind (CC2). Mixed representation, or SC-type, equations are solved as well as SS equations fully in spectral space.

Computational results are presented for scattering from various periodic surfaces. The results include examples with grazing incidence, a very rough surface and a highly oscillatory surface. The examples vary over a parameter set which includes the geometrical optics regime, physical optics or resonance regime, and a renormalization regime.

The objective of this study was to determine the best computational method for these problems. Briefly, the SC method was the fastest, but it did not converge for large slopes or very rough surfaces for reasons we explain. The SS method was slower and had the same convergence difficulties as SC. The CC methods were extremely slow but always converged. The simplest approach is to try the SC method first. Convergence, when the method works, is very fast. If convergence does not occur with SC, then SS should be used, and failing that CC.     

Theoretical and computational aspects of scattering from rough surfaces: one-dimensional perfectly reflecting surfaces

Abstract

We discuss the scattering of acoustic or electromagnetic waves from one-dimensional rough surfaces. We restrict the discussion in this report to perfectly reflecting Dirichlet surfaces (TE polarization). The theoretical development is for both infinite and periodic surfaces, the latter equations being derived from the former. We include both derivations for completeness of notation. Several theoretical developments are presented. They are characterized by integral equation solutions for the surface current or normal derivative of the total field. All the equations are discretized to a matrix system and further characterized by the sampling of the rows and columns of the matrix which is accomplished in either coordinate space (C) or spectral space (S). The standard equations are referred to here as CC equations of either the first (CC1) or second kind (CC2). Mixed representation, or SC-type, equations are solved as well as SS equations fully in spectral space.

Computational results are presented for scattering from various periodic surfaces. The results include examples with grazing incidence, a very rough surface and a highly oscillatory surface. The examples vary over a parameter set which includes the geometrical optics regime, physical optics or resonance regime, and a renormalization regime.

The objective of this study was to determine the best computational method for these problems. Briefly, the SC method was the fastest, but it did not converge for large slopes or very rough surfaces for reasons we explain. The SS method was slower and had the same convergence difficulties as SC. The CC methods were extremely slow but always converged. The simplest approach is to try the SC method first. Convergence, when the method works, is very fast. If convergence does not occur with SC, then SS should be used, and failing that CC.     

Monostatic and bistatic statistical shadowing functions from a one-dimensional stationary randomly rough surface: II. Multiple scattering

Abstract

The integral equation model (IEM) has been developed over the last decade and it has become one of the most widely used theoretical models for rough-surface scattering in microwave remote sensing. In the IEM model the shadowing function is typically either omitted or a form based on geometric optics with single reflection is used. In this paper, a shadowing function for one-dimensional rough surfaces which incorporates multiple scattering, finite surface length and both monostatic and bistatic configurations is developed. For any uncorrelated process, the resulting equation can be expressed in terms of the monostatic statistical shadowing function with single reflection, derived in the preceding companion paper. The effect of correlation between the surface slopes and heights for a Gaussian surface is studied to illuminate the range over which such correlations can be ignored. It is found that while the correlation between surface slopes and heights in the monostatic statistical shadowing function with single reflection can be ignored, when calculating the average shadowing function with double reflection the correlation between slopes and heights between points must be incorporated.     

Monostatic and bistatic statistical shadowing functions from a one-dimensional stationary randomly rough surface: II. Multiple scattering

The integral equation model (IEM) has been developed over the last decade and it has become one of the most widely used theoretical models for rough-surface scattering in microwave remote sensing. In the IEM model the shadowing function is typically either omitted or a form based on geometric optics with single reflection is used. In this paper, a shadowing function for one-dimensional rough surfaces which incorporates multiple scattering, finite surface length and both monostatic and bistatic configurations is developed. For any uncorrelated process, the resulting equation can be expressed in terms of the monostatic statistical shadowing function with single reflection, derived in the preceding companion paper. The effect of correlation between the surface slopes and heights for a Gaussian surface is studied to illuminate the range over which such correlations can be ignored. It is found that while the correlation between surface slopes and heights in the monostatic statistical shadowing function with single reflection can be ignored, when calculating the average shadowing function with double reflection the correlation between slopes and heights between points must be incorporated.     

Monostatic and bistatic statistical shadowing functions from a one-dimensional stationary randomly rough surface according to the observation length: I. Single scattering

When solving electromagnetic rough-surface scattering problems, the effect of shadowing by the surface roughness often needs to be considered, especially as the illumination angle approaches grazing incidence. This paper presents the Ricciardi-Sato, as well as the Wagner and the Smith formulations for calculating the monostatic and bistatic statistical shadowing functions from a one-dimensional rough stationary surface, which are valid for an uncorrelated Gaussian process with an infinite surface length. In this paper, these formulations are extended to include a finite surface length and any uncorrelated process. The inclusion of a finite surface length is needed to extend the single-reflection shadowing function to the more general multiple-reflection case, presented in the following companion paper. Comparisons of these shadowing functions with the exact numerical solution for the shadowing (using surfaces with Gaussian and Lorentzian autocorrelation functions for a Gaussian process) shows that the Smith formulation without correlation is a good approximation, and that including correlation only weakly improves the model. This paper also presents a method to include the shadowing effect in the electromagnetic scattering problem.     

Bistatic scattering and depolarization by randomly rough surfaces: application to the natural rough surfaces in X-band

The scattering of waves by random rough surfaces has important applications in the remote sensing of oceans and land. The problem of developing a model for rough surfaces is very difficult since, at best, the scattering coefficient σ⁰ is dependent upon (at least) the radar frequency, geometrical and physical parameters, incident and observation angles, and polarization. The problem of electromagnetic scattering from a randomly rough surface is analysed using the Kirchhoff approximation (stationary phase, scalar approximation), the small-perturbation model and the two-scale models. A first major new consideration in this paper is the polarimetric signature calculations as a function of the transmitter location and receiver location for a bistatic radio-link. We calculate the like- and cross-polarized received power directly using the scattering coefficients, without calculating the Mueller matrix. Next, a study of the regions of validity of the Kirchhoff and small-perturbation rough surface scattering models (in the bistatic case) is presented. Comparisons between the numerical calculations and the models are made for various surface rms heights and correlation lengths both normalized to the incident wavenumber (denoted by σ and L, respectively). By using these two parameters to form a two-dimensional space, the approximate regions of validity are then established. The second major new consideration is the development of a theoretical two-scale model describing bistatic reflectivity as well as the numerical results computed for the bistatic radar cross section from rough surfaces especially from the sea and snow-covered surfaces. The results are used to show the Brewster angle effect on near-grazing angle scattering.     

Bistatic scattering and depolarization by randomly rough surfaces: application to the natural rough surfaces in X-band

Abstract

The scattering of waves by random rough surfaces has important applications in the remote sensing of oceans and land. The problem of developing a model for rough surfaces is very difficult since, at best, the scattering coefficient σ⁰ is dependent upon (at least) the radar frequency, geometrical and physical parameters, incident and observation angles, and polarization. The problem of electromagnetic scattering from a randomly rough surface is analysed using the Kirchhoff approximation (stationary phase, scalar approximation), the small-perturbation model and the two-scale models. A first major new consideration in this paper is the polarimetric signature calculations as a function of the transmitter location and receiver location for a bistatic radio-link. We calculate the like- and cross-polarized received power directly using the scattering coefficients, without calculating the Mueller matrix. Next, a study of the regions of validity of the Kirchhoff and small-perturbation rough surface scattering models (in the bistatic case) is presented. Comparisons between the numerical calculations and the models are made for various surface rms heights and correlation lengths both normalized to the incident wavenumber (denoted by σ and L, respectively). By using these two parameters to form a two-dimensional space, the approximate regions of validity are then established. The second major new consideration is the development of a theoretical two-scale model describing bistatic reflectivity as well as the numerical results computed for the bistatic radar cross section from rough surfaces especially from the sea and snow-covered surfaces. The results are used to show the Brewster angle effect on near-grazing angle scattering.     

Implementation and validation of the curvilinear coordinate method for the scattering of electromagnetic waves from two-dimensional dielectric random rough surfaces

The curvilinear coordinate method is applied for analysing 2-D dielectric random rough surfaces. The theory is based on Maxwell's equations written in a non-orthogonal coordinate system. For each medium, this method leads to an eigenvalue system. The scattered fields within two media are expanded as linear combinations of eigensolutions satisfying the outgoing wave condition. The boundary conditions allow the scattering amplitudes to be determined. The coherent and incoherent intensities are estimated by averaging the scattering amplitudes over several realizations. The theory is verified by comparison with results obtained by other exact method. A discussion on the C-method and the Sparse-Matrix CAnonical Grid method is proposed in terms of accuracy and computation time.     

Theoretical study of the Kirchhoff integral from a two-dimensional randomly rough surface with shadowing effect: application to the backscattering coefficient for a perfectly-conducting surface

In this paper, the backscattering coefficient of a two-dimensional randomly rough perfectly-conducting surface is investigated using the Kirchhoff approach with a shadowing function. The rough surface height/slope correlations assumed to be Gaussian are accounted for in this analysis. The scattering coefficient is then formulated in terms of a characteristic function for the integrations over the surface heights, in terms of expected values for the integrations over the surface slopes. Numerical comparisons of Kirchhoff's approach (KA) with the stationary-phase (SP) approximation are made with respect to the choice of the one-dimensional surface height autocorrelation function and the shadowing effect. For an isotropic surface the results show that SP underestimated the incoherent backscattering coefficient compared with KA. Moreover, when the correlation between the slopes and the heights is neglected, the shadowing effect may be ignored.     

Theoretical study of the Kirchhoff integral from a two-dimensional randomly rough surface with shadowing effect: application to the backscattering coefficient for a perfectly-conducting surface

Abstract

In this paper, the backscattering coefficient of a two-dimensional randomly rough perfectly-conducting surface is investigated using the Kirchhoff approach with a shadowing function. The rough surface height/slope correlations assumed to be Gaussian are accounted for in this analysis. The scattering coefficient is then formulated in terms of a characteristic function for the integrations over the surface heights, in terms of expected values for the integrations over the surface slopes. Numerical comparisons of Kirchhoff's approach (KA) with the stationary-phase (SP) approximation are made with respect to the choice of the one-dimensional surface height autocorrelation function and the shadowing effect. For an isotropic surface the results show that SP underestimated the incoherent backscattering coefficient compared with KA. Moreover, when the correlation between the slopes and the heights is neglected, the shadowing effect may be ignored.     

The curvilinear coordinate method associated with the short-coupling-range approximation for the study of scattering from one-dimensional random rough surfaces

This paper deals with the scattering of an incident plane wave from a mono-dimensional rough surface. The analysis of the scattering phenomenon is performed using the curvilinear coordinate method (C method) associated with the short-coupling-range approximation (SCRA) and the Huygens principle. The C method is based on the solving of Maxwell’s equations under their covariant form written in a nonorthogonal coordinate system which fits the scattering surface profile. This leads to eigenvalue problems. The scattered surface fields are expressed as linear combinations of eigensolutions and the combination coefficients are determined using the boundary conditions. Electromagnetic fields are obtained using the surface fields and the Huygens principle. The short-coupling-range approximation (SCRA) applied with the C method allows a significant saving in computation time. We confirm the efficiency and the validity of this new approach with respect to the C method.