Vegetation/Forest Effects in Microwave Remote Sensing of Soil Moisture

This thesis includes (1) the distorted Born approximation (DBA) and an improved coherent model for vegetation-covered surfaces at L-band for data-cube based soil moisture retrieval; (2) a unified approach for combined active and passive remote sensing of vegetation-covered surfaces with the same input physical parameters; (3) Numerical Maxwell Model in 3D (NMM3D) simulations of a vegetation canopy comprising randomly distributed dielectric cylinders; and (4) a hybrid method based on the generalized T matrix of single objects and Foldy-Lax equations for NMM3D full-wave simulations of the realistic vegetation/forest with vector spherical, spheroidal and cylindrical wave expansoins. The main contributions and novelty of this thesis are NMM3D full-wave simulations of vegetation/forest canopy using the generalized T matrix of the single object and Foldy-Lax equations of multiple scattering among many objects. Before this work, the large-scale full-wave simulations of vegetation/forst such as many tree trunks were deemed very difficult. The NMM3D full-wave simulation results showed that the results of past models significantly overestimate attenuation in a vegetation/forest canopy. The NMM3D full-wave models predict transmissions that are several times greater than that of past models. A much greater microwave transmission means the microwave can better penetrate a vegetation/forest canopy and thus it can be used to retrieve soil moisture. The thesis starts with the DBA to compute the backscattering coefficients for various kinds of vegetation-covered surfaces such as pasture, wheat and canola fields. For the soybean fields, an improved coherent branching model is used. The novel feature of the analytic coherent model consists of conditional probability functions to eliminate the overlapping effects of branches in the former branching models. In order to make use of complex physical models for real time retrieval for satellite missions, the outputs of the physical model are provided as lookup-tables (data-cubes). By inverting the lookup-tables, time-series retrieval of soil moisture is performed. Next, the DBA is extended to calculate the bistatic scattering coefficients. Emissivities are calculated by integrating the bistatic scattering coefficients over the hemispherical solid angle. The backscattering coefficients and emissivities calculated using this approach form a consistent model for combined active and passive microwave remote sensing. In the analytical physical models mentioned above, as well as in another commonly used approach of the radiative transfer equation (RTE), the attenuation of the wave is accounted for by the attenuation rate per unit distance, which originates from the concept of an “effective medium”. Such a model is unsuitable for a vegetation canopy. Because of these issues, NMM3D full-wave simulations of vegetation are pursued. Firstly, the scattering of a vegetation canopy consisting of cylindrical scatterers is calculated. The approach for solving Maxwell’s equations is based on the Foldy-Lax multiple scattering equations (FL) combined with the body of revolution (BOR). For a layer of extended-cylinders distributed in clusters, the NMM3D simulations at C-band show very different results from DBA/RTE. The method FL-BOR is limited for rotationally symmetric objects such as cylinders and circular disks. To perform NMM3D full-wave simulations for realistic vegetation/forests, a hybrid method is used, which is a hybrid of the off-the-shelf techniques and newly developed techniques. The newly developed techniques are the three key steps of the hybrid method: (1) extracting the generalized T matrix of each single object using vector spheroidal/cylindrical waves, (2) vector wave transformations, and (3) solving FL for all the objects.PHDElectrical EngineeringUniversity of Michigan, Horace H. Rackham School of Graduate Studieshttps://deepblue.lib.umich.edu/bitstream/2027.42/153347/1/huanght_1.pd

Review Article: Global Monitoring of Snow Water Equivalent Using High-Frequency Radar Remote Sensing

Seasonal snow cover is the largest single component of the cryosphere in areal extent, covering an average of 46 × 106 km2 of Earth\u27s surface (31 % of the land area) each year, and is thus an important expression and driver of the Earth\u27s climate. In recent years, Northern Hemisphere spring snow cover has been declining at about the same rate (∼ −13 % per decade) as Arctic summer sea ice. More than one-sixth of the world\u27s population relies on seasonal snowpack and glaciers for a water supply that is likely to decrease this century. Snow is also a critical component of Earth\u27s cold regions\u27 ecosystems, in which wildlife, vegetation, and snow are strongly interconnected. Snow water equivalent (SWE) describes the quantity of water stored as snow on the land surface and is of fundamental importance to water, energy, and geochemical cycles. Quality global SWE estimates are lacking. Given the vast seasonal extent combined with the spatially variable nature of snow distribution at regional and local scales, surface observations are not able to provide sufficient SWE information. Satellite observations presently cannot provide SWE information at the spatial and temporal resolutions required to address science and high-socio-economic-value applications such as water resource management and streamflow forecasting. In this paper, we review the potential contribution of X- and Ku-band synthetic aperture radar (SAR) for global monitoring of SWE. SAR can image the surface during both day and night regardless of cloud cover, allowing high-frequency revisit at high spatial resolution as demonstrated by missions such as Sentinel-1. The physical basis for estimating SWE from X- and Ku-band radar measurements at local scales is volume scattering by millimeter-scale snow grains. Inference of global snow properties from SAR requires an interdisciplinary approach based on field observations of snow microstructure, physical snow modeling, electromagnetic theory, and retrieval strategies over a range of scales. New field measurement capabilities have enabled significant advances in understanding snow microstructure such as grain size, density, and layering. We describe radar interactions with snow-covered landscapes, the small but rapidly growing number of field datasets used to evaluate retrieval algorithms, the characterization of snowpack properties using radar measurements, and the refinement of retrieval algorithms via synergy with other microwave remote sensing approaches. This review serves to inform the broader snow research, monitoring, and application communities on progress made in recent decades and sets the stage for a new era in SWE remote sensing from SAR measurements

Multiple Volume Scattering in Random Media and Periodic Structures with Applications in Microwave Remote Sensing and Wave Functional Materials

The objective of my research is two-fold: to study wave scattering phenomena in dense volumetric random media and in periodic wave functional materials. For the first part, the goal is to use the microwave remote sensing technique to monitor water resources and global climate change. Towards this goal, I study the microwave scattering behavior of snow and ice sheet. For snowpack scattering, I have extended the traditional dense media radiative transfer (DMRT) approach to include cyclical corrections that give rise to backscattering enhancements, enabling the theory to model combined active and passive observations of snowpack using the same set of physical parameters. Besides DMRT, a fully coherent approach is also developed by solving Maxwell’s equations directly over the entire snowpack including a bottom half space. This revolutionary new approach produces consistent scattering and emission results, and demonstrates backscattering enhancements and coherent layer effects. The birefringence in anisotropic snow layers is also analyzed by numerically solving Maxwell’s equation directly. The effects of rapid density fluctuations in polar ice sheet emission in the 0.5~2.0 GHz spectrum are examined using both fully coherent and partially coherent layered media emission theories that agree with each other and distinct from incoherent approaches. For the second part, the goal is to develop integral equation based methods to solve wave scattering in periodic structures such as photonic crystals and metamaterials that can be used for broadband simulations. Set upon the concept of modal expansion of the periodic Green’s function, we have developed the method of broadband Green’s function with low wavenumber extraction (BBGFL), where a low wavenumber component is extracted and results a non-singular and fast-converging remaining part with simple wavenumber dependence. We’ve applied the technique to simulate band diagrams and modal solutions of periodic structures, and to construct broadband Green’s functions including periodic scatterers.PHDElectrical EngineeringUniversity of Michigan, Horace H. Rackham School of Graduate Studieshttp://deepblue.lib.umich.edu/bitstream/2027.42/135885/1/srtan_1.pd

Multiple Volume Scattering in Random Media and Periodic Structures with Applications in Microwave Remote Sensing and Wave Functional Materials

The objective of my research is two-fold: to study wave scattering phenomena in dense volumetric random media and in periodic wave functional materials. For the first part, the goal is to use the microwave remote sensing technique to monitor water resources and global climate change. Towards this goal, I study the microwave scattering behavior of snow and ice sheet. For snowpack scattering, I have extended the traditional dense media radiative transfer (DMRT) approach to include cyclical corrections that give rise to backscattering enhancements, enabling the theory to model combined active and passive observations of snowpack using the same set of physical parameters. Besides DMRT, a fully coherent approach is also developed by solving Maxwell’s equations directly over the entire snowpack including a bottom half space. This revolutionary new approach produces consistent scattering and emission results, and demonstrates backscattering enhancements and coherent layer effects. The birefringence in anisotropic snow layers is also analyzed by numerically solving Maxwell’s equation directly. The effects of rapid density fluctuations in polar ice sheet emission in the 0.5~2.0 GHz spectrum are examined using both fully coherent and partially coherent layered media emission theories that agree with each other and distinct from incoherent approaches. For the second part, the goal is to develop integral equation based methods to solve wave scattering in periodic structures such as photonic crystals and metamaterials that can be used for broadband simulations. Set upon the concept of modal expansion of the periodic Green’s function, we have developed the method of broadband Green’s function with low wavenumber extraction (BBGFL), where a low wavenumber component is extracted and results a non-singular and fast-converging remaining part with simple wavenumber dependence. We’ve applied the technique to simulate band diagrams and modal solutions of periodic structures, and to construct broadband Green’s functions including periodic scatterers.PHDElectrical EngineeringUniversity of Michigan, Horace H. Rackham School of Graduate Studieshttps://deepblue.lib.umich.edu/bitstream/2027.42/137141/1/srtan_1.pd

Three-Dimensional Electromagnetic Scattering from Layered Media with Rough Interfaces for Subsurface Radar Remote Sensing

The objective of this dissertation is to develop forward scattering models for active microwave remote sensing of natural features represented by layered media with rough interfaces. In particular, soil profiles are considered, for which a model of electromagnetic scattering from multilayer rough surfaces with/without buried random media is constructed. Starting from a single rough surface, radar scattering is modeled using the stabilized extended boundary condition method (SEBCM). This method solves the long-standing instability issue of the classical EBCM, and gives three-dimensional full wave solutions over large ranges of surface roughnesses with higher computational e±ciency than pure numerical solutions, e.g., method of moments (MoM). Based on this single surface solution, multilayer rough surface scattering is modeled using the scattering matrix approach and the model is used for a comprehensive sensitivity analysis of the total ground scattering as a function of layer separation, subsurface statistics, and sublayer dielectric properties. The buried inhomogeneities such as rocks and vegetation roots are considered for the first time in the forward scattering model. Radar scattering from buried random media is modeled by the aggregate transition matrix using either the recursive transition matrix approach for spherical or short-length cylindrical scatterers, or the generalized iterative extended boundary condition method we developed for long cylinders or root-like cylindrical clusters. These approaches take the field interactions among scatterers into account with high computational efficiency. The aggregate transition matrix is transformed to a scattering matrix for the full solution to the layered-medium problem. This step is based on the near-to-far field transformation of the numerical plane wave expansion of the spherical harmonics and the multipole expansion of plane waves. This transformation consolidates volume scattering from the buried random medium with the scattering from layered structure in general. Combined with scattering from multilayer rough surfaces, scattering contributions from subsurfaces and vegetation roots can be then simulated. Solutions of both the rough surface scattering and random media scattering are validated numerically, experimentally, or both. The experimental validations have been carried out using a laboratory-based transmit-receive system for scattering from random media and a new bistatic tower-mounted radar system for field-based surface scattering measurements.Ph.D.Electrical EngineeringUniversity of Michigan, Horace H. Rackham School of Graduate Studieshttp://deepblue.lib.umich.edu/bitstream/2027.42/91459/1/xduan_1.pd

Passive Microwave Remote Sensing of Snow Layers Using Novel Wideband Radiometer Systems and RFI Mitigation

Climate change can reduce the availability of water resources in many regions, and it will affect agriculture, industry, and energy supply. Snowpack monitoring is important in water resource management as well as flood and avalanche protection. The rapid melting process due to global warming changes the snowpacks' annual statistics, including the extent, and the snow water equivalent (SWE) of seasonal snowpacks, which results in non-stationary annual statistics that should be monitored in nearly daily intervals. The development of advanced radiometric sensors capable of accurately measuring the snowpack thickness and SWE is needed for the long-term study of the snowpack parameters' statistical changes. Passive microwave radiometry provides a means for measuring the microwave emission from a scene of snow and ice. A Wideband Autocorrelation Radiometer (ac{WiBAR}) operating from 1-2~GHz measures spontaneous emission from snowpack at long wavelengths where the scattering is minimized, but the snow layer coherent effects are preserved. By using a wide bandwidth to measure the spacing between frequencies of constructive and destructive interference of the emission from the soil under the snow, it can reveal the microwave travel time through the snow, and thus the snow depth. However, narrowband radio frequency interference (RFI) in the WiBAR's frequency of operations reduces the ability of the WiBAR to measure the thickness accurately. In addition, the current WiBAR system is a frequency domain, FD-WiBAR, system that uses a field-portable spectrum analyzer to collect the data and suffers from high data acquisition time which limits its applications for spaceborne and airborne technologies. In this work, a novel frequency tunable microwave comb filter is proposed for RFI mitigation. The frequency response of the proposed filter has a pattern with many frequencies band-pass and band rejection that preserves the frequency span while reducing the RFI. Moreover, we demonstrate time-domain WiBAR, TD-WiBAR, which presented as an alternative method for FD-WiBAR, and is capable of providing faster data acquisition. A new time-domain calibration is also developed for TD-WiBAR and evaluated with the frequency domain calibration. To validate the TD-WiBAR method, simulated laboratory measurements are performed using a microwave scene simulator circuit. Then the WiBAR instrument is enhanced with the proposed comb filter and showed the RFI mitigation in time-domain mode on an instrument bench test. Furthermore, we analyze the effects of an above snow vegetation layer on brightness temperature spectra, particularly the possible decay of wave coherence arising from volume scattering in the vegetation canopy. In our analysis, the snow layer is assumed to be flat, and its upward emission and surface reflectivities are modeled by a fully coherent model, while an incoherent radiative transfer model describes the volume scattering from the vegetation layer. We proposed a unified framework of vegetation scattering using radiative transfer (RT) theory for passive and active remote sensing of vegetated land surfaces, especially those associated with moderate-to-large vegetation water contents (VWCs), e.g., forest field. The framework allows for modeling passive and active microwave signatures of the vegetated field with the same physical parameters describing the vegetation structure. The proposed model is validated with the passive and active L-band sensor (PALS) acquired in SMAPVEX12 measurements in 2012, demonstrating the applicability of this model.PHDElectrical and Computer EngineeringUniversity of Michigan, Horace H. Rackham School of Graduate Studieshttp://deepblue.lib.umich.edu/bitstream/2027.42/169653/1/maryamsa_1.pd

Exploring bistatic scattering modeling for land surface applications using radio spectrum recycling in the Signal of Opportunity Coherent Bistatic Simulator

The potential for high spatio-temporal resolution microwave measurements has urged the adoption of the signals of opportunity (SoOp) passive radar technique for use in remote sensing. Recent trends in particular target highly complex remote sensing problems such as root-zone soil moisture and snow water equivalent. This dissertation explores the continued open-sourcing of the SoOp coherent bistatic scattering model (SCoBi) and its use in soil moisture sensing applications. Starting from ground-based applications, the feasibility of root-zone soil moisture remote sensing is assessed using available SoOp resources below L-band. A modularized, spaceborne model is then developed to simulate land-surface scattering and delay-Doppler maps over the available spectrum of SoOp resources. The simulation tools are intended to provide insights for future spaceborne modeling pursuits

Modélisation de l’émission micro-onde hivernale en forêt boréale canadienne

La caractérisation du couvert nival en forêt boréale est un élément important pour la compréhension des régimes climatiques et hydrologiques. Depuis plusieurs années, l’utilisation des micro-ondes passives est étudiée pour l’estimation de l’équivalent en eau de la neige (SWE : Snow Water Equivalent) à partir de capteurs satellitaires. Les algorithmes empiriques traditionnels étant limités en forêt boréale, le couplage d’un modèle de transfert radiatif (MTR) micro-onde passive (qui prend en compte les contributions du sol, de la neige, de la végétation et de l’atmosphère) avec un modèle de neige pour l’inversion du SWE semble une avenue prometteuse. La thèse vise donc à coupler un MTR avec le schéma de surface du modèle climatique canadien (CLASS) dans une perspective d’application opérationnelle pour les estimations de SWE à partir de données satellitaires micro-onde à 10.7, 19 et 37 GHz. Dans ce contexte, certains aspects centraux du MTR, dont l’effet de la taille des grains ainsi que la contribution de la végétation sont développés et quantifiés. Le premier aspect étudié dans la thèse concerne l’adaptation du modèle d’émission micro-onde passive DMRT-ML (Dense media radiative transfer theory – multi layer) pour l’intégration d’une nouvelle métrique représentant la taille des grains (surface spécifique des grains de neige: SSA). L’étude basée sur des mesures radiométriques et de neige in situ, montre la pertinence de l’utilisation de la SSA dans DMRT-ML et permet d’analyser le sens physique de l’adaptation nécessaire pour amener le modèle à simuler les températures de brillance (T[indice inférieur B) de la neige avec une erreur quadratique moyenne minimale de l’ordre de 13 K. Dans un contexte du couplage entre le modèle de neige de CLASS et DMRT-ML, un modèle d’évolution de la SSA est ensuite implémenté dans CLASS. Les SSA simulées par le module développé sont validées avec des données in situ basées sur la réflectance de la neige dans l’infrarouge à courte longueur d’onde pour différents types d’environnement. Au niveau de la contribution de la végétation, le modèle γ-ω a été étudié à partir de différentes bases de données (satellite, avion et au sol) en forêt boréale dense. L’étude montre l’importance de la considération de la diffusion (ω) pour l’estimation de l’émission de la végétation, paramètre auparavant généralement négligé aux hautes fréquences. Ensuite, des relations entre les transmissivités et certains paramètres structuraux de la forêt, dont l’indice de surface foliaire (LAI), ont été établies pour des forêts boréales en été. Des valeurs d’albédo de diffusion (ω) ainsi que les paramètres définissant la réflectivité du sol (QH) en forêt boréale ont aussi été inversées. Finalement, les simulations de T [indice inférieur] B issues du couplage du MTR (DMRT-ML, modèle γ-ω, et modèle atmosphérique) avec CLASS (dont les SSA simulées) ont été comparées avec les données AMSR-E sur une série temporelle continue de sept ans. Les premières comparaisons montrent une différence entre les paramètres de végétation (γ-ω) d’été et d’hiver, ainsi qu’une importante contribution des croûtes de glace dans la neige au signal. Les simulations du modèle ajusté montrent une bonne correspondance avec les observations d’AMSR-E (de l’ordre de 3 à 7 K selon la fréquence et la polarisation). Des tests de sensibilité montrent par contre une faible sensibilité du MTR/CLASS au SWE pour des forêts denses et des couverts nivaux épais. Le MTR-CLASS développé pourrait permettre l’assimilation de températures de brillance satellitaires en forêt boréale dans des systèmes opérationnels pour l’amélioration de paramètres de surface, dont la neige, dans les modèles météorologiques et climatiques