Investigation of Sea Ice Using Multiple Synthetic Aperture Radar Acquisitions

The papers of this thesis are not available in Munin. Paper I: Yitayew, T. G., Ferro-Famil, L., Eltoft, T. & Tebaldini, S. (2017). Tomographic imaging of fjord ice using a very high resolution ground-based SAR system. Available in IEEE Transactions on Geoscience and Remote Sensing, 55 (2):698-714. Paper II: Yitayew, T. G., Ferro-Famil, L., Eltoft, T. & Tebaldini, S. (2017). Lake and fjord ice imaging using a multifrequency ground-based tomographic SAR system. Available in IEEE Journal of Selected Topics in Applied Earth Observations and Remote Sensing, 10(10):4457-4468. Paper III: Yitayew, T. G., Divine, D. V., Dierking, W., Eltoft, T., Ferro-Famil, L., Rosel, A. & Negrel, J. Validation of Sea ice Topographic Heights Derived from TanDEMX Interferometric SAR Data with Results from Laser Profiler and Photogrammetry. (Manuscript).The thesis investigates imaging in the vertical direction of different types of ice in the arctic using synthetic aperture radar (SAR) tomography and SAR interferometry. In the first part, the magnitude and the positions of the dominant scattering contributions within snow covered fjord and lake ice layers are effectively identified by using a very high resolution ground-based tomographic SAR system. Datasets collected at multiple frequencies and polarizations over two test sites in Tromsø area, northern Norway, are used for characterizing the three-dimensional response of snow and ice. The presented experimental results helped to improve our understanding of the interaction between radar waves and snow and ice layers. The reconstructed radar responses are also used for estimating the refractive indices and the vertical positions of the different sub-layers of snow and ice. The second part of the thesis deals with the retrieval of the surface topography of multi-year sea ice using SAR interferometry. Satellite acquisitions from TanDEM-X over the Svalbard area are used for analysis. The retrieved surface height is validated by using overlapping helicopter-based stereo camera and laser profiler measurements, and a very good agreement has been found. The work contributes to an improved understanding regarding the potential of SAR tomography for imaging the vertical scattering distribution of snow and ice layers, and for studying the influence of both sensor parameters such as its frequency and polarization and scene properties such as layer stratification, air bubbles and small-scale roughness of the interfaces on snow and ice backscattered signal. Moreover, the presented results reveal the potential of SAR interferometry for retrieving the surface topography of sea ice

Modeling of Subsurface Scattering from Ice Sheets for Pol-InSAR Applications

Remote sensing is a fundamental tool to measure the dynamics of ice sheets and provides valuable information for ice sheet projections under a changing climate. There is, however, the potential to further reduce the uncertainties in these projections by developing innovative remote sensing methods. One of these remote sensing techniques, the polarimetric synthetic aperture radar interferometry (Pol-InSAR), is known since decades to have the potential to assess the geophysical properties below the surface of ice sheets, because of the penetration of microwave signals into dry snow, firn, and ice. Despite this, only very few studies have addressed this topic and the development of robust Pol-InSAR applications is at an early stage. Two potential Pol-InSAR applications are identified as the motivation for this thesis. First, the estimation and compensation of the penetration bias in digital elevation models derived with SAR interferometry. This bias can lead to errors of several meters or even tens of meters in surface elevation measurements. Second, the estimation of geophysical properties of the subsurface of glaciers and ice sheets using Pol-InSAR techniques. There is indeed potential to derive information about melt-refreeze processes within the firn, which are related to density and affect the mass balance. Such Pol-InSAR applications can be a valuable information source with the potential for monthly ice sheet wide coverage and high spatial resolution provided by the next generation of SAR satellites. However, the required models to link the Pol-InSAR measurements to the subsurface properties are not yet established. The aim of this thesis is to improve the modeling of the vertical backscattering distribution in the subsurface of ice sheets and its effect on polarimetric interferometric SAR measurements at different frequencies. In order to achieve this, polarimetric interferometric multi-baseline SAR data at different frequencies and from two different test sites on the Greenland ice sheet are investigated. This thesis contributes with three concepts to a better understanding and to a more accurate modeling of the vertical backscattering distribution in the subsurface of ice sheets. First, the integration of scattering from distinct subsurface layers. These are formed by refrozen melt water in the upper percolation zone and cause an interesting coherence undulation pattern, which cannot be explained with previously existing models. This represents a first link between Pol-InSAR data and geophysical subsurface properties. The second step is the improved modeling of the general vertical backscattering distribution of the subsurface volume. The advantages of more flexible volume models are demonstrated, but interestingly, the simple modification of a previously existing model with a vertical shift parameter lead to the best agreement between model and data. The third contribution is the model based compensation of the penetration bias, which is experimentally validated. At the investigated test sites, it becomes evident that the model based estimates of the surface elevations are more accurate than the interferometric phase center locations, which are conventionally used to derive surface elevations of ice sheets. This thesis therefore improves the state of the art of subsurface scattering modeling for Pol-InSAR applications, demonstrates the model-based penetration bias compensation, and makes a further research step towards the retrieval of geophysical subsurface information with Pol-InSAR

Sea ice local surface topography from single-pass satellite InSAR measurements: a feasibility study

Quantitative parameters characterizing the sea ice surface topography are needed in geophysical investigations such as studies on atmosphere–ice interactions or sea ice mechanics.Recently, the use of space-borne single-pass interferometric synthetic aperture radar (InSAR) for retrieving the ice surface topography has attracted notice among geophysicists. In this paper the potential of InSAR measurements is examined for several satellite configurations and radar frequencies, considering statistics of heights and widths of ice ridges as well as possible magnitudes of ice drift. It is shown that, theoretically, surface height variations can be retrievedwith relative errors < 0.5 m. In practice, however, the sea ice drift and open water leads may contribute significantly to the measured interferometric phase. Another essential factor is the dependence of the achievable interferometric baseline on the satellite orbit configurations. Possibilities to assess the influence of different factors on the measurement accuracy are demonstrated: signal-to-noise ratio, presence of a snow layer, and the penetration depth into the ice. Practical examples of sea surface height retrievals from bistatic SAR images collectedduring the TanDEM-X Science Phase are presented

Ice Shelf Melt Rates and 3D Imaging

Ice shelves are sensitive indicators of climate change and play a critical role in the stability of ice sheets and oceanic currents. Basal melting of ice shelves plays an important role in both the mass balance of the ice sheet and the global climate system. Airborne- and satellite based remote sensing systems can perform thickness measurements of ice shelves. Time separated repeat flight tracks over ice shelves of interest generate data sets that can be used to derive basal melt rates using traditional glaciological techniques. Many previous melt rate studies have relied on surface elevation data gathered by airborne- and satellite based altimeters. These systems infer melt rates by assuming hydrostatic equilibrium, an assumption that may not be accurate, especially near an ice shelf’s grounding line. Moderate bandwidth, VHF, ice penetrating radar has been used to measure ice shelf profiles with relatively coarse resolution. This study presents the application of an ultra wide bandwidth (UWB), UHF, ice penetrating radar to obtain finer resolution data on the ice shelves. These data reveal significant details about the basal interface, including the locations and depth of bottom crevasses and deviations from hydrostatic equilibrium. While our single channel radar provides new insight into ice shelf structure, it only images a small swatch of the shelf, which is assumed to be an average of the total shelf behavior. This study takes an additional step by investigating the application of a 3D imaging technique to a data set collected using a ground based multi channel version of the UWB radar. The intent is to show that the UWB radar could be capable of providing a wider swath 3D image of an ice shelf. The 3D images can then be used to obtain a more complete estimate of the bottom melt rates of ice shelves

Five decades of radioglaciology

Radar sounding is a powerful geophysical approach for characterizing the subsurface conditions of terrestrial and planetary ice masses at local to global scales. As a result, a wide array of orbital, airborne, ground-based, and in situ instruments, platforms and data analysis approaches for radioglaciology have been developed, applied or proposed. Terrestrially, airborne radar sounding has been used in glaciology to observe ice thickness, basal topography and englacial layers for five decades. More recently, radar sounding data have also been exploited to estimate the extent and configuration of subglacial water, the geometry of subglacial bedforms and the subglacial and englacial thermal states of ice sheets. Planetary radar sounders have observed, or are planned to observe, the subsurfaces and near-surfaces of Mars, Earth's Moon, comets and the icy moons of Jupiter. In this review paper, and the thematic issue of the Annals of Glaciology on ‘Five decades of radioglaciology’ to which it belongs, we present recent advances in the fields of radar systems, missions, signal processing, data analysis, modeling and scientific interpretation. Our review presents progress in these fields since the last radio-glaciological Annals of Glaciology issue of 2014, the context of their history and future prospects

Ice in Norwegian subarctic fjords and coastal regions: An examination of ice formation, properties, and trends based on remote sensing and in situ data

Larger fjords along the Norwegian coast are ice free all year due to the influence of warm Atlantic water, however, sea ice can form in the inner parts of fjords and in smaller fjord branches. While a wide breadth of work exists examining mainland Norwegian fjords often focused on water mass dynamics and their link to biological processes, little research has been conducted on specifically sea ice in these locations. The overarching goal of the work presented here is to address this gap in research and knowledge. To do so, first an assessment of ice extent in fjords and other coastal areas along the coast of mainland Norway from 2001 to 2019 between February through May is completed. Through the development of an automated method to estimate quantitatively ice extent in optical satellite imagery, variability in ice extent between years as well as within single seasons is highlighted for regions and specific fjords. An investigation into the factors potentially driving ice formation is conducted with focused placed on air temperature, snowfall, and rainfall plus snowmelt. Next, seven fjords located in northern Norway are studied in-depth through the collection of ice and ocean measurements over the course of three winter seasons spanning 2017 into 2020. Ice stratigraphy, bulk ice salinity and oxygen isotopic composition, ocean temperature, salinity and oxygen isotopic composition, and river isotopic composition are analyzed. Possible connections to temperature and snowfall are revisited in addition to river runoff illuminating further variations not only in ice extent between seasons but also ice properties and the factors driving ice formation. Lastly, bulk ice salinity and oxygen isotopic composition are given a closer look through the inversion of previously derived, separate, relationships of each to growth rate and interface water composition. Results provide a promising method to deduce ocean and weather conditions during ice growth in fjords and coastal areas where continuous measurement through winter is not possible

Probabilistic graphical techniques for automated ice-bottom tracking and comparison between state-of-the-art solutions

We present improvements to existing frameworks for automated extraction of ice interfaces applied to two-dimensional and three-dimensional radar echograms of polar ice sheets. These improvements consist of novel image pre-processing steps and empirically-derived cost functions that allow for the integration of further domain-specific knowledge into the models employed. Along with an explanation of our modifications, we demonstrate the results obtained by our proposed models and algorithms, such as a 43% decrease in mean tracking error in the case of three-dimensional imagery. We also present the results obtained by several state-of-the-art ice-interface tracking solutions, and compare all automated results with manually-corrected ground-truth data. Furthermore, we perform a self-assessment of tracking results by analyzing the differences found between the automatically extracted ice-layers in cases where two separate radar measurements have been made at the same location