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

GNSS-IR Model of Sea Level Height Estimation Combining Variational Mode Decomposition

The Global Navigation Satellite System-Reflections (GNSS-R) signal has been confirmed to be used to retrieve sea level height. At present, the GNSS-Interferometric Reflectometry (GNSS-IR) technology based on the least square method to process signal-to-noise ratio (SNR) data is restricted by the satellite elevation angle in terms of accuracy and stability. This paper proposes a new GNSS-IR model combining variational mode decomposition (VMD) for sea level height estimation. VMD is used to decompose the SNR data into intrinsic mode functions (IMF) of layers with different frequencies, remove the IMF representing the trend item of the SNR data, and reconstruct the remaining IMF components to obtain the SNR oscillation item. In order to verify the validity of the new GNSS-IR model, the measurement data provided by the Onsala Space Observatory in Sweden is used to evaluate the performance of the algorithm and its stability in high elevation range. The experimental results show that the VMD method has good results in terms of accuracy and stability, and has advantages compared to other methods. For the half-year GNSS SNR data, the root mean square error (RMSE) and correlation coefficient of the new model based on the VMD method are 4.86 cm and 0.97, respectively

Advanced GNSS-R instruments for altimetric and scatterometric applications

This work is the result of more than eight years during a bachelor thesis, a master thesis, and the Ph.D. thesis dedicated to the development of the Microwave Interferometric Reflectometer (MIR) instrument. It summarizes all the knowledge acquired during this time, and describes the MIR instrument as detailed as possible. MIR is a Global Navigation Satellite System - Reflectometer (GNSS-R), that is, an instrument that uses Global Navigation Satellite System (GNSS) signals scattered on the Earth's surface to retrieve geophysical parameters. These signals are received below the noise level, but since they have been spread in the frequency domain using spread-spectrum techniques, and in particular using the so-called Pseudo Random Noise (PRN) codes, it is still possible to retrieve them because of the large correlation gain achieved. In GNSS-R, two main techniques are used for this purpose: the conventional technique cGNSS-R and the interferometric one iGNSS-R, each with its pros and cons. In the former technique, the reflected signal is cross-correlated against a locally generated clean-replica of the transmitted signal. In the latter technique the reflected signal is cross-correlated with the direct one. Nowadays multiple GNSS systems coexist, transmitting narrow and wide, open and private signals. A comparison between systems, signals, and techniques in fair conditions is necessary. The MIR instrument has been designed as an airborne instrument for that purpose: the instrument has two arrays, an up-looking one, and a down-looking one, each with 19 dual-band antennas in a hexagonal distribution. The instrument is able to form 2 beams at each frequency band (L1/E1, and L5/E5A), which are pointing continuously to the desired satellites taking into account their position, as well as the instrument's position and attitude. The data is sampled and stored for later post-processing. Last but not least, MIR is auto-calibrated using similar signals to the ones transmitted by the GNSS satellites. During the instrument development, the Distance Measurement Equipment/TACtical Air Navigation (DME/TACAN) signals from the Barcelona airport threatened to disrupt the interferometric technique. These signals were also studied, and it was concluded that the use of a mitigation systems were as strongly recommended. The interferometric technique was also affected by the unwanted contribution of other satellites. The impact of these contributions was studied using real data gathered during this Ph.D. thesis. During these 8 years, the instrument was designed, built, tested, and calibrated. A field campaign was carried out in Australia between May 2018 and June 2018 to determine the instrument's accuracy in sensing soil moisture and sea altimetry. This work describes each of these steps in detail and aims to be helpful for those who decide to continue the legacy of this instrument.Este trabajo es el resultado de más de 8 años de doctorado dedicados al desarrollo del instrumento Microwave Interferometric Reflectometer (MIR). Esta tesis resume todo el conocimiento adquirido durante este tiempo, y describe el MIR lo más detalladamente posible. El MIR es un Reflectómetro de señales de Sistemas Globales de Navegación por Satélite (GNSS-R), es decir, es un instrumento que usa señales de GNSS reflejadas en la superficie de la tierra para obtener parámetros geofísicos. Estas señales son recibidas bajo el nivel de ruido, pero dado que han sido ensanchadas en el dominio frecuencial usando técnicas de espectro ensanchado, y en particular usando códigos Pseudo Random Noise (PRN), es todavía posible recibirlas debido a la elevada ganancia de correlación. En GNSS-R existen dos técnicas para este propósito: la convencional (cGNSS-R), y la interferométrica (iGNSS-R), cada una con sus pros y sus contras. En la primera se calcula la correlación cruzada de la señal reflejada y de una réplica generada del código transmitido. En la segunda técnica se calcula la correlación cruzada de la señal reflejada y de la señal directa. Hoy en día muchos sistemas GNSS coexisten, transmitiendo señales de distintos anchos de banda, algunas públicas y otras privadas. Una comparación entre sistemas, señales, y técnicas en condiciones justas es necesaria. El MIR es un instrumento aerotransportado diseñado como para ese propósito: el instrumento tiene dos arrays de antenas, uno apuntando al cielo, y otro apuntando al suelo, cada uno con 19 antenas doble banda en una distribución hexagonal. El instrumento puede formar 2 haces en cada banda frecuencial (L1/E1 y L5/E5A) que apuntan continuamente a los satélites deseados teniendo en cuenta su posición, y la posición y actitud del instrumento. Los datos son guardados para ser procesados posteriormente. Por último pero no menos importante, el MIR se calibra usando señales similares a las transmitidas por los satélites de GNSS. Durante el desarrollo del instrumento, señales del sistema Distance Measuremt Equi Distance Measurement Equipment/TACtical Air Navigation (DME/TACAN) del aeropuerto de Barcelona mostraron ser una amenaza para la técnica interferométrica. Estas señales fueron estudiadas y se concluyó que era encarecidamente recomendado el uso de sistemas de mitigación de interferencias. La técnica interferométrica también se ve afectada por las contribuciones no deseadas de otros satélites, llamado cross-talk. El impacto del cross-talk fue estudiado usando datos reales tomados durante esta tesis doctoral. A lo largo de estos 8 años el instrumento ha sido diseñado, construido, testeado y calibrado. Una campaña de medidas fue llevada a cabo en Australia entre Mayo de 2018 y Junio de 2018 para determinar la capacidad del instrumento para estimar la humedad del terreno y la altura del mar. Este documento describe cada uno de estos pasos al detalle y espera resultar útil para aquellos que decidan continuar con el legado de este instrumento.Postprint (published version

Contributions to land, sea, and sea ice remote sensing using GNSS-reflectometry

This PhD thesis researches the use of passive remote sensing techniques using signals transmitted from the navigation satellites (GNSS) in order to retrieve different geophysical parameters. The thesis consists of two different parts without taking into account the introduction, the state of the art and the conclusions. The first part analyzes the Interference Pattern Technique (IPT), which was previously used in another PhD thesis, and proposes some key improvements. First, the addition of horizontal polarization to the existing vertical polarization is proposed. Then, the retrieval of soil moisture is studied using the horizontal polarization only and combining both polarizations to correct for the surface roughness effects. It is also demonstrated that the phase difference between the two interference patterns is directly related to soil moisture content. A field campaign was conducted in Australia to test empirically all the theoretical developments and algorithms. Secondly, the possibility of measuring Significant Wave Height (SWH) and Mean Sea Surface Level (MSSL) using the IPT is studied. A three month field campaign over coastal sea is devoted to that study. The SWH retrieval is a new estimation algorithm based on measuring the point where the interference pattern loses coherence. The MSSL retrieval is based on the estimation of the IPT oscillation frequency, testing different spectral estimators to improve the accuracy. Since the IPT is limited in coverage due to its static requirements, the research conducted in this thesis migrated to scatterometric GNSS-R techniques. The main goal that migration was to increase coverage of the different GNSS-R instruments. Therefore, the second part of this thesis analyzes the applicability of a scatterometric technique from different platforms: ground-based (mobile and fixed), airborne, and spaceborne. The ground-based still platforms have allowed to develop a soil moisture retrieval algorithm. The ground-based moving platforms have extended the validity of that algorithm. Airborne platforms have been used to study the reflected electric field statistics when the surface reflecting surface is varying (smooth or rough land, and sea). They have also been used to develop different algorithms to measure the coherent and incoherent scattered components depending on the data structure (real-data or complex data). Coherent reflectivity measured from airborne platforms has been compared to other techniques such microwave radiometry, which is highly used in the soil moisture retrieval from spaceborne sensors, and other sensors using optical, multispectral and thermal frequency bands. These relationships between microwave radiometry and GNSS-R measurements suggests the potential synergy of both techniques. A sea ice detection algorithm is also developed using scatterometric GNSS-R data from the UK TDS-1 mission. This algorithm is based on measuring the degree of coherence of the reflected waveform. Finally, a field campaign was conducted to study the effect of vegetation on the GNSS signals that pass through it in order to take into account and correct the effect of vegetation in the GNSS-R data and in the soil moisture retrieval algorithms.Aquesta tesi doctoral aprofundeix en el coneixement de les tècniques de teledetecció passives utilitzant senyals emesos pels satèl·lits de navegació (GNSS) amb l'objectiu de recuperar diferents paràmetres geofísics del terreny. La tesi conté dues parts ben diferenciades a banda de la introducció, estat de l'art i conclusions. La primera part analitza la tècnica coneguda com a patró d'interferències, utilitzada prèviament en una altra tesi doctoral, i proposa certes millores per la seva aplicabilitat. En primer lloc es decideix afegir polarització horitzontal a la ja existent polarització vertical, i s'estudia la recuperació d'humitat del sòl utilitzant només polarització horitzontal i combinant les dues polaritzacions per corregir els efectes de la rugositat del terreny. A continuació es demostra que la mesura de desfasament entre els dos patrons d'interferència està directament relacionada amb la humitat del terreny. Es va realitzar una campanya de mesures a Austràlia per provar empíricament tots els desenvolupaments teòrics i algorismes proposats. En segon lloc s'analitza l'aplicabilitat del patró d'interferències en la mesura de l'altura de les onades (SWH) i del nivell del mar (MSSL), tots dos de forma precisa. L'estimació de l'alçada de les onades és un procés totalment nou basat en mesurar el punt on el patró d'interferències perd la coherència. L'estimació del nivell del mar es basa en l'anàlisi espectral del patró d'interferències provant diferents estimadors espectrals. Atès que la tècnica del patró d'interferència està limitada en cobertura per les seves característiques estàtiques, la investigació duta a terme en aquesta tesi doctoral va migrar cap a tècniques GNSS-R escateromètriques. El principal objectiu a assolir va ser el d'augmentar la cobertura dels diferents instruments GNSS-R de mesura. En conseqüència, la segona part d'aquesta tesi analitza l'aplicabilitat d'aquestes tècniques des de diferents plataformes terrestres (mòbils i fixes), aerotransportades i satèl·lit. Les plataformes terrestres fixes han permès derivar algoritmes de recuperació d'humitat i les mòbils estendre la validació d'aquests. Les plataformes aerotransportades s'han utilitzat per mirar l'estadística del camp elèctric reflectit quan la superfície on es reflecteixen els senyals GNSS va variant (terra plana o terra rugosa, i mar). També han servit per desenvolupar diferents algorismes amb l'objectiu de determinar les components coherent i incoherent del senyal reflectit. De la mateixa manera, dades de reflectivitat coherent mesurades des d'aquestes plataformes han estat comparades amb altres tècniques de teledetecció passiva com la radiometria de microones, altament utilitzada en la mesura d'humitat de terreny, i altres sensors òptics, multi-espectrals, i tèrmics. Aquests resultats han permès suggerir la possible sinergia de dades d'ambdues tecnologies. Un algorisme per detectar la presència de gel sobre el mar també ha estat desenvolupat mitjançant l'ús de dades GNSS-R escateromètriques satel·litals de la missió UK TDS-1. Aquest algorisme es basa en mesurar el grau de coherència de la forma d'ona reflectida. Finalment, s'ha realitzat un estudi de l'efecte de la vegetació en els senyals GNSS que la travessen, per tal de poder corregir aquest efecte en els algoritmes de recuperació d'humitat del terreny

Flaring and pollution detection in the Niger Delta using Remote Sensing

Merged with duplicate record 10026.1/6553 on 28.02.2017 by CS (TIS)Abstract Through the Global Gas Flaring Reduction (GGFR) initiative a substantial amount of effort and international attention has been focused on the reduction of gas flaring since 2002 (Elvidge et al., 2009). Nigeria is rated as the second country in the world for gas flaring, after Russia. In an attempt to reduce and eliminate gas flaring the federal government of Nigeria has implemented a number of gas flaring reduction projects, but poor governmental regulatory policies have been mostly unsuccessful in phasing it out. This study examines the effects of pollution from gas flaring using multiple satellite based sensors (Landsat 5 TM and Landsat 7 ETM+) with a focus on vegetation health in the Niger Delta. Over 131 flaring sites in all 9 states (Abia, Akwa Ibom, Bayelsa, Cross Rivers, Delta, Edo, Imo, Ondo and Rivers) of the Niger Delta region have been identified, out of which 11 sites in Rivers State were examined using a case study approach. Land Surface Temperature data were derived using a novel procedure drawing in visible band information to mask out clouds and identify appropriate emissivity values for different land cover types. In 2503 out of 3001 Landsat subscenes analysed, Land Surface Temperature was elevated by at least 1 ℃ within 450 m of the flare. The results from fieldwork, carried out at the Eleme Refinery II Petroleum Company and Onne Flow Station, are compared to the Landsat 5 TM and Landsat 7 ETM+ data. Results indicate that Landsat data can detect gas flares and their associated pollution on vegetation health with acceptable accuracy for both Land Surface Temperature (range: 0.120 to 1.907 K) and Normalized Differential Vegetation Index (sd ± 0.004). Available environmental factors such as size of facility, height of stack, and time were considered. Finally, the assessment of the impact of pollution on a time series analysis (1984 to 2013) of vegetation health shows a decrease in NDVI annually within 120 m from the flare and that the spatio-temporal variability of NDVI for each site is influenced by local factors. This research demonstrated that only 5 % of the variability in δLST and only 12 % of the variability in δNDVI, with distance from the flare stack, could be accounted for by the available variables considered in this study. This suggests that other missing factors (the gas flaring volume and vegetation speciation) play a significant role in the variability in δLST and δNDVI respectively

Numerical modeling of thermal bar and stratification pattern in Lake Ontario using the EFDC model

Thermal bar is an important phenomenon in large, temperate lakes like Lake Ontario. Spring thermal bar formation reduces horizontal mixing, which in turn, inhibits the exchange of nutrients. Evolution of the spring thermal bar through Lake Ontario is simulated using the 3D hydrodynamic model Environmental Fluid Dynamics Code (EFDC). The model is forced with the hourly meteorological data from weather stations around the lake, flow data for Niagara and St. Lawrence rivers, and lake bathymetry. The simulation is performed from April to July, 2011; on a 2-km grid. The numerical model has been calibrated by specifying: appropriate initial temperature and solar radiation attenuation coefficients. The existing evaporation algorithm in EFDC is updated to modified mass transfer approach to ensure correct simulation of evaporation rate and latent heatflux. Reasonable values for mixing coefficients are specified based on sensitivity analyses. The model simulates overall surface temperature profiles well (RMSEs between 1-2°C). The vertical temperature profiles during the lake mixed phase are captured well (RMSEs < 0.5°C), indicating that the model sufficiently replicates the thermal bar evolution process. An update of vertical mixing coefficients is under investigation to improve the summer thermal stratification pattern. Keywords: Hydrodynamics, Thermal BAR, Lake Ontario, GIS