Multiple-scattering in radar systems: a review

Although extensively studied within the lidar community, the multiple scattering phenomenon has always been considered a rare curiosity by radar meteorologists. Up to few years ago its appearance has only been associated with two- or three-body-scattering features (e.g. hail flares and mirror images) involving highly reflective surfaces. Recent atmospheric research aimed at better understanding of the water cycle and the role played by clouds and precipitation in affecting the Earth's climate has driven the deployment of high frequency radars in space. Examples are the TRMM 13.5 GHz, the CloudSat 94 GHz, the upcoming EarthCARE 94 GHz, and the GPM dual 13-35 GHz radars. These systems are able to detect the vertical distribution of hydrometeors and thus provide crucial feedbacks for radiation and climate studies. The shift towards higher frequencies increases the sensitivity to hydrometeors, improves the spatial resolution and reduces the size and weight of the radar systems. On the other hand, higher frequency radars are affected by stronger extinction, especially in the presence of large precipitating particles (e.g. raindrops or hail particles), which may eventually drive the signal below the minimum detection threshold. In such circumstances the interpretation of the radar equation via the single scattering approximation may be problematic. Errors will be large when the radiation emitted from the radar after interacting more than once with the medium still contributes substantially to the received power. This is the case if the transport mean-free-path becomes comparable with the instrument footprint (determined by the antenna beam-width and the platform altitude). This situation resembles to what has already been experienced in lidar observations, but with a predominance of wide- versus small-angle scattering events. At millimeter wavelengths, hydrometeors diffuse radiation rather isotropically compared to the visible or near infrared region where scattering is predominantly in the forward direction. A complete understanding of radiation transport modeling and data analysis methods under wide-angle multiple scattering conditions is mandatory for a correct interpretation of echoes observed by space-borne millimeter radars. This paper reviews the status of research in this field. Different numerical techniques currently implemented to account for higher order scattering are reviewed and their weaknesses and strengths highlighted. Examples of simulated radar backscattering profiles are provided with particular emphasis given to situations in which the multiple scattering contributions become comparable or overwhelm the single scattering signal. We show evidences of multiple scattering effects from air-borne and from CloudSat observations, i.e. unique signatures which cannot be explained by single scattering theory. Ideas how to identify and tackle the multiple scattering effects are discussed. Finally perspectives and suggestions for future work are outlined. This work represents a reference-guide for studies focused at modeling the radiation transport and at interpreting data from high frequency space-borne radar systems that probe highly opaque scattering media such as thick ice clouds or precipitating clouds.

Understanding three-dimensional effects in polarized observations with the ground-based ADMIRARI radiometer during the CHUVA campaign

Measurements of down-welling microwave radiation from raining clouds performed with the Advanced Microwave Radiometer for Rain Identification (ADMIRARI) radiometer at 10.7-21-36.5 GHz during the Global Precipitation Measurement Ground Validation ""Cloud processes of the main precipitation systems in Brazil: A contribution to cloud resolving modeling and to the Global Precipitation Measurement"" (CHUVA) campaign held in Brazil in March 2010 represent a unique test bed for understanding three-dimensional (3D) effects in microwave radiative transfer processes. While the necessity of accounting for geometric effects is trivial given the slant observation geometry (ADMIRARI was pointing at a fixed 30 elevation angle), the polarization signal (i.e., the difference between the vertical and horizontal brightness temperatures) shows ubiquitousness of positive values both at 21.0 and 36.5 GHz in coincidence with high brightness temperatures. This signature is a genuine and unique microwave signature of radiation side leakage which cannot be explained in a 1D radiative transfer frame but necessitates the inclusion of three-dimensional scattering effects. We demonstrate these effects and interdependencies by analyzing two campaign case studies and by exploiting a sophisticated 3D radiative transfer suited for dichroic media like precipitating clouds.

Cooperative scattering by cold atoms

We have studied the interplay between disorder and cooperative scattering for the single scattering limit in the presence of a driving laser. Analytical results have been derived and we have observed cooperative scattering effects in a variety of experiments, ranging from thermal atoms in an optical dipole trap, atoms released from a dark MOT and atoms in a BEC, consistent with our theoretical predictions.

Shower approach in the simulation of ion scattering from solids

An efficient approach for the simulation of ion scattering from solids is proposed. For every encountered atom, we take multiple samples of its thermal displacements among those which result in scattering with high probability to finally reach the detector. As a result, the detector is illuminated by intensive “showers,” where each event of detection must be weighted according to the actual probability of the atom displacement. The computational cost of such simulation is orders of magnitude lower than in the direct approach, and a comprehensive analysis of multiple and plural scattering effects becomes possible. We use this method for two purposes. First, the accuracy of the approximate approaches, developed mainly for ion-beam structural analysis, is verified. Second, the possibility to reproduce a wide class of experimental conditions is used to analyze some basic features of ion-solid collisions: the role of double violent collisions in low-energy ion scattering; the origin of the “surface peak” in scattering from amorphous samples; the low-energy tail in the energy spectra of scattered medium-energy ions due to plural scattering; and the degradation of blocking patterns in two-dimensional angular distributions with increasing depth of scattering. As an example of simulation for ions of MeV energies, we verify the time reversibility for channeling and blocking of 1-MeV protons in a W crystal. The possibilities of analysis that our approach offers may be very useful for various applications, in particular, for structural analysis with atomic resolution.

Ultra-long raman laser with a feedback based on the Rayleigh scattering

This paper reports the Rayleigh scattering effects in ultra-long Raman fibre laser. It has been found that in a long fibre cavity (-100 km) the distributed feedback due to Rayleigh back scattering at propagation of light between fibre Bragg grating reflectors may be comparable with the lumped feedback provided by the FBG itself. As a result, Raman lasing in the fibre span limited by lumped (FBG) reflector at one side only appears possible due to significant reflection from the RS-based "random" distributed mirror at the other side. Thus, it concludes that a distributed Rayleigh scattering "random" mirror can form a cavity together with a single FBG spliced to the opposite cavity end.

Elastic and inelastic scattering of fast neutrons in fusion reactor materials

In this work, the angular distributions for elastic and. inelastic scattering of fast neutrons in fusion .reactor materials have been studied. Lithium and lead material are likely to be common components of fusion reactor wall configuration design. The measurements were performed using an associated particle time-of- flight technique. The 14 and 14.44 Mev neutrons were produced by the T(d,n} 4He reaction with deuterons being accelerated in a 150kev SAMES type J accelerator at ASTON and in.the 3. Mev DYNAMITRON at the Joint Radiation Centre, Birmingham respectively. The associated alpha-particles and fast. neutrons were detected.by means of a plastic scintillator mounted on a fast focused photomultiplier tube. The samples used were extended flat plates of thicknesses up to 0.9 mean-free-path for Lithium and 1.562 mean-free-path for Lead. The differential elastic scattering cross-sections were measured for 14 Mev neutrons for various thicknesses of Lithium and Lead in the angular range from zero to; 90º. In addition, the angular distributions of elastically scattered 14,.44 Mev .neutrons from Lithium samples were studied in the same angular range. Inelastic scattering to the 4.63 Mev state in 7Li and the 2.6 Mev state, and 4.1 Mev state in 208Pb have:been :measured.The results are compared to ENDF/B-IV data files and to previous measurements. For the Lead samples the differential neutron scattering:cross-sections for discrete 3 Mev ranges and the angular distributions were measured. The increase in effective cross-section due to multiple scattering effects,as the sample thickness increased:was found to be predicted by the empirical .relation ....... A good fit to the exoerimental data was obtained using the universal constant............ The differential elastic scattering cross-section data for thin samples of Lithium and Lead were analyzed in terms of optical model calculations using the. computer code. RAROMP. Parameter search procedures produced good fits to the·cross-sections. For the case of thick samples of Lithium and Lead, the measured angular distributions of :the scattered neutrons were compared to the predictions of the continuous slowing down model.

Gamma rays associated with neutron scattering in fusion reactor materials

The gamma-rays produced by the inelastic scattering of 14 MeV neutrons. in fusion reactor materials have been studied using a gamma-ray spectrometer employing a sodium iodide scintillation detector. The source neutrons are produced by the T(d,n)4He reaction using the SAMES accelerator at the University of Aston in Birmingham. In order to eliminate the large gamma-ray background and neutron signal due to the sensitivity of the sodium iodide detector to neutrons, the gamma-ray detector is heavily shielded and is used together with a particle time of flight discrimination system based on the associated particle time of flight method. The instant of production of a source neutron is determined by detecting the associated alpha-particle enabling discrimination between the neutrons and gamma-rays by their different time of flight times. The electronic system used for measuring the time of flight of the neutrons and gamrna-rays over the fixed flight path is described. The materials studied in this work were Lithium and Lead because of their importance as fuel breeding and shielding materials in conceptual fusion reactor designs. Several sample thicknesses were studied to determine the multiple scattering effects. The observed gamma-ray spectra from each sample at several scattering angles in the angular range Oº - 90° enabled absolute differential gamma-ray production cross-sections and angular distributions of the resolved gamma-rays from Lithium to be measured and compared with published data. For the Lead sample, the absolute differential gamma-ray production cross-sections for discrete 1 MeV ranges and the angular distributions were measured. The measured angular distributions of the present work and those on Iron from previous work are compared to the predictions of the Monte Carlo programme M.O.R.S.E. Good agreement was obtained between the experimental results and the theoretical predictions. In addition an empirical relation has been constructed which describes the multiple scattering effects by a single parameter and is capable of predicting the gamma-ray production cross-sections for the materials to an accuracy of ± 25%.

Ultra-long raman laser with a feedback based on the Rayleigh scattering

This paper reports the Rayleigh scattering effects in ultra-long Raman fibre laser. It has been found that in a long fibre cavity (-100 km) the distributed feedback due to Rayleigh back scattering at propagation of light between fibre Bragg grating reflectors may be comparable with the lumped feedback provided by the FBG itself. As a result, Raman lasing in the fibre span limited by lumped (FBG) reflector at one side only appears possible due to significant reflection from the RS-based "random" distributed mirror at the other side. Thus, it concludes that a distributed Rayleigh scattering "random" mirror can form a cavity together with a single FBG spliced to the opposite cavity end.

Nonlinear mixture model for hyperspectral unmixing

Proceedings of International Conference - SPIE 7477, Image and Signal Processing for Remote Sensing XV - 28 September 2009

Extracting temporal and spatial distributions information about marine mucilage phenomenon based on Modis satellite images; a case study of the Tyrrhenian and the Adriatic Sea, 2010-2012

Dissertation submitted in partial fulfillment of the requirements for the Degree of Master of Science in Geospatial Technologies.

On the way to solve lighting problems in underwater imaging

A major obstacle to processing images of the ocean floor comes from the absorption and scattering effects of the light in the aquatic environment. Due to the absorption of the natural light, underwater vehicles often require artificial light sources attached to them to provide the adequate illumination. Unfortunately, these flashlights tend to illuminate the scene in a nonuniform fashion, and, as the vehicle moves, induce shadows in the scene. For this reason, the first step towards application of standard computer vision techniques to underwater imaging requires dealing first with these lighting problems. This paper analyses and compares existing methodologies to deal with low-contrast, nonuniform illumination in underwater image sequences. The reviewed techniques include: (i) study of the illumination-reflectance model, (ii) local histogram equalization, (iii) homomorphic filtering, and, (iv) subtraction of the illumination field. Several experiments on real data have been conducted to compare the different approaches

Inversion of crosshole seismic data in heterogenous environments : Comparison of waveform and ray-based approaches

High-resolution tomographic imaging of the shallow subsurface is becoming increasingly important for a wide range of environmental, hydrological and engineering applications. Because of their superior resolution power, their sensitivity to pertinent petrophysical parameters, and their far reaching complementarities, both seismic and georadar crosshole imaging are of particular importance. To date, corresponding approaches have largely relied on asymptotic, ray-based approaches, which only account for a very small part of the observed wavefields, inherently suffer from a limited resolution, and in complex environments may prove to be inadequate. These problems can potentially be alleviated through waveform inversion. We have developed an acoustic waveform inversion approach for crosshole seismic data whose kernel is based on a finite-difference time-domain (FDTD) solution of the 2-D acoustic wave equations. This algorithm is tested on and applied to synthetic data from seismic velocity models of increasing complexity and realism and the results are compared to those obtained using state-of-the-art ray-based traveltime tomography. Regardless of the heterogeneity of the underlying models, the waveform inversion approach has the potential of reliably resolving both the geometry and the acoustic properties of features of the size of less than half a dominant wavelength. Our results do, however, also indicate that, within their inherent resolution limits, ray-based approaches provide an effective and efficient means to obtain satisfactory tomographic reconstructions of the seismic velocity structure in the presence of mild to moderate heterogeneity and in absence of strong scattering. Conversely, the excess effort of waveform inversion provides the greatest benefits for the most heterogeneous, and arguably most realistic, environments where multiple scattering effects tend to be prevalent and ray-based methods lose most of their effectiveness.

Seismic characterization of heterogeneous sedimentary and fractured rocks based on Biot's theory of poroelasticity

Understanding and quantifying seismic energy dissipation, which manifests itself in terms of velocity dispersion and attenuation, in fluid-saturated porous rocks is of considerable interest, since it offers the perspective of extracting information with regard to the elastic and hydraulic rock properties. There is increasing evidence to suggest that wave-induced fluid flow, or simply WIFF, is the dominant underlying physical mechanism governing these phenomena throughout the seismic, sonic, and ultrasonic frequency ranges. This mechanism, which can prevail at the microscopic, mesoscopic, and macroscopic scale ranges, operates through viscous energy dissipation in response to fluid pressure gradients and inertial effects induced by the passing wavefield. In the first part of this thesis, we present an analysis of broad-band multi-frequency sonic log data from a borehole penetrating water-saturated unconsolidated glacio-fluvial sediments. An inherent complication arising in the interpretation of the observed P-wave attenuation and velocity dispersion is, however, that the relative importance of WIFF at the various scales is unknown and difficult to unravel. An important generic result of our work is that the levels of attenuation and velocity dispersion due to the presence of mesoscopic heterogeneities in water-saturated unconsolidated clastic sediments are expected to be largely negligible. Conversely, WIFF at the macroscopic scale allows for explaining most of the considered data while refinements provided by including WIFF at the microscopic scale in the analysis are locally meaningful. Using a Monte-Carlo-type inversion approach, we compare the capability of the different models describing WIFF at the macroscopic and microscopic scales with regard to their ability to constrain the dry frame elastic moduli and the permeability as well as their local probability distribution. In the second part of this thesis, we explore the issue of determining the size of a representative elementary volume (REV) arising in the numerical upscaling procedures of effective seismic velocity dispersion and attenuation of heterogeneous media. To this end, we focus on a set of idealized synthetic rock samples characterized by the presence of layers, fractures or patchy saturation in the mesocopic scale range. These scenarios are highly pertinent because they tend to be associated with very high levels of velocity dispersion and attenuation caused by WIFF in the mesoscopic scale range. The problem of determining the REV size for generic heterogeneous rocks is extremely complex and entirely unexplored in the given context. In this pilot study, we have therefore focused on periodic media, which assures the inherent self- similarity of the considered samples regardless of their size and thus simplifies the problem to a systematic analysis of the dependence of the REV size on the applied boundary conditions in the numerical simulations. Our results demonstrate that boundary condition effects are absent for layered media and negligible in the presence of patchy saturation, thus resulting in minimum REV sizes. Conversely, strong boundary condition effects arise in the presence of a periodic distribution of finite-length fractures, thus leading to large REV sizes. In the third part of the thesis, we propose a novel effective poroelastic model for periodic media characterized by mesoscopic layering, which accounts for WIFF at both the macroscopic and mesoscopic scales as well as for the anisotropy associated with the layering. Correspondingly, this model correctly predicts the existence of the fast and slow P-waves as well as quasi and pure S-waves for any direction of wave propagation as long as the corresponding wavelengths are much larger than the layer thicknesses. The primary motivation for this work is that, for formations of intermediate to high permeability, such as, for example, unconsolidated sediments, clean sandstones, or fractured rocks, these two WIFF mechanisms may prevail at similar frequencies. This scenario, which can be expected rather common, cannot be accounted for by existing models for layered porous media. Comparisons of analytical solutions of the P- and S-wave phase velocities and inverse quality factors for wave propagation perpendicular to the layering with those obtained from numerical simulations based on a ID finite-element solution of the poroelastic equations of motion show very good agreement as long as the assumption of long wavelengths remains valid. A limitation of the proposed model is its inability to account for inertial effects in mesoscopic WIFF when both WIFF mechanisms prevail at similar frequencies. Our results do, however, also indicate that the associated error is likely to be relatively small, as, even at frequencies at which both inertial and scattering effects are expected to be at play, the proposed model provides a solution that is remarkably close to its numerical benchmark. -- Comprendre et pouvoir quantifier la dissipation d'énergie sismique qui se traduit par la dispersion et l'atténuation des vitesses dans les roches poreuses et saturées en fluide est un intérêt primordial pour obtenir des informations à propos des propriétés élastique et hydraulique des roches en question. De plus en plus d'études montrent que le déplacement relatif du fluide par rapport au solide induit par le passage de l'onde (wave induced fluid flow en anglais, dont on gardera ici l'abréviation largement utilisée, WIFF), représente le principal mécanisme physique qui régit ces phénomènes, pour la gamme des fréquences sismiques, sonique et jusqu'à l'ultrasonique. Ce mécanisme, qui prédomine aux échelles microscopique, mésoscopique et macroscopique, est lié à la dissipation d'énergie visqueuse résultant des gradients de pression de fluide et des effets inertiels induits par le passage du champ d'onde. Dans la première partie de cette thèse, nous présentons une analyse de données de diagraphie acoustique à large bande et multifréquences, issues d'un forage réalisé dans des sédiments glaciaux-fluviaux, non-consolidés et saturés en eau. La difficulté inhérente à l'interprétation de l'atténuation et de la dispersion des vitesses des ondes P observées, est que l'importance des WIFF aux différentes échelles est inconnue et difficile à quantifier. Notre étude montre que l'on peut négliger le taux d'atténuation et de dispersion des vitesses dû à la présence d'hétérogénéités à l'échelle mésoscopique dans des sédiments clastiques, non- consolidés et saturés en eau. A l'inverse, les WIFF à l'échelle macroscopique expliquent la plupart des données, tandis que les précisions apportées par les WIFF à l'échelle microscopique sont localement significatives. En utilisant une méthode d'inversion du type Monte-Carlo, nous avons comparé, pour les deux modèles WIFF aux échelles macroscopique et microscopique, leur capacité à contraindre les modules élastiques de la matrice sèche et la perméabilité ainsi que leur distribution de probabilité locale. Dans une seconde partie de cette thèse, nous cherchons une solution pour déterminer la dimension d'un volume élémentaire représentatif (noté VER). Cette problématique se pose dans les procédures numériques de changement d'échelle pour déterminer l'atténuation effective et la dispersion effective de la vitesse sismique dans un milieu hétérogène. Pour ce faire, nous nous concentrons sur un ensemble d'échantillons de roches synthétiques idéalisés incluant des strates, des fissures, ou une saturation partielle à l'échelle mésoscopique. Ces scénarios sont hautement pertinents, car ils sont associés à un taux très élevé d'atténuation et de dispersion des vitesses causé par les WIFF à l'échelle mésoscopique. L'enjeu de déterminer la dimension d'un VER pour une roche hétérogène est très complexe et encore inexploré dans le contexte actuel. Dans cette étude-pilote, nous nous focalisons sur des milieux périodiques, qui assurent l'autosimilarité des échantillons considérés indépendamment de leur taille. Ainsi, nous simplifions le problème à une analyse systématique de la dépendance de la dimension des VER aux conditions aux limites appliquées. Nos résultats indiquent que les effets des conditions aux limites sont absents pour un milieu stratifié, et négligeables pour un milieu à saturation partielle : cela résultant à des dimensions petites des VER. Au contraire, de forts effets des conditions aux limites apparaissent dans les milieux présentant une distribution périodique de fissures de taille finie : cela conduisant à de grandes dimensions des VER. Dans la troisième partie de cette thèse, nous proposons un nouveau modèle poro- élastique effectif, pour les milieux périodiques caractérisés par une stratification mésoscopique, qui prendra en compte les WIFF à la fois aux échelles mésoscopique et macroscopique, ainsi que l'anisotropie associée à ces strates. Ce modèle prédit alors avec exactitude l'existence des ondes P rapides et lentes ainsi que les quasis et pures ondes S, pour toutes les directions de propagation de l'onde, tant que la longueur d'onde correspondante est bien plus grande que l'épaisseur de la strate. L'intérêt principal de ce travail est que, pour les formations à perméabilité moyenne à élevée, comme, par exemple, les sédiments non- consolidés, les grès ou encore les roches fissurées, ces deux mécanismes d'WIFF peuvent avoir lieu à des fréquences similaires. Or, ce scénario, qui est assez commun, n'est pas décrit par les modèles existants pour les milieux poreux stratifiés. Les comparaisons des solutions analytiques des vitesses des ondes P et S et de l'atténuation de la propagation des ondes perpendiculaires à la stratification, avec les solutions obtenues à partir de simulations numériques en éléments finis, fondées sur une solution obtenue en 1D des équations poro- élastiques, montrent un très bon accord, tant que l'hypothèse des grandes longueurs d'onde reste valable. Il y a cependant une limitation de ce modèle qui est liée à son incapacité à prendre en compte les effets inertiels dans les WIFF mésoscopiques quand les deux mécanismes d'WIFF prédominent à des fréquences similaires. Néanmoins, nos résultats montrent aussi que l'erreur associée est relativement faible, même à des fréquences à laquelle sont attendus les deux effets d'inertie et de diffusion, indiquant que le modèle proposé fournit une solution qui est remarquablement proche de sa référence numérique.

An experimental investigation of dechanneling in copper single crystals

This investigation comprises a comparison of experimental and theoretical dechanneling of MeV protons in copper single crystals. Dechanneling results when an ion's transverse energy increases to the value where the ion can undergo small impact parameter collisions with individual atoms. Depth dependent dechanneling rates were determined as functions of lattice temperature, ion beam energy and crystal axis orientation. Ion beam energies were IMeV and 2MeV,temperatures ranged from 35 K to 280 K and the experiment was carried out along both the (lOa) and <110) axes. Experimental data took the form of aligned and random Rutherford backscattered energy spectra. Dechanneling rates were extracted from these spectra using a single scattering theory that took explicit account of the different stopping powers experienced by channeled and dechanneled ions and also included a correction factor to take into account multiple scattering effects along the ion's trajectory. The assumption of statistical equilibrium and small angle scattering of the channeled ions allows a description of dechanneling in terms of the solution of a diffusion like equation which contains a so called diffusion function. The diffusion function is shown to be related to the increase in average transverse energy. Theoretical treatments of increase in average transverse energy due to collisions of projectiles with channel electrons and thermal perturbations in the lattice potential are reviewed. Using the diffusion equation and the electron density in the channel centre as a fitting parameter dechanneling rates are extracted. Excellent agreement between theory and experiment has been demonstrated. Electron densities determined in the fitting procedure appear to be realistic. The surface parameters show themselves to be good indicators of the quality of the crystal.

On the way to solve lighting problems in underwater imaging

A major obstacle to processing images of the ocean floor comes from the absorption and scattering effects of the light in the aquatic environment. Due to the absorption of the natural light, underwater vehicles often require artificial light sources attached to them to provide the adequate illumination. Unfortunately, these flashlights tend to illuminate the scene in a nonuniform fashion, and, as the vehicle moves, induce shadows in the scene. For this reason, the first step towards application of standard computer vision techniques to underwater imaging requires dealing first with these lighting problems. This paper analyses and compares existing methodologies to deal with low-contrast, nonuniform illumination in underwater image sequences. The reviewed techniques include: (i) study of the illumination-reflectance model, (ii) local histogram equalization, (iii) homomorphic filtering, and, (iv) subtraction of the illumination field. Several experiments on real data have been conducted to compare the different approaches