High resolution 3-D imaging via multi-pass SAR

Spaceborne/airborne synthetic aperture radar (SAR) systems provide high resolution two-dimensional terrain imagery. The paper proposes a technique for combining multiple SAR images, acquired on flight paths slightly separated in the elevation direction, to generate high resolution three-dimensional imagery. The technique could be viewed as an extension to interferometric SAR (InSAR) in that it generates topographic imagery with an additional dimension of resolution. The 3-D multi-pass SAR imaging system is typically characterised by a relatively short ambiguity length in the elevation direction. To minimise the associated ambiguities we exploit the relative phase information within the set of images to track the terrain landscape. The SAR images are then coherently combined, via a nonuniform DFT, over a narrow (in elevation) volume centred on the 'dominant' terrain ground plane. The paper includes a detailed description of the technique, background theory, including achievable resolution, and the results of an experimental study.

A Small Baseline DIFSAR Approach for Investigating Deformations on Full Resolution SAR Interferograms

This paper presents a differential synthetic apertureradar (SAR) interferometry (DIFSAR) approach for investigatingdeformation phenomena on full-resolution DIFSAR interferograms.In particular, our algorithm extends the capabilityof the small-baseline subset (SBAS) technique that relies onsmall-baseline DIFSAR interferograms only and is mainly focusedon investigating large-scale deformations with spatial resolutionsof about 100 100 m. The proposed technique is implemented byusing two different sets of data generated at low (multilook data)and full (single-look data) spatial resolution, respectively. Theformer is used to identify and estimate, via the conventional SBAStechnique, large spatial scale deformation patterns, topographicerrors in the available digital elevation model, and possibleatmospheric phase artifacts; the latter allows us to detect, onthe full-resolution residual phase components, structures highlycoherent over time (buildings, rocks, lava, structures, etc.), as wellas their height and displacements. In particular, the estimation ofthe temporal evolution of these local deformations is easily implementedby applying the singular value decomposition technique.The proposed algorithm has been tested with data acquired by theEuropean Remote Sensing satellites relative to the Campania area(Italy) and validated by using geodetic measurements.

The magma budget of Volcan Arenal, Costa Rica from 1980 to 2004

Using topographic data collected by radar interferometry, stereo-photogrammetry, and field survey we have measured the changing surface of Volcan Arenal in Costa Rica over the period from 1980 to 2004. During this time this young volcano has mainly effused basaltic andesite lava, continuing the activity that began in 1968. Explosive products form only a few percent of the volumetric output. We have calculated digital elevation models for the years 1961, 1988 and 1997 and modified existing models for 2000 and 2004. From these we have estimated the volume of lava effused and coupled this with the data presented by an earlier study for 1968-1980. We find that a dense rock equivalent volume of 551 M m(3) was effused from 1968 to 2004. The dense rock equivalent effusion rate fell from about 2 m(3) s(-1) to about 0.1-0.2 m(3) s(-1) over the same period, with an average rate of about 0.5 m(3) s(-1). Between 1980 and 2004, the average effusion rate was 0.36 m(3) s(-1), a similar rate to that measured between 1974 and 1980. There have been two significant deviations from this long-term rate. The effusion rate increased from 1984 to 1991, at the same time as explosivity increased. After a period of moderate effusion rates in the 1990s, the rate fell to lower levels around 1999. (c) 2006 Elsevier B.V. All rights reserved.

Glaciers and climate in northern Sweden during the 19th and 20th century

Our understanding of the climate of northern Sweden during the late Holocene is largely dependent on proxy-data series. These datasets remain spatially and temporally sparse and instrumental series are rare prior to the mid 19th century. Nevertheless, the glaciology and paleo-glaciology of the region has a strong potential significance for the exploration of climate change scenarios, past and future. The aim of this thesis is to investigate the 19th and 20th century climate in the northern Swedish mountain range. This provides a good opportunity to analyse the natural variability of the climate before the onset of the industrial epoch. Developing a temporal understanding of fluctuations in glacier front positions and glacier mass balance that is linked to a better understanding of their interaction and relative significance to climate is fundamental in the assessment of past climate. I have chosen to investigate previously unexplored temperature data from northern Sweden from between 1802 and 1860 and combined it with a temperature series from a synoptic station in Haparanda, which began operation in 1859, in order to create a reliable long temperature series for the period 1802 to 2002. I have also investigated two different glaciers, Pårteglaciären and Salajekna, which are located in different climatic environments. These glaciers have, from a Swedish perspective, long observational records. Furthermore, I have investigated a recurring jökulhlaup at the glacier Sälkaglaciären in order to analyse glacier-climate relationships with respect to the jökulhlaups. A number of datasets are presented, including: glacier frontal changes, in situ and photogrammetric mass balance data, in situ and satellite radar interferometry measurements of surface velocity, radar measurements, ice volume data and a temperature series. All these datasets are analysed in order to investigate the response of the glaciers to climatic stimuli, to attribute specific behaviour to particular climates and to analyse the 19th and 20th century glacier/climate relationships in northern Sweden. The 19th century was characterized by cold conditions in northern Sweden, particularly in winter. Significant changes in the amplitude of the annual temperature cycle are evident. Through the 19th century there is a marked decreasing trend in the amplitude of the data, suggesting a change towards a prevalence of maritime (westerly) air masses, something which has characterised the 20th century. The investigations on Salajekna support the conclusion that the major part of the 19th century was cold and dry. The 19th century advance of Salajekna was probably caused by colder climate in the late 18th and early 19th centuries, coupled with a weakening of the westerly airflow. The investigations on Pårteglaciären show that the glacier has a response time of ~200 years. It also suggests that there was a relatively high frequency of easterly winds providing the glacier with winter precipitation during the 19th century. Glaciers have very different response times and are sensitive to different climatic parameters. Glaciers in rather continental areas of the Subarctic and Arctic can have very long response times because of mass balance considerations and not primarily the glacier dynamics. This is of vital importance for analyzing Arctic and Subarctic glacier behaviour in a global change perspective. It is far from evident that the behaviour of the glacier fronts today reflects the present climate.

Interferometria SAR per lo studio di movimenti e generazione di modelli digitali del terreno in Antartide

Satellite SAR (Synthetic Aperture Radar) interferometry represents a valid technique for digital elevation models (DEM) generation, providing metric accuracy even without ancillary data of good quality. Depending on the situations the interferometric phase could be interpreted both as topography and as a displacement eventually occurred between the two acquisitions. Once that these two components have been separated it is possible to produce a DEM from the first one or a displacement map from the second one. InSAR DEM (Digital Elevation Model) generation in the cryosphere is not a straightforward operation because almost every interferometric pair contains also a displacement component, which, even if small, when interpreted as topography during the phase to height conversion step could introduce huge errors in the final product. Considering a glacier, assuming the linearity of its velocity flux, it is therefore necessary to differentiate at least two pairs in order to isolate the topographic residue only. In case of an ice shelf the displacement component in the interferometric phase is determined not only by the flux of the glacier but also by the different heights of the two tides. As a matter of fact even if the two scenes of the interferometric pair are acquired at the same time of the day only the main terms of the tide disappear in the interferogram, while the other ones, smaller, do not elide themselves completely and so correspond to displacement fringes. Allowing for the availability of tidal gauges (or as an alternative of an accurate tidal model) it is possible to calculate a tidal correction to be applied to the differential interferogram. It is important to be aware that the tidal correction is applicable only knowing the position of the grounding line, which is often a controversial matter. In this thesis it is described the methodology applied for the generation of the DEM of the Drygalski ice tongue in Northern Victoria Land, Antarctica. The displacement has been determined both in an interferometric way and considering the coregistration offsets of the two scenes. A particular attention has been devoted to investigate the importance of the role of some parameters, such as timing annotations and orbits reliability. Results have been validated in a GIS environment by comparison with GPS displacement vectors (displacement map and InSAR DEM) and ICEsat GLAS points (InSAR DEM).

Satellite image map of King George Island, Antarctica

A new topographic database for King George Island, one of the most visited areas in Antarctica, is presented. Data from differential GPS surveys, gained during the summers 1997/98 and 1999/2000, were combined with up to date coastlines from a SPOT satellite image mosaic, and topographic information from maps as well as from the Antarctic Digital Database. A digital terrain model (DTM) was generated using ARC/INFO GIS. From contour lines derived from the DTM and the satellite image mosaic a satellite image map was assembled. Extensive information on data accuracy, the database as well as on the criteria applied to select place names is given in the multilingual map. A lack of accurate topographic information in the eastern part of the island was identified. It was concluded that additional topographic surveying or radar interferometry should be conducted to improve the data quality in this area. In three case studies, the potential applications of the improved topographic database are demonstrated. The first two examples comprise the verification of glacier velocities and the study of glacier retreat from the various input data-sets as well as the use of the DTM for climatological modelling. The last case study focuses on the use of the new digital database as a basic GIS (Geographic Information System) layer for environmental monitoring and management on King George Island.

Desarrollo de Técnicas Interferométricas para Radares de Alta Resolución en Milimétricas

En esta tesis doctoral se aborda el desarrollo de técnicas interferométricas para radares de alta resolución en milimétricas. Específicamente, se centra en el desarrollo de técnicas para radares terrestres y en concreto aquellas relacionadas con el seguimiento, clasificación y generación de imágenes de blancos móviles. Aparte del desarrollo teórico y simulación de las diferentes técnicas, uno de los grandes retos de esta tesis es la verificación experimental de las mismas. Para este propósito, se empleará un prototipo de sensor radar interferométrico de muy alta resolución y de ondas milimétricas. El primer capítulo de la tesis realiza una pequeña introducción a las técnicas interferométricas y presenta los objetivos, motivación y organización del presente trabajo. El segundo capítulo hace una pequeña introducción a los radares interferométricos de alta resolución en milimétricas y presenta el sensor radar experimental con el que se llevarán a cabo la validación de las distintas técnicas presentadas. El tercer capítulo recoge las distintas técnicas interferométricas desarrolladas para el seguimiento y cálculo de la altura de blancos móviles en radares de alta resolución. Entre ellas se encuentran: el seguimiento de blancos móviles en sucesión de imágenes ISAR, el cálculo de altura de blancos en imágenes radar bidimensionales distancia-tiempo y su equivalente en imágenes distancia-Doppler. El cuarto capítulo presenta una aplicación de las técnicas interferométricas para la vigilancia de tráfico en carreteras. Se describirán dos configuraciones radar que permiten calcular la velocidad de todos los blancos iluminados por el radar identificando de manera unívoca a los blancos en función del carril por el que circulen. El quinto capítulo presenta técnicas interferométricas aplicadas a la vigilancia en entornos marítimos basadas en la generación de imágenes interferométricas. Para demostrar la viabilidad del uso de estas imágenes se ha desarrollado un simulador de blancos móviles extensos realistas. Abstract This Ph. D thesis deals with the development of radar interferometry techniques for high-resolution millimeter wave sensors. It focuses on the development of techniques for ground-based radars and specifically those related to monitoring, classification and imaging of moving targets. Apart from the theoretical development and simulation, another major technical challenge of this thesis is the experimental verification of the different techniques. For that purpuse, a very high resolution interferometric millimeter wave radar sensor is used. The first chapter of the thesis makes a brief introduction to the interferometric techniques and shows the goals, motivation and organization of this work. The second chapter provides a brief introduction to high resolution interferometric radars in millimeter waves and presents the experimental radar sensor which will be used for the validation of the various techniques presented. The third chapter presents the different interferometric techniques developed for monitoring and obtaining the height of moving targets in high resolution radars. Among them are: tracking of moving targets in a succession of ISAR images, targets height calculation using bidimensional range-time radar images and the equivalente technique using range-Doppler images. The fourth chapter presents the application of interferometric techniques for road traffic monitoring. Two radar configurations are described. Both of them are able to obtain the speed for simultaneuslly illuminated targets, and univocally identify each target based on the detected road lane. The fifth chapter presents the application of interferometric techniques to maritime surveillance based on interferometric imaging. To demonstrate the feasibility of the presented techniques a realistic simulator of extended moving targets has been developed.

Hydrologic Dynamics of the Ground-Water-Dependent Sian Ka’an Wetlands, Mexico, Derived from InSAR and SAR Data

The 5,280 km2 Sian Ka’an Biosphere Reserve includes pristine wetlands fed by ground water from the karst aquifer of the Yucatan Peninsula, Mexico. The inflow through underground karst structures is hard to observe making it difficult to understand, quantify, and predict the wetland dynamics. Remotely sensed Synthetic Aperture Radar (SAR) amplitude and phase observations offer new opportunities to obtain information on hydrologic dynamics useful for wetland management. Backscatter amplitude of SAR data can be used to map flooding extent. Interferometric processing of the backscattered SAR phase data (InSAR) produces temporal phase-changes that can be related to relative water level changes in vegetated wetlands. We used 56 RADARSAT-1 SAR acquisitions to calculate 38 interferograms and 13 flooding maps with 24 day and 48 day time intervals covering July 2006 to March 2008. Flooding extent varied between 1,067 km2 and 2,588 km2 during the study period, and main water input was seen to take place in sloughs during October–December. We propose that main water input areas are associated with water-filled faults that transport ground water from the catchment to the wetlands. InSAR and Landsat data revealed local-scale water divides and surface water flow directions within the wetlands.

Mass balances of different sectors in Antarctica in 1996, 2000 and 2006

Large uncertainties remain in the current and future contribution to sea level rise from Antarctica. Climate warming may increase snowfall in the continent's interior, but enhance glacier discharge at the coast where warmer air and ocean temperatures erode the buttressing ice shelves. Here, we use satellite interferometric synthetic-aperture radar observations from 1992 to 2006 covering 85% of Antarctica's coastline to estimate the total mass flux into the ocean. We compare the mass fluxes from large drainage basin units with interior snow accumulation calculated from a regional atmospheric climate model for 1980 to 2004. In East Antarctica, small glacier losses in Wilkes Land and glacier gains at the mouths of the Filchner and Ross ice shelves combine to a near-zero loss of 4 ± 61 Gt/yr. In West Antarctica, widespread losses along the Bellingshausen and Amundsen seas increased the ice sheet loss by 59% in 10 years to reach 132 ± 60 Gt/yr in 2006. In the Peninsula, losses increased by 140% to reach 60 ± 46 Gt/yr in 2006. Losses are concentrated along narrow channels occupied by outlet glaciers and are caused by ongoing and past glacier acceleration. Changes in glacier flow therefore have a significant, if not dominant impact on ice sheet mass balance.

(Table 1) Mass loss of the Greenland Ice Sheet

Using satellite radar interferometry observations of Greenland, we detected widespread glacier acceleration below 66° north between 1996 and 2000, which rapidly expanded to 70° north in 2005. Accelerated ice discharge in the west and particularly in the east doubled the ice sheet mass deficit in the last decade from 90 to 220 cubic kilometers per year. As more glaciers accelerate farther north, the contribution of Greenland to sea-level rise will continue to increase.

Three-dimensional space-borne synthetic aperture radar (SAR) imaging with multiple pass processing

Three-dimensional (3D) synthetic aperture radar (SAR) imaging via multiple-pass processing is an extension of interferometric SAR imaging. It exploits more than two flight passes to achieve a desired resolution in elevation. In this paper, a novel approach is developed to reconstruct a 3D space-borne SAR image with multiple-pass processing. It involves image registration, phase correction and elevational imaging. An image model matching is developed for multiple image registration, an eigenvector method is proposed for the phase correction and the elevational imaging is conducted using a Fourier transform or a super-resolution method for enhancement of elevational resolution. 3D SAR images are obtained by processing simulated data and real data from the first European Remote Sensing satellite (ERS-1) with the proposed approaches.

Sistemas radar aplicados a la investigación de subsidencia y movimientos de ladera

En este trabajo se presentan los sistemas radar satélite y terrestres, así como los métodos de análisis de imágenes radar clásicos y avanzados para la investigación de los movimientos del terreno, haciendo énfasis en la subsidencia y los movimientos de ladera. Para ello en primer lugar se describen los distintos sensores radar disponibles así como las principales características de las imágenes radar generadas. A continuación se detallan los aspectos fundamentales de la interferometría diferencial, de los distintos métodos de interferometría diferencial avanzada y del radar terrestre. Finalmente se presentan los resultados obtenidos en distintas zonas de estudio: la subsidencia por explotación del acuífero en el área metropolitana de Murcia, la subsidencia minera y los movimientos de ladera de la Sierra de Cartagena, los movimientos de ladera de la cuenca del río Gállego y el deslizamiento del Portalet.