Information retrieval from spaceborne GNSS Reflectometry observations using physics- and learning-based techniques

This dissertation proposes a learning-based, physics-aware soil moisture (SM) retrieval algorithm for NASA’s Cyclone Global Navigation Satellite System (CYGNSS) mission. The proposed methodology has been built upon the literature review, analyses, and findings from a number of published studies throughout the dissertation research. Namely, a Sig- nals of Opportunity Coherent Bistatic scattering model (SCoBi) has been first developed at MSU and then its simulator has been open-sourced. Simulated GNSS-Reflectometry (GNSS-R) analyses have been conducted by using SCoBi. Significant findings have been noted such that (1) Although the dominance of either the coherent reflections or incoher- ent scattering over land is a debate, we demonstrated that coherent reflections are stronger for flat and smooth surfaces covered by low-to-moderate vegetation canopy; (2) The influ- ence of several land geophysical parameters such as SM, vegetation water content (VWC), and surface roughness on the bistatic reflectivity was quantified, the dynamic ranges of reflectivity changes due to SM and VWC are much higher than the changes due to the surface roughness. Such findings of these analyses, combined with a comprehensive lit- erature survey, have led to the present inversion algorithm: Physics- and learning-based retrieval of soil moisture information from space-borne GNSS-R measurements that are taken by NASA’s CYGNSS mission. The study is the first work that proposes a machine learning-based, non-parametric, and non-linear regression algorithm for CYGNSS-based soil moisture estimation. The results over point-scale soil moisture observations demon- strate promising performance for applicability to large scales. Potential future work will be extension of the methodology to global scales by training the model with larger and diverse data sets

GNSS Reflectometry for land surface monitoring and buried object detection

Global Navigation Satellite System Reflectometry (GNSS-R) is attracting growing interest nowadays for several remote sensing applications. As a bistatic radar, the transmitter and the receiver are not co-located and in the special case of GNSS-R, the GNSS satellites are acting as transmitters and the receiver can be mounted either in a static position or onboard a aircraft or low orbit satellite. Various information about the surface from where the GNSS signals are reflected or scattered can be extracted by means of reflected signal strength, code delay, carrier phase delay, interference with direct GNSS signals and so on. Possible applications cover soil moisture retrieval, ice topography and thickness detection, snow depth estimation, vegetation coverage, sea state monitoring such as sea wind and surface roughness, sea salinity… In this work, soil moisture retrieval was mostly focused on. Hardware including antennas and receivers was studied and designed. Our first strategy of soil moisture retrieval is to apply a single Left Hand Circular Polarization (LHCP) antenna for reflected signal reception. Therefore multiple types of antennas such as the helix antenna, the patch antenna and several commercial antennas were designed, simulated or tested in the anechoic chamber. Two receiver solutions were used in our group and both of them apply the SiGe GPS frontend. The first solution is a PC based one: the collection and store of the raw incoming reflected GPS signals were done by the NGrab software (designed by NAVSAS Group of Politecnico di Torino) installed in a standard PC. The other solution was developed in our group and it is operated by a single Hackberry board, which consists of power supply, storage subsystem and customized Linux Debian operating system. The light weight and small size enable this compact receiver to perform flight measurement onboard UAVs. Both of the above mentioned receivers only store raw sampled data and no real time signal processing is performed on board. Post processing is done by Matlab program which makes correlations in both time and frequency domain with incoming signals using the local generated GPS C/A code replica. The so-called Delay Doppler Map (DDM) is therefore generated through this correlation. Signal to Noise Ratio (SNR) can be calculated through Delay Waveform (DW) which is extracted from DDM at the Doppler frequency where the correlation peak exists. Received signal power can be obtained knowing the noise power which is given in a standard equation. In order to better plan a static measurement and to georeference specular points on the surface, programs for georeferencing specular points on either Google Maps or an x-y plane centered at the receiver position were developed. Fly dynamics in terms of roll, pitch and yaw influencing the antenna gain due to the variation of incident angles were also studied in order to compensate the gain to the received signal. Two soil moisture retrieval algorithms were derived corresponding to two receiving schemes. The first one is for the receiving of only LHCP reflected signals. In this case, the surface is assumed to be perfectly smooth and the received signal is seen to consist of only coherent component caused by specular reflection. Dielectric constant can be retrieved from the processed SNR. Two measurement campaigns were carried out using this single LHCP system. The first campaign is a flight measurement overflown a big portion of rice fields when most of the fields were flooded. It was a test measurement on the SNR sensitivity to water/no-water surfaces and an attempt of dielectric constant retrieval was also performed. SNR showed good sensitivity to the surface water content and dielectric constant was also checked to be reasonable. The second campaign is in static positions and it includes two experiments. This campaign initially aimed at testing the sensitivity of the compact receiver to different surface moisture. Results of both SNR and retrieved dielectric constant showed to be coherent with the surface moisture changes. The other retrieval algorithm is for the receiving of both LHCP and RHCP reflected signals concurrently. The cross polarization power ratio (LHCP/RHCP) is believed to be independent of surface roughness by several previous studies and this idea was also verified during the deriving process for either specular reflection case (only coherent component) or diffuse scattering condition (incoherent component). For diffuse scattering, three well known models were applied which are the Kirchhoff Approximation in stationary-phase approximation (Kirchhoff Geometrical Optics, KGO), Kirchhoff Approximation in Physical Optics Approximation (KPO) and Small Perturbation Method (SPM). These three models cover different roughness surfaces from very rough (KGO) to slightly rough surfaces (SPM). All the derived results of cross polarization ratio for the three models were verified to be independent of surface properties and depend on only dielectric constant of soil and incident angle. A new application of GNSS-R technique for the possibility of detection of buried objects was firstly investigated by our group. It has the potential use for man-made mines detection in the military field. Two measurement campaigns were carried out and the variation of the SNR level due to the presence of a metallic object was investigated. The first measurement campaign was performed in a static condition on a sandy terrain to check the functionality of the system. And the presence of the metallic object was detected also in the case of wet terrain. In the second measurement campaign, the antenna was moving along a given path and the possibility of detecting the object dimensions was highlighted. The results show the possibility of adopting this technique on board a remotely controlled UAV for metal object and even its dimension detection. A measurement of snow depth attempting to relate it to reflected LHCP SNR is briefly presented and discussed in Chapter 7

Temporal Characteristics of Boreal Forest Radar Measurements

Radar observations of forests are sensitive to seasonal changes, meteorological variables and variations in soil and tree water content. These phenomena cause temporal variations in radar measurements, limiting the accuracy of tree height and biomass estimates using radar data. The temporal characteristics of radar measurements of forests, especially boreal forests, are not well understood. To fill this knowledge gap, a tower-based radar experiment was established for studying temporal variations in radar measurements of a boreal forest site in southern Sweden. The work in this thesis involves the design and implementation of the experiment and the analysis of data acquired. The instrument allowed radar signatures from the forest to be monitored over timescales ranging from less than a second to years. A purpose-built, 50 m high tower was equipped with 30 antennas for tomographic imaging at microwave frequencies of P-band (420-450 MHz), L-band (1240-1375 MHz) and C-band (5250-5570 MHz) for multiple polarisation combinations. Parallel measurements using a 20-port vector network analyser resulted in significantly shorter measurement times and better tomographic image quality than previous tower-based radars. A new method was developed for suppressing mutual antenna coupling without affecting the range resolution. Algorithms were developed for compensating for phase errors using an array radar and for correcting for pixel-variant impulse responses in tomographic images. Time series results showed large freeze/thaw backscatter variations due to freezing moisture in trees. P-band canopy backscatter variations of up to 10 dB occurred near instantaneously as the air temperature crossed 0⁰C, with ground backscatter responding over longer timescales. During nonfrozen conditions, the canopy backscatter was very stable with time. Evidence of backscatter variations due to tree water content were observed during hot summer periods only. A high vapour pressure deficit and strong winds increased the rate of transpiration fast enough to reduce the tree water content, which was visible as 0.5-2 dB backscatter drops during the day. Ground backscatter for cross-polarised observations increased during strong winds due to bending tree stems. Significant temporal decorrelation was only seen at P-band during freezing, thawing and strong winds. Suitable conditions for repeat-pass L-band interferometry were only seen during the summer. C-band temporal coherence was high over timescales of seconds and occasionally for several hours for night-time observations during the summer. Decorrelation coinciding with high transpiration rates was observed at L- and C-band, suggesting sensitivity to tree water dynamics.The observations from this experiment are important for understanding, modelling and mitigating temporal variations in radar observables in forest parameter estimation algorithms. The results also are also useful in the design of spaceborne synthetic aperture radar missions with interferometric and tomographic capabilities. The results motivate the implementation of single-pass interferometric synthetic aperture radars for forest applications at P-, L- and C-band

GNSS reflectometry for land remote sensing applications

Soil moisture and vegetation biomass are two essential parameters from a scienti c and economical point of view. On one hand, they are key for the understanding of the hydrological and carbon cycle. On the other hand, soil moisture is essential for agricultural applications and water management, and vegetation biomass is crucial for regional development programs. Several remote sensing techniques have been used to measure these two parameters. However, retrieving soil moisture and vegetation biomass with the required accuracy, and the appropriate spatial and temporal resolutions still remains a major challenge. The use of Global Navigation Satellite Systems (GNSS) reflected signals as sources of opportunity for measuring soil moisture and vegetation biomass is assessed in this PhD Thesis. This technique, commonly known as GNSS-Reflectometry (GNSS-R), has gained increasing interest among the scienti c community during the last two decades due to its unique characteristics. Previous experimental works have already shown the capabilities of GNSS-R to sense small reflectivity changes on the surface. The use of the co- and cross-polarized reflected signals was also proposed to mitigate nuisance parameters, such as soil surface roughness, in the determination of soil moisture. However, experimental evidence of the suitability of that technique could not be demonstrated. This work analyses from a theoretical and an experimental point of view the capabilities of polarimetric observations of GNSS reflected signals for monitoring soil moisture and vegetation biomass. The Thesis is structured in four main parts. The fi rst part examines the fundamental aspects of the technique and provides a detailed review of the GNSS-R state of the art for soil moisture and vegetation monitoring. The second part deals with the scattering models from land surfaces. A comprehensive description of the formation of scattered signals from rough surfaces is provided. Simulations with current state of the art models for bare and vegetated soils were performed in order to analyze the scattering components of GNSS reflected signals. A simpli ed scattering model was also developed in order to relate in a straightforward way experimental measurements to soil bio-geophysical parameters. The third part reviews the experimental work performed within this research. The development of a GNSS-R instrument for land applications is described, together with the three experimental campaigns carried out in the frame of this PhD Thesis. The analysis of the GNSS-R and ground truth data is also discussed within this part. As predicted by models, it was observed that GNSS scattered signals from natural surfaces are a combination of a coherent and an incoherent scattering components. A data analysis technique was proposed to separate both scattering contributions. The use of polarimetric observations for the determination of soil moisture was demonstrated to be useful under most soil conditions. It was also observed that forests with high levels of biomass could be observed with GNSS reflected signals. The fourth and last part of the Thesis provides an analysis of the technology perspectives. A GNSS-R End-to-End simulator was used to determine the capabilities of the technique to observe di erent soil reflectivity conditions from a low Earth orbiting satellite. It was determined that high accuracy in the estimation of reflectivity could be achieved within reasonable on-ground resolution, as the coherent scattering component is expected to be the predominant one in a spaceborne scenario. The results obtained in this PhD Thesis show the promising potential of GNSS-R measurements for land remote sensing applications, which could represent an excellent complementary observation for a wide range of Earth Observation missions such as SMOS, SMAP, and the recently approved ESA Earth Explorer Mission Biomass.La humedad del suelo y la biomasa de la vegetaci on son dos parametros clave desde un punto de vista tanto cient co como econ omico. Por una parte son esenciales para el estudio del ciclo del agua y del carbono. Por otra parte, la humedad del suelo es esencial para la gesti on de las cosechas y los recursos h dricos, mientras que la biomasa es un par ametro fundamental para ciertos programas de desarrollo. Varias formas de teledetección se han utilizado para la observaci on remota de estos par ametros, sin embargo, su monitorizaci on con la precisi on y resoluci on necesarias es todav a un importante reto tecnol ogico. Esta Tesis evalua la capacidad de medir humedad del suelo y biomasa de la vegetaci on con señales de Sistemas Satelitales de Posicionamiento Global (GNSS, en sus siglas en ingl es) reflejadas sobre la Tierra. La t ecnica se conoce como Reflectometr í a GNSS (GNSS-R), la cual ha ganado un creciente inter es dentro de la comunidad científ ca durante las dos ultimas d ecadas. Experimentos previos a este trabajo ya demostraron la capacidad de observar cambios en la reflectividad del terreno con GNSS-R. El uso de la componente copolar y contrapolar de la señal reflejada fue propuesto para independizar la medida de humedad del suelo de otros par ametros como la rugosidad del terreno. Sin embargo, no se pudo demostrar una evidencia experimental de la viabilidad de la t ecnica. En este trabajo se analiza desde un punto de vista te orico y experimental el uso de la informaci on polarim etrica de la señales GNSS reflejadas sobre el suelo para la determinaci on de humedad y biomasa de la vegetaci on. La Tesis se estructura en cuatro partes principales. En la primera parte se eval uan los aspectos fundamentales de la t ecnica y se da una revisi on detallada del estado del arte para la observaci on de humedad y vegetaci on. En la segunda parte se discuten los modelos de dispersi on electromagn etica sobre el suelo. Simulaciones con estos modelos fueron realizadas para analizar las componentes coherente e incoherente de la dispersi on de la señal reflejada sobre distintos tipos de terreno. Durante este trabajo se desarroll o un modelo de reflexi on simpli cado para poder relacionar de forma directa las observaciones con los par ametros geof sicos del suelo. La tercera parte describe las campañas experimentales realizadas durante este trabajo y discute el an alisis y la comparaci on de los datos GNSS-R con las mediciones in-situ. Como se predice por los modelos, se comprob o experimentalmente que la señal reflejada est a formada por una componente coherente y otra incoherente. Una t ecnica de an alisis de datos se propuso para la separacióon de estas dos contribuciones. Con los datos de las campañas experimentales se demonstr o el bene cio del uso de la informaci on polarim etrica en las señales GNSS reflejadas para la medici on de humedad del suelo, para la mayor a de las condiciones de rugosidad observadas. Tambi en se demostr o la capacidad de este tipo de observaciones para medir zonas boscosas densamente pobladas. La cuarta parte de la tesis analiza la capacidad de la t ecnica para observar cambios en la reflectividad del suelo desde un sat elite en orbita baja. Los resultados obtenidos muestran que la reflectividad del terreno podr a medirse con gran precisi on ya que la componente coherente del scattering ser a la predominante en ese tipo de escenarios. En este trabajo de doctorado se muestran la potencialidades de la t ecnica GNSS-R para observar remotamente par ametros del suelo tan importantes como la humedad del suelo y la biomasa de la vegetaci on. Este tipo de medidas pueden complementar un amplio rango de misiones de observaci on de la Tierra como SMOS, SMAP, y Biomass, esta ultima recientemente aprobada para la siguiente misi on Earth Explorer de la ESA

Bistatic Radar Configuration for Soil Moisture Retrieval: Analysis of the Spatial Coverage

Some outcomes of a feasibility analysis of a spaceborne bistatic radar mission for soil moisture retrieval are presented in this paper. The study starts from the orbital design of the configuration suitable for soil moisture estimation identified in a previous study. This configuration is refined according to the results of an analysis of the spatial resolution. The paper focuses on the assessment of the spatial coverage i.e., on the verification that an adequate overlap between the footprints of the antennas is ensured and on the duty cycle, that is the fraction of orbital period during which the bistatic data are acquired. A non-cooperating system is considered, in which the transmitter is the C-band Advanced Synthetic Aperture Radar aboard Envisat. The best performances in terms of duty cycle are achieved if the transmitter operates in Wide Swath Mode. The higher resolution Image Swath Modes that comply with the selected configuration have a duty cycle that is never less than 12% and can exceed 21%. When Envisat operates in Wide Swath Mode, the bistatic system covers a wide latitude range across the equator, while in some of the Image Swath Modes, the bistatic measurements, collected from the same orbit, cover mid-latitude areas. In the latter case, it might be possible to achieve full coverage in an Envisat orbit repeat cycle, while, for a very large latitude range such as that covered in Wide Swath Mode, bistatic acquisitions could be obtained over about 65% of the area