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.     

A monostatic illumination function with surface reflections from one-dimensional rough surfaces

When solving scattering or emissivity problems for rough surfaces, the shadowing effect is often taken into account. Furthermore, for rough surfaces with large root mean square slope, surface reflections of the incidence or emission ray should not be neglected, especially at large observation angles. In this paper, a model of the monostatic statistical illumination function for one-dimensional rough surfaces with single surface reflection is developed, which is based on the Smith illumination function. A Monte Carlo ray-tracing algorithm is used to evaluate the accuracy of the present model. It is shown that, when neglecting the correlation between heights and slopes of the surface, the present model agrees quite well with the Monte Carlo result. Moreover, the result is improved if the correlation between heights and slopes is taken into account. For practical purposes, an empirical factor is introduced to improve the performance of the uncorrelated first-order illumination function to avoid computing the correlated one, which takes a long computation time. Besides, the first-order illumination function is significant at large observation angles, which could be promising to overcome problems in models of surface infrared emissivity where underestimation occurs compared with experimental measurements.     

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.     

Shadowing function with single reflection from anisotropic Gaussian rough surface. Application to Gaussian, Lorentzian and sea correlations

In this paper, the monostatic (transmitter and receiver are located at the same place) and bistatic (transmitter and receiver are distinct) statistical shadowing functions from an anisotropic two-dimensional randomly rough surface are presented. This parameter is especially important in the case of grazing angles for computing the bistatic scattering coefficient in optical and microwave frequencies. The objective of this paper is to extend the previous work (Bourlier C, Berginc G and Saillard J 2002 Waves Random Media 12 145-74), valid for a one-dimensional surface, to a two-dimensional anistropic surface by considering a joint Gaussian process of surface slopes and heights separating two points of the surface. The monostatic average (statistical shadowing function average over the statistical variables) shadowing function is then performed in polar coordinates with respect to the incidence angle, the azimuthal direction and the surface height two-dimensional autocorrelation function. In addition, for a bistatic configuration, it depends on the incidence angle and azimuthal direction of the receiver. For Gaussian and Lorentzian correlation profiles and practically important power-type spectra such as the Pierson-Moskowitz sea roughness spectrum, the numerical solution, obtained from generating the surface Gaussian elevations (Monte Carlo method), is compared with the uncorrelated and correlated models. The results show that the correlation underestimates the shadow slightly, whereas the uncorrelated results weakly overpredict the shadow and are close to the numerical solution.     

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.     

A new bistatic model for electromagnetic scattering from randomly rough surfaces

In this paper we propose a bistatic model for electromagnetic scattering from a Gaussian rough surface with small to moderate heights. It is based on the integral equation formulation where the spectral representations of the Green's function and its gradient are in complete forms, a general approach similar to those used in the advanced integral equation model (AIEM) and the integral equation model for second-order multiple scattering (IEM2M). Yet this new model can be regarded as an extension to these two models on two accounts: first it has made fewer and less restrictive assumptions in evaluating the complementary scattering coefficient for single scattering, and second it contains a more rigorous analysis by the inclusion of the error function related terms for the cross- and complementary scattering coefficients, which stems from the absolute phase term in the spectral representation of the Green's function. It is expected that our result for the complementary scattering coefficient is more accurate and more general, even when the effect of the error function related terms is neglected. As a result, the proposed model is expected to have wider applicability with a better accuracy. Numerical simulations are provided to demonstrate the validity of the proposed model.     

Three dimensional analyses of scattering by pressure-release sinusoidal surfaces

Scattering by pressure-release sinusoidal surfaces in three dimensions is analyzed using the Fresnel phase approximation and realistic source and receiver directivity approximations. Geometrical shadowing and second-order scattering are explicitly included to explore the validity of the Kirchhoff approximation. No restrictions on the surface heights and slopes are made. The "goodness" of the resulting expressions is verified by requiring the pressure scattered by a sinusoidal surface to reduce to the image solution as the surface amplitude goes to zero. The first-order scattered pressure achieves a very good approximation to the image solution, and the second-order scattered pressure goes to zero, as expected, under this requirement. The theory is compared with available experimental scattering measurements, and the agreement is good. Because the slopes on the experimental surface are very steep, shadowing corrections are indispensible to achieving accurate first and second order scattering results. Shadowing has a greater impact on the scattering prediction than the second-order scattering contribution. This suggests that the Kirchhoff approximation may be more robust when incorporated into a theory including a detailed shadowing treatment as well as the Fresnel phase approximation and a good directivity approximation.     

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

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.     

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.     

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.     

Scattering of a Narrow Pulsed Wave Beam by a Randomly Rough,Locally Lambertian Surface under Conditions of Strong Shadowing

We consider the problem of scattering of a narrow pulsed wave beam irradiating a randomly rough surface in the atmosphere under conditions of strong shadowing of one surface element by another. Expressions for the average power recorded by a receiver are obtained for the case where a locally Lambertian surface with Gaussian distribution of heights and slopes is irradiated by a delta impulse. It is shown that the shadowing and the atmospheric turbulence give rise to a considerable distortion of the received optical signal. The obtained analytical expressions for the received power are in good agreement with the results of numerical calculations.     

一维粗糙面散射的射线追踪法研究

在考虑遮蔽效应的情况下,应用射线追踪法对一维导体和介质的高斯粗糙面的散射系数进行了研究,并分别计算了考虑一次反射和二次反射的散射系数。同时,利用蒙特卡罗法生成一维高斯粗糙面,计算了考虑一次、二次反射和遮蔽时不同均方根斜率粗糙面的平均散射系数。     

风驱粗糙海面背景下海雾对可见和红外光多次散射特性研究

通过测量空间辐射反演海雾气溶胶的微观特性是一种重要的遥感方法,但是海雾对辐射的反射会受到海面背景的影响。该工作研究了海雾与海面的耦合多次散射。利用Mie理论研究了海雾气溶胶的单次散射特性,利用辐射传输理论研究太阳光在无下垫面海雾中传输多次散射特性;在基尔霍夫近似下研究了风驱粗糙海面的散射特性,并考虑了海面的遮蔽效应,得到了风驱海面随风速的反射函数。根据累加法研究了海雾和海面的耦合散射,计算结果表明海雾在有粗糙海面作为背景的情况下,其整体反射有较大增强。     

Scattering of a Narrow Wave Beam on a Rough,Locally-Specular Surface in the Case of Pulsed Illumination

We consider scattering of a pulsed narrow wave beam on a rough surface with a locally-specular indicatrix. Analytical expressions for the average received power are obtained for normal distributions of heights and slopes of the rough surface in two cases in which the direction to the receiver is close to or strongly different from the direction of specular reflection. It is shown that in these cases, the received echo pulses have drastically different profiles determined by the source and receiver parameters, the scheme of sounding, and the variance of the heights and slopes of the rough surface. The obtained analytical expressions for the average received power agree well with the results of numerical simulations.     

One method of solving rough surface wave diffraction problems

The problems of wave diffraction on a rough surface is considered by reducing surface scattering to volume scattering and then using methods for study of multiple scattering developed in the theory of wave propagation in random-inhomogeneous media. Using this approach, the Dyson equation is solved in the Bourret approximation as is the Bethe-Salpeter equation in the ladder approximation for the correlation function of the scalar field scattered on an infinite ideally reflective statistically homogeneous surface. Solutions of these equations automatically consider shadowing of some portions of the surface by others.Translated from Izvestiya Vysshikh Uchebnykh Zavedenii, Radiofizika, Vol. 27, No. 1, pp. 65–70, January, 1984.The author is greatly indebted to A. B. Shmelev and A. G. Vinogradov for their valuable remarks and evaluation of the study.     

Effective probability density function of rough surface slopes when strong shadowing is present

It is shown that at low grazing angles, the slope probability density function (PDF) of the nonshadowed part of a rough surface can differ significantly from the slope PDF of the overall surface, if surface heights and slopes are functionally dependent. If the surface steepness has a tendency to increase with height, the effective slopes of the illuminated part of the surface can be significantly steeper than the average slope of the surface as a whole. This fact can play a crucial role in any theoretical interpretation of experimental results concerning radar scattering by the sea surface at low grazing angles.     

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.     

Three-dimensional radiative transfer in an anisotropically scattering, semi-infinite medium: generalized reflection function

Three-dimensional radiative transfer in an anisotropic scattering medium exposed to spatially varying, collimated radiation is studied. The generalized reflection function for a semi-infinite medium with a very general scattering phase function is the focus of this investigation. An integral transform is used to reduce the three-dimensional transport equation to a one-dimensional form, and a modified Ambarzumian's method is applied to formulate a nonlinear integral equation for the generalized reflection function. The integration is over both the polar and azimuthal angles; hence, the integral equation is said to be in the double-integral form. The double-integral, reflection function formulation can handle a variety of anisotropic phase functions and does not require an expansion of the phase function in a Legendre polynomial series. Complicated kernel transformations of previous single-integral studies are eliminated. Single and double scattering approximations are developed. Numerical results are presented for a Rayleigh phase function to illustrate the computational characteristics of the method and are compared to results obtained with the single-integral method. Agreement between the two approaches is excellent; however, as the transform variable increases beyond five the number of quadrature points required for the double-integral method to produce accurate solutions significantly increases. A new interpolation scheme produces accurate results when the transform variable is large.