Monitoring Soil Moisture and Freeze/Thaw State Using C-band Imaging Radar

Soil moisture is an important state variable in many hydrological and meteorological applications. This thesis explores the use of the C-band synthetic aperture radar (SAR) parameters to monitor soil moisture and freeze/thaw state in a cold-season hydrologic environment. The circular-linear compact polarimetric (CP) configuration is considered as a possible alternative of the quad polarimetric (QP) system because it acquires images with wider swath and reduced complexity, cost and energy requirement of the radar system while maintaining the information content of the acquired imagery. In this study, 15 RADARSAT-2 QP images were acquired from October 2013 to June 2014 and CP images were simulated from each RADARSAT-2 QP imagery acquired. Field measurements of soil properties were collected along with the radar imagery acquisitions. The backscattering coefficients in all polarizations were able to discriminate frozen and unfrozen soils. But their correlations with soil moisture content were weak if examining frozen or unfrozen soils separately. The Oh et al. (1992) model was implemented in this study to compare with acquired RADARSAT-2 data. A good agreement was found between the linear polarimetric backscattering coefficients simulated by the Oh model and the RADARSAT-2 data, indicating that the study site even with 10 cm tall standing hay was consistent with a bare soil site at C-band and the Oh model can be applied to frozen soils. With respect to CP parameters, the first and fourth Stokes parameters and m-δ surface and volume scattering components can detect soil freeze/thaw state and have potential for frozen/unfrozen soils mapping. The influence of vegetation on selected CP parameters was also evident in this study. Results demonstrated the utility of C-band radar in detecting soil freeze/thaw state rather than monitoring the changes of soil moisture content. More image acquisitions during the freezing and thawing periods, continuous field measurements of soil moisture and state, and ground measurements collected over wider study area can help further develop understanding of the CP parameters and facilitate future use of the CP mode. The contribution of this thesis is to provide better understanding of the CP parameters at a specific site and to demonstrate that CP parameters can replicate QP SAR variables to detect surface soil conditions

Détection des cycles de gel/dégel de la couche active du sol en toundra arctique à partir d’imageries radar à synthèse d’ouverture (RSO) multicapteur en bande C

L’augmentation de la température de l’air moyenne annuelle, chiffrée à +2,3 °C pour les régions de l’arctique Canadien entre 1948 et 2016, a des impacts considérables sur le couvert nival arctique et sur la végétation en place. Ces deux paramètres influencent le régime thermique du sol et donc, les cycles de gel/dégel de sa couche active dans l’écosystème arctique. L’importance du suivi de ces cycles réside dans leur influence sur plusieurs paramètres de la cryosphère tels que le cycle hydrologique et du carbone, la saison de croissance de la végétation, l’état du pergélisol sous-jacent ainsi que l’épaisseur de sa couche active. L’utilisation de données ponctuelles ou provenant de capteurs micro-onde passive à basse résolution présente un enjeu pour le suivi spatial et temporel de ces cycles. Le projet vise à développer un algorithme de détection des cycles de gel/dégel du sol en toundra arctique à partir d’imageries RSO multicapteur (i.e., Sentinel-1 et RADARSAT-2) ayant une couverture temporelle quasi journalière en bande C, afin d’évaluer l’impact de la variabilité spatiale et temporelle des paramètres influençant le régime thermique du sol tel que, les écosystèmes terrestres (i.e., écotype) et la présence de neige. L’étude se concentre sur une zone à l’intérieur du bassin versant du lac Greiner à proximité de la ville de Cambridge Bay au Nunavut. La normalisation de l’angle d’incidence a permis de diminuer le bruit dans les séries temporelles ainsi que de rendre possible l’utilisation d'images acquises à l'intérieur de plusieurs orbites d’observation. Cela a aussi permis d’uniformiser les données des deux capteurs pour les combiner en une seule série temporelle. Deux algorithmes de détections ont été utilisés, soit un algorithme de seuil saisonnier (STA) ainsi qu’un algorithme de détection de changement (CPD). La validation s’est faite à partir des données spatialement distribuées de température du sol et de l’air indépendamment sous forme de précision (%) et de délai (#jours) de détection. Les deux algorithmes ont permis d’obtenir une précision de détection de plus de 97% sur les sites de référence. Une spatialisation, pixel par pixel, de la méthode STA a permis la création de cartes de jour de gel/dégel pour le site d’étude. La combinaison des cartes de jour de transition avec la carte d’écotype a permis de modéliser l’impact des caractéristiques des écotypes sur le jour de transition. Les résultats obtenus dans ce projet démontrent clairement le potentiel de l’utilisation des données RSO en bande C pour la détection des cycles de gel/dégel, ce qui constitue un résultat important en raison de la quantité grandissante de données à cette fréquence (e.g., RCM, Sentinel-1A-C-D). La méthode présentée dans ce projet pourrait permettre de créer des cartes de transition pour tout le bassin versant du lac Greiner à partir de données RSO en bande C.Abstract : The observed average annual surface temperature increase of 2.3°C in the Canadian Arctic regions between 1948 and 2016 has significant effects on the Arctic snow cover and on the vegetation in place. Those two parameters influence the thermal regime of the ground and therefore the freeze and thaw (F/T) cycles of the soil active layer in the Arctic tundra ecosystem. The importance of monitoring these cycles lies in their influence on several parameters of the cryosphere such as the hydrological and carbon cycle, the vegetation growing season, the state of the underlying permafrost and the thickness of its active layer. The use of punctual data or low-resolution passive microwave sensors presents a challenge for the spatial and temporal monitoring of these cycles. The project aims to develop an algorithm for soil freeze/thaw cycles detection in arctic tundra from multisensor C-band imagery (i.e., Sentinel-1 and RADARSAT-2) to assess the impact of the spatial and temporal variability of the parameters influencing the thermal regime of the ground, such as the terrestrial ecosystems (i.e., ecotype) and the snow cover. The study focused on a region of the Greiner lake watershed on Victoria Island in Nunavut. An incidence angle normalization was applied to the backscatter time series to remove influence of the acquisition angle on backscatter and to allow for the use of images acquired within several orbits of observation. This also standardized the data from the two sensors to combine them into a single time series. Two detection algorithms were used on the normalized backscatter coefficient data, namely a seasonal threshold algorithm (STA) and a change point detection algorithm (CPD). A spatially distributed network of soil and air temperature were used for validation in the form of accuracy (%) and delay (#days) of detection. Both algorithms achieved a detection accuracy of more than 97% for the entire analysis period on the reference sites. A pixel-by-pixel spatialization of the STA method allowed to create F/T transition maps for the extended study site. The combination of the transition maps with the ecotype data made it possible to model the impact of ecotype characteristics on the day of transition. The results obtained in this project clearly demonstrate the potential of using C-band for the detection of F/T cycles, which is an important aspect due to the increasing number of data at this frequency (e.g., RCM, Sentinel -1A-C-D). The method presented in this project could then make it possible to create transition maps for the entire Greiner Lake watershed from C-band SAR data and thus improve the integration of this parameter in climate models

Evaluation of the potential of ALOS PALSAR L-band quadpol radar data for the retrieval of growing stock volume in Siberia

Because of the massive wood trade, illegal logging and severe damages due to fires, insects and pollution, it is necessary to monitor Siberian forests on a large-scale, frequently and accurately. One possible solution is to use synthetic aperture radar (SAR) remote sensing technique, in particular by combining polarimetric technique. In order to evaluate the potentiality of ALOS PALSAR L-band full polarimetric radar for estimation of GSV, a number of polarimetric parameters are investigated to characterise the polarisation response of forest cover. Regardless of the weather conditions, a high correlation (R=-0.87) is achieved between polarimetric coherence and GSV. The coherence in sparse forest is always higher than in dense forest. The coherence level and the dynamic range strongly depends on the weather conditions. The four-component polarimetric decomposition method has been applied to the ALOS PALSAR L-band data to compare the decomposition powers with forest growing stock volume (GSV). Double-bounce and volume scattering powers show significant correlation with GSV. The correlation between polarimetric decomposition parameters and GSV is enhanced if the ratio of ground-to-volume scattering is used instead of considering polarimetric decomposition powers separately. Two empirical models have been developed that describe the ALOS PALSAR L-band polarimetric coherence and ground-to-volume scattering ratio as a function of GSV. The models are inverted to retrieve the GSV for Siberian forests. The best RMSE of 38 m³/ha and R²=0.73 is obtained based on polarimetric coherence. On the other hand, using the ratio of ground-to-volume scattering the best retrieval accuracy of 44 m³/ha and R²=0.62 is achieved. The best retrieval results for both cases are observed under unfrozen condition. Saturation effects for estimated GSV versus ground-truth GSV are not observed up to 250 m³/ha

Understanding of crop lodging induced changes in scattering mechanisms using RADERSAT-2 and Sentinel-1 derived metrics

Abstract. Crop lodging – the bending of crop stems from the vertical – is a major yield-reducing factor in cereal crops and causes deterioration in grain quality. Accurate assessment of crop lodging is important for improving estimates of crop yield losses, informing insurance loss adjusters and influencing management decisions for subsequent seasons. The role of remote sensing data, particularly synthetic aperture radar (SAR) data has been emphasized in the recent literature for crop lodging assessment. However, the effect of lodging on SAR scattering mechanisms is still unknown. Therefore, this research aims to understand the possible change in scattering mechanisms due to lodging by investigating SAR image pairs before and after lodging. We conducted the study in 26 wheat fields in the Bonifiche Ferraresi farm, located in Jolanda di Savoia, Ferrara, Italy. We measured temporal crop biophysical (e.g. crop angle) parameters and acquired multi-incidence angle RADARSAT-2 (R-2 FQ8-27° and R-2 FQ21-41°) and Sentinel-1 (S-1 40°) images corresponding to the time of field observations. We extracted metrics of SAR scattering mechanisms from RADARSAT-2 and Sentinel-1 image pairs in different zones using the unsupervised H/α decomposition algorithm and Wishart classifier. Contrasting results were obtained at different incidence angles. Bragg surface scattering increased in the case of S-1 (6.8%), R-2 FQ8 (1.8%) while at R-2 FQ21, it decreased (8%) after lodging. The change in double bounce scattering was more prominent at low incidence angle. These observations can guide future use of SAR-based information for operational crop lodging assessment in particular, and sustainable agriculture in general

Detection of soil permittivity and soil freezing using satellite microwave radars

Remote sensing of soil permittivity and soil freezing was investigated using two different satellite based microwave radars: ASCAT and ASAR. ASCAT is a scatterometer with a good temporal resolution but coarse spatial resolution. ASAR is a synthetic aperture radar and has fine spatial resolution, but lacks good temporal coverage. Soil permittivity is related to soil moisture, which is considered an essential climate vari- able since it has an effect on both weather and climate. Soil freezing affects hydrological and carbon cycles, surface energy balance, photosynthesis of vegetation and the activity of soil microbes. A semi-empirical model for backscattering of forested land was used to acquire soil permittivity retrievals from satellite measurements using the method of least squares. The onset of soil freezing was determined from the permittivity retrievals using a simple threshold method. A five year time series of satellite observations from July 2007 to June 2012 (April 2012 for ASAR) was investigated in Sodankylä in Northern Finland. The satellite based retrievals were compared against in situ measurements of soil permittivity, soil temperature, soil frost and snow depth. According to the results the satellite permittivity retrievals correlate with each other, but not with in situ permittivity measurements. ASCAT retrieval shows some correlation with in situ temperature measurements, which could impair its correlation with in situ permittivity. The explanation for this phenomenon needs further research. Comparison of soil freezing onset dates from satellite retrievals with in situ soil temperature and soil frost measurements showed quite good agreement for most years, and did not seem to be affected by first snowfall, even though the permittivity retrievals appeared to react in a similar way to snow cover and soil freezing. This indicates that with better calibration of the permittivity threshold limit this method could be used for soil freeze detection. Auxiliary information about air temperature and snow cover could also be used to filter out possible false estimates before freezing and after the snow cover starts to affect the satellite retrievals

Quantitative Estimation of Surface Soil Moisture in Agricultural Landscapes using Spaceborne Synthetic Aperture Radar Imaging at Different Frequencies and Polarizations

Soil moisture and its distribution in space and time plays an important role in the surface energy balance at the soil-atmosphere interface. It is a key variable influencing the partitioning of solar energy into latent and sensible heat flux as well as the partitioning of precipitation into runoff and percolation. Due to their large spatial variability, estimation of spatial patterns of soil moisture from field measurements is difficult and not feasible for large scale analyses. In the past decades, Synthetic Aperture Radar (SAR) remote sensing has proven its potential to quantitatively estimate near surface soil moisture at high spatial resolutions. Since the knowledge of the basic SAR concepts is important to understand the impact of different natural terrain features on the quantitative estimation of soil moisture and other surface parameters, the fundamental principles of synthetic aperture radar imaging are discussed. Also the two spaceborne SAR missions whose data was used in this study, the ENVISAT of the European Space Agency (ESA) and the ALOS of the Japanese Aerospace Exploration Agency (JAXA), are introduced. Subsequently, the two essential surface properties in the field of radar remote sensing, surface soil moisture and surface roughness are defined, and the established methods of their measurement are described. The in situ data used in this study, as well as the research area, the River Rur catchment, with the individual test sites where the data was collected between 2007 and 2010, are specified. On this basis, the important scattering theories in radar polarimetry are discussed and their application is demonstrated using novel polarimetric ALOS/PALSAR data. A critical review of different classical approaches to invert soil moisture from SAR imaging is provided. Five prevalent models have been chosen with the aim to provide an overview of the evolution of ideas and techniques in the field of soil moisture estimation from active microwave data. As the core of this work, a new semi-empirical model for the inversion of surface soil moisture from dual polarimetric L-band SAR data is introduced. This novel approach utilizes advanced polarimetric decomposition techniques to correct for the disturbing effects from surface roughness and vegetation on the soil moisture retrieval without the use of a priori knowledge. The land use specific algorithms for bare soil, grassland, sugar beet, and winter wheat allow quantitative estimations with accuracies in the order of 4 Vol.-%. Application of remotely sensed soil moisture patterns is demonstrated on the basis of mesoscale SAR data by investigating the variability of soil moisture patterns at different spatial scales ranging from field scale to catchment scale. The results show that the variability of surface soil moisture decreases with increasing wetness states at all scales. Finally, the conclusions from this dissertational research are summarized and future perspectives on how to extend the proposed model by means of improved ground based measurements and upcoming advances in sensor technology are discussed. The results obtained in this thesis lead to the conclusion that state-of-the-art spaceborne dual polarimetric L-band SAR systems are not only suitable to accurately retrieve surface soil moisture contents of bare as well as of vegetated agricultural fields and grassland, but for the first time also allow investigating within-field spatial heterogeneities from space

Study of the speckle noise effects over the eigen decomposition of polarimetric SAR data: a review

This paper is focused on considering the effects of speckle noise on the eigen decomposition of the co- herency matrix. Based on a perturbation analysis of the matrix, it is possible to obtain an analytical expression for the mean value of the eigenvalues and the eigenvectors, as well as for the Entropy, the Anisotroopy and the dif- ferent a angles. The analytical expressions are compared against simulated polarimetric SAR data, demonstrating the correctness of the different expressions.Peer ReviewedPostprint (published version

Comparison of Target Detectors to Identify Icebergs in Quad-Polarimetric L-Band Synthetic Aperture Radar Data

Icebergs represent hazards to ships and maritime activities and therefore their detection is essential. Synthetic Aperture Radar (SAR) satellites are very useful for this, due to their capability to acquire data under cloud cover and during day and night passes. In this work, we compared six state-of-the-art polarimetric target detectors to test their performance and ability to detect small-sized icebergs 120 m, as they are easier to detect). However, the differences between quad- and dual- or single-polarimetric detectors became much more evident when the PF value was fixed to low detection probabilities 10−6 (i.e., smaller icebergs). In the single-polarimetric mode, the HV channel showed PD values of 0.62 for open ocean and 0.26 for sea ice, compared to values of 0.81 (open ocean) and 0.77 (sea ice) obtained with quad-polarimetric detectors