Open access

Electromagnetic Models for Remote Sensing of Layered Rough Media

Written By

Pasquale Imperatore, Antonio Iodice and Daniele Riccio

Published: February 1st, 2010

DOI: 10.5772/9108

Chapter metrics overview

2,230 Chapter Downloads

View Full Metrics

1. Introduction

Each region of the Earth’s crust can be morphologically modeled as a suitable layered structure, in which some amount of roughness is presented by every interface. Actually, propagation in stratified soil, sand cover of arid regions, forest canopies, urban buildings, snow blanket, snow cover ice, sea ice and glaciers, oil flood on sea surface, and other natural scenes can be modeled referring to most likely discrete (piecewise-constant) systems, rather than continuous, with some amount of roughness presented by every interface. Moreover, a key issue in remote sensing of other Planets is to reveal the content under the surface illuminated by the sensors: also in this case a layered model is usually employed.

The aim of this chapter is to provide a structured presentation of the main theoretical and conceptual foundations for the problem of the electromagnetic wave interaction with layered rough media. In the first part, special emphasis is on the analytical models obtainable in powerful framework of the perturbation approach. The comprehensive scattering model based on the Boundary Perturbation Theory(BPT), which permits to systematically analyze the bi-static scattering patterns of 3D multilayered rough media, is then presented highlighting the formal connections with all the previously existing simplified perturbative models, as well as its wide relevance in the remote sensing applications scenario. The polarimetric Scattering Matrixof a multilayered medium with an arbitrary number of rough interfaces is also provided. The second part is devoted to a mathematical description which connects the concepts of localscattering and globalscattering. Consequently, a functional decomposition of the BPT global scattering solution in terms of basic single-scattering local processes is rigorously established. The scattering decompositiongives insight into the BPT analytical results, so enabling a relevant physical-revealing interpretation involving ray-series representation. Accordingly, in first-order limit, the way in which the character of the local scattering processes emerges is dictated by the nature of the structural filter action, which is inherently governed by the series of coherent interactions with the medium boundaries. As a result, the phenomenologically successful BPT model opens the way toward new techniques for solving the inverse problem, for designing SAR processing algorithms, and for modelling the time-domain response of layered structures.

Figure 1.

Geometry for an N-rough boundaries layered medium.


2. Problem definition

When stratified media with rough interfaces are concerned, the possible approaches to cope with the EM scattering problem fall within three main categories. First, the numerical approaches do not permit to attain a comprehensive understanding of the general functional dependence of the scattering response on the structure parameters, as well as do not allow capturing the physics of the involved scattering mechanisms. Layered structures with rough interfaces have been also treated resorting to radiative transfer theory(RT). However, coherent effects are not accounted for in RT theory and could not be contemplated without employing full wave analysis, which preserves phase information. Another approach relies on the full-wave methods. Although, to deal with the electromagnetic propagation and scattering in complex random layered media, several analytical formulation involving some idealized cases and suitable approximations have been conducted in last decades, the relevant solutions usually turn out to be too complicated to be generally useful in the remote sensing scenario, even if simplified geometries are accounted for. The proliferation of the proposed methods for the simulation of wave propagation and scattering in a natural stratified medium and the continuous interest in this topic are indicative of the need of appropriate modelling and interpretation of the complex physical phenomena that take place in realistic environmental structures. Indeed, the availability of accurate, sound physical and manageable models turns out still to be a strong necessity, in perspective to apply them in retrieving of add-valued information from the data acquired by microwave sensors. For instance, such models are high desirable for dealing with the inversion problem as well as for the effective design of processing algorithms and simulation of Synthetic Aperture Radar signals. Generally speaking, an exact analytical solution of Maxwellequations can be found only for a few idealized problems. Subsequently, appropriate approximation methods are needed. Regarding the perturbative approaches, noticeable progress has been attained in the investigation on the extension of the classical SPM (small perturbation method) solution for the scattering from rough surface to specific layered configurations. Most of previous existing works analyze different layered configurations in the first-order limit, using procedures, formalisms and final solutions that can appear of difficult comparison (Yarovoy et al., 2000), (Azadegan and Sarabandi, 2003), (Fuks, 2001). All these formulations, which refer to the case of a single rough interface, have been recently unified in (Franceschetti et al, 2008). On the other hand, solution for the case of two rough boundaries has also been proposed in (Tabatabaeenejad and Moghaddam, 2006).

Methodologically, we underline that all the previously mentioned existing perturbative approaches, followed by different authors in analyzing scattering from simplified geometry, imply an inherent analytical complexity, which precludes the treatment to structures with more than one (Fuks, 2001) (Azadegan et al., 2003) (Yarovoy et al., 2000) or two (Tabatabaeenejad er al., 2006) rough interfaces.

The general problem we intend to deal with here refers to the analytical evaluation of the electromagnetic scattering by layered structure with an arbitrary number of rough interfaces (see Figure 1). As schematically shown in Figure 1, an arbitrary polarized monochromatic plane wave


is considered to be incident on the layered medium at an angleθ0irelative to thez^direction from the upper half-space, where in the field expression a time factorexp(-jt)is understood, and where, using a spherical frame representation, the incident vector wave direction is individuated byθ0i,ϕ0i


whereki=kxix^+kyiy^is the two-dimensional projection of incident wave-number vector on the plane z=0. The parameters pertaining to layer m with boundaries –d m-1and -d mare distinguished by a subscript m. Each layer is assumed to be homogeneous and characterized by arbitrary and deterministic parameters: the dielectric relative permittivity m, the magnetic relative permeability μ mand the thickness m=d m-d m-1. With reference to Figure 1, it has been assumed that in particular, d 0=0. In the following, the symbol denotes the projection of the corresponding vector on the plan z=0. Herer=(r,z), so we distinguish the transverse spatial coordinatesr=(x,y)and the longitudinal coordinate z. In addition, each mth rough interface is assumed to be characterized by a zero-mean two-dimensional stochastic processζm=ζm(r)with normal vectorn^m. No constraints are imposed on the degree to which the rough interfaces are correlated.

A general methodology has been developed by Imperatore et al. to analytically treat EM bistatic scattering from this class of layered structures that can be described by small changes with respect to an idealized (unperturbed) structure, whose associated problem is exactly solvable. A thorough analysis of the results of this theoretical investigation (BPT), which is based on perturbation of the boundary condition, will be presented in the following, methodologically emphasizing the development of the several inherent aspects.

Figure 2.

Geometry for a flat boundaries layered medium


3. Preliminary notation and definitions

This section is devoted preliminary to introduce the formalism used in the following of this chapter. The Flat Boundaries layered medium (unperturbed structure) is defined as a stack of parallel slabs (Figure 2), sandwiched in between two half-spaces, whose structure is shift invariant in the direction of x and y (infinite lateral extent in x-y directions). With the notationsTm1|mpandRm1|mp, respectively, we indicate the ordinary transmission and reflection coefficients at the interface between the regions (m-1) and m, with the superscript p{v, h} indicating the polarization state for the incident wave and may stand for horizontal (h) or vertical (v) polarization (Tsang et al., 1985) (Imperatore et al. 2009a). In addition, we stress that:


The generalized reflection coefficientsm1|mp, for the p-polarization (TE or TM), at the interface between the regions (m-1) and m are defined as the ratio of the amplitudes of upward- and downward-propagating waves immediately above the interface, respectively. They can be expressed by recursive relations as in (Chew W. C., 1997) (Imperatore et al. 2009a):


Likewise, at the interface between the regions (m+1) and m,m+1|mpis given by:




wherekm=k0μmεmis the wave number for the electromagnetic medium in the mth layer, withk0=ω/c=2π/λ, and wherek=kxx^+kyy^is the two-dimensional projection of vector wave-number on the plane z=0. It should be noted that the factors


take into account the multiple reflections in the mth layer. On the other hand, the generalized transmission coefficients in downward direction can be defined as:


where p {v, h}. The generalized transmission coefficients in upward direction are then given by:


which formally express the reciprocity of the generalized transmission coefficients for an arbitrary flat-boundaries layered structure (Imperatore et al. 2009b). In addition, with reference to a layered slab sandwiched between two half-space, we consider the generalized transmission coefficients in upward direction for the layered slab between two half-spaces (m,0), which are defined as


Note also that


The generalized transmission coefficients in downward direction for the layered slab between two half-spaces (m,0), can be defined as


On the other hand, it should be noted that the0|mpare distinct from the coefficients0|mp(slab), because in the evaluation of0|mpthe effect of all the layers under the layer m is taken into account, whereas0|mp(slab)are evaluated referring to a different configuration in which the intermediate layers 1...m are bounded by the half-spaces 0 and m. In the following, we shown how the employing the generalized reflection/transmission coefficient notions not only is crucial in obtaining a compact closed-form perturbation solution, but it also permit us to completely elucidate the obtained analytical expressions from a physical point of view, highlighting the role played by the equivalent reflecting interfaces and by the equivalent slabs, so providing the inherent connection between local and global scattering responses.


4. Spectral Representation of the Stochastic Geometry Description

In this section, the focus is on stochastic description for the geometry of the investigated structure, and the notion of wide-sense stationary process is detailed. First of all, when the description of a rough interface by means of deterministic functionζm(r)is concerned, the corresponding ordinary 2-D Fourier Transform pair can be defined as


Let us assume now thatζm(r), which describes the generic (mth) rough interface, is a 2-D stochastic process satisfying the conditions


where the angular bracket denotes statistical ensemble averaging, and whereBζm(ρ)is the interface autocorrelation function, which quantifies the similarity of the spatial fluctuations with a displacement . Equations (17)-(18) constitute the basic assumptions defining a wide sense stationary (WSS) stochastic process: the statistical properties of the process under consideration are invariant to a spatial shift. Similarly, concerning two mutually correlated random rough interfacesζmandζn, we also assume that they are jointly WSS, i.e.


whereBζmζn(ρ)is the corresponding cross-correlation function of the two random processes.

It can be readily derived that


The integral in (15) is a Riemann integral representation forζm(r), and it exists ifζm(r)is piecewise continuous and absolutely integrable. On the other hand, when the spectral analysis of a stationary random process is concerned, the integral (15) does not in general exist in the framework of theory of the ordinary functions. Indeed, a WSS process describing an interfaceζm(r)of infinite lateral extension, for its proper nature, is not absolutely integrable, so the conditions for the existence of the Fourier Transform are not satisfied. In order to obtain a spectral representation for a WSS random process, this difficulty can be circumvented by resorting to the more general Fourier-Stieltjes integral (Ishimaru, 1978); otherwise one can define space-truncated functions. When a finite patch of the rough interface with area A is concerned, the space-truncated version of (15) can be introduced as


subsequently,ζ˜m(k)=limAζ˜m(kA)is not an ordinary function. Nevertheless, we will use again the (15)-(16), regarding them as symbolic formulas, which hold a rigorous mathematical meaning beyond the ordinary function theory (generalized Fourier Transform). We underline that by virtue of the condition (17) directly follows also thatζ˜m(k)=0. Let us consider


where the asterisk denotes the complex conjugated, and where the operations of average and integration have been interchanged. When jointly WSS processesζmandζnare concerned, accordingly to (19), the RHS of (22) must be a function ofrronly; therefore, it is required that


where (∙) is the Dirac delta function, and whereWmn(κ)is called the (spatial) cross power spectral density of two interfacesζmandζn, for the spatial frequencies of the roughness. Equation (23) states that the different spectral components of the two considered interfaces must be uncorrelated. This is to say that the (generalized) Fourier transform of jointly WSS processes are jointly non stationary white noise with average powerWmn(k). Indeed, by using (23) into (22), we obtain


where the RHS of (24) involves an (ordinary) 2D Fourier Transform. Note also that as a direct consequence of the fact thatζn(r)is real we have the relationζ˜n(k)=ζ˜n*(k). Therefore, settingρ=rrin (24), we have


The cross-correlation functionBζmζn(ρ)of two interfacesζmandζn, is then given by the (inverse) 2D Fourier Transform of their (spatial) cross power spectral density, and Equation (25) together with its Fourier inverse


may be regarded as the (generalized) Wiener-Khinchin theorem. In particular, when n=m, (23) reduces to


whereWm(κ)is called the (spatial) power spectral density of nth corrugated interfaceζmand can be expressed as the (ordinary) 2D Fourier transform of n-corrugated interface autocorrelation function, i.e., satisfying the transform pair:


which is the statement of the classical Wiener-Khinchin theorem. We emphasize the physical meaning ofWm(κ)dκ=Wm(κx,κ)ydκxdκy: it represents the power of the spectral components of the mth rough interface having spatial wave number between xand x+d xand yand y+d y, respectively, in x and y direction. Furthermore, from (20) and (26) it follows that


This is to say that, unlike the power spectral density, the cross power spectral density is, in general, neither real nor necessarily positive. Furthermore, it should be noted that the Dirac’s delta function can be defined by the integral representation


By using in (27) and (23) the relationδ(0;A)=A/(2π)2, we have, respectively, that the (spatial) power spectral density of nth corrugated interface can be also expressed as


and the (spatial) cross power spectral density of two interfacesζmandζnis given by


It should be noted that the domain of a rough interface is physically limited by the illumination beamwidth. Note also that the different definitions of the Fourier transform are available and used in the literature: the sign of the complex exponential function are sometimes exchanged and a multiplicative constant(2π)2may appear in front of either integral or its square root in front of each expression (15)-(16). Finally, we recall that the theory of random process predicts only the averages over many realizations.


5. Perturbative Field Formulation

With reference to the geometry of Figure 1, in order to obtain a solution valid in each region of the structure, we have to enforce the continuity of the tangential fields:


whereΔEm=Em+1EmΔHm=Hm+1Hm, and the surface normal vector is given by:


with the slope vectorγm


and whereis the nabla operator in the x-y plane. In order to study the fieldsEmandHmwithin the generic mth layer of the structure, we assume then that, for each mth rough interface, the deviations and slopes of surface with respect to the reference mean plane z=-d mare small enough in the sense of (Ulaby et al, 1982) (Tsang et al., 1985), so that the fields can be expanded about the reference mean plane. The fields expansion can be then injected into the boundary conditions (34)-(35), so that, retaining only up to the first-order terms, the following nonuniform boundary conditions can be obtained (Imperatore et al. 2009a)


where the field solution has been formally represented as:


Therefore, the boundary conditions from each mth rough interface can be transferred to the associated equivalent flat interface. In addition, the right-hand sides of Eqs. (38) and (39) can be interpreted as effective magnetic (JHmp(1)) and electric (JEmp(1)) surface current densities, respectively, with p denoting the incident polarization; so that we can identify the first-order fluctuation fields as being excited by these effective surface current densities imposed on the unperturbed interfaces. Accordingly, the geometry randomness of each corrugated interfaces is then translated in random current sheets imposed on each reference mean plane (z=-d m), which radiate in an unperturbed (flat boundaries) layered medium. As a result, within the first-order approximation, the field can be than represented as the sum of an unperturbed partEn(0),Hn(0)and a random part, so thatEn(rz)En(0)+En(1)Hn(rz)Hn(0)+Hn(1). The first is the primary field, which exists in absence of surface boundaries roughness (flat-boundaries stratification), detailed in (Imperatore et al. 2009a); whereasEn(1),Hn(1)can be interpreted as the superposition of single-scatter fields from each rough interface. In order to perform the evaluation of perturbative development, the scattered field is then represented as the sum of up- and down-going waves, and the first-order scattered field in each region of the layered structure can be represented in the form:


Therefore, a solution valid in each region of the structure can be obtained from (42)-(45) taking into account the non uniform boundary conditions (38)-(39). In order to solve the scattering problem in terms of the unknown expansion coefficientsSm±q(1)(k), we arrange their amplitudes in a single vector according to the notation:


Subsequently, the nonuniform boundary conditions (38)-(39) can be formulated by employing a suitable matrix notation, so that for the (q=h) horizontal polarized scattered wave we have (Imperatore et al. 2008a) (Imperatore et al. 2009a):




is the term associated with the effective source distribution, wherein the expressions of the effective currentsJ˜Emp(1)andJ˜Hmp(1), imposed on the (flat) unperturbed boundary z = −d m, for an incident polarization p {v, h} are detailed in (Imperatore et al. 2009a); and where Z 0is the intrinsic impedance of the vacuum. Furthermore, the fundamental transfer matrix operator is given by:


with the superscripts q {v, h} denoting the polarization. Moreover, it should be noted that on a (kth) flat interface Eq. (47) reduces to the uniform boundary conditions, thus getting:


We emphasize that Eqs. (47) states in a simpler form the problem originally set by Eqs. (38)-(39): as matter of fact, solving Eq. (47) m implies dealing with the determination of unknown scalar amplitudesSm±q(1)(k)instead of working with the corresponding vector unknownsEm(1),Hm(1). Therefore, the scattering problem in each mth layer is reduced to the algebraic calculation of the unknown expansion scattering coefficients vector (46). As a result, when a structure with rough interfaces is considered, the enforcement of the original non uniform boundary conditions through the stratification (m=0,..., N-1) can be addressed by writing down a linear system of equations with the aid of the matrix formalism (47)-(48) with m=0,..., N-1. As a result, the formulation of non-uniform boundary conditions in matrix notation (47)-(48) enables a systematic method for solving the scattering problem: For the N-layer stratification of Figure 1, we have to find 2N unknown expansion coefficients, using N vectorial equations (47), i.e., 2N scalar equations. It should be noted that, for the considered configuration, the relevant scattering coefficientsSN+q(1)(k),S0q(1)(k)are obviously supposed to be zero. The scattering problem, therefore, results to be reduced to a formal solution of a linear equation system. We finally emphasize that here we are interested in the scattering from the stratification; therefore, the determination of the unknown expansion coefficientsS0+q(1)(k)of the scattered wave into the upper half-space is required only. Full expressions for these coefficients are reported in (Imperatore et al. 2009a).


6. BPT Closed-form Solution

The field scattered upward in the upper half-space in the first-order limit can be written in the form (see (42)-(45)):


By employing the method of stationary phase, we evaluate the integral (51) in the far field zone, obtaining:


with q {v, h} is the polarization of the scattered field. The scattering cross section of a generic (nth) rough interface embedded in the layered structure can be then defined as


where < > denotes ensemble averaging, where q {v, h} and p {v, h} denote, respectively, the polarization of scattered field and the polarization of incident field, and where A is the illuminated surface area. The estimate of the mean power density can be obtained by averaging over an ensemble of statistically identical interfaces. Therefore, considering that the (spatial) power spectral densityWn(κ)of nth corrugated interface is defined as in (32), the scattering cross section relative to the contribution of the nth corrugated interface, according to the formalism used in [Franceschetti et al. 2008], can be expressed as:


with p, q {v, h} denoting, respectively, the incident and the scattered polarization states, which may stand for horizontal polarization (h) or vertical polarization (v); the coefficientsα˜qpm,m+1are relative to the p-polarized incident wave impinging on the structure from upper half space 0 and to the q-polarized scattering contribution from structure into the upper half space, originated from the rough interface between the layers m, m+1:


wherem|0p(k)are the generalized transmission coefficients in upward direction (see (11)). Furthermore, we stress when the backscattering case (k^s×k^i=0) is concerned, our cross-polarized scattering coefficients (55)-(58) evaluated in the plane of incidence vanish, in full accordance with the classical first-order SPM method for a rough surface between two different media (Ulaby et al, 1982) (Tsang et al., 1985).

We now show that the solution, given by the expression (55)-(58), is susceptible of a straightforward generalization to the case of arbitrary stratification with N-rough boundaries. Taking into account the contribution of each nth corrugated interface, the global bi-static scattering cross section of the N-rough interface layered media can be expressed as:


with p, q {v, h}, where the asterisk denotes the complex conjugated, whereα˜qpi,i+1are given by (55)-(58), and where the cross power spectral densityWij, between the interfaces i and j, for the spatial frequencies of the roughness is given by (33). As a result, the scattering from the rough layered media is sensitive to the correlation between rough profiles of different interfaces. In fact, a real layered structure will have interfaces cross-correlation somewhere between two limiting situations: perfectly correlated and uncorrelated roughness. Consequently, the degree of correlation affects the phase relation between the fields scattered from each rough interface. Obviously, when the interfaces are supposed to be uncorrelated, the second term in (59) vanishes and, in the first-order approximation, the total scattering from the structure arises from the incoherent superposition of radiation scattered from each interface.

As a result, an elegant closed form solution is established, which takes into account parametrically the dependence of scattering properties on structure (geometric and electromagnetic) parameters. In addition, as it will be shown in the following, the proposed global solution turns out to be completely interpretable with basic physical concepts, clearly discerning the physics of the involved scattering mechanisms.


7. Generalized Scattering Matrix

In this section, to emphasize the polarimetric character of the BPT solution, we introduce the generalized bistatic scattering matrix of the layered rough media, which can be then formally expressed by:


In particular,


characterizes the polarimetric response of the generic (mth) rough interface of the layered structure, for a plane wave incident in directionkiand for a given observation directionks, with


and wherein


where we have introduced the notation


It should be noted that (63)-(67) are obtained directly by (55)-(58) by making use of (11). Note also that the coefficientsα˜qpm,m+1are expressible in a direct closed-form and depend parametrically on the unperturbed structure parameters. We also emphasize that the scattering configuration we have considered is compliant with the classical Forward Scattering Alignment (FSA) convention adopted in radar polarimetry.

Denoting with the superscript T the transpose, it can be verified that the scattering matrix satisfies the following relationship (Imperatore et al. 2009b)


which concisely expresses the reciprocity principle of the electromagnetic theory, as expected. This is to say that the scattering experiment is invariant for interchanging the role of transmitter and receiver. Note that the inversion of the projections on the z=0 planekiandksare directly related to the inversion of the incident and scattered vector waveki=kikz0iz^andks=ks+kz0sz^, respectively.

As a result, the presented closed-form solution permits the polarimetric evaluation of the scattering for a bi-static configuration, once the three-dimensional layered structure’s parameters (shape of the roughness spectra, layers thickness and complex permittivities), the incident field parameters (frequency, polarization and direction) and the observation direction are been specified. Therefore, our formulation leads to a direct functional dependence (no integral evaluation is required) and, subsequently, allows us to show that the scattered field can be parametrically evaluated considering a set of parameters: some of them refer to an unperturbed structure configuration, i.e. intrinsically the physical parameters of the smooth boundary structure, and others which are determined exclusively by (random) deviations of the corrugated boundaries from their reference position. Procedurally, once the generalized reflection/transmission coefficients are recursively evaluated, the (63-67) can be than directly computed, so that the scattering cross section (59) or the generalized scattering matrix (60) of a structure with rough interfaces can be finally predicted. Finally, it should be noted that the method to be applied needs only the classical gently-roughness assumption, without any further approximation.


8. Unifying Perspective on Perturbation Approaches

In this section, we first summarize and discuss the previous existing scattering models introduced to cope with simplified layered geometry with only one (or two) rough interface, whose derivation methods belong to the class of perturbative methods. In the perspective of providing a unifying insight for the different perturbative formulations, the aim is to reconsider the state of art in an organized mathematical framework, analytically demonstrating the formal consistency of BPT general scattering solution, which permits to deal with layered media with an arbitrary number of rough interfaces, with the previous existing perturbative models, whose relevant first-order solutions can appear already of difficult mutual comparison (Franceschetti et al 2003) (Franceschetti et a. 2008).

In (Fuks, 2001) a model to calculate scattering from a rough surface on top of a stratified medium (see the geometry of Figure 3a) has been proposed. Expressions for scattering bi-static cross section were obtained by using the plane wave expansion of scattered EM fields and an equivalent current method, without using to the Green’s function formalism. With reference to the geometry pictured in Figure 3b, an analytical small-perturbation-based model was developed to deal with a slightly rough interface boundary covered with a homogeneous dielectric layer (Azadegan et al., 2003) ( Sarabandi et al, 2000). Starting from a perturbation series expansion and by employing the Green’s function formalism, a solution to predict the first order bi-static scattering coefficients was found. On the other hand, the backscattering problem from the two-middle layer structure with one embedded corrugation, as schematized in Figure 3c, was investigated in (Yarovoy et al., 2000) in the first-order approximation, by using the small perturbation method combined with the Green’s function approach. This approach leads to some cumbersome analytical expressions for backscattering coefficients.

As matter of fact, all these models, which refer to different simplified geometry, employ different perturbative procedures and different notations in the relative analytical derivation, so that the resulting solutions turn out mutually of difficult comparison. Besides, the finding of the connection between these existing functional forms is not a trivial task. With regard to these models, in (Franceschetti et al., 2008) it was essentially demonstrated the equivalence of the relevant analytical procedures and the consistency of the respective solutions. It should be mentioned that the models in (Fuks, 2001) and (Azadegan et al., 2003) (Sarabandi et al, 2000) are derived for a bi-static configuration. Conversely, the solution derived in (Yarovoy et al., 2000) with reference to the geometry of Figure 3c, which is even relatively more general since contemplates flat-boundaries stratification above and under the roughness, is given only in backscattering case. On the other hand, none of the pertinent configurations of these simplistic considered models is directly applicable to an actual remote sensing scenario. In fact, the natural stratified media are definitely constituted by corrugated interfaces, each one exhibiting a certain amount of roughness, whereas the flatness is an idealization which does not occur in natural media. More specifically, it can occur that, for a given roughness, one might consider an operational EM wavelength for which the roughness itself can be reasonably neglected. However, in principle, there is no defensible motivation, beyond the relevant limitation of the involved analytical difficulties, for considering the effect of only one interfacial roughness, neglecting the other relevant ones. This poses not only a conceptual limitation. In fact, in the applications perspective of retrieving geo-physical parameters by scattering measurements, whether there is a dominant interfacial roughness, and, in case, which the dominant one is, should be established after the remote sensing data are analyzed and, conversely, they cannot constitute a priori assumptions.

Each of the existing first-order models referring to a simplified geometry with one (Fuks, 2001) (Azadegan et al., 2003) (Sarabandi et al, 2000) (Yarovoy et al., 2000) or two (Tabatabaeenejad et al., 2006) rough interfaces, can be rigorously regarded as a particular case of our general model. Indeed, it can be analytically demonstrated that when the general geometry reduces to each simplified one, the consistency of the relevant solutions formally holds. In fact, when the (63)-(66) are specialized for the case of Figure 3a, the factors0|mp(k)exp(jkzmΔm)turn out to be unitary and the general solution formally reduces to the one in (Fuks, 2001). Similarly, by specializing the solution to the configuration of Figure 3b (Azadegan et al., 2003), the computation is reduced to onlyα˜qp1,2, in whichξ01±p(k)=T0|1p(k)[1+R0|1p(k)R1|2p(k)ej2kz1Δ1]1ejkz1Δ1[1±R1|2p(k)], so the equivalent solution in (Azadegan et al., 2003) (Franceschetti et al., 2008) is formally obtained. Finally, specializing the general solution to the geometry of Figure 3c, and considering the backscattering case (ki=ks), the computation is reduced to onlyα˜qp1,2in whichξ01±p(ki)=T0|1p(ki)[1+R0|1p(ki)1|2p(ki)ej2kz1iΔ1]1ejkz1Δ1[1±1|2p(ki)], so we formally obtain the equivalent solution in (Yarovoy et al., 2000) (Franceschetti et al., 2008).

Analytically speaking, this allows us to obtain, in a unitary formal framework, a comprehensive insight into the first-order perturbation solutions formalism for scattering from stratified structure with rough interfaces.

Finally, the Boundary Perturbation Theory results can be also regarded as a generalization of the classical SPM for rough surface (Ulaby et al., 1982) (Tsang et al., 1985) to layered media with rough interfaces.

Figure 3.

Simplified geometry considered by other Authors.


9. Wave Scattering Decomposition

In this section, the focus is on the intrinsic significance of the global BPT scattering solution, getting more concrete insight into the physics of the problem of the scattering from rough interfaces of a layered media. In order to be able to express the solution in terms of readable basic physical phenomena, a key point is to exploit the local scattering concept, differently from (Yarovoy et al., 2000) and (Fuks, 1998) wherein the authors resort to the radar contrast.

It should be noted that the exact analytic decomposition of the solution in terms of local interactions is rigorously feasible, since, in the first-order perturbative approximation, the scattering amplitude can be written as a single space integral with a kernel that depends only on the rough interface height and on its first-order derivatives at a given point. Moreover, since in the limit of first-order BPT solution the global response of a structure with all rough interfaces can be directly obtained considering the superposition of the response from each interface, we firstly focus our attention to a generic embedded rough interface. Afterwards, the general interpretation for a layered structure with an arbitrary number of rough interfaces can be addressed.

To focus formally on the relations among local and global scattering concepts, the identified Wave Scattering Decomposition (Imperatore et al 2008c) (Imperatore et al. 2009c), for the global scattering response of the structure in terms of the four types of local interactions, can be expressed with a compact notation as:




captures the local response of the mth rough interface between two layer of permittivityεm,εm+1respectively, and the transfer vectorPmqpis related to the coherent propagation inside the stratification (Imperatore et al. 2009c). Specifically, four distinct types of local interaction with an embedded rough interface can be distinguished: two of them identifiable as local scattering through the relevant interfacial roughness and other ones as local scattering from the roughness. We emphasize that the corresponding coefficientsαqpm,m+1andαqpm+1mrefer to cases in which both the observation and incidence directions are, respectively, above and under the roughness; whereasβqpm,m+1andβqpm+1mconcern the local scattering contribution that cross the roughness in opposite directions. In addition, we stress that the local scattering coefficients are formally identical to the classical ones relative to a rough surface between two half-spaces (Ulaby et al., 1982) (Tsang et al., 1985). On the other hand, we emphasize that the transfer vector, which measures the influence of the stratification on the local scattering, whatever the roughness is, can be expressed in terms of the generalized transmission/reflection coefficients (Imperatore et al. 2009c). Once the local nature has been recognized, the solution can be suitably expanded, so that it can be expressed as a ray series or optical geometric series, where each term of the series is susceptible of a powerful physical interpretation (Imperatore et al. 2008c) as illustrated in the next section.


10. Physical Interpretation

In this section, we show in detailed the physical meaning of the wave scattering decomposition obtained in the previous section. The analytical solution (69), after suitable expansion of the elements of the transfer vector, is then susceptible of an expression in terms of an infinite sum of contributions (geometric power series). Consequently, the suitably expanded solution can be expressed as an optical geometric series, where each term of the series is susceptible of a direct physical interpretation. In particular, each individual term of the absolutely summable infinite series can be physically identified as a wave propagating in the structure that experiences a single-scattering local interaction with the roughness.

To this purpose, we introduce the following useful notations


and recognize that these factors correspond to a complete roundtrip in the intermediate layer with coherent reflections at the layer boundaries. Moreover, in order to provide a symmetrical expansion, it is possible to explicit the factorMmh(ki), which is associated with the multiple round trip in the mth layer and is included in0|mh(ki); so we can write:


It should be noted that the0|mpare distinct from the coefficients0|mpslab, because in the evaluation of the former the effect of all the layers under the layer m is taken into account, whereas the latter is evaluated referring to a different configuration in which the intermediate layers 1...m are bounded by the half-spaces 0 and m.

Furthermore, we focus our attention on the two layers just above (m) and under (m+1) the generic roughness. In Figure 4, the remaining part of the structure is visualized condensed in two equivalent slabs constituted, respectively, by the intermediate layers m+2,...,N-1 (under the (m+1)th layer) and 1,...,m-1 (above the mth layer). Without loss of generality, since analogous considerations hold for the other polarization combinations, the analysis can be conducted for the hh case only. Consequently, four distinct families of rays can be recognized; each one associated to one type of local interaction, so that each term of the expansion of the (69) can be readily identified as follows (Imperatore et al. 2009c):

a) Local upward scattered waves from rough interface: each of these waves (see Figure 4.a) undergoes a coherent transmission into mth layer (0|mh(ki)exp(jkzmiΔm)through the intermediate layers 1,…,m-1, then j 1complete round-trips ([Λm(ki)]j1) in the mth layer with coherent reflections at the incident angle (kithen an incoherent local scattering from the rough interface (αhhmm+1(ks,ki)upward within the observation planekssubsequently j 2complete round-trips ([Λmh(ks)]j2in the mth layer with coherent reflections at the scattering angle (ksand finally a coherent transmission (exp(jkzmsΔm)m|0h(slab)(ks)in the upper half-space through the intermediate layers m-1,…,1.

b) Local upward scattered waves through rough interface: each of these waves (see Figure 4.b) undergoes a coherent transmission into the mth layer (0|mh(ki)exp(jkzmiΔm)through the intermediate layers 1,…,m-1, then j 1complete round-trips[Λm(ki)]j1in the mth layer with coherent reflections at the incident angle (kisubsequently a coherent transmission (Tm|m+1h(ki)) followed by n 1complete round-trips in the (m+1)th layer ([Λm+1h(ki)]n1and by further bounce on the (m+1)th flat interface (m+1|m+2h(ki)exp(j2kz(m+1)iΔm+1)at the incident angle (kiand after that an incoherent local scattering through the rough interface (βhhm+1,m(kski)upward within the observation plane (kssubsequently j 2complete round-trips ([Λmh(ks)]j2in the mth layer with coherent reflections at the scattering angle (ksand finally a coherent transmission (exp(jkzmsΔm)m|0h(slab)(ks)) in the upper half-space through the intermediate layers m-1,…,1.

c) Local downward scattered waves through rough interface: each of these waves (see Figure 4.c) undergoes a coherent transmission into the mth layer (0|mh(ki)exp(jkzmiΔm)), through the intermediate layers 1,…,m-1, then j 1complete round-trips ([Λm(ki)]j1) in the mth layer with coherent reflections at the incident angle (ki), then an incoherent local scattering through the rough interface (βhhm,m+1(ks,ki)) downward in the observation plane (ks) followed by further bounce on the (m+1)th flat interface (m+1|m+2h(ks)exp(j2kz(m+1)sΔm+1)) with subsequently n 2complete round-trips in the (m+1)th layer ([Λm+1h(ks)]n2) at the scattering angle (ks), next a coherent transmission (Tm+1|mh(ks)) followed by subsequently j 2complete round-trips ([Λmh(ks)]j2) in the mth layer with coherent reflections at the scattering angle (ks), and finally a coherent transmission (exp(jkzmsΔm)m|0h(slab)(ks)) in the upper half-space through the intermediate layers m-1,…,1.

d) Local downward scattered waves from rough interface: each of these waves (see Figure 4.d) undergoes a coherent transmission into the mth layer (0|mh(ki)exp(jkzmiΔm)), through the intermediate layers 1,…,m-1, then j 1complete round-trips ([Λm(ki)]j1) in the mth layer with coherent reflections at the incident angle (ki), next a coherent transmission (Tm|m+1h(ki)) followed by n 1complete round-trips in the (m+1)th layer ([Λm+1h(ki)]n1) and by further bounce on the (m+1)th flat interface (m+1|m+2h(ki)exp(j2kz(m+1)iΔm+1)) at the incident angle (ki), and after that an incoherent local scattering from the rough interface (αhhm+1,m(ks,ki)), downward in the observation plane (ks), followed by further bounce on the (m+1)th flat interface (m+1|m+2h(ks)exp(j2kz(m+1)sΔm+1)) with subsequently n 2complete round-trips in the (m+1)th layer ([Λm+1h(ks)]n2) at the scattering angle (ks), then a coherent transmission (Tm+1|mh(ks)) followed by subsequently j 2complete round-trips ([Λmh(ks)]j2) in the mth layer with coherent reflections at the scattering angle (ks), and finally a coherent transmission (exp(jkzmsΔm)m|0h(slab)(ks)) in the upper half-space through the intermediate layers m-1,…,1.

Note also that when an arbitrary layered structure with all rough interfaces is concerned, since in the first-order limit the multiple scattering contributions are neglected, the relative physical interpretation can be obtained effortless by superposition of the several ray contributions obtained considering separately each rough interface.

Figure 4.

Physical interpretation for the scattering from an arbitrary layered structure with an embedded rough interface.

Figure 5.

Physical Interpretation: bistatic configuration.

Figure 6.

Physical Interpretation: mono-static configuration.

As a result, the obtained interpretation (Figure 4) enables the global scattering phenomenon involved to be visualized as a superposition of local interactions, emphasizing the role of the interference effects in the structure as well (Imperatore et al. 2009c). It should be also noted that, despite the expansion is attained rigorously without any further approximation with respect to the solution proposed (see (54)-(58)), the resulting interpretation turns out to be extremely intuitive and surprisingly simple. In particular, when the configuration reduces to a rough interface covered by a dielectric layer, as the reader can easily verify, we obtain the interpretation (Franceschetti et al., 2008) for the bistatic and monostatic configuration illustrated in Figure 5 and Figure 6, respectively.

11. Scattering patterns computation

In this section, we present some numerical examples aimed at studying scattering coefficients (59). To this purpose, we consider the canonical layered media with three rough interfaces pictured in Figure 8, which is representative of several situations of interest. In common with classical theoretical studies of the scattering of waves from random surfaces, we assume that the interfaces constitute Gaussian 2D random processes with Gaussian correlations, whose spectral representation is given by


where, with regard to the nth interface, nand l nare the surface height standard deviation and correlation length, respectively. In order to perform a consistent comparison, we refer to interfaces with the same roughness. In addition, we suppose no correlation between the interfaces. For instance, we analyze the layered medium with three rough interfaces schematized in Figure 7, which can be parametric characterized as follows. We assume k 0l n=1.5, k 0σ n=0.15 for n= 0, 1, 2. In addition, the considered vertical profile is characterized by the following parameters: 0=1, 1=3.0+j0.0, 2=5.5+j0.00055, 3=10.5+j1.55; Δ 1/=1.50, Δ 2/=2.80 (see Table 1). Once this reference structure has been characterized, we study the scattering cross section of the structure as a function of the scattering direction in the upper half-space, assuming fixed the incident direction. It should be noted that, also considering a limited number of layers, the number of parameters involved by the model makes difficult the jointly visualization of the multi-variables dependency. As matter of fact, once the structure has been parametrically defined and incident direction has been fixed, it is possible to visualize the scattering cross section of the structure as a function of the scattering direction in the upper half-space (Imperatore et al. 2008c). Therefore, to characterize the re-irradiation pattern of the structure in three-dimensional space, scattering cross-section distributions are represented (Figure 8) as function of zenithal and azimuthal angles and are treated as three-dimensional surfaces. To save space, only the case hh is considered. In addition, we assume fixed the incidence angleθ0i= 45 (k^i=x^). Therefore, to evaluate the effect on the global response of each rough interface, the several single contributions are shown in Figure 8a, Figure 8b, and Figure 8c, respectively. In addition, the total contribution is also pictured (Figure 8d). It should be noted that to visualize the patterns an offset of +40dB has been considered for the radial amplitude, so that scattering coefficients less than -40dB are represented by the axes origin.

Figure 7.

Three rough interfaces layered media.

Figure 8.

Bi-static scattering coefficients hh for a three rough interfaces layered media: 0contribution (a), 1contribution (b), 2contribution (c), total contribution (d) (note that scattering coefficients values less than -40 dB are represented by the axes origin).

f1.0 GHzk0l01.5
25.5+j 0.00055k0l21.5
310.5+j 1.55

Table 1.

Parametric characterization of the layered media of Figure 7.

12. Conclusion

A quantitative mathematical analysis of wave propagation in three-dimensional layered rough media is fundamental in understanding intriguing scattering phenomena in such structures, especially in the perspective of remote sensing applications. The results of the Boundary Perturbation Theory (BPT), as introduced by P. Imperatore and his coauthors in many different papers, essentially constitutes the content of this chapter in which the theoretical body of results is presented in organized manner, emphasizing the applications perspective. These formally symmetric and physically revealing analytical results are crucial and will contribute to innovatory applications in microwave remote sensing. For instance, they open the way toward new techniques for solving the inverse problem, for designing SAR processing algorithms, and for modelling the time-domain response of complex layered structures.


  1. 1. FranceschettiG.ImperatoreP.IodiceA.RiccioD.RuelloG.2003Scattering from a Layered Medium with One Rough Interface: Comparison and Physical Interpretation of Different Methods, Proceedings of IEEE International Geoscience and Remote Sensing Symposium (IGARSS 2003), Toulouse, France,29122914, Jul. 2003.
  2. 2. FranceschettiG.ImperatoreP.IodiceA.RiccioD.RuelloG.2008Scattering form Layered Strucutures with One rough Interface: A Unified Formulation of Perturbative Solutions.IEEE Trans. on Geoscience and Remote Sensing,466June 2008,16341643.
  3. 3. ImperatoreP.IodiceA.RiccioD.2008aPerturbative Solution for the Scattering from Multilayered Structure with Rough Boundaries,Proceedingsof Microwave Radiometry and Remote Sensing of the Environment (MICRORAD 2008),14, 11-14 March 2008, Florence, Italy.
  4. 4. ImperatoreP.IodiceA.RiccioD.2008bSmall Perturbation Method for Scattering From Rough Multilayers, Proceedings of IEEE International Geoscience and Remote Sensing Symposium, 2008 (IGARSS 2008),5pp.V-271-V-274, Boston 7-11 July 2008.
  5. 5. ImperatoreP.IodiceA.RiccioD.2008cInterpretation of Perturbative Solution for the Scattering From Layered Structure With Rough Interfaces, Proceedings of IEEE International Geoscience and Remote Sensing Symposium, 2008 (IGARSS 2008),4pp. IV-1141-IV-1144, Boston 7-11 July 2008.
  6. 6. ImperatoreP.IodiceA.RiccioD.2009aElectromagnetic Wave Scattering From Layered Structures with an Arbitrary Number of Rough Interfaces.IEEE Trans. on Geoscience and Remote Sensing,47410561072, April 2009.
  7. 7. ImperatoreP.IodiceA.RiccioD.2009bTransmission Through Layered Media With Rough Boundaries: First-Order Perturbative Solution.IEEE Transaction on Antennas and Propagation,57514811494, May 2009.
  8. 8. ImperatoreP.IodiceA.RiccioD.2009cPhysical Meaning of Perturbative Solutions for Scattering from and Through Multilayered Structure with Rough Interfaces.submitted to IEEE Trans. on Antennas and Propagation.
  9. 9. SarabandiK.AzadeganR.2000Simulation of interferometric SAR response to rough surfaces covered with a dielectric layer, Proceedings IEEE IGARSS, 24-28 July 2000173917414
  10. 10. AzadeganR.SarabandiK.2003Analytical formulation of the scattering by a slightly rough dielectric boundary covered with a homogeneous dielectric layer, Proceedings of IEEE AP-S Int. Symp., Columbus, OH, Jun. 2003,420423.
  11. 11. FuksI. M.1998Radar contrast polarization dependence on subsurface sensing, Proceedings of IEEE IGARSS’98,3Seattle, WA, USA, July 6-10, 1998,14551459.
  12. 12. FuksI.M.2001Wave diffraction by a rough boundary of an arbitrary plane-layered medium.IEEE Transaction on Antennas and Propagation,494Apr. 2001.630639.
  13. 13. YarovoyA. JonghR. V.LigthardL. P.2000Scattering properties of a statistically rough interface inside a multilayered medium.Radio Science,35n.2, 2000,445462.
  14. 14. TabatabaeenejadA.MoghaddamM.2006Bistatic scattering from three-dimensional layered rough surfaces.IEEE Trans. on Geoscience and Remote Sensing,44821022114, Aug. 2006.
  15. 15. ChewW. C.1997Waves and Fields in Inhomogeneous Media. IEEE Press, 1997.
  16. 16. IshimaruA.1978Wave Propagation and Scattering in Random Media,2Academic Press, New York, 1978.
  17. 17. TsangL.KongJ. A.ShinR. T.1985Theory of Microwave Remote Sensing. New York: Wiley, 1985.
  18. 18. UlabyF. T.MooreR. K.FungA. K.1982Microwave Remote Sensing, vol. I, II, II. Reading, MA: Addison-Wesley, 1982.

Written By

Pasquale Imperatore, Antonio Iodice and Daniele Riccio

Published: February 1st, 2010