When a beam emitted from an active monostatic sensor system sweeps across a volume, the echoes from scatterers present will fluctuate from ping to ping due to various interference phenomena and statistical processes. Observations of these fluctuations can be used, in combination with models, to infer properties of the scatterers such as numerical density. Modeling the fluctuations can also help predict system performance and associated uncertainties in expected echoes. This tutorial focuses on “physics-based statistics,” which is a predictive form of modeling the fluctuations. The modeling is based principally on the physics of the scattering by individual scatterers, addition of echoes from randomized multiple scatterers, system effects involving the beampattern and signal type, and signal theory including matched filter processing. Some consideration is also given to environment-specific effects such as the presence of boundaries and heterogeneities in the medium. Although the modeling was inspired by applications of sonar in the field of underwater acoustics, the material is presented in a general form, and involving only scalar fields. Therefore, it is broadly applicable to other areas such as medical ultrasound, non-destructive acoustic testing, in-air acoustics, as well as radar and lasers.

## LINKS TO SOFTWARE USED TO CREATE FIGURES

The software used to produce all predictions in the figures is in the supplementary material at https://doi.org/10.1121/1.5052255. The software is also stored online (Lee and Baik, 2018), where it is subject to future revisions.

## I. INTRODUCTION

Echoes, as measured through the receiver of an active monostatic sensor system, will typically fluctuate from ping to ping as the beam emitted from the system scans across a volume containing scatterers or as the scatterers in that volume move through the beam (Fig. 1). It is essential to understand the echo statistics for accurate interpretation of the scattering and for modeling system performance. Toward that goal, understanding the underlying physical processes that give rise to the fluctuations allows one to accurately predict and interpret the echo statistics. In the simplest case in which the propagation medium is homogeneous and there are no boundaries present (i.e., a “direct path” geometry), the fluctuations are due to a combination of several statistical processes: the interference between overlapping echoes when multiple scatterers are present, the random nature of the echoes from individual scatterers (not including beampattern effects), and the modulation of the echo due to the random location of the scatterer in the beam. Once the geometry is further complicated by heterogeneities in the medium and/or the presence of boundaries, the propagated signals will become refracted and/or rescattered, giving rise to more propagation paths that are potentially random and, in turn, also contributing to the fluctuations. This tutorial presents key concepts and formulations associated with predicting the echo fluctuations over a wide range of scenarios in terms of the physics of the scattering, system parameters, and signal theory.

The statistical behavior of echoes is important across a diverse range of active sensor systems and applications involving the use of either acoustic waves (such as with sonar, medical ultrasonics, or non-destructive testing) or electromagnetic waves (such as with radar or light) to study individual discrete scatterers, assemblages of scatterers, or rough interfaces that cause scattering. Understanding echo statistics has been integral in interpreting radar clutter (Watts and Ward, 2010) and sonar reverberation and clutter (Ol'shevskii, 1978; Gallaudet and de Moustier, 2003; Abraham and Lyons, 2010), sonar classification of marine life and objects on the seafloor (Stanton and Clay, 1986; Medwin and Clay, 1998), medical ultrasound classification of human tissue (Eltoft, 2006; Destrempes and Cloutier, 2010, 2013; Oelze and Mamou, 2016), and non-destructive ultrasound testing of materials (Li , 1992). Within each of these areas, echo statistics are used in the detection and classification of scatterers, discriminating between scatterers of interest from clutter (i.e., unwanted echoes), estimating numerical density of scatterers, and determining the performance of sensor systems for use in detection and classification of scatterers. The studies of clutter involve characterizing the statistics of unwanted echoes from features or objects in the environment that have properties resembling the target of interest. The clutter may be due to the sea surface (radar/sonar/laser), marine life or seafloor (sonar), human tissue (medical ultrasound), or grains in solids (non-destructive testing).

The above applications have many elements of statistical theory in common. Those common elements are treated formally in Goodman (1985) and Jakeman and Ridley (2006). There are also notable differences between interpreting echoes from acoustic and electromagnetic systems, such as the presence of shear waves and polarization, respectively, which are summarized in Le Chevalier (2002). Differences in echo statistics between scalar fields (both acoustics and electromagnetic, ignoring shear waves and polarization, respectively) and those fields with polarization effects (electromagnetic only) are summarized in Jakeman and Ridley (2006).

While there has been much work conducted in the area of echo statistics, the focus has generally involved describing the variability of echoes through use of generic statistical functions whose parameters need to be determined from experimental data. Here, “generic” refers to those functions generally devoid of a physical basis and derived solely from random signal theory, such as a Gaussian process. These generic functions include the Rayleigh, Rice, K, Weibull, log normal, and Nakagami-m distributions (App). Only under a narrow range of conditions can some of the generic functions be connected to the physics. There has been a much smaller set of studies in which the statistics were modeled from first principles, and hence being predictive, through a combination of the physics of scattering, sensor system parameters, and signal theory. To date, there has been no text or tutorial with a comprehensive coverage of physics-based phenomena.

This tutorial presents the fundamentals of echo statistics associated with an active monostatic sensor system, written in a manner that is readable by a wide audience. While much of the material is motivated by the use of sonar in underwater applications, the material is broadly applicable to other types of active acoustic systems used in various media including air, solids, and biological tissue, as well as to radar and lasers. Because of the generality of the treatment, the system used to transmit waves and receive the echoes will heretofore be referred to as a “sensor system” regardless of application and whether it involves acoustic or electromagnetic waves. The formulations describe scalar fields, which can be directly applied to either acoustic or electromagnetic signals in the absence of the presence of shear waves or polarization effects, respectively. Only brief reference is made to those latter effects.

The emphasis in this tutorial is on describing the statistics of the magnitude of the complex echo (or echo envelope) in terms of the physics of the scattering, sensor system parameters, and signal processing. More specifically, the emphasis is on *physics-based* methods, in contrast to the generic approaches referenced above. The methods are focused solely on “first-order statistics” which concern statistical properties of the signal at a single instant in time. Also, the material can apply to single beam (fixed or scanning) and multi-beam systems where the scatterer(s) in each case are randomly located in each beam.

Since there are too many combinations of types of systems, signal processing, and environments to adequately describe in a single paper, the material is focused on the fundamental aspects of echo statistics due to discrete scatterers that are not specific to any particular system, signal processing, or environment. The range of beamwidths, types of scatterers, and number of scatterers in the predictions illustrate the corresponding range of statistical behavior of the echoes from a wide range of system parameters and distributions of scatterers. However, the formulations are completely limited to (1) first-order statistics as described above, and are mostly limited to (2) narrowband signals that are long enough so that the echoes from the scatterers overlap significantly, (3) direct path geometries where there are no boundaries present, and (4) a homogeneous medium. Beyond those limited scenarios, examples are presented of more complex cases involving pulsed signals in which the echoes from the scatterers only partially overlap and the presence of boundaries and/or heterogeneities, including waveguide effects. There are also many more important cases not covered explicitly, but that can be described using these formulations as a basis, including advanced signal processing and complex environments.

A common theme of the material involves predicting the degree to which the statistics of the echo magnitude deviate from the classic Rayleigh distribution (described in more detail later), with a focus on the “tail” of the echo distributions. Experimental observations of echo statistics, particularly in the limit of a large number of scatterers, tend toward the Rayleigh distribution and deviations from that distribution contain information on the scatterers. The tail corresponds to the highest values of echo magnitude and will typically be the part of the echo that is detected above the background noise or reverberation. This tail is shown to contain valuable information as it is a function of the numerical density of scatterers, type of scatterer, bandwidth, and beamwidth.

The tutorial is organized as follows: Secs. II–IV provide qualitative descriptions and illustrations of concepts that are important to echo statistics, a summary of the range of applications that exploit echo statistics, and equations specific to scattering and random processes (not specific to scattering) that are common to many formulations in echo statistics. Sections V–VIII draw from the material in Secs. II–IV in a presentation of echo statistics formulations for a wide range of important physical scenarios. These scenarios include beampattern effects associated with main lobes of various width, narrowband and broadband signals, completely overlapping echoes (long signal) and partially overlapping echoes (short signal), single scatterers and mixed assemblages of scatterers, elongated and randomized scatterers, and geometries involving either a direct path and a homogeneous medium or ones involving the presence of boundaries and/or heterogeneities. Sections VI and VII include a progression of cases, all involving a direct path geometry with a homogeneous medium using narrowband signals that are long enough so that the echoes from the scatterers overlap significantly. The two sections begin with the case involving no beampattern effects (omni-directional beam) (Sec. VI), and then later incorporating beampattern effects (directional beam) (Sec. VII). Section V, besides providing an overview of the material of the rest of the tutorial, outlines how the results in Secs. VI and VII can be extended to more complex cases, such as the ones given in Sec. VIII. Section VIII draws from Secs. VI and VII for selected complex cases involving partially overlapping echoes and the presence of boundaries and/or heterogeneities. Although many details of the derivations are not given in Secs. VI–VIII, references to the concepts and equations in Secs. II–IV are given, as well as references to previously published papers.

## II. KEY ELEMENTS OF ECHO STATISTICS

Fundamental to the echo statistics are the (1) various types of interference between overlapping echoes from multiple scattering features from a single scatterer, multiple scatterers, and/or multiple propagation paths in a heterogeneous medium or a medium containing boundaries, (2) stochastic nature of scattering by individual scatterers (not including sensor system effects), and (3) random modulation due to random location of the scatterers in the sensor beampattern (Fig. 1). These elements are first discussed qualitatively with important aspects illustrated. Formal mathematical treatments based on the physical processes, sensor system, and signal types are presented later. All processes are described quantitatively in terms of the probability density function (PDF), which is the probability of occurrence of a random variable for each of its values. The PDF and related quantities are defined formally in Sec. IV B 1.

### A. Interfering signals

Echoes from one or more scatterers can be decomposed into the sum of the echoes from multiple scattering highlights from each individual scatterer, echoes from multiple scatterers when more than one scatterer is present, and echoes from multiple propagation paths when the medium is heterogeneous and/or there are one or more boundaries present. Regardless of scenario, because of randomness generally associated with the scatterer or environment, the echoes will tend to interfere randomly with each other, causing fluctuations from ping to ping.

In order to understand this interference phenomenon, we first examine the simple case of a single sinusoidal signal whose amplitude (i.e., envelope) does not vary in time. This signal corresponds to the echo from a point scatterer fixed in the beam of a narrowband system. In this case, the amplitude of the signal is constant [Fig. 2(a)]. The PDF describing this single-valued quantity is the delta function [Fig. 2(a)]. Note that when expressing the signal as a complex variable (which is the case for most of this tutorial), then it is the *magnitude* of the signal that is of interest which is equal to the amplitude of the sine wave in this example.

The statistics change dramatically once another signal with the same magnitude but of random phase is added. This second signal corresponds to the echo from a second identical scatterer (or scattering feature or ray path) whose echo overlaps with that of the first. Now, the PDF of the magnitude of the sum is spread over a range of values, depending on whether the signals added constructively, destructively, or something in between [Fig. 2(b)]. As more signals of equal magnitude and random phase are added, the PDF continues to change and, in the limit of an infinite number of signals of random phase, the PDF converges to the Rayleigh PDF [Fig. 2(c)].

This example involved single-frequency signals of infinite extent. This idealized signal can be used to approximate pulsed (or “gated”) sine wave signals used in many sensor systems where echoes from scatterers significantly overlap. The same principles also apply, although with more complexity, to pulsed signals (narrowband and broadband) of much shorter extent where the multiple signals generally overlap only partially as presented in Sec. VIII A.

### B. Stochastic scattering from a single scatterer

The simplest of scatterers is a point scatterer, as it gives the same echo regardless of orientation. In this case in which effects from the sensor system are ignored (i.e., omnidirectional beam) and when using a sinusoidal signal, the echo magnitude from a single point scatterer is constant with a delta function PDF [Fig. 3(a)]. For more realistic scatterers of finite extent, such as elongated ones whose orientation and/or shape may change in time, the echo magnitude will vary from ping to ping. Now, the PDF is significantly broadened and the echo statistics are more complex [Fig. 3(b)].

There is some correspondence between this and the former case involving sinusoidal signals. The echo from a point scatterer corresponds to a single sine wave with a constant magnitude. In some cases, an echo from an irregular finite-sized scatterer can be decomposed into echoes from different parts of the body of the scatterer, which corresponds to the addition of multiple random-phased sinusoidal signals. Also, there are conditions under which the echo PDF from a complex scatterer can be the Rayleigh PDF, which corresponds to the sum of an infinite number of randomly phased sine waves.

### C. Beampattern effects

The echo measured through the receiver of the sensor system depends, in part, on the location of the scatterer in the two-way beampattern. The beampattern will modulate the echo according to the location, with stronger values being associated with the scatterer being near the center of the main lobe and weaker values corresponding to locations well away from the center. For a randomly located scatterer, the modulation will correspondingly be random, adding to the variability of the echo. The distribution associated with this random modulation is referred to as the “beampattern PDF.”

The randomizing effects of the beampattern on the echo from a randomly located object can be illustrated by varying the solid angle within which the object is allowed to move (Fig. 4). In this simple example, a circular piston transducer is used to both send and receive the signal and the solid angle is centered about the axis of symmetry of the transducer (which is the center of the main lobe). A point scatterer is allowed to move randomly within the solid angle at a fixed distance from the transducer. In the first case, the solid angle is 0 and the scatterer is restricted to remain fixed in the center of the beam where the response is maximum [Fig. 4(a)]. The corresponding echo is single valued from ping to ping and its PDF is the delta function [right panel of Fig. 4(a)]. As the solid angle increases from zero allowing the scatterer to randomly move a greater amount across the beam, the ping-to-ping variability in the echo increases correspondingly [Figs. 4(b)–4(d)]. When only the portion of the main lobe above the highest sidelobe is involved in the motion [Figs. 4(a)–4(c)], then the echo PDF is either single valued or approximately power-law-distributed. However, once sidelobes are involved, then the echo PDF is more complex, as a non-monotonic characteristic is present due to the sidelobe structure [“spiky” section of curve in Fig. 4(d), where values of echo are lower].

These effects will be shown in Sec. VII A 4 to be qualitatively similar over a wide range of beamwidths. In essence, regardless of beamwidth (several degrees or several tens of degrees), if the scatterer is randomly located within a solid angle containing at least most of the main lobe of the beampattern, then the echo will correspondingly be modulated over a wide range of values of the beampattern.

## III. EXPLOITING ECHO STATISTICS FOR VARIOUS APPLICATIONS

Understanding and quantifying the statistical variability of a signal is useful in a diverse range of applications. Below are described two common uses of echo statistics, one in which the variability of the observed signal is used to infer important information regarding the scatterers, and the other in which the degree to which the signal varies for a given scenario is predicted to understand the error or uncertainty regarding the expected signal. Specific examples are given in Sec. III A where echo statistics are used as an inference tool spanning use of sonar, radar, and medical ultrasound.

### A. Inferring information on scatterers

Regardless of the field and type of sensor system, one can exploit properties of the variability of the echo to make inferences of important quantities such as number of scatterers, scatterer characteristics, discriminating between and classifying different types of scatterers, and probability of false alarm when detecting a scatterer of interest that is interspersed with other unwanted scatterers. Key to the success of the inference is using a model of the echo statistics whose parameters can be linked to the physical properties of interest. The description of such models is the focus of this tutorial.

Once the statistics model is determined for the particular application, there are numerous methods to infer parameters of modeled echo statistics from data including use of least squares, maximum likelihood estimators (MLE) (Azzalini, 1996), and method of moment estimators (MME) (Joughin , 1993). Performance of an unbiased estimator can be evaluated formally through the Cramer-Rao lower bound (CRLB) method which estimates a lower bound of the variance of the estimator (Hogg and Craig, 1978). Abraham and Lyons (2002b) provides an example of use of MLE, MME, and CRLB for inference of modeled parameters from data (Appendix B of that paper).

Many applications using echo statistics as an inference tool are described in citations given in Sec. I. Selected examples spanning three types of sensor systems—sonar, radar, and medical ultrasound—are briefly summarized below.

#### 1. Inferring number of scatterers

The shape of the echo PDF can be used to infer the number of scatterers. This topic has especially seen much attention spanning many applications, with attempts to connect parameters of either physics-based or generic PDFs, such as the K PDF, to numerical density of scatterers. Stanton (2015) used physics-based methods (such as described in this tutorial in Sec. VII B) to relate the shape of the PDF of the magnitude of the echo to the number of scatterers in a sonar beam in the ocean. The inferences were conducted for several different types of scatterers (“bottom-like,” “compact stationary,” and “compact non-stationary”) that ranged in densities from resolved to unresolved. Lee and Stanton (2015) used similar physics-based methods with a sonar beam to infer the numerical density of fish (more details on this broadband method given in Sec. VIII A). Other studies have related a parameter of a generic PDF to numerical density of scatterers. For example, Abraham and Lyons (2002b) analytically related the shape parameter of the K PDF to numerical density of scatterers (for certain conditions; general and not specific to sonar) and applied the results to using sonar to estimate numerical density of scatterers on the seafloor. Tunis (2005) performed a controlled laboratory experiment with various dilute solutions containing cancerous cells (acute myeloid leukemia and prostate adenocarcinoma) to empirically relate a parameter of the gamma distribution of the medical ultrasound echo to numerical density of the cells.

For the general case in which there is a mixture of different sized scatterers, Lee and Stanton (2014) have formulated the echo PDF based on various physical parameters—sensor beampattern, scattering amplitudes, and numerical density of each type of scatterer. This method is not specific to any type of sensor system (and is described in Sec. VII C) and the simulations demonstrated conditions under which the number of scatterers (especially within the type that dominated the echo) could be inferred (see Figs. 6 and 7 of that paper).

#### 2. Discriminating between echo from scatterer of interest and background

For a single scatterer of interest interspersed within an aggregation of other unwanted scatterers (in a volume or on a surface), echoes from the scatterer of interest can be discriminated from those containing only the surrounding unwanted or “background” scatterers. Here, the echo containing the scatterer of interest will also be contaminated with echoes from the background. Ferrara (2011) conducted an experimental study on the ocean involving use of radar to detect and classify echoes from ships and oil rigs. In their data, there were echoes from two types of regions containing (1) both the scatterers of interest (ships or oil rigs; one at a time) and the sea surface and (2) only the sea surface. The experimental radar-echo statistics data were fit to the generalized K PDF. They demonstrated that combinations of generalized K parameters (such as ratios of the parameters) could be used to unambiguously distinguish between echoes involving the scatterers of interest (ships and oil rigs) and the background (sea surface) scattering from those due to the surrounding sea surface alone. There was only one false alarm out of 229 detected targets.

#### 3. Removing beampattern effects to isolate properties of resolved scatterer

While the previous case in Sec. III A 2 involves classifying a scatterer of interest when its echo is confounded with that of the surrounding background, the case simplifies once the echo is resolved from the background. This can occur when the scatterer is suspended in the volume away from any boundary or other neighboring scatterer. In this case, through inversion methods, variability due to the scatterer's random location in the beampattern can be removed to isolate echo variability due to the target alone. The target echo variability (with no beampattern effects) contains information on the scatterer, such as its size and degree to which it is elongated or rough.

For example, in Clay (1983), the PDF of the echo (as described in Sec. VII A 1 of this tutorial), which is a function of both the beampattern PDF and the PDF of the scattering amplitude of the scatterer, is formulated in terms of a convolution involving each of those individual PDFs. Through knowledge of the beampattern properties, the effects of the beampattern are removed from the echoes through deconvolution, leaving only the PDF of the scattering amplitude. The statistics of the scattering amplitude are then used to classify the scatterer in terms of its size and type. Although the paper was intended for use of sonar to detect and classify echoes from fish, the equations are general and applicable to any sensor system. Clay's method was refined in Stanton and Clay (1986).

#### 4. Discriminating between different types of aggregations of scatterers

Oelze and Mamou (2016) reviewed a number of previously published studies concerning use of medical ultrasound to detect and classify *in vivo* soft tissue. One such study (Fig. 3 of that paper) combined the spectral content of the broadband echo [to estimate effective scatterer diameter (ESD)] with echo statistics data fit to the homodyned K distribution. The feature analysis plot of ESD versus two parameters of the homodyned K distribution provided strong discrimination (with no overlap) between echoes from three different types of cancerous tumors (mammary) in rodents (rat fibroadenomas, mouse carcinomas, and mouse sarcomas).

#### 5. Further considerations

In order for the echo statistics to best be exploited, there are further considerations in using echo statistics as a tool with a practical sensor system.

##### a. Directionality of sensor beam improves inference techniques.

The degree to which the echoes are non-Rayleigh contains valuable information such as numbers of scatterers that can be inferred, as summarized above. Once the echoes are Rayleigh distributed, the amount of information to be inferred is limited. The beampattern causes echoes that may be otherwise Rayleigh-like to be non-Rayleigh through reducing the effective number of scatterers that are “seen” by the narrow mainlobe of the system. Furthermore, echoes that are non-Rayleigh before beampattern effects will deviate even more from the Rayleigh PDF once the beampattern is accounted for. The process of inferring information from the field of scatterers can therefore be enhanced or even optimized through varying, when possible, the width of the main lobe of the beampattern so that the echoes will be non-Rayleigh over the range of expected values.

##### b. Noise and “tail” of echo PDF.

In any practical sensor system, there will be various sources of noise (including electronic noise) and background reverberation that will normally dominate the lower values of echoes. Thus, any information to be inferred must generally come from the higher portion of the echo PDF—that is, the “tail.” Many analyses focus solely on properties of the tail that is above the detection threshold of the system.

### B. Predicting error or uncertainty in signal magnitude

In many applications, it is desired to predict the magnitude of the signal after it has propagated through an environment. The applications can deal with ones such as those described above in which one wishes to infer properties of the environment or scatterers from the signal or other types of applications in which the performance of the sensor system is being studied. Ideally, one wants to know all key parameters of the environment and scatterers so that the properties of the signal that propagates through the environment can be known to 100% certainty. However, due to a combination of the uncertainties of the values of the many parameters of the environment [e.g., roughness of surface(s), material properties, etc.] and scatterers (e.g., number, type, orientation, etc.) as well as random variations of those same quantities due to naturally occurring processes, the predicted signal in a natural environment is normally a random variable.

To quantify the degree to which the signal varies (that is, predict what is commonly called its error or uncertainty), the random variables associated with parameters of the environment and scatterers must first be estimated. Using these random variables, the variability of the signal can then be predicted through use of physics-based equations as described below. Regardless of the approach used, there are different degrees of variability. First, predicting the variability of the signal due to random variations in the environment and scatterers from the expected mean properties. A higher order variability would be due to errors in those estimates of the random variations in the environment and scatterers and associated deviations from the predicted variability of the signal (that is, an error in a prediction of a signal PDF).

Beyond the above considerations in which there is no noise in the system, there is also uncertainty in the predictions associated with random noise. This noise can be due to a combination of ambient noise in the environment and electrical system noise. Also, diffuse background reverberation is sometimes considered to be part of the noise. Although this topic is outside the scope of this tutorial, it is an important consideration in modeling practical systems. The noise can be accounted for through various methods which include (1) modeling the noise as an additive random complex variable to the (complex) echo signal (Jones , 2017) and (2) modeling the echo magnitude PDF as a “mixture” PDF, which is the sum of the noise-free echo magnitude PDF and the noise-only PDF (Stanton and Chu, 2010; Abraham , 2011).

## IV. KEY EQUATIONS FOR SCATTERING AND ASSOCIATED STATISTICAL PROCESSES

Fundamental equations are first given describing deterministic scattering such as for a fixed location or orientation of the scatterer. In order to describe the scattering for a randomized case such as when the location or orientation vary randomly, general equations describing statistical processes that occur in such scenarios are given, which are then incorporated into the equations for deterministic scattering to describe randomized scattering. All of the below equations relate the echo and its fluctuations directly to the physical properties of the scatterers and sensor system. This type of treatment is referred to as “physics based.”

### A. Deterministic scattering

The scattering described in this section is from a single ping or statistical realization. This deterministic description will then be randomized for use in predicting the statistics of the echoes. The scattering geometry involves use of an active monostatic sensor system where the transmitter and receiver are collocated (Fig. 1). The scattering can involve one or multiple scatterers. Also, all analyses involve scatterers distributed within a thin shell of constant radius with the sensor system at the center. The thickness of the shell is much smaller than its radius.

This (shell) geometry is for modeling a long narrowband signal in which echoes from all of the scatterers within the shell completely overlap. For a sensor system that uses short pulses to detect a more broadly distributed set of scatterers (such as in many practical applications), the thickness of this thin shell can be related approximately to the duration of the pulse. A rigorous approach to modeling pulsed systems where the echoes generally only partially overlap is given in Sec. VIII A.

#### 1. Single scatterer

*V*received by the sensor system due to the echo from a single scatterer is

_{s}*j*= $ \u2212 1$, $\omega $ is the angular frequency of the sinusoidal signal, $k$ is the acoustic or electromagnetic wavenumber ( = 2 $\pi $/

*λ*, where

*λ*is the wavelength),

*α*is the absorption coefficient of the medium so that $ e \u2212 2 \alpha r$ is the two-way loss due to absorption, and $ b ( \theta , \varphi )$ is the two-way beampattern of the sensor system whose values lie in the range [0,1]. The term $ b ( \theta , \varphi )$ is the product of the beampatterns of the transmitter and receiver and the terms $\theta $ and $ \u2009 \varphi $ are the angular coordinates of the scatterer. Specifically, $ b ( \theta , \varphi ) = b T ( \theta , \varphi ) b r ( \theta , \varphi )$, where $ b T$ and $ b r$ are the transmitter and receiver beampatterns, respectively. The term $ f b s$ is the backscattering amplitude of the scatterer and is a complex variable.

In this formulation, the signal is assumed to be at a single frequency (i.e., narrowband) of infinite temporal extent. Also, the acoustic and electromagnetic fields associated with $ G T$ and $ G R$ are assumed to be scalar quantities. For acoustics, this scalar quantity is pressure, the compressional component of the field, and assumes no shear component in the medium (although conversion of compression to shear wave can take place within the scatterer). For the electromagnetic field, the scalar quantity is one polarization component of the electric or magnetic field. This assumption in this latter case treats each polarization component independently and ignores coupling between the components [Chap. 4 of Goodman (1985); Chap. 4 of Jakeman and Ridley (2006)].

^{2}) are suppressed. Note that the term $ \sigma b s$ should not be confused with a similar representation for backscattering cross section, $ \sigma \u2009$, that is commonly used where $ \sigma = 4 \pi \sigma b s$ and TS = $ 10 \u2009 log ( \sigma \u2009 / 4 \pi ) .$

*α*in Eq. (1) have been suppressed and $ r 0 = 1$

*m*. Here, the magnitude of the signal is based simply on the absolute value of the signal.

#### 2. Multiple scatterers

*N*scatterers is

*i*th scatterer, respectively. As with Eq. (1), the signal is assumed to be at a single frequency of infinite temporal extent. With this type of signal, there is 100% overlap between the echoes from all individual scatterers. The simple summation of echoes from individuals reflects the assumption that only single-order scattering is being considered and higher-order scattering (i.e., re-scattering of echoes between individuals) is assumed to be negligible.

*N*scatterers as received through the sensor system is given by

*i*th scatterer as received through the sensor system is

*i*th scatterer are in the term $ \Delta i$, which includes phase shifts due to differences in range within the shell (i.e.,

*2kr*) and due to the scatterer and beamformer. For acoustic/electromagnetic wavelengths that are small compared with the shell thickness, $ \Delta i$ will generally vary randomly and uniformly in the range [0 2 $\pi $] for randomly distributed scatterers.

_{i}### B. Probability density function and related quantities

#### 1. Definitions and equations

The statistics described herein involve the probability of occurrence of random variables. This is in contrast to other types of statistics such as statistical tests (e.g., the t- and Kolmogorov-Smirnov tests). The principal statistical quantities used are the probability density function (PDF or *p*), cumulative distribution function (CDF), and probability of false alarm (PFA), which are interrelated from the following expressions for one-dimensional continuous random variables (Ol'shevskii, 1978; Papoulis, 1991; and Goodman, 1985). Expressions involving multi-dimensional random variables will appear later in context. While the below expressions are general, the integration limits associated with echo magnitude statistics later will reflect the fact that the magnitude (*x*) is always positive and $ p = 0$ for $ x < 0$.

Note that the CDF is rigorously referred to as the “distribution” or “distribution function.” However, the PDF is also commonly referred to as a “distribution” as well as “frequency function” throughout the literature. While “PDF,” “distribution,” and “distribution function” are commonly interchanged with no change in meaning (in context), the term “PDF” will be principally used herein.

*dP*of a random variable

_{X}*X*occurring in the differential interval [

*x, x + dx*] is expressed in terms of the probability density function $ p X ( x )$ of

*X*,

*X*is equal to unity.

*X*up to an arbitrary point

*x*is determined from the integral in Eq. (9) and is referred to as the cumulative distribution function,

*X*for all values above

*x*can be related to the integral of the probability density function as

*not*from the target of interest. Conversely, when the PDF is used to describe the anticipated echo from a scatterer (or “target”) of interest, the PD gives a measure of the probability of detecting the target of interest. The limit,

*x*, in Eq. (12) is the threshold above which the PFA and PD are calculated (Chap. 7 of Ol'shevskii, 1978).

#### 2. Calculating PDFs: Directly and from Monte Carlo simulations

Some PDFs are, conveniently, closed-form analytical solutions that can be calculated directly. For example, the Rayleigh and K PDFs given later, as well as the Gaussian PDF are in a simple analytical form and are straight forward to calculate. Other analytical solutions, such as those given in integral form, are also straight forward to calculate through numerical integration. However, PDFs for many realistic cases are generally not in closed form and require numerical simulations involving the scattering equations to estimate them. For example, once parameters of the sensor system and scatterers are specified, models of the echo at the *signal* level (not PDF level) such as in Eq. (6) are simulated to create the statistics of the echo.

A common method to estimate PDFs is through use of Monte Carlo simulations which generally involves making calculations of many statistical realizations of the signal to form an ensemble from which the PDFs are estimated. For example, in Eqs. (6) and (7), each *i*th scatterer will have a distribution of scattering amplitudes $ f b s ( i )$ and locations ( $ \theta i , \u2009 \varphi i$). Each realization involves randomly selecting one of those scattering amplitudes and locations for each scatterer, then calculating the signal using Eq. (6). The process is repeated many times, with each selection of random values being statistically independent (when independence is required) from the other selections.

Use of Monte Carlo simulations tends to be more general than analytical methods as there are fewer assumptions and, hence, limitations. For example, the scatterers can be arbitrarily correlated in space (such as being in a thin layer or small patch) and can have arbitrary phase-shift distributions. While many analytical models are restricted to cases involving independently distributed scatterers whose phases are uniformly distributed over [0 2 $\pi $], use of the Monte Carlo simulations do not have such restrictions.

Once the many simulations have been completed, the PDF is commonly estimated by putting the realizations of the signal into “bins” to form a histogram. For an analysis of echo magnitude statistics in the example described above, the result of each calculation is put into a magnitude bin (i.e., quantized value of echo magnitude) so that a histogram can be formed. These simulations require many realizations so that there can be correspondingly many narrow bins in order to produce a histogram that is an accurate representation of the actual PDF.

This binning approach is intuitive and is a method used in this paper, when appropriate, due to its simplicity. Conditions under which this method are used depend upon a combination of the type of structure in the PDF and corresponding number of realizations required to model that structure. In some cases, such as for a smoothly varying PDF, the computation time is reasonable. Calculations of PDFs for other applications where there is structure such as the presence of narrow peaks or nulls in the PDF curve may require more realizations and correspondingly significant computer time. When making too few calculations in this latter case, there can be artifacts in the result, such as smoothed or completely missed peaks or nulls. Thus, when there is the presence of narrow peaks or nulls in the PDF, a closed-form analytical method is used in this paper to determine the PDF, when possible. Beyond these approaches, the kernel density estimation (KDE) method was used to reduce the number of realizations needed to produce a reliable estimate of the echo PDF (Botev , 2010; Lee and Stanton, 2015; Scott, 1992). The calculations illustrated in this paper typically involve 10^{7} realizations.

Finally, for applications that extract information from the tail of the PDF, estimation methods such as importance sampling can be used to reduce the variance in the estimate and to increase the efficiency of the Monte Carlo process through selectively sampling the more desired (tail) samples (Agapiou , 2017).

#### 3. Non-uniform spacing of bins

Depending upon the types of features one is investigating in a PDF, the curves will either be plotted on a linear-linear or logarithmic-logarithmic scale. While the former scale may be more intuitive, the latter is especially useful when examining the tail of the PDF which typically has low values relative to the maximum. Choice of type of scale influences how the bins are determined. For linear-linear plots of PDFs, equally-spaced bins for the horizontal axis are normally used. However, when plotting PDFs on a logarithmic-logarithmic scale, the width of the bins should be equal on a log scale, which is non-uniform on a linear scale. Otherwise, if the bins were equal on a linear scale, but plotted on a log scale, the density of points on the plots would increase throughout the plot, and not fully characterize the shape of the PDF.

#### 4. Normalization

##### a. Vertical scale.

The probability of a variable occurring over *any* of the values of the random value *x* over the entire range is, by definition, unity. Therefore, the integral over *x* of any PDF over all values is unity. PDFs are commonly derived with a constant factor introduced that is determined through normalizing the area under the PDF curve to unity. From this property, it follows that the CDF will begin at a value of 0 for the smallest value of *x* and reach its maximum value of unity at the largest value of *x.* Similarly, the PFA will begin at unity and decrease to the value of 0 for the corresponding smallest and largest values of *x*, respectively.

##### b. Horizontal scale.

In some applications, it is also important to normalize the (horizontal) scale associated with the random variable. This can be the case when the calibration of the system is not known accurately, the propagation loss of the signal in the medium is not known accurately, or when only the shape of the PDF, CDF, and PFA are of interest regardless of the echo strength. Through normalization, only the relative values of the random variable will be considered. One convenient approach is to normalize the random variable by its root-mean-square (*rms*) value $ \u27e8 x 2 \u27e9 1 / 2$ and plot the PDF, CDF, and PFA versus the random variable divided by $ \u27e8 x 2 \u27e9 1 / 2$, where $ \u27e8 \cdots \u27e9$ is the average over a statistical ensemble of values. In this case, the area under the PDF curve (with an argument normalized by $ \u27e8 x 2 \u27e9 1 / 2$) is preserved under the transformation and is unity.

Regardless of the type of scale (linear-linear or logarithmic-logarithmic) or uniformity of spacing of bins, all normalizations are first calculated on a linear-linear scale. For example, an equation such as Eq. (9) is on a linear-linear scale and can be used to normalize the PDF to unity while accounting for non-uniform spacing in the integral.

### C. Fundamental statistical processes relevant to echo statistics

#### 1. Randomizing the deterministic scattering equations

Random fluctuations of echoes involve several fundamental statistical processes. For example, in Eq. (1), the beampattern is shown to be a function of the angular coordinates. In general, the scatterer will be randomly located in the beam, making the angular coordinates of the scatterer random variables. Since the beampattern is a function of the random angular coordinates, then the beampattern function is, in turn, a random variable for a randomly located scatterer. The scattering amplitude in Eq. (1) is also generally a randomly variable due to the random nature of the scatterer. Since the echo $ e \u0303$ in Eq. (4) as measured through the receiver of the sensor system is the product of the two random variables, the beampattern function and scattering amplitude, then $ e \u0303$ is also a random variable. Finally, once there are multiple scatterers in the beam of the sensor system, the resultant echo $ e \u0303$ [Eq. (6)], as measured through the receiver, will be the sum of the random individual echoes [Eq. (7)] and will, in turn, be a random variable.

These statistical processes—function of a random variable(s), multiplication of two random variables, and addition of random variables—are of wide applicability, are not specific to sensor systems or scattering, and appear in standard textbooks on statistics. Formulas summarizing these general processes are given below for later reference in the scattering problem. While only the simplest of cases involving one or two random variables are given, formulas involving more random variables are given in the references and/or later in context of the application.

#### 2. Function of a single random variable

If the function *Z* is a function of the random variable *X*, then *Z*(*X*) is also a random variable. The formulations relating the PDF of *Z* to the PDF of *X* are based on the fundamental principle that the probability of occurrence of an event in one space (*X* in this case) is the same as that in the transformed space (*Z* in this case). The resultant PDF $ p Z$(z) is then given by one of two equations depending upon whether *Z*(*X*) varies monotonically or non-monotonically with respect to *X*. Specifically, *Z*(*X*) is monotonic with *X* if it either solely increases or solely decreases over the range of *X* such that for any value of *Z*, there is only one (unique) value of *X*. Conversely, for the non-monotonic case, *Z*(*X*) both increases and decreases over the range of *X* so that there can be multiple values of *X* for a given value of *Z* [both of these cases and the below equations are described on pp. 23–27 of Goodman (1985)].

*Z*(

*X*) varies monotonically with respect to

*X*over the entire range of

*X*, then the following expressions can be written where the differential probabilities in the two spaces are equated to each other,

*X*and

*Z*,

*Z*in terms of the PDF of

*X*for this monotonic case,

*Z*(

*X*) varies

*non-monotonically*with respect to

*X*over the range of

*X*, the PDF is described by a similar equation, but summed over

*M*contiguous segments where

*Z*(

*X*) varies monotonically within each segment,

*x*(

*z*) and

*x*(

_{m}*z*) in Eqs. (16) and (17) are the inverse functions

*z*

^{−}^{1}(

*x*) and

*z*

^{−}^{1}(

*x*), respectively. In practice, these inverse functions can be determined numerically from the forward analytical function, plots, or tables of

_{m}*z*(

*x*) and

*z*(

*x*).

_{m}#### 3. Function of two random variables

The above analysis involving a function of one random variable is extended to the case of a function of two random variables. In this case, if the function *Z* is a function of the random variables *X* and *Y*, then *Z*(*X*, *Y*) is also a random variable. Relating the PDF of *Z* to the PDF(s) of *X* and *Y* involves the same process as in the previous case of one random variable in which the probability of occurrence of an event in one space is set equal to that of the other space. This process generally involves first determining the Jacobian of the transformation relating the two spaces, although that will not be shown explicitly below (Papoulis, 1991).

*Z*occurring for any value below

*z*is given in terms of

*x*and

*y*as

*P*( $ Z \u2264 z$) is also the

_{Z}*CDF*. Here, $ p X , Y ( x , y )$ is the joint probability density function of the random variables

_{Z}*X*and

*Y*, and

*D*is the region or regions in the

_{Z}*xy*plane containing values of

*x*and

*y*where

*Z*(

*X,Y*) $ \u2264 z$ (

*D*is illustrated in Fig. 6-7 of Papoulis, 1991). This equation for

_{Z}*P*( $ Z \u2264 z$) is a two-dimensional form of Eq. (10). The PDF of

_{Z}*z*can be expressed by taking the differential of $ P Z$ ( $ Z \u2264 z$) above,

*dD*is now the differential region or regions(s) in the

_{z}*xy*plane whereby the values of

*x*and

*y*are bounded by the differential area determined by the range $ z \u2264 Z \u2264 z + d z$ [Eq. (6-36) and Fig. 6-7 of Papoulis, 1991].

This equation is complex to solve and depends upon the characteristics and form of $ p X , Y ( x , y ) .$ For simple forms such as *Z = XY* and *Z = X + Y*, where *X* and *Y* are independent random variables, then the solution for *p _{Z}(z*) in each case is in closed form. Since those two cases are used throughout this tutorial, they are treated separately in Secs. IV C 4 and IV C 5. The more complex case in which

*Z*is a general function of

*X*and

*Y*is used only once later in the tutorial and the solution will be given in the context of that application (Sec. VII A 7).

#### 4. Product of two random variables

*X*and

*Y*are random variables, then the product

*Z = XY*is also a random variable as described in Sec. IV C 3. If

*X*and

*Y*are independent of each other, then $ p X , Y ( x , y )$

*= p*(

_{X}*x)p*(

_{Y}*y*) in the integrand in Eq. (19), where

*p*(

_{X}*x*) and

*p*(

_{Y}*y*) are the PDFs of the random variables

*X*and

*Y*, respectively. Inserting this product of the two PDFs into the integrand in Eq. (19), the PDF $ p Z ( z )$ of the product

*Z*=

*XY*can be shown to be

*product*of two random variables (

*Z = XY*), the change in variables

*y = z/x*is used, and the area

*dxdy*for

*dD*is mapped to the area (

_{z}*1/|x|*)

*dxdz*. Through this mapping, the double integral for

*dD*is replaced with a single integral over

_{z}*x*. Replacing

*dxdy*in Eq. (19) with (

*1/|x|*)

*dxdz*, the

*dz*drops out of both sides of the equation and the integral is only over

*dx*as shown. The absolute value sign is used for the variable

*x*so that the differential area will be positive for all values of

*x*. Note also that the term

*|x|*in the factor (

*1/|x|*) for the area is equal to the absolute value of the Jacobian of the transformation in the derivation of Eq. (6-74) in Papoulis (1991).

Once Eq. (20) is used in physical applications, the range over which one or more of the physical parameters may be constrained and its corresponding PDF will be zero outside of that range. The integration limit(s) may reflect that constraint by only spanning the range over which the PDF is non-zero as shown later.

#### 5. Sum of random variables

There is a variety of methods to evaluate the PDF of the sum of independent random variables, ranging from purely analytical to purely numerical. Sometimes, a “purely” analytical method still requires numerical evaluation, such as when an integral or series summation are involved and numerical integration or summation are performed, respectively. Two common methods are discussed below: the method of characteristic functions and Monte Carlo simulations.

##### a. Method of characteristic functions.

A commonly used analytical method involves use of characteristic functions (CFs) where the CF of a random variable is the Fourier transform of its PDF (Goodman, 1985). Addition of an arbitrary number of independent random variables involves first taking the product of their corresponding CFs. This product is the CF of the sum of the random variables. The PDF of that sum is then the inverse Fourier transform of the CF product.

This CF approach can be derived from Eq. (18) for the case of *Z = X + Y* where *Z*, *X*, and *Y* are all random variables. Since the random variables, *X* and *Y*, are independent of each other, then $ p X , Y ( x , y )$* = p _{X}(x)p_{Y}(y*). Using this relationship in the integrand of Eq. (18), the PDF of

*Z*, $ p Z ( z )$, can be shown to be the convolution of the two functions,

*p*(

_{X}*x*) and

*p*(

_{Y}*y*). This convolution can then be shown to be equivalent to the product of the Fourier transforms of

*p*(

_{X}*x*) and

*p*(

_{Y}*y).*Since these Fourier transforms are, as defined above, the CFs of the two functions, then the method of characteristic functions follows as described above (Goodman, 1985; Papoulis, 1991). The method is extendable to the sum of an arbitrary number (

*N*) of independent random variables by first expressing the convolution integral by formulating the sum of two random variables where one of the random variables is the sum of

*N*−1 random variables and the other random variable is the remaining variable. The PDF of the sum of the

*N−*1 random variables is determined through a similar process involving the sum of

*N−*2 random variables, and so on. After completing this iterative process, the PDF of the summed

*N*random variables is related to the product of the Fourier transforms of the PDFs of the

*N*random variables.

Acoustic and electromagnetic signals are complex and normally constructed of a real and imaginary term, making them two dimensional. Since the method of characteristic functions is extendable to multi-dimensional variables, this method can be applied to determine the PDF of the sum of complex signals. For the case in which the phase of the summed signal is uniformly distributed [0 2 $\pi $] and each component of the signal has a zero mean (i.e., a “circularly symmetric signal” in the complex plane), then the CF and PDF of the signal magnitude are a Hankel transform pair. Application of the CF to calculate the PDF of the magnitude and magnitude squared of complex signals is summarized in Jakeman and Ridley (2006, Chap. 4), including an extension to signals where the phase is not uniformly distributed. Methods to numerically evaluate the PDF for circularly symmetric signals (via the CF and Hankel transform) are given in Drumheller (1999).

Barakat used a broadly similar approach to the Hankel transform method by extending the 1D CF method to circularly symmetric complex signals through constructing an orthogonal component of the sum, equal in magnitude to the original single component of the sum, resulting in an exact, analytical expression for the PDF of the magnitude of the sum of complex random variables (Barakat, 1974). However, our experience in applying the Barakat approach has resulted in convergence issues due to truncation of the infinite series summation that must be evaluated (Chu and Stanton, 2010; Lee and Stanton, 2014).

##### b. Monte Carlo simulations.

The method of Monte Carlo simulations is discussed in more general terms in Sec. IV B 2 and will only be briefly summarized here in the context of adding random variables. This is a commonly used numerical approach that involves simulating a statistical ensemble of a large number of realizations of the process of interest so that the PDF can be formed. Performing Monte Carlo simulations to evaluate the sum of random variables in predictions of scattering will provide a stable and accurate solution (see, for example, Stanton , 2015). When using this method to simulate signals associated with a sensor system, random variables are normally added through (1) a phasor addition in the frequency domain of single frequency or narrowband signals of long enough extent that the echoes are completely overlapping [such as with Eq. (6)] or (2) addition in the time domain of short signals when the echoes are only partially overlapping and/or when broadband signals of any duration are used (Sec. VIII A). In the case of phasor addition, the signals are first represented in complex form and then the real and imaginary components are added separately before being recombined to calculate the signal magnitude.

#### 6. Sum of infinite number of random variables (central limit theorem; Rayleigh PDF)

The Rayleigh PDF is widely used in describing echo statistics. It is commonly used as the “starting point” in describing the statistics of the echo magnitude, especially when there are many scatterers or many highlights from an individual scatterer contributing to the echo. When the statistics do not follow the Rayleigh PDF, deviations of the statistics from the Rayleigh PDF are frequently described. The deviations in the higher values of the echo magnitude, i.e., the “tail” of the distributions, are of particular interest. The term “non-Rayleigh statistics” is commonly associated with those distributions that deviate from the Rayleigh PDF. The Rayleigh PDF and associated CDF and PFA are illustrated (Fig. 5).

## V. IN-DEPTH TREATMENT OF ECHO STATISTICS: OVERVIEW

As discussed in Sec. I, various important aspects of echo statistics will now be examined in detail. The treatment will draw from the concepts and equations given in Secs. II–IV. Generally, deterministic equations for the echo magnitude [Eqs. (4) and (6)], which are based on solutions to the wave equation, are randomized with respect to various physical quantities such as random location in beampattern and random orientation of scatterer. They are randomized using fundamental statistics equations given in Sec. IV C. This approach differs from other approaches such as first randomizing parameters of the governing differential equation before solving the equation [see, for example, the summary in Sec. 12.6 of Jakeman and Ridley (2006), and references therein].

In Sec. VI, the treatment begins with the simplest of cases—single-frequency signals of infinite extent in which echoes from all scatterers completely overlap and direct path propagation in a homogeneous medium without the complicating effects due to reflections from boundaries. These idealized signals approximate gated sine waves where the echoes overlap significantly. In addition, given the importance of the random modulation effect of the beampattern on the echo when a scatterer is randomly located in the beam, the echo statistics are first described without those effects. This is equivalent to the sensor system having an omnidirectional beam. In Sec. VII, beampattern effects are then incorporated into the analysis. In Sec. VIII, the work is further extended into more realistic and complex cases involving pulsed signals (narrowband and broadband) in which echoes may only partially overlap, the presence of a boundary near a scatterer, and propagation and scattering in a waveguide with a heterogeneous medium.

### A. How to use this material for realistic signals/environments, and advanced signal/beam processing

The number of combinations of types of systems, signals, signal processing and beamforming algorithms, and environments is limitless and cannot be adequately described within this tutorial. The material is therefore aimed toward the more fundamental aspects of echo statistics that are not specific to any particular system or environment, but that many applications either have in common or could use as a basis. For example, most material involves the following: (1) *Signal type*: Long, narrowband signals are used in which the echoes from all scatterers completely overlap (in all cases in Secs. VI and VII, and some cases in Sec. VIII). (2) *Signal processing*: The magnitude of the signal as measured at a single instant in time is measured—that is, “first-order statistics” is modeled (in all cases in Secs. VI–VIII). The instant in time may be fixed or randomly selected, but it is not adaptively chosen according to a particular echo magnitude. (3) *Processing of beam data and/or beamforming*: Echoes from a single beam are modeled (fixed or scanning; or one selected from a multi-beam system) in which the scatterers are randomly distributed in space (in all cases in Secs. VII and VIII). The echo is sampled from a single beam for a random spatial distribution of scatterers and the beam is not steered adaptively to select or focus on a particular scatterer. (4) *Environment*: direct path geometries in which the medium is homogeneous and there are no reflecting boundaries (in Secs. VI and VII; and one case in Sec. VIII).

The following cases for systems, signals, and environments of greater complexity are examined in Sec. VIII: (1) *Signal type and signal processing*: pulsed signals are modeled in one example in which the echoes from the scatterers only partially overlap (Sec. VIII A). In Sec. VIII A, the pulse is further shortened through use of matched filter processing. (2) *Environment*: in two examples, geometries in which there are one or more boundaries present and there are heterogeneities in the medium (Secs. VIII B and VIII C, which cover boundary interaction and waveguide effects).

The several cases modeled in Sec. VIII, while far from spanning the many possible complex scenarios, provide examples for how the fundamental formulations involving the more simple cases can be applied to the more complex cases. For example, the case involving a pulsed signal shows how a time series can be constructed due to the interference between the partially overlapping echoes from the scatterers (Sec. VIII A). The examples involving the presence of one or more boundaries and/or heterogeneities in the medium show the different types of effects associated with the boundaries and medium heterogeneities. For a single boundary near a scatterer, it can provide an added source of interference due to the interaction of the incident signal and the boundary and scatterer (Sec. VIII B). For two parallel boundaries and/or a medium with a local minimum in wave speed, a waveguide is formed and the signal can propagate along multiple paths that are guided by the boundaries or local minimum (Sec. VIII C). These multiple paths represent additional sources of fluctuations in the echoes.

### B. Peak sampling, pulsed signals with boundaries, and beyond

There are many more important cases not covered explicitly in this tutorial, but that can be described using these formulations. For example, there are systems that use peak sampling signal processing, such as recording the maximum echo magnitude in a time window or adaptively steering a beam toward the scatterer with the largest echo. When multiple scatterers are present, the process of peak sampling in both cases (time and angle/space) will bias the statistics toward higher values than the magnitudes modeled with first-order statistics and a fixed beam. This process involves “extremal” statistics (Stanton, 1985), which is outside the scope of this tutorial. However, for the time-based peak sampling, the method in Sec. VIII A that produces a time series could then serve as a basis of the time series with which a peak sampling algorithm could be applied. As shown in Stanton (1985), the bias in this case increases with the ratio of window duration to ping duration of the signal. For the case of a scanning beam or multi-beam system adaptively focusing on the peak echo in a field of multiple scatterers, this extremal statistics formulation can be adapted from a time series to a space series of echoes scanned across angles. The bias here will increase as the angular window is increased.

Another important case involves use of pulsed signals in the presence of boundaries in which the echoes from the scatterers are only partially overlapping. Here, the method given in Sec. VIII A to produce a time series can be incorporated into the formulations in Secs. VIII B or VIII C that involve the presence of boundaries. Furthermore, advanced signal processing such as peak sampling can also be incorporated into this case as described above.

## VI. IN-DEPTH TREATMENT OF ECHO STATISTICS: NO BEAMPATTERN EFFECTS

Given the complex effects of the beampattern on the echo statistics (Fig. 4), fluctuations of the echoes without the influence of the beampattern are first examined separately. This is equivalent to a sensor system with an omnidirectional beam so that the echo value is the same regardless of angular location (at constant range) in the beam. In this simplified case, fluctuations can be due to a combination of interference between overlapping echoes and scattering effects. Fluctuations due to those effects are treated separately below.

### A. Addition of random signals (generic signals; not specific to scattering)

This initial treatment solely involves analyses of generic random signals, devoid of explicit representation of the sensor system and field of scatterers. This analysis forms the foundation for the more complex analyses given later that first involve scattering from objects and, eventually, system effects. The case of an arbitrary number of arbitrary signals is first presented, which is then followed by some commonly used special cases, including the Rice and K PDFs.

#### 1. Arbitrary cases

*N*arbitrary complex signals completely overlap, the resultant signal

*A*is the coherent sum

Since $ a i$ and $ \Delta i$ are random variables, then so is *A*. The fluctuations of *A* from realization to realization depend strongly on the statistical properties of $ a i$ and $ \Delta i$. The phase shifts $ \Delta i$ play a major role in the fluctuations. For example, for the simple case in which $ a i$ is constant for all *i* (i.e., $ a i = a$) and $ \Delta i$ is randomly and uniformly distributed over the range 0–2 $\pi $, *A* will fluctuate greatly from realization to realization due to variability in constructive and destructive interference effects associated with phase variability alone (Figs. 2 and 6). In one realization, there may be complete constructive interference and *A* is at a maximum. In another realization, there may be complete destructive interference and *A* is at a minimum. And, generally, *A* will take on intermediate values due to partial interference.

The characteristics of the fluctuations in this case also depend greatly on *N* as illustrated. When there is only one signal (*N* = 1), the signal is single valued and the PDF of the signal magnitude $ | A |$ is the delta function (Fig. 2). In the other extreme in which there are an infinite number of signals (*N* = $\u221e$), the PDF of $ | A |$ is the Rayleigh PDF as given in Eq. (21) (Fig. 2). The PDF of $ | A |$ takes on other shapes for intermediate values of *N* (Fig. 6). In the more general case in which $ a i$ is a random variable (not equal), the curves will fluctuate in a similar fashion and, in the limit of *N* = $\u221e$, the PDF of $ | A |$ will become Rayleigh via the CLT.

Equation (24) is broadly applicable to the scattering problem as it could represent the summation of scattering highlights from within a single scatterer or the summation of echoes from multiple scatterers. Characteristics of the scatterer, sensor system, and scattering geometry can be incorporated into *N*, $ a i$, and $ \Delta i$.

#### 2. Sine wave plus noise (Rice PDF)

*N*must be sufficiently large so that the magnitude of the summation converges to a Rayleigh random variable. The PDF of $ | A |$ is the Rice PDF (Rice, 1954),

*rms*value of the noise term, $ \u27e8 | A | 2 \u27e9$ is the mean square of the sine wave plus noise [i.e., mean square of Eq. (25)], and $ I 0$ is the zeroth-order modified Bessel function of the first kind. Note that in the original derivation by Rice, the noise term in Eq. (25) can also involve the summation of signals of different frequencies as well (i.e., more general than the single frequency case shown here).

The shape of the Rice PDF depends strongly on $\gamma $ (Fig. 7). For example, in the limit as $\gamma $ approaches infinity, the PDF is close to a Gaussian PDF. This corresponds to the limit of high SNR where the signal is dominated by the constant sine wave. In the other extreme, as $\gamma $ approaches zero, the PDF approaches the Rayleigh PDF. This latter case corresponds to the limit of low SNR where the signal is dominated by the noise. The shape of the Rice PDF changes smoothly for all intermediate values of $\gamma $.

In addition to this formula being widely applicable in modeling noisy signals, it can also be applied to scattering problems where the “signal” [sine wave in Eq. (25)] is the mean scattered field and the “noise” [summation term in Eq. (25)] is the component of the scattered signal that fluctuates about the mean (Stanton and Clay, 1986; Stanton and Chu, 1992). For example, it could be used to model the fluctuations of the echo from a sphere near a rough interface where the individual echoes from the sphere and rough interface are delta-function- and Rayleigh-distributed, respectively. Or, rather than a rough interface, the sphere could be surrounded by a cloud of smaller scatterers whose individual echoes are of random phase. If the constant sine wave in Eq. (25) is considered to represent a single scatterer of interest (such as the sphere) and the noise term represents the background reverberation (such as from the rough interface or cloud of smaller scatterers), then the term $\gamma $ is the signal-to-reverberation ratio.

#### 3. Special distributions of N or a_{i} (K PDF)

*N*and/or $ a i$ in Eq. (24) are random variables. This can be divided into two categories—one in which $ \Delta i$ is randomly and uniformly distributed over the range 0–2 $\pi $, and the other in which $ \Delta i$ is non-uniformly distributed. The former category will first be discussed. If

*N*follows a negative binomial PDF and its average value tends to infinity, the statistics of $ | A |$ in Eq. (24) for arbitrary $ a i$ is described by the K PDF (Jakeman and Pusey, 1976; Abraham and Lyons, 2002b),

*K*is the modified Bessel function of the second kind (and served in the naming of the PDF), Γ is the gamma function, $ \alpha K$ is the shape parameter, and $ \lambda K$ is a scale parameter equal to the mean square of the signal divided by $ \alpha K$.

The K PDF has a single mode and varies smoothly with $ | A |$ (Fig. 8). The distribution tends to the Rayleigh PDF in the limit of high $ \alpha K $. Note that the function *K* also has other, less common names, including the Basset function and the modified Bessel function of the *third* kind.

The K PDF can also be derived in several other ways. For example, it has been shown that for finite *N* and if $ a i$ follows an exponential PDF, then $ | A |$ is K-distributed (Abraham and Lyons, 2002b). Beyond methods involving summing sinusoidal signals, the K PDF has been shown to be due to the product of two independent random variables: a Rayleigh-distributed term and one that is chi distributed (Ward, 1981). This product was later written in extended form as the product of a Rayleigh-distributed term and the square root of a term that is gamma distributed (Abraham and Lyons, 2002b). Here, the square root of the gamma-distributed term is related to the chi-distributed term through analytical continuation of the integer number of summed terms in the chi distribution to a non-integer number in the gamma distribution. In another derivation, the K PDF has also been shown to result from a Rayleigh PDF whose mean-square value is gamma distributed (Jakeman and Tough, 1987). Both of these latter derivations are referred to as a “compound representation.”

Equation (29) has been widely used to describe echo statistics in both acoustic and electromagnetic applications. While there has generally not been a direct connection to the physics of the scattering, the tail of the distribution has generally followed those from experimental data after an empirical fit. As discussed above, through the interpretation of Abraham and Lyons (2002b), the number of scatterers has been related to the shape parameter of the K PDF for the specific case in which the amplitudes $ a i$ of the individual echoes are exponentially distributed. Since $ a i$ in this case is observed through the receiver of the sensor system, the exponential PDF includes the effects of both fluctuations from the stochastic nature of the scatterer and the variability due to the scatterer being randomly located in the beam. Details of those effects are given in Sec. VII. Also, in his expression of a K-distributed magnitude being due to the product of Rayleigh- and chi-distributed random variables, Ward (1981) attributed the Rayleigh term as being due to quickly varying interference between scatterers (such as from phase shift differences within a patch of scatterers) and the chi term being due to slowly varying changes in the echo from larger-scale variations in the “bunching” or patchiness of scatterers.

The above K PDF involves signals who phases are uniformly and randomly distributed over [0 2 $\pi $]. However, when the distribution of phases is non-uniform, then $ | A |$ can be described by the *generalized* K-distribution (not shown) (Jakeman and Tough, 1987). This distribution is described by three parameters— $ \lambda K$ and $ \alpha K$ as given above, plus a third that describes the non-uniform phase distribution. In addition, the generalized K PDF in Jakeman and Tough (1987) can also describe an n-dimensional random walk (generalized from the two-dimensional walk for the standard K PDF). For the case of a non-uniform phase distribution, this is a random walk with directional bias. Similar to the compound representation of the K PDF above being a Rayleigh PDF with its mean-square value being gamma distributed, with this directional bias, the generalized K PDF can be (compound) represented by a Rice PDF with the mean square noise and constant amplitude signal components each being gamma distributed in a correlated way (Jakeman and Tough, 1987). The generalized K PDF can be applicable to the case in which one or several scatterers dominate the scattering from within a field of many scatterers, thus skewing the distribution of phases into one that is non-uniform (Ferrara , 2011).

Another generalization of the K PDF is the *homodyned* K-distribution (not shown). Like the generalized K PDF, it is also a three-parameter distribution. And, like the generalized K PDF, it can also be (compound) represented by a Rice PDF, but with only the mean square noise component (but not the constant amplitude component) being gamma distributed (Jakeman and Tough, 1987).

The properties and potential relations to scattering of the K PDF, generalizations of the K PDF, and other generic PDFs are summarized in Destrempes and Cloutier (2010).

#### 4. Adding independent realizations of the complex signal A

As discussed in Secs. IV C 6 and VI A 1, in the limiting case of the sum of an infinite number of random complex variables, the PDF of the magnitude of the sum is Rayleigh distributed (as per the CLT). Specifically, for Eq. (24) where the phase shifts $ \Delta i$ are randomly and uniformly distributed over the range 0–2 $\pi $, in the limit of *N* = $\u221e$, $ | A |$ is Rayleigh distributed. Now, consider the case in which there is an ensemble of statistically independent realizations of the complex signal *A* (in the limit of large *N*) where the magnitude of each realization of *A* is Rayleigh distributed and with the same mean square value. Then by extension of the CLT, the magnitude of the sum of those realizations of *A* is likewise Rayleigh distributed, even for a finite number of realizations. However, if each realization of *A* is modulated by a multiplicative term with a magnitude that is a random variable, the resultant magnitude of the sum of a finite number of these (complex) signals can possibly be strongly non-Rayleigh. As described in Sec. VI A 3, such a modulation can be caused by patchiness of the scatterers, resulting in echo PDFs that can be derived through a compound representation. Also, when the scatterers are uniformly distributed (i.e., no patchiness effects), there can also be modulation caused by the beampattern where the echo from the scatterer is randomly modulated by its random location in the beam, causing strongly non-Rayleigh echoes as described in Secs. VII B and VII C.

### B. Complex scatterers with stochastic properties

Scatterers can range in complexity from the simplest of form, a point scatterer, in which the scattering amplitude is constant for all orientations, to an arbitrarily shaped object whose echo varies from orientation to orientation. In this section, the statistics of scattering by an individual are examined in a progression of complexity. The point scatterer, Rayleigh scatterer (defined below), and smooth and rough prolate spheroids with both fixed and random orientation are modeled. A summary of scatterers with other complexities are given at the end.

#### 1. Point scatterer

The simplest case is the point scatterer or, more precisely, very small scatterer. The dimensions of this scatterer are sufficiently small compared with an acoustic or electromagnetic wavelength that the echo is constant for all orientations. It is constant because there is only one scattering highlight from this object and nothing else with which to interfere. The PDF of the echo magnitude is the delta function [Figs. 2(a) and 3(a)].

#### 2. Rayleigh scatterer

At the other extreme from a point scatterer is a scatterer whose echo is Rayleigh distributed. This so-called “Rayleigh scatterer” can be in many forms. For example, it could be a small patch of many point scatterers, each with an echo whose phase is randomly and uniformly distributed over the range [0 2 $\pi $]. Or, it could be a single spherical scatterer whose surface is randomly rough and can be described as a collection of many scattering highlights bounded by the surface. The echo from each highlight has a phase that is randomly and uniformly distributed over the range [0 2 $\pi $]. In each case, from the central limit theorem, the echo from the patch or rough sphere is Rayleigh distributed [Fig. 2(c); Eq. (21)]. The fluctuations occur from ping to ping as the patch or rough sphere are rotated, or from realization to realization of a randomized spatial distribution of scatterers in the patch or randomized roughness of the rough sphere.

#### 3. Randomized prolate spheroid

A sequence of formulations is presented, beginning with the deterministic description of the scattering by a smooth prolate spheroid at fixed orientation, then randomizing its orientation, and then further randomizing the spheroid by roughening its boundary. In contrast to the Rayleigh scatterer, the echoes from the randomized prolate spheroid are generally non-Rayleigh because of the elongated shape of the prolate spheroid. The degree to which the echoes are non-Rayleigh can be connected to various parameters of the scatterer through the physics-based formulas given below. The formulas given below are adapted from Bhatia (2015).

##### a. Smooth boundary, fixed orientation.

*β*is the angle between the direction of the incident acoustic or electromagnetic wave and the plane that is normal to the lengthwise axis of the prolate spheroid [Fig. 9(a)]. Note that this equation is equivalent to Eq. (7) of Bhatia (2015), but is in a form in which no terms have singularities.

*β*= 0 and

*β*= $ \pi / 2$ correspond to broadside and end-on incidence, respectively, relative to the incident wave. The terms $c$ and $ b 1$ are the lengths of the semimajor and semiminor axes of the prolate spheroid, respectively (the length and width of the spheroid are $ 2 c$ and $ 2 b 1$, respectively). The spheroid is axisymmetric about the length-wise axis, leading to only one unique value of semiminor axis. The term $ b 1$ is not to be confused with the beampattern

*b*. The aspect ratio of the spheroid is defined to be the ratio $ c / b 1$ (or, equivalently, length/width). The scattering is a strong function of the aspect ratio and orientation [Fig. 9(b)]. At broadside incidence, the above formula reduces to

##### b. Smooth boundary, random orientation.

We now randomize the scattering by first randomizing the orientation. This is done by making *β* a random variable with an associated PDF $ p \beta ( \beta )$. For simplicity, the prolate spheroid will only rotate in a single plane about its minor axis in a plane containing the direction vector of the incident plane wave [Fig. 9(a)].

*β*, then $ f s s$ is a random variable as well. Inserting $ | f S S |$ and $ p \beta ( \beta )$ into Eq. (16), the PDF of the magnitude of the scattering amplitude of a randomly oriented smooth prolate spheroid is

*β*is

*β*only varies over the range 0– $\pi $/2 because of the symmetry of the scattering over the other angles.

The echo PDF for the randomly oriented smooth prolate spheroid is calculated for a range of aspect ratios (not including the case of a sphere where the aspect ratio is unity) by inserting Eq. (34) into Eq. (33) [Fig. 10(a)]. For each (non-unity) aspect ratio, the echo PDF is characterized by a smoothly varying function for most magnitudes, but with strong narrow peaks at the end points. These peaks are attributed to the fact that the backscattering near broadside and end-on incidence varies slowly with orientation angle, which increases the probability of occurrence at those corresponding echo values. The range of echo values as well as probability of occurrence are both shown to be a strong function of aspect ratio of the spheroid. For the case of a sphere, the echo is constant for all *β* and the echo PDF is a delta function [Fig. 10(a)].

##### c. Rough boundary, random orientation.

*rms*and, for each orientation, the term is statistically independent from its value at all other orientations. For a randomly oriented prolate spheroid, $ | f s s |$ is also a random variable. Thus, the magnitude of the scattering amplitude for the randomly rough, randomly oriented prolate spheroid is the product of two random variables, $ | f s s |$ and the Rayleigh-distributed modulation term. The statistics of the echoes from randomly rough, randomly oriented prolate spheroids can be described using Eq. (20),

_{min}and $ | f s s |$

_{max}are the minimum and maximum values of the magnitude of the scattering amplitude of the smooth prolate spheroid, respectively, which correspond to end-on and broadside incidence. Those two terms replaced the limits − $\u221e$ and + $\u221e$, respectively, in the integral in Eq. (20) to reflect the fact that the magnitude of the scattering amplitude of the smooth prolate spheroid is within the range [ $ | f s s |$

_{min}, $ | f s s |$

_{max}] and that the corresponding PDF $ \u2009 p s s$ is zero for values of its argument outside that range. Without the above constraint, the integral still would have been constrained by the Rayleigh PDF whose argument is limited to only positive values, which would have led to integral limits of zero to infinity.

The resultant echo PDFs of the scattering by the randomly rough, randomly oriented prolate spheroid are significantly different from the smooth counterpart [Fig. 10(b)]. The curves for each aspect ratio are now smoothly varying, have a mode, and do not have any singularities.

#### 4. Randomly rough objects with full range of roughness and circumferential/internal waves

There are many more important complexities of the scattering and associated echo statistics than described above and they are summarized briefly below. The above example involved the simplified case in which the boundary of the prolate spheroid was impenetrable, only the echo from the front interface was accounted for, and its roughness was sufficiently large so that the echo was Rayleigh distributed for any fixed orientation. However, it is not uncommon for an object to have circumferential and internally refracted/reflected waves and/or a boundary roughness that spans small-to-intermediate levels relative to a wavelength of the incident signal. Each of these latter cases can lead to statistics of the echo that are significantly different than that of the very rough impenetrable prolate spheroid discussed above and must be accounted for.

Given the complexity of predicting the scattering by objects with those various types of waves, the echo statistics will be summarized below according to a progression of complexity. In Sec. VI B 4 a, effects on the scattering due to circumferential and internal waves for smooth objects are briefly reviewed. In Sec. VI B 4 b, the object is roughened, and the statistics of the echoes from only the front interface are considered. This is similar to the case of the rough impenetrable prolate spheroid above but with the added complexity of accounting for small-to-intermediate roughness where the echoes are non-Rayleigh. In Sec. VI B 4 c, circumferential/internal waves associated with penetrable materials are incorporated with boundary roughness included.

##### a. Circumferential and internal waves (smooth boundary).

Objects with common material properties not only produce a significant echo from the front interface, but they also support various classes of circumferential and internally refracted/reflected waves due to incident signals (*acoustic*: Hackman, 1993; Marston, 1992; Numrich and Uberall, 1992; Uberall, 1973; *electromagnetic*: Chen, 1998; Crispin and Siegel, 1968; Moser , 1981; Murphy , 1980; Newton, 1982; *acoustic and electromagnetic*: Bowman , 1987; Nussenzveig, 1969). For the case of impenetrable objects, such as infinitely dense materials for acoustics applications and perfectly conducting materials for electromagnetic applications, there is a weak signal that diffracts or “creeps” around the boundary and reradiates in all directions including back to the receiver. For the case of penetrable objects, such as elastic materials for acoustics applications and dielectric materials for electromagnetic applications, those diffracted waves also exist. In addition, with elastic objects, there can be strong surface elastic waves that are generated from the incident acoustic signal and will travel around the boundary (in addition to the diffraction boundary wave). Furthermore, other processes are also involved for the penetrable objects such as waves that refract into the interior of the object and reflect internally. Because of these various effects, the echo from a smooth symmetrical object such as a sphere or cylinder will generally be composed of the sum of echoes associated with the front-facing interface (sometimes called the “specular” echo), internal transmission and reflections, and circumferential waves.

The total echo from these objects will vary according to the constructive and destructive interference between the individual echoes associated with each of the different scattering phenomena described above. For a smooth object such as a sphere or cylinder, the phase of the echo associated with each phenomenon will vary strongly with signal frequency (or more precisely, *ka*, where *a* is the radius of the object). Because of these frequency dependences, which also vary with each phenomenon, the pattern of echo vs frequency will contain a series of peaks and strong nulls associated with the constructive and destructive interference, respectively, between the different phenomena.

##### b. Front interface only: Small-to-intermediate roughness, Rice PDF.

The following discussion concerns the echo from the front interface only, in isolation from the circumferential and internal waves.

When the roughness of the surface is smaller than the wavelength of the incident signal or, more precisely, when *k* $ \sigma B$ < 1, where $ \sigma B$ is the standard deviation of the boundary (or, equivalently, the *rms* deviation from the mean boundary), the magnitude of the echo from the front-facing surface of a rough object at normal incidence to that surface is generally not Rayleigh distributed. This is because the phases of the echoes from the individual scattering features of the surface are relatively narrowly distributed, in contrast to being uniformly distributed [0 2 $\pi $] such as in the above case of the (very) rough spheroid. This effect has been studied mostly in the context of scattering by rough planar interfaces and, to a much lesser extent, for rough bounded objects. In either case, in the limit of small roughness (i.e., *k* $ \sigma B$ ≪ 1), the echo PDF will tend to the delta function. In the opposite limit (i.e., *k* $ \sigma B$ ≫ 1), the echo PDF will tend to the Rayleigh PDF (or at least be Rayleigh-like) as in the above example.

It has been shown that the echo statistics from the full range of roughness at normal incidence of a planar boundary can be described by the Rice PDF, as given above in Eq. (26) (Stanton and Clay, 1986). In this formulation, the scattered signal is decomposed into the sum of two components, the coherent mean (constant amplitude) and fluctuation component. The PDFs of the magnitudes of the two components are the delta function and Rayleigh PDF, respectively. The mean component is related to the reflection coefficient of the interface modified by the term $ e \u2212 2 k 2 \sigma B 2$—a term originally derived for rough planar interfaces by Eckart (1953). The fluctuation term is related to the scattering phenomena that cause deviations in the echo from this mean. These two scattering terms, the constant amplitude and fluctuation components, correspond to the sinusoidal signal and noise, respectively, in the original Rice formulation.

The Rice PDF shape parameter has also been explicitly connected to parameters of the roughness and sensor system for rough planar interfaces (normal incidence), as summarized in Stanton and Clay (1986). Parameters of the roughness are the *rms* deviations ( $ \sigma B$) from the normal to the mean surface and the two-dimensional autocorrelation function (along the surface) of the surface. Parameters of the sensor system are the frequency and beampattern. Although the studies summarized in Stanton and Clay (1986) are specific to underwater acoustic signals, the analysis involving the Rice PDF is formulated for scalar fields and could be applied to other types of sensor systems, such as radar, when scalar field representations are appropriate.

These formulations connecting the Rice PDF to fluctuations due to randomized scattering by a rough planar surface were extended to the case of acoustic scattering by rough elastic cylinders immersed in a fluid in Stanton and Chu (1992) and Gurley and Stanton (1993). Because of the complexity of the scattering by the elastic objects, the analysis in Stanton and Chu (1992) was divided into two formulations—a general one which described all echoes (i.e., including those associated with the front interface, the circumferential waves, and internally refracted waves) and a simple analysis involving only the front interface. The fluctuations of the echo from the rough front (curved) interface were related to the Rice PDF in the same manner as with the case of the rough planar interface involving a scalar field, but taking into account curvature of the surface. This simpler approach is clearly an approximation to the other, more general, case where the other waves played a major role in the fluctuations, especially near nulls due to interference effects (as discussed below). However, the Rice PDF, when used to model fluctuations of the echo from the front interface only, was demonstrated in these studies to predict the general trend of the fluctuations as parameters such as *ka* and *k* $ \sigma B$ were varied, where *a* is the mean radius of the rough cylinder.

An important element in the modeling of both the above simple and general cases was the fact that the Rice PDF assumes a noise term whose quadrature (i.e., real and imaginary) phase components have the same mean square values. However, the quadrature components of the random component of the scattered (scalar) signal from rough interfaces are generally not equal. Furthermore, the fluctuations of the scattered signal are sensitive only to the component of the “noise” (i.e., fluctuation component of scattered signal) that is *in phase* with the mean scattered field. These effects are quantified and accounted for in Stanton and Chu (1992) where the Rice PDF shape parameter is formulated in terms of the in-phase component of the fluctuation component of the scattered field.

##### c. Randomized circumferential and internal waves added to front interface echo; nulls and attenuation effects.

The complexity is now increased by accounting for circumferential and internally refracting waves so that the echo consists of *all* waves—due to the front interface and circumferential/internal waves. As discussed above, in the case of the smooth penetrable object, these waves will interfere with each other and cause deep nulls in the pattern of echo magnitude versus frequency. When the penetrable object is roughened, the phases of the echoes associated with the various scattering phenomena will correspondingly vary, each in a different manner specific to the respective phenomenon. The total echo (sum of all components) will fluctuate from realization to realization due to the random roughness in a manner broadly similar to that from an impenetrable rough object but with important differences. For example, the center frequency of the nulls from the destructive interference will vary from realization to realization. Because of the steepness of the null (typically 30 dB variation within a narrow band of frequencies), the echo at frequencies near that of the null will fluctuate significantly.

The fluctuations near the null for acoustic scattering by randomly rough elastic cylinders immersed in a fluid have been observed in both numerical (Stanton and Chu, 1992) and experimental (Gurley and Stanton, 1993) studies. In these studies, the shape parameter of the Rice PDF, when fit with the simulations or experimental data, is shown to vary dramatically at frequencies near each null. The shape parameter decreases by as much as two orders of magnitude near the null which corresponds to a similar increase in the degree to which the echo fluctuates.

The roughness not only affects the phase shifts of circumferential and internal waves, but also their magnitude. In one study, the dominant acoustic Lamb-wave (or “plate wave”) of a spherical elastic shell along all meridional paths was randomized due to the roughness. The variability in path length and, hence, phase of this circumferential wave was related to the roughness. The total Lamb wave echo summed from all paths was shown to be attenuated exponentially due to the decrease in coherence of the signal (Stanton , 1998).

Although the above two examples involved scattering of acoustic waves by elastic objects immersed in a fluid, the same principles apply to other sensor systems such as for medical ultrasound or radar applications. In general, when the phases of these different types of scattered signals vary randomly due to roughness, then the shape parameter of the echo PDF will vary significantly near any null in the echo versus frequency pattern. Similarly, any type of circumnavigating or internally refracting wave can experience attenuation due to the different ray paths adding incoherently.

## VII. IN-DEPTH TREATMENT OF ECHO STATISTICS: WITH BEAMPATTERN EFFECTS

Beampattern effects are now added to the above treatment which involved the statistics of echoes from scatterers in the absence of beampattern effects (i.e., equivalent to an omnidirectional beam). Once the scatterers are placed in a directional beam of the sensor system, the echo received by the system becomes modulated by the beampattern. If the location of the scatterer is random, then the modulation is correspondingly random. In this case, the beampattern function is a random variable with a PDF referred to as the beampattern PDF.

The effect of the beampattern on echo statistics can be profound, as discussed in Sec. II C and illustrated (Fig. 4), and accounting for it can be complex. For example, in the simplest case of a point scatterer whose scattering amplitude is delta function distributed, the resultant echo PDF due to the scatterer being randomly located in a directional beam has a trend that is approximately power-law with strong structure superimposed [Fig. 4(d)].

The below treatment begins with a general formula for the echo PDF due to a single scatterer randomly distributed in a beam. The properties of the scatterer, spatial distribution, and beampattern are all arbitrary in this formula. Following the general formula is a progression of examples beginning with the simplest of cases—a point scatterer randomly and uniformly distributed within only the mainlobe of the beam of an axisymmetric transducer (i.e., excluding sidelobes). The analysis then extends to complex scatterers and the entire beam (including sidelobes) and, finally, to the case involving a beampattern from an arbitrary transducer in which the beampattern is not axisymmetric and a spatial distribution of scatterers that is not uniform. After the treatment of single scatterers, the echo statistics associated with multiple scatterers in the beam is presented. In the first set of examples, all scatterers are identical, which is then followed by the more general case of assemblages of scatterers of varying types.

The formulations can be applied to single beam (fixed or scanning) and multi-beam systems, provided that the scatterer(s) are randomly located in the beam. The formulations are general and are not specific to any particular system. Although specific types of signal processing or adaptive beamforming that some systems incorporate are not modeled, the formulations presented herein can serve as a basis for the modeling as discussed in Secs. V A and V B.

### A. Single scatterer randomly located in beam

#### 1. Accounting for beampattern effects in echo PDF

*beampattern*PDF, $ p b ( b )$. Consider now a scatterer whose scattering properties are random (such as a randomly rough and/or randomly oriented elongated scatterer). The scattering amplitude is now a random variable and its magnitude can be described by the PDF $ p s ( | f b s | )$ (Secs. II B and VI B). The magnitude of the echo $ e \u0303$ received by the system is equal to the product of the magnitude of the scattering amplitude $ | f b s |$ and the beampattern $ b ( \theta , \u2009 \varphi )$ [Eq. (4)]. With both of these latter two terms being random variables, then $ e \u0303$ is also a random variable (Sec. IV C 4). Using Eq. (20), the PDF of $ e \u0303$ can be written in terms of $ p S$ and $ p b$ as [Ehrenberg (1972)]

*x*is used to denote $ | f b s |$. The term b ( = $ \u2009 e \u0303$/

*x*) is implicitly the argument of $ p b$. Using the same procedure, $ p e ( e \u0303 )$ can be expressed in an alternate, but equivalent form

*b*) is implicitly the argument of $ p S$.

Equations (36a) and (36b) are equivalent because of the commutative nature of the product of the two random variables in Eq. (20). The limits of ± $\u221e$ in the integral in Eq. (20) are reduced to the ranges [ $ e \u0303$ $\u221e$] and [0 1] in Eqs. (36a) and (36b), respectively, since the values of the beampattern *b* only span the range [0 1] and the corresponding beampattern PDF $ p b$ is zero outside that range. Finally, the range $ 0 \u2264 b \u2264 1$ is used in each equation as they apply to the *entire* beampattern, whereas some applications later will involve only portions of the beampattern, such as for values of the mainlobe only above the highest sidelobe. In those cases, the limits in the integrals are modified accordingly.

These expressions, originally derived by Ehrenberg (1972) [using Eq. (20)] for echo intensity with identical form, are given for echo magnitude (i.e., not intensity) in Ehrenberg (1981) and are also described in reviews in Stanton and Clay (1986) and Ehrenberg (1989). While use of one form over the other [Eq. (36a) vs Eq. (36b)] was not explained in the early papers, it is possible that one form may be more conducive for evaluation, such as in numerical integration (Bhatia , 2015).

The above integral relationship in Eqs. (36a) and (36b) between the echo PDF $ p e ( e \u0303 )$ and the PDFs of the magnitude of the scattering amplitude and beampattern function is completely general, as it applies to an arbitrary stochastic scatterer that is randomly located (at approximately constant range) in an arbitrary beampattern over an arbitrary spatial distribution. The constraint of the scatterers being at approximately constant range is consistent with Eq. (4) as these equations apply to the scatterers distributed within a thin shell at a nearly constant distance from the sensor system. This eliminates the range-dependent effects in analysis of the echo fluctuations.

#### 2. PDF of spatial distribution of scatterer

The beampattern PDF depends not only on the beampattern function $ b ( \theta , \u2009 \varphi )$, but also the PDF, $ p \theta , \varphi ( \theta \u2009 , \varphi ) ,$ of the angular location of the scatterer in the beam. This probability, in combination with the beampattern function, determines the degree to which the echo is randomly modulated by the beampattern. For example, if the scatterer were fixed in the center of the beam, then $ p \theta , \varphi $ is a delta function peaked at $ ( \theta \u2009 , \varphi )$ = (0, 0) and the echo is multiplied by unity for all realizations [Fig. 4(a)]. In the other extreme, if the scatterer were randomly located throughout the entire half-space, then $ p \theta , \varphi $ is finite for all $ \theta \u2009 \u2009 and \u2009 \varphi .$ In this latter case, the echo is randomly modulated by all values of the beampattern resulting in a wide range of echo values, even for a point scatterer [Fig. 4(d)].

Two simple examples are treated here involving the scatterer being randomly and uniformly distributed in a half-plane (2D) and half-space (3D) at approximately a constant range. The 2D and 3D cases apply to geometries in which the sensor system is detecting scatterers that are distributed throughout a thin semicircular arc and a thin hemispherical shell, respectively. The 2D case may apply to geometries where (1) the transducer is a line with a beampattern that only varies in one plane, (2) the sensor system is located within a thin layer of scatterers and is looking along the layer from within that layer, or (3) the system is in a waveguide and long-range echoes only vary with respect to one dimension (such as in the plane parallel to the waveguide boundaries as described in Sec. VIII C). In all of these cases, only a half-plane or half-space are considered as the transducer is assumed to be baffled sufficiently so that there is no “back” radiation.

In each example, $ p \theta , \varphi $ (which reduces to $ p \theta $ in some cases) is calculated using Eq. (8), where $ p X ( x )$ in that equation is used to represent the probability of occurrence of the scatterer per unit volume at angular location $ ( \theta \u2009 , \varphi )$ and *x* represents the volume *V* (Chap. 10 of Medwin and Clay, 1998). The term $ d P X$ $ ( = d P V )$ in Eq. (8) is the differential probability of occurrence of the scatterer in the differential volume *dV* at location $ ( \theta \u2009 , \varphi )$. For the case of the scatterers being located within a thin hemispherical shell of constant radius, which is the 3D geometry in most examples in this tutorial, then *dx = dV* = *r ^{2}* $ sin \u2009 \theta d \theta d \varphi $

*Δr*, where

*r*is the radius of the hemispherical shell, $\theta $ is the spherical polar angle, $\varphi $ is the spherical azimuthal angle, and

*Δr*is the thickness of the shell. For the case of the scatterers being located within a thin-shelled semi-circular arc of constant radius

*r*, which is the 2D geometry sometimes used in this tutorial, then

*dx = dV = r*$ d \theta \Delta w \Delta r$, where $ \Delta w$ is the (thin) width of the circular arc (strip). The total volume of the thin shell in each of these two cases is

*2*$\pi $

*r*

^{2}

*Δr*and $\pi $

*r*$ \Delta w$Δ

*r*, respectively.

For a scatterer uniformly distributed within the volume, then the differential probability of occurrence $ d P V$ per unit differential volume $ d V$ is held constant (i.e., $ d P V / d V$ = *p _{v}* = constant). Using that constraint, and the fact that the integral of $ d P V$ over the total volume that the scatterer can occupy is unity [and, hence, $ d P V / d V$ =

*p*= (total volume)

_{v}^{−1}], then

*p*= ( $\pi $

_{v}*r*$ \Delta w$

*Δr*)

^{−1}and

*p*= (

_{v}*2*$\pi $

*r*)

^{2}Δr^{−1}for the 2D and 3D cases, respectively. Through these changes in variables, Eq. (8) becomes $ d P V$ = $ p \theta d \theta $ and $ d P V$ = $ p \theta , \varphi d \theta d \varphi $ for these two cases, respectively, where expressions for $ p \theta $ and $ p \theta , \varphi $ are given below.

*d*$\theta $ and spanning all values of $ \varphi .$ Since the scatterer is randomly and uniformly distributed across all values of $\varphi $ within each ring, then the probability $ p \theta , \varphi $ only depends upon the area of each ring, which is proportional to $ sin \u2009 \theta $ (which appears in the expression for

*dV*above). Accounting for the uniform distribution across all $\varphi $ [0, $ 2 \pi $] for a given value of $\theta $ yields a factor ( $ 2 \pi $)

^{−1}in Eq. (38).

#### 3. Beampattern PDF for mainlobe only (axisymmetric transducer, uniformly distributed scatterer)

##### a. Exact solution.

##### b. Power law approximation for beampattern PDF.

*k*ensures the integral of Eq. (43) over

_{0}*b*is unity and is given by

*b*approaching unity or $ k a T \u2009 sin \u2009 \theta $ approaching zero (that is, for angles near the center of the beam), the exponent in Eq. (43) is (Chu and Stanton, 2010)

*b*under these limiting conditions. This equation was derived under the assumption that the scatterers are uniformly distributed at approximately constant range within the mainlobe for polar angles in the range $ 0 \u2264 \theta \u2264 \theta S L$ (all $\varphi $) (i.e., a 3D case). Note that the power-law form in Eq. (43) applies to both an intensity-based analysis (Ehrenberg, 1972; Ehrenberg , 1981) and magnitude-based analysis such as in this tutorial and in Chu and Stanton (2010).

#### 4. Beampattern PDF for entire beam (axisymmetric transducer, uniformly distributed scatterer)

*m*in Eq. (46), each value of

*m*corresponds to a portion of the beampattern that is monotonic and is associated with the structure of the mainlobe and sidelobes (Fig. 11).

With the sidelobes accounted for in Eq. (47), the beampattern PDF has significant structure involving singularities (Fig. 12). Each sharp peak in the PDF is associated with the peak of a sidelobe, while the smoothly varying portion with a nearly constant slope at the higher values of *b* is associated with the portion of the mainlobe above the highest sidelobe as discussed above. The beampattern PDF is also shown to vary with beamwidth (i.e., different $ k a T$). The narrower the beam, the more sidelobes are present, which leads to correspondingly more structure in the PDF. There is also some similarity in the occurrence of the sharp peaks as beamwidth is varied.

#### 5. Beampattern PDF for 2D and 3D distribution of scatterers

For any distribution of scatterers containing at least most of the main lobe of the beampattern, the beampattern PDF will generally be qualitatively similar for all distributions of scatterers. Specifically, the PDF will generally have a downward trend, such as with the (approximately) power law illustrated in Figs. 4 and 12. Naturally, there will be some differences associated with the different distributions. As shown in Fig. 4, if the scatterers are only in the main lobe of the beam and do not encounter sidelobes, then the beampattern PDF decreases monotonically with amplitude. Once the sidelobes are encountered, then there is structure superimposed on this decreasing trend. Furthermore, the beampattern PDF will also vary depending upon whether or not the scatterers reside in a half-plane containing the MRA or the entire half space (Fig. 13). While the beampattern PDF is shown to be generally similar between the two cases, there is a singularity associated with the 2D case at the maximum value of amplitude.

Beampattern PDFs, such as those illustrated in Fig. 13, for 2D (half plane) and 3D (half space) distributions of scatterers at approximately constant range can be predicted analytically using Eqs. (37) and (39), respectively, with Eq. (41) used in Eq. (46). However, for the purpose of illustration and comparison with the analytical predictions in Fig. 12, the curves in Fig. 13 were produced through Monte Carlo simulations of Eq. (41) in which random draws of $\theta $ (fixed $\varphi $) and ( $\theta $, $\varphi $) were made for the 2D and 3D cases, respectively.

#### 6. Echo PDF for different types of individual scatterers in axisymmetric beam

The above expressions for the beampattern PDF in Secs. VII A 2–VII A 4 are now incorporated into the general expression in Eq. (36) to calculate the effects of the beampattern on the echo PDF. The examples below give a progression of simple to complex scatterers, located either in the main lobe of the beam only or anywhere in the entire beam (including sidelobes).

##### a. Point scatterer; main lobe only.

*δ*( $ | f b s |$− $ | f p s |$) for the PDF of the scattering amplitude into Eq. (36a) (where

*x*denotes $ | f b s |$) with the upper limit of the integral $ e \u0303 / b S L$, reflecting that the integral only involves values of the beampattern above the highest sidelobe. The inequality involving

*b*given after Eq. (48) explicitly indicates the range of

*b*over which the equation applies, which also corresponds to the range in polar angles, $ 0 \u2264 \theta \u2264 \theta S L$ (all $\varphi $).

This simple example illustrates the significance of the beampattern in echo statistics. With the simplest case of a point scatterer with a constant scattering amplitude, the echo received through the receiver of a directional sensor system involving the main lobe only is approximately *power law distributed* as shown in Eq. (48) [Figs. 4(c) and 12(a)]. As demonstrated below, the echo PDF is further complicated once a more complex scatterer and the entire beampattern (including sidelobes) are used.

##### b. Rayleigh scatterer; mainlobe only.

*b*given after Eq. (49) explicitly indicates the range of

*b*over which the equation applies, which also corresponds to the range in polar angles, $ 0 \u2264 \theta \u2264 \theta S L$ (all $\varphi $).

##### c. Complex scatterer; entire beam.

The analysis is further extended to more complex cases involving a point scatterer, Rayleigh scatterer, randomly oriented smooth prolate spheroid, and randomly oriented rough prolate spheroid, each randomly and uniformly located (one at a time) in a half-plane at approximately constant range involving the entire beam [mainlobe and all sidelobes; $ \u2212 \pi / 2 \u2264 \theta \u2264 \pi / 2$ (fixed $\varphi $)] (Fig. 14). For each of those cases, the PDF of the magnitude of the scattering amplitude *p _{s}*( $ | f b s |$) using

*δ*( $ | f b s |$− $ | f p s |$), Eq. (21), Eq. (33), and Eq. (35), respectively, is inserted into Eq. (36a), where

*x*denotes $ | f b s |$. For the case of using an axisymmetric beampattern due to a circular planar piston transducer, the beampattern PDF given in Eq. (46) is used in Eq. (36a) in which the entire beam is accounted for (0 $ \u2264 b \u2264 1$), while using Eq. (41) for the beampattern function and Eq. (37) for the 2D distribution of scatterers. Also, the prolate spheroids rotate in the plane and Eq. (34) is used for their orientation distribution.

As in Sec. VI in which beampattern effects are not included, predictions using these equations that now incorporate beampattern effects show dependence of echo PDF [and now also PFA through use of Eq. (12)] with type of scatterer (Fig. 14). However, the beampattern significantly alters the shape of the echo PDF over the counterpart cases not involving the beampattern: i.e., the delta function PDF for a point scatterer, Rayleigh PDF for a Rayleigh scatterer, and the PDFs associated with a randomized prolate spheroid as illustrated in Fig. 10. Furthermore, the degree to which the PDFs deviate from the Rayleigh PDF increases once the beampattern is included. As with the examples excluding beampattern effects, the more elongated the scatterers become, the greater the degree to which the echo PDF is non-Rayleigh. As with the PDFs, the slope of the tail of all PFAs depends upon scatterer type. Note that these examples in Fig. 14 which involve a 2D distribution of scatterers are qualitatively similar to the corresponding examples involving 3D distributions in Sec. VII B for the single scatterer (*N* = 1) cases.

#### 7. Beampattern PDF for non-axisymmetric beampattern, non-uniform distribution of scatterer

The above cases involve the simpler examples in which the beampattern is axisymmetric and the scatterer is randomly and uniformly distributed in a half-plane or half-space at approximately constant range. Those examples apply to many important scenarios. However, there are also important cases in which the beampattern may not be axisymmetric (such as a rectangular transducer or mills cross array) and where the location of the scatterer is non-uniformly distributed. In this more general scenario, the beampattern is now a function of two random variables— $ \theta \u2009 and \u2009 \varphi $. Furthermore, the PDF of the angular location of the scatterer is now a function of both of those variables. Below, the beampattern PDF for the most general case of non-axisymmetric beampattern and non-uniform distribution of scatterer is first given, which is followed by the simplified case of a non-axisymmetric beampattern with a uniformly distributed scatterer.

##### a. Non-uniformly distributed scatterer.

*z = b, P*(where “

_{Z}= P_{B}*B*” is the random variable for the beampattern function

*b*), $ x = \theta , \u2009 y = \varphi $, and $ p X , Y ( x , y )$ = $ p \theta , \varphi ( \theta \u2009 , \varphi )$. Differentiating

*P*with respect to $\theta $ and rearranging terms gives the following expression for the beampattern PDF (Ehrenberg, 1972):

_{B}*m*accounts for the fact that the beampattern function is not monotonic. Each segment within the integrand associated with a value of

*m*is monotonic. These segments are related to the regions

*D*in Eq. (18) where the function

_{Z}*b*has the same value for multiple values of $\theta $. This is illustrated specifically for the case of the beampattern function in Fig. 11.

##### b. Uniformly distributed scatterer.

Given the complexity associated with the asymmetry of the sensor beam, both of the above two equations must generally be evaluated numerically. Equation (51) has been evaluated to predict the echo statistics associated with a rectangular transducer in which the one-way beamwidths in the two orthogonal planes were 5° and 20° (Stanton , 2015). The beampattern PDF of the two-way beampattern, as illustrated in Fig. 2 in that paper, is qualitatively similar to the ones illustrated in this tutorial in that the PDF trends toward smaller values as the beampattern value increases. However, the structure in the beampattern PDF associated with the rectangular transducer is much different than that associated with the circular transducer in this tutorial because of the lack of axial symmetry in the beampattern. Specifically, the sharp spikes shown in Fig. 12(a) of this tutorial that are associated with the sidelobes for a circular transducer are much larger in magnitude, but fewer in number, than the corresponding spikes in the beampattern PDF associated with the rectangular transducer. The beampattern PDF for the rectangular transducer was used in Stanton (2015) in interpreting experimental data, as summarized in Sec. III A 1 of this tutorial.

### B. Multiple identical scatterers randomly located in beam

In this case, there are two or more scatterers present at the same time, each randomly, uniformly, and independently distributed in the beam at approximately constant range. The transmit signals are long enough so that the echoes from all scatterers are assumed to completely overlap. The scatterers are “identical” in that they possess the same statistical properties. For example, the magnitude of the scattering amplitude of each scatterer could be Rayleigh distributed with the same mean scattering cross section. Or, each scatterer could be a randomly rough, randomly oriented prolate spheroid with the same mean dimensions (and, hence, the same mean scattering cross section). Although the statistical properties are the same, the scattering amplitudes of the scatterers are generally different from each other for any given ping or realization since they are statistically independent of each other.

As discussed in Sec. IV C 5, there are various methods to calculate the echo statistics in this case in which the sum of multiple random variables (i.e., echoes from multiple scatterers) is calculated. The methods range from closed-form analytical to pure numerical approaches involving Monte Carlo simulations of summations of phasors. Because of its generality, the latter case is used in the below analysis.

In this simple “phasor summation” method, Eq. (6) is used in which a phasor (the summand) is calculated for each scatterer and each realization. For each realization, the phasors are added together coherently to form the total echo as measured by the sensor system. The echo PDF is estimated through forming a histogram of the total echo magnitude through the binning method or using the kernel density estimation (KDE) method described in Sec. IV B 2. The phase shift term $ \Delta i$ varies randomly and uniformly in the range [0 2 $\pi $], reflecting the random location of the scatterers (range-wise) and high frequencies (short wavelengths) of the signal. All scatterers are distributed within a thin hemispherical shell so that there are no significant differences in the range-dependence of the losses due to spreading and absorption. The magnitude of the echo from each individual scatterer is given by Eq. (7). As discussed above, since the scattering amplitude and location of the scatterer in the beam are random variables (leading to the beampattern function being a random variable), then this individual echo magnitude is a random variable as well.

Three sets of examples are investigated using three scatterer types from above—a point scatterer, Rayleigh scatterer, and randomly oriented rough prolate spheroid. For each phasor, the randomized terms (scattering amplitude, beampattern function, and phase) are randomly drawn from numerically generated statistical distributions. For the point scatterer, the magnitude of the scattering amplitude is simply a constant in Eq. (7). For the other two scatterers, the magnitude of the scattering amplitude is determined using random draws from numerically generated random variables. The magnitude of the scattering amplitude of the Rayleigh scatterer for each realization was randomly drawn from a numerically generated Rayleigh random variable [whose PDF is given in Eq. (21)]. Generating the magnitude of the scattering amplitude of the randomly oriented rough prolate spheroid began by randomly drawing an orientation angle *β*, which was then used to calculate the magnitude of the scattering amplitude of the smooth prolate spheroid from Eq. (30). That value, in turn, was multiplied by a Rayleigh distributed random variable (which was randomly drawn using the same method as for the Rayleigh scatterer). Note that the prolate spheroid calculations could have used Eq. (35) directly to describe the statistics of the random variable scattering amplitude, from which a random draw could have been made. However, for the purposes of illustration, the scattering was described from the beginning (random draw of orientation angle), which would be the process used for a more general scatterer for which there is not a closed-form solution.

For each realization, the (axisymmetric) beampattern is calculated for a random location (polar angle $\theta $) with the scatterer angular (location) distribution PDF of $ sin \u2009 \theta \u2009$ from Eq. (39), and the phase shift is sampled from a uniform distribution [0 2 $\pi $]. Also, for each realization, each of the above random variables is generated by employing inverse transform sampling, in which samples of any probability distribution is generated at random through its CDF (Devroye, 1986). Using these terms, the summand in Eq. (6) is calculated and summed over the *N* scatterers for each realization. The process is repeated for millions of realizations (typically 10^{7}) and a histogram is formed representing the statistics of the magnitude of the echo as received by the sensor system.

In each example, the echo PDFs are shown to vary significantly with type and number of scatterers (Figs. 15–17). As with the case of a single scatterer in the beam, the PDF is significantly different than when beampattern effects are not accounted for. As expected, as the number of scatterers increases, the PDFs approach the Rayleigh PDF. In general, for a small number of scatterers, the echo PDF deviates significantly from the Rayleigh PDF, both in the small and large echo magnitude regions of the PDF. Echo PFAs are also illustrated, which also show significant dependences upon type and number of scatterer. The degree to which the PDF deviates from a Rayleigh PDF also varies with beamwidth (Fig. 18). For a fixed number of scatterers, the narrower the beam, the greater the deviation from a Rayleigh PDF. Note that the “noisy” characteristic in portions of some of the plots of PDF in Figs. 15–18 is due to the relatively low number of realizations in the Monte Carlo simulations that were within those particular log-spaced magnitude bins (i.e., when both $ p e$ and $ e \u0303 / \u27e8 e \u0303 2 \u27e9 1 / 2$ are low).

It is important to emphasize that the number of scatterers *N* is distributed uniformly in a thin hemispherical shell in a half-space. For a narrow beamwidth, the mainlobe occupies just a small fraction of that half-space. Since the mainlobe weights the echoes from the individual scatterers much more than the sidelobes, the scatterers within the mainlobe will generally contribute the most to the signal. In essence, the signal will be dominated by the *effective* number of scatterers—those that are within the mainlobe. The average number of scatterers within the mainlobe is indicated in parentheses after the value of *N* in the legends to these and later figures and are also given in Table I.

Beamwidth (−3 to −3 dB) (deg.) . | ka
. _{T} | N (half space)
. | Average number within mainlobe . |
---|---|---|---|

1 | 132.74 | 100 | 0.0417 |

3 | 44.2511 | 1 | 0.00375 |

3 | 44.2511 | 10 | 0.0375 |

3 | 44.2511 | 25 | 0.0937 |

3 | 44.2511 | 100 | 0.375 |

3 | 44.2511 | 250 | 0.937 |

3 | 44.2511 | 1000 | 3.75 |

3 | 44.2511 | 2500 | 9.37 |

10 | 13.2907 | 100 | 4.13 |

20 | 6.6707 | 100 | 16.0 |

Beamwidth (−3 to −3 dB) (deg.) . | ka
. _{T} | N (half space)
. | Average number within mainlobe . |
---|---|---|---|

1 | 132.74 | 100 | 0.0417 |

3 | 44.2511 | 1 | 0.00375 |

3 | 44.2511 | 10 | 0.0375 |

3 | 44.2511 | 25 | 0.0937 |

3 | 44.2511 | 100 | 0.375 |

3 | 44.2511 | 250 | 0.937 |

3 | 44.2511 | 1000 | 3.75 |

3 | 44.2511 | 2500 | 9.37 |

10 | 13.2907 | 100 | 4.13 |

20 | 6.6707 | 100 | 16.0 |

The echo statistics vary with beamwidth even though *N* is constant because the effective number of scatterers varies with beamwidth—that is, for fixed *N*, the number of scatterers within the mainlobe that contributes to the higher values of the scattering decreases with decreasing beamwidth (Fig. 18). For example, for *N* = 100, the number of scatterers within a 1°-wide mainlobe is, on average, equal to 0.0417 (Table I). Before beampattern effects are included, the sum of 100 independent, random phase signals would be very close to Rayleigh distributed (Fig. 6). However, with such a narrow beam with only a small fraction of the 100 scatterers within the mainlobe, the resultant echo PDF will be strongly non-Rayleigh.

### C. Multiple scatterers of different types and sizes

Once there is more than one type or size of scatterer present, the echo PDF will not only depend upon the type and number of scatterers, and sensor beamwidth, as in the above examples in Sec. VII B which involved one type of scatterer at a time, but will also depend upon the variability of values of mean scattering cross section and variability of numbers of scatterers across the different types of scatterers. Furthermore, the echo PDF will depend upon the differences in spatial distributions of the different types of scatterers. For example, consider two distinct cases: (1) a *split aggregation* in which scatterers of only one type reside in “monotype” patches and each echo contains contributions from only one type of scatterer and (2) an *interspersed aggregation* in which the different types of scatterers are uniformly interspersed with each other so that each echo contains contributions from all of the types of scatterers (Fig. 19). In each case, the beam is much narrower than both the patch and analysis window containing the scatterers and the echoes from many pings or independent realizations are observed within each window (Fig. 19). The echo statistics will generally be significantly different for the two distributions, even when the same set of scatterers of varying types are involved in each distribution (Lee and Stanton, 2014). And, for each type of distribution, the echo PDF will vary as the number and mean scattering cross section of each scatterer type are varied.

Below, the echo PDF is formulated for each type of spatial distribution—first for cases involving two types of scatterers, then generalizing to an arbitrary number of scatterer types. For the split aggregation, a mixture PDF (to be defined below) involving separate phasor summations for the monotype patches is used and, for the interspersed aggregation, the phasor summation method is used to form the echo PDF from the echoes from all types of scatterers within the analysis window. Examples of echo PDFs from each type of distribution are given in Sec. VII C 3 for the simple case involving two different sizes of scatterers, “weak” (type A) and “strong” (type B), for a range of numbers and sizes of scatterers. Both types of scatterers are Rayleigh scatterers, but with different mean scattering cross sections.

#### 1. Split aggregation of type A and B scatterers—mixture PDF

#### 2. Interspersed aggregation of type A and B scatterers—coherent phasor sum

*i*th scatterer of the

*k*th type as received through the sensor system is

*i*th scatterer of the

*k*th type, respectively. There are

*N*and

*M*type A and B scatterers, respectively. The term

*k*corresponds to type A or type B scatterer in this example. Once $ e \u0303$ is calculated from Eq. (54) for many realizations, the echo PDF is formed.

#### 3. Comparisons between echo PDFs from split and interspersed two-type aggregations

*rms*magnitude of the scattering amplitudes of the strong (“

*S*”) to weak (“

*W*”) scatterers is given by

*rms*magnitude of the scattering amplitude of the denoted scatterer type [this notation is chosen to be consistent with that of Lee and Stanton (2014) and is not to be confused with $ \lambda R$ of Eqs. (21)–(23) of this tutorial].

With the number of weak scatterers fixed at 2500, the number of strong scatterers is varied over the range 25–2500 for two values of $ r S W$ (5 and 20). The resultant echo PDFs are shown to vary in shape over all combinations of these parameters (Figs. 20 and 21). All PDFs deviate significantly from the Rayleigh PDF. In each example, the tail of the PDF is elevated above the Rayleigh PDF. The degree to which the tail is elevated is especially pronounced for the larger ratio of strong-to-weak scattering amplitude ( $ r S W$* = *20) in both types of aggregations and for the cases involving fewer numbers of strong scatterers in the interspersed aggregations. This is consistent with the intuition that the strong scatterers dominate the echo and can cause the echo to be non-Rayleigh when they are small in number. Also, very importantly, the PDFs with the same parameters, but different spatial distribution (split- and interspersed aggregation), are significantly different from each other. Note that, as with some of the previous simulations, the “noisy” characteristic in portions of some of the plots of PDF in Figs. 20 and 21 are due to the relatively low number of realizations in the Monte Carlo simulations that were within those particular log-spaced magnitude bins (i.e., when both $ p e$ and $ e \u0303 / \u27e8 e \u0303 2 \u27e9 1 / 2$ are low).