Analytical comparison between the surface current integral equation and the second-order small-slope approximation

This paper is the third in a series discussing a new approximate bistatic model for electromagnetic scattering from perfectly conducting rough surfaces. Our previous approach supplemented the Kirchhoff model through the addition of new terms involving linear orders in slope and surface elevation differences that arise naturally from a second iteration of the surface current integral equation. This completion of the Kirchhoff was shown to provide the correct first-order small perturbation method (SPM-1) in the general bistatic context. The agreement with SPM-1 was achieved because differences of surface heights are no longer expanded in powers of surface slope. While consistent with SPM, our previous formulation fails to reconverge toward the Kirchhoff model, at some incidence and scattered angles, when the illuminated surface satisfies the high frequency roughness condition. This weakness is also shared with the first-order small slope approximation (SSA-1) which is structurally equivalent to our previous formulation where the polarization is independent of surface roughness. The second-order small slope approximation (SSA-2), which satisfies the SPM-1 and second-order small perturbation method (SPM-2) limits by construction, was shown by Voronovich to converge toward the tangent plane approximation of the Kirchhoff model under high frequency conditions. In the present paper, we show that, in addition to the linear orders in our previous model, one must now include cross-terms between slope and surface elevation to ensure convergence toward both high frequency and small perturbation limits. With the inclusion of these terms, our new formulation becomes comparable to the SSA-2 (second-order kernel) without the need to evaluate all the quadratic order slope and elevations terms. SSA-2 is more complete, however, in the sense that it guarantees convergence toward the second-order Bragg limit (SPM-2) in the fully dielectric case in addition to both SPM-1 and Kirchhoff. Our new generalization is shown to explain correctly extra depolarization in specular conditions to be caused by surface curvature and surface autocorrelation for incoherent and coherent scattering, respectively. This result will have large repercussions on the interpretation of bistatically reflected signals such as those from GPS.     

Local and non-local curvature approximation: a new asymptotic theory for wave scattering

Abstract

We present a new asymptotic theory for scalar and vector wave scattering from rough surfaces which federates an extended Kirchhoff approximation (EKA), such as the integral equation method (IEM), with the first and second order small slope approximations (SSA). The new development stems from the fact that any improvement of the ‘high frequency’ Kirchhoff or tangent plane approximation (KA) must come through surface curvature and higher order derivatives. Hence, this condition requires that the second order kernel be quadratic in its lowest order with respect to its Fourier variable or formally the gradient operator. A second important constraint which must be met is that both the Kirchhoff approximation (KA) and the first order small perturbation method (SPM-1 or Bragg) be dynamically reached, depending on the surface conditions. We derive herein this new kernel from a formal inclusion of the derivative operator in the difference between the polarization coefficients of KA and SPM-1. This new kernel is as simple as the expressions for both Kirchhoff and SPM-1 coefficients. This formal difference has the same curvature order as SSA-1 + SSA-2. It is acknowledged that even though the second order small perturbation method (SPM-2) is not enforced, as opposed to the SSA, our model should reproduce a reasonable approximation of the SPM-2 function at least up to the curvature or quadratic order. We provide three different versions of this new asymptotic theory under the local, non-local, and weighted curvature approximations. Each of these three models is demonstrated to be tilt invariant through first order in the tilting vector.     

A practical second-order electromagnetic model in the quasi-specular regime based on the curvature of a 'good-conducting' scattering surface

This letter presents an approximate second-order electromagnetic model where polarization coefficients are surface dependent up to the curvature order in the quasi-specular regime. The scattering surface is considered 'good-conducting' as opposed to the case for our previous derivation where perfect conductivity was assumed. The model reproduces dynamically, depending on the properties of the scattering surface, the tangent-plane (Kirchhoff) or the first-order small-perturbation (Bragg) limits. The convergence is assumed to be ensured by the surface curvature alone. This second-order model is shown to be consistent with the small-slope approximation of Voronovich (SSA-1+SSA-2) for perfectly conducting surfaces. Our model differs from SSA-1 + SSA-2 in its dielectric expression, to correct for a full convergence toward the tangent-plane limit under the 'good-conducting' approximation. This new second-order formulation is simple because it involves a single integral over the scattering surface and therefore it is suitable for a vast array of analytical and numerical applications in quasi-specular applications.     

A practical second-order electromagnetic model in the quasi-specular regime based on the curvature of a ‘good-conducting’ scattering surface

Abstract

This letter presents an approximate second-order electromagnetic model where polarization coefficients are surface dependent up to the curvature order in the quasi-specular regime. The scattering surface is considered ‘good-conducting’ as opposed to the case for our previous derivation where perfect conductivity was assumed. The model reproduces dynamically, depending on the properties of the scattering surface, the tangent-plane (Kirchhoff) or the first-order small-perturbation (Bragg) limits. The convergence is assumed to be ensured by the surface curvature alone. This second-order model is shown to be consistent with the small-slope approximation of Voronovich (SSA-1+SSA-2) for perfectly conducting surfaces. Our model differs from SSA-1 + SSA-2 in its dielectric expression, to correct for a full convergence toward the tangent-plane limit under the ‘good-conducting’ approximation. This new second-order formulation is simple because it involves a single integral over the scattering surface and therefore it is suitable for a vast array of analytical and numerical applications in quasi-specular applications.     

Bistatic scattering from an anisotropic sea surface: Numerical comparison between the first-order SSA and the TSM models

The first-order small-slope approximation (SSA-1) model is used for numerical predictions of the normalized radar cross section (NRCS) of an anisotropic ocean surface in bistatic configurations for the K_u-band radar frequency. The calculations were made by assuming the Elfouhaily et al. surface-height spectrum for fully developed seas. In the forward-backward case, the SSA-1 presents an agreement with the geometric optics limit of the Kirchhoff approximation results in the near-specular directions where it is well known that the last model works well. In the fully bistatic case, SSA-1 numerical results are compared with those of the two-scale model in several configurations as a function of wind speed, wind direction, incident/scattering angles and for co-and cross-polarization states. Good agreement between the two models is noted in the co-and cross-polarization case with a small difference of about 1-2 dB. But in certain configurations, the SSA-1 model tends to overestimate the radar cross section peak behaviour. This irregularity is discussed and interpreted. The main purpose of this paper is to analyse NRCS predictions based on the SSA-1 model in a fully bistatic configuration.     

A new bistatic model for electromagnetic scattering from perfectly conducting random surfaces

In this paper, we extend the Kirchhoff approach, which is widely used for near-nadir backscattering calculations, to include the proper polarization sensitivity for general bistatic scattering from gently sloping, perfectly conducting surfaces. Previously, Holliday has shown how the inclusion of terms from the second iteration of the surface-current integral equation is required to obtain agreement with the small perturbation method for backscattering conditions. Here we employ a similar approach by retaining all terms in this iterative expansion through first order in the surface slope to derive expressions for the standard Kirchhoff field as well as for a supplementary field that contains the polarization sensitivity. A polarization vector notation is introduced to simplify the inclusion of tilting effects from larger-scale features on the scattering surface. In connection with this latter development, we provide a clarification of the earlier work by Valenzuela on this topic together with an extension to the bistatic problem. These extensions to the standard Kirchhoff approach form the basis for our composite bistatic scattering model which should provide a convenient and powerful tool for calculations involving passive as well as active microwave scattering from random surfaces.     

Extension of the local curvature approximation to third order and improved tilt invariance

The second-order local curvature approximation (LCA2) is a theory of rough surface scattering that reproduces fundamental low and high frequency limits in a tilted frame of reference. Although the existing LCA2 model provides agreement with the first order small perturbation method up to the first order in surface tilt, results reported in this paper produce a new formulation of the model that achieves consistency with perturbation theory to first order in surface height and arbitrary order in surface tilt. In addition, extension of the modified LCA to third order is presented, and allows the theory to match the second-order small perturbation method to arbitrary order in surface tilt. Crucial to the development of the theory are a set of identities involving relationships among the small perturbation method (i.e. low frequency) and Kirchhoff approximation (i.e. high frequency) kernels; a set of new identities obtained in our derivations is also presented. Sample results involving 3D electromagnetic scattering from penetrable rough surfaces, as well as 2D scattering from Dirichlet sinusoidal gratings, are provided to compare the new results with the existing LCA2 model and with other rough surface scattering theories.     

Extension of the local curvature approximation to third order and improved tilt invariance

The second-order local curvature approximation (LCA2) is a theory of rough surface scattering that reproduces fundamental low and high frequency limits in a tilted frame of reference. Although the existing LCA2 model provides agreement with the first order small perturbation method up to the first order in surface tilt, results reported in this paper produce a new formulation of the model that achieves consistency with perturbation theory to first order in surface height and arbitrary order in surface tilt. In addition, extension of the modified LCA to third order is presented, and allows the theory to match the second-order small perturbation method to arbitrary order in surface tilt. Crucial to the development of the theory are a set of identities involving relationships among the small perturbation method (i.e. low frequency) and Kirchhoff approximation (i.e. high frequency) kernels; a set of new identities obtained in our derivations is also presented. Sample results involving 3D electromagnetic scattering from penetrable rough surfaces, as well as 2D scattering from Dirichlet sinusoidal gratings, are provided to compare the new results with the existing LCA2 model and with other rough surface scattering theories.     

Study of the Electromagnetic Scattering from the Very Rough Fractal Sea Surface

In this paper, based on the fundamental formulae of the first-order and second-order Kirchhoff approx-imation mad with consideration of the shadowing effect, the backscattering enhancement of the one-dimensional very rough fractal sea surface with Pierson-Moskowitz spectrum is studied under the second-order Kirchhoff approximation at microwave frequency. The numerical results are compared with those of the first-order Kirchhoff approximation and integral equation method. The dependencies of the bistatic scattering cross section and the backscattering enhancement on the incident angle, fractal dimension, and windspeed over the sea surface are analyzed in detail.     

Tilt-invariant theory of rough-surface scattering: I

Abstract

The small-slope approximation (SSA) in rough-surface scattering theory uses the surface slope as a small parameter of expansion. But, from the physical point of view, the slope may not be a restrictive parameter because we can change the slope of a surface simply by tilting the coordinate system. We present the theory of rough-surface scattering in a coordinate-invariant form. The new method, tilt-invariant approximation (TIA), leads to a different expansion that does not require that the slope of a surface be small. For a small Rayleigh parameter this approximate solution provides the correct perturbation theory, for a large Rayleigh parameter it provides the Kirchhoff approximation with several correcting terms.     

Study of the Electromagnetic Scattering from the Very Rough Fractal Sea Surface

In this paper, based on the fundamental formulae of the first-order and second-order Kirchhoff approx-imation and with consideration of the shadowing effect, the backscattering enhancement of the one-dimensional veryrough fractal sea surface with Pierson-Moskowitz spectrum is studied under the second-order Kirchhoff approximationat microwave frequency. The numerical results are compared with those of the first-order Kirchhoff approximation andintegral equation method. The dependencies of the bistatic scattering cross section and the backscattering enhancementon the incident angle, fractal dimension, and windspeed over the sea surface are analyzed in detail.     

Tilt-invariant theory of rough-surface scattering: I

The small-slope approximation (SSA) in rough-surface scattering theory uses the surface slope as a small parameter of expansion. But, from the physical point of view, the slope may not be a restrictive parameter because we can change the slope of a surface simply by tilting the coordinate system. We present the theory of rough-surface scattering in a coordinate-invariant form. The new method, tilt-invariant approximation (TIA), leads to a different expansion that does not require that the slope of a surface be small. For a small Rayleigh parameter this approximate solution provides the correct perturbation theory, for a large Rayleigh parameter it provides the Kirchhoff approximation with several correcting terms.     

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

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.     

The IEM2M rough-surface scattering model for complex-permittivity scattering media

The integral equation model (IEM) was developed in the late 1980s and arguably became the most cited and implemented rough-surface scattering model in the field of radar remote sensing for Earth observation. It was derived by applying a second-order iteration in the incident electromagnetic field to the integral equations of the surface fields as given by Poggio and Miller. It is thus an extension of the first-order, Born approximation of these equations that produce the classical Kirchhoff approximation. The IEM has been subject to numerous amendments and variations over the last 20 years due to the imperfect introduction and handling of the Weyl representation of the spherical wave in its first version. The work presented here is a further development of the contribution made by the same author in 2001 (IEM2M), which was the first version of IEM able to blend analytically both the Kirchhoff and the small-perturbation approximations for the bistatic case. The improvement reported in this article is concerned with the inclusion of evanescent waves in the formulation of the model and the extension of the range of applicability of the second-order scattering terms to interfaces with complex-permittivity scattering media.     

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.     

The weighted curvature approximation in scattering from sea surfaces

A family of unified models in scattering from rough surfaces is based on local corrections of the tangent plane approximation through higher-order derivatives of the surface. We revisit these methods in a common framework when the correction is limited to the curvature, that is essentially the second-order derivative. The resulting expression is formally identical to the weighted curvature approximation, with several admissible kernels, however. For sea surfaces under the Gaussian assumption, we show that the weighted curvature approximation reduces to a universal and simple expression for the off-specular normalized radar cross-section (NRCS), regardless of the chosen kernel. The formula involves merely the sum of the NRCS in the classical Kirchhoff approximation and the NRCS in the small perturbation method, except that the Bragg kernel in the latter has to be replaced by the difference of a Bragg and a Kirchhoff kernel. This result is consistently compared with the resonant curvature approximation. Some numerical comparisons with the method of moments and other classical approximate methods are performed at various bands and sea states. For the copolarized components, the weighted curvature approximation is found numerically very close to the cut-off invariant two-scale model, while bringing substantial improvement to both the Kirchhoff and small-slope approximation. However, the model is unable to predict cross-polarization in the plane of incidence. The simplicity of the formulation opens new perspectives in sea state inversion from remote sensing data.     

A simplified asymptotic theory for ocean surface electromagnetic wave scattering

The normalized radar cross-section (NRCS) expression of the Local Curvature Approximation (LCA-1) is derived to first order. The polarization sensitivity of this model is compared to the Kirchhoff Approximation (KA), Two-Scale Model (TSM), Small Slope Approximation (SSA-1) and Small Perturbation Method (SPM-1) to first order in the backscattering configuration. Analytical comparisons and numerical simulations show that LCA-1 and TSM could be rewritten with the same formulation and that their polarization sensitivities are comparable. Comparisons with experimental data acquired in C- and Ku-band reveal that the polarization sensitivities of these models are not adequate. However, the NRCS azimuth modulation predicted by LCA-1 is found to be dependent on polarization and sea surface roughness. This property of the LCA-1 model yields to an azimuth modulation for the polarization ratio. Based on the surface curvature correction concept, a simplified electromagnetic model is proposed. The curvature correction is restricted to the resonant wave-number of the sea roughness spectrum. This is found to reproduce the polarization ratio given by experimental data versus incidence angle and wind speed.     

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.     

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 coefficient from one- and two-dimensional random surfaces using the stationary phase and scalar approximation with shadowing effect: comparisons with experiments and application to the sea surface

Abstract

In this paper, the bistatic scattering coefficient from one- and two-dimensional random surfaces using the stationary phase method and scalar approximation with shadowing effect is investigated. Both of these approaches use the Kirchhoff integral. With the stationary phase, the bistatic cross section is formulated in terms of the surface height joint characteristic function where the shadowing effect is investigated. In the case of the scalar approximation, the scattering function is computed from the previous characteristic function and in terms of expected values for the integrations over the slopes, where the shadowing effect is analysed analytically. Both of these formulations are compared with experimental data obtained from a Gaussian one-dimensional randomly rough perfectly-conducting surface. With the stationary-phase method, the results are applied to a two-dimensional sea surface.